应物11-1变摆长单摆运动的模拟选读.doc

东 北 石 油 课 程 设 计 5年 月 日  东北石油大学课程设计任务书 课程 计算物理和MATLAB课程设计 题目 变摆长单摆运动的模拟 专业 应用物理学 姓名 王士愚 学号 110901340107 主要内容、基本要求、主要参考资料等 主要内容： 设计基于MATLAB的变摆长单摆运动的模拟。单摆是一种理想的物理模型，它由理想化的摆球和摆线组成．摆线由质量不计、不可伸缩的细线提供；摆球密度较大，而且球的半径比摆线的长度小得多，这样才可以将摆球看做质点，由摆线和摆球构成单摆。 基本要求： 2.功能明确，操作方便；可以增加其他功能或修饰。 3.代码应适当缩进，并给出必要的注释，以增强程序的可读性。 主要参考资[1] Steven E. Koonin, 秦克诚译. 计算物理学. 北京：高等教育出版社1993. [2] 马文淦等. 计算物理学. 合肥：中国科学技术大学出版社1992. [3] 张志涌. 精通MATLAB6.5. 北京：北京航空航天大学出版社，2003  目 录 第1章 绪论 1 1.1基本概念 1 1.2 设计分析 1 第2章 基本设计 2 2.1 完全理想条件下单摆的数学模型 2 2.2 现实情况下单摆的数学模型 4 2.3理想条件下特殊摆动的数学模型 5 第3章 结论 7 参考文献 8 附录 9 第1章 绪论 1.1基本概念 单摆是一种理想的物理模型，它由理想化的摆球和摆线组成．摆线由质量不计、不可伸缩的细线提供；摆球密度较大，而且球的半径比摆线的长度小得多，这样才可以将摆球看做质点，由摆线和摆球构成单摆。在满足偏角10°的条件下，单摆的周期： 1-1 用一根绝对挠性且长度不变、质量可忽略不计的线悬挂一个质点，在重力作用下在铅垂平面内作周期运动，就成为单摆。模拟单摆运动，主要是使单摆能够在最低点左右摆动动，用来在物理课的时候让同学能够更直接的感受到单摆的运行轨迹以及运行周期与初始角。 对于单摆运动，在大学物理和理论力学教材中有很多例子，主要研究了包括摆长和摆角的不同以及存在的阻尼对单摆运动的影响。 1.2 设计分析 首先根据理想条件,摆线质量忽略不计,空气阻力忽略不计。 设摆线长度为l,摆球质量为m,重力加速度为g,系统的初始时刻为t=0,在任意t≥0时刻摆球的线速度为v(t),角速度为ω(t),角位移(t),以单摆的固定位置为坐标原点建立直角坐标系,水平方向为x轴方向。示意图1-1所示： 图 1-1 在 t 时刻，摆锤所受切向力是重力mg在其运动圆弧切线方向上的分力，即 1-1第2章 基本设计 2.1 完全理想条件下单摆的数学模型 根据牛顿第二运动定律，切向加速度为： 2-1 因此得到单摆的运动微分方程组： 2-2 2-3 使用欧拉算法求解：将dv(t) ??v(t ??dt) ??v(t)和d(t) ??(t ??dt) ?(t)代入 式2-2及式2-3中，并以仿真步进量Δ作为dt 的近似，得到基于时间的递推方程： 2-4 2-5 注:本递推方程仅适合于摆角，也就是要求无论初始速度多少，摆角的最大幅度不能超过90度，如果超过90度比并且初始速度为0时放手小球会自由下落一段时间才能摆动，本递推方程无法描述。 据此编写仿真程序：（程序见附录） 在以上假设条件下得到仿真图形如图2-1： 图 2-1完全理想条件下单摆的运动模拟 在其他条件不变的情况下，仅仅改变细线的长度L=3，再次进行仿真，仿真图象如图2-2： 图 2-2完全理想条件下L=3单摆的运动模拟 对比两幅图象可以看出，在理想条件下，同样摆角下，单摆的摆臂变化，影响单摆的最大线速度以及单摆的周期，当摆臂增加时,最大线速度增加，同时单