参考

舰船冲击谱若干计算方法比较研究道客巴巴

崔济东的反应谱计算软件SPECTR

问题

用Abaqus进行仿真计算的加速度冲击响应的结果与SPECTR部分不一致，尝试找一下原因。

SPECTR计算结果

matlab simulink响应谱计算结果

结果与Abaqus仿真结果基本一致，加速度响应从14HZ到18Hz中途有下降，而SPECTR则是一直上升。

matlab控制代码及simulink模型

%RSget%
fn=[0,2.2,6,10,14,16,18,20];%单自由度系统固有频率
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;%获取当前激活模型名
stoptime=1;
Num=length(k1);%响应谱求解点数
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));%针对每组输入，分别返回输出 % 仿真结束后，你可以利用simout获取仿真结果
    %获取所有的输出结果
    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])

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34