matlab汽车座椅脉冲振动冲击仿真

1、内容简介
略220-座椅振动

2、内容说明

1、隔振系统微分方程为：

该隔振系统为非线性隔振系统，当，系统称为准零刚度隔振系统。

具体参数为：，，，

2、需要解决的问题

在给定频率下，求解上述微分方程所确定的系统的时域响应。画出响应的响应曲线。如1Hz，2Hz等。

3、改变参数求系统时域响应

如：激励信号变为正弦激励，多频激励、随机激励或冲击激励。

也是分别画出响应曲线。

响应曲线为横轴时间，纵轴位移的曲线。

 

 

 

（1）多频激励为

（2）冲击激励为

加速度激励，，冲击激励总时长t=0.5s。

（3）随机激励比较麻烦，需要用MATLAB将加速度功率谱，应用谐波叠加法得到时域激励信号，然后求解。

随机激励的公式：

路面不平度功率谱：

这里，的值为路面硬度测量值。参考空间频率为，，。

如果车辆以固定的水平速度行驶，对于给定路面来说，由于路面不平度产生的力可以用下面的级数来描述(即为时域激励信号)：

3、仿真分析

clc;
clear all;
close all;
%% bͼ

t0_1 = 0.2;
zmax_1 = 0.1;
sigma_1 = 26;
[t_11,e_11,t_21,u_11]=impluse1(t0_1,zmax_1,sigma_1);
figure
plot(t_11,e_11*1000,t_21,u_11)
xlabel('Time(s)');
ylabel('Displacement(mm)');
title zmax=0.1
legend('Excitation','Response of the system');

t0_2 = 0.2;
zmax_2 = 1;
[t_12,e_12,t_22,u_12]=impluse1(t0_2,zmax_2,sigma_1);
figure
plot(t_12,e_12*1000,t_22,u_12)
xlabel('Time(s)');
ylabel('Displacement(mm)');
title zmax=1
legend('Excitation','Response of the system');

t0_3 = 0.2;
zmax_3 = 5;
[t_13,e_13,t_23,u_13]=impluse1(t0_3,zmax_3,sigma_1);
figure
plot(t_13,e_13*1000,t_23,u_13)
xlabel('Time(s)');
ylabel('Displacement(mm)');
title zmax=5
legend('Excitation','Response of the system');


figure
plot(t_21,u_11,t_22,u_12,t_23,u_13)
xlabel('Time(s)');
ylabel('Displacement(mm)');
title sigma=26
legend('zmax=0.1','zmax=1','zmax=5');

 

4、参考论文
A vibration isolation system in low frequency excitation region using negative stiffness structure for vehicle seat