7自由度整车振动模型附带matlab计算程序

kyu16866

当一只路面四个车轮的不平度功率谱后，容易求得输出响应的频率响应函数矩阵，和响应的频谱矩阵，如果求车身质心加速度频率响应是不是abs(H(1,1)+H(1,2)+H(1,3)+H(1,4))*(2*pi*f)^2，四个车轮输入对车身的线性叠加？
%汽车7自由度振动响应计算
clear all
m=1235;Jx=472;Jy=1731;v=80*1000/60/60;
csp=4e-7;

mt1=49;mt2=49;mt3=49;mt4=49;
M=[m 0 0 0 0 0 0;0 Jx 0 0 0 0 0;0 0 Jy 0 0 0 0;0 0 0 mt1 0 0 0;0 0 0 0 mt2 0
0;...
    0 0 0 0 0 mt3 0;0 0 0 0 0 0 mt4];
k1=22741;k2=22741;k3=26144;k4=26144;
kt1=302342;kt2=302342;kt3=492982;kt4=492982;
c1=1228;c2=1228;c3=1210;c4=1210;
ct1=430;ct2=430;ct3=430;ct4=430;
l=2.655;%轴距
l1=1.1135;l2=1.5415;l3=0.7075;l4=0.7075;%l1,l2前轴距与后轴距，l3,l4：1/2轮距
C=[c1+c2+c3+c4 -l1*c1+l2*c2-l1*c3+l2*c4 -l3*c1-l3*c2+l4*c3+l4*c4 -c1 -c2 -c3 -
c4;...
    -l1*c1+l2*c2-l1*c3+l2*c4 l1^2*c1+l2^2*c2+l1^2*c3+l2^2*c4 l1*l3*c1-l2*l3*c2-
l1*l4*c3+l2*l4*c4 l1*c1 l2*c2 l1*c3 l2*c4;...
    -l3*c1-l3*c2+l4*c3+l4*c4 l1*l3*c1-l2*l3*c2-l1*l4*c3+l2*l4*c4...
    l3^2*c1+l3^2*c2+l4^2*c3+l4^2*c4 l3*c1 l3*c2 -l4*c3 -l4*c4;...
    -c1 l1*c1 l3*c1 c1+ct1 0 0 0;...
    -c2 -l2*c2 l3*c2 0 c2+ct2 0 0;...
    -c3 l1*c3 -l4*c3 0 0 c3+ct3 0;...
    -c4 -l2*c4 -l4*c4 0 0 0 c4+ct4];
K=[k1+k2+k3+k4 -l1*k1+l2*k2-l1*k3+l2*k4 -l3*k1-l3*k2+l4*k3+l4*k4 -k1 -k2 -k3 -
k4;...
    -l1*k1+l2*k2-l1*k3+l2*k4 l1^2*k1+l2^2*k2+l1^2*k3+l2^2*k4 l1*l3*k1-l2*l3*k2-
l1*l4*k3+l2*l4*k4 l1*k1 l2*k2 l1*k3 l2*k4;...
    -l3*k1-l3*k2+l4*k3+l4*k4 l1*l3*k1-l2*l3*k2-l1*l4*k3+l2*l4*k4...
    l3^2*k1+l3^2*k2+l4^2*k3+l4^2*k4 l3*k1 l3*k2 -l4*k3 -l4*k4;...
    -k1 l1*k1 l3*k1 k1+kt1 0 0 0;...
    -k2 -l2*k2 l3*k2 0 k2+kt2 0 0;...
    -k3 l1*k3 -l4*k3 0 0 k3+kt3 0;...
    -k4 -l2*k4 -l4*k4 0 0 0 k4+kt4];
Ct=[0 0 0 0;0 0 0 0;0 0 0 0;ct1 0 0 0;0 ct2 0 0;0 0 ct3 0;0 0 0 ct4];%轮胎阻尼矩阵
Kt=[0 0 0 0;0 0 0 0;0 0 0 0;kt1 0 0 0;0 kt2 0 0;0 0 kt3 0;0 0 0 kt4];%轮胎刚度矩阵
for n=1:1000
     f=0.1*n;
Kq=Kt+2*pi*f*j*Ct;
H=inv(K-(2*pi*f)^2*m+(2*pi*f)*j*C)*Kq;%(:,:,n)
四轮路面不平度功率谱矩阵
sq=csp*v/f/f;
tao=1/v;
tt=2*pi*f*tao*j;
Sq=sq*[1 exp(-tt) 0.18 0.18*exp(-tt);exp(-tt) 1 0.18*exp(-tt) 0.18;...
     0.18 0.18*exp(-tt) 1 exp(-tt);0.18*exp(-tt) 0.18 exp(-tt) 1];
Sz1=conj(H)*Sq*H.';%响应功率谱密度矩阵
% Sz=abs(Sz1);
% Szq(n)=abs(Sz(1,1))*(2*pi*f)^4;%车身质心加速度功率谱
% Sd1(n)=Sz(1,1)-l1*Sz(2,2)-l3*Sz(3,3)-Sz(4,4);
% Sd2(n)=Sz(1,1)-l2*Sz(2,2)-l3*Sz(3,3)-Sz(4,4);%Sz(5,5);
% Sd3(n)=Sz(1,1)-l1*Sz(2,2)+l4*Sz(3,3)-Sz(4,4);%Sz(6,6);
% Sd4(n)=Sz(1,1)+l2*Sz(2,2)+l4*Sz(3,3)-Sz(4,4);%Sz(7,7);
% Sf1(n)=sqrt((k1*Sd1)^2+(c1*2*pi*f*Sd1)^2);
% Sf2(n)=sqrt((k2*Sd2)^2+(c2*2*pi*f*Sd2)^2);
% Sf3(n)=sqrt((k3*Sd3)^2+(c3*2*pi*f*Sd3)^2);
% Sf4(n)=sqrt((k4*Sd4)^2+(c4*2*pi*f*Sd4)^2);
% Df(n)=abs(H(1,1)-H(4,1));%-l1*H(2,1)-l3*H(3,1)
Dj(n)=abs(H(1,1)+H(1,2)+H(1,3)+H(1,4))*(2*pi*f)^2;%车身质心加速度频率响应
end
f1=0.1:0.1:100;
% plot(f1,Szq)
plot(f1,Dj)

画出来的Dj是车身质心加速度频率响应吗？

kyu16866

求出的频率响应函数H是7x4的矩阵，代表4个输入对7个输出Hij (i=1-7,j=1-4)可我想要4个车轮输入对一个输出的影响  请高手指导，，，，，，，

kyu16866

哈哈，终于理解  这个应该怎么求四个输入对某一个响应的功率谱密度啦！！！！

若无其事

学长，求联系方式

lzplzp

同样的问题

zxx4477

csp=4e-7  是什么， 不平度系数？