iypt2020年第6題下沉的碗動力學分析以及數(shù)值模擬
以下是matlab源代碼：
global M g p R r
M=0.5;%千克
g=9.8;
p=1000;%水的密度
R=0.03;
r=0.005;
t1=[0,75];
H1=[0;0;0];
[t,H] = ode45(@dH,t1,H1);
plot(t,H);
function dHdt=dH(t,H)
global M g p R r
dHdt=zeros(3,1);
dHdt(1)=H(2);
dHdt(2)=M*g-p*g*pi*(H(1)-H(3))*(R^2-r^2);
dHdt(3)=sqrt(-4*g*r^4/R^4*(H(3)-H(1)));
end
標簽：
