论文部分内容阅读
【摘要】通过对涟钢冷轧主厂房边坡的地质情况调查研究，用ANSYS有限元软件对边坡的稳定性情况进行了模拟分析，根据分析结果和边坡具体情况，设计采用预应力锚索加网格梁加固方案。然后对其进行边坡加固，并对加固效果和锚索的受力情况进行了模拟分析，分析结果证明加固效果良好。
【关键词】锚杆（索）；ANSYS有限元；锚固系统；影响因素
LIANGANG cold rolling plant main slope stability analysis and evaluation of the effect of anchoring systems
Jiang Zhi-xin，Lv Ping-bo
（Nuclear Industry Geological Bureau of Hunan Province， 303 Brigade Changsha Hunan 410000）
【Abstract】LIANGANG rolled through the main plant Slope geology research， using ANSYS finite element software for slope stability conditions were simulated and analyzed， according to the results of the analysis and the specific circumstances of the slope， designing prestressed anchor Mississauga net grid beam reinforcement program. Then its slope reinforcement， and reinforcement effect and anchor the force carried out a simulation analysis and the results prove that the reinforcement effect is good.
【Key words】Bolt （Cable）；ANSYS finite element；Anchoring system；Factors
1. 引言
涟钢冷轧主厂房为湖南省重点建设项目，该主厂房工程建筑面积约90000m2 ，由12 个布置合理的纵横跨间组成，均采用全钢刚架结构，最大跨度 36m，基本柱距 12m，最大柱距 36 m，厂房高度 20m～42m；柱子采用阶形柱，上柱为焊接“H”型钢，下柱为钢管砼格构柱；吊车梁为实腹式焊接“H”型钢梁， 高度1200mm～3000mm，最高轨面标高17.5m。厂房沿线地形起伏较大，坡角陡，一般50°左右，出露以泥灰岩为主，岩体风化程度普便较大。所以，边坡的稳定性对厂房的安全至关重要。因此，研究此类边坡的稳定性，对于设计和施工有非常重要的现实意义[1]。
2. 风化岩体边坡的特点
（1） 风化残积层。
该层厚度不大，大多数地段基岩裸露，出露覆盖层主要为沉积物，岩性上部主要是 （含砾）低液限粘（粉）土，褐黄色、土黄色、灰黄色等，可塑～硬塑状，砾质成份主要为硅化岩角砾，粒径0.5～2cm，次棱角～次圆状，含量约占25～35%左右，层厚0～5.85m，分布于山坡坡脚地段。
（2）泥灰岩。
青灰色、灰色，泥中厚层薄层状，岩质较软，易风化。
全风化层：岩性已全风化成碎石土组成，层厚0～2.5m。
强风化层：节理裂隙较发育，岩芯破碎，呈碎块状，层厚0～2.5m。
中风化层：节理裂隙发育，岩石相对较硬，钻孔岩芯较完整，多呈短柱状。
经原位测试、室内试验、并参考有关规范规程结合当地工程经验，推荐各工程地质分区的主要岩土力学指标值如表1所示[1]。
3. 风化岩体对边坡稳定性的影响
3.1 由于风化岩体结构的破坏，从很多方面降低了边坡的稳定性，具体可从以下几个方面来分析[2～3]。
（1） 透水性能增加，抗剪强度降低，从而降低了边坡的稳定性。
在泥灰岩中，由于各矿物颗粒的膨胀率各异。当风化作用导致其不均匀的膨胀和收缩；或者由于水的溶蚀、水解岩中某些可溶物质，都将削中和破坏岩石中矿物颗粒之间的联结，导致岩石力学性能降低，特别是岩石的透水性能增大，抗剪强度降低。
（2）形成和加剧的岩体裂隙，导致软弱结构面的产生，从而影响边坡的稳定性。
由于风化作用，使岩体沿着已有的联结比较软弱的地方，如未裂开的层理、片理、劈理方向上，矿物结晶颗粒之间的结合面上，以及在矿物结晶颗粒解理方向上，形成新的裂隙。或对原有裂隙进一步加深、增宽、延伸和扩大。这种形成或加剧岩体裂隙的作用，主要是由于水的楔入作用和冻胀作用。
3.2 在边坡开挖后，通过对现场工程地质情况的调查，发现在强风化带、中风化带中岩体裂隙较为发育，裂隙面倾角较陡，裂面多呈闭合状，且裂隙中多有泥质胶结的充填。幸运的是没有一组节理倾向与边坡的临空面倾向一致[4]。
4. 边坡开挖有限元计算
（1）本文利用ANSYS有限元软件对开挖后的边坡的坡面进行稳定性验算，通过验算结果了解开挖后边坡的稳定性情况。图1为边坡开挖前和开挖后的整体模型和网格划分。
（2）模型采用D-P屈服准则。建立模型时，先生成加固锚索所在位置的关键点和节点，然后再根据节点建立锚索单元，采用beam3二维梁单元来模拟锚索，这样可以得到锚索所受到的轴力和剪力，整个支护结构共有12根锚索。然后再建立台阶处的关键点，生成台阶处的网格梁单元，根据加固设计的情况整个模型共建立了5个台阶，因为网格梁可以看作是柔性支护结构，所以也可以采用beam3二维梁单元来进行模拟。接下来，再建立边坡的坡面曲线和坡体的面，在此过程中生成开挖时需要杀死的单元的面，这样就可以方便的通过生死单元来模拟边坡的逐级开挖的过程了[5]。 （3）计算模型经过分台阶开挖荷载逐步释放的计算后，可以看到在未进行任何有效的支护的情况下边坡开挖后的位移和受力情况，图2为边坡在分台阶开挖完成之后的竖向位移图。从图中可以看出，在开挖之后边坡有明显的向下滑动趋势，边坡的竖向位移最大值达到7.13cm。从图中还可以清楚地看出，如果步对开挖的边坡进行相应地支护措施，则由于边坡的滑移还会引起了左侧向上隆起，隆起量可以达到将近20cm。
（5）从图4边坡的总的位移矢量图中我们也能很清楚的看出，边坡在没有有效的支护的情况下，即便是分多台阶进行开挖也会对整个边坡造成很大的扰动，使边坡产生明显的滑移。另外还可以从位移的云图中看出边坡单元的运动趋势呈圆弧状破坏。
5. 边坡加固设计与稳定性分析
（1）坡面每10m高设一台阶，台阶宽2m，上设截水沟。在边坡稳定性系数不足1.2的断面之间以3.5～5m间距，沿坡面布置规格400×400mm的地梁， 每根地梁安装3根预应力锚索。第一级坡面的地梁之间浆砌片石窗式护坡，窗内浆砌片石造景。从第二级至第三级坡面的地梁之间浆砌片石窗式护坡，岩层内安装泄水管，窗内坡面栽植灌木和植草。第四级至第五级坡面强风化层较厚，为了防止地表水渗入岩层，造成岩层软化，采取先在岩层内注浆封闭透水层，再在坡面栽植灌木和植草。其布置如表2所示[6]。
（3）如图5边坡加固后的竖向位移云图所示，在对边坡进行的相应的加固措施之后，边坡的竖向位移有显著的减小，边坡的竖向滑移最大值减小到了2.92cm，左侧向上隆起位移最大值也减小到了8.57cm。对边坡的水平位移而言，加固后与加固前相比，也可以清楚的看出，加固后的边坡水平向位移也有很明显的减小。边坡水平位移的最大值由加固前的21.7cm减小到了6.90cm，而且边坡左侧的水平位移最大值也由边坡加固前的10.5cm减小到了加固后的4.35cm。由这两个数据的对比我们可以清楚的看出，边坡的滑动在加固后得到了明显的控制（边坡加固后的水平向位移云图见图6）。
（4）图7为加固后的边坡的总的位移矢量图，从图中我们也能很清楚的看出，边坡在进行了有效的支护的情况下，整个边坡的变形得到了很好的控制。综合边坡的竖向位移和水平位移以及总的位移，可以说这次加固就本边坡而言还是比较成功的。
（5）对于锚索（杆）的具体受力情况，我们同样在进行每步的开挖运算之后，给出锚索在工作状态时的轴力和剪力的分布图，如图8a和图8b所示。对于网格梁的具体受力情况，我们同样在进行每步的开挖运算之后，给出网格梁在工作状态时的轴力和剪力的分布图，在本文的分析中是选用ANSYS的二维梁单元beam3来对网格梁进行模拟分析的，网格梁的受力情况如图8c所示。
6. 结论
本章通过对边坡的地质情况调查，参考相关的资料，对其边坡的情况进行了描述，同时选取了相关的一些参数。然后通过ANSYS有限元程序对其边坡的稳定性情况进行了模拟分析，分析边坡在开挖后的位移情况，说明边坡处于不稳定状态。根据边坡稳定性的分析结果，确定了预应力锚索加地格梁的加固方案，并通过有限元方法验证了此加固方案是合理和有效的。随后又对影响边坡锚固效应和锚固支护效果的因素进行了分析研究，并对支护系统的受力进行了分析。提出理论一些对工程有指导意义的方法和措施，但是岩土锚固问题对于不同的工程还需要因地制宜进行合理有效的分析和评价。
（6）从上图的边坡支护系统的受力情况可以清楚地看出，锚杆对边坡的滑动起到了很好的抑制作用，锚索（杆）的应力最大值为0.55MPa，位于第9根锚杆的后端位置。从锚索（杆）的受力特点看，锚索（杆）的受力集中在锚杆的中部，呈现中部大，两端小的情况。另外，在下部的几根锚索（杆）的里端，出现了受拉的情况，受拉的拉应力值不大。对整个边坡的稳定性没有影响。从图中还可以看出，锚索的轴力和剪力沿锚索全长呈非线性分布，但由于它随着岩体变形而变化，所以要给出它的轴力和剪力的简单计算公式或实测值是很困难的[9]。文献[10]中也给出了锚杆在工作时和拉拔试验时受力的分布情况，与本次模拟的结果比较一致，也说明了本边坡的分析是比较符合实际的。
参考文献
[1] 汤雷. 锚固力作用机理及锚固技术研究. [中国矿业大学博士论文]，2010.
[2] 张发明，邵蔚侠. 岩质高边坡预应力锚固问题研究. 河海大学学报，2009年 第27卷 第6期 73～76.
[3] 陈卫忠，朱维申. 节理岩体加固效果及其在边坡工程中的应用. 勘察科学技术，2011. 1： 3～6.
[4] 郝文化. ANSYS土木工程应用实例[M]. 北京： 中国水利水电出版社， 2009.
[5] 交通部第二公路勘察设计院. 路基. 北京：人民交通出版社，2006.5： 40～60.
[6] 张宏博，李英勇，宋修广. 边坡锚固工程中锚索预应力的变化研究. 山东工学学报（工学版），2002.12： 575～578.
[7] 冯文学. 预应力锚索在路基边坡加固施工中的应用. 山西建筑， 2003.5： 233～234.
[8] 孙玉科，牟会宠，姚宝魁. 边坡岩体稳定性分析. 科学出版社，2008： 95～120.
[9] 陆士良，汤雷，杨新安. 锚杆锚固力与锚固技术. 煤炭工业出版社，2008.11：138～1.
【关键词】锚杆（索）；ANSYS有限元；锚固系统；影响因素
LIANGANG cold rolling plant main slope stability analysis and evaluation of the effect of anchoring systems
Jiang Zhi-xin，Lv Ping-bo
（Nuclear Industry Geological Bureau of Hunan Province， 303 Brigade Changsha Hunan 410000）
【Abstract】LIANGANG rolled through the main plant Slope geology research， using ANSYS finite element software for slope stability conditions were simulated and analyzed， according to the results of the analysis and the specific circumstances of the slope， designing prestressed anchor Mississauga net grid beam reinforcement program. Then its slope reinforcement， and reinforcement effect and anchor the force carried out a simulation analysis and the results prove that the reinforcement effect is good.
【Key words】Bolt （Cable）；ANSYS finite element；Anchoring system；Factors
1. 引言
涟钢冷轧主厂房为湖南省重点建设项目，该主厂房工程建筑面积约90000m2 ，由12 个布置合理的纵横跨间组成，均采用全钢刚架结构，最大跨度 36m，基本柱距 12m，最大柱距 36 m，厂房高度 20m～42m；柱子采用阶形柱，上柱为焊接“H”型钢，下柱为钢管砼格构柱；吊车梁为实腹式焊接“H”型钢梁， 高度1200mm～3000mm，最高轨面标高17.5m。厂房沿线地形起伏较大，坡角陡，一般50°左右，出露以泥灰岩为主，岩体风化程度普便较大。所以，边坡的稳定性对厂房的安全至关重要。因此，研究此类边坡的稳定性，对于设计和施工有非常重要的现实意义[1]。
2. 风化岩体边坡的特点
（1） 风化残积层。
该层厚度不大，大多数地段基岩裸露，出露覆盖层主要为沉积物，岩性上部主要是 （含砾）低液限粘（粉）土，褐黄色、土黄色、灰黄色等，可塑～硬塑状，砾质成份主要为硅化岩角砾，粒径0.5～2cm，次棱角～次圆状，含量约占25～35%左右，层厚0～5.85m，分布于山坡坡脚地段。
（2）泥灰岩。
青灰色、灰色，泥中厚层薄层状，岩质较软，易风化。
全风化层：岩性已全风化成碎石土组成，层厚0～2.5m。
强风化层：节理裂隙较发育，岩芯破碎，呈碎块状，层厚0～2.5m。
中风化层：节理裂隙发育，岩石相对较硬，钻孔岩芯较完整，多呈短柱状。
经原位测试、室内试验、并参考有关规范规程结合当地工程经验，推荐各工程地质分区的主要岩土力学指标值如表1所示[1]。
3. 风化岩体对边坡稳定性的影响
3.1 由于风化岩体结构的破坏，从很多方面降低了边坡的稳定性，具体可从以下几个方面来分析[2～3]。
（1） 透水性能增加，抗剪强度降低，从而降低了边坡的稳定性。
在泥灰岩中，由于各矿物颗粒的膨胀率各异。当风化作用导致其不均匀的膨胀和收缩；或者由于水的溶蚀、水解岩中某些可溶物质，都将削中和破坏岩石中矿物颗粒之间的联结，导致岩石力学性能降低，特别是岩石的透水性能增大，抗剪强度降低。
（2）形成和加剧的岩体裂隙，导致软弱结构面的产生，从而影响边坡的稳定性。
由于风化作用，使岩体沿着已有的联结比较软弱的地方，如未裂开的层理、片理、劈理方向上，矿物结晶颗粒之间的结合面上，以及在矿物结晶颗粒解理方向上，形成新的裂隙。或对原有裂隙进一步加深、增宽、延伸和扩大。这种形成或加剧岩体裂隙的作用，主要是由于水的楔入作用和冻胀作用。
3.2 在边坡开挖后，通过对现场工程地质情况的调查，发现在强风化带、中风化带中岩体裂隙较为发育，裂隙面倾角较陡，裂面多呈闭合状，且裂隙中多有泥质胶结的充填。幸运的是没有一组节理倾向与边坡的临空面倾向一致[4]。
4. 边坡开挖有限元计算
（1）本文利用ANSYS有限元软件对开挖后的边坡的坡面进行稳定性验算，通过验算结果了解开挖后边坡的稳定性情况。图1为边坡开挖前和开挖后的整体模型和网格划分。
（2）模型采用D-P屈服准则。建立模型时，先生成加固锚索所在位置的关键点和节点，然后再根据节点建立锚索单元，采用beam3二维梁单元来模拟锚索，这样可以得到锚索所受到的轴力和剪力，整个支护结构共有12根锚索。然后再建立台阶处的关键点，生成台阶处的网格梁单元，根据加固设计的情况整个模型共建立了5个台阶，因为网格梁可以看作是柔性支护结构，所以也可以采用beam3二维梁单元来进行模拟。接下来，再建立边坡的坡面曲线和坡体的面，在此过程中生成开挖时需要杀死的单元的面，这样就可以方便的通过生死单元来模拟边坡的逐级开挖的过程了[5]。 （3）计算模型经过分台阶开挖荷载逐步释放的计算后，可以看到在未进行任何有效的支护的情况下边坡开挖后的位移和受力情况，图2为边坡在分台阶开挖完成之后的竖向位移图。从图中可以看出，在开挖之后边坡有明显的向下滑动趋势，边坡的竖向位移最大值达到7.13cm。从图中还可以清楚地看出，如果步对开挖的边坡进行相应地支护措施，则由于边坡的滑移还会引起了左侧向上隆起，隆起量可以达到将近20cm。
（5）从图4边坡的总的位移矢量图中我们也能很清楚的看出，边坡在没有有效的支护的情况下，即便是分多台阶进行开挖也会对整个边坡造成很大的扰动，使边坡产生明显的滑移。另外还可以从位移的云图中看出边坡单元的运动趋势呈圆弧状破坏。
5. 边坡加固设计与稳定性分析
（1）坡面每10m高设一台阶，台阶宽2m，上设截水沟。在边坡稳定性系数不足1.2的断面之间以3.5～5m间距，沿坡面布置规格400×400mm的地梁， 每根地梁安装3根预应力锚索。第一级坡面的地梁之间浆砌片石窗式护坡，窗内浆砌片石造景。从第二级至第三级坡面的地梁之间浆砌片石窗式护坡，岩层内安装泄水管，窗内坡面栽植灌木和植草。第四级至第五级坡面强风化层较厚，为了防止地表水渗入岩层，造成岩层软化，采取先在岩层内注浆封闭透水层，再在坡面栽植灌木和植草。其布置如表2所示[6]。
（3）如图5边坡加固后的竖向位移云图所示，在对边坡进行的相应的加固措施之后，边坡的竖向位移有显著的减小，边坡的竖向滑移最大值减小到了2.92cm，左侧向上隆起位移最大值也减小到了8.57cm。对边坡的水平位移而言，加固后与加固前相比，也可以清楚的看出，加固后的边坡水平向位移也有很明显的减小。边坡水平位移的最大值由加固前的21.7cm减小到了6.90cm，而且边坡左侧的水平位移最大值也由边坡加固前的10.5cm减小到了加固后的4.35cm。由这两个数据的对比我们可以清楚的看出，边坡的滑动在加固后得到了明显的控制（边坡加固后的水平向位移云图见图6）。
（4）图7为加固后的边坡的总的位移矢量图，从图中我们也能很清楚的看出，边坡在进行了有效的支护的情况下，整个边坡的变形得到了很好的控制。综合边坡的竖向位移和水平位移以及总的位移，可以说这次加固就本边坡而言还是比较成功的。
（5）对于锚索（杆）的具体受力情况，我们同样在进行每步的开挖运算之后，给出锚索在工作状态时的轴力和剪力的分布图，如图8a和图8b所示。对于网格梁的具体受力情况，我们同样在进行每步的开挖运算之后，给出网格梁在工作状态时的轴力和剪力的分布图，在本文的分析中是选用ANSYS的二维梁单元beam3来对网格梁进行模拟分析的，网格梁的受力情况如图8c所示。
6. 结论
本章通过对边坡的地质情况调查，参考相关的资料，对其边坡的情况进行了描述，同时选取了相关的一些参数。然后通过ANSYS有限元程序对其边坡的稳定性情况进行了模拟分析，分析边坡在开挖后的位移情况，说明边坡处于不稳定状态。根据边坡稳定性的分析结果，确定了预应力锚索加地格梁的加固方案，并通过有限元方法验证了此加固方案是合理和有效的。随后又对影响边坡锚固效应和锚固支护效果的因素进行了分析研究，并对支护系统的受力进行了分析。提出理论一些对工程有指导意义的方法和措施，但是岩土锚固问题对于不同的工程还需要因地制宜进行合理有效的分析和评价。
（6）从上图的边坡支护系统的受力情况可以清楚地看出，锚杆对边坡的滑动起到了很好的抑制作用，锚索（杆）的应力最大值为0.55MPa，位于第9根锚杆的后端位置。从锚索（杆）的受力特点看，锚索（杆）的受力集中在锚杆的中部，呈现中部大，两端小的情况。另外，在下部的几根锚索（杆）的里端，出现了受拉的情况，受拉的拉应力值不大。对整个边坡的稳定性没有影响。从图中还可以看出，锚索的轴力和剪力沿锚索全长呈非线性分布，但由于它随着岩体变形而变化，所以要给出它的轴力和剪力的简单计算公式或实测值是很困难的[9]。文献[10]中也给出了锚杆在工作时和拉拔试验时受力的分布情况，与本次模拟的结果比较一致，也说明了本边坡的分析是比较符合实际的。
参考文献
[1] 汤雷. 锚固力作用机理及锚固技术研究. [中国矿业大学博士论文]，2010.
[2] 张发明，邵蔚侠. 岩质高边坡预应力锚固问题研究. 河海大学学报，2009年 第27卷 第6期 73～76.
[3] 陈卫忠，朱维申. 节理岩体加固效果及其在边坡工程中的应用. 勘察科学技术，2011. 1： 3～6.
[4] 郝文化. ANSYS土木工程应用实例[M]. 北京： 中国水利水电出版社， 2009.
[5] 交通部第二公路勘察设计院. 路基. 北京：人民交通出版社，2006.5： 40～60.
[6] 张宏博，李英勇，宋修广. 边坡锚固工程中锚索预应力的变化研究. 山东工学学报（工学版），2002.12： 575～578.
[7] 冯文学. 预应力锚索在路基边坡加固施工中的应用. 山西建筑， 2003.5： 233～234.
[8] 孙玉科，牟会宠，姚宝魁. 边坡岩体稳定性分析. 科学出版社，2008： 95～120.
[9] 陆士良，汤雷，杨新安. 锚杆锚固力与锚固技术. 煤炭工业出版社，2008.11：138～1.