基于非线性有限元分析的边坡失稳过程模拟.docx

基于非线性有限元分析的边坡失稳过程模拟 1 干扰能量法 在工程结构的稳定分析中，以往的方法是静力分析。评估物体是否处于天平状态的依据是通过测量物体的阻力和滑动力之比的大小来确定的。由于力是矢量,应用静力分析法在推导公式和计算时,必须遵循矢量法则,其过程相当繁琐,尤其对三维问题更为复杂。此外,该法还需事先人为确定最危险滑动面,而最危险滑动面的确定具有相当大的难度。干扰能量法是根据系统的能量准则,以考察体受扰动后产生的干扰能量值作为判据去确定考察体是否稳定的一种方法。由于能量是标量,无需考虑方向性问题,所以推导和运算极为简便。此外,该法还可应用代数叠加方法迅速得到考察体内任意一部分,甚至整体的干扰能量值,从而判断考察体是否稳定,并找出考察体内局部稳定性薄弱部位。 结合长江下游某岸坡工程,对边坡逐级进行了土体非线性有限元计算,得到边坡的应力应变状态,并在此基础上,分别进行了相应的干扰能量计算,对边坡的失稳过程进行了探讨。 2 干扰能量法的基本原则 2.1 干扰位移引起的稳定性 工程稳定问题的力学机理与变形体的稳定性类似,对于任一变形体,假设系统总势能为Πp,由文献知,可用总势能的二阶变分来判别系统的稳定状态:当δ2Πp0时,稳定状态;当δ2Πp=0时,临界状态;当δ2Πp0时,失去稳定。 设Πp(ui0)为考察体在原有平衡位置的总势能,Πp(ui0+δui)为考察体受到δui干扰位移后的总势能,ΔΠp则为相应于干扰位移引起的干扰能量,由台劳基数展开: 因扰动前的总势能具有驻值,即δΠp=0,故有 ΔΠp≈21δ2Πp(2) ΔΠp为受干扰位移后系统势能的改变量,所以为一种干扰能量。假设存在某种形式的干扰位移使得ΔΠp0,则可以得出系统处于不稳定状态的结论;但是若存在某种形式干扰位移使得ΔΠp0,并不能得到系统处于稳定状态的结论。须寻找最不利干扰位移(使ΔΠp最小),然后计算出最小干扰能量,对系统进行稳定判别。 2.2 单元能量干扰能量计算 考虑计入几何非线性效应,可推得干扰能量的有限元格式为 分别为原有状态的应力和荷载;和分别为干扰位移增量、应变增量、小变形应变增量及应力增量;ΔU表示干扰后贮存在物体内的变形能增量,它使系统恢复原位,称之为干扰内能;ΔW表示外力在干扰位移中所作的功,它将减少系统的能量。 当ΔWΔU时,系统将出现动能而导致失稳;当ΔW=ΔU时,系统处于临界状态;当ΔWΔU时,系统除去干扰后可恢复原位,系统稳定。 基于几何非线性有限元格式及虚功原理,可得干扰能量的计算式: 式中为单元弹塑性刚度矩阵;为单元几何刚度矩阵或初应力刚度矩阵。 2.3 干扰能量极值 在干扰能量的求解过程中,寻找最不利干扰位移是非常关键的。考虑到这一位移必须使干扰能量ΔΠp得到极小值,并且满足位移协调条件,选取具有待定系数的振型函数的叠加作为干扰能量试函数,然后根据极值原理求出待定系数,最终确定最不利干扰位移。 3 干扰能量法在工程稳定分析中的应用 3.1 局部和整体稳定性 采用能量的观点可得如下的稳定安全系数: 由于干扰能量是标量,可方便地进行局部和整体稳定性的分析。若以一个单元为考察对象,计入相邻单元的作用力,可得ΔUe和ΔWe,则该单元的稳定安全系数若考虑某一部分的稳定性,只需在上述基础上进行代数相加,便可得到相应的稳定安全系数。 3.2 干扰能量法 进行工程稳定分析时,可基于有限元的计算结果用干扰能量法求出每个单元的干扰能量值、安全系数值,然后分别作出相应的等值线图,则最小值曲线经过的区域即为工程薄弱部位。 4 局部破坏的影响 长江下游某岸坡高19.8 m,坡角约为50°,主要由淤泥质黏土、淤泥质亚黏土等组成,1998年2月进行观测时,发现了滑塌现象。 一般地,边坡的失稳破坏并不是瞬间发生的,而是一个由局部破坏逐渐扩展,最终贯通形成滑面的渐进过程。由于材料本身及外部因素的影响,坡体中的应力分布往往出现局部的应力集中现象,当应力超过材料允许强度时,其状态不稳定,必然导致局部破坏,而一旦发生局部破坏,必然发生应力释放、应力转移和应力重新调整,而破坏区的邻近区域所受影响最大,该部位可能由原来的没有超过强度值转变为超过强度值而发生破坏,并进行应力释放,又把多余荷载转移到其他部位,这一过程反复进行,破坏面会不断延伸。考虑到本例岸坡所处地段的河流冲刷作用较强,在发生局部破坏后,如果破坏单元和水流接触,则很容易被水流冲走,从而导致结构重新调整,岸坡的应力状态也会发生相应变化。 为此,首先采用非线性有限元对边坡进行数值分析,然后利用其计算结果,用干扰能量法得到边坡各单元的干扰能量及安全系数值,进而根据干扰能量及安全系数的等值线图确定边坡的薄弱部位,然后挖除掉薄弱部位单元,对边坡重新进行数值分析,这种过程反复进行3次,从而对边坡的失稳过程进行了数值模拟,