当前位置: 首页 > news >正文

ROCKET PROPULSION ELEMENTS——DEFINITIONS AND FUNDAMENTALS笔记

火箭推进理论包括力学、热力学、化学

单位换算

磅力 l b f lbf lbf
磅,磅质量 l b m lbm lbm l b lb lb p o u n d s pounds pounds l b s lbs lbs

1 s l u g = 1 l b f 1 f t / s 2 1 slug = \frac{1 lbf}{1 ft/s^2} 1slug=1ft/s21lbf
1 k g = 1 N 1 m / s 2 1 kg = \frac{1 N}{1 m/s^2} 1kg=1m/s21N

1 l b f = 4.45 N 1 lbf = 4.45 N 1lbf=4.45N
1 f t = 30.48 c m 1 ft = 30.48 cm 1ft=30.48cm
1 s l u g = 14.593904 k g 1 slug = 14.593904kg 1slug=14.593904kg

1 c a l = 4.184 J 1 cal = 4.184 J 1cal=4.184J
1 e r g = 1 0 − 7 J 1 erg = 10^{-7} J 1erg=107J
1 e V = 1.60218 × 1 0 − 19 J 1 eV = 1.60218\times 10^{-19} J 1eV=1.60218×1019J
1 B t u = 1055.05585 J 1 Btu = 1055.05585 J 1Btu=1055.05585J

1 l b = 0.454 k g 1 lb = 0.454kg 1lb=0.454kg
1 k g = 2.20462262 l b 1kg = 2.20462262 lb 1kg=2.20462262lb

定义

总冲(量) I t = ∫ 0 t F d t I_t=\int_{0}^{t} F dt It=0tFdt

质量比冲 I s = ∫ 0 t F d t ∫ m ˙ d t I_s = \frac{\int_{0}^{t} F dt}{\int \dot{m} dt} Is=m˙dt0tFdt,单位 m / s m/s m/s
重量比冲 I s = ∫ 0 t F d t ∫ m ˙ g d t I_s = \frac{\int_{0}^{t} F dt}{\int \dot{m} g dt} Is=m˙gdt0tFdt,单位 s s s
书中的比冲指重量比冲,对于固发,重量比冲往往根据总冲 与 发动机初始与最终重量之差 来计算
推进剂总流量 m ˙ \dot{m} m˙

等效推进剂总重量 w = I t / I s w = I_t / I_s w=It/Is
推进剂重量流量 w ˙ = m ˙ g \dot{w} = \dot{m} g w˙=m˙g

体积流量是指流经管道横截面的流体的体积与该体积通过截面所用的时间之比。$m^3/t$
质量流量是指流经管道横截面的流体的质量与该体积通过截面所用的时间之比。$kg/s$

等效排气速度 c = F / m ˙ c = F/\dot{m} c=F/m˙

推进剂比耗量为比冲倒数

质量比为最终质量与初始质量之比 M R = m f / m o MR = m_f / m_o MR=mf/mo
总的质量比为各级质量比之积

举例:三级火箭质量分别是 m 0 , m 1 , m 2 m_0, m_1, m_2 m0,m1,m2
各级质量比为: ( m 1 + m 2 ) / ( m 0 + m 1 + m 2 ) (m_1+m_2)/(m_0+m_1+m_2) (m1+m2)/(m0+m1+m2) m 2 / ( m 1 + m 2 ) m_2/(m_1+m_2) m2/(m1+m2)
总质量比为: m 2 / ( m 0 + m 1 + m 2 ) m_2/(m_0+m_1+m_2) m2/(m0+m1+m2)

推进剂质量分数 为推进剂质量与初始质量之比 ζ = ( m o − m f ) / m 0 = m p / m o \zeta = (m_o - m_f)/m_0 = m_p / m_o ζ=(momf)/m0=mp/mo
推进剂质量 m p m_p mp

推进系统死重包括:贮存推进剂、燃烧推进剂所需的硬件

冲重比为总冲与初始重量之比
推重比为推力与初始重量之比
初始重量 w 0 = ( m f + m p ) g w_0 = (m_f+m_p)g w0=(mf+mp)g

补充

重量(weight):由于引力(Gravity )而作用到物体上的力
W = m g W=mg W=mg g g g为重力(Gravity )加速度

测量 measurements
标准 standard
国际单位制 International System of Units (SI)
公制 metric system;SI只是公制中的一种
英制 British Imperial Units,美国采用英制

视重,表观重量(apparent weight):体重计的读出的重量(由于体重计与人整体作变速运动,所以体重计对人的支撑力是变化)
人静止时, F N = m g F_N = mg FN=mg
人在下蹲时,经历先加速后减速最后静止,三个阶段
(以向下为正)
向下加速,速度向下,加速度方向向下, m g − F N = m a mg - F_N = ma mgFN=ma F N < m g F_N < mg FN<mg,失重weightlessness
向下减速,速度向下,加速度方向向上, m g − F N = − m a mg - F_N = -ma mgFN=ma F N > m g F_N > mg FN>mg,超重overweight

地球上物体所受的重力只是引力的一个分量

质量有三个属性:惯性质量、主动引力质量、被动引力质量

惯性质量(inertial mass): F = m a F=ma F=ma

引力质量(gravitational mass)主要可以分为主动引力质量和被动引力质量两种,前者决定物体产生的引力场的强弱,后者决定物体在处于其他引力场时受到的引力大小
F = G M 1 M 2 r 2 F=G\frac{M_1 M_2}{r^2} F=Gr2M1M2

狭义相对论提出,物体质量随其运动速度的增加而增加
m = m 0 1 − ( v / c ) 2 m = \frac{m_0}{\sqrt{1-(v/c)^2}} m=1(v/c)2 m0,静止质量 m 0 m_0 m0
狭义相对论提出,质能方程
E = m c 2 E=mc^2 E=mc2
由上面两个方程可以得出
m 0 c 2 = E 2 − p 2 c 2 m_0 c^2 = \sqrt{E^2 - p^2 c^2} m0c2=E2p2c2 p = m v p=mv p=mv

广义相对论中的等效原理:惯性质量和引力质量是等价

推力

假设推力与质量流量恒定,气体出口速度均匀并指向轴向,总推力为:
F = m ˙ v 2 + ( p 2 − p 3 ) A 2 F= \dot{m} v_2 + (p_2 - p_3)A_2 F=m˙v2+(p2p3)A2
p 2 = p 3 p_2=p_3 p2=p3,则称发动机喷管具有最佳膨胀比
第一项为动量推力,是推进剂质量流量与排气相对飞行器的速度之积
第二项为压差推力,喷管出口面积乘以排气压力与环境压力之差

燃烧室内表面的压力 p 1 p_1 p1是最高的
喷管喉部压力 p t p_t pt
喷管出口压力 p 2 p_2 p2
大气压 p 3 p_3 p3随高度增加而降低

补充—拉瓦尔喷管

拉瓦尔喷管,渐缩渐阔喷管(de Laval nozzle,convergent-divergent nozzle,CD nozzle,con-di nozzle)

拉瓦尔喷管三段中的流体为:亚音速、音速、超音速

等熵流动方程
d p p = γ d ρ ρ \frac{dp}{p} = \gamma \frac{d \rho}{\rho} pdp=γρdρ
d p = c 2 d ρ dp = c^2 d\rho dp=c2dρ
其中, γ = C p / C v , c = γ p ρ = γ R T \gamma = C_p / C_v,c = \sqrt{\gamma \frac{p}{\rho}}=\sqrt{\gamma RT} γ=Cp/Cvc=γρp =γRT c c c为音速

动量守恒方程
u ′ = u + ( Δ v Δ x ) Δ x u' = u + (\frac{\Delta v}{\Delta x})\Delta x u=u+(ΔxΔv)Δx
p ′ = p + ( Δ p Δ x ) Δ x p' = p + (\frac{\Delta p}{\Delta x})\Delta x p=p+(ΔxΔp)Δx
A ′ = A A' = A A=A
Δ x Δ t = v \frac{\Delta x}{\Delta t} = v ΔtΔx=v
m = ρ Δ x A m = \rho \Delta x A m=ρΔxA

F = m a = m Δ v Δ t F = ma = m \frac{\Delta v}{\Delta t} F=ma=mΔtΔv
p A − p ′ A ′ = m v ′ − v Δ t pA - p'A' = m \frac{v' - v}{\Delta t} pApA=mΔtvv,将上面前三个方程代入得
p A − ( p + ( Δ p Δ x ) Δ x ) A = m ( Δ v Δ x ) Δ x Δ t pA - (p + (\frac{\Delta p}{\Delta x})\Delta x) A = m \frac{(\frac{\Delta v}{\Delta x})\Delta x}{\Delta t} pA(p+(ΔxΔp)Δx)A=mΔt(ΔxΔv)Δx,再代入后两个方程得
− Δ p Δ x = ρ v Δ v Δ x - \frac{\Delta p}{\Delta x} = \rho v \frac{\Delta v}{\Delta x} ΔxΔp=ρvΔxΔv
− d p d x = ρ v d v d x - \frac{d p}{d x} = \rho v \frac{d v}{d x} dxdp=ρvdxdv
ρ v d v = − d p \rho v dv = - dp ρvdv=dp

连续性方程(质量守恒方程)
m ˙ = ρ v A = C \dot{m} = \rho v A = C m˙=ρvA=C m ˙ \dot{m} m˙为质量流量
l n ρ + l n v + l n A = l n C ln \rho + ln v + ln A = ln C lnρ+lnv+lnA=lnC
d ( l n ρ + l n v + l n A ) = 0 d(ln \rho + ln v + ln A) = 0 d(lnρ+lnv+lnA)=0
d ρ ρ + d v v + d A A = 0 \frac{d \rho}{\rho} + \frac{d v}{v} + \frac{d A}{A} = 0 ρdρ+vdv+AdA=0

由动量守恒方程、等熵流动方程联立得
d p = c 2 d ρ , ρ v d v = − d p dp = c^2 d\rho, \rho v dv = - dp dp=c2dρ,ρvdv=dp
− ρ v d v = c 2 d ρ -\rho v dv = c^2 d\rho ρvdv=c2dρ
− ρ M 2 d v = v d ρ -\rho M^2 dv = v d\rho ρM2dv=vdρ
− M 2 d v / v = d ρ / ρ -M^2 dv / v = d \rho / \rho M2dv/v=dρ/ρ
再将上面方程代入质量守恒方程得
− M 2 d v / v = d ρ / ρ , d ρ ρ + d v v + d A A = 0 -M^2 dv / v = d \rho / \rho, \frac{d \rho}{\rho} + \frac{d v}{v} + \frac{d A}{A} = 0 M2dv/v=dρ/ρ,ρdρ+vdv+AdA=0
− M 2 d v / v + d v v + d A A = 0 -M^2 dv / v + \frac{d v}{v} + \frac{d A}{A} = 0 M2dv/v+vdv+AdA=0
( 1 − M 2 ) d v / v = − d A / A (1 - M^2) dv/v = -dA/A (1M2)dv/v=dA/A

( 1 − M 2 ) d v / v = − d A / A (1 - M^2) dv/v = -dA/A (1M2)dv/v=dA/A
如果流动为超音速( M > 1 M>1 M>1),若截面积增大( d A > 0 dA>0 dA>0),则 d v > 0 dv>0 dv>0,即速度增加
如果流动为超音速( M > 1 M>1 M>1),若截面积减小( d A < 0 dA<0 dA<0),则 d v < 0 dv<0 dv<0,即速度减小
如果流动为亚音速( M < 1 M<1 M<1),若截面积增大( d A > 0 dA>0 dA>0),则 d v < 0 dv<0 dv<0,即速度减小
如果流动为亚音速( M < 1 M<1 M<1),若截面积减小( d A < 0 dA<0 dA<0),则 d v > 0 dv>0 dv>0,即速度增加
动量守恒 ρ v d v = − d p \rho v dv = - dp ρvdv=dp,由于 ρ v > 0 \rho v > 0 ρv>0,因此 d v dv dv d p dp dp反号,即速度增加,压力减小;速度减小,压力增加

亚音速流可视为不可压缩流,密度恒定,而质量守恒 ρ v A = C \rho v A = C ρvA=C,因此面积减小时,速度增大;面积增大时,速度减小。

超音速流为可压缩流,密度变化,质量守恒 ρ v A = C \rho v A = C ρvA=C,因此面积减小时, ρ v \rho v ρv增大;面积增大时, ρ v \rho v ρv减小。
如果流动为超音速( M > 1 M>1 M>1 − M 2 d v / v = d ρ / ρ -M^2 dv / v = d \rho / \rho M2dv/v=dρ/ρ,则 − ρ / v > d ρ / d v - \rho / v > d \rho / dv ρ/v>dρ/dv 0 > d ρ / d v 0 > d \rho / dv 0>dρ/dv,因此 d v dv dv d ρ d\rho dρ反号。
如果流动为超音速( M > 1 M>1 M>1), − M 2 d v / v = d ρ / ρ -M^2 dv / v = d \rho / \rho M2dv/v=dρ/ρ,则 ∣ d ( l n ρ ) ∣ > ∣ d ( l n v ) ∣ |d(ln\rho)| > |d(lnv)| d(lnρ)>d(lnv),即速度变化量远小于密度变化量, ρ v \rho v ρv的变化主要取决于 d ρ d \rho dρ
因此,面积增大, ρ v \rho v ρv增大,则 d ρ > 0 , d v < 0 d \rho >0, dv < 0 dρ>0,dv<0;面积减小, ρ v \rho v ρv减小,则 d ρ < 0 , d v > 0 d \rho < 0, dv > 0 dρ<0,dv>0

https://www.grc.nasa.gov/WWW/k-12/airplane/conmo.html
https://www.grc.nasa.gov/WWW/k-12/airplane/isentrop.html
https://www.grc.nasa.gov/WWW/k-12/airplane/tunnozt.html

排气速度

等效排气速度 c = v 2 + ( p 2 − p 3 ) A 2 / m ˙ c = v_2 + (p_2 - p_3)A_2/\dot{m} c=v2+(p2p3)A2/m˙

特征排气速度 c ∗ = p 1 A t / m ˙ c^* = p_1 A_t / \dot{m} c=p1At/m˙

能量与效率

射流功率为单位时间内输出的喷射物质动能
对于恒定的气体喷射速度 v v v,射流功率 P j e t = m ˙ v 2 2 P_{jet} = \frac{\dot{m} v^2}{2} Pjet=2m˙v2

比功率为射流功率除以推进剂系统装填质量

化学火箭发动机输入功率 P c h e m = m ˙ Q R J P_{chem} = \dot{m} Q_R J Pchem=m˙QRJ
单位质量化学推进剂能获得的最大能量为其燃烧反应热 Q R Q_R QR

化学火箭发动机燃烧效率是单位质量推进剂的实际反应热与理论反应热之比
燃烧效率乘以 P c h e m P_{chem} Pchem即为可用功率,可用功率被转换为排气射流的动能

化学火箭发动机内效率 η i n t = P j e t / ( η c o m b P c h e m ) \eta_{int} = P_{jet}/(\eta_{comb} P_{chem}) ηint=Pjet/(ηcombPchem)

推进剂热量 = 燃烧室可用能量 + 燃烧损失(不良混合、不完全燃烧)
燃烧室可用能量 = 排气射流总能量 + 壁面热量损失
排气射流总能量 = 排气射流动能 + 排气射流不可用热能
排气射流动能 = 可用于推进飞行器的能量 + 排气残余动能

飞行器功率 P v e h i c l e = F v P_{vehicle}=Fv Pvehicle=Fv
推进效率 η p = P v e h i c l e / ( P v e h i c l e + ( w ˙ / g 0 ) ( c − u ) 2 / 2 ) \eta_p = P_{vehicle}/(P_{vehicle} + (\dot{w}/g_0)(c-u)^2/2) ηp=Pvehicle/(Pvehicle+(w˙/g0)(cu)2/2)

相关文章:

  • Jumpserver界面设置及界面功能
  • Docker12:Docker网络
  • S7-200SMART案例分析——运动控制编程三轴控制(一)
  • 四十二、路由层
  • python-(4-1)数据类型的应用(int、float、bool、字符串)
  • [Interview]Java 面试宝典系列之 Java 多线程
  • 纷享销客《2022新增长系列之快消行业橙皮书》重磅发布
  • 2023秋招面试准备
  • 【学习笔记】(数学)线性代数-矩阵的概念和特殊矩阵
  • 用ARM进行汇编语言编程(2)算数指令,CPSR寄存器与逻辑运算
  • 计算机毕业设计ssm趣评美食管理评论系统lrt3w系统+程序+源码+lw+远程部署
  • 《JavaScript从入门到精通》|变量作用域|垃圾回收|闭包【函数进阶篇】
  • 工程项目管理概述
  • 【Q-Learning】TD算法的一种
  • 【QT】Qt调用OCX控件详解
  • 网络传输文件的问题
  • 2017届校招提前批面试回顾
  • git 常用命令
  • js操作时间(持续更新)
  • React16时代,该用什么姿势写 React ?
  • storm drpc实例
  • 订阅Forge Viewer所有的事件
  • 如何设计一个微型分布式架构?
  • 找一份好的前端工作,起点很重要
  • Linux权限管理(week1_day5)--技术流ken
  • ​七周四次课(5月9日)iptables filter表案例、iptables nat表应用
  • #传输# #传输数据判断#
  • (003)SlickEdit Unity的补全
  • (03)光刻——半导体电路的绘制
  • (11)工业界推荐系统-小红书推荐场景及内部实践【粗排三塔模型】
  • (Arcgis)Python编程批量将HDF5文件转换为TIFF格式并应用地理转换和投影信息
  • (Redis使用系列) Springboot 使用redis的List数据结构实现简单的排队功能场景 九
  • (附源码)ssm考生评分系统 毕业设计 071114
  • (三) diretfbrc详解
  • (万字长文)Spring的核心知识尽揽其中
  • (转)C#调用WebService 基础
  • (转载)虚函数剖析
  • .NET Framework .NET Core与 .NET 的区别
  • .net 逐行读取大文本文件_如何使用 Java 灵活读取 Excel 内容 ?
  • .Net(C#)自定义WinForm控件之小结篇
  • .NET/C# 利用 Walterlv.WeakEvents 高性能地定义和使用弱事件
  • .NET和.COM和.CN域名区别
  • .NET牛人应该知道些什么(2):中级.NET开发人员
  • .NET序列化 serializable,反序列化
  • .NET中使用Redis (二)
  • @for /l %i in (1,1,10) do md %i 批处理自动建立目录
  • [ CTF ]【天格】战队WriteUp- 2022年第三届“网鼎杯”网络安全大赛(青龙组)
  • [ web基础篇 ] Burp Suite 爆破 Basic 认证密码
  • [2017][note]基于空间交叉相位调制的两个连续波在few layer铋Bi中的全光switch——
  • [2544]最短路 (两种算法)(HDU)
  • [android] 手机卫士黑名单功能(ListView优化)
  • [Angularjs]asp.net mvc+angularjs+web api单页应用
  • [AR Foundation] 人脸检测的流程
  • [BZOJ 1040] 骑士
  • [bzoj4240] 有趣的家庭菜园