clear?all;
clc;
%%%定義Bouc-Wen模型參數%%%
k0=7.1;
a=0.114;
A=1.0;
beta=
gamma=
z（1）=0;

v=4.0;??????????????%定義單向剪切速率4?1/s
gammaeqmax=2;
T=gammaeqmax/v;
num=500;
N=num+1;
dt=T/num;
t=0:dt:T;
%%%%%%%%%%%
for?i=1:N
????u（i）=v*t（i）;
????z（i+1）=z（i）+dz;
????dz=
????dgammaeq（i）=v;
????taoeq（i）=2*gammaeq（i）*（C1eq+2*C2eq*gammaeq（i）^2+3*C3eq*gammaeq（i）^4）;
????gammaov（1）=0;
????taoov（i）=2*gammaov（i）*（C1ov+2*C2ov*gammaov（i）^2+3*C3ov*gammaov（i）^4）;
????dgammaov（i）=dgammaeq（i）-2*taoov（i）*（2/3*gammaov（i）^2+1）/eta（i）;
????gammao

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1120??2016-01-12?19:49??BoucWen.m