第三章隧道施工过程数值模拟方法与ANSYS实现.docVIP

第三章隧道施工过程数值模拟方法与ANSYS实现 2.1隧道施工过程数值模拟方法 2.1.1开挖（卸载）的模拟 ①基本的模拟思想 隧道开挖时破坏了岩体内有的应力平衡，围岩内的各点在地应力的作用下，在一定范围内围岩产生位移，形成松弛，与此同时也会使围岩的物理力学性质恶化，也就是我们所说的“二次应力场”。隧道的开挖导致围岩应力场及位移场的变化，一般都是通过卸载过程来实现的。在对卸载过程进行模拟时，通常有采用的就是在已知边界初始地应力作用下，沿预定开挖线进行的“开挖卸载模拟方法”。这种方法的位移场真实地反应了开挖所引起的位移变化，是工程需要了解的重要部分。 ②实现卸载的具体方法 正确模拟卸载过程的效果是地下工程数值模拟的一个重要课题。开挖卸载之前，沿开挖边界上的各点都处于一定的初始应力状态，开挖使这些边界的应力解除，也就是我们所说的卸载，从而引起围岩变形和应力场的变化。对上述过程的模拟通常所采用的方法有两种：“反转应力释放法”和“地应力自动释放法”。 “反转应力释放法”是把沿开挖边界上的初始地应力反向后转换成等价的“释放荷载”，施加与开挖边界，在不考虑初始地应力的情况下进行有限元分析，将由此得到的围岩位移作为由于工程开挖卸载产生的岩体位移，由此得到的应力场与初始应力场叠加即为开挖后的应力场。对于大型的地下工程或者复杂的施工方法，应力场多次叠加，使得分析过程过于繁杂，另外，进行弹塑性分析时，由于应力场需要叠加，对围岩屈服的判断需做特殊的处理，增加了分析的复杂度，降低了分析的准确性。 “地应力自动释放法”则是认为隧道的开挖打破了开挖边界上各点的初始应力平衡状态，开挖边界上的节点受力不平衡，为获得新的力学平衡，围岩就要产生相应的变形，引起应力的重分布，从而直接得到开挖后围岩的应力场和位移场。分部开挖时，对于每一步的开挖，将这一步被开挖的单元变为“空单元”，即在开挖边界产生新的力学边界条件，然后直接进行计算就可以得到工况开挖后的结果，接着可用同样的方法进行下一步的开挖分析。 该方法更符合隧道开挖后围岩应力重分布的真实过程，反映了开挖后围岩卸载的机理，可以实现连续的开挖分析。对于弹塑性分析计算只需建立弹塑性模型，其余计算过程同线弹性，相比“反转应力释放法”有着极大的便。 2.1.2施工过程的模拟 在对隧道工程进行结构分析时，我们不但关注建成后隧道结构和围岩的稳定性，而且关注各个施工阶段中围岩和尚未完成的结构的受力和变形情况。根据新奥法施工原理，隧道开挖后，围岩从变形到破坏有一个时间历程。因此要想真实地模拟隧道开挖与支护的整个施工作业流程，不仅要考虑岩体的复杂形态，施工方法，支护施作时机等因素，当我们把岩体的变形看作是弹性或者弹塑性问题，建立平面应变模型进行隧道结构分析时，为了比较真实地模拟施工过程，常采用“应力逐步释放”的方法来模拟隧道开挖与支护的时空效应，具体实现的方法也常用“施加虚拟支撑力逐步释放法”。 “施加虚拟支撑力逐步释放法”是在“地应力自动释放法”的基础上，通过在 开挖边界施加虚拟支撑力的方法，来模拟围岩的逐步卸载，如图所示。初始应力 状态（a）为初始地应力状态，在（b）阶段，隧道的开挖引起开挖边界上的释放 节点荷载 ii f=α?f 1 。为实现这一过程，在初始应力场中开挖隧道单元的同时，在 开挖边界上各相应的节点施加虚拟支撑力(1)() 1i 1i P=?α??f，则产生新的荷载边 界条件，继续进行计算，就直接得到开挖后的围岩位移场和应力场，在阶段（c）， 这时只需要将虚拟支撑力减小为(1)() 2 i12i P=?α?α??f，继续进行计算即得到初 期支护后围岩和支护的位移和应力，在阶段（d），二次衬砌施作后，剩余节点荷 载被完全释放，只需要除去虚拟支撑力，继续计算就可以得到竣工后的围岩和衬 砌的位移和应力。其中 1 α， 2 α， 3 α为此阶段的地应力释放率，可根据量测资料加 以确定，通常近似地将它定位本阶段隧道控制测点的变形值与施工完毕稳定以后 该控制测点的总变形值的比值。在缺乏实测变形资料的情况下，亦可按工程类比 法加以确定，并根据试算结果予以修正。 i f的计算详见下节。 2.1.3释放荷载的计算 在模拟隧道施工过程中，都要进行释放荷载的计算。释放荷载的确定也有两种 方法，一种是将释放边界一侧单元的初始应力转换成相应的等效节点荷载，然后 通过叠加，计算开挖边界上各节点总的等效节点荷载 ∫ Ω Ω F e=BTd 00 σ =[]=?∑ e Te iixiy fffF 0 其中 0 σ为单元初始应力向量；BT为应变矩阵的转置；Ω为积分区间，平面问题为 单元面积，空间问题为单元体积。 另一种确定释放荷载的方法是：根据预计开挖边界量测单元的初始应力通过插 值球的节点节点上的应力，然后假定两相邻节点节点之间应力变化为线性分布