大型地下钢筋混凝土岔管结构优化分析

1、大型地下钢筋混凝土岔管构造优化分析摘要：本文根据岩体和混凝土衬砌结合受力特点，系统地提出了大型地下钢筋混凝土岔管构造优化分析方法。针对岔管空间几何特征，结合三维有限元的原理，提出了采用程序分析判断，自动生成各种不同类型的三维有限元地下钢筋混凝土岔管构造网格的计算方法。该方法只须输入几个关键的控制参数，便能自动剖分所须岔管网格，为大型地下钢筋混凝土岔管构造优化分析提供一种快速网格生成方法。根据岩体损伤开裂和弹塑性破坏特征，从岩体开挖应力扰动的塑性耗散能扩散理论和混凝土衬砌构造开裂损伤原理，提出了大型地下钢筋混凝土岔管构造型式优化的评估方法。通过对工程实例分析，证明该方法能综合反响地下钢筋混凝土岔

2、管构造型式和参数选择以及岔管所处的地质环境和条件的影响，为优化地下钢筋混凝土岔管构造合理型式提供了一种有效的分析方法。关键词：地下岔管三维有限元构造优化地下钢筋混凝土岔管构造是与围岩结合承载受力构造，采用常规的构造力学方法分析其衬砌型式，已不能充分考虑围岩和衬砌结合工作特性1。讨论大PD值(P表示压力管道承受的水压力，D表示压力管道的直径)地下钢筋混凝土岔管衬砌构造与围岩结合受力的优化分析方法，对合理选择地下钢筋混凝土岔管构造形式，进步岔管衬砌设计程度，缩短设计周期，加快施工进度，节省工程投资，简化施工工艺，具有重要经济和社会效应。本文根据岩体弹塑性破坏特征，采用岩体塑性耗散能评估地下钢筋混凝

3、土岔管构造开挖的合理型式；根据混凝土损伤开裂破坏特征，采用损伤应变能评估地下钢筋混凝土岔管衬砌受力的优化型式；根据三维有限元的根本原理和空间几何特征，采用交互式图形显示技术，对各种不同类型的地下钢筋混凝土岔管构造的三维有限元网格进展自动生成，为大PD值地下钢筋混凝土岔管构造优化分析提供一种快速和有效评估计算方法。1地下岔管构造优化分析有限元网格自动生成方法地下钢筋混凝土岔管构造是一种复杂的多重介质空间构造，采用有限元数值方法进展优化分析时，首先须对地下钢筋混凝土岔管构造进展离散化。复杂空间构造的三维有限元网格离散剖分工作是一项繁杂费时的工作，而且数据输入时容易出现人为过失。虽然国内有一些著名软

4、件如ViziAD，AutAD具有较好的交互式绘图功能，但是这些软件往往需要有限元使用者有一定的建模技巧，同样须花相当多的时间。当进展岔管构造优化分析时，每分析一个方案，都必须重新建模，对修改构造方案进展重新离散和重新组织数据，所以对优化参数的选择和优化构造的判断非常不便2。因此编制一个适应于各种不同类型的地下钢筋混凝土岔管有限元网格自动生成程序，对岔管进展优化分析是非常重要的。地下岔管钢筋混凝土岔管，衬砌与围岩共同承载，因此网格剖分应包括衬砌和围岩。为了适用于各种地下岔管衬砌与围岩结合承载的三维有限元网格自动生成，本文利用映射法原理编制岔管构造网格的自动生成程序。只要用户给出构造轮廓的某些特征

5、值，程序根据这些特征值，调用相关算法，即可进展网格的自动生成，不再需要用户计算单元信息(包括各种约束信息、锚杆信息、受载面信息).因此大大减小了前处理工作量，节省了时间，降低了出错率。使得地下岔管构造优化分析计算变得简捷而快速。程序生成的有限元数据可以通过图形显示程序显示出来，并可通过人机交互交式，对生成的网格进展修改、添加、删除，保证有限元分析数据的准确性。地下钢筋混凝土岔管有限元网格自动剖分程序主要分3步进展。(1)根据钢筋混凝土岔管构造的几何特征和有限单元信息规定，首先对规那么岔管构造进展离散。如图1(a)对称岔管，只需输入主、支管直径D1，D2；衬砌厚度d1,d2；主管和总管长度S1和

6、S；岔管的半岔角、半锥角和两支管的间距L，那么可分别建立主管和支管的轮廓方程：(1)式中：X，Y，Z为以岔管相贯线交点为圆点的三维坐标；0，1，2为沿X轴线所取的不同段面常数。对于岔管段可根据岔角参数(半岔角、半锥角)建立以腰梁和U梁为部分坐标的相关方程，通过数学变换那么可求所需相贯线的坐标。再根据几何关系，那么可推出岔管锥段斜截面的交线参数方程3。如图1(b)所示，任意倾角为1的斜截面与母线D的交点为P，由图中的几何关系那么可得到P点在截面上的向径P和夹角的表达式为：p=s(+1)s(-)/s(+1-)s，tg=(1-s)sins/s2+ssin2(2)式中：为母线D在对称面AB上的投影与A

7、的夹角(A)；角为母线D的旋转角，即从A开场旋转到D的角度AED的角度。根据点的坐标和精度要求控制，选取一定间距截面和不同的角，按照式(2)那么计算出斜截面上的任意交点P的坐标。对于卜型岔采用类似方法可求得岔管轮廓线尺寸和相贯线坐标。图1对称岔管构造参数图2对称岔和卜型岔衬砌单元图3对称岔整体有限元单元网格(2)对岔管外围的围岩进展网格剖分。首先可根据已离散的岔管衬砌单元和节点，按其规那么将围岩向外扩大34层单元，为此可采用对岔管衬砌外边界进展自动标识，而后将其向围岩拓展。然后根据岩体所要计算的边界范围以及网格离散的计算精度确定计算模型边界的单元编号和结点坐标。再采用插值方法对中间岩体网格进展

8、逐步加密过渡，那么可确定岔管外围的全部岩体单元网格。这一部分关键在于对岔管岔角处网格剖分，因为岔角处应力较为集中，往往是有限元分析的重点，要求网格较密，同时又要求与第一部分网格相协调。因此在对岔管外围规那么部分岩体网格进展离散时，可采用管道截面交线逐渐由圆或椭圆过渡到矩形，从而使得岔管岔角处围岩形成为一个大的三棱柱，对于较为规那么的三棱柱，可以较为容易地进展六面体网格剖分。(3)对岔管分岔角处进展修圆处理。地下钢筋混凝土衬砌岔管分岔角处的应力集中现象十清楚显，在开挖时，常常对分岔角处进展修圆处理。岔角修圆后过渡段不再是一锥管，但截面交线参数方程仍为椭圆。本文采用在分岔角处进展滚动圆滑处理，即根

9、据修圆程度要求，选用一个适宜半径的圆球，将岔管过渡段母线滚动圆滑，求得母线与斜截面的交点。然后根据修圆的控制范围选定不修圆的边界，再对修圆后的母线和不修圆的边界之间的网格节点坐标进展插值计算，那么可得到较均匀的修圆岔角网格。采用上述方法，只需修改初始的几个输入参数，那么可离散如图2、图3所示各种岔管参数网格。转贴于论文联盟.ll.2地下岔管构造形式优化评估方法地下钢筋混凝土岔管构造形式的优劣不但受岔管本身的形态的影响，还与所处的地质环境、岩体特性、原岩应力场有极大的关系。岔管构造本身的形态直接影响钢筋混凝土岔管衬砌的受力特性，而地下岔管所处的围岩稳定特性是地下岔管构造稳定的关键。因此对地下钢筋

10、混凝土岔管构造形式的优化评估，可从岔管衬砌整体受力特点和岔管洞室开挖围岩稳定特性两方面来评估。2.1岔管洞室开挖围岩稳定特性评估岩体是一种弹塑性介质，其塑性变形是不可恢复的非线性变形。当地下岔管所处的围岩特性不同、原岩应力状态相异时，地下岔管洞室开挖后，其围岩的破坏状态和应力扰动情况也就不同。岩体发生塑性破坏时，岩体的塑性应变能增量将全部耗散掉。塑性应变能耗散越多，说明地下岔管洞室开挖对围岩的扰动破坏越大，围岩的稳定性态也就越差。地下岔管洞室构造形式不同，导致的岔管洞室的应力分布和破坏特性也就不同。因此对地下岔管洞室构造形式进展优化评估时，首先可综合分析岔管洞室开挖后围岩的变形、应力分布规律和

11、围岩的张裂破坏区、剪切破坏区、压裂破坏区的大小，然后根据开挖终了围岩的塑性散耗能大小综合评估各种岔管形式的优劣。岩体的张裂破坏区是指围岩受张力而产生张性破坏的区域。当岔管洞室开挖后，单元的第一主应变1满足式(3)条件，那么认为该单元产生张裂破坏1=/kE(3)式中：为现场直剪试验推测出的岩体单轴抗压强度，E为现场岩体载重试验测得的岩体弹模，k是根据试验和设计要求所取的平安系数，为岩体的临界应变值4。岩体的剪切破坏区是指洞室开挖后围岩应力状态到达临界应力的屈从岩体区域。通常可根据有限单元计算的应力状态按所选择的屈从函数F进展判断。当单元的屈从函数F()0时，那么认为该单元岩体已产生屈从破坏。岩体

12、的压裂破坏区是指洞室开挖后，围岩的最大压应力值超过岩体抗压强度的区域。在有限元计算中，当单元的最大压应力3(为岩体的容许抗压强度)时，那么将该单元计入压裂破坏。在有限元计算迭代过程中，每迭代一次首先根据应力状态计算单元的屈从函数值F，对所有F0的单元计入塑性破坏的标识，并算出单元体积vpi。然后根据单元累积的应变状态计算单元的第一主应变1,对所有1的单元计入张裂破坏的标识，并算出单元体积vki。根据单元累积的应力状态计算单元的第三主应力3，对所有3的单元计入压裂破坏的标识，并算出单元体积vyi。在每期开挖迭代完后，按下式求出总的张裂破坏岩体vk、剪切破坏岩体Vp和压裂破坏岩体体积Vy：，岩体塑

13、性耗散能也就是岩体进入破坏后，其应力增量在对应的应变增量上所做的功。当岩体单元进入裂破坏后，该单元单位体积的塑性耗散能增量为5：=(ij-*ij)dPij=pijdpij(5)式中：ij是按加载直接计算超出于屈从面或破坏面外的计算应力状态，dpij是进入破坏后的应变增量，*ij是将ij沿垂直于屈从面法向拉到屈从面上所对应的应力，即F(*ij)=0.根据微分和矩阵原理，可求得塑性应力增量为：pij=-F(ij)/T(6)每迭代一次，对所有进入破坏状态单元的塑性耗散能进展求和，那么可求出岩体总的塑性耗散能Ep：(7)从上式可看出岩体开挖应力扰动越大，围岩的塑性功也就越大，最终导致围岩的塑性耗散能E

14、p也就越多，其围岩的稳定就越差。所以塑性耗散能Ep值综合反响了岩体分期开挖的稳定特性。2.2岔管衬砌受力状态特性评估地下岔管衬砌构造形式不同，充水运行后衬砌的受力特性也就不同。根据弹性力学原理，构造受外力作用后，其外力所做的功将全部转变成构造的应变能。有限单元的单位应变能由体变能U0v和形变能U0d两部分组成，即5：U0=1/2ijij=U0v+U0d=1/18KI21+1/2GJ2(8)式中：K=E/3(1-2),G=E/2(1+)，I1为第一应力不变量，J2为第二应力偏量不变量。衬砌构造体变能是由于体积变化引起的内部能量变化，不引起构造的破坏。而形变能引起构造形状变化，形变能越大越容易引起

15、构造破坏。当岩体采用非线性介质，而衬砌构造采用线弹性介质分析时，那么可用衬砌构造总的形变能大小来判别衬砌构造的受力好坏。即(9)形变能越大(I=1,2,n表示对所有衬砌单元求和)，说明衬砌构造受力条件越差。当岩体和衬砌构造都采用非线性介质分析时，那么可根据衬砌构造开裂后的损伤系数D，按开裂后衬砌构造应变能的总改变量的大小来判别衬砌构造的受力好坏。即：(10)式中：k=1,2,表示对所有开裂衬砌单元求和。D越大，衬砌单元开裂程度越大，最终累计的开裂应变能的改变量E也就越大，说明衬砌构造破坏越大。损伤系数D计算、混凝土衬砌的钢筋计算参见文献6。转贴于论文联盟.ll.3工程实例分析本文结合桐柏抽水蓄

16、能地下钢筋混凝土岔管构造型式，对地下岔管构造型式进展了三维有限元数值模拟优化分析。该电站总装机1200，利用已建成的桐柏水库作为上、下库，发电的平均水头300.地下岔管的主管直径为9.0，两支管的直径为5.5，衬砌厚度0.6.岔管处的岩体为较完好的中细粒花岗岩，岩体的物理力学参数为：弹模E=17GPa，泊松比=0.25，粘结力=1.5Pa，摩擦系数f=1.2，抗拉强度Re=1.5Pa，抗压强度R=130Pa，容重r=27kN/3.根据实测地应力值和山体地形特征，采用三维有限元进展反演分析计算，得出地下岔管处岩体最大的初始地应力值为1112Pa.表1岔管构造参数(单位：)工对称岔管参数卜型岔管参

17、数况主管直径支管直径衬砌厚度岔角/锥角/修圆系数作用水头主管直径支管直径衬砌厚度岔角/锥角/修圆系数作用水头19.05.50.640141.03759.05.50.64061.037529.05.50.653141.03759.05.50.65381.037539.05.50.660141.03759.05.50.660101.037549.05.50.653141.23759.05.50.65381.45375注：修圆系数大于1表示修圆程度。表2各种岔管构造不同计算工况下围岩和衬砌破坏指标工况岔管洞室开挖围岩破坏指标充水衬砌弹性计算充水衬砌非线性计算衬砌破坏指标岔方是塑性体积/3开裂体积/3

18、总破坏体积量/3塑性耗散能/103NT岔角最大应力/Pa岔角位移/衬砌形变能/103NT岔角最大应力/Pa岔角位移/塑性体积/3开裂体积/3压裂体积/3总破坏体积量/3塑性耗散能/103NT岔角位移/管否型修式案圆1否927186111322-0.560.391018.80.092729621234320.13对2否846172101821-0.670.40998.350.092809031183280.13称3否82216098220-0.670.40968.150.092728601132260.14岔4是917169108622-0.680.351018.20.0826990111702

19、80.121否66412979313-0.630.361018.70.0818311751358580.14卜2否70312682913-0.700.361008.30.0920811051313520.15型3否71612984513-0.700.41998.10.0924010661306500.16岔4是80414795115-0.710.341029.90.07242108011323520.15根据地下钢筋混凝土岔管构造的布置情况，对卜型岔管和对称岔管构造分别拟定了4种构造型式进展分析，每种岔管构造型式又分别计算了洞室开挖、隧洞充水运行衬砌按线弹性分析、隧洞充水运行衬砌按非线性分析3

20、种计算工况，岔管的构造参数见表1.不管是卜型岔还是对称岔，岔管构造均剖分了17616个空间等参单元，20220个节点，其中地下岔管衬砌构造剖分了4032个单元，其单元网格见图23.通过对各种岔管构造型式的分析计算，将岔管洞室和衬砌的破坏指标列入表2，从以下几个方面可以看出不同的岔管型式的受力特点和优化特性。(1)从岔管洞室构造开挖围岩稳定特性来看，对称岔的岔角越大，岔管洞室开挖后的应力扰动就越大，围岩的破坏体积也就越大；而卜型岔的岔角越大，岔管洞室开挖后的应力扰动就越小，围岩的破坏体积相应也就越校由此说明地下岔管洞室构造的优化型式不但与岔管的构造型式有关，还与岔管所处的地质环境和地质条件有关。

21、(2)比拟对称岔和卜型岔的洞室开挖稳定特性，卜型岔各种方案的岩体破坏指标比对称岔各种方案的破坏指标都要小(见图4和图5).例如方案2，卜型岔岩体总破坏体积8293比对称岔10183小1893，占22.8%，岩体的塑性耗散能前者比后者小53.8%.由此说明该工程卜型岔的洞室构造比对称岔的洞室构造受力条件要优裕一些。(3)当将岔角处的岩体进展修圆后，由于岔角处的跨度加大，不管对称岔，还是卜型岔洞室围岩的破坏指标均有所加大。说明岔管的岔角修圆对洞室围岩的整体稳定是不利的。(4)当岔管充水运行后，从衬砌的受力特点来看，不管是对称岔，还是卜型岔，岔角越大，衬砌构造的破坏指标越小，既衬砌构造受力条件越好。

22、说明衬砌构造的优化型式主要受岔管的构造参数影响，围岩地质条件影响相对较校(5)比拟对称岔和卜型岔衬砌构造的受力特性可看出，对称岔在岔角处的部分最大拉应力值比卜型岔略大，但衬砌开裂破坏的体积前者比后者平均小12.1%，塑性耗散能小86%.说明对称岔由于受力对称，衬砌总体应力扰动要小一些，但岔角处破口稍大，所以导致岔角的部分应力增大。(6)当对称岔采用1.2的修圆系数修圆岔角，岔角处的部分应力由8.35Pa降到8.2Pa，说明对岔角修圆能有效的缓解应力集中。但对卜型岔采用1.45的修圆系数修圆岔角，岔角处的部分应力由8.3Pa增至9.9Pa.说明修圆系数过大时，岔角处的相贯线太长，导致岔角处受力状态恶化，使得部分应力和衬砌的破坏体积都会增加。综上分析比拟说明，地下钢筋混凝土岔管