使用Matlab提取骨架曲线与耗能能力-Toy模板网

这篇具有很好参考价值的文章主要介绍了使用Matlab提取骨架曲线与耗能能力。希望对大家有所帮助。如果存在错误或未考虑完全的地方，请大家不吝赐教，您也可以点击"举报违法"按钮提交疑问。

骨架曲线的提取

代码网上很多，但是却没有教程，所以代码拿来一跑就出错，为了以后方便复习及后面的师弟师妹好使用，在此做个手把手教程。在此感谢133办公室曾师兄，没有他的指导代码仍是跑不动的。

先上代码，点击matlab中的新建脚本，将此代码复制进去修改即可。

使用Matlab提取骨架曲线与耗能能力

clear all
clc;
hyster=importdata('D:\matlab\bin\xjx.mat');%导入数据
LineNum=size(hyster,1);
EndCircle=1;
%--------------------找各级循环分割点位置----------------------负->正or负->零
for i=3:LineNum
if (hyster(i,1)*hyster(i-1,1)<=0)&&(hyster(i-1,1)<0);%--默认初始加载从正半周开始
EndCirclePoints(EndCircle)=i; %分割点的行号EndCirclePoints
EndCircle=EndCircle+1; %中间分割点个数为EndCircle-1
end;
end;

%--------------------检查最后是否存在不完整滞回环----------------------
if EndCirclePoints(EndCircle-1)<LineNum
LoopNum=length(EndCirclePoints)+1; %LoopNum仅为完整循环的个数，应还有最后一个不完整循环
else
LoopNum=length(EndCirclePoints);
end;
hyster=hyster(1:EndCirclePoints(EndCircle-1),:);
deflection=hyster(:,1);%位移
Force=hyster(:,2);%力
% figure()
% plot(deflection,Force,'r') %%绘图
% hold on
Vmp=max(Force(:,1));%%最大荷载(正向）
Vmn=min(Force(:,1));%%最大荷载(负向）
dmp=deflection(find(hyster(:,2)==Vmp),1);%%最大位移（正向）
dmn=max(deflection(find(hyster(:,2)==Vmn),1));%%最大位移（负向）
% plot(dmp,Vmp,'*',dmn,Vmn,'*')
% hold on
Vcrack_p=0.85*Vmp;%%峰值85%
Vcrack_n=0.85*Vmn;%%峰值85%
% plot([0,40],[Vmp,Vmp],'b',[0,40],[Vp,Vp],'k',[-40,0],[Vmn,Vmn],'b',[-40,0],[Vn,Vn],'k')
Vp=0.7*Vmp;%%峰值70%
Vn=0.7*Vmn;%%峰值70%
%% 正向