动力有限元方程显式逐步积分格式稳定性条件.pdfVIP

动力有限元方程显式逐步积分格式稳定性条件 景立平 (中国地震局工程力学研究所，黑龙江哈尔滨15008m ■翼t利用单元的最大频率是系统的特征频率的上界的概念，结合模态分析方法，对常用的动力宵限元方程的2 种显式逐步积分格式的稳定性进行了分析．给jnr便于实际应用的稳定性条件．并利用算例对给出的稳定性条件 进行了验证。 关t调t数值模拟；有限元；显式积分格式；稳定性 中圈分类号，P135．Ol；TB115 文■标讽码。A 文章螬号l CONDITIoNoFSTABILITYoFEXPLICITFINITEELEMENTSTEP ⅨTEGRALPATTERNFoRⅥ炯MoTl0N JING Li·ping AbstractiThe oftwo fiIlite are inthis useof血e stabilityexplicitst印integralpattemsanalyzedpapeLBy 血atthema】【imum of the of el锄emis 1irnitmecharacteristiconcof tlle f忙quency supper systcm，andcombining modal condj巧on0f isobtained．Atmesallledme，￡he∞ndid册is捃s拄矗。d a衄lysjs瑚et}10d，t1】es扭bilily by出e numerical calcul撕on． words：numerical Key simul撕on；finite elemem；explicitst印integralp舭m；s诅bility 式(1)在时间域上的逐步求解一般有隐式和显 1 引 言 式2种逐步积分途径。常用的隐式逐步积分格式有 有限元分析是近场波动数值模拟常用的方法。 法每一时步计算均需要求解线性方程组。当空间自 对于一个实际空间区域的波动数值模拟，利用有限 由度较少时可以用这类方法计算，而当自由度较多 元进行空间离散处理之后，根据虚位移原理可以得 时，如三维波动问题，将耗费很长时间，有些问题 甚至难以完成计算。因此，在实际波动数值分析中 到如下～般的动力有限元求解方程： 常采用显式逐步积分方法。显式积分方法不需要求 【M】{☆(f))十[c】li(f))+【蜘fⅣ(f))={，(f))(1) 解耦联的方程组，需要的计算机内存相对较小，但 式中：【_jlf]，[C】，【捌分别为求解体系的质量矩阵、一般是条件稳定的。此外，根据波传播的物理性质 阻尼矩阵和刚度矩阵；{n(t)1，{d(t)1，【“(t)1分别为知道波速是有限的，因此，波动具有时空局部特征。 单元节点处介质运动的加速度、速度和位移矢量； 也就是介质内某一点在某一时刻的运动只能影响邻 {F(f)}为作用于单元节点处的外力矢量。 近区域相邻时刻的运动。这一局部特征意昧着任一 2004一11—23 收■日期l2004一lo一21；●国日期I 薹土疆目，地震学联台摹金课题(1