参考

A comparative study of several calculation methods for ship shock spectrum道客巴巴

Cui Jidong's response spectrum calculation softwareSPECTR

问题

用AbaqusThe results of the acceleration shock response for simulation calculations are compared with SPECTRPartially inconsistent,Try to find out why.

SPECTR计算结果

matlab simulinkResponse spectrum calculation results

结果与Abaqus仿真结果基本一致,Acceleration response from14HZ到18HzThere is a drop in the middle,而SPECTRIt keeps going up.

matlabControl codes and simulink模型

%RSget%
fn=[0,2.2,6,10,14,16,18,20];%Natural frequencies of a single degree of freedom system
wn=2*pi*fn;
k1=wn.^2;
m=1;
s_=0.05;
wd=sqrt(1-s_^2)*wn;
c1=2*s_*wn*m;
modelname=bdroot;%Get the currently active model name
stoptime=1;
Num=length(k1);%Response spectrum solution points
Color=[0.5,0,0;
    1,0,0;
    0,0.5,0;
    0,1,0;
    0,0,0.5;
    0,0,1;
    0.5,0.5,0;
    0.5,0,0.5
    0,0.5,0.5
    0.5,05,0.5];
for i=1:Num
    k=k1(i);
    c=c1(i);
    simout{
    i}=sim(modelname,'StopTime',num2str(stoptime));%针对每组输入,Return the output separately % 仿真结束后,你可以利用simoutGet simulation results
    %Get all output results
    plot(simout{
    i}.get('tout'),simout{
    i}.get('yout'),'color',Color(i,:));
    hold on;
    RS_min(i)=min(simout{
    i}.get('yout'));
    RS_max(i)=max(simout{
    i}.get('yout'));
end
figure;
plot(fn,[RS_min;RS_max;-RS_min])