群英水库浆砌石重力拱坝三维有限元温度应力计算剖析.docVIP

群英水库浆砌石重力拱坝三维有限元温度应力计算剖析 　　摘要：对于坝高超过100 m的砌石拱坝，结构受力复杂，工程经验相对较少，特别是考虑到北方地区冬夏季节温差很大，将产生不可忽略的温度应力，因此有必要对大坝及相关部位进行三维有限元计算复核，并特别考虑了温度应力对结构受力的影响。大型通用有限元分析软件ANSYS温度计算模块为三维空间结构温度应力计算提供了极大的方便。通过对群英水库浆砌石重力拱坝各种工况的计算分析，认为在内外温差较大时，温度场在拱坝结构设计时应引起充分的重视。 　　关键词：温差；拱坝；温度应力；温度场；有限元 　　中图分类号：TV641文献标志码：A文章编号2015）002004604 　　群英水库建于河南省焦作市西北部22 km的修武县西村乡境内，[JP2]所在河流为海河流域卫河上游的大沙河上。水库主要建筑物有大坝、输水洞、电站厂房等。大坝为定圆心、变半径的浆砌石重力拱坝，坝顶弧长15428 m，坝顶中心角为80°，外半径1105 m，坝顶厚45 m，坝底厚52 m，厚高比为052，弧高比为154，坝顶高程49050 m，最大坝高1005 m，为亚洲同类型第一高坝。大坝采用坝顶溢流，溢洪道布置在坝顶中部，[JP]无闸控制，采用克奥型非真空溢流堰，挑流消能。 　　按照《砌石坝设计规范》，国内已建坝高超过100 m以上的砌石坝仅有3座，因此对于坝高超过100 m的砌石坝， 特别是考虑到北方地区冬夏季节温差很大，将产生不可忽略的温度应力。夏季受太阳辐射影响，坝体表面升温较高，而太阳辐射对水库水体的影响则不大。反之冬季将导致坝体表面温度较低，而水体的温度相对稳定。因此有必要对大坝及相关部位进行三维有限元计算复核，并特别考虑了温度应力对结构受力的影响。 　　1ANSYS温度应力计算理论基础 　　1.1稳态温度场有限元原理 　　根据模型边界条件，由热传导方程推导可得三维稳态温度场的导热方程为 　　由上述方程可知，温度应力的有限元计算与外荷载下普遍意义的有限元计算是相似的，温度变化引起的相当体力和相当面力可作为外荷载那样，由虚位移原理移置到各个节点上，即由温度场分析所得的节点温度输入到结构相应节点上，并输入初始温度场的各节点温度值，就完成了不均匀温度场的变化值输入，以此可进行温度应力的计算[3]。 　　2实例计算分析 　　2.1计算工况 　　拱坝设计荷载组合分为基本荷载组合和特殊荷载组合两类。本次结构安全分析根据群英水库大坝的实际情况和《砌石坝设计规范》（SL 25-2006）的要求，分析了以下4 种基本组合和4 种特殊组合工况。 　　（1）基本组合。 　　基本组合1：正常蓄水位（47700 m）+正常温降+坝体自重+扬压力； 　　基本组合2：设计洪水位（4817500 m）+正常温升+坝体自重+扬压力； 　　基本组合3：[JP2]死水位（42500 m）+正常温降+坝体自重+扬压力；[JP] 　　基本组合4：[JP2]死水位（42500 m）+正常温升+坝体自重+扬压力。[JP] 　　（2）特殊组合。 　　特殊组合1：校核洪水位（48520 m）+正常温升+坝体自重+扬压力； 　　特殊组合2：正常蓄水位（47700 m）+正常温降+坝体自重+扬压力+地震； 　　[JP2]特殊组合3：死水位（42500 m）+正常温降+坝体自重+[JP]扬压力+地震； 　　特殊组合4：死水位（42500 m）+正常温升+坝体自重+扬压力+地震。 　　2.2计算模型 　　根据相关规范和类似工程，利用有限元法对拱坝进行结构分析时，单元的剖分应尽量达到设计所要求的精度，单元形式结合拱坝体型合理选用。建模地基范围尽可能考虑大坝的特点和基础地质条件，满足计算要求。计算坐标系采用如下：X轴为顺河向，以向下游为正、向上游为负；Y轴为垂直向，与高程一致，以向上为正、向下为负；Z轴为坝轴线方向（横河向），以向右岸为正、向左岸为负。根据群英水库浆砌石重力拱坝的实际及坝基岩层分布情况，三维整体计算模型的范围如下：以位于河床最大坝高坝段中部的坝踵处为基准，顺河向分别沿坝踵、坝趾处往上下游方向取一倍坝高作为坝基长度；竖直方向沿坝底面处向下截取一倍坝高作为坝基深度；坝轴线方向以大坝坝肩两端为基准，向左右岸方向各延伸一倍坝高作为坝基宽度。同时，考虑计算的需要，该三维有限元模型忽略了廊道、输水洞等附属结构的影响。建模时，尽量模拟坝体的实际形态，使之反映实际情况。划分网格时，单元选用精确度较高的8节点六面体等参单元SOLID 45，单元节点总数为187 983个，单元总数为172 252个，群英水库拱坝的三维有限元模型图见图1。 　　2.3参数选取 　　根据大坝设计情况以及地质分布情况，其中坝体材料主要