一种广义精细积分法.pdfVIP

第39卷第5期 力 学 学 报 VoL39．No．5 ChineseJournalofTheoreticaland Mechanics Sep．，2007 2007年9月 Applied ·种广义精细积分法D 富明慧2)刘祚秋林敬华 (中山大学应用力学与工程系，广东省近岸海洋工程重点实验室，广州510275) 摘要提出了求解非齐次动力方程特解的一种精细数值积分法，该方法与通解精细积分法具有相同精度．首先 选取一个积分形式的非齐次方程特解，将积分区域划分为2“份，并对之进行精细的数值积分；然后针对载茼 为多项式，指数函数及三角函数的情况，将积分求和转化为一个递推过程，按此只需n次矩阵乘法就能计算 出积分和，从而得到非齐次方程的特解．该方法的优点是能与通解的精细积分过程有机地结合起来，具有极高 的精度和效率。同时还具有较广泛的适用范围．算侧结果证明了该方法的有救性． 关键词结构动力分析。精细积分，指数矩阵运算，数值积分，递推算法 中国分类号：0321文献标识码： A 文章编号z 0459-1879(2007)05-0672—06 引 言 同时还具有高精度(与通解精细积分的精度相同)和 高效率(归纳出一个递推算法)，以及较广的适用范 精细积分法自1994年由钟万勰ll硎提出后，在 围．由于本文方法将通解和特解这两个精细积分过 结构动力分析中获得了日益广泛的应用．精细积分 程有机地结合起来，是精细积分法的进一步发展， 法在求解动力方程的齐次方程时具有突出的优势， 故称之为广义精细积分法(GPTSIM)． 一方面它既有极高的精度，可认为是计算机上的精 确解，另一方面该求解过程稳定性好，因此具有较 1通解的精细积分 高的效率．文献[3】的工作使精细积分法的性能进一 结构动力学方程为 步得到提高．但对于非齐次方程的特解，由于经常会 涉及矩阵求逆问题14,51，这种矩阵求逆的运算不仅 MX+c支+KX=F(t) 计算量大，有时甚至是不稳定的或者不可能的，这 也在一定程度上限制了精细积分法的应用．针对这 其中M，c和jr分别为质量阵、阻尼阵和刚度阵， 一问题，文献【6—8】提出了增维精细积分的方法，该 x为位移向量，F为载荷向量．方程(1)的初始条 法的特点是先将非齐次方程转化为齐次方程，然后 件为 再进行精细积分，从而避免矩阵求逆，但这种转化 (2) 是有条件的．文献[9，10】提出了一种利用数值积分 x(o)=xo，x(0)=Xo 求解非齐次方程特解的方法，这种方法具有简便易 其中xo，支。分别是初位移向量和初速度向量．方 行的特点，其精度会随选取积分点个数的增加而逐 程(1)可改写为 渐提高，但与此同时计算量也随之急剧增加，实际 上特解要想达到或者接近通解的精度几乎是不可能 o=Hv+r(t1 的． 综上所述，到目前为止尚缺乏一种普遍适用、且