水利工程论文-高拱坝水垫塘反拱底板衬砌结构的非线性静力分析.doc

水利工程论文-高拱坝水垫塘反拱底板衬砌结构的非线性静力分析摘要：本文利用Ansys程序对高拱坝水垫塘衬砌结构进行了全过程计算。在止水触单元模拟水垫塘衬砌块结构之间的接缝以及衬砌与基岩的接触，并对衬砌块与拱端支座的连接形式进行了研究；在止水结构破坏后，衬砌块结构下部受到方向向上的脉动压力和时均压力，拱的作用表现出来，并且衬砌块与基岩之间的锚固钢筋开始作用，这时的结构可看作是一种“反吊拱”。本文利用弹簧单元模拟衬砌块与基岩之间的钢筋进行计算，得到一些用于指导水垫塘设计的结果。关键词：水垫塘衬砌接触元弹簧元弹性地基梁(板)峡谷地段修建的高拱坝，其泄洪消能布置往往要设置水垫塘，水垫塘衬砌的底板有平底板和反拱底板两种，平底板的稳定问题，现在已研究得比较清楚，反拱底板在近几年得到广泛应用，主要是因为底板拱结构抵抗破坏的能力强，稳定性要优于平底板。但是由于反拱底板面积比较大，在施工过程中要设置温度缝和施工缝，因此反拱底板被分成一系列相互独立又相互联系的板块，其受力过程表现为一种高度的非线性，特别是在止水破坏后，各板块相互撞击、滑动表现为典型的接触行为，计算相当复杂。崔广涛等从理论分析和模型实验两方面论证了反拱型底板的受力条件好，其稳定性优于平底板；刘沛清等也对反拱型底板的稳定性进行了研究，并提出了相应的稳定计算模式。但是他们只是把衬砌块作为刚体，把块间连接看作“铰”(其实是一种机构)，或者把整个水垫塘结构看作一个三铰拱或无铰拱，和实际的块间结构有很大的区别，本文利用接触单元来模拟相邻板块间的接触，接触单元可以模拟块间的接触、咬合、摩擦、分离、撞击等不同状况，特别是在止水结构破坏后，衬砌块结构下部受到方向向上的脉动压力和时均压力，拱的作用表现出来，并且衬砌块与基岩之间的锚固钢筋开始作用，这时的反拱底板整体可以看作是一个“反吊拱”。利用弹簧单元来模拟锚固到基岩中的钢筋，进行有限元计算，得到比较理想的结果。1有限元单元形式的选择以及非线性接触理论接触是一种高度的非线性行为，由于考虑摩擦使问题变得困难起来。接触单元是覆盖在模型接触面上的一层单元，有限元模型通过指定接触单元来识别可能的接触匹对，本文所使用的接触单元为三节点单元，如图1所示。这种单元形式不受接触形状的限制，使用范围广，这种单元实际上是一种点对线接触，k点为接触点，ij为目标线。不过可以用它来模拟线对线的接触，因为可以把线指定为一组节点。对于水垫塘衬砌块，由于其相邻块之间地位相同，可以用对称接触来模拟，即先选择a线作为目标线，b线上的点作为接触点，形成接触对，然后再反过来利用b线作为目标线，a线上的点作为接触点，形成另一组接触对。图1接触单元示意图2单元坐标单元坐标如图2所示，单元法向量和单元切向量分别为：n=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方向的运动受到目标面摩擦力的限制，本文中的摩擦类型选为库仑摩擦类型，程序中提供了一个不管接触压力的值，而人为提供等效剪应力的选项，如果剪应力达到之，滑动发生。在求摩擦力之前，需要求出接触点相对于目标面的切向位移us=1/2(s*-s0*)l(5)式中：s0*上步求解过程中的接触位置。其实这个位移是由粘性位移和滑动位移两部分组成，即为粘性位移，为滑动位移。摩擦力为：在粘性阶段：fs=ktFfs在滑动阶段：fs=fs(6)式中：kt:粘结刚度，fs库仑摩擦类型的静态摩擦限值(fs=-fn)，为摩擦系数，F为静态到动态转化摩擦因子。向力可以通过罚函数法或罚函数和拉格朗日联合法来实现，这时力将加到节点上，直接触点穿透进入目标面。如图3所示。罚函数法：(7)(8)式中：kn为法向接触刚度。罚函数和拉格朗日联合法：fn=min(0,kng+i+1)(9)式中：i+1为节点在i+1处的拉格朗日增值因子，当g时，i+1=i。为容限因子，为材料因子(0)。为方便求单元的刚度矩阵和荷载向量，引入两个向量(10)(11)由图3在法向上列平衡方程是：fn,k=fn,i+fn,j=fn(12)同理，由切向得：fs,k=fs,i+fs,j=fs(13)写成向量的形式为：(14)单元的刚度矩阵为：在粘性阶段：Ke=KnNnNnT+KsNsNsT(15)在滑动阶段：Ke=KnNnNnT(16)在未接触阶段：Ke=0(17)式中：Ks为切向接触刚度。弹簧单元能在轴向上承受拉压力和扭矩，本文所用的是二维弹簧元，每个单元由两个节点构成，每个节点有两个自由度。2计算实例本文结合溪洛渡水垫塘进行计算，