flac边坡计算方法.docxVIP

AC法简述 各种数值计算方法在岩土、边坡稳定分析中应用虽近几年有很大的发展，但 这些方法其理论本身以及采用的算法都有着各自的局限性。 有限元法和边界元法 都有小变形的假设。近几年发展起来的快速拉格朗日法（Fast Lagrangian Analysis of Continua，简称FLAC法）则是较好的吸取了其他数值方法的优点并克服其缺点而形 成的一种新型数值计算方法。 FLAC法基本原理类同于离散单元法，但它却能像有限元那样适用于多种材料模式 与边界条件的非规则区域的连续问题求解。在求解过程中，FLAC又采用了离散元的 动态松弛法，不需求解大型联立方程组，便于计算。另外，FLAC法不但能处理一 般的大变形问题，而且能模拟岩体沿某一软弱面产生滑动的变形。 FLAC法能针对 不同的材料特性，使用相应的本构方程来比较真实地反映实际材料的动态行为。FLAC 法还可考虑锚杆、挡墙、抗滑桩等支护结构与围岩的相互作用。 FLAC法基本原理 下面以显式拉格朗日有限差分为例，介绍FLAC法的基本原理。 对于函数F，由高斯定理， OF - J Fn dS = J —— dV （2.130） 式中:V——函数求解域（或单元）的体积； B——V的边界； ni——V的单位外法线矢量。 dF 定义梯度丁的平均值为: o xi （OF） = :JvOFdV（2.131） ox V V ox ii 式中，＜.＞表示求平均值。 对于一个具有N条边的多边形，上式可写成对N条边求和的形式。 ,OF、 1“—… （云）=V￡ F n AS . （2.132） i 式中：^Sj—多边形的边长； F ——F在△、「上的平均值。 ou浇-（2.133） oxj假定以速度气代替式（2.132）中的％，且气取各条边两端的结点（即差分网络的角 点） ou 浇-（2.133） ox j （°t） = _L￡「（u a + Ua）n AS - ox 2V 「i i j j」 j 对于四边形单元（2.16）： （°^_ 1（ A ） ox 2Vj（2.134）U ⑴ + U（2）n AS（ a）+ （°^_ 1 （ A ） ox 2V j （2.134） i i j i i i j i +U（3）+ U（4）n AS（c）+ U（4）+ U（i）n AS（d） i i j i i i j i 同理可求出当］值。则由几何方程可求得单元的平均应变增 a」 同理可求出 量: 量: ,a 、 11,au，au、| a 、 (Ae ) = I ( Q + ( i) I At (2.135) ij 21 ox ex I j i 由广义虎克定律，各向同性材料的本构方程 (2.136)为 (2.136) aij = 2vejj + 人第弓 其中，入、v为拉梅常数:。=￡订=￡ 11+￡ 22+S 33为体积应变。 t l 3=J) 0 (一 j) 因此，单元的平均应力增量可表达 成： (5)= 3%(谯)+ (5)= 3%(谯)+ 2/以气) (2.137) 同时，若以应力表示应变，则其本构关系 为 (Ae ) = 1—^(Ab ) + — 18 (2.138) ij E ij E i ij 式中：p——泊松比； E——弹性模量； I1——应力第一不变量。 这样，通过上述各式的迭代求解，便可求出每一迭代时步相应各单元的应力和应变值。 上面已求出了各域(单元)的应力，下面来求各结点的平衡力。 由结点的运动方程： * + Pg. = Pii.(2.139) i 式中：气——总加速度； g. 重力加速度。 对式(2.139)沿积分路径积分(图2.17)得 Pii = — ￡(。)n AS +pg (2.140) i V . j i i  其中，￡ ( a ..)n^S.为某结点周围单元作用在该结点上的集中力。 .. 1「 u = — F + g (2.141) 式中：F——指作用在结点上的合力(净 力)。 利用中心差分，得某结点加速度和速 度： ..u (t + At / 2) 一 u (t — At / 2)叩)=―H―-(2.142) u (t + At/2) = u (t —At/2) + u (t)At(2.143) i i i 其中，如u,(t — At/2)为结点上一时步的速度，而％［《］△ t也已求出。 进一步得结点位移： u (t + At) = u (t) + u (t + At/2)At (2.144) i i i §2.8.3 FLAG法计算过程 按照上述思路，通过迭代求解，便可求出各个时步边坡上各单元(或结点)的应力、位 移值，进而可模拟出整个边坡变形破坏的全过程。其迭代求解过程如图2.18所示。 结点不平衡力结点速度单元应力单元应变结点位移应力与应变关系 结点不平