两相介质弹塑性地震反应分析方法研究-期刊.pdfVIP

两相介质弹塑性地震反应分析方法研究 on MethodtoElasto—PlasticSeismic of Study Response Analysis Porous Fluid．SaturatedMedia 李亮1杜修力1赵成刚2李立云1 (1．北京工业大学建筑工程学院．北京．100022z 2北京交通大学土木建筑工程学院，北京，100044) 摘要：针对流体饱和两相多孔介质的弹塑性波动方程组，运用伽辽金方法对该波动方程组进行 空间离散，得到两相介质弹塑性波动方程组的伽辽金弱式，然后应用中心差分法与Newmark常平 均加速度法相结合的时域积分方法，对上述波动方程组进行时间离散，构造求解两相介质弹塑性 波动方程组的显式时间积分歹4式，从而形成流伴饱和两相多孔介质弹塑性动力反应计算分析的时 域显式有限元方法。该方法采用了解耦技术，不需要求解联立方程组，具有节省计算机内存空间 和能够提高计算速度等优点。数值算例表明，该方法能够对两相介质在输^地震波作用下的弹塑 性动力反应进行较为准确的计算和分析。 关键词：流体饱和两相多孔介质弹塑性地震反应时域显式有限元方法 1 引言 流体饱和两相多孔介质在地震波作用下的弹塑性动力反应的计算分析是土动力学中的重 要研究课题之一。由于该问题的复杂性，描述该问题的方程为固相与液相耦联的二阶偏微分 方程组，因此在一艘隋况下无法得到解析解，而是必须采用数值方法进行求解。求解该问题 所采用的数值计算方法将显著地影响计算时所占用的计算机内存资源量以及计算所需要消耗 的机时，从而决定了计算求解的效率。现阶段求解两相多孔介质弹塑性动力反应问题的数值 方法以有限元方法为主，而现有的有限元方法基本都是隐式方法w或者隐-显式方法w。当这 些方法应用于时域问题的求解时，存在如下缺点：即在每一个计算的时间分步中都要求解一 个联立的线性方程组。当所求解的问题自由度很多时，其计算工作量将很大，所需要的计算 时间和计算机内存空间将极为可观，因此极大地限制了它们的应用和发展。为克服上述缺点， 赵成刚等人…应用中心差分法和Newmark常平均加速度法相结合的方法，提出了两相多孔介 质弹性动力反应计算分析的时域显式有限元方法。该方法的特点：采用了解耦技术，只需要 求单元刚度矩阵而不需要组集系统的总体刚度矩阵，不需要求解联立方程组，因而具有节省 计算机内存空间和能够提高计算速度等优点。通过将应用该方法得到的计算结果与相应的解 析解计算结果进行对比，从而验证了该方法的正确性。本文将把该方法推广应用于两相多孔 介质的弹塑性波动问题之中，形成两相多iL介质弹塑性动力反应计算分析的时域显式有限元 方法，这是对两相多孔介质时域有限元方法的丰富和发展。 2 两相多孔介质弹塑性波动方程组 针对平面问题，增量形式的流体饱和两相多孔介质的弹塑性波动方程组为’ ·车亮．巍津饱和两相多孔介质动力反应计算分析．北京交通大学博士学位敝，2004 两相介质弹塑性地震反应分析方法研究 103 白。警地^。)等十k等嘞警地：岷)等 +‰等叶叫降+等]+华降+等]} 一。降一掣]卸叫见学 mt， ’‰等讹^，，等心等岷学也悔，等‘ 一叫掣一下O(dUy)1Ph)岛学 m2) ‘ ‰{n降+等卜n)[警+等]) +吨掣一掣h等 口，， 日{《訾+等卜，陪+等]} 十叫掣一掣卜等