水电站蜗壳保压浇混凝土结构的三维仿真分析2行政.docVIP

水电站蜗壳保压浇混凝土结构的三维仿真分析2行政论文范文大全 水电站蜗壳保压浇混凝土结构的三维仿真分析 2 计算方法 2.1 温度场计算 引水压力管道合拢前，考虑到坝体混凝土已强迫冷却到稳定温度厨行纵缝灌浆，因此先计算在边界气温作用下的坝体准稳定温度场，在此基础上，进行模拟厂房蜗壳保压浇混凝土施工过程的温度场仿真计算，厂房混凝土按自然冷却考虑。瞬态温度场的计算采用在空间上用有限元离散、在时间上用向后差分的隐式差分格式 。计算采用自行开发的温度场仿真计算程序3dustpcg. 2.2 应力与变形仿真计算 在获得大坝与厂房的温度场后进行应力与变形分析，将大坝和厂房作为整体进行仿真计算，从厂房混凝土浇筑开始至运行期，模拟厂房混凝土保压浇筑过程和卸压、运行期加压过程以及由此引起的钢蜗壳与外围混凝土间的接触问题和水库蓄水过程等，考虑了自重、温度和徐变的影响。其中，采用初应变法[2，3]考虑徐变的影响，接触面单元物理方程见文献 。计算采用自行开发的可同时考虑混凝土温度徐变影响和缝面接触问题的结构仿真分析程序3dcrcpcg. 2.3 有限元方程组的快速解法 计算范围包括坝体、厂房水下结构、钢管和部分岩基。为适应仿真计算的需要，厂房基本上按分层分区划分网格，各浇筑层分两层单元，整个计算模型共划分结点50927个，单元45632个。应力分析时总自由度最后将近14.5万个，不仅计算规模大、计算时段多，而且还要进行接触问题非线性迭代，计算工作量相当庞大。为此，在3dustpcg和3dcrcpcg程序中采用笔者提出的对称逐步超松驰预处理共轭梯度迭代法(ssor-pcg)的改进迭代格式 作为求解器，与常用的大型有限元方程组的一维变带宽存储的三角形分解直接解法相比，在存储量和计算工作量方面都降低一个数量级以上，使得在微机上快速求解大型问题成为可能。方程组的规模越大，其效率更加显著。由于计算工作量少，计算舍入误差也小。 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.45～0.49mpa，约占内水压力增量(0.62mpa)的76%.运行水位相同时，夏季运行期(7月份，下同)比冬季运行期的平均传压大0.30～0.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 传压/mpa a点 b点 c点 平均 a点 b点 c点 平均 a点 b点 c点 平均 a点 b点 c点 平均 卸压 前夕 - - - - 0.11 0.10 0.16 0.12 -