数学模型(研究思路).docVIP

本数学模型中的火箭，理想化为精确地沿重力反方向竖直上升，大多数情况下考虑空气阻力。空气阻力大小与速度的关系理想化为f=kv2.飞行过程中火箭的质量变化，重心移动等均忽略不计。 较为简单的数据采集及计算： 自由落体、冲量计、功能转换计算。 如自由落体，可以通过记录火箭飞行至最高点的时间坐标t1,以及落地点的时间坐标t2,在不考虑空气阻力的情况下，可以求得 冲量计。火箭头锥部安装冲量计，测出落地瞬间冲量为Ix。忽略空气阻力， 则可以求得 .则 功能转换：让火箭在水面上起飞，落入水中，测得深度为h ,忽略空气、水的阻力。则 可以求得 以上三种方法都极为简单地忽略了空气阻力，甚至水的阻力。显然，要想更高精度地得到火箭的各种数据，这些方法都是不可取的，下面提供几种较为精确地计算方法。 加速度计装置的应用。 简易的加速度计，可以理解为如图的装置。其中转轴半径r,周期为T。 使用时在火箭内部放置该装置即可记录下加速度与时间的图像。由于火箭到达最高点后会头朝下下落，因此会有一个转体的过程，所以，得到的图像需要我们进行转录。 我们得到的原图大概会是这样的： 因为物块和火箭都受重力加速度 所以图中显示出的是发动机和空气 阻力产生的加速度，不包括重力 加速度。 我们记 。 在图中t1时刻时，an达到最大，t2时刻时推力与空气阻力相等。 在t2—t3时刻阻力大于推力，火箭减速上升。到t3速度减为0，箭体转体变为头朝下。在t3—t4时刻火箭开始做加速运动，空气阻力逐渐增大。火箭在t4时刻落地。 因为火箭在t3时刻有一转体动作，为了便于理解，我们将图像拆分为两部分，并将符号代入。则Y轴为加速度，X轴即加速度计转轴经过的路程 ，则。因此X轴即为时间轴。 再计入重力加速度：将t轴移动一个单位g.. 上图中函数t0时未经过远点是 因为没有考虑起飞前地面支持产生的 加速度。 由图可以看出， 加速度计采集到的数据主要是各加速度与施加的函数，可转化为力与时间的函数欲用此数据计算位移，虽可行但是比较烦琐。可采用三角测量法。 二．三角测量法的应用 三角测量法需要一个普通的摄象机，一把刻度尺或量角器。相机摆放距离需远大于火箭的飞行高度，使火箭飞行全过程不超出相机视野并且相机不得移动。 记火箭飞行高度为h,影象中于标尺上测得的高度为h’，仰角θ。 相机为n帧每秒即fps=n。刻度尺距离相机x,距离火箭起点999x。 则h=1000h’=1000xtanθ 若在一幅图上，即一帧的画面中，火箭位移在刻度尺上为△h’, 则瞬时速度 将瞬时速度近似看作该帧中间时刻速度， 则瞬时加速度 这种方法更直观地获取了高度数据，也可借此画出 v-t,a-t图象。 多普勒效应的应用 要应用多普勒效应，就得先确定一个声源。由于火箭发动机本身频率范围过大，而且频率会改变，经过多次试验后仍无法判断频率变化，故放弃以发动机声音作为声源的方案。 经过考虑，我们认为在火箭上安装一个频率稳定，频率范围在发动机声音频率范围之为的声源的方案具有一定的可行性，故在此进行详细讨论。 我们记测得的频率为δ’,记声源本身频率为δ。V为声源移动速度，v源 为声音传播速度，可以取344m/s.记声音采样频率为n,可以取44.1kHz根据多普勒效应公式 将录音机置于火箭正下方，可以推得 根据此公式，可以绘出v随时间变化的函数。 则函数曲线下半部分面积即为火箭垂直高度。 另， 因为采样频率为n,则最高精度的瞬时加速度 若n取44.1kHz，则求出的加速度曲线精确度是十分可观的。 因为，利用多普勒效应，我们也得到了v-t , a-t 的图象。 四．求空气阻力公式中k的值（f=kv2） 由于没有优良的设备条件，比如风洞等。我们无法得知火箭在流体中的受力情况。 在此提供一种方法，准确度虽不高，且f=kv2的公式较为理想化，但思路还是有借鉴价值的，在这里提供，仅仅作为参考。 在加速度计得出的图象中，下落部分所测得加速度，实为火箭壳体空气阻力造成的（加速度计在火箭内部，几乎无空气阻力）。从an-t 图我们可以得出f-t图象 将函数记为f(t). 根据公式f=kv2. 再根据 a合-t 图，我们可以描绘出v随t的变化函数即V（t）。 将V（t）代入公式，可得出，k 将为一个常数。 至此，我们求得了k的值，从而可以进行更复杂的计算。 火箭发动机总冲量及比冲量的计算。 根据多普勒效应的数据，我们已经能够描绘出火箭上升阶段