关于伽辽金法的一点注记1)

2022-04-28李映辉李翔宇

力学与实践 2022年2期

陈波李映辉李翔宇

(西南交通大学力学与航空航天学院,成都 610031)

伽辽金法类属于加权残值法,在《有限元方法》和《计算结构力学》教科书中都会讲授[1-2]。该方法在具体应用时会涉及试函数的选取,通常要求试函数为一组线性无关的函数,并满足结构的所有位移和力边界条件。一些文献指出在难于选取满足所有边界条件的试函数时,可以选取仅满足位移边界条件不必满足力边界条件的试函数,通过增加试函数的个数弥补不满足力边界条件引起的误差。然而,这种简化与加权残值法的基本原理相违背,获得结果可能无法收敛于正确解。下面以计算同时受轴向力和横向分布力悬臂梁的横向位移和固有频率为例进行说明。

1 轴向力下悬臂梁横向位移和固有频率

图1(a) 所示为端部受轴向力P作用的均匀悬臂梁,长度L,抗弯刚度EI,单位长质量密度ρA,横向分布载荷q(x,t),横向位移w(x,t)；其微段dx受力如图1(b) 示,图中M为弯矩,V为剪力,θ为截面转角。梁微幅振动方程为

图1

边界条件为

下面用伽辽金法分别求解悬臂梁横向位移和固有频率。

1.1 横向位移

式(1) 中去掉时间导数项,假设外激励与时间t无关,得横向弯曲的微分方程为

式(2)中用ws(x)和qs(x)分别表示梁横向位移和静态横向分布载荷。假设

其中,aj为待定系数

为满足位移边界的试函数,式(4) 中βj是超越方程cos(βjL)cosh(βjL)+1 = 0 的根。式(3) 代入到式(2),两边同乘ϕi(x) (i= 1,2,··· ,n),并在[0,L] 积分得到矩阵方程为

其中

由式(5) 可求得

将式(7) 代入式(3) 得悬臂梁横向位移。

1.2 固有频率

略去式(1) 中外激励项q(x,t),并设w(x,t) =W(x)sin(ωt+ϕ) (W(x) 为振型函数,ω为系统的固有频率),得到

取仅满足位移边界的试函数

其中,bj为待定系数,ϕj(x) 见式(4)。参照上述横向位移求解过程,得矩阵齐次方程为

其中

为使式(10) 有非平凡解,则令其系数矩阵的行列式为零,即

由式(12) 可求出悬臂梁固有频率。

2 验证和讨论

考虑物理参数：L=1 m,EI=4.0 MNm2,ρA=100 kg/m[3]。不失一般性,引入无量纲量

其中,Pcr=EIπ2/(4L2) 为悬臂梁临界屈曲载荷[4],为了使梁不发生屈曲失稳,应保证1。

为检验伽辽金法的正确性和收敛性,图2(a) 和图2(b) 分别给出了=0 时,试函数数目对悬臂梁端部横向位移和前三阶固有频率的影响。注意,悬臂梁(= 0) 在均布载荷作用下横向位移及其固有频率的精确解已在材料力学和振动力学教材给出。可见随着试函数数目增加,悬臂梁端部横向位移和固有频率的伽辽金解迅速收敛于精确解。值得指出的是,由于=0,此时伽辽金的试函数(式(4)) 同时满足悬臂梁的位移和力边界条件,因此随着试函数数目的增加,伽辽金解收敛于精确解。为减小试函数数目的影响,下面计算皆取n=10。

图2 悬臂梁横向位移的解析解与不同试函数个数下的伽辽金解比较

下面分析试函数(式(4)) 仅满足位移边界条件时(右端部不满足剪力边界条件),伽辽金法求解轴力作用下悬臂梁的横向位移和固有频率的正确性。图3给出了不同轴向压力下悬臂梁横向位移的精确解和伽辽金解,其中悬臂梁横向位移的精确解采用格林函数法求得,具体计算见文献[4]。可见随着轴向压力的增加,横向位移的精确解逐渐增大,该结果是预料之中的,因为轴向压力使梁软化增加其位移；然而横向位移的伽辽金解竟呈现截然相反的趋势,它表明轴向压力使梁硬化减小其位移,显然这个结论是错误的。进一步,表1 给出了不同轴压下悬臂梁前三阶固有频率。从表中可见,随着轴向压力增大,前三阶固有频率的精确解减小,一阶固有频率变化最为剧烈,而一阶固有频率的伽辽金解增加,二、三阶固有频率仅微弱减小,再次得到错误的结果。