%第二个M文件
function dx=weifenzu(t,x)
u=1; P=1;L=1;E=1;B=1;D=1652;F=1; %先假设各参量都为常数1,后期再进行修正
Re=x(4)*abs(x(5)-x(2))*x(9)/u; %雷诺数的表达式
q=(12+2.754*(x(4)*abs(x(5)-x(2))*x(9)/u)^0.55*P^0.33)*(x(7)-x(8))*L*x(1)/x(9)^2+6*x(1)*(x(7)^4-x(8)^4)*E*B/x(9);
dx=[(-f(x(1),x(4),Cd(Re),x(5),x(2),x(9))-J(x(8),x(1),Re)*x(2))/x(2)^2
f(x(1),x(4),Cd(Re),x(5),x(2),x(9))/(x(1)*x(2))
q/(x(1)*x(2))
((-0.4*D^3+2.2*x(5)*D^2+36.256*x(7)*D-10623*D-3.2*D*x(5)^2+90.64*x(5)*x(7)+10496.112*x(5)+1.4*x(5)^3)*0.046*x(4)^2+(-2.4*J(x(8),x(1),Re)*x(5)^2+1.4*J(x(8),x(1),Re)*x(2)*x(5)-0.4*J(x(8),x(1),Re)*x(3)-1.4*f(x(1),x(4),Cd(Re))*x(5)+0.4*f(x(1),x(4),Cd(Re))*x(2)+0.4*q)*x(4)+1.4*J(x(8),x(1),Re)*x(6))/(1.4*x(5)*x(6)-x(4)*x(5)^3)
(-0.018*x(4)*D^3+0.101*x(4)*x(5)*D^2-488.658*x(4)*D+1.668*x(4)*x(7)*D-0.147*x(4)*D*x(5)^2+488.658*x(4)*x(5)-1.668*x(4)*x(5)*x(7)+0.064*x(4)*x(5)^3+0.4*q-1.4*J(x(8),x(1),Re)*x(5)^2-0.4*J(x(8),x(1),Re)*x(3)+1.4*J(x(8),x(1),Re)*x(2)*x(5)+0.4*f(x(1),x(4),Cd(Re))*x(2)-1.4*f(x(1),x(4),Cd(Re))*x(5))/(x(4)*x(5)^2-1.4*x(6))
(0.018*x(5)^4-0.055*D*x(5)^3-(1.668*x(7)-0.055*D^2-488.658)*x(5)^2-488.658*D*x(5)-0.018*x(5)*D^3+1.668*x(5)*x(7)*D)*x(4)^2+(-0.4*J(x(8),x(1),Re)*x(5)^3+0.064*x(6)*x(5)^2+0.4*J(x(8),x(1),Re)*x(2)*x(5)^2-0.4*f(x(1),x(4),Cd(Re),x(5),x(2),x(9))*x(5)^2-0.129*D*x(5)*x(6)-0.4*J(x(8),x(1),Re)*x(3)*x(5)+0.4*f(x(1),x(4),Cd(Re),x(5),x(2),x(9))*x(2)*x(5)+0.4*q*x(5)+0.064*x(6)*D^2)*x(4)-1.4*J(x(8),x(1),Re)*x(5)*x(6)-1.4*f(x(1),x(4),Cd(Re),x(5),x(2),x(9))*x(6)+1.4*J(x(8),x(1),Re)*x(2)*x(6)
(0.002*x(4)^2*x(5)^5-(0.007*D*x(4)^2+0.048*J(x(8),x(1),Re)*x(4))*x(5)^4+(58.775*x(4)^2-0.265*x(7)*x(4)^2+0.007*D^2*x(4)^2-0.048*f(x(1),x(4),Cd(Re),x(5),x(2),x(9))*x(4)+0.048*J(x(8),x(1),Re)*x(2)*x(4)+0.008*x(6)*x(4))*x(5)^3+((-205.708*D+0.7*x(7)*D+0.147*x(7)*D-0.002*D^3+146.924*D-0.5*D*x(7))*x(4)^2+0.048*q*x(4)-0.048*J(x(8),x(1),Re)*x(3)*x(4)-0.015*D*x(6)*x(4)+2.4*J(x(8),x(1),Re)*x(4)*x(7)+0.048*f(x(1),x(4),Cd(Re),x(5),x(2),x(9))*x(2)*x(4)-0.168*J(x(8),x(1),Re)*x(6))*x(5)^2+(-488.658*x(7)*x(4)^2+1.668*x(4)^2*x(7)^2-0.101*x(7)*x(4)^2*D^2+0.008*x(6)*x(4)*D^2+1.4*f(x(1),x(4),Cd(Re),x(5),x(2),x(9))*x(7)*x(4)-1.4*J(x(8),x(1),Re)*x(2)*x(4)*x(7)+0.168*J(x(8),x(1),Re)*x(2)*x(6)-0.168*f(x(1),x(4),Cd(Re),x(5),x(2),x(9))*x(6))*x(5)+(1710.304*x(7)*D-5.837*D*x(7)^2+0.018*x(7)*D^3-1221.646*x(7)*D+4.169*D*x(7)^2)*x(4)^2+(0.4*J(x(8),x(1),Re)*x(3)*x(7)-0.4*q*x(7)-0.4*f(x(1),x(4),Cd(Re),x(5),x(2),x(9))*x(2)*x(7))*x(4)-1.4*J(x(8),x(1),Re)*x(6)*x(7))/(1.4*x(4)*x(5)*x(6)-x(4)^2*x(5)^3)
((x(7)-x(8))*(x(1)+x(4))^0.5)/(0.014*F^0.9*x(6)^0.5);
47.619*F^0.9+13.143*F^0.9*(x(4)*abs(x(5)-x(2))*x(9)/u)^0.5
];
end |