宋志伟， 李 威， 渠鸿飞

(华中科技大学 船舶与海洋工程学院， 湖北 武汉 430074)

结构的动力稳定性一直受到人们广泛的关注[1～5]。解析法[1]和数值方法[2～4]已经被成功应用于结构动力稳定性分析中。本文提出了求解梁动力稳定性的拟小波方法[6]，利用该方法计算了两端简支和固支梁的动力失稳区，并讨论了周期性轴向力中恒定项对动力失稳区的影响，与解析解对比验证了采用拟小波法求解梁动力稳定性的可行性和有效性。

1 振动方程

当直梁受到周期性轴向力PS+PDcos(θt)作用时，根据欧拉梁理论，忽略转动惯量和剪切效应，其振动方程为[1]：

(1)

其中,E为杨氏模量,I为截面惯性矩,m为单位长度质量,PS+PDcos(θt)为周期性轴向力，PS为恒定项,PD和θ分别为幅值和激励圆频率，u为横向位移,t为时间。

为了简化，引入以下无量纲量

则,(1)式可以化为:

(2)

将(2)式改写为:

(3)

两端简支边界条件为:

U(0,τ)=U(1,τ)=0，

两端固支边界条件为:

U(0,τ)=U(1,τ)=0，

2 数值方法和离散式子

2.1 拟小波数值方法

拟小波方法是Wei[6]等人提出来的一种新型的数值方法，已经被成功应用于科学和工程多领域中[7～10]。关于该方法的数学理论和应用，可参考Wei等[8～10]的文献。根据Shannom定理,引进拟尺度函数[7]

δΔ,σ(x-xk)=

(4)

(5)

由于拟小波具有良好的局域特性,实际计算只需要在网格点x附近取2W个计算点即可达到计算精度。(5)式对空间坐标x的n阶导数为

(n=0,1,2, …)

(6)

在(1)式中需要求空间的二阶和四阶导数，其表达式见文献[8]。

2.2 离散式子

本文利用拟小波数值离散格式(6)式离散(3)式的空间导数, 四阶Runge - Kutta (RK4)法离散时间导数。具体如下[7]: 将空间X坐标进行均匀等分,单元网格大小记为ΔX=(1-0)/N(N为单元网格总数)。……