查看: 134|回复: 0|关注: 0

[未答复] 关于方程求解,为什么得到的结果有问题?

[复制链接]

新手

5 麦片

财富积分


050


3

主题

3

帖子

0

最佳答案
发表于 6 天前 | 显示全部楼层 |阅读模式
请各位大佬救救孩子吧!
这是文章的方程和求解结果
这是我的代码和结果,麻烦看看哪出问题了啊,参数都没有错,我检查很多遍了
function []=zhouniu
clear;
clc;
clf;
tspan=0:0.01:100;
y=zeros(4,1);
[t,y]=ode15s(@dfunction1_2,tspan,y);

plot(y(:,1),y(:,2),'k')

xlabel('位移','FontName','Times New Roman','FontSize',12);
ylabel('速度','FontName','Times New Roman','FontSize',12);

end

function f=dfunction1_2(t,y)
f=zeros(4,1);
E=2.1*10^11;G=7.8*10^10;rou=7850;C1=1.35*10^(-8);
C2=-1.9*10^(-4);
b=0.1;
Kc=25000000;
Ca=4000;Cam=1.7;
s0=0.001;A=10^4;B=10.9;nb=1;rh=0.2;e=10^(-6);
roufl=1497.83;miufl=0.2;Kf=25*10^6;miue=0.35;miud=0.3;

Fb=67144.5;wd=3.14159;l=200;L=4000;di=0.0702056;de=0.0889;
Di=0.0714375;De=0.1651;
% 参数设置
Mb=(rou*pi*(De^2-Di^2)*l)/4;
Md=(roufl*pi*(Di^2+Cam*De^2)*l)/4;
Ma=Mb+Md+((rou*pi*(de^2-di^2)*L)/12);

I=(rou*pi*(De^4-Di^4)*l)/32+(rou*pi*(de^4-di^4)*L)/96;
Ct=(pi*miufl*l*De^3)/(2*(2*rh-De));
Ka=E*pi*(de^2-di^2)/(4*L);   
Kt=G*pi*(de^4-di^4)/(32*L);
u=-sign(y(4))*(2/pi)*atan(A*y(4))*((miue-miud)/(1+B*abs(y(4)))+miud);
ROP=C1*Fb*sqrt(wd)+C2;
c=2*pi*ROP/wd;
F0=Fb*(1+Ka/Kc);

if (y(1))>=s0*sin(nb*y(3))+e
    F=Kc*(y(1)-s0*sin(nb*(y(3))));
    T=Kc*(y(1)-s0*sin(nb*(y(3))))*((2/3)*rh*u+b*sqrt(rh*c));
elseif (y(1))<=s0*sin(nb*y(3))-e
    F=0;
    T=0;
elseif (y(1))>s0*sin(nb*y(3))-e && (y(1))<s0*sin(nb*y(3))+e
    F=(Kc*(y(1)-s0*sin(nb*y(3))+e))/2;
    T=Kc*((2/3)*rh*e*u+(b*sqrt(rh*c)*(y(1)-s0*sin(nb*y(3))+e))/2);
end

f(1)=y(2);
f(2)=(F0-F-Ka.*y(1)-Ca.*y(2))/Ma;
f(3)=y(4);
f(4)=(-T-Kt.*(y(3)-wd*t)-Ct.*y(4))/I;
end

主要方程

主要方程

文章结果

文章结果

我的结果

我的结果
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /3 下一条

快速回复 返回顶部 返回列表