任意面内荷载作用下薄圆盘的自由振动与屈曲分析.docxVIP

?

?

任意面内荷载作用下薄圆盘的自由振动与屈曲分析

?

?

李国荣周叮李雪红霍瑞丽

摘要：基于单集中力作用下半无限平面的应力分布公式，利用外载荷叠加原理，得到自平衡面内集中力系作用下薄圆盘的应力分布公式，通过积分计算进一步获得自平衡面内分布力系作用下薄圆盘的应力分布表达式。取切比雪夫多项式与边界函数的乘积作为容许函数，应用里兹法分别导出薄圆盘在任意面内静力荷载作用下的横向自由振动与屈曲的特征值方程，数值求解特征值方程得固有频率和屈曲荷载。与取幂级数和傅里叶级数乘积作为容许函数以及有限元结果对比验证了方法的快速收敛性和高精度。

关键词：横向振动;屈曲;薄圆盘;面内荷载;里兹法

：O326;TU311.3文献标志码：A：1004-4523（2021）05-1001-08

DOI：10.16385/j.cnki.issn.1004-4523.2021.05.014

引言

薄壁结构具有重量轻、经济性好等优点，被广泛应用于航空、桥梁、机械等领域，在使用过程中人们发现其较易发生共振。而结构边界约束等亦会引起薄壁结构的面内初应力，导致其原有的动力特性发生改变，使得实际使用中出现的共振难以控制。

圆盘在经典边界条件下的振动研究已十分成熟，一些学者研究了圆盘在各种边界条件下的振动特性。石先杰等[1]采用谱几何法分析了弹性边界条件下圆盘横向自由振动特性。武兰河等[2]采用微分容积法求解任意边界条件下中厚圆盘的轴对称自由振动。李秋红等[3]采用改进的Fourier?Bessel级数方法和Rayleigh?Ritz法对任意弹性边界条件下圆盘的自由振动进行分析。Shi等[4]提出了一种求解任意边界条件下圆盘自由振动的统一方法，将位移解用一种简单的三角级数展开形式表示。Zhang等[5]基于简化板理论和改进的二维傅里叶?里兹法，建立了圆盘在各种弹性边界条件下振动特性的统一分析模型。对于面内荷载作用下的振动特性也有学者进行了研究，杜国君[6]分析了均布载荷作用下圆盘在非线性弯曲静平衡构形附近的微幅自由振动，使用修正迭代法求解，得到了固有频率?载荷的特征关系。Maretic等[7]以贝塞尔函数形式给出了等压作用下圆盘横向振动的解析解。上述研究表明：初始荷载显著影响结构的动力学性能，而固有频率是反映结构力学特性的一个重要参数。

关于薄圆盘的屈曲问题，早在19世纪末，Bryan[8]就指出：沿周边受压圆盘的屈曲荷载与径向对称屈曲模态相对应。而对于复杂荷载作用下薄圆盘的屈曲问题，此时的挠度将不再是径向对称，求其解析解是十分困难的，这方面的研究文章还不多见。Aung等[9]研究了薄圆盘在中间和边缘径向载荷作用下的弹性屈曲问题。马连生等[10]根据圆盘轴对称特征值问题数学上的相似性，研究了不同理论间圆盘特征值的解析关系。邱文彪[11]利用圆盘的应变能构造圆盘屈曲问题的哈密顿体系，将圆盘临界载荷和对应屈曲模态的求解归结为广义辛本征值和辛本征解问题。丁华建等[12]利用有限元软件计算了圆盘在不同计算单元模式下的屈曲问题，结果表明，常用的板和壳单元对薄板屈曲问题求解有一定的局限性。Ayoub等[13]研究了圆盘在一对径向集中力作用下的自由振动与屈曲特性。张海军等[14]使用Ritz法求解了薄圆盘在径向面内荷载下的振动和屈曲特性。

本文基于弹性力学基本公式，引入参数h为力的作用线到圆心的距离，推导出圆盘在任意自平衡面内荷载作用下的应力场解析解。将应力场方程代入圆盘能量方程，得到任意面内荷载下振动与屈曲的特征方程，采用切比雪夫多项式和边界函数的乘积作为位移试函数，通过求解特征方程可以对任意面内荷载作用下圆盘的振动与屈曲特性进行分析。研究结果弥补了分析初应力下圆盘动力特性及稳定性时，面内载荷只能为径向加载的不足。

1任意自平衡面內荷载作用下薄圆盘的应力场

1.1自平衡集中力系作用下圆盘的应力场

考察半无限空间中如图1所示的薄圆盘区域，集中力P沿弦AB方向作用于A点，设其力的作用线AB至圆心O的距离为弦心距h，与r轴的夹角为。

由集中力作用下半无限平面应力公式可知[15]，在集中力P作用下，点M处沿r2方向的应力为

为维持与集中力P的平衡，沿圆盘边界需存在合力指向A点的分布载荷。当存在多个集中力时，每个集中力都会在圆盘边界引起对应的分布荷载。当圆盘边界的集中力自平衡时，则边界载荷径向均布，满足

把圆盘从半空间取出，将式（2）反向施加于圆盘边界即可满足边界处分布荷载为0。因此，自平衡集中荷载作用下圆盘内任一点的应力为式（1）与（2）的反向叠加。

过圆心O点做CD与力P作用线平行，将点M沿r2的应力转化为以点O为原点，r为正方向的极坐标应力：

1.2自平衡分布力系作用下圆盘的应力场

如图2所示，假设一个面内分布力函数表达式为，其中为分布力