iypt2020年第6题下沉的碗动力学分析以及数值模拟

以下是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