内容正文:
第一章 解常微分方程组
第6节 小球弹簧类碰撞模型
创意图:
物理模型描述:一动碰一静;即 v1=v0;v2=0
碰撞后的速度满足 v1’=(m_1-m_2)/(m_1+m_2)*v_1
V2’=2m_1/(m_1+m_2)*v_1
步骤1:分别列出两个小球整个运动过程的微分方程
V1’(t,x1,v1,x2,v2)=k (x_2 - x_1 - l) / m_1 (sgn(l - (x_2 - x_1)) + 1) / 2
V2’(t,x1,v1,x2,v2)=-k (x_2 - x_1 - l) / m_2 (sgn(l - (x_2 - x_1)) + 1) / 2
x1’(t,x1,v1,x2,v2)=v1
X2’(t,x1,v1,x2,v2)=v2
Sgn(x)为符号函数,这里实现两小球间距离小于弹簧原长时相互作用力满足胡克定律
当两小球间距离大于弹簧原长时相互作用力为0。
然后解这个常微分方程组即可完成制作。
本例中的n与d为设置弹簧的两个参数。
学科网(北京)股份有限公司
学科网(北京)股份有限公司
$