考虑畸变的起重机薄壁箱梁的自由振动分析

2022-11-03谭敏尧程文明李杭飞臧付连

西南交通大学学报 2022年5期

谭敏尧，程文明，李杭飞，臧付连

（西南交通大学机械工程学院，四川成都 610031）

在起重机械中薄壁箱形梁在工程中得到了广泛的使用，在这些应用中减轻重量是非常重要的.同时，由于这种结构的特点是高刚度重量比，且在静载荷和动载荷的作用下薄壁箱形梁对扭转和不稳定性都非常敏感. 薄壁箱形梁在反对称荷载作用下，除引起截面扭转外，还产生截面畸变. 在薄壁箱形结构设计中，传统方法仅考虑延伸、弯曲、扭转、翘曲4 种位移形变，在高阶模态下，若不考虑畸变效应，将大大影响自由振动频率预测的准确度. 为了提高薄壁箱形梁结构的计算精度，根据结构力学特征建立广义位移模型，在原有基础上，增加畸变形变，这对薄壁箱形梁振动特性的研究具有重要的意义.

Gere 等[1]采用Rayleigh-Ritz 方法研究了各种边界条件下均布梁的弯扭耦合自由振动. 根据圣维南理论，Dokumaci 等[2-3]在静载荷和动载荷作用下建立了一些存有扭转的模型. Ohga 等[4]利用传递矩阵法研究了开、闭截面薄壁构件的自由振动. 此外，Adam[5]分析了单对称梁的耦合弯扭振动. Arpaci等[6-7]研究了任意截面梁的自由振动. 对于任意开口截面薄壁梁的振动，Jrad 等[8]推导了不同边界条件下弯扭耦合模态的高自由振动模态的解析解.Gökdağ等[9]研究了单对称开截面梁的自由振动和受迫振动的弯扭耦合. 在模型中，由虚功原理导出了包含质量和刚度特性的动力刚度矩阵. Rafezy 等[10]建立了三维非对称截面的复合材料梁弯曲运动的精确动力刚度矩阵. 李海超等[11]提出一种半解析法来分析圆柱壳结构自由振动特性. 蒲育等[12]基于扩展的N阶广义剪切变形梁理论获得弹性地基简支梁的自由振动精确解. Ovesy 等[13]通过考虑不同的本构方程假设，对薄壁组合梁的自由振动特性进行了研究. Vo 等[14]利用精剪变形理论对任意叠合梁进行了振动和屈曲分析. 该理论解释了剪切应变随梁高的抛物线变化，产生的耦合称为三耦合振动和屈曲. Filippi 等[15]提出了一类用于组合结构静力和动力分析的新型精细化梁理论. 这些梁模型是通过对梁截面上位移场的3 个分量执行切比雪夫多项式的高阶展开得到的. 考虑轴向翘曲效应和材料各向异性对壳壁中表面剪切应变的影响，Wang 等[16]提出了一种精度较高的薄壁封闭截面组合梁结构动力学模型.

虽然一些文献对薄壁闭合梁进行了动力分析，但很少有文献考虑弯曲、扭转和畸变的耦合来求解薄壁闭合梁. 分析发现，薄壁箱梁的截面畸变对固有频率有显著影响. 因此，要正确地预测薄壁箱梁的动力特性，必须同时考虑弯曲、扭转和畸变形变.Kim 等[17]就薄壁箱梁静力和动力进行了翘曲、扭转和畸变耦合变形的静力和动力特性分析. Kaya 等[18]采用微分变换方法（DTM）分析轴向加载的闭口截面复合材料Timoshenko 梁的自由振动响应，该梁的材料耦合特性是由于铺层方向的影响而发生弯曲和扭转振动. Kim[19]对在轴向力作用下组合梁的空间耦合进行了自由振动分析，利用能量原理推导出运动方程和力-位移关系方程，并基于位移分量的幂级数展开式给出位移参数的显式表达式. 基于能量法，Prokić等[20]利用Benscoter 理论推导了薄壁梁的微分运动方程，研究了简支薄壁梁的自由振动行为. 但仅考虑到横向剪切变形、转动惯量及翘曲影响，假设了薄壁箱形梁板面完全刚性，即忽略了箱梁横截面的畸变形变. Dancila 等[21]提出了除非考虑薄壁箱形梁横截面畸变形变，否则无法准确预测更高的固有频率.

本文基于弹性力学理论，用能量法分析了简支薄壁箱形梁的自由振动. 考虑到弯曲、扭转和畸变的耦合作用，预测了薄壁箱形梁自由振动的固有频率和模态振型. 根据Benscoter 理论，在广义位移的基础上，将各变形模态用广义坐标位移表示出来，进而求出薄壁箱形梁的固有频率和模态振型的精确显示公式. 该表达式简单明了，适合于一般的简支薄壁箱形梁理论的实际使用.

1 动力学基本理论

1.1 基本假设

为了确定薄壁箱形梁的几何形状，使用了两个坐标系. 如图1 所示. 图中：θ为扭转角. (x，y，z)为笛卡尔直角坐标系，其中，x轴与薄壁箱形梁的梁长重合，y轴和z轴与横截面的主轴重合. (n，s，x)为局部坐标系，其中，n为沿等值线（截面中线）测量的法向坐标，s为任意起始点沿等值线测量的剖面坐标. 点M(w，v，u)为横截面上任意一点，其中，v、w、u分别为y、z、x轴方向的位移. 薄壁箱形梁的理论基于以下假设：

图1 薄壁箱形梁坐标系Fig. 1 Coordinate systems of thin-walled box beams

1）在各板中，中线平面沿轮廓线方向的线剪切应变分布与圣维南扭转的线剪切应变分布相同；

2）假设薄壁箱形梁截面的中面轮廓不可延伸；

3）变形量足够小，可采用线弹性理论.

1.2 运动学方程

根据薄壁箱形梁的静力屈曲分析[22]可以得出，薄壁箱形梁上以任意一点M为初始点的位移分量在小扭转的情况下，使用Vlasov 模型产生的经典位移关系为

式中：χ为畸变角；i= 1 ～ 4 是箱梁各板的编号；ωi为扭转翘曲函数，如式（4）；Ωi为畸变翘曲函数，如式（5）；ηi、ψni分别为横截面y轴和z轴方向的畸变位移，如式（6）、（8）；t为时间；Zi(s)、Yi(s)、hi(s)均为横截面系数，分别如式（9）～（11）； ξ1、 ξ2分别如式（12）、（13）.

式（4）～（13）中：b和h分别为箱梁横截面的宽和高；(z0,y0)为横截面的剪切中心；Lside为箱梁横截面各边相对应的长度.

根据Kim 等[17]理论研究，横截面畸变形变的几何系数ψsi(s)可满足主要的静、动平衡条件. 此外，当每个板是不等厚时，κ是不连续的.

本文重点研究薄壁箱形梁在动载荷作用下的行为. 因此，在自由振动和线性动力学框架下，用式（1）～（3）作为运动学方程是充分、有效的. 为此需要对应变、动能和外载荷进行分析，进而推导运动控制方程. 应变能U的微分为

式中：σxi和τsi为薄壁箱形梁的正应力和切应力分量；εxi和 γsi分别为正应变和剪切应变，分别如式（15）、（16）；Fi为各板的面积；L为梁长.

式中：u0、w0、v0分别为任意点M在u、w、v方向的初始位移.

应变能的变化可以用薄壁单元截面上的内力函数表示. 在动力学系统变分原理的基础上，推导动力学运动方程. 根据式（14）对横截面进行积分后，应变能可表示为

式中：N为轴向力；My和Mz分别为关于y轴和z轴的弯矩；Bωθ为双力矩；Bωχ为畸变引起的双力矩；Msv为圣维南扭矩；Msvχ为畸变引起的扭矩. 具体如式（18）所示.

2 运动控制方程

在能量原理的基础上，应用哈密顿理论推导运动控制方程. 箱梁总能量的微分可由应变能U、动能T和外载荷功W表示，如式（19）.

式中：t1和t2分别为起、止时间.

利用虚位移原理可以得到薄壁箱形梁的运动方程. 图2 中考虑了横截面x1=x和x2=x+dx之间的一个小单元，在中面上承受单位面积上的外载荷p(px,py,pz)，其中：px、py、pz分别为x、y、z单位面积上的外载荷. 图2 中：ρ为单位体积质量密度. 在横截面上作意一点的应力为

图2 薄壁箱形梁的微单元Fig. 2 A micro-element of a thin-walled box beam

式中：α为轮廓中线的法线与x轴之间的夹角；τxs为切应力.

虚位移δu满足必要的连续性条件和位移边界条件，可以采用与实位移向量相同的形式. 根据虚位移原理，虚位移可表示为

加速度矢量可以表示为

等质量密度薄壁箱形梁单元的动能为

动能的微分表示为

根据能量原理，薄壁箱形梁单元单位长度外载荷和惯性力的虚功为

经整理计算，应用虚位移原理，可获得关于虚位移参数δu0、δv0、δw0、δθ、δχ的方程，对于任意虚位移参数均能恒等地满足这些方程，即虚位移参数的系数表达式必须消失. 将式（17）代入式（18），可以得到控制方程和相关的边界条件，则关于延伸、弯曲、扭转和畸变耦合的方程为

式中：mz和my分别为单位长度外力引起的关于z轴和y轴的弯矩；mp和mω分别为单位长度外力引起扭矩和双力矩；mpχ和mΩ分别为单位长度外力引起畸变扭矩和畸变双力矩.

式（26）～（30）是薄壁箱形梁延伸、弯曲、扭转和畸变的自由振动的最一般形式，变量u、v、w、θ、χ是完全耦合的.

对于具有连续质量分布且结构形式简单的结构可采用广义位移法. 对于简支梁，设任意时间t的任一截面的位移可用三角函数表示为

式中：V(x)、W(x)、Θ(x)、X(x)分别为横向位移、纵向位移、扭转和畸变的振幅，如式（31）.

式中： s in λk为满足位移边界条件的振动模态函数，其中， λk=kπx/L，为形状函数,k=1，2，3，…；CV、CW、CΘ、CX均为常数.

在简支梁的情况下，各支点处的支撑为防止旋转和可自由弯曲，其边界条件为

为了求出CV、CW、CΘ、CX为非零解，式（30）～（32）的导数必须为0，即

式中：ω*=，E为弹性模量；G为剪切模量；Iy和Iz分别为关于y、z轴的惯性矩；It和Iωθ分别为圣维南扭转和翘曲常数；Id和Iωχ分别为畸变引起的扭转和翘曲常数；Ip为扭转和翘曲耦合惯性矩；Ipχ为畸变与翘曲耦合惯性矩；A为横截面面积；A1、A2分别为纵向、横向弯曲与畸变耦合惯性矩参数；A3、A4分别为纵向、横向弯曲与翘曲耦合惯性矩参数.

经整理可得四阶频率代数方程为

式中：a、c、d、e和f均为常数.

3 算例分析

3.1 畸变模态下薄壁箱形梁的自由振动（算例1）

为了验证理论公式（式（34））的正确性，求解薄壁箱形梁在畸变模态下自由振动时的频率和振型，并与有限元分析（ANSYS）结果进行对比. 图3 为简支薄壁箱形梁的几何参数. 通过对薄壁箱梁的模态分析，验证该方法的有效性和准确性. 该梁的材料参数和几何参数为：E= 210 GPa，G= 80.7 GPa，ρ=7 850 kg/m3，b= 150 mm，h= 200 mm. 薄壁箱形梁的板厚为：m1=m3= 4 mm，m2= 8 mm，m4= 5 mm. 在有限元分析过程中，ANSYS 软件中采用4node181的壳单元进行分析，在网格划分后生成1 932 个节点和1 890 个单元.

图3 薄壁箱形梁的几何参数Fig. 3 Geometrical parameters of thin-walled box beams

表1 为考虑横截面畸变形变的理论方法和有限元分析方法下薄壁箱形梁其中的4 阶模态（2、4、6、8 阶）的固有频率对比. 为了对结果进行评估，将有限元法的固有频率表示为ωFEA，本方法的固有频率表示为ω1,Δ为本理论的固有频率与有限元固有频率的相对误差.

由表1 可以看出：固有频率结果具有较强的一致性，其中最大误差为0.89%（在第8 阶时），验证了本文分析方法可以较准确地进行薄壁箱形梁固有频率的预测；伴随着模态阶数的升高，固有频率的误差逐渐增加.

表1 薄壁箱形梁其中4 阶的固有频率比较Tab. 1 Comparison of 2-, 4-, 6-, and 8-order natural frequencies of thin-walled box beams

图4 为2、4、6 阶和8 阶的弯曲模态振型.

图4 弯曲模态振型Fig. 4 Bending mode shapes

由图4 可以看出：薄壁箱形梁的所有振型均为傅立叶变换，所得到的弯曲模态振型与ANSYS 弯曲模态振型结果吻合良好.

3.2 不同长度对自由振动的影响（算例2）

为了研究畸变形变对薄壁箱梁自由振动特性的影响，考虑了3 种不同长度的结构. 将本文的计算结果与Prokić等[20]的弯扭耦合（不考虑畸变形变）情况下求解的自由振动的精确解进行对比分析. 薄壁箱形梁的材料参数和几何参数与算例1 相同. 表2 为箱梁长度分别为3、4、5 m 时薄壁箱形梁的前4 阶自由振动频率及其相对误差. Prokić方法计算的弯扭组合结果用ω1表示；考虑畸变形变的相关结果用ω2表示；Δ1为Prokić 结果与有限元分析法的自由振动频率的相对误差；Δ2为考虑畸变形变情况下的自由振动频率与有限元分析法的自由振动频率的相对误差.

表2 简支薄壁箱形梁的固有频率Tab. 2 Natural frequencies of simply supported thin-walled box beams

在表2 中，将ω1与ωFEA进行比较，可以看出：本方法的解与有限元方法的结果有较强的一致性，且准确性和有效性较高，尤其是在第4 阶模态情况下，最大相对误差Δ2仅为0.38%（当L=3 m）；与Prokić 方法相比，自由振动频率ω1＜ω2，即考虑畸变形变可有效降低自由振动频率，进而提高了自由振动数值精确解的有效值；长细比越大，自振频率越小；在高阶模态下考虑畸变时的固有频率更接近于有限元法的解，因此，在高阶模态上考虑畸变效应的振动分析更占优势.