MATLAB 编程 关于刚度退化模型 clough双线性

function [tag,tag_2,Pp,F_present,kk]=clough_model(tag1,X_present,X_velocity,original_stiff,X_yeild) %% 1 tag=tag1 ; if tag==0 %从0点加载 kk=original_stiff(:,1); plastic=0 ; F_present=kk*X_present+plastic; if X_present>X_yeild %%正向位移大于屈服位移 tag=1; tag_2=1; Pp=(X_yeild-X_previous)/(X_present-X_previous); %对第一类拐点处理 kk=original_stiff(:,2); plastic=(original_stiff(:,1)-original_stiff(:,2))*X_yeild; F_present=kk*X_present+plastic; elseif X_present<-X_yeild %%反向位移大于屈服位移 tag=-1; tag_2=1; Pp=(-X_yeild-X_previous)/(X_present-X_previous); kk=original_stiff(:,2); plastic=-(original_stiff(:,1)-original_stiff(:,2))*X_yeild; F_present=kk*X_present+plastic; end end %% 2这里是第二类拐点可不做处理 if tag==1 %正向进入强化阶段 kk=original_stiff(:,2); plastic=(original_stiff(:,1)-original_stiff(:,2))*X_yeild; F_present=kk*X_present+plastic; if X_velocity<0 %进入屈服后的卸载阶段 X_ultimate_1=X_present; F_ultimate_1=F_present; tag=2; kk=original_stiff(:,1); plastic=F_ultimate_1-kk*X_ultimate_1; F_present=kk*X_present+plastic ; end end %% 3 if tag==-1 %方向进入强化阶段 kk=original_stiff(:,2); plastic=-(original_stiff(:,1)-original_stiff(:,2))*X_yeild; F_present=kk*X_present+plastic; if X_velocity >0 NX_ultimate_1=X_present; NF_ultimate_1=F_present; tag=-2; kk=original_stiff(:,1); plastic=NF_ultimate_1-kk*NX_ultimate_1; F_present=kk*X_present+plastic; end end %% 4 if tag==2 %正向屈服后的卸载阶段 kk=original_stiff(:,1) ; plastic=F_ultimate_1-kk*X_ultimate_1; F_present=kk*X_present+plastic; if F_present<0 %卸载至反向加载阶段 tag=-3; tag_2=1; Pp=F_previous/(F_previous-F_present); %第三类拐点处理 X_Fzero_1=X_previous+Pp*(X_present-X_previous); F_zero_1=0; kk=(NF_ultimate_1-F_zero_1)/(NX_ultimate_1-X_Fzero_1); plastic=F_zero_1-kk*X_Fzero_1; F_present=kk* X_present+plastic; elseif X_present>X_ultimate_1 %加载至强化阶段 tag=1; tag_2=1; Pp=(X_ultimate_1-X_previous)/(X_present-X_previous); kk=original_stiff(:,2); plastic=(original_stiff(:,1)-original_stiff(:,2))*X_yeild; F_present=kk*X_present+plastic; end end %% 5 if tag==-2 %屈服后反方向卸载阶段 kk=original_stiff(:,1); plastic=NF_ultimate_1-kk*NX_ultimate_1; F_present=kk*X_present+plastic; if F_present>0 %加载阶段 tag=3; tag_2=1; Pp=F_previous/(F_previous-F_present); NX_Fzero_1=X_previous+Pp*(X_present-X_previous); NF_zero_1=0; kk=(F_ultimate_1-NF_zero_1)/(X_ultimate_1- NX_Fzero_1) ; plastic=NF_zero_1-kk*NX_Fzero_1; F_present=kk*X_present+plastic; end end %% 6 if tag==-3 kk=(NF_ultimate_1-F_zero_1)/(NX_ultimate_1-X_Fzero_1); plastic=F_zero_1-kk*X_Fzero_1; F_present=kk*X_present+plastic; if X_velocity>0 tag=-2; kk=original_stiff(:,1); plastic=NF_ultimate_1-kk* NX_ultimate_1; F_present=kk*X_present+plastic; end end %% 7 if tag==3 %反向-正向加载 kk=(F_ultimate_1-NF_zero_1)/(X_ultimate_1-NX_Fzero_1); plastic=NF_zero_1-kk*NX_Fzero_1; F_present=kk*X_present+plastic ; if X_velocity<0 tag=2; kk=original_stiff(:,1); plastic=F_ultimate_1-kk*X_ultimate_1; F_present=kk*X_present+plastic; elseif X_present>X_ultimate_1 tag =1; tag_2=1 ; Pp=(X_ultimate_1-X_previous)/(X_present-X_previous); kk=original_stiff(:,2); plastic=(original_stiff(:,1)-original_stiff(:,2))*X_yeild; F_present=kk*X_present+plastic; end end