matlab-miner,修正miner,Manson,Corten-Dolan的疲劳寿命预测算法

clear; A=[ 0.925 0.295 1991 1.005 0.375 30733 1.080 0.300 1991 1.000 0.220 16451 0.925 0.455 1938 1.080 0.450 1944 1.000 0.530 1940 1.160 0.530 121 1.005 0.685 340 0.850 0.530 133 1.080 0.610 278 1.235 0.605 10 0.925 0.605 263 0.845 0.685 21 0.770 0.610 5 1.315 0.685 1 1.080 0.760 56 0.925 0.765 56 1.155 0.685 12 0.850 0.840 4 1.235 0.765 2 1.000 0.840 57 0.770 0.760 2 1.160 0.840 8 0.925 0.915 14 1.005 0.995 26 1.230 1.070 1 0.770 1.070 1 0.845 0.995 2 1.080 0.920 11 1.080 1.070 6 0.770 0.920 1 1.000 1.150 11 0.925 1.075 4 1.155 0.995 1 1.000 1.300 4 1.000 1.500 2 1.000 1.600 1]; jizhun=156; % 156 130 % 182 156 130 % 156 130 104 % 164 148 % 180 164 148 % 148 117 102 B(:,2)=A(:,1)*jizhun; B(:,1)=A(:,2)*jizhun; B(:,3)=A(:,3); Sa=1*jizhun; Sm=1.6*jizhun; % N1=1.93e10/(430/(430-Sm)*Sa)^2.55; N1=1.08e12/(430/(430-Sm)*Sa)^3.26; % N1=1.53e8/(430/(430-Sm)*Sa-11)^2.03; % N1=2.14e8/(488/(488-Sm)*Sa-57)^2.04; % N1=3.92e8/(488/(488-Sm)*Sa-36)^1.91; % N1=1.33e9/(488/(488-Sm)*Sa-8.4)^2.4; N2=426991;%% % 189271 380307 % 186140 426991 1042064 % 20404 38856 83375 % 144183 362211 % 113330 174402 286725 % 47261 114516 189823 Fai=1/log(N1/N2)*log(log(0.35*((N1/N2)^0.25))/log(1-0.65*((N1/N2)^0.25))); Z=log(0.35*((N1/N2)^0.25))/(N1^(Fai)); for i=1:38 % B(i,4)=((B(i,1)*430./(430-B(i,2)))).^2.55; %156 130 B(i,4)=((B(i,1)*430./(430-B(i,2)))).^3.26; % 182 156 130 % B(i,4)=((B(i,1)*430./(430-B(i,2)))-11).^2.03;% 156 130 104 % B(i,4)=((B(i,1)*488./(488-B(i,2)))-57).^2.04;% 164 148 % B(i,4)=((B(i,1)*488./(488-B(i,2)))-36).^1.91; % % B(i,4)=((B(i,1)*488./(488-B(i,2)))-8.4).^2.4;%101 % B(i,5)=(1.93*10^10)/B(i,4); B(i,5)=(1.08*10^12)/B(i,4); % B(i,5)=(1.53*10^8)/B(i,4); % B(i,5)=(2.14*10^8)/B(i,4); % B(i,5)=(3.92*10^8)/B(i,4); % B(i,5)=(1.33*10^9)/B(i,4); B(i,6)=B(i,5).*exp(Z.*(B(i,5).^Fai)); B(i,7)=B(i,5)-B(i,6); B(i,8)=B(i,3)./B(i,6); B(i,9)=B(i,3)./B(i,7); end T=1/sum(B(:,8))+1/sum(B(:,9)) shiyan=6.92; % 4.19 4.87 % 1.68 6.92 15.56 % 1.52 3.21 6.76 % 4.12 4.6 % 2.49 3.88 6.40 % 2.18 2.97 5.96 Bai=abs(shiyan-T)/shiyan