用矩阵形式表示方程:M*X''+C*X'+K*X=F(t)
方程实际上是六自由度系统的振动响应方程,就是说X是一个6×1的矩阵,M,C,K为质量,阻尼,刚度矩阵。而且C,K不是对角阵,各个自由度之间有耦合。下面是一些参量代码:
function dy=odefun(t,y)
%六自由度常微分非齐次方程
m=87.7;
I_x=2.57;I_y=4.94;I_z=3.04;
a=0.2;b=0.15;h=0.3;h_1=h/2;h_2=h/2;
k_x=1.5e5;k_y=8.75e5;k_z=1e5;
k_ox=5e4;k_oy=5e4;k_oz=5e4;
c_x=0.12*(2*sqrt(m*k_x));c_y=0.12*(2*sqrt(m*k_y));c_z=0.12*(2*sqrt(m*k_z));
c_ox=0.12*(2*sqrt(m*k_ox));c_oy=0.12*(2*sqrt(m*k_oy));c_oz=0.12*(2*sqrt(m*k_oz));
%系统参数
c_11=4*c_z+2*c_oz;
c_22=(4*c_x+2*c_ox)*b^2+(4*c_y+2*c_oy)*a^2;
c_33=4*c_x+2*c_ox;
c_44=4*c_x*h_1^2+2*c_oy*h_2^2+(4*c_z+2*c_oz)*b^2;
c_55=4*c_