Bouc Wen_Bouc_Wen_bouc_bouc-wen模型_Bouc-wen_Boucwen

clear all;close all;clc; ratio=0.025; k=500000; exp=5; Fy=250000; Dy=0.5; x=0:0.01:5.005; dt=0.005; nt=5.005/0.005+1; load ('weiyi.txt') d=weiyi(:,2,1); f=x; z0=0; z=0; for i=2:(length(d)) u=d(i); du=(d(i)-d(i-1))/10; for j=1:10 sign=0; if(du*z>0) sign=1; end dz=k/Fy*du*(1-sign*power(abs(z),exp)); z=z+dz; end f(i)=ratio*k*u+(1-ratio)*k*Dy*z; end plot(d,f);