《第二章 初等数学模型.doc》由会员分享，可在线阅读，更多相关《第二章 初等数学模型.doc（13页珍藏版）》请在金锄头文库上搜索。
1、第三章 初等数学模型所谓初等数学模型主要是指建立模型所用的数学知识和方法主要是初等的，而不是高等的。在解决实际问题的过程中，往往主要是是看解决问题的效果和应用的结果如何，而不在于用了初等的方法还是高等的方法，对于数学建模也是这样。本章介绍了量纲分析法、比例与函数建模法，并给出了相应的一些模型。第一节 量纲分析法量纲分析提出于20世纪初，是物理学中常用的一种定性分析方法，也是在物理领域中建立数学模型的一个有力工具。它是在经验和实验的基础上, 利用物理定律的量纲齐次原则，确定各物理量之间的关系。1.1 量纲齐次原则许多物理量是有量纲的，有些物理量的量纲是基本的，另一些物理量的量纲则可以由基本量纲根
2、据其定义或某些物理定律推导出来。例如在动力学中，把长度， 质量和时间的量纲作为基本量纲，记为； 而速度的量纲可表示为. 在国际单位制中，有7个基本量：长度、质量、时间、电流、温度、光强度和物质的量，它们的量纲分别为L、M、T、I、J、和N，称为基本量纲。任一个物理量q的量纲都可以表成基本量纲的幂次之积， 量纲齐次性原则 用数学公式表示一个物理定律时，等式两端必须保持量纲一致。量纲分析就是在保证量纲一致的原则下，分析和探求物理量之间关系。先看一个具体的例子，再给出。1.2量纲分析的一般方法例1 （单摆运动）质量为的小球系在长度为的线的一端，线的另一端固定，小球偏离平衡位置后，在重力作用下做往复摆
3、动，忽略阻力，求摆动周期的表达式。解：在这个问题中有关的物理量有设它们之间有关系式 -（1.1）其中为待定常数，入为无量纲的比例系数，取（1.1）式的量纲表达式有 整理得： -（1.2）由量纲齐次原则应有 -（1.3）解得： 代入（31）得 -（1.4）（1.4）式与单摆的周期公式是一致的1.3 Buckingham Pi定理下面我们给出用于量纲分析建模的 Buckingham Pi定理，定理（Buckingham Pi定理） 设n个物理量之间存在一个函数关系 -（1.5）为基本量纲，。的量纲可表示为 矩阵称为量纲矩阵，若的秩则（1.5）式与下式等价，其中F为一个未定的函数关系，为无量纲量，且
4、可表示为 -(1.6)而为线性齐次方程组的基本解向量. 利用Pi定理建模，关键是确定与该问题相关的几个基本量纲的无量纲量。第二节 航船的阻力 长吃水深度的船以速度航行，若不考虑风的影响，航船受到的阻力除依赖于船的诸变量以外，还与水的参数密度P，粘性系数, 以及重力加速度g有关。我们利用pi定理分析f和上述物理之间的关系1 航船问题中涉及的物理量及其量纲为我们要寻求的关系式为 -（2.1）这些物理量中涉及到的基本量纲为L、M、T2 写出量纲距阵 3 解齐次线性方程组 , 可得 个基本解向量 由（1.6）式，可给出4个无量量纲 -（2.2）由Pi 定理，（2.1）等价于下列方程 -（2.3）这里是
5、未定的函数由（2.2），（2.2）可得阻力f的显式表达式， -（2.4）其中是由（2.3）可得到的未知的函数关系,在力学上, 称为Froude数，记为Fr ;称为Reynold数，记为Re , 因此（2.4）又可写为 -（2.5）4. 下面我们利用物理模拟进一步确定航船在水中的阻力。 设:分别表示模型和原型中的各物理量，由(2.5)有 当无量纲量 -（2.6）成立时 , 可得 -（2.7）则此时由模型船的阻力，及其它的；可确定原型船的阻力。下面,我们讨论一下(2.6)成立的条件,如果在实验中采用跟实际同样的水质,则 又故可得 ： -（2.8）要使得(2.8)成立 , 必有；也即模型船与原型船一
6、样大，这显然排除了物理模拟的可行性。若考虑选用不同的水质，仍设 则(2.6)式化为 -（2.7）由(2.7)可得 ，若按1：20的比例，显然无法找到如此小的粘性系数的液体。实际上的一种近似处理方法是，在一定条件下Re数的影响很小，这样可近似得到， 类似地分析，只要 即有 -（2.8）由（2.8）式就容易确定原型船的阻力第三节 抛射问题 下面我们通过一个例子，介绍如何使用无量纲化方法简化模型。抛射问题：在某星球表面以初速度竖直向上发射火箭，记星球半径为，星球表面重力加速度为，忽略阻力，讨论发射高度随时间T的变化规律。模型建立：设轴竖直向上， 时 ，火箭和星球质量分别记为和，由牛顿第二定律和万有引
7、力定律可得: -（3.1）以 代入（2.17）可得 故得如下初值问题 -（3.2）(3.2)式的解可以表为 也即发射高度是以 为参数的的函数,下面我们采用无量纲化方法化简方程（3.2）， 显然抛射问题中的基本量纲为 ， 而，所谓无量纲化是指，对(3.2)式中的和分别构造且有相同的参数组合和，使得新变量 为无量纲量，其中 ， 称为特征尺度或参考尺度；把方程（3.2）化为 对 的微分方程，即可简化模型，如何寻找特征尺度，这里我们以为例，首先写出参数的量纲距阵 t的量纲向量为 记为 : 求解线性方程组 得通解: 任取，即得到一种特征尺度，例如 得 ，得 得 同理可得的几种特征尺度等以下，我们利用不同
8、的和化简（3.2） 1. 令 ; 则由 , 代入(3.2)可得: -（3.3）(3.3)式的解可表为 ,含一个独立参数且为无量纲量.2. 令 , 类似地可将(3.2)化为 : -（3.4）3. 令 可将(3.3)化为 -（3.4）按照现有科技能力, ,在（3.4）中令，则有 -（3.5）(2.22) 的解为: ，代回原变量 得: , -（3.6）(3.6)式恰为假定火箭运动过程中所受星球引力不变的运动方程。无量纲化是用数学工具研究物理问题时常用的方法，恰当地选择特征尺度不仅可以减少参数的个数，而且可以帮助人们决定舍弃哪些次要因素。虽然量纲分析法可以得到一些重要有用的结论，但它也有很大的局限性，
9、它要用于对实际问题的定性分析。思考题： 在（3.3）和(3.4)式中,令会得到什么结果?第四节 比例与函数建模 本节介绍的几个模型,都是利用基本的比例关系与函数建立起的数学模型。4.1 动物体型问题 问题： 某生猪收购站，需要研究如何根据生猪的体长（不包括头尾）估计其体重？模型假设：1 将四足动物的躯干（不含头尾）视为质量为m的圆柱体，长度为，截面面积，直径为d，图4-1 图4-12 把圆柱体的躯干看作一根支撑在四肢上的弹性梁，动物在体重f作用下的最大下垂为，即梁的最大弯曲，根据弹性力学弯曲度理论有： -（4.1） 3 以生物进化学的角度，可认为动物的相对下垂度已达到一个最合适的数值，也即为常
10、数。模型建立： ； - -（4.2） 由（1）式 可令 为比例常数由（2）式 -（4.3） 令 由假设3，为常数 -（4.4）因此生猪的体重与体长的四次方成正比，在实际工作中，工作人员可由实际经验及统计数据找出常数K，则可近似地由生猪的体长估计它的体重。4.2 双重玻璃的功效问题： 房间居室的窗户有的是双层的，即在窗户上装两层玻璃，且中间留有一定的空隙，试比较双层玻璃窗与单层玻璃窗的热量流失？模型假设：1。设双层玻璃窗的两玻璃的厚度都为d，两玻璃的间距为L；单层玻璃窗的玻璃厚度为2d，所用玻璃材料相同，如图4-22假设窗户的封闭性能很好，两层玻璃之间的空气不流动，即忽略热量的对流，只考虑热量的传导.3. 室内温度和室外温度保持不变，热传导过程处于稳定状态，即单位时间通过单位面积的热量为常数 4. 玻璃材料均匀，热传导系数为常数. 图4-2模型建立： 对于厚度为d的均匀介质，两侧温度差为，则单位时间由温度高的一侧向温度低的一侧通过单位面积的热量Q满足 k为热传导系数 设玻璃的热传导系数为，空气的热传导系数为 （1） 先考虑单层玻璃的单位时间，单位面积的热量传导 -（4.5）（2） 考虑双层玻璃情形 此时热量先通过厚度为d的玻璃
