土石坝地震反应分析可视化研究.pdfVIP

蔓至旦全里圭垫查兰堂查垒望堡苎叁 盔鎏：!!!!：! 土石坝地震反应分析的可视化研究 』 L 宿晓辉李俊杰 (大连理工大学土木工程系．大连116024) 摘要土石坝在地震作用下的动力反应是一个连续而又复杂的过程．本文采用可视化技术，将土石 坝的地震反应过程通过屏幕适时反映，为全面进行坝体地震反应分析给出了一种新的途径和方法． I前言 inScientific 随着计算机及图形学技术的进步，科学计算可视化(、，isualization 用计算机图形学的理论和方法，将：[程计算所产生的大规模数据信息转化为直观的、易 于为理解的、用图形和图像形式表j示的静态的或动态的画面，便于研究人员加深对被模 拟对象的物理规律及其反应过程的认识，从而摆脱处理冗长的抽象数据的被动局面，丰 富了科学发现的方法【11。 可视化技术在用于固体、流体形态的连续变化方面，具有独特的优势，对土石坝 而言，当地震发生时，坝体的反应是一连续而复杂的过程，采用VISC技术可以清楚地 分析坝体的地震反应。近十几年来，国内外基于等价线性化方法的土石坝动力反应分析 已比较成熟，并可适用于各种土石坝坝型的计算。虽然此类程序可以得出坝体地震位移、 加速度和应力，但同时需要研究人员在程序中加八大量的比较和输出模块，还需整理大 量的输出数据，浪费时间与人力，并在程序运行过程中看不到每一瞬时的反应情况，缺 少一个整体、连续的概念。基于上述情况，本文结台大连理工大学研制的二维土石坝有 限元程序，编制了采用VISC技术的分析系统，初步解决了上述问题阱。 2土石坝动力有限元程序的几点说明 2．1土的动本构关系 土的动应力一应变关系是±工结构物动力反应分析的重要基础。土在周期荷载作用 下的应力一应变关系有两个特点，一个是非线性，另一个是滞后性。而土的骨干曲线反 映了动应变的非线性，滞回曲线反映了应变对应力的滞后性。 9(k：)。。(以)“5(t／m2) (1) G一==69 式中(t：)一为与材料的密度、粒径等有关的常数，口：为有效平均主应力。 土的动剪切模量、动阻尼比与动剪应变的关系由动三轴试验给出。 ——628— 2．2单元类型 土体单元采用四节点等参单元；在计算混凝土心墙土石坝、沥青混凝土心墙土石 坝和钢筋混凝土面板堆石坝时，心墙和面板采用两节点梁单元，梁单元转动方向的自由 度可以考虑弯矩的作用：对两种不同力学性质材料的接触面加入接触面单元，用来模拟 接触面的滑移或开裂；在进行库水与坝体动力耦合计算时，水体采用四边形八节点可压 缩等参单元。 2,3方程解法与迭代方法 土石坝体系动力有限元方程，采用逐步积分法求解，将输入地震过程分为几个时 段，在每一时段内，求出各单元的最大动反应，对非线性土体采用等价线性化迭代计算， 土体单元要满足 l警卜，f等卜％ ∽ 式中Gi为第i次迭代土体的动剪切模量：Dj为第i次迭代土体的动阻尼比。一般迭代扯5次 即可达到收敛要求。最后给出坝体各节点的位移、速度、加速度和各单元的最大应力值。 3土石坝地震反应可视化系统的组成及特点 根据土石坝动力分析有限元程序的流程，土石坝动力有限元可视化分析系统分为 如下三个模块： 3．1输入数据的前处理模块 依据输入数据，以交互方式返回输入控制数据，确定计算模型，绘出计算网格。 该模块采用人机交互方式输入，在帮助命令的提示下，提高了输入数据的准确性和准确 性，同时方便检查输入数据，以图形形式显示的有限元网格，可及时修改土石坝结构有 限元网格的错误，提高了效率。 3．2动力分析模块 在土石坝的动力计算中，坝体的位移、速度、加速度、应力指标是判别坝体是否 发