理论力学第二次大作业论文——自由落体及有初速度条件下落体偏东问题的研究.docxVIP

自由落体及有初速度条件下落体偏东问题的研究 问题背景 落体偏东是指在北半球，当物体从高处自由下落时落点会向东偏离的现象。随着科技的飞速发展和人类对天空的不断探索，落体偏东现象成为了我们必须深入探讨的课题。早期，在国外曾做过几个落体偏东的实验（见表1）与理论值y0t=g2ω(t-sin2ωt/2ω)cosλ（或y0h=23hωcosλ2hg）进行比较。式中ω是地球自传的角速度，g是重力加速度（设与纬度无关），λ是落体所处的纬度，h是质量为m的物体由静止落下的高度，t是落体下落的时间。 如今落体偏东理论已渐趋成熟，在精确制导及载人航天等许多领域都得到了应用和拓展。本文将对此问题进行简单分析。 问题假设 落体为质点，不考虑空气阻力。 不考虑除落体和地球组成的系统外的其他外力。 地球为球体，质心与地心重合，自转速度恒定。 模型建立与求解 3.1 自由落体的偏东问题 模型建立 地球自转速度为 ω=0.72722×10-4rad/s 地球的半径为 R=6378164m 质量为m的物体位于地球赤道时的牵连惯性力FIe达到最大，即 FIe=m×6378164×ω2=0.03373m 与物体自重相比 FIemg=0.00344 可见物体的牵连惯性力远小于其重力，且与运动无关。故可忽略前因惯性力或将其合并入重力，而只考虑科氏惯性力。 以自由质点为例。以质点初始位置为坐标原点，Oz轴沿地球半径方向，Oy轴沿纬线切线向东，Ox轴沿经线切线向南。如图： G （1） (1) 由质点相对运动微分方程 设质点所在位置纬度为λ，则有 式(1)投影到x,y,z轴 (2) (3) (4) 初始条件 方程(4)中ω、y的一阶导数均为小量，与g相比，略去右端第二项 积分并考虑初始条件 z=-12gt2 方程(3)中x的一阶导数相对z的一阶导数为小量，可以略去。得到 积分并考虑初始条件 带入方程(1)，积分并考虑初始条件 在北半球sin20,物体在下落过程中有少量偏东，更少量偏南。南半球则是少量偏东，更少量偏北。 模型求解 基于以上的研究，物体在无初速度、不考虑空气阻力的条件下做自由落体运动。注意到偏南量为，其中的，可知落体的偏南量为一数量级很小的数值，考虑落体的偏南在生活中并无实际价值，所以在接下来对以上建立的模型的求解过程中，我们不考虑落体偏南，而只考虑落体偏东效应。 观察落体的偏东量表达公式，影响落体偏东量的两个因素为下落时间及物体所在纬度，而自由落体的下落时间又直接与物体所在高度相关。于是，我们得到了需要考虑的两个变量：高度及纬度。 上述求得的表达式虽然可以直观地表达落体偏东量与影响因素的函数关系，但是有一定的近似处理，在实际的求解中，为了得到更为精确的结果，我们使用了数值计算中的四五阶龙格库塔方法求解最初建立的常微分方程组，在给定的初始条件下，依据方程可以得到各个时刻的物体的状态。 首先，我们考虑同一高度下，不同纬度下的落体偏东。 如图，由于物体在南北半球的落体偏东值对称，我们选取了物体在北半球的情况进行研究。该图可以近似看作物体从高处自由落体的下落轨迹图。在赤道处，落体有最大的偏东量，随着纬度变大，落体偏东量逐渐减小，且减小的幅度越来越大。注意到，在南北两极处，物体的偏东量为0，这是因为物体下落的速度与平行，所以物体不受科氏力的作用，不发生偏东。 接下来，我们考虑不同高度，物体偏东量大小随纬度变化的曲线，用Matlab做图如下： 观察图形，我们发现，落体的偏东量随高度的增大而增大。 结果分析 总结上述过程，我们得到以下结论： 自由落体偏东量与高度和纬度的函数关系为： 自由落体偏东随纬度增大而减小，在两极处偏东量为0； 自由落体偏东随高度增大而增大。 具有x轴方向初速度的落体偏东问题 模型建立 假定质点从有限的高度处且在x轴方向上以初速度为下落，那么我们可认为不变，且除重力以外的,其初始条件变为 (5) 把式（5）代入(2)得 所以： 把式（5）代入（3）得 所以： 把式（5）代入（4）得 所以： 再将计算出来的代入(2)，(3)，（4）得 略去后得 再次积分得 消去后得 (6) 如上，我们得到了描述具有x轴方向初速度下的物体下落轨迹方程。 模型求解 为了研究有x轴方向初速度的落体偏量，我们将z=0代入到公式（6）中，得到下公式： 其中的分别为落体的偏南（偏北）量和偏东（偏西）量，注意到落体的偏南量为一数量级很小的值，在模型建立的近似处理中，我们将落体的偏南（偏北）量忽略，所以接下来我们只研究落体的偏东。 先考虑速度的变化对落体的偏东（偏西）量的影响。为简便起见，我们选择特定高度及纬度下研究速度的影响。先考虑高度为200m，纬度为时的情况，用Matlab做图如下： 观察图像，我们发现存在