三峡下游坝面钢衬钢筋混凝土管道有限元分析.docx

三峡下游坝面钢衬钢筋混凝土管道有限元分析工程概况三峡水电站为巨型水电站，坝型为混凝土重力坝。坝高145m，坝顶高程185m，正常蓄水位175m，坝段宽24m。采用坝后式厂房，装机26台，采用单机单管形式，单机容量700MW。电站压力管道管径为12.4m，HD值为1730m2。压力管道采用坝后背管布置形式，在下游坝面预留浅槽，管道埋入坝面以内约1/2管径。压力管道背管部分采用钢衬钢筋混凝土管形式，管道规模大，结构新颖，技术复杂，管道工作条件十分复杂。材料及参数管道钢衬采用16MnR钢板，板厚30mm，E=210 GPa，μ=0.3，ρ=0.0078Mkg/m3。管道外包混凝土厚度为2.0m，标号采用C30，E=25.5GPa，μ=0.167，ρ=0.0025Mkg/m3。坝体混凝土标号C20，E=20.0GPa，μ=0.167，ρ=0.0024Mkg/m3。计算模型3.1计算断面计算中取12m管道作为研究对象，计算时不考虑大坝下游坡度，按照管道轴线处于水平位置计算。计算断面尺寸（单位：m）见图1。图1 三峡下游坝面钢衬钢筋混凝土管道尺寸约束条件假定混凝土未开裂，钢衬与混凝土之间没有缝隙，坝段纵缝可允许坝段混凝土自由变形，所以模型底面位移全约束，管道轴线方向位移被约束，其他不做约束。单元类型及网格划分管道内衬用Shell Elastic 4node 63单元（四面体平面板壳单元），Real Constants厚度设置为0.03m。材料参数E=210 GPa，μ=0.3，ρ=0.0078Mkg/m3。管道外包混凝土和坝体混凝土用Solid Brick 8node 185单元（8结点六面体块体单元），但其材料参数不同。管道外包混凝土E=25.5GPa，μ=0.167，ρ=0.0025Mkg/m3。坝体混凝土E=22.5GPa，μ=0.167，ρ=0.0025Mkg/m3。计算模型共12365个结点，11380个单元。计算模型网格见图2，钢衬网格如图3所示。图2计算模型网格图3钢衬网格计算方案在钢衬内表面施加大小为0.5MPa的内水压力，施加在shell单元节点上，如图4所示。图4 约束及荷载计算结果分析5.1变形分析5.1.1运行工况变形混凝土管道在施加内水压力以后的变形情况如图5所示。其中最大变形发生在管道内表面顶部稍偏左右两侧区域，其最大值为0.568mm。图5 管道及坝体混凝土变形5.1.2 位移云图混凝土管道X、Y、Z方向位移及总位移云图见图6~图9。图6 X方向位移云图图7 Y方向位移云图图8 Z方向位移云图图9 总位移云图5.1.3 位移分析说明由以上位移云图可知，X方向最大位移发生在钢衬混凝土腰部附近区域，其值为0.457mm；Y方向最大位移发生在管道顶部附近区域，其值为0.526mm；Z方向位移总体表现为在管道整个断面上基本为0；管道总位移的最大值位于管道内表面顶部稍偏左右两侧区域，其值为0.568mm。另外，原结构以面YOZ为对称面，所施加的荷载同样关于面YOZ对称，因此推断其各向位移也应该关于面YOZ对称。从位移云图可以看出，各向位移均具有关于面YOZ的对称性，据此可以认为位移成果是可信的。5.2 应力分析5.2.1混凝土应力等值线图混凝土应力分析分别考虑X、Y、Z方向的正应力，以及第一、二、三主应力。具体混凝土应力情况见图10~图15。图10 X方向正应力云图图11 Y方向正应力云图图12 Z方向正应力云图图13 第一主应力云图图14 第二主应力云图图15 第三主应力云图5.2.2钢衬Mises应力云图钢衬X、Y、Z方向正应力和第一、二、三主应力云图从略，此处给出钢衬Mises应力云图，见图16所示。图16 钢衬Mises应力云图5.2.3应力分析说明　数值（MPa）位置　拉压拉压X方向正应力1.7390.365管道顶部内表面管道腰部Y方向正应力1.6510.388钢衬混凝土左右两边腰部附近区域管道底部内表面Z方向正应力0.2770.046钢衬混凝土左右两边腰部附近区域计算模型下部左右两个拐角处第一主应力1.7390.037管道顶部内表面计算模型下部左右两个拐角处第二主应力0.2770.046钢衬混凝土左右两边腰部附近区域计算模型下部左右两个拐角处第三主应力0.0320.388坝体底部管道内表面腰部断面内外力校核6.1计算方法进行腰部断面的内外力平衡校核，计算方法如下：沿管道腰部取截面Ⅰ—Ⅰ，见图17。图17 计算模型示意图计算参数：P—钢衬受到的内水压力，A—表示腰部断面。理论上，应在ANSYS整体坐标系下对截面应力σy积分，与自重、水压力在Y方向分力比较。此处为排除重力的影响，不对混凝土、钢衬的自重进行计算，只考虑加上水压力（管道中考虑水的自重）后在钢衬和混凝土中引起的附加应力。又由于钢衬