第5章 地下结构连续介质设计模型数值计算法（4）：有限元数值分析.ppt

5.3.3结构体系离散化地下结构有限单元法的基础是将地下结构和计算范围内的岩土体离散为有限个仅在节点处相互连接的单元体的组合，对于平面应变问题常用的单元包括线单元和面单元。一般情况下往往采用两节点或三节点杆单元来模拟锚杆，采用两节点或三节点梁单元来模拟喷射混凝土，采用三节点或六节点三角形常应变元、四节点或八节点四边形等参单元来模拟围岩和二次衬砌等，因此，地下结构体系离散化后往往是各种类型单元的组合。由于四边形等参单元具有应力变化连续、精度较高、便于网格划分等优点，在模拟围岩和二次衬砌时采用四边形等参单元比较合适。单元划分的疏密、大小和形状都会影响有限单元法计算精度。根据误差分折，应力的误差与单元尺寸的一次方成正比，位移的误差与单元尺寸的二次方成正比。因此，从理论上讲，单元划分得越密越小、形状越规则，计算精度则越高。但是，单元划分越密，单元数就越多，对计算设备的要求就越高，耗费的计算时间也越长。为解决这一矛盾，一般在地下结构周围区域、地质构造区域等应力、位移交化梯度较大以及荷载有突变的区域，单元划分可以加密，而其他区域则可以稀疏一些。疏密区单元大小相差不宜过大，应尽可能均匀过渡。另外，在地下结构体系离散化时尚需注意以下几点：（1）单元各边长相差不能过大，两边夹角不能过小，各夹角最好尽量相等。（2）单元边界应当划分在材料的分界面上和开挖的分界线上，一个单元不能包含两种材料。（3）集中荷载作用点、荷载突变处及锚杆的端点处必须布置节点。（4）地下结构和围岩在几何形状和材料特性等方面都具有对称性时，应充分利用其对称性取部分计算区域进行离散，以减少计算工作量。（5）单元的划分要考虑到分部开挖的分界线和部分开挖区域的分界线。5.3.4边界条件和初始应力由于地下工程都是在具有初始应力场的岩土体中开挖的，因而数值计算中一般采用施加反向释放荷载来模拟开挖过程，即由于开挖而在洞周形成释放荷载，其值等于沿开挖边界上原先的应力并以原来相反的方向作用于开挖边界上。计算范围的外边界可以采取两种方式处理：（1）位移边界条件，即一般假定边界点位移为零（也有假定为弹性支座或给定位移的，但地下工程分析中很少采用）；（2）力边界条件，由岩土体中的初始应力场确定，包括自由边界（p=0）条件。（3）定混合边界条件，即节点的一个自由度给定位移，另一个自由度给定节点力（二维问题）。图计算范围的边界条件当然，无论采用哪种边界条件都会存在一定的误差，且这种误差随计算范围的减小而增大，靠近边界处误差最大，这种现象称为“边界效应”。边界效应在动力分析中影响更为显著，需妥善处理。当地下结构为浅埋时，上部为自由边界，考虑重力作用，两侧作用三角形分布初始地应力，侧压力系数为；图浅埋隧洞的初始应力图浅埋隧洞地层分层时的初始应力当地下结构为深埋时，上部及侧向边界上均作用均布的初始地应力，可用隧道研究范围的中心点处的初始地应力作为整个研究范围的初始地应力，侧压力系数以实测为准或根据经验确定。图深埋隧洞的初始应力构造应力引起的初始地应力将不再是主应力方向，此时应根据构造应力的方向进行转化。图构造应力与坐标系XOY的关系假设构造应力的大小及方向与坐标系XOY的关系如左图所示，则根据弹性力学知识有：初始应力场是自重应力场和构造应力场的叠加。5.5地下结构有限元分析过程在选定计算范围，建立计算力学模型，进行结构体系离散后，可根据连续介质力学的识进行单元的力学特性分析和整体平衡方程的建立，在整体平衡方程中引入边界条件，即可计算单元上各节点的位移。（1）单元力学特性的分析单元力学特性分析是有限单元分析的核心部分。①形函数的选择：将每个单元内部任意一点的位移用节点位移表示。②利用几何方程建立应变与节点位移的关系；算单元应力与应变；③利用物理（本构）方程建立应力与节点位移的关系；④利用虚功原理建立节点力和节点位移的关系；（2）建立体系的平衡方程，并求解单元的应力和应变①由单元刚度矩阵组装成总刚度矩阵；②将节点荷载按照力的平衡条件汇集成体系的总荷载，得到体系平衡方向；