248土质边坡穗定分析一原理,方法.程序 [Bw] asato N 0 N2 0 N3 N4 2x+2x,1x+g m0 m2 0 3 0 m4 0m210m20 Im21 I11 722 /712 23 713 m24 /714 其中 mi=guNis+21 (1=1,2,3,4) (9.69) m2i=o2INis+222 Nit (2)单元矩阵的积分。将式(955)、式(956)、式(957)、式(9.58)分别代入式(925) 式(9.33),可以得到相应各单元矩阵。注意到dxdy=det[/dtds,可得到一般的表达式 ∫1F(s)d。用高斯积分法来计算式(925)式(93)各单元矩阵系数的值。一般用四点 法,此时 F(s,1) a, a F(S,,t,) (9.71) 式中:s1=1=0.57735;52=12=-0.57735;a1=a2=10。 924本构关系 进行固结和结构分析的一个重要工作是确定土的应力应变关系,即矩阵[O。这里,涉 及到土的本构关系。已有许多这方面的研究成果,现结合土体的应力应变分析,作一简要的 1.弹性模型 建立在广义定律基础上的弹性理论对式(6.32)中的[C]=[C]的表达式为 E (1+v)(1-2v) 式中:E为弹性模量;为泊桑比。 2.非线性弹性模型 土的变形特征至少与周围应力a有直接的关系,因此,相应不同的a采用不同的E,v值 是比较符合实际的。248 土质边坡稳定分析 原理 ⋅ 方法 ⋅ 程序                                 ∂ ∂ + ∂ ∂ ∂ ∂ + ∂ ∂ ∂ ∂ + ∂ ∂ ∂ ∂ + ∂ ∂ = 21 11 22 12 23 13 24 14 21 22 23 24 11 12 13 14 1 2 3 4 1 2 3 4 21 22 11 12 21 22 11 12 0 0 0 0 0 0 0 0 = 0 0 0 0 0 0 0 0 0 0 [ ] m m m m m m m m m m m m m m m m N N N N N N N N t Q s Q t Q s Q t Q s Q t Q s Q BW (9.68) 其中 m1i = Q11Nis + Q12Nit (i =1,2,3,4) (9.69) m2i = Q21Nis + Q22Nit (i = 1,2,3,4) (9.70) (2) 单元矩阵的积分 将式(9.55) 式(9.56) 式(9.57) 式(9.58)分别代入式(9.25)~ 式(9.33) 可以得到相应各单元矩阵 注意到 dxdy = det[J]dtds 可得到一般的表达式 ∫ ∫ − − 1 1 1 1 F(s,t)dsdt 用高斯积分法来计算式(9.25)~式(9.33)各单元矩阵系数的值 一般用四点 法 此时 (9.71) ∫ ∫ − − ∑∑ = = ≈ 1 1 1 1 2 1 2 1 ( , ) ( , ) i j i j i j F s t dsdt α α F s t 式中 s1 = t1 = 0.57735 ; s2 = t2 = −0.57735 ;α1 =α 2 =1.0 9. 2. 4 本构关系 进行固结和结构分析的一个重要工作是确定土的应力应变关系 即矩阵[C] 这里 涉 及到土的本构关系 已有许多这方面的研究成果 现结合土体的应力应变分析 作一简要的 回顾 1. 弹性模型 建立在广义定律基础上的弹性理论对式(6.32)中的[C]= [Ce ]的表达式为             − − − + − = ν ν ν ν ν ν ν 2 1 0 0 1 0 1 0 (1 )(1 2 ) [ ] E Ce (9.72) 式中 E 为弹性模量 ν为泊桑比 2. 非线性弹性模型 土的变形特征至少与周围应力σ3有直接的关系 因此 相应不同的σ3采用不同的 E ν值 是比较符合实际的