system: [SYSTEM]

user: 3. 人字形支护结构平面有限元模型 3.1. 基本假定 1) 忽略排桩与土的竖向摩擦效应及空间效应，假设人字形支护结构为平面杆件体系； 2) 假设排桩桩底及斜向加筋水泥土地下连续墙墙底为铰支或固支；斜向加筋水泥土地下连续墙加筋体与桩顶的连接为铰接或刚接； 3) 将坑底以下排桩及斜向加筋水泥土地下连续墙的嵌固段视为弹性地基梁(如图2所示)。 图2. 人字形支护结构平面杆系有限元力学模型 3.2. 土压力确定 1) 墙后土体作用在墙上的土压力 根据朗金土压力理论可得作用在坑底以上假想墙背BC上(如图1所示)的土压力强度及总的土压力分别为： σ a = γ z K a − 2 c K a (1) E a ， = 1 2 ( γ H 1 K a − 2 c K a ) H 1 (2) 式中： K a = tan 2 ( 45 − φ / 2 ) 为主动压力系数；z为计算点深度(m)； c ， φ 分别为坑底以上墙后土的黏聚力(kPa)和内摩擦角(˚)， γ 为墙后土的重度(kN/m 3 )； H 1 为基坑开挖深度(m)。 为求得坑底以上作用在斜向加筋水泥土地下连续墙的土压力，偏于安全地假设斜向加筋水泥土地下连续墙与土的接触面光滑，则土压力计算简图如图1所示。取墙后土体ABC为研究对象，由平衡条件可得坑底以上作用在斜向加筋水泥土地下连续墙上总的土压力为： E a = 1 2 sin α ( γ H 1 K a − 2 c K a ) H 1 (3) 式中： α 为斜向加筋水泥土地下连续墙与水平面的夹角(˚)。 随着基坑开挖深度的增加，人字形支护结构前后不平衡土压力增加，人字形支护结构及桩墙之间的土体将类似于俯斜式挡土墙向基坑开挖侧移动。结合相关土压力分布规律研究，作用在人字形支护结构的整体土压力可按以下分布形式确定，坑底以上的土压力 σ a 1 + σ b 1 呈三角形分布，最大值在坑底，坑底至墙底土压力 σ b 2 呈矩形分布。由(3)可推得坑底以上及以下作用在斜向加筋水泥土地下连续墙的法向土压力强度分别为： σ a 1 + σ b 1 = γ z K a − 2 c K a (4) σ b 2 = γ H 1 K a − 2 c K a (5) 式中： σ a 1 为墙后土压力通过桩墙土体传递于坑底以上排桩上的土压力(kPa)； σ b 1 为作用于坑底以上斜向加筋水泥土地下连续墙上的法向土压力(kPa)； σ b 2 为作用于坑底以下斜向加筋水泥土地下连续墙上的法向土压力(kPa)，如图2所示。 基坑开挖后，在斜向加筋水泥土地下连续墙上的土压力作用下，由于坑底以下桩墙之间距离较大，水泥土墙后土压力通过桩墙之间土体的传递作用弱，因此，坑底以下不考虑水泥土墙后土压力传递给排桩。由于坑底以上桩墙之间的距离较近，桩墙间的土体会将水泥土墙后土压力传递给排桩，对排桩而言是倾覆力矩，对斜向加筋水泥土地下连续墙而言是抗侧移土压力。假设排桩上的土压力 σ a 1 及斜向加筋水泥土地下连续墙上的土压力 σ b 1 按下式对整体主动土压力( γ z K a − 2 c K a )进行分配 [ 8 ] σ a 1 = K 1 ( γ z K a − 2 c K a ) σ b 1 = K 2 ( γ z K a − 2 c K a ) K 1 + K 2 = 1 } (6) 式中：K 1 、K 2 分别为坑底以上考虑土压力传递效应的排桩桩背及斜向加筋水泥土地下连续墙墙背土压力分配系数。 分配系数K 1 按如下方法确定 [ 8 ]。见图3，桩墙间土体以排桩E点为滑动起点，坑底面以上桩墙间土体在重力作用下对滑动点E产生的抗滑力矩为M E ；由于受后排桩阻挡，墙后滑动面内土体以C点为滑动起点，墙后滑动面内土体在重力作用下对滑动点C产生的抗滑力矩为M C 。M E 与M C 一起组成土体的总抗滑动力矩M。分配系数K 1 定义为抗滑力矩M E 与总抗滑力矩M之比值，整理得分配系数K 1 为： K 1 = M E M = b 2 H 1 2 tan 2 ϕ D ϕ D = 2 { 45 ∘ − arc [ tan ( 45 ∘ − ϕ 2 ) ] − 2 c γ H 1 } } (7) 式中： ϕ D 为水泥土墙后土体等效内摩擦角(˚)，对于黏性土，按式(7)中第2式计算。 图3. 分配系数计算简图 2) 桩墙之间土体作用在排桩上的土压力 基于排桩和斜向加筋水泥土地下连续墙对桩墙之间土体的约束，近似将桩墙之间的土体视为侧限条件下无限长弹性土体，则作用在排桩上的土压力为 [ 9 ] σ = μ γ 1 z / ( 1 − μ ) 。式中： μ 、 γ 1 分别为桩墙之间土的泊松比及重度(kN/m 3 )；z为计算点深度(m)。同时，假设桩墙之间的土体作用于基坑开挖面以上排桩上的土压力为三角形分布，即坑底处的土压力强度为 p 1 = μ γ 1 H 1 / ( 1 − μ ) ；桩墙之间的土体作用在坑底以下排桩上的土压力为矩形分布，其土压力强度为 p 1 = μ γ 1 H 1 / ( 1 − μ ) (如图2所示)。 3.3. 斜向水泥土墙与土的侧摩擦阻力确定 斜向加筋水泥土地下连续墙侧摩阻力采用桩土界面传递函数法来考虑，将斜向水泥土墙划分成若干单元，每一单元与土体之间用非线性弹簧联系，以模拟斜向水泥土墙与土之间的摩擦特性。这些非线性弹簧的应力-应变关系就是墙侧摩阻力 τ 与剪切位移s间的关系。本模型的传递函数采用Kezdi形式，其表达式为 [ 10 ] τ ( z ) = K γ z tan φ [ 1 − exp ( − k s s u − s ) ] (8) 式中： K = 1 − sin φ ′ 为土的侧压力系数； 分别为土的有效内摩擦角(˚)及土的重度(kN/m 3 )；k为与土的类别及密实度有关的系数。s u 为墙侧摩擦阻力充分发挥时的临界位移，根据墙侧土的情况，可取3~6 mm。 3.4. 人字形支护结构平面杆系有限元力学模型及计算实现 根据前述的分析可得作用在坑底以上及以下斜向加筋水泥土地下连续墙上的法向土压力分别为 σ b 1 和 σ b 2 ；作用在坑底以上及以下排桩上的土压力分别为 σ a 1 + σ 及 。基于前述的假定可得人字形支护结构平面杆系有限元力学模型如图2所示。 1) 排桩及水泥土墙的刚度矩阵 建立刚度矩阵时，不考虑排桩桩侧摩擦作用，只考虑斜向水泥土墙的侧摩擦阻力。将排桩及斜向水泥土墙划分为m个单元，基于结构力学中的矩阵位移法，可获得排桩及斜向水泥土墙的单元刚度矩阵分别为： k p e = [ 12 E I / l 3 6 E I / l 2 − 12 E I / l 3 6 E I / l 2 6 E I / l 2 4 E I / l − 6 E I / l 2 2 E I / l − 12 E I / l 3 − 6 E I / l 2 12 E I / l 3 − 6 E I / l 2 6 E I / l 2 2 E I / l − 6 E I / l 2 4 E I / l ] (9) k w e = [ E A / l 0 0 − E A / l 0 0 0 12 E I / l 3 6 E I / l 2 0 − 12 E I / l 3 6 E I / l 2 0 6 E I / l 2 4 E I / l 0 − 6 E I / l 2 2 E I / l − E A / l 0 0 E A / l 0 0 0 − 12 E I / l 3 − 6 E I / l 2 0 12 E I / l 3 − 6 E I / l 2 0 6 E I / l 2 2 E I / l 0 − 6 E I / l 2 4 E I / l ] (10) 式中：E为桩身或墙身材料弹性模量(MPa)；I为桩或墙的截面惯性矩(m 4 )；l为单元长度(m)，A为水泥土墙单位长度的面积(m 2 /m)。 由结构力学可知，单元节点力与节点位移之间的关系为： { f } i = [ k b e ] i { w } i (11) 式中：对于排桩： { f } i = { V j − 1 M j − 1 V j M j } i T 为第i单元的节点力列阵； { w } i = { w j − 1 θ j − 1 w j θ j } i T 为第i单元节点位移列阵； 对于水泥土墙： { f } i = { U j − 1 V j − 1 M j − 1 U j V j M j } i T 为第i单元的节点力列阵； { w } i = { u j − 1 w j − 1 θ j − 1 u j w j θ j } i T 为第i单元节点位移列阵。 [ k b e ] 对应排桩或水泥土墙的单元刚度矩阵。 对坑底以下排桩及斜向加筋水泥土地下连续墙的嵌固段视为弹性地基梁，将排桩或水泥土墙与地基的接触面也相应划分为m个子域，各子域的地基反力可视为均匀分布，则每个单元上的地基反力的合力为： R j = K j w j (12) 式中： K j = l b 为排桩或水泥土墙与地基接触面积上的集中变基床系数；b为排桩或墙的计算宽度。 利用单元节点的位移连续性条件及节点力平衡条件，由式(11)可得排桩或水泥土墙的整体刚度方程为： { F } = [ K B ] { w } (13) 式中：对于排桩： { F } = { F 1 F M 1 ⋯ F m F M m } T ， { w } i = { w 1 θ 1 ⋯ w m θ m } T ； 对于水泥土墙： { F } = { F N 1 F 1 F M 1 ⋯ F N m F m F M m } T ， { w } = { u 1 w 1 θ 1 ⋯ u m w m θ m } T ； [ K B ] 对应排桩或水泥土墙的整体刚度矩阵。 2) 地基刚度方程 在土压力作用下排桩或墙节点有位移和转角，根据接触条件，地基中相应的节点也应考虑沉降和转角。但弹簧与排桩或墙接触处只能承担反力 R j = K j s j ，而不能抵抗转动，因此，基底反力偶 R m j = 0 。由此可得地基的整体刚度方程为： { R } = [ K S ] { s } (14) 式中： [ K s ] 对应于排桩或水泥墙的地基刚度矩阵。 对于排桩的地基：其节点力和节点位移列阵分别为： { R } = { R 1 ， R M 1 ， ⋯ ， R m ， R M m } T ，     { s } = { s 1 ， θ 1 ， ⋯ ， s m ， θ m } T ； 对水泥土墙的地基：其节点力和节点位移列阵分别为： { R } = { N 1 ， R 1 ， R M 1 ， ⋯ ， N m ， R m ， R M m } T ，     { s } = { u 1 ， s 1 ， θ 1 ， ⋯ ， u m ， s m ， θ m } T 3) 排桩或斜向加筋水泥土地下连续墙与地基相互作用的整体刚度方程 对于坑底以下视为弹性地基梁的排桩及斜向加筋水泥土地下连续墙的嵌固段，根据各节点的位移连续性条件及静力平衡条件，由式(13) (14)可得排桩或斜向加筋水泥土地下连续墙和地基相互作用的整体刚度方程为： { P } = [ A ] { w } (15) 式中： [ A ] 为排桩或斜向加筋水泥土地下连续墙的嵌固段和地基相互作用的整体刚度矩阵； { P } 为嵌固段的节点力列阵； 为嵌固段的节点位移列阵。 4) 人字形支护结构整体刚度方程 利用坑底以上与坑底以下的排桩和斜向加筋水泥土地下连续墙连接节点、排桩桩顶和斜向加筋水泥土地下连续墙墙顶连接节点的位移连续性条件及静力平衡条件，由式(13) (15)可得排桩与斜向加筋水泥土地下连续墙组成的人字形支护结构整体刚度方程为： { F } = [ K ] { δ } (16) 式中： { K } 为人字形支护结构整体刚度矩阵。由式(16)可看出，如果已知人字形支护结构的节点力，就可求出相应的节点位移。 问题的边界条件为： 桩 底 为 铰 支 时 : M n = 0 ， w n = 0 桩 底 为 固 支 时 : w n = 0 ， φ n = 0 墙 底 为 铰 支 时 : M n = 0 ， u n = w n = 0 桩 底 为 固 支 时 : u n = w n = 0 ， φ n = 0 } 。 (17)