施工期坝体温度场及温度应力全过程仿真的研究.pdfVIP

计算机科学2008Voi．35№．4B 施工期坝体温度场及温度应力全过程仿真研究 周华1 熊忠明2郝鹏2彭成佳2周艳国3 (武汉大学水资源与水电工程科学国家重点实验室 武汉430072)1 (贵阳勘测设计研究院 贵阳550002)2(武汉大学土建学院 武汉430072)s 摘要本文针对碾压混凝土坝体施工特点，研究了混凝土的热学参数及边界条件随时间的变化规律、混凝土浇筑过程 等，应用瞬态热传导三维有限元分析方法对重力坝施工期的温度场和徐变应力进行了动态仿真计算，分析了坝体在施工 期温度和应力的时空分布规律。本文的研究成果为坝体施工过程中温控措施、施工方案的优化提供了重要的依据。 关键词碾压混凝土，施工期，温度场，徐变应力，三雏仿真 式中：卢为表面放热系数(埘／(m2．h℃))；A为混凝土的导热 1 引言 系数(kJ／(rrL}L℃))；L为环境温度(℃)。 碾压混凝土是近几十年来发展起来的筑坝新技术，尤其 计算中初始条件和第一、第二类边界条件均可用温度边 是碾压混凝土重力坝的发展速度非常快，具有快速施工、节省 界条件来表示，第三类边界条件表示了固体与流体接触时的 水泥、降低造价、缩短工期等特点。在混凝土碾压施工期间， 传热条件，可用对流边界条件来表示。 由于碾压混凝土坝采用整体通仓铺筑，施工速度快，层面间歇 2．3有限元控制方程 期短，因此产生的水化热就积蓄在坝体内不易散发，混凝土内 根据热传导基本方程和边界条件，利用变分原理可导出 的最高温度并不比常态混凝土低多少，同时大量水化热的散 有限元控制方程为 发持续较长时间，利用层面散发效果不大，致使温度成为其主 [H]{T)+{F)一O (6) 要荷载之一，也会使坝体产生较大的温度应力，严重时还会产 式中：[H]为系数矩阵；{F)为边界温度荷载列阵。 生裂缝，因此控制温度应力是碾压混凝土筑坝技术中的重要 2．4徐变应力计算原理 问题‘11。 采用增量初应变方法计算坝体混凝土由于水化热、变温、 本文采用通用有限元分析软件ANSYS中的热分析模块外荷载引起的徐变应力变化规律。在计算时域内把时间离散 (THERMAI，)对西南地区某重力坝溢流坝段进行了全过程 成行个时段，在任一时段Atl内的增量有限元方程为： 仿真分析，计算了施工期坝体温度场和徐变应力的时空分布 K；△盈=△Pi (7) 规律，得到了一些有益的结论，为今后我国碾压混凝土坝的设 式中△Pi一△P5+△只+△Pi，△P{为自重、面力等外荷载引起 计和施工提供参考依据，具有一定的工程参考价值。 的等效结点力，△P；为变温引起的等效结点力，△B为徐变 2仿真分析基本原理 引起的等效结点力；Kr为整体刚度矩阵，△鑫为整体结点位移 增量。 2．1热传导方程 如果采用指数形式的徐变度[2]，则在计算徐变应变增量 按各向同性材料考虑时，非稳定温度场的热传导方程为 时具有递推公式。本文计算中采用8参数指数形式的徐变度 磊2n(孬十万十再)十一娶_n(磐+努+磐)+塑3r (1)(1) 公式如下： 式中：口为导温系数(mz／h)；T为温度(℃)；r为时间(h)；口为 P—D2“叫] (8) 混凝土的绝热温升(℃)。 式中At，B1，C1，D1，Az，Bz，c2，D2采用复形法进行拟合。 2．2初始条件及边界条件 (1)初始条件