Boun-Wen.rar

function Ydot=boucwen(t,Y,A,b,r,n,m,c,k,a,F) %该函数表示Bouc-Wen模型 %t为时间，Y为结构的状态； %A,b,r,n为Bouc-Wen模型的参数； %m,c,k,a为结构的质量矩阵、阻尼矩阵、弹性刚度矩阵、屈服后与屈服前刚度比； %F为外荷载； %Ydot为Y的一阶导数； cn=length(m); Ydot=[Y([cn+1:2*cn]); -inv(m)*k*a*Y(1:cn)-inv(m)*c*Y((cn+1):2*cn)-inv(m)*k*(1-a)*Y([(2*cn+1):3*cn])-inv(m)*F; %(3.25)第二式F(t)项错误，此处已更正 A.*Y([(cn+1):2*cn])-b.*abs(Y([(cn+1):2*cn])).*abs(Y([(2*cn+1):3*cn])).^(n-1).*Y([(2*cn+1):3*cn])-r.*Y([(cn+1):2*cn]).*abs(Y([(2*cn+1):3*cn])).^n]; %(3.24),(3.25)