第六章 代数方程模型.pptVIP

作业与实践： P213 1，2,3 6.3原子弹爆炸的能量估计与量纲分析 第六章 代数方程和差分方程模型 凯里学院数学科学学院数学建模课件 6.3 原子弹爆炸的能量估计 1945年7月16日美国科学家在新墨西哥州阿拉莫戈多沙漠试爆了全球第一颗原子弹, 震惊世界! 当时资料是保密的, 无法准确估计爆炸的威力. 英国物理学家泰勒研究了两年后美国公开的录像带, 利用数学模型估计这次爆炸释放的能量为19.2千吨. 后来公布爆炸实际释放的能量21千吨 185.0 62.0 67.3 4.61 59.0 3.26 42.8 1.36 31.9 0.66 175.0 53.0 65.6 4.34 48.7 1.93 41.0 1.22 28.8 0.52 145.0 34.0 64.3 4.07 46.9 1.79 38.9 1.08 25.4 0.38 130.0 25.0 62.9 3.80 46.0 1.65 36.3 0.94 19.9 0.24 106.5 15.0 61.1 3.53 44.4 1.50 34.2 0.80 11.1 0.10 r(m) t(ms) r(m) t(ms) r(m) t(ms) r(m) t(ms) r(m) t(ms) 泰勒测量： 时刻t 所对应的“蘑菇云”的半径r 原子弹爆炸的能量估计 爆炸产生的冲击波以爆炸点为中心呈球面向四周传播,爆炸的能量越大，在一定时刻冲击波传播得越远. 冲击波由爆炸形成的“蘑菇云”反映出来. 泰勒用量纲分析方法建立数学模型, 辅以小型试验,又利用测量数据对爆炸的能量进行估计. 物理量的量纲 长度 l 的量纲记 L=[l] 质量 m的量纲记 M=[m] 时间 t 的量纲记 T=[t] 动力学中基本量纲 L, M, T 速度 v 的量纲 [v]=LT-1 导出量纲 加速度 a 的量纲 [a]=LT-2 力 f 的量纲 [f]=LMT-2 引力常数 k 的量纲 [k] 对无量纲量?，[?]=1(=L0M0T0) 量纲齐次原则 =[f][l]2[m]-2=L3M-1T-2 在经验和实验的基础上利用物理定律的量纲齐次原则,确定各物理量之间的关系. 量纲齐次原则 等式两端的量纲一致 量纲分析~利用量纲齐次原则寻求物理量之间的关系. 例：单摆运动 l mg m 求摆动周期 t 的表达式 设物理量 t, m, l, g 之间有关系式 ?1, ?2, ?3 为待定系数，?为无量纲量 (1)的量纲表达式 与 对比 对 x,y,z的两组量测值x1,y1,z1 和x2,y2,z2, p1 = f( x1,y1,z1), p2 = f( x2, y2,z2 ) 为什么假设这种形式? 设p= f(x,y,z) x,y,z的量纲单位缩小a,b,c倍 p= f(x,y,z)的形式为 量纲齐次原则 单摆运动 单摆运动中 t, m, l, g 的一般表达式 y1~y4 为待定常数, ?为无量纲量 基本解 设 f(q1, q2, ?, qm) = 0 ys = (ys1, ys2, …,ysm)T , s = 1,2,…, m-r F(? 1, ?2,…, ?m-r ) = 0 与 f (q1, q2, ?, qm) =0 等价, F未定. Pi定理 (Buckingham) 是与量纲单位无关的物理定律，X1,X2, ?Xn 是基本量纲, n?m, q1, q2, ?qm 的量纲可表为 量纲矩阵记作 线性齐次方程组 有 m-r 个基本解，记作 为m-r 个相互独立的无量纲量, 且 则 记爆炸能量为E，将“蘑菇云”近似看成一个球形. 时刻 t 球的半径为 r t, E 空气密度ρ, 大气压强P 基本量纲：L, M, T 原子弹爆炸能量估计的量纲分析方法建模 r与哪些因素有关？ r t E ρ P L M T 量纲矩阵 y=(1,-2/5,-1/5,1/5,0) y=(0,6/5,-2/5,-3/5,1)T 原子弹爆炸能量估计的量纲分析方法建模 有2个基本解 两个无量纲量 原子弹爆炸能量估计的数值计算 时间 t 非常短 能量 E 非常大 泰勒根据一些小型爆炸试验的数据建议 用r, t 的实际数据做平均 空气密度 ?=1.25 (kg/m3) 1千吨(TNT能量) = 4.184*1012焦尔 E=19.7957 (千吨) E=8.2825×1013(焦耳) 实际值21千吨 泰勒的计算 t r 最小二乘法拟合 r=atb E=8.0276×1013 (焦耳)即19.2千吨 取y平均值得c=6.9038 模型检验 b=0.4058 ~2/5 量纲分析法的评