不要意思 程序有点长,  只看前面和后面就可以了哈。就是 根据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