复摆运动分析研究.pdfVIP

理论力学实验报告之复摆运动分析 姓名：桂馨 班级：机械 0803 学号： 200804000306 成绩： 复摆运动分析研究 一、实验目的 1、进一步掌握动力学基本理论，掌握复摆运动得规律； 2、掌握利用理论力学理论知识解决复杂力学问题的能力； 3、提高利用计算机进行辅助分析的能力 二、实验内容 工程中对于几何形状复杂的物体，常用实验得方法测定其转动惯量。欲求其对轴 O 的 转动惯量，可将物体在轴 O悬挂起来，并使其作微幅摆动，利用其摆动周期得等时性计算 出转动惯量。当摆动角度增大时，复摆是否还具有等时性呢？ 本实验对复摆得运动规律进行详细分析研究：建立复摆得运动微分方程，利用 Manlab 对复摆进行仿真计算，研究复摆的摆角对运动周期得影响。 三、实验原理 1、动量矩定理或刚体定轴转动微分方程 2、运动微分方程得 Matlab 数值求解 在生产和科研中，所建立的微分方程往往很复杂，且大多得不出解析解。而在实际 中一般是要求得到解在若干点上满足规定精确度的近似值。对于常微分方程 y f (t , y ) y(t0 ) y0 其数值解是指由初始点 t0 开始的若干离散的 t 值处， 即对 t0t1t2 , tn, 求出准确 值 y(to),y(t1),y(t2) , y(tn) 的近似值 y0,y1,y2 , yn 。因此研究常微分方程的数值解法是 十分必要的。 在求常微分方程数值解方面， Matlab 具有丰富的函数， 将其统称为 solver ，其一般格 式为： ［T,Y ］ =solver(odefun,tspan,) 上面的调用格式中，各参数的具体含义如下： 参数“ odefun ”表示 ODE函数的名称； 参数“ tspan ”，当 tspan 便是二元向量［ to, tn ］时， tspan 是用来定义求解数值解 的时间区间的；当 tspan 表示多元向量 [to, t1, t2 ,, tn] 时，命令将会在 tspan 定义的 时间序列进行数值求解，此时 tspan 的元素必须按照单调次序排列； 参数 y0 表示为微分方程的初始数值； 参数 T 是所求得的自变量数据列向量； 参数 Y 表示所求微分方程的因变量数据矩阵。 Solver 为命令 ode45,ode23,ode113,ode15s,ode23s,ode23t,ode23tb 之一。其中 ode45,ode23,ode113 用于求解非刚性微分方程， ode15s,ode23s,ode23t,ode23tb 用于求解 刚性微分方程。 ode45 是大部分场合的首选算法。 Matlab 中常微分方程的数值求解命令直接求解的是一阶常微分方程。 而动力学微分方 第 1 页共 4 页 理论力学实验报告之复摆运动分析 程一般是两阶的。利用数值方法求解动力学微分方程需要首先对其进行降阶增维处理。 四、实验过程 1、力学模型建立、描述 所研究复摆如下图所示： 2、数学模型建立 根据上图列出复摆的运动微分方程，其中复摆质量为 m，质心为 C，质心到悬挂点 O 的距离为 a； mga