水工隧洞围岩支护非线性特性分析与衬砌配筋计算方法研究.pdfVIP

DoctorForum

水工隧洞围岩支护非线性特性分析与衬砌配筋

计算方法研究

孔词1,2

(1.长江设计集团有限公司，湖北武汉430010；2.武汉大学水资源工程与调度全国重点实验室，

湖北武汉430072)

为满足国家经济发展、能源以及民生需求，中国已修建了众多水力发电工程和引调水工程，并有若干类似工程

处于在建或拟建中。受限于地形地貌条件，水电和引调水工程难免包含众多复杂水工隧洞。为维持隧洞的稳定性

和安全性，围岩、初期支护以及二次衬砌结构形成联合作用系统共同承担外部荷载。当遭遇不良地质条件时，围

岩与支护结构将表现出诸多显著的非线性变形行为，合理模拟围岩与支护结构的非线性力学行为，揭示其非线性

力学特性是对水工隧洞围岩与支护结构进行分析和评价的基础，而高效计算衬砌配筋是保证衬砌运行安全的关

键。围绕水工隧洞围岩与支护结构非线性特性分析以及衬砌配筋计算方法中的几个关键问题，展开系统性研究，

主要研究内容及结论如下。

(1)形成了围岩与结构面非线性数值分析与安全评价方法。在经典弹塑性力学理论框架的基础上，考虑岩体屈

服后偏应力的损伤，建立了基于Zienkiewicz-Panda屈服准则和Drucker-Prager塑性势函数的岩体弹塑性损伤本构模

型，并给出了详细的本构方程积分算法，通过岩石圆柱体试样压缩试验，充分验证了该本构模型及其积分算法的

正确性。其次，通过严格的力学推导，建立了结构面峰值抗剪强度模型，并与薄层接触单元模型相结合，形成了

结构面非线性数值分析方法。最后，在单元屈服接近度概念的基础上，提出了单元危险系数评价指标，该指标能

对单元整个变形过程中的危险程度进行定量评价。

(2)形成了有无钢垫板全长黏结式锚杆非线性数值模拟方法，研究了锚杆非线性力学特性、破坏模式和支护性

能影响因素。引入隐式锚杆柱单元模拟锚杆的物理加固作用，在分析有无钢垫板锚杆不同支护机理的基础上，提

出了相应的支护反力计算公式，以此模拟锚杆力学约束效应。为求解锚杆支护反力，根据微元体的局部平衡条

件，建立锚杆-围岩接触面控制方程的增量格式，根据有无钢垫板锚杆的不同边界条件，给出相应的定解格式，并

考虑锚杆-围岩接触面的黏结、软化、滑移现象以及锚杆本身的屈服、硬化行为，建立了有无钢垫板锚杆的非线性

数值模型。通过若干数值算例，充分验证了该数值模型的正确性，并研究了有无钢垫板锚杆的非线性力学特性和

破坏模式。以圆形隧洞为例，详细研究了钢垫板、锚杆安装时机、锚杆长度以及锚杆半径对锚杆受力特性、支护

效果的影响。

(3)提出了一个新的钢拱架结构单元模型，形成了钢拱架+喷射混凝土联合支护结构非线性数值模拟方法。采用刚

度等效原理，对钢拱架结构的复杂截面型式进行适当简化，建立了钢拱架结构的数学模型，并基于钢材二次塑性流

动模型和金属材料混合硬化模型，建立了钢拱架结构的混合硬化弹塑性本构模型。基于有限元的变形等效、节点荷

载等效以及单元平衡方程，推导钢拱架结构的刚度叠加公式，以此模拟钢拱架结构的力学作用，形成了钢拱架结构

的隐式非线性有限元模拟方法，并通过圆形隧洞算例进行验证。然后，基于应力张量的谱分解原理和弹性损伤力学

理论，建立喷射混凝土弹性损伤本构模型，最终形成了钢拱架+喷射混凝土联合支护结构非线性数值分析方法。最

后，将所提方法应用于圆形和马蹄形隧洞中，系统研究了联合支护结构的非线性承载特性。

(4)建立了混凝土塑性损伤本构模型，给出了该模型参数的校准方法和流程，并结合钢筋组合式模型模拟方

法，形成了钢筋混凝土衬砌非线性数值分析方法。基于有效应力空间塑性理论、应力张量的谱分解原理和基于热

力学框架的连续损伤力学理论，建立了混凝土塑性损伤本构模型，并给出了详细的本构方程积分方法，通过若干

数值实验充分验证了该模型的合理性。通过单轴循环拉伸/压缩试验，对影响塑性损伤本构模型的核心参数进行了

详尽的参数敏感性分析，给出了模型参数的校准方法和流程。采用组合式模型模拟钢筋结构的力学作用，并结合

混凝土的塑性损伤本构模型，形成钢筋混凝土衬砌非线性数值分析方法。考虑混凝土损伤以后的内力重分布，以

内力平衡方程为基础，推导了钢筋应力和裂缝宽度计算公式。

(5)提出了满足承载能力极限状态计算规定以及裂缝宽度限制要求的钢筋混凝土衬砌配筋计算方法，并开发了

相