求解非线性动力方程边值问题的加权残数打靶法.docVIP

求解非线性动力方程边值问题的加权残数打靶法 85 工租实际中,由于动力微分方程的非线性性质,大多数只能使用数值方法米求解.样条加权 残数法是一种求解微分方程有效方法,它具有收敛快,稳定性好,便于计算机实施等优点[1]. 对于非线性动力问题,样条加权残数法有明显的优点[]. 1求解非线性动力问题的样条子域法 1.1非线性结构增量动力平衡方程 对于非线性结构,方程(1)是一个非线性微分方程,山于,均为时间f的函数,所以 方程(1)亦可写作 m+,D(f)+Is(f)=,(f)(3) 要求解微分方程(3),可采用增:肇法,即将时间离散为着千段,在每个时间间隔Af内, 把非线性问题当作线性问题处理,集合所有分段q的线性解便可得到非线性动力问题的近似 解.这是一种解决非线性问题常用的分段线性化方法. 方程(3)的增量动力平衡方程为 ,nA”()+c(t)△”(t)+(f)AH(f)=△,(f)(4) 式中 c(f): 0” k(t=+ Af(t)=,(f+At)一,() 方程(4)可作线性方程处理,是非线性动力平衡方程线性化的结果. 1.2三次B样条函数 在术文中,直接使用三次B样条砸数()作为时域的分段试函数..()的定义为 (f+2).r∈[一2,一1] I(f,J-2).一4(f+1).f∈[一1,0] 3(f)=1/6I(2一).一4(1一f)∈[o,1](5) I(2一r).f∈[1,2] ogt;2 1.3求解增量动力方程的样条子域法[2] 根据加权残数法的概念,在时间区间[ff]内设试函数”()为 “(f)=[a]{(f)} 式中 [a]=[a～laoa…a.…] 是待定系数矩阵. {(f)}:[(f一.)(f0):(r)….(r.).?.] △f=f+l—f. e6 =!= Af 是时间区间叭f…]上无量纲化的时间变量,有0≤≤】. 因为也()是用四个区间的曲线表达的分段多项式,在任一时间区间[f 只有四条曲线非零,这四条曲线方程为: 1/6(1一f)., 1/6[Z—r).一4(1一f).]l 116[(1+r).一4r.]f l/6r.j 对应于这四条曲线的待定系数分别为a,a,,a.+.,a. 所以,在Et.f…]内,试函数”(≠)为 )=a÷(1一r)3+a,-16-[(2一一4(1一]+a;+’i1(1+一川+ z ÷r. (),()可由上式求导分别得到. 在=.时,r有=0,所以在时刻有 1(a+a ?+a. 1 (1f)=刍(_aI_【+ai.t): ()=1(a i-l-2a.+a)J 根据增量方程(4),在[f?]内取 hu()=“(f)一”(.){ A”(f)=“(f)一”(f.)l △(f)=(f)一(f;)f A/(t)=,(f)一,(f;)j 把(7),(8),(9)式代入方程(4),形成残差(r)为 ){一寿m+c托)+÷r(11_『).1]R(/al- + {墨m_△3r2-厂4rc(fi)+_E(2(1_r)一4]k(t))a. + {一麦m+堑珏-3r.∞_)+÷[(1+r)r3-1) + 锰m+c(ff)+(fI))aI+_,(卅,(f 利用加权残数法的公式 f(r)(r)dr:oJD (6) (7) (8) (9) (10) 在子域法中,取权函数矿(t)=1,可导出样条函数——子域法系数的递推公式 去,(t)dt-_La l砉D--18--k(t,C(tk(tf一)l一 击一m+古k(t) L—1 m+ 6~tC(t)+R(t. (12) 当f:.时,=0,a一,a,a由初始条件定fl},l!lI把f=.时的初j!ff条件代入方程(8), 则可解出 a+”u..AtZ ～一— A — tz . 13) 口”.一”.i’ at=t~o+”:,,f+ 式中,.是把=.时的初始条件代入非线性动力平衡方程而得到. 在f:+l时刻有 u(t…)=一1(a r+姐…+am) 0I :(cI=j(a.十aIT”)}㈨) (t)=壶(a.一2a一十a一)I 利用初始条件求出初始待定系数a—,a,a后,可按递推公式(11)计淖出其余待定系 数,则由公式(14)计算任一时刻+,的位移,速度,加速度.但在非线性动力rg题中,因C (t.),k(t)是与位移及速度有关的变嚣,所以在每一步计算中,除了计算出待定系数ar下外, 还应计算出”,,把所求出的蜥,.代入c(),fl(t㈩)才能进行下一步计算,系数递 掖不能单独进行. 在文献(2)中,作者按文献[6]讨论了样条函数子域法计算公式的稳定性,所得结论是: 样条函数子域法是条件稳定的,其稳定条件是t ≤o.3～o.6(15) 式中,T.为在时间间隔[}J.+]内分段线性化方程的周期. 对于非线性变化大I地方.可按(15)式改变步长△f,得到较准确的动力响直.