怎么免费申请网站/武汉seo排名优化公司
不要意思 程序有点长, 只看前面和后面就可以了哈。就是 根据epsilon0(i) 来求最后面那个式子的值。。。拜托了哈
for i=2:(500)
(省略)
epsilon0(i)=(DN-YN)/L;
if epsilon0(i)>=0
if epsilon0(i)>epsilon0(i-1)
if 0<=epsilon0(i) && epsilon0(i)<=0.0023
if epsilon0(i-1)==0
sigma0(i)=Ea*epsilon0(i);
GANG=Ea;
elseif epsilon0(i)~=0
aa1=epsilon0min;bb1=epsilontmin;cc1=sigma0min; % this is the point
epsilont(i)=bb1;
EF=Em*Ea/(bb1*(Ea-Em)/H+Em);
sigma0(i)=cc1+EF*(epsilon0(i)-aa1);
GANG=EF;
end
elseif 0.0023<=epsilon0(i) && epsilon0(i)<=0.0596
aa1=epsilon0min;bb1=epsilontmin;cc1=sigma0min; % this is the point
epsilont(i)=bb1;
% this is the condition for minor loop judgement
sigmatp1=sigmams1+epsilontmin*(sigmamf1-sigmams1)/H;
EF=Em*Ea/(bb1*(Ea-Em)/H+Em);
sigma0(i)=cc1+EF*(epsilon0(i)-aa1);
GANG=EF;
if sigma0(i)>sigmatp1
epsilont(i)=H*((epsilon0(i)-epsilonms1)/(epsilonmf1-epsilonms1));
sigma0(i)=sigmams1+epsilont(i)*(sigmamf1-sigmams1)/H;
a0=epsilon0(i);b0=epsilont(i);c0=sigma0(i); % this is the point
GANG=(sigmamf1-sigmams1)/(epsilonmf1-epsilonms1);
end
elseif epsilon0(i)>0.0596
sigma0(i)=sigmamf1+Em*(epsilon0(i)-epsilonmf1);
GANG=Em;
end
elseif epsilon0(i)-epsilon0(i-1)<0
epsilon0max=a0;epsilontmax=b0;sigma0max=c0;
epsilont(i)=b0;
epsilontp3=epsilonaf1+epsilontmax*(epsilonas1-epsilonaf1)/H;
sigmatp3=sigmaaf1+epsilontmax*(sigmaas1-sigmaaf1)/H;
if epsilon0(i)<=epsilon0max && epsilon0(i)>=epsilontp3;
ER=Em*Ea/(epsilontmax*(Ea-Em)/H+Em);
sigma0(i)=sigma0max+ER*(epsilon0(i)-epsilon0max);
GANG=ER;
elseif epsilon0(i)=0
epsilont(i)=epsilontmax-epsilontmax*((epsilontp3-epsilon0(i))/(epsilontp3-epsilonaf1));
sigma0(i)=sigmaaf1+epsilont(i)*(sigmatp3-sigmaaf1)/epsilontmax;
epsilon0min=epsilon0(i);
epsilontmin=epsilont(i);
sigma0min=sigma0(i); % this is the point
GANG=(sigmaas1-sigmaaf1)/(epsilonas1-epsilonaf1);
elseif epsilon0(i)>epsilontmax
sigma0(i)=sigmamf1+Em*(epsilon0(i)-epsilonmf1);
GANG=Em;
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
elseif epsilon0(i)<0
if epsilon0(i)
if 0>=epsilon0(i) && epsilon0(i)>=-0.0023
if epsilon0(i-1)==0
sigma0(i)=Ea*epsilon0(i);
GANG=Ea;
elseif epsilon0(i)~=0
aa11=epsilon00min;bb11=epsilonttmin;cc11=sigma00min; % this is the point
epsilont(i)=bb11;
EFF=Em*Ea/(bb11*(Ea-Em)/HH+Em);
sigma0(i)=cc11+EFF*(epsilon0(i)-aa11);
GANG=EFF;
end
elseif -0.0023>=epsilon0(i) && epsilon0(i)>=-0.0596
aa11=epsilon00min;bb11=epsilonttmin;cc11=sigma00min; % this is the point
epsilont(i)=bb11;
% this is the condition for minor loop judgement
sigmatp11=sigmams2+epsilonttmin*(sigmamf2-sigmams2)/HH;
EFF=Em*Ea/(bb11*(Ea-Em)/HH+Em);
sigma0(i)=cc11+EFF*(epsilon0(i)-aa11);
GANG=EFF;
if sigma0(i)
epsilont(i)=HH*((epsilon0(i)-epsilonms2)/(epsilonmf2-epsilonms2));
sigma0(i)=sigmams2+epsilont(i)*(sigmamf2-sigmams2)/HH;
a00=epsilon0(i);b00=epsilont(i);c00=sigma0(i); % this is the point
GANG=(sigmamf2-sigmams2)/(epsilonmf2-sigmamf2);
end
elseif epsilon0(i)
sigma0(i)=sigmamf2+Em*(epsilon0(i)-epsilonmf2);
GANG=Em;
end
elseif epsilon0(i)-epsilon0(i-1)>0
epsilon00max=a00;epsilonttmax=b00;sigma00max=c00;
epsilont(i)=b00;
epsilontp33=epsilonaf2+epsilonttmax*(epsilonas2-epsilonaf2)/HH;
sigmatp33=sigmaaf2+epsilonttmax*(sigmaas2-sigmaaf2)/HH;
if epsilon0(i)>=epsilon00max && epsilon0(i)<=epsilontp33;
ERR=Em*Ea/(epsilonttmax*(Ea-Em)/HH+Em);
sigma0(i)=sigma00max+ERR*(epsilon0(i)-epsilon00max);
GANG=ERR;
elseif epsilon0(i)>epsilontp33 && epsilon0(i)<=0
epsilont(i)=epsilonttmax-epsilonttmax*((epsilontp33-epsilon0(i))/(epsilontp33-epsilonaf2));
sigma0(i)=sigmaaf2+epsilont(i)*(sigmatp33-sigmaaf2)/epsilonttmax;
epsilon00min=epsilon0(i);
epsilonttmin=epsilont(i);
sigma00min=sigma0(i); % this is the point
GANG=(sigmaaf2-sigmaas2)/(epsilonaf2-epsilonas2);
elseif epsilon0(i)
sigma0(i)=sigmamf2+Em*(epsilon0(i)-epsilonmf2);
GANG=Em;
end
end
EPSILONT=epsilont(i); ELAST=GANG;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
(省略)
end
end