clear;clc;
% 框架参数 ================================================================
m1=4.7835e5;
m2=4.7835e5;
m3=5.1779e5;
k1=1.9947e8;
k2=2.3635e8;
k3=1.3931e8;
M=[...
m1 0 0
0 m2 0
0 0 m3
];
K=[...
k1+k2 -k2 0
-k2 k2+k3 -k3
0 -k3 k3
];
C=[...
2.7421 -0.8976 -0.1596
-0.8976 2.4274 -0.7353
-0.1596 -0.7353 1.5075
].*1e5;
invM=inv(M);
E=[...
1
1
1];
% =========================================================================
t_end=20;
step = 0.02;
ts(:,1)= 0 : step : t_end; % 地面加速采样时间
ag = sin(ts); % 地面加速度
Ft = [ts, -ag*(M*E)'];
% t_end = ts(end);
opts=simset('FixedStep',step);
sim('frame',t_end,opts);
评论1