水利工程论文-水电站蜗壳保压浇混凝土结构的三维仿真分析.doc_第1页
水利工程论文-水电站蜗壳保压浇混凝土结构的三维仿真分析.doc_第2页
水利工程论文-水电站蜗壳保压浇混凝土结构的三维仿真分析.doc_第3页
水利工程论文-水电站蜗壳保压浇混凝土结构的三维仿真分析.doc_第4页
水利工程论文-水电站蜗壳保压浇混凝土结构的三维仿真分析.doc_第5页
已阅读5页,还剩2页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

水利工程论文-水电站蜗壳保压浇混凝土结构的三维仿真分析摘要:三峡工程水电站厂房蜗壳采用保压浇外围混凝土的结构形式。为研究钢蜗壳与外围混凝土交界面的接触性态,分别对冬季和夏季浇混凝土情况进行了模拟施工过程的三维有限元仿真计算,给出了交界面在不同季节不同水位运行期的传力和间隙,结果表明温度对传力的影响显著。对冬季浇筑情况,研究了通过提高保压水温来减小高温季节高水位运行期的传力;对夏季浇筑情况,研究了通过降低保压水头来减小蜗壳混凝土在低温季节低水位运行期的间隙。关键词:蜗壳保压浇混凝土温度效应仿真分析三维有限元法温度变化引起的混凝土和钢蜗壳间的不协调温度变形会影响两者间的相互作用,这种相互作用是随运行季节变化的,并与浇筑季节和保压水温有关,因此选择合适的浇筑季节和保压水温也同样重要。若再考虑施工过程及混凝土弹模、徐变和水泥水化热随时间变化的影响,两者的相互作用实际上是十分复杂的,只有通过模拟施工过程的仿真计算才能分析清楚。然而,通常都没有考虑温度的影响1,模拟施工过程的仿真分析更是少见。三峡水电站蜗壳采用保压浇混凝土的结构形式,装机高程57m.正常设计水位175m(考虑水击力的影响,相应内水压力为1.395MPa),初期运行水位135m,两者相差较大。通过有关研究工作,设计采用的保压水头为70m.由于机组台数多和工期要求,可能在不同季节都会有机组施工。一方面,由于保压水头较小,低温季节浇外围混凝土的机组在高温季节高水位运行时,混凝土结构承担的荷载较大;另一方面,由于初期运行水头较小,高温季节浇外围混凝土的机组在低温季节低水位运行时,钢蜗壳与混凝土将不能完全贴紧运行。这两方面可分别采取提高保压水温和降低保压水头的措施。为此,本文通过三维有限元仿真计算,分析冬、夏季浇筑外围混凝土情况,在不同季节、不同水位运行期钢蜗壳与外围混凝土交界面的传力或可能存在的间隙,并研究合适的保压水头和保压水温。1计算条件1.1结构计算模型三峡工程挡水坝为混凝土重力坝,电站厂房为坝后式。本文取一左岸岸坡坝段机组为研究对象,由于厂坝联接段压力钢管用垫层管取代了伸缩节,且垫层管下游端未设止推环,大坝的变形将通过垫层管影响到钢蜗壳的位移以及钢蜗壳与外围混凝土间的相互作用,因此计算模型中包括了大坝和垫层管,坝体包括钢管坝块及实体坝块,两坝块间为永久横缝。厂坝间分缝51m以下岩坡进行接缝灌浆。主厂房模拟至67m,上、下游副厂房只模拟下部实体部分,见图1、图2.垫层管长10m,其中坝内长5.8m,厂内长4.2m.厂内段垫层管可分为2段,即一期混凝土2.2m段和三期混凝土2.0m段。钢蜗壳进口位于垫层管下游端的下游侧1.1m处。图1大坝及厂房立体(含部分基础)图2厂房横剖面1.2材料参数(1)钢管、座环和固定导叶等钢材:E=210GPa,=0.30,=1.210-5/,=78.0kN/m3;蜗壳钢板厚度2064mm.(2)基岩:厂房基岩E=26GPa;大坝基岩E=1026GPa;=0.23,=0.8510-5/,导温系数a=0.083m2/d.(3)大坝混凝土:E=26GPa,=0.167,=0.8510-5/,a=0.083m2/d.(4)钢蜗壳与外围混凝土交界面:摩擦系数f=0.5.(5)厂房混凝土:=0.167,=0.8510-5/,a=0.083m2/d,=24.5kN/m3;弹模、绝热温升和徐变度分别见式(1)式(3),不计混凝土自身体积变形。E(t)=33.0t/(5.12+t)(单位:GPa)(1)Q(t)=24.2(1-e-0.837t0.849)(单位:)(2)C(t,)=C1()(1-e-0.3(t-)+C2()(1-e-0.005(t-)(单位:10-6/MPa)(3)式中:t混凝土龄期(d);加荷龄期(d);C1()=7.58+183.1/;C2()=12.4+35.3/。1.3主要边界温度曲线T(t)及表面放热系数(1)与大气接触的结构外表面T(t)=17.35+11.55sin(t-108.8)(),=15.1w/m2(4)式中:t时间(d),以1月1日为原点;=2/365,下同。厂房混凝土浇筑期间,厂房左右两侧边界也按式(4)考虑,以后按绝热边界考虑。(2)运行期上、下游副厂房内表面和主厂房67m表面T(t)=22.0+8.0sin(t-120.0)(),=5.13w/m2(5)(3)运行期引水压力管道和蜗壳内表面T(t)=17.53+8.5sin(t-130.0)(),=2326w/m2(6)(4)运行期大坝上游表面。运行期库水位以下的大坝上游表面采用库水水温边界条件,库水温度曲线随高程变化。1.4施工过程厂房混凝土浇筑层厚1.53.0m,分4区,对角两区同时浇筑,间歇期78d,15d浇一层。分别模拟了冬季和夏季浇外围混凝土情况,施工过程见表1.厂房混凝土入仓温度见表2.1.5保压水温与保压水头对冬季浇筑外围混凝土情况,在保压水头70m条件下,研究了保压水温控制措施,共计算了3种方案:(1)无保温措施:蜗壳内水体取当时河水,水体初温11.7,钢蜗壳外表面取气温边界条件(4);(2)采取保温措施:水体初温16,钢蜗壳外表面采用泡沫塑料保温,=2.0w/m2。(3)采取加温措施:水体温度低于22时就加温,计算时钢蜗壳内表面取恒温22,外表面按绝热考虑。在(1)、(2)方案中,保压水体同样划分单元参加计算,水的导温系数为0.01238m2/d.以上钢蜗壳外表面边界条件仅适用于尚未被混凝土覆盖的区域。表1厂房施工过程施工项目冬季浇外围混凝土情况(年.月)夏季浇外围混凝土情况(年.月)一期混凝土1999.102000.41999.92000.4二期(外围)混凝土2000.122001.22000.62000.9卸压2001.32000.9压力钢管合拢2001.42000.10三期混凝土2001.52001.82000.112001.2表2厂房混凝土入仓温度()区域11月3月4月10月5月9月强约束区自然2016弱约束区自然2018脱离约束区自然2020对夏季浇筑外围混凝土情况,蜗壳内水体取当时河水,水体初温25.5,在无保温措施情况下,计算了两种保压水头方案:70m和62m.2计算方法2.1温度场计算引水压力管道合拢前,考虑到坝体混凝土已强迫冷却到稳定温度厨行纵缝灌浆,因此先计算在边界气温作用下的坝体准稳定温度场,在此基础上,进行模拟厂房蜗壳保压浇混凝土施工过程的温度场仿真计算,厂房混凝土按自然冷却考虑。瞬态温度场的计算采用在空间上用有限元离散、在时间上用向后差分的隐式差分格式2。计算采用自行开发的温度场仿真计算程序3DUSTPCG.2.2应力与变形仿真计算在获得大坝与厂房的温度场后进行应力与变形分析,将大坝和厂房作为整体进行仿真计算,从厂房混凝土浇筑开始至运行期,模拟厂房混凝土保压浇筑过程和卸压、运行期加压过程以及由此引起的钢蜗壳与外围混凝土间的接触问题和水库蓄水过程等,考虑了自重、温度和徐变的影响。其中,采用初应变法2,3考虑徐变的影响,接触面单元物理方程见文献4。计算采用自行开发的可同时考虑混凝土温度徐变影响和缝面接触问题的结构仿真分析程序3DCRCPCG.2.3有限元方程组的快速解法计算范围包括坝体、厂房水下结构、钢管和部分岩基。为适应仿真计算的需要,厂房基本上按分层分区划分网格,各浇筑层分两层单元,整个计算模型共划分结点50927个,单元45632个。应力分析时总自由度最后将近14.5万个,不仅计算规模大、计算时段多,而且还要进行接触问题非线性迭代,计算工作量相当庞大。为此,在3DUSTPCG和3DCRCPCG程序中采用笔者提出的对称逐步超松驰预处理共轭梯度迭代法(SSOR-PCG)的改进迭代格式5作为求解器,与常用的大型有限元方程组的一维变带宽存储的三角形分解直接解法相比,在存储量和计算工作量方面都降低一个数量级以上,使得在微机上快速求解大型问题成为可能。方程组的规模越大,其效率更加显著。由于计算工作量少,计算舍入误差也小。3计算成果2003年7月水库蓄水至135m,同年10月机组投入运行。2009年汛后蓄水至175m,下游设计洪水位76.4m.仿真计算是从厂房混凝土浇筑开始至2020年止,几个主要特征时刻的计算成果见表3表4.3个特征点A、B、C的位置见图2.初期水位运行期的计算成果为2004年的计算结果,且未考虑水击力的影响;正常水位运行期的计算成果是2020年的计算结果,考虑了水击力的影响。计算中考虑了管内水重的影响。表中结点的传压(即钢蜗壳与外围混凝土间交界面的法向接触应力)或间隙是由接触面单元形心的值按单元面积绕节点加权平均获得,而平均传压或间隙是由蜗壳段所有接触面单元形心的值按单元面积加权平均获得。由这些成果表中可见,间隙或传压的分布是不均匀的。3.1冬季浇外围混凝土情况3.1.1无保温措施从表3可见,卸压前夕(2月份),平均传压0.12MPa.到初期水位运行期冬季(1月份,下同),平均传压只有0.10MPa,说明此时的内水压力由钢蜗壳承担。运行季节相同时,正常水位运行期比初期水位运行期的平均传压大0.450.49MPa,约占内水压力增量(0.62MPa)的76%.运行水位相同时,夏季运行期(7月份,下同)比冬季运行期的平均传压大0.300.34MPa,约为剩余水压的46%.可见,温度对传压的影响是显著的。利用夏冬两季3个特征点处钢板的温差及其外围3m厚混凝土内的平均温差以及平面应变问题的组合圆环轴对称温度接触应力公式,这种影响得到了验证。正常水位运行期夏季的平均传压为0.89MPa,超过了剩余水压。3.1.2采取保温措施采取保温措施后,平均传压虽有所减小,但幅度很小,仅0.02MPa,可见采取保温措施并不能有效地减小传压。从温度场来看,由于泡沫塑料保温主要对短期的温度骤变(如寒潮袭击)有效,而且上半圆蜗壳长时间处于低温空气中,浇筑到高程58m时,采取保温措施情况的上半圆蜗壳钢板温度只有7.47,比无保温措施情况只高1.37,时间越长,两者的温度越接近。由于上覆混凝土水泥水化热温升的影响,卸压时上半圆蜗壳钢板温度回升到16.4。3.1.3采取加温措施从表3可见,采取加温措施后,正常水位运行期夏季的平均传压为0.73MPa,比无保温措施情况减小0.16MPa,说明加温措施可比较有效地减小传压。初期水位运行期冬季尚存在间隙,平均间隙为0.53mm,接触部位主要在下半圆和上座环附近。表3冬季浇筑外围混凝土情况蜗壳与混凝土界面的间隙和传压时间无保温措施加温措施间隙/mm传压/MPa间隙/mm传压/MPaA点B点C点平均A点B点C点平均A点B点C点平均A点B点C点平均卸压前夕-0.110.100.160.12-0.010.040.110.07初期运行期冬季-0.050.190.160.100.440.76-0.53-0.040.06初期运行期夏季-0.300.640.610.44-0.130.510.210.27正常运行期冬季-0.520.820.630.59-0.310.680.330.43正常运行期夏季-0.761.091.130.89-0.530.950.840.733.2夏季浇外围混凝土情况3.2.1保压水头70m正常水位运行期夏季的平均传压为0.68MPa,比冬季浇筑情况且无保温措施情况(表3)小0.21MPa.初期水位运行期冬季,相当一部分区域尚存在间隙,平均间隙为0.66mm,比冬季浇外围混凝土但采取加温措施情况的0.53mm大24

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论