matlab 加权残值法,基于加权残值法的高阶直接积分算法

第s期第2蜷 浙 ：}工 大 学 学 报 (自燕科学版) m VO1

．24

1 990年9月 Journal of：Zhejiang Unive~sity (Natural ScienCe) Septeml r|1990

基于加权残值法的高阶直接积分算法

段继伟 蔡承文

(力 学 系 )

提 要

本文以三次的Hermite插值多项式逼近位移响应，应用加权残值法推导了一个直接积 分递

推算{去。在分析算法性能的基础上，推荐了一族无条件稳定的三阶算法 通过算铘和几种类似算

珐作了比较。

关键词：结构动力学，振动，动响监，数值积分法

＼

l 算法的推导

设结构被离散化为N自由度线性系统，运动方程为

．

MU +CU+KU=P(t) (1)

式中M，c、K为质量，阻尼，刚度矩阵，U P为位移、 载荷列阵。崂 时间步硖，在时间间

隔[t·， -]中，用三次的Hermit,插值多项式来近似位移和载荷变化：

U(f)=U。矸o( )+hUoHo(~)+U。+1Hl(f)+hU叶 IH I(f) (2)

P(f)=P。H。(f)+hP H 0(f)+Pl+1H I(f)+hP" l H1(f) (3)

一

H口 卜 。 2 ， 一2 。 (d)

／4I=3f 一2f。， ／41= 一f +f。

式中 f是无量纲时间}f=t +fh。近似式(2)、(3)代入运动方程(1)一般不能 满 足， 记

其残值为R(f)。取两个权函数 -(f)，w：(f)，命残值的加权积分为零：

． ． 。

．

r1

- ljoⅣ，( )R《f)df 0 ( 1，2) (5)

-

， 引进参数 ．，B， ．(t=1， 2)；

m ，

J )州

本文于1989年8月收

口

(

、●-rl、，--●●，

。

／ f

"

．

“ “

●●，J r●●●J

维普资讯 http://www.cqvip.com

幕5南 ， 段缝伟等：基于加肛残值法的高阶直接积分算法 745

因为在两个加权方程(5)中只涉及上述积分，所以可用 6个参数 ，p 的某种组合来代

表权函数，而不必计较权函数的具体形式。文献[1—4]的方法都可作为本文算法 的 特殊 情

况。整理(5)式得递推格式如下

BI L： ： ⋯C⋯,s，1 州s I：=+ 啦“】

(7)

B。1=(6—12a )M+(6a。一6pf)hc+(3pi一21' )h。K I

(8)

Bf z=(一2+6a )M+(一2口j+3pi)hC+(一pi+ j)h K j

c 1=(6—12a。)M+(6a 一6pf)hc+(一1 4-3pi一2 )h K 1

， (9)

c．2=(4—6a。)M+(一l+4口i一3pf) c+(一口 +2p — f)h K J

f(3p1—2 I)I(1—3BI+2 1)ll

S=I l，

(3月z一．2 z)I(1—3月2+2 t)I J

f(一pj+ 1)I( 1—2p1+ I)I1 ‘

T=1 I (10)

l(一 ≈+ i)I( z一 B2+ ≈)I l

单步的递推方程(7)的系数矩阵是2N×2N阶的。但在经典阻尼(满足cM’ K=KM’-c)的

假定下，它可以解耦为两个N×N阶的矩阵方程

AoU口+I=DI1U +Dlt(hU-)+Ell( P +1)+E1 z(h P．)

+Fl1(h。P．+1)+Fl 2(h。P ) ．

A0(hUt+1)=D2IU +D2：(Ⅲ