matlab四分之一模型车+人体和座椅

Posted 2022-05-04 studyer_domi

1、内容简介
略

2、内容说明

求四分之一车模型的传递率曲线：

1.：随机路面作为输入

2：四分之一整车模型

3.：在模型上有四自由度人-座椅模型

4.：输出是：座椅的传递率曲线

 

下图是座椅模型

3、仿真分析

clc
clear all
close all
%% ---------基本参数------
m0_1=65;%非簧载质量
m0_2=325;%簧载质量
k0_1=19e4;%轮胎刚度
c0_1=4e2;%轮胎阻尼
k0_2=18e3;%悬架刚度
c0_2=16e2;%悬架阻尼

%% Specify the model name
model = 'model';

%% Specify the analysis I/Os
% Specify block name as the analysis I/Os
% to linearize the block model/model
io = 'model/model';

%% Specify the operating point
% Use the model initial condition
op = operpoint(model);

%% Linearize the model
sys = linearize(model,io,op);

%% Plot the resulting linearization
step(sys)
sys_new = tf(sys);
f = 0.1:0.1:1000;
s = 1i*2*pi*f;
num = 6.452e07*s.^5 + 5.61e10*s.^4 + 1.351e13*s.^3+ 1.168e15*s.^2 + 5.285e16*s + 4.651e17;
den = s.^12 + 385.3*s.^11 + 1.235e05*s.^10 + 2.307e07*s.^9+ 3.2e09*s.^8 + 3.073e11*s.^7 + 2.131e13*s.^6+ ...
    1.023e15*s.^5 + 3.404e16*s.^4+ 6.541e17*s.^3 + 3.284e18*s.^2+ 2.338e19*s + 4.651e17;
ss = num./den;
semilogy(f,abs(ss))
grid;
xlabel 频率/Hz
ylabel 相对路面传递率
%% 
f = 0.1:0.01:20;
s = 1i*2*pi*f;
num = 4980*s.^4 + 4.275e06*s.^3 + 9.945e08*s.^2 + 7.893e10*s + 3.191e12;
den = s.^8 + 350.1*s.^7 + 1.079e05*s.^6 + 1.811e07*s.^5 + 2.204e09*s.^4+ 1.688e11*s.^3 + 7.885e12*s.^2 + 1.602e14*s + 3.191e12;
ss = num./den;
semilogy(f,abs(ss))
grid;
xlabel 频率/Hz
ylabel 相对车体传递率

 

 

 

4、参考论文
略