粘弹性偏心梁的复模态解耦解耦.docxVIP

粘弹性偏心梁的复模态解耦解耦 当梁截面有两个对称轴时，截面的切割中心与形状中心重叠，梁的弯曲运动和旋转运动相互独立。然而对于许多工程应用中的梁结构,横截面只有一个或没有对称轴,剪切中心和形心不重合,弯曲运动和扭转运动之间存在耦合。在这种情况下,即使对于线弹性结构,也不能单纯使用弯曲运动和扭转运动的叠加来计算梁的动力响应。近些年许多学者研究了弯扭耦合梁的自由振动和阻尼作用下动力响应问题,但大多数研究仅考虑简单阻尼并使用实模态法进行求解,到目前为止还没有关于使用复模态法求解粘弹性偏心梁弯扭耦合问题的报道。 粘弹性材料在工业界有着广泛的用途,因此粘弹性系统的振动问题日益受到重视。粘弹性系统的自由振动和简谐受迫振动通常利用待定系数法进行分析,但该方法应用于任意激励的响应有一定的困难。模态分析法是处理振动问题的一种有效方法。实模态法的应用受到阻尼可正交化条件的局限,而复模态法可以求解一般性阻尼问题,因而复模态法可有效用于解决粘弹性系统的振动问题。复模态法在有阻尼离散系统和连续系统振动问题中应用已较广泛,但在连续系统弯扭耦合振动问题中的应用还较少。 本文使用Kelvin-Voigt模型表示材料的粘弹性性质,对横截面只有一个对称轴的单向偏心粘弹性梁建立运动方程,利用复模态的正交性进行解耦,求解结构在外部激励下的一般动力响应表达式。通过跟踪边界条件行列式零点的方法求解了结构复频率和复模态。通过算例分析了阻尼和弯扭耦合作用对复频率、复模态和外部激励的动力响应等方面的影响。 1 状态方程形式 图1所示薄壁梁结构截面只有一个对称轴,点s和点c分别为横截面的剪切中心和形心,两者之间的距离为xc。设轴向坐标y与梁的弹性轴(经过梁每个横截面的剪切中心)重合,x轴与对称轴重合。EIx和GJ分别是z方向侧向弯曲刚度和绕y轴的扭转刚度,μ是单位长度梁的质量,Is是单位长度梁绕y轴的极质量惯性矩。剪切中心在z方向的弯曲位移是v(y,t),绕y轴的扭转位移是ψ(y,t)。 考虑粘弹性阻尼影响,设应力应变满足Kelvin-Voigt关系 σ=Eε+Cε˙ετ=Gγ+Cγ˙γ(1)σ=Eε+Cεε˙τ=Gγ+Cγγ˙(1) 式中σ、τ分别为正应力和剪应力;ε、γ分别为正应变和剪应变;Cε、Cγ分别为与正应变和剪应变相关的阻尼系数。上标圆点表示对时间t的微分。 采用伯努利-欧拉梁模型,由Hamilton原理推导出粘弹性偏心梁单向弯扭耦合运动方程为 EΙzu?″+CεΙz˙u?″+μ¨u-μxc¨ψ=f(y,t)-GJψ″-CγJ˙ψ″+ΙS¨ψ-μxc¨u=Τ(y,t)(2) 式中上标撇号表示对坐标y的微分。 为化简运动方程方程,令 C1=CεIz?4/?y4,C2=CγJ?2/?y2,K1=EIz?4/?y4,K2=GJ?2/?y2(3) 为将弯扭耦合运动方程转化为状态方程形式,补充两个恒等式 μ˙u-μ˙u-μxc˙ψ+μxx˙ψ=0-μxc˙u+μxc˙u+ΙS˙ψ-Ιs˙ψ=0(4) 由此,状态方程形式的运动方程为 A˙z+Bz=F(5) 其中F={f0Τ0}Τz={u˙uψ˙ψ}Τ A=[C1μ0-μxcμ0-μxc00-μxc-C2ΙS-μxc0ΙS0]B=[Κ10000-μ0μxc00-Κ200μxc0-ΙS] 状态方程的A、B两矩阵中包含了K1、-K2、C1和-C2四个微分因子。在适当边界条件下(固支、简支、自由),利用分部积分可以证明其均为自伴随因子,A、B两个矩阵对称,具有如下性质 ∫L0zT1Az2dy=∫L0zT2Az1dy∫L0zT1Bz2dy=∫L0zT2Bz1dy(6) 设方程(5)的齐次解为z(y,t)=?(y)eλt,得到特征方程 λA?+B?=0(7) 求解特征方程(7)可得特征值λn和特征向量?n(y),由于A、B两矩阵元素均为实数,所以小阻尼情况下特征值λn和特征向量?n(y)为共轭成对复数,其表达式如下 λn=-δn+jωn(8) ?n(y)=?Rn(y)+j?In(y) (9) 其中j=√-1?δn、ωn为实数,δn是粘弹性梁的衰减系数,ωn为振动频率。 利用特征方程(7)和A、B两矩阵的对称性(6),可以证明得到特征向量关于A和B两矩阵的加权正交关系 ∫L0?ΤmA?ndy={anm=n0m≠n(10) ∫L0?ΤmB?ndy={bnm=n0m≠n(11) 其中由特征方程(7)可以得到an和bn关系为 λnan+bn=0(12) 利用状态空间特征向量对A和B两矩阵的加权正交关系,可对运动方程进行解耦,转化为常微分方程组。由于特征向量为复数,所以称为复模态法。把状态空间向量z对复模态进行展开 z(y,t)=∞∑n=1?n(y)qn(t) (13) 把式(13)代入式(5),所得方程等号两边同时乘以