|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
fun.m文件
function F=fun(dei,Cd,Fd,beta,omiga,Fti)
F = Cd.*dei+Fd.*(1-exp(-beta.*dei./omiga)./(1+exp(-beta.*dei./omiga)))-Fti;
sigmoid.m文件
%sigmoid模型%---------力计算位移
%F(t)=Cd*de+Fd*sgn(de)-F0=Kd*(x-e)-F0
%F(t)=Fd*(1-exp(-beta*de/omiga))/(1+exp(-beta*de/omiga))+Cd*de
clear all
Cd=8.684*1e3;%模型参数
Fd=1400;
beta=541;
omiga=12.56;
Hz=2;%sinwave frequency
dt=0.002;%荷载时间步长
dF=25;%荷载步长
Fmax=1500;%荷载范围
Ft=1500.*sin(12.56.*(0:dt:10));
de=zeros(1,length(Ft)+1);%预定义变量---内速度
e=de;%预定义变量---内位移
for i=1:length(Ft)
de(i+1)=fzero(@de(i) fun(de(i),Cd,Fd,beta,omiga,Fti),de(i)); %Cd.*dei+Fd.*(1-exp(-beta.*dei./omiga)./(1+exp(-beta.*dei./omiga)))-Fti
e(i+1)=e(i)+de(i+1).*dt;%积分
end
comet(e(1,2:end),Ft)
%comet(de(1,2:end),Ft)
但是最后运行不了 |
|