某水电站溢洪道闸室正常运行期三维有限元分析.docVIP

某水电站溢洪道闸室正常运行期三维有限元分析 - 工程设计 简介： 为了验证溢洪道闸室结构在正常运行期的安全性与合理性，本文运用三维有限元方法，对闸室结构的位移、应力、抗滑稳定等情况进行了计算分析。计算结果表明：闸室的位移、应力、抗滑稳定等情况均满足相关要求，闸室结构安全、体型合理。关键字：溢洪道闸室 三维有限元 结构设计1、概述 　　某溢洪道闸室长50 m，堰上游宽81.45 m，下游宽72 m。堰顶高程209 m，由4孔15m宽 ×21.09 m高的设闸溢流堰组成，堰体上游坡度为1：0.667，堰面曲线为Y=0.04285X1.85。闸墩末端宽度为4 m，最宽处约6.5 m，闸墩采用预应力混凝土结构。闸室设弧形工作门和钢叠梁检修门，门库设在左侧。基础齿槽高程190.00 m，在齿槽内设帷幕灌浆检查排水廊道。闸墩顶设交通桥。溢洪道闸室构筑在岩石地基上，正常运行期水位为228.00m。　　该闸室存在结构复杂、闸墩较高等特点，而传统计算方法难以反映截面突变、刚度变化等因素对力学性能的影响，也不能准确描述关键部位的应力状态和变形情况，因此，需要在传统计算方法外辅以有限元法进行校核补充。本文采用三维有限元方法研究了溢洪道闸室在正常运行期的位移、应力、抗滑稳定等情况，为闸室的结构设计提供了参考依据。　　2、计算模型　　2.1 三维有限元模型　　溢洪道闸室三维有限元计算模型的计算范围为：闸室上游侧取1.5倍闸室高度，下游侧取2.0倍闸室高度，左右两侧和基础分别取一倍闸室高度，闸室高度取46.5m。基岩与闸室混凝土按固结处理。　　采用空间六面体和四面体等参单元对整体结构进行网格剖分。网格剖分时主要参照以下四个原则进行[1]：　　（1） 在现有计算机内存和硬盘等外部条件限制下，尽可能多地增加单元和节点数量，以提高计算精度。　　（2）溢流堰、中间三个闸墩、牛腿等部位基本采用全六面体网格剖分；左右边墩及挡水坝结构复杂，采用四面体网格剖分；基础岩体采用四面体网格剖分。　　（3） 基础岩体的单元尺寸由构筑物边界到模型边界处由小到大过渡，以提高构筑物与岩体接触面附近的计算精度。　　（4）网格剖分应尽量反映构筑物的轮廓形状、材料分区、荷载分布等情况。　　按照上述网格剖分原则，溢洪道闸室结构三维有限元模型的总单元数为723955个，节点数为209184个。溢洪道闸室结构三维有限元网格剖分图见图1。采用右手坐标系，坐标原点位于闸室轴线竖直面、溢洪道中心线竖直面、V190.00m高程水平截面的交点处。X轴正向指向顺水流方向，Y轴正向指向竖直方向向上，Z轴正向指向横河流方向。　　2.2 混凝土分区及材料特性　　计算分析时，闸室混凝土主要分为三部分，溢流堰及左右挡水坝采用C20混凝土，闸墩采用C25混凝土，弧形闸门支座牛腿采用C40混凝土。各区混凝土静态物理及力学参数见表1。闸室构筑在新鲜砂岩上，闸室下游为新鲜页岩，砂岩和页岩静态物理及力学参数见表2。　　　　图1溢洪道闸室结构三维有限元网格剖分图　　表1 混凝土静态物理及力学参数标号容重（N/m3）弹性模量E（Pa）泊松比ν轴心抗压强度标准值（Pa）轴心抗拉强度标准值（Pa）C202.40×1042.55×10100.16713.50×1061.50×106C252.40×1042.80×10100.16717.00×1061.75×106C402.40×1043.25×10100.16727.00×1062.45×106　　表2 基岩静态物理及力学参数岩体种类比重（kg/m3）弹性模量E（Pa）泊松比ν饱和抗压强度（Pa）砂岩2.60×1038.00×1090.18100.00×106页岩2.60×1037.00×1090.2045.00×106　　2.3 荷载及边界条件　　闸室正常运行期的荷载主要有结构自重、静水压力、弧形闸门上游水重、扬压力、预应力、设备自重等。模型边界采用刚性链杆约束[2，3]。　　3、计算结果及分析　　采用三维有限元方法，对闸室结构的位移、应力、抗滑稳定等情况进行了计算分析。计算分析内容包括整体结构的位移结果、关键部位的位移结果、门槽部位路径位移结果；整体结构的应力结果、关键部位的应力结果、闸墩根部的路径应力结果、关键截面的八种应力结果等。受篇幅限制，本文仅给出有代表性的重要结果。　　3.1 闸墩及门槽标号说明 　　由于闸室结构较为复杂，闸墩和门槽数量较多，为了在计算结果及分析中方便描述，各部位分别用相应的符号代替。从左至右的五个闸墩编号为：D1、D2、D3、D4、D5；从左至右的八个检修闸门门槽编号为：J1、J2、J3、J4、J5、J6、J7、J8；从左至右的八个弧形闸门门槽编号为：H1、H2、H3、H4、H5、H6、H7、H8。闸墩及门槽编号说明图见图2。　　　　图2 闸墩及门槽编号