有限元法在岩土工程中的实际应用.pptVIP

如果分层足够多， 也可以接近曲线 也可以不作处理 修正：v’=v*s’/s 位移修正是以位移瞬时完成为前提 固结问题不必修正！！！ 逐级施工，对每一级荷载增量，荷载是一次性施加；因此，对每级荷载下位移增量修正后累加； 如果分层无穷多，顶面位移为0； 逐级施工两种模拟方法： 1、假定初始应力，形成网格，单元数增加，求[D]； 2、不形成网格，只作为荷载； 三、开挖 开挖，就是使开挖面上的应力解除， 1．假定开挖面应力为挖去土的自重应力 第Ⅰ层， 在斜面上，由бz和бx两分量叠加，求开挖面上与自重应力等效的结点荷载。 将等大反向的等效结点荷载作用于开挖面上个结点，作有限元计算。 第Ⅱ、Ⅲ层解除的自重应力，分别为γh2、γh3。 2．除第一层外，其它各层用有限元算得的应力:界面应力内插法 第一层开挖后（荷载近似为rz），有限元可算得第二开挖面上的应力。以结点M为例，周围单元1、2、3、4的应力由有限元求得。 假定： 用四个单元形心坐标代入，分别得出四单元的计算应力，有四个方程可得a1、a2、a3、a4，代入上式，用M点坐标代入，求得M的б，即由单元应力内插得开挖面上的应力，进而求得结点力。 由于第一层开挖，第二开挖面上的应力已不是rz，用有限元算得的应力更合理。 如бx 对边界点C或D，可由其下面的结点（E、M点）线性外插： 图2－21（a） 图2－21（b） 3．由有限元解得的位移求开挖面上的结点力 （1）由各开挖面上的初始自重应力求相应结点力 {FⅠ}、{FⅡ}、{FⅢ} {RI} {RⅡ} {RⅡI} （实际上求得所有结点的结点力，但仅取用开挖面结点力保留后用） （2）挖第一层 {RⅠ}由{FⅠ}有累加而得 [K]{Δδ}=－{RⅠ} 由第II开挖面以上的各单元，计算相应由第一层开挖引起的结点力 {ΔFⅡ}1= 实际求得所用结点的结点力 仅对开挖面上的结点求结点荷载{ΔRⅡ}1={ΔFⅡ} 1 同理求得{ΔFⅢ}1 则： {RⅡ}＇={RⅡ}+Δ{RⅡ}1 {RⅢ}＇={RⅢ}+Δ{RⅢ}1 1、切线迭代； 首先，根据假定的初始应力求Et、Vt 或者直接由弹塑性模型求[D] [K] 位移、应变、应力 L次迭代 （1）{б}L-1→Et、Vt→[D]L （2）[D]L→[K]L （3）{б}L-1→{F}L-1→{R}L-1 （4）{DR}L={R}L-{R}L-1 （5）[K]{Dδ}L={DR}L （6）{Dδ}L→{Dε}L→εL=εL-1+ DεL （7） εL 由非线性关系 бL （二）、余量迭代法 2、常劲度迭代 直接用初始模量、或直接假定的弹性常数 L次迭代 （1）{б}L-1→Et、Vt→[D]L （2）[D]L→[K]L （3）{б}L-1→{F}L-1→{R}L-1 （4）{DR}L={R}L-{R}L-1 （5）[K]{Dδ}L={DR}L （6）{Dδ}L→{Dε}L→εL=εL-1+ DεL （7） εL 由非线性关系 бL （二）、余量迭代法 将非线性的σ～ε关系用初应力{σ0}的变化来实现。 （三）、初应力迭代法 用初应力法迭代求解非线性方程组实际上就是设想有一个初应力{σ0}，使其非线性的应力应变关系 与线性的应力应变关系 ＋{σ0}等效（如图）。其中， 为弹性应力。事实上，前一种关系式是真实的，而后一种线性关系式是假想的，只有将{σ0}调整到适当值时，两者才能等效。调整{σ0}使得非线性方程组得到求解的方法即初应力迭代法。 （三）、初应力迭代法 将非线性的σ～