八级载荷谱计算：软件助力精准分析

封面文献中提到8级载荷谱的划分方法：

相应步骤：

1、通过用Ncode的进行雨流计数 的方法，统计出极值，

2、求出幅值和均值，由于载荷谱时域数据是应变数据，所以参考弹性模量 的参数值，转换成应力。

3、去除小载荷

4、设定分组间隔

5、求出每组数据频数

ps： 下图中的分组用了matlab 自带函数hist，这个函数用于等间隔分组没问题，用在非等间隔分组就会有问题。后续会给出一个针对该论文的八级非等间隔分组的自己编写的函数 hist_8 .

Matlab代码

function [p n N]=f_j(t,m,i)
% p 为估计参数

% t为极值数据；m截断下限值；e为应力应变换算系数，一般取0.2; i为样本量

f=0.2*(t(:,1)-t(:,2));j=0.2*(t(:,1)+t(:,2))/2;

f_f=find(f<m); % 需要提出数据的位置标示
f(f_f)=[];j(f_f)=[]; %截断数据
size_f=size(f)
size_j=size(j)

f_g=1200*f./(1200-j);

N=[f,j]; %提取截断后的数据


  lev=max(f_g)*[0.0625	0.2	0.35	0.5	0.65	0.7875	0.9	0.975];
     f_g_hist=hist(f_g,lev)
     figure(2*i-1);hist(f_g,lev); %幅度直方图;
   saveas(gcf,int2str(2*i-1),'jpg');
figure(2*i);wblplot(f_g);   %幅值wpp图
  saveas(gcf,int2str(2*i),'jpg');
     max_f_g=max(f_g);
n=length(f);

p1=wblfit(f_g);

p=[p1 f_g_hist max_f_g];

for j=2*i-1:2*i; % 关闭图片
     figure(j);close
end
end 

武汉格发信息技术有限公司，格发许可优化管理系统可以帮你评估贵公司软件许可的真实需求，再低成本合规性管理软件许可,帮助贵司提高软件投资回报率，为软件采购、使用提供科学决策依据。支持的软件有: CAD,CAE,PDM,PLM,Catia,Ugnx, AutoCAD, Pro/E, Solidworks ,Hyperworks, Protel,CAXA,OpenWorks LandMark,MATLAB,Enovia,Winchill,TeamCenter,MathCAD,Ansys, Abaqus,ls-dyna, Fluent, MSC,Bentley,License,UG,ug,catia,Dassault Systèmes,AutoDesk,Altair,autocad,PTC,SolidWorks,Ansys,Siemens PLM Software,Paradigm,Mathworks,Borland,AVEVA,ESRI,hP,Solibri,Progman,Leica,Cadence,IBM,SIMULIA,Citrix,Sybase,Schlumberger,MSC Products...