泄洪底孔应力变形三维有限元剖析.docVIP

泄洪底孔应力变形三维有限元剖析 　　摘要:本文对重力坝应力变形计算的主要方法进行调研和分析对比,采用有限元方法,结合黄河炳灵水电站重力坝,利用大型通用有限元软件ANSYS,对泄洪底孔坝段建立了三维有限元模型,完成了应力变形计算分析。分别对6种工况下泄洪底孔坝段、弧门支承梁、不同剖面上的应力变形分布状况进行了对比分析,为评价泄洪底孔的结构适应性和设计方案的合理性、安全可靠性提供了理论依据。 　　关键词:重力坝泄洪底孔三维有限元应力变形 　　 　　 　　1.引言 　　结合黄河炳灵水电站工程,利用大型通用有限元软件ANSYS,对该水电站重力坝泄水底孔坝段的应力变形进行三维有限元分析计算,得出不同剖面上的应力变形分布状况,综合评价泄洪底孔的结构适应性和设计方案的合理性、安全可靠性,为设计提供可靠依据。 　　构建底孔坝段有限元模型,分别计算不同工况下泄洪底孔坝段、弧门支承梁等部分的应力变形分布规律。对计算结果进行对比分析,评价泄洪底孔在重力坝中的结构适应性及设计方案的安全可靠性。 　　2.三维有限元分析的基本原理 　　有限元法的基本原理是首先把一个连续体离散为有限个单元,这些单元仅在结点处连接,组成单元的集合体,依次替代原来的结构。然后,对各单元选择一个位移模式,使单元内部各点的位移可以近似用结点位移来表达,接着,根据总势能最小原则导出表示结点力和结点位移关系的单元刚度矩阵,再将单元刚度矩阵叠加起来,形成该系统的总刚度矩阵,根据给定的问题与应当满足的微分方程及边界条件,组成整个系统的平衡方程,解之,求出结点值,最后算出单元的应力与应变。 　　3.泄洪底孔三维有限元分析 　　有限元计算模型施加的结构边界约束为:在坝基底部施加三向约束,坝基上、下游边界面施加相应法向链杆约束,坝基和坝段左、右侧边界面施加相应法向链杆约束。 　　计算模型选取泄水底孔坝段作为计算对象,对坝段进行应力变形分析。计算域为坝基岩体取1倍的坝高,上游自坝踵向上游延伸1倍的坝高,下游自坝趾向下游延伸1倍的坝高。 　　整体坐标系的坐标原点在坝段左侧坝踵处。坝轴线指向右岸为X轴正向,下游方向为Y轴正向,铅直向上为z轴正向。三维有限元模型坐标系为笛卡尔直角坐标系,X轴顺水流方向,指向下游为正;Y轴铅直向上为正;Z轴垂直水流方向,指向右岸为正。泄洪底孔坝段三维有限元计算分析模型(未选取地基部分)如图-1所示。 　　闸墩等大体积混凝土用8节点棱柱体SOLID45单元模拟,板用壳SHELL93单元模拟,梁柱用BEAM4三维梁单元模拟,地基接触面用CONTA174和TARGE170单元模拟。模型单元总数为169626;节点总数为47851。泄洪底孔坝段三维有限元模型网格离散图见图-2所示。 　　 　　图-1泄洪底孔坝段三维计算分析模型 　　 　　图-2泄洪底孔坝段三维有限元模型网格离散图 　　4.应力变形计算结果 　　为了能清楚叙述泄洪底孔模型的应力变形分布情况,把整个结构分成泄洪底孔坝段、弧门 　　支承梁2部分,各部分的体形详见图-3和图-4。 　　 　　图-3泄洪底孔坝段体形图 　　图-4弧门支承梁体形图 　　 　　5.泄洪底孔坝段应力变形分析 　　 　　泄洪底孔坝段应力极值统计见表-1,轴向位移极值统计见表-2。 　　表-1泄洪底孔坝段应力极值统计表(MPa) 　　 　　 　　 　　表-2泄洪底孔坝段位移极值统计表(mm) 　　 　　 　　从泄洪底孔的坝段的应力极值的分布情况来看。在X、Y、Z方向的拉应力最大值都是由工况1来控制。由应力等值线图可以看出(由于篇幅限制此处略去剖面应力等值线图),各工况下泄洪底孔坝段的应力分布比较相似,且最大拉应力都出现在楔形块的扭面处,即在楔形块靠近基岩部位,可能是因为这个部位是混凝土楔形块与基岩相接触且形状有突变,从而引起了这一部位的应力集中。 　　从轴向位移极值统计表可以看出,坝段混凝土整体结构在X方向上的变形最大值1.114mm,出现在弧门支承梁上,这主要是因为弧门推力的作用。Y方向上的变形最大值2.846mm,出现在灌浆排水廊道顶部。Z方向上的变形最大值2.523mm,出现在坝轴线下10～14米的楔形块部位。 　　 　　6.弧门支承梁应力分析 　　 　　弧门支承梁应力极值统计见表-3,轴向位移极值统计见表-4。 　　可以看出工况二、工况三以及工况五这三种工况弧门大梁都施加了弧门推力,因此这三种工况下的弧门大梁的应力也比较大,弧门大梁的受力和配筋主要由这几种工矿来控制。而在这三种工况中,弧门大梁的正应力和剪应力相差不大,其中工况三的应力稍大。 　　由应力等值线图可以看出(由于篇幅限制此处略去剖面应力等值线图),X方向的最大正拉应力x发生在弧门推力作用的节