机械模态分析作业

机械模态分析作业:如图1所示是一个单自由系统附件一个减振器形成的的两自由振动系统,已知m1=105kg，m2=7kg，k1=10000N/m，k2=410N/m，c2=1.15N·m-1·s，F1(t)=F1ejωt。求：（简化为粘性比例阻尼进行实模态分析）1．物理坐标下的振动微分方程；2．频响函数矩阵；3．频响函数的模态展式矩阵；4．脉冲相应函数；5．画出H11（ω）的幅频特性曲线，相频特性曲线，实频特性曲线，虚频特性曲线，Nyquist图，Bode图；6．固有频率，阻尼固有频率；7．画出振型图；8．模态坐标系下的振动微分方程；9．模态参数：复模态质量，复模态刚度，复模态阻尼。10.按实模态系统，给出灵敏度分析。11.集全班同学的数据（必要的话再补做不同m2，k2，c2参数下的数据，画出x1的最大振幅与m2，k2，c2，的变化曲线，从而分析出减振器的最佳参数。解：1．振动微分方程对质量m1、m2绘分离体图（如图1-1），用牛二定律列分离体在铅垂方向的力平衡方程得1221221111122122122()()()()Fcxxkxxkxmxcxxkxxmx（1.1）将（1.1）整理可得：112211221122222222000mxccxkkkxFmcckkxxx（1.2）且m1=105、m2=7、k1=10000、k2=410、c2=1.15，代入（1.2）得：0410410-410-104101.151.15-1.15-15.17001051212121Fxxxxxx（1.3）可以得出此二自由度系统振动微分方程为：()MxCxKxft其中M=700105；C=1.151.15-1.15-15.1；K=410410-410-10410；f(t)=01F图1-1、系统的分离体图2．频响函数矩阵由书P25（1.4-58）公式可知，此二自由度系统频响函数矩阵为一2×2方阵，其表达式为：21()()HKMjC，其中M=700105；C=1.151.15-1.15-15.1；K=410410-410-10410；图1两自由度振动系统（2.1）写成矩阵形式：（2.2）3．频响函数的模态展式矩阵1）求解瑞利阻尼矩阵由于粘性阻尼矩阵C无法进行正交性对角化，故不能直接应用坐标变换将（1.3）解耦。由于在该题中，粘性阻尼相对很小，对于小阻尼振动系统，可以利用瑞利比例阻尼来代替粘性阻尼，以获得可对角化的阻尼矩阵。（1）瑞利比例阻尼系数的确定瑞利比例阻尼：C，其中410410-410-10410；a、为瑞利比例阻尼系数瑞利比例阻尼系数存在以下关系：1112222222，其中i为圆频率iif（if为系统固有频率，书中表示为i）；i为阻尼比iii将上式写为矩阵形式：111222122122可得：1111222122122，其中、iif，iii（3.1）由此可知，只要我们确定了一个系统任意两阶的固有频率及其阻尼比，就可以确定出瑞利比例阻尼系数，从而得到瑞利比例阻尼矩阵。（2）求该二阶系统的一、二阶固有频率及其阻尼比利用求解该系统振动微分方程()MxCxKxft的特征值i来确定固有频率及其阻尼比。由书P23（1.4-43）-（1.4-46）公式为求解步骤，下面利用Matlab来计算固有频率i和阻尼比i：编写Matlab程序polynomial.m求特征方程，程序如下：symsx;m1=105;m2=7;k1=10000;k2=410;c2=1.15;M=[m10;0m2];C=[c2-c2;-c2c2];K=[k1+k2-k2;-k2k2];y=det(M*x^2+C*x+K)解以上求得的多项式：p=[735644115920115004100000];x0=roots(p)由特征值可得：2650.72593.72872.02211、0396.02593.72872.0112882.102626.107253.022、0707.02626.107253.0222（3）求瑞利比例阻尼系数及瑞利比例阻尼矩阵根据公式（3.1）编写Matlab程序rayleigh.m求解特征方程，程序如下：functionCr=rayleigh()%--计算瑞利阻尼系数alpha和beta--xi1=0.0396;xi2=0.0707;f1=7.2650;f2=10.2882;omega1=2*pi*f1;omega2=2*pi*f2;A=[1/(2*omega1)omega1/2;1/(2*omega2)omega2/2];xi=[xi1;xi2];x=inv(A)*xi;alpha=x(1,1)beta=x(2,1)%--计算瑞利阻尼矩阵Cr(2*2)--alpha=-1.8801;beta=0.0026;m1=105;m2=7;k1=10000;k2=410;M=[m10;0m2];K=[k1+k2-k2;-k2k2];Cr=alpha*M+beta*K;可知：瑞利比例阻尼系数、0026.0瑞利比例阻尼矩阵12.0947-1.0660-1.0660-3445.170-C2）求解模态矩阵（及特征矢量矩阵）书P23已说明根据粘性比例阻尼振动系统的微分方程所求得的特征矢量与该系统无阻尼振动下求得的特征矢量相等。因此，我们可以利用求此二阶系统在无阻尼振动下的微分方程的特征矢量更简单的得出模态矩阵改写Matlab程序polynomial.m求解此二阶系统在无阻尼振动下的微分方程的特征方程，程序如下：symsx;m1=100;m2=5;k1=10000;k2=500;M=[m10;0m2];K=[k1+k2-k2;-k2k2];y=det(K-x^2*M)解以上求得的多项式：p=[735-1159204100000];x0=roots(p)可知：554653201.、1597.104。将其分别代入回0)(MK，可得：035.1178410-410-767.47862111、0319.1179-410-410-7685.526-2212求得模态矩阵1.2848-11.67501122223）求解频响函数的模态展式矩阵（1）求模态质量矩阵、模态刚度矩阵和模态阻尼矩阵116.5550001394.10591.2848-11.6750117001051.2848-111.67501MmdiagTi121340006721.806251.2848-11.675011410410-410-104101.2848-111.67501KkdiagTi187.6-001843.8-1.2848-11.67501112.0947-1.0660-1.0660-3445.1701.2848-111.67501CcdiagTi（2）由此可得频响函数的模态展式为：221()TiiiiiiHkmjc（3.3）写成矩阵形式为：22111221112212222211122211122211122221222111221221222222111222111222()()()()()kmjckmjckmjckmjcHHHHHkmjckmjckmjckmjc将所求、[]idiagm、[]idiagk、[]idiagc代入：1.65072848.1-1.2848-1)6.1875550.116-121340(305613667501167501118.18431059.1394-8062.56721)()()()()(1122211211j...jHHHHH）（4．脉冲响应函数对（3.3）作傅立叶逆变换，得到脉冲响应函数矩阵：21()siniTtiidiiidihtetm（4.1）5．11()H的幅频、相频、实频、虚频特性曲线以及导纳图和博德图1）11()H的幅频特性曲线：11()H与的关系11222222221111222211()(1)4(1)4Hkk，其中01，i为阻尼比。代入可得：8471)x(0.07074)1(1213401)x(0.03964)7802.521(8062.567211)(22211H（5.1）编写Matlab程序figure1.m画图，程序如下：omega=5:.01:15;o1=omega/sqrt(52.7802);o2=omega/sqrt(105.8471);k1=56721.8062;k2=121340;xi1=0.0396;xi2=0.0707;y11=1./(k1.*sqrt((1-o1.^2).^2+(2.*xi1.*o1).^2))+1./(k2.*sqrt((1-o2.^2).^2+(2.*xi2.*o2).^2));plot(omega,y11,'LineWidth',2);gridonxlabel('频率Hz')ylabel('幅值mm')title('m1的一阶幅频特性')输出图形：2）11()H的相频特性曲线：11与的关系112211221222arctanarctan11。代入可得：-110.2882x0.07072arctan1x0.03962arctan（5.2）编写Matlab程序figure2.m画图，程序如下：omega=5:.01:15;o1=omega/sqrt(52.7802);o2=omega/sqrt(105.8471);k1=56721.8063;k2=121340;xi1=0.0396;xi2=0.0707;y11=atan((-2.*xi1.*o1)./(1-o1.^2))+atan((-2.*xi2.*o2)./(1-o2.^2));plot(omega,y11,'LineWidth',2);axis([515-22])gridonxlabel('频率Hz')ylabel('相位角')title('m1的一阶相频特性')输出图形：3）11()H的实频特性曲线：11()RH与的关系221211222222221111222211()1414RHkk。代入可得：8471)x(0.07074)1(12134084711)x(0.03964)7802.