1、基于水电站钢衬钢筋混凝土坝后背管有限元分析【摘要】采用非线性有限元法对水电站钢衬钢筋混凝土压力管道进行了分析，分析外包钢筋混凝土的裂缝开裂次序以及钢衬的承载比例等，给出了钢衬和外包钢筋混凝土的应力、位移分布规律。为水电站钢衬混凝土压力管道的设计与施工提供了一定的参考依据。【关键词】水电站;有限元法;钢衬钢筋混凝土管;承载比例系数一、钢筋混凝土单元ANSYS 软件的 SOLID65 单元是专为混凝土、岩石等抗压能力远大于抗拉能力的非均匀材料开发的单元，它可以模拟混凝土中的加强钢筋，以及材料的拉裂和压碎现象。（一）单元中混凝土模型的本构关系SOLID65 单元是一种 8 节点的块体单元，它是在普通
2、 8 节点三维等参单元的基础上增加了针对混凝土材料参数和整体式钢筋模型。在 SOLID65 单元中使用弹塑性本构关系来描述混凝土受拉时的应力应变关系，采用 Drucker-Prager 屈服准则。其塑性流动为关联流动，混凝土单轴受压时的应力-应变关系采用 Saenz 提出的公式：=E01+E0Es-20+02JZ）JY（1）式中：E0 为初始弹性模量;Es=0/0 为应力达到峰值时的割线弹性模量;0，0分别为应力-达到峰值时的应力，应变。（二）压碎与开裂模拟使用 SOLID65 单元进行非线性计算，当应力组合达到破坏面时，则单元进入压碎或开裂状态。如果在单轴、双轴或三轴压力作用下，某个积分点上
3、的材料失效了，就意味着这个点上的材料压碎了。单元进入压碎状态，则单元刚度为零，且应力完全释放。通过修正的应力-应变关系，引入垂直与裂缝表面方向上的一个缺陷平面来表示在某个积分点上出现了裂缝。当裂缝张开时，后继荷载产生了在裂缝表面的滑动或剪切时，引入一个剪切力传递系数t 来模拟剪切力的损失。当裂缝闭合时，所有垂直于裂缝面的压应力都能传递到裂缝上，但是剪切力只传递原来的c 倍。在抗裂计算过程中，裂缝闭合的判据为：开裂应变ckck0。其中，开裂应变ckck 的定义为：JZ（ckck=ck1+1-（ck2+ck3）只有一条裂缝ck1+ck2 有两条裂缝ck1 有三条裂缝JZ）JY（3）其中，ck=Tc
4、k，式中：Tck为坐标转换矩阵。JZ（=eln-1+n-thn-plnJZ）JY（4）式中：n 为荷载步数;eln-1 为前一步弹性应变;n为应变增量;thn为热应变增量;pln为塑性应变增量。如果ckck0，则认为裂缝是张开的，在某个积分点上出现了裂缝之后，则认为在下一步迭代中裂缝是张开的。二、结构计算模型（一）工程概况混凝土坝后的钢衬钢筋混凝土压力管道是近年来发展起来的水电站压力管道的新形式。这种压力管道是一种按钢衬与鋼筋混凝土联合承载设计的复合型结构。鉴于此种结构的复杂性，将采用非线性有限元法对水电站坝后钢衬钢筋混凝土压力管道进行结构分析，给出钢衬和外包钢筋混凝土的应力、位移分布规律，分
5、析外包钢筋混凝土的裂缝开裂次序以及钢衬的承载比例等。参与本内容研究的水电站位于贵州省清镇市和贵州省黔西县交界处，是利用乌江干流水电阶梯开发的第二级东风水电站已有建筑物兴建的。水电站最小水头 95m，最大水头132m，设计水头 117m，装机容量 125MW。电站拦河大坝为双曲率抛物线薄拱坝，最大坝高162m，坝后背管采用钢衬钢筋混凝土压力管道，钢衬内半径 r=2.50m，钢衬厚度=20mm，外包钢筋混凝土厚度 1m。（二）计算方案在对背管结构进行非线性有限元分析时，考虑了两种计算方案：方案 1，采用钢衬钢筋混凝土压力管道结构形式;方案 2，采用预应力钢衬钢筋混凝土压力管道结构形式。管道设计内水
6、压力1.17MPa。（三）基本参数水电站坝后钢衬钢筋混凝土压力管道钢衬材料为 16MnR 钢，泊松比=0.3，弹性模量E=206GPa，抗拉强度设计值s=300MPa。外包钢筋混凝土厚为 1.0m，强度等级为 C25，泊松比=0.167，弹性模量 E=28GPa，轴心抗压强度设计值 fc=12.5MPa，轴心抗拉强度设计值 ft=1.30MPa。混凝土管环向配置两层钢筋，钢筋强度等级为级，每层为40200。钢管和外包混凝土之间初始缝隙为 0.5mm，用以模拟施工、温度变化和混凝土徐变的影响。（四）计算模型在坝后背管斜直段切出 1m 长的管段进行计算分析。对计算单元划分时，钢衬采用 4节点的壳体
7、单元来模拟。外包混凝土及坝体采用 8 节点的块体单元模拟。普通钢筋和预应力钢绞线采用杆单元模拟。钢衬和外包混凝土留有的缝隙用弹簧单元连接，管坝之间的垫层也用弹簧单元来模拟。三、计算结果分析（一）钢管应力从图 1 可以看出，钢管在 90处环向应力随内水压力变化曲线分为 3 部分，内水压力在 0.7MPa 以内变化时，钢管在 90处环向应力基本呈线性变化。内水压力在 0.7MPa 至1.7MPa 范围内变化时，钢管在 90处环向应力变化较为平缓，这是由于此时内水压力由钢衬和钢筋混凝土管共同承担，使钢管环向应力增加较慢。内水压力在 1.7MPa 至 2.4MPa范围内变化时，钢管在 90处环向应力变
8、化较为剧烈，这是由于此时外包混凝土管已经出现贯通性裂缝，使混凝土管的承载能力降低，钢管分担荷载增加较快。预应力钢衬钢筋混凝土压力管道方案比普通钢衬钢筋混凝土压力管道方案钢管的环向应力小得多，并且预应力钢衬钢筋混凝土压力管道裂缝出现的较晚。（二）混凝土应力在设计内水压力 1.17MPa 下，从图 2 可以看出，混凝土管在 90处环向应力内缘较小，外缘较大，在 0至 45范围内环向应力内缘较大，外缘较小，环向应力变化较大。但在 0至 270范围内的混凝土环向应力较小，并且环向应力变化较为平缓。（三）钢筋应力钢衬钢筋混凝土压力管道内、外层环钢筋在 90处环向应力随内水压力变化曲线如图3、图 4 所示
9、，得出内、外层环钢筋在 90处环向应力随内水压力变化规律同钢管环向应力随内水压力变化规律类似。但在预应力钢衬钢筋混凝土压力管道方案中，内、外层环钢筋在内水压力为零的情况下，钢筋在 90处环向应力为压应力。并且内层环钢筋比外层环钢筋较早进入塑性状态。（四）混凝土开裂规律当内水压力达到设计内压 1.17MPa 时，钢衬的最大环向拉应力为 91.05MPa，外裹混凝土管的最大环向拉应力为 1.12MPa。当内水压力达到设计内压 1.77MPa 时，外裹混凝土管出现第一条初裂缝，初裂缝出现的位置大约在=35和 145的混凝土管内缘上。当内水压力达到设计内压 1.96MPa 时，外裹混凝土管出现第一条贯
10、穿性裂缝，贯穿性裂缝出现在大约=50和 130的混凝土管上。随着内水压力的增加，贯穿性裂缝条数迅速增加，当内水压力达到 2.08MPa 时，上半圆混凝土基本上都已裂穿。当内水压力达到设计内压2.12MPa 时，外裹混凝土管在-50230 范围内都已开裂。（五）钢管承载比例钢衬钢筋混凝土压力管道钢衬承载比例系数如表 1 所示，钢衬承载比例系數定义为钢衬承担的内水压力与总的内水压力的比值，钢衬承载比例系数先减小后略微增加并趋于平缓。四、结语分析外包钢筋混凝土的裂缝开裂次序以及钢衬的承载比例，采用非线性有限元法对水电站坝后钢衬钢筋混凝土压力管道进行结构分析，给出钢衬和外包钢筋混凝土的应力、位移分布规
11、律，通过提高压力管道的安全系数，对钢管具有保护作用。钢衬钢筋混凝土压力管道是一种受力较为合理的结构形式，外包钢筋混凝土和钢衬联合承载，预应力钢衬钢筋混凝土压力管道可以有效的延缓裂缝的出现，具有广阔的应用前景，为水电站钢衬混凝土压力管道的设计与施工提供了一定的参考依据。参考文献：1姬栋宇，童丽萍.鸭池河水电站钢衬混凝土坝后背管结构分析J.人民黄河，2009（12）.2中华人民共和国水利部.SL/T191-96 水工混凝土结构设计规范S.北京：中国水利水电出版社，1996.3国茂华，侯建国，张长青.坝后背管外包混凝土裂缝间距计算分析J.人民黄河，2009（11）.4马文亮，王清云，张建华.预应力钢衬钢筋混凝土坝后背管非线性分析J.长江科学院院报，2012（9）.作者简介：谢静仪（1991-），男，湖南益阳人，助教，湖南城建职业技术学院土木工程系，研究方向：现场施工管理及土木工程教学。