双单摆摆动的模拟.pdfVIP

计算物理 课程设计说明书 题目： 双单摆摆动的模拟 学生姓名 学 号 院 （系） 理学院 专 业 应用物理 指导教师 2013 年12 月30 日 目 录 一．正文 1 1.选题背景： 1 2.方案论证： 1 3.过 (设计或实验)论述： 2 （a）.简述 2 （b）.方程分析 2 （c）．计算机模拟 3 4.结果分析： 3 5.结论或总结： 5 二．设计体会及今后的改进意见 5 参考文献资料 一．正文 1.选题背景： 双单摆在理论力学中是一个经典的模型，在对其研究中，通过运用拉格朗日力学和小 振动系统的处理方法可以近似的求解出双单摆的解析解，但是却无法让我们直观的观察双 单摆的运动情况，而借助计算机，运用龙格-库塔法可以模拟双单摆的运动情况，将解析 解转换为数字解进而进行模拟。 2.方案论证： 依据拉格朗日力学进行双单摆的求解可以有效的了解两个相互联系的单摆的运动情 况，根据 matlab 的数字模拟更加直观的展现了双单摆的运动情况，对动量守恒的研究有 着很明显的帮助。在实际生活中双单摆的运动情况比较复杂，我们不能够直观的了解双单 摆的运动情况，原因有许多，因为有许多干扰因素导致了理想的双单摆模型并不存在，而 通过计算机模拟我们可以精确的去学习和论证双单摆。 1 3.过程(设计或实验)论述： （a）.简述 如上图所示，为双单摆的基本模型，通过观察发现两个小球之间是相互关联的，要想了解 两个小球的运动情况，先要写出它们之间满足的方程，以横轴为X轴，竖轴为Y轴建立平 面直角坐标系，取θ 和θ 为两个广义坐标。 0 （b）.方程分析 由拉格朗日力学可求 d2 1 d2 d d 2 0 2 2 0   T ml  ml ml cos(  ) dt2 2 dt2 dt dt 0   Vmgl(2cos cos ) 0 LTV 带入由拉格朗日方程中去 d d  0  令y  ,y  ,y  ,y  1 0 2 dt 3 4 dt dy 1 y dt 2 dy 1 2 2 2  [ly sin(y y )cos(y y )gsiny cos(y y )ly sin(y y )2gsiny ] dt 2llcos (y y2 2 3 1 3 1 3 3 1 4 3 1 1 3 1) dy 3 y dt 4 dy 1 4  [ly sin(y y )cos(y y )2gsiny cos(y y )2ly sin(y y )2gsiny ]2 2