《高拱坝水垫塘反拱底板衬砌结构的非线性静力分析》由会员分享，可在线阅读，更多相关《高拱坝水垫塘反拱底板衬砌结构的非线性静力分析（12页珍藏版）》请在金锄头文库上搜索。
1、1高拱坝水垫塘反拱底板衬砌结构的非线 性静力分析摘要：本文利用 Ansys 程序对高拱坝水垫塘衬砌结构进行了全过程计算。在止水结构未破坏以前，衬砌块结构可以看作是作用在弹性地基上的板，利用点 线接触单元模拟水垫塘衬砌块结构之间的接缝以及衬砌与基岩的接触，并对衬砌块与拱端支座的连接形式进行了研究；在止水结构破坏后，衬砌块结构下部受到方向向上的脉动压力和时均压力，拱的作用表现出来，并且衬砌块与基岩之间的锚固钢筋开始作用，这时的结构可看作是一种“反吊拱”。本文利用弹簧单元模拟衬砌块与基岩之间的钢筋进行计算，得到一些用于指导水垫塘设计的结果。 关键词：水垫塘 衬砌 接触元 弹簧元 弹性地基梁(板) 峡
2、谷地段修建的高拱坝，其泄洪消能布置往往要设置水垫塘，水垫塘衬砌的底板有平底板和反拱底板两种，平底板的稳定问题，现在已研究得比较清楚，反拱底板在近几年得到广泛应用，主要是因为底板拱结构抵抗破坏的能力强，稳定性要优于平底板。但是由于反拱底板面积比较大，在施工过程中要设置温度缝和施工缝，因此反拱底板被分成一系列相互独立又相互联系的板块，其受力过程表现为一种高度的非线性，特别是在止水破坏后，各板块相互撞击、滑动表现为典型的接触行为，计算相当复杂。崔广涛等从理论分析和模型实验两方面论证了反拱型底板的受力条件好，其稳定性优于平底板；刘沛清等也对反拱型底板的稳定性进行了研究，并提出了相应的稳定计算模式。但是
3、他们只是把衬砌块作为刚体，把块间连接看作“铰”(其实是一种机构)，或者把整个水垫塘结构看作一个三铰拱或无铰拱，和2实际的块间结构有很大的区别，本文利用接触单元来模拟相邻板块间的接触，接触单元可以模拟块间的接触、咬合、摩擦、分离、撞击等不同状况，特别是在止水结构破坏后，衬砌块结构下部受到方向向上的脉动压力和时均压力，拱的作用表现出来，并且衬砌块与基岩之间的锚固钢筋开始作用，这时的反拱底板整体可以看作是一个“反吊拱”。利用弹簧单元来模拟锚固到基岩中的钢筋，进行有限元计算，得到比较理想的结果。1 有限元单元形式的选择以及非线性接触理论接触是一种高度的非线性行为，由于考虑摩擦使问题变得困难起来。接触单
4、元是覆盖在模型接触面上的一层单元，有限元模型通过指定接触单元来识别可能的接触匹对，本文所使用的接触单元为三节点单元，如图 1 所示。这种单元形式不受接触形状的限制，使用范围广，这种单元实际上是一种点对线接触，k 点为接触点，ij为目标线。不过可以用它来模拟线对线的接触，因为可以把线指定为一组节点。对于水垫塘衬砌块，由于其相邻块之间地位相同，可以用对称接触来模拟，即先选择 a 线作为目标线，b 线上的点作为接触点，形成接触对，然后再反过来利用 b 线作为目标线，a 线上的点作为接触点，形成另一组接触对。图 1 接触单元示意3图 2 单元坐标单元坐标如图 2 所示，单元法向量和单元切向量分别为：n
5、=vs(1)s=(xj-xi)/l(2)式中：v:ij 与整体坐标 x-y 平面的方向角；xi、xj：i、j 点的位置向量；l:目标单元的长度。接触点与目标面的距离为：g=(xk-xi)Tn(3)s*=-1+2(xk-xi)Ts/l(4) 式中：xk:k 点的位置向量。当两板块接触时 k 点沿 n 方向运动，受到目标面法向应力的限制，沿 s 方向的运动受到目标面摩擦力的限制，本文中的摩擦类型选为库仑摩擦类型，程序中提供了一个不管接触压力的值，而人为提供等效剪应力的选项，如果剪应力达到之，滑动发生。在求摩擦力之前，需要求出接触点相对于目标面的切向位移4us=1/2(s*-s0*)l(5)式中：s
6、0*上步求解过程中的接触位置。其实这个位移是由粘性位移和滑动位移两部分组成，即为粘性位移，为滑动位移。摩擦力为：在粘性阶段：fs=ktFfs在滑动阶段：fs=fs(6)式中：kt:粘结刚度，fs 库仑摩擦类型的静态摩擦限值(fs=-fn)， 为摩擦系数，F 为静态到动态转化摩擦因子。向力可以通过罚函数法或罚函数和拉格朗日联合法来实现，这时力将加到节点上，直接触点穿透进入目标面。如图 3 所示。罚函数法：5(7)(8)式中：kn 为法向接触刚度。罚函数和拉格朗日联合法：fn=min(0,kng+i+1)(9)式中：i+1 为节点在 i+1 处的拉格朗日增值因子，当g 时，i+1=ii+kng；当
7、g 时，i+1=i。 为容限因子， 为材料因子(0)。为方便求单元的刚度矩阵和荷载向量，引入两个向量(10)(11)由图 3 在法向上列平衡方程是：6fn,k=fn,i+fn,j=fn(12) 同理，由切向得：fs,k=fs,i+fs,j=fs(13)写成向量的形式为：(14)单元的刚度矩阵为：在粘性阶段：Ke=KnNnNnT+KsNsNsT(15)在滑动阶段：Ke=KnNnNnT7(16)在未接触阶段：Ke=0 (17)式中：Ks 为切向接触刚度。弹簧单元能在轴向上承受拉压力和扭矩，本文所用的是二维弹簧元，每个单元由两个节点构成，每个节点有两个自由度。 2 计算实例本文结合溪洛渡水垫塘进行计
8、算，工程概况为：水库的正常蓄水位 600.00m，相应库容为 115.7 亿 m3，最大坝高 278m，装机容量初步拟定 12600MW。拦河大坝为混凝土双曲拱坝，最大坝高 278m；泄水建筑物有坝身 7 个表孔(12.511m)和八个深孔(66.7m)，具有“高水头、大泄量、窄河谷”的特点，最大泄洪功率高达 100 000MW，60%的洪水通过坝身渲泄，坝身泄量达 30000m3/s，坝身的表孔和深孔的最大单宽流量超过 200m3/s。水垫塘衬砌每块园心角 8.64,半径 81.25m，分为 9 块，厚度为 4m。在止水结构未破坏以前，衬砌块承受向下的水荷载作用，其实可看作是作用在弹性地基上
9、的梁(或板)，利用接触单元模拟板块之间的连接，用四边形单元模拟板块，得到有限元模型如图 4，当止水结构破坏后，水垫塘底板受到向上的脉动压力和时均压力，拱的作用表现出来，并且衬砌块与基岩之间的锚固钢筋开始作用，利用四边形单元模拟板块，接触单元模拟板块之间的连接，弹簧单元模拟锚固锚筋，得到有限元模型如图 5。计算结果如表 1、2、3、4 和图 6 所示。8图 3 接触单元受力分析图 4 止水破坏前有限元图 5 止水破坏后有限元图 6 止水破坏后反吊拱在 40 米水压力下的变位表 1 拱端设缝情况下的拱端位移(止水破坏前)表 2 拱端未设缝情况下的拱端推力(止水破坏前)9工况 位移/mm温升(-54
10、0) 0.86温降(40-5) -0.86静水压-0.85静水压+10m 冲击水压-0.96工况 推力/(T/m)温升(-540) 10-70.2180温降(40-5) 70.2176静水压71.3792静水压+10m 冲击水压 85.164640m 脉动压力 -284.7316表 3 止水破坏前拱端设缝与未设缝情况下的主应力大小位置比较设缝未设缝11计算工况第一主应力第三主应力第一主应力第三主应力大小/(MPa) 大小/(MPa) 位置 大小/(MPa) 位置 大小/(MPa) 位置温升(-540) 0.012 -0.434 5 0.300 1 -0.270125温降(40-5) 0.012
11、 -0.434 5 0.300 1 -0.270 5静水压0.016 -0.435 5 0.349 1 -0.280 5静水压+10m 冲击水压0.031 -0.438 5 0.551 1 -0.379 5注：表中位置数是指距拱端第几块表 4 反吊拱衬砌块需要提供的锚固力13块数1 2 3 4 5锚固力/(T/m)301.54 343.53 336.74 322.22 319.98锚筋剪力/(N/m)773 424.13 117.49 32.77 0板块位移/mm1.878 5.085 4.825 4.692 4.644平底板锚固力/T407.92 407.9314407.93 407.93
12、407.93注：表中块数是指距拱端所说，锚筋 2810003 结论与讨论通过本文的计算可以得到如下结论：(1)由表 3 知，对于底板混凝土强度是没有问题的。(2)由表 2 知，在温升和止水结构破坏后，拱端承受的推力是相当大的，这和试验数据项一致2。因此应该在拱端块和拱端支座之间设缝，这将使拱端力减小，但是缝宽不能太大，太大反拱起作用太慢，并且由表 1 得，缝宽在 1cm 即可。(3)由图 6 可知，对于止水完全破坏的反吊拱，最大位移出现在距坝端的第二块。这已由模型实验所验证2。所以这两块应是设计中锚固的重点。(4)反拱底板的锚固力比平底板要少 20%以上，反拱底板是一种安全经济合理的结构形式。(5)用接触元模拟水垫塘衬砌块之间的接缝是可行的。参 考 文 献：1 朱伯芳.有限元原理与应用(第二版)M.北京：中国水利水电出版社，1998，10.2 练继建，杨令强，等.溪洛渡水垫塘水弹性模型实验(中间报告)R.天津：天津大学，2000，12.3 刘沛清.泄水建筑物消能防冲论文);“论文选集M.北京：中国科学技术出版社，2000. 4 Johnson K L. Contact MechanicsM.Cambridge University Press,1985. 5 Kohnke Peter.Ansys Theory Reference Release5，6R.1999.
