龙格库塔求解车辆振动模型_多自由度振动系统的龙格库塔法

1、内容简介
略

2、内容说明
略

3、仿真分析

clear
clc
close all
t0=0; % 开始时间
t1=10;%结束时间
h=0.01;%步长
u = 1; %输入幅值
x0=[0;u;0;1;0;1;0;1];  %初始状态
T=t0:h:t1; %时间序列
X=zeros(length(x0),length(T)); %
F = u*sin(T); %激励源
Fdt = u*cos(T); %激励源速度
Fdtdt = -u*sin(T);%激励加速度
X(:,1)=x0;
m1 = 50;% 轮胎重量
m2 = 1000;% 车体重量
m3 = 80;% 人体重量
k1 = 1800;c1 = 200;     %轮胎刚度阻尼
k2 = 14000;c2 = 5000;   %悬挂刚度阻尼
k3 = 1000; c3 = 500;    %座椅刚度和阻尼

for j=1:length(T)-1
    X(:,j+1)=RK45(T(j),X(:,j),@(t,x)f(t,x),h); %调用龙格库塔函数
end
%输出信息
figure
plot(T,F);
title Z1位移
figure
plot(T,Fdt);
title Z1速度
figure
plot(T,Fdtdt);
title Z1加速度

4、参考论文
略
企q鹅q号：1762016542