课程设计说明书
院(系)名 称: 宇 航 学 院 学 生 姓 名: 李 东 来 学 号:12151075
专 业 名 称: 飞行器动力工程(航天) 指 导 教 师:张黎辉
2016.1.22
课程设计任务书
一、课程设计题目:
设计实验用液体火箭发动机推力室
二、课程设计题目的原始数据及设计技术要求
推力:500N
燃料:气氧+75%酒精 余氧系数:α=0.8 燃烧室压力:2MPa 出口压力:0.1MPa
三、课程设计任务:
1进行热力计算、推力室结构参数计算:确定圆柱形燃烧室直径、 长度,喉部直径,喷管收敛段、扩张段长度,喷管出口直径。 2进行喷嘴设计、推力室水冷却计算。 3 详细设计并绘制推力室部件总图。 4 零件设计: 5 撰写设计说明书。
四、课程设计日期:自2015年12月14日至2016年1月22日 学生:李东来指导教师:张黎辉 班级:121516教研室主任:
目录
1. 设计参数 ............................................................................................................ 1 2. 推力室参数计算结果 ........................................................................................ 1 3. 推力室结构参数计算 ........................................................................................ 1 4. 推力室头部设计 ................................................................................................ 3 4.1 燃料喷嘴设计 ................................................................................................ 3 4.2 氧化剂喷嘴: ................................................................................................ 3 5. 推力室身部设计 ................................................................................................ 4
5.1 推力室圆筒段冷却计算 ........................................................................ 4
5.1.1 燃气的气动参数 ......................................................................... 4 5.1.2 计算燃气与内壁面的对流换热密度 ......................................... 4 5.1.3 计算燃气与内壁面的辐射热流密度 ......................................... 5 5.1.4 计算总热流密度、总热流量及冷却剂流量 ............................. 6 5.1.5 确定冷却通道参数 ..................................................................... 6 5.1.6 计算内壁面和外壁面温度 ......................................................... 6 5.2 推力室喉部冷却计算 ............................................................................ 7
5.2.1 燃气的气动参数 ......................................................................... 7 5.2.2 计算燃气与内壁面的对流换热密度 ......................................... 7 5.2.3 计算燃气与内壁面的辐射热流密度 ......................................... 8 5.2.4 计算总热流密度、总热流量及冷却剂流量 ............................. 8 5.2.5 确定冷却通道参数 ..................................................................... 9 5.2.6 计算内壁面和外壁面温度 ......................................................... 9
6. 推力室强度校核 .............................................................................................. 10
6.1推力室圆筒段强度校核 ....................................................................... 10 6.2喷管强度校核 ....................................................................................... 10 7. 点火器设计 ...................................................................................................... 11 8. 螺栓强度校核 .................................................................................................. 12 9. 整体结构分析 .................................................................................................. 12
9.1头部结构 ............................................................................................... 12 9.2喷注器 ................................................................................................... 13 9.3点火器 ................................................................................................... 13 9.4推力室 ................................................................................................... 13 9.6密封结构 ............................................................................................... 13 10. 感悟 ................................................................................................................ 14 参考文献 ............................................................................................................. 15
1. 设计参数
推力:Ftc=500N 推进剂:气氧+75%酒精 余氧系数:α=0.8 燃烧室压力:pc=2MPa 出口压力:pe=1atm
2. 推力室参数计算结果
化学当量混合比rm0=1.485,实际混合比rmc=1.188,地面理论比冲Istct =2317.87m/s,特征速度c∗=1649.35m/s,扩张比εe=3.635。
3. 推力室结构参数计算
取燃烧室效率ηC=0.96;喷管效率ηn=0.96。 推力室总质量流量为:qmc=从而得出:
推力室氧化剂质量流量:qmoc=qmc×r
rmc+1
mc
FtcIstct ηcηn
=0.234kg/s
=0.127kg/s
推力室燃料质量流量:qmfc=qmc−qmoc=0.107kg/s 喷管喉部面积:At=
c∗⋅qmcpc
=1.930×10−4m2
喉部直径Dt=
4Atπ
=15.68mm,圆整取Dt=16mm ,则At=2.01×10−4m2
取燃烧室的特征长度L=2.4m
−43
燃烧室容积Vc=L×At=4.8255×10m
利用燃烧室收缩比求燃烧室直径
取燃烧室收缩比为εc=16则燃烧室横截面积为:Ac=εcAt=3.217 ×10−3m2燃烧室直径为:Dc= c⋅Dt=64mm 设计推力室喷管双圆弧收敛段型面
选择R1=1.5Rt=11.76mm,圆整取R1=12mm 选择R 2,取ρ=2.5,则R2=ρRc=ρ cRt=80mm 则收敛段长度为:Lc2
=Rt k+ρ c − ρ−1 c+k+1
2
2
式中,k=1.5,ρ=2.5,εc=16,则计算得Lc2=61.96mm。圆整取Lc2=62mm
以R1和R2所作圆弧切点的位置为 =c2=8.09mm
H=Lc2− =53.91mm
2
y=kRt+Rt− k2Rt− 2=11.14mm
收敛段容积为:Vc2=1.17982385×10−4m3 燃烧室圆柱段的长度为:Lc1=
Vc−Vc2Ac
=113.33mm,圆整取114mm
喷管扩张段直径De= e⋅Dt=30.505mm,圆整取32mm 喷管扩张段与喉部截面之间用半径R3=Dt=16mm的圆弧过渡
n=1.7930,求大圆弧相对半径 喷管出口角取2βe=150查得喷管相对长度L
0=R
De
2Ln+ 1.5−2D−1
t
2
nsinβe− 1.5− cosβe 2 1−L
2D
t
e
=4.56
0Dt=72.97mm,圆整取73mm 则:R0=R
nDt=28.688mm,圆整取29mm Ln=L
X0=Ln+R0sinβe=38.528mm,圆整取39mm
Y0=R0cosβe−βm=sin −1
De2
=56.375mm,圆整取57mm =22.46°。
Ln+R0sinβe
R0+Dt
4. 推力室头部设计
采用带切向孔的直流-离心式喷嘴,燃料采用切向式离心喷嘴,氧化剂采用直流式喷嘴;排布方式:中间1个喷嘴,外圈均布3个喷嘴,燃料和氧化剂喷嘴数量为nf=n0=4。
4.1燃料喷嘴设计
已知:75%酒精密度:ρf=877.3kg∕m3 酒精喷嘴压降∆p0=0.2pc=0.4MPa
根据经验数据确定流量系数。取l d=3,得到μ=0.8 则有
Anf==5.049×10−6m2
df= n
nf
fπ
=1.035mm,圆整取1.1mm
lf=3df=3.3mm
该喷嘴为自击式,故αs始终为0,取喷嘴偏转角度αf=45°。
4.2氧化剂喷嘴:
氧化剂喷嘴质量流量
qmo =qmoc nf=
0.127
kg/s=0.042kg/s 选取喷嘴压降Δpo=0.4MPa,取流量系数μo=0.8。 气氧压力
pin=pc+Δp=2.4MPa
标准状况下,氧气多变指数k=1.4,ρst=1.43kg∕m3。由状态方程
ppst
= ρ得
st
ρ
k
pin1∕k
ρo= ρst=12.17kg∕m3
st
气氧的喷出速度为
2kpc
k−1
2×1.42
1.4−1
wo=
RTin 1− in
=
×259.8×273.15× 1− =52.94m/s
由气体直流喷嘴的质量流量方程得
A=qno
mo o=0.032×106mm2=81.49mm2oo 喷孔直径d4A
o= non0
π=5.88mm,圆整取do=6mm,lo=18mm
由几何关系可知壁厚b=
2
lf+d =3.1mm,圆整取b=4mm
5. 推力室身部设计
5.1推力室圆筒段冷却计算
5.1.1燃气的气动参数
圆筒段燃气温度T st =3006K 燃气多变指数k=1.15
燃气定压比热容cp=2.296kJ kg∙K 燃气粘度μ=1.0×10−4Pa⋅s 燃气普朗特数Pr=0.5881
5.1.2计算燃气与内壁面的对流换热密度
圆筒段横截面积
12A=πDc=3.217×10−3m2
喷管喉部过渡平均半径
R =0.5 R 1+R 3 =14mm
假设内壁温度:T wg =500K 。
利用巴兹法计算燃气与内壁面的对流换热系数: 根据T
T wg
st
=0.166,查表得到考虑附面层内燃气性能变化的修正系数σ=1.52。
燃气与内壁面的对流换热系数
0.026μ0.2cppc0.8Dt0.1At0.9
g= ς=1067W/(m2⋅K)
Dt
ns燃气与内壁面的对流换热密度qk= g Tst−Twg =2.675×106W∕m2
5.1.3计算燃气与内壁面的辐射热流密度
c1
根据D=1.78,查得气体对整个壁面辐射的平均射线长
c
L
l=0.85Dc=54.4mm
水蒸气分压pH2O=pc⋅
nHOnn
=1.12MPa =0.34MPa
二氧化碳分压pCO2=pc⋅
nCO计算得到pH2Ol=0.061MPa⋅mpCO2l=0.018MPa⋅m
查图得水蒸气发射率ε0H 2O =0.014,指数关系n=1+kH2OpH2O=1.6,则水蒸气的实际发射率为
εH2O=1− 1−ε0H2O =1− 1−0.014 1.6=0.0223
查图得二氧化碳发射率εC O2=0,则总的发射率为
εg=εH2O+εCO2−εH2OεCO2=0.0223
壁面发射率一般取为ε0w =0.8,则实际有效壁面发射率为
εW=ε0W 1+ 1−ε0W 1−εg =0.956
由于壁面温度较低,故壁面对燃气的辐射可以忽略,因此燃气辐射热流密度
n
为
Tg4Twg4Tg4
qr=5.67⋅εw⋅ εg −αg =5.67⋅εw⋅εg =9.870×104W∕m2
5.1.4计算总热流密度、总热流量及冷却剂流量
总热流密度q=qk+qr=2.774×106W∕m2
取推力室圆筒段壁厚δ=1mm,则外壁面温度 TWf=TWg−总热流量Φ=q⋅2πRclc=6.36×104W
若要求冷却水通过冷却通道时的温升为40K ,则冷却水流量为
Φ
q co ==0.381kg/s
pco co 冷却水的温度可以取为
T co =
20+602
qδλw
=492K
+273 K =313K
5.1.5确定冷却通道参数
推力室内壁面及肋条材料为1Cr18Ni9Ti ,导热系数λW =23.26W ∕ m ⋅K ,内壁厚δ=1mm ,冷却通道高度h=3mm,肋条厚b=1.5mm,冷却通道宽a=2mm,则当量直径ⅆe=a+b=2.4mm,冷却通道数n=57。
2a
5.1.6计算内壁面和外壁面温度
冷却通道面积 冷却剂流速
e=
A=na =3.42×10−4m2 Vt=
μco
qcoA⋅ρco
=1.114m/s
1000×1.25×2.4×10−3
0.658×10冷却剂雷诺数 冷却剂普朗特数 冷却剂努塞尔数
ρcoVtde
==4063
Pr=
Cpcoμcoλco
=4.31
0.4
Nu f =0.023R 0.8=31.8 e Pr
冷却剂和外壁面的对流换热系数
f=
Nufλcode
=8321W∕m2
液体壁面温度 气体壁面温度
TWf=η
q
p f
+Tf=480K
Twg=λ+Twf=488K
w
qδ
由计算结果可知,推力室圆筒段气体壁面Twg=484K小于假定的温度500K, 相差2.4%小于5%,符合冷却要求。
5.2推力室喉部冷却计算
5.2.1燃气的气动参数
喉部燃气温度 燃气多变指数 燃气定压比热容 燃气粘度
T st =2845K 。 k=1.1478。 C p =2.28kJ ∕kg ⋅K μ=0.0001Pa ⋅s Pr=0.5967 pt=1.1552MPa
燃气普朗特数 喉部压力
5.2.2计算燃气与内壁面的对流换热密度
喉部横截面积
A=4πDt2=2.01×10−4m2
1
喷管喉部过渡半径 R =0.5 R 1+R 3 =14mm 假设内壁温度:
T wg =740K 。
利用巴兹法计算燃气与内壁面的对流换热系数 根据T
T wg
st
=0.26,查表得到考虑附面层内燃气性能变化的修正系数σ=1.35。
燃气与内壁面的对流换热系数
0.026μ0.2Cppc0.8Dt0.1At0.9
g= ς=11320W/(m2⋅K)
t
ns燃气与内壁面的对流换热密度
qk= g Tst−Twg =2.38×107W∕m2
5.2.3计算燃气与内壁面的辐射热流密度
c1
根据D=1.78,查得气体对整个壁面辐射的平均射线长
c
L
l=0.85Dc=54.4mm
水蒸气分压
pH2O=pt⋅pCO2=pt⋅
nHOnnCOn
=0.670MPa =0.208MPa
pCO2l=0.011MPa⋅m
二氧化碳分压 计算得到
pH2Ol=0.036MPa⋅m
查图得水蒸气发射率ε0H 2O =0.01,指数关系n=1+kH2OpH2O=1.4,则水蒸气的实际发射率为
εH2O=1− 1−ε0H2O =0.014
查图得二氧化碳发射率εC O2=0,则总的发射率为
εg=εH2O+εCO2−εH2OεCO2=0.014
壁面发射率一般取为ε0w =0.8,则实际有效壁面发射率为
εW=ε0W 1+ 1−ε0W 1−εg =0.958
由于壁面温度较低,故壁面对燃气的辐射可以忽略,因此燃气辐射热流密度为
Tg4Twg4Tg4qr=5.67⋅εw⋅ εg −αg =5.67⋅εw⋅εg =49820W∕m2
n
5.2.4计算总热流密度、总热流量及冷却剂流量
总热流密度
q=qk+qr=2.39×107W∕m2
取推力室圆喉部壁厚δ=1mm,则外壁面温度
TWf=TWg
qδ2.39×107×1×10−3−= 740−K=672K w取喉部截面长度10mm 总热流量
Φ=q⋅S=1.15×104W
Φ1.15×104ΔT ===6.5K
co co 取定性温度313K
冷却水流量为q co =0.421kg/s ,则喉部内冷却水温升:
5.2.5确定冷却通道参数
推力室内壁面及肋条材料为1Cr18Ni9Ti ,导热系数λW =23.26W ∕ m ⋅K ,内壁厚δ=1mm ,冷却通道高度h=3mm,肋条厚b=1.5mm,冷却通道宽a=2mm,则当量直径ⅆe=a+b=2.4mm,冷却通道数n=14。
2a
5.2.6计算内壁面和外壁面温度
冷却通道面积 冷却剂流速
A=14a =8.4×10−5m2
co
Vt=A⋅ρ=5.01m/s
co
q
冷却剂雷诺数 冷却剂普朗特数 冷却剂努塞尔数
Re=Pr=
ρcoVtdeμcoCpcoμcoλco
=18357 =4.31
0.4Nu f =0.023R 0.8=106 e Pr
冷却剂和外壁面的对流换热系数
Nufλco
f==28223W∕m2
e
液体壁面温度 气体壁面温度
TWf=η
q
p f
+Tf=683K
Twg=λ+Twf=727K
w
qδ
由计算结果可知,推力室喉部气体壁面Twg=727K大于假定的温度740K, 相差1.7%小于5%,,且当Twg=727K时1Cr18Ni9Ti 不会失效,故符合冷却要求。
6. 推力室强度校核
6.1推力室圆筒段强度校核
圆筒段所选材料为1Cr18Ni9Ti ,查机械手册得:抗拉强度为ςb=550MPa,屈服极限为ςs=200MPa
将推力室圆筒段视为承压薄壁圆筒,其周向应力为
pcDc
ς1=ςt==32MPa
1轴向应力为
ς2=ςx=
径向应力为
ς3=ςr=pc=2 MPa
推力室圆筒段内壁温度约400K,温度较高,材料的屈服极限有所下降,故采用较大的安全系数,取n=2.5
许用应力为
ςs =
由第三强度理论有
ςr3=ς1−ς3=30MPa
由第四强度理论有
ςr4
(ς1−ς2)+(ς2−ς3)+(ς3−ς1)= =30.3MPa
2
2
2
F
=1.22MPa ςs
=80MPa 由此推力室圆柱段满足强度条件
6.2喷管强度校核
喷管所选材料为1Cr18Ni9Ti ,查机械手册得:抗拉强度为ςb=550MPa,屈服极限为ςs=200MPa
喷管喉部周向应力
ptDt
ς1=ςt==4.58MPa
1
喷管喉部轴向应力为
ς2=ςx=
喷管喉部径向应力为
ς3=ςr=pt=1.207MPa
F514==4.72MPa 推力室圆筒段内壁温度约640K,温度更高,材料的屈服极限下降更多,故采用更大的安全系数,取n=5
许用应力为
ςs =
由第三强度理论有
ςr3=ς1−ς3=3.483MPa
由第四强度理论有
(ς1−ς2)+(ς2−ς3)+(ς3−ς1)
ςr4= =3.41MPa
2
2
2
ςs
=40MPa 由此推力室喷管段满足强度条件
综上,推力室各段在保证发动机性能和冷却要求的前提下均满足强度条件。
7. 点火器设计
点火器采用棒料加工而成,在棒料的一侧边钻孔,用来安装酒精管嘴和电点火塞,气氧管嘴安装在上部平面。为保证点火腔有一定的压力,应设计燃气通道喉部,具体计算如下:
At(点火)=
C∗∙q点火点火
=2.831×10−5m2
dt(点火)=6mm
8. 螺栓强度校核
燃烧室压力为2MPa,酒精腔压力为2.4MPa,气氧腔压力为2.5MPa,喷管出口压强为1atm 。本地面试验用气氧酒精发动机采用强度等级为4.6级受拉螺栓,螺栓材料为钢,其强度极限为
ςb=400MPa
则其屈服极限为
ςs=0.6×400=240MPa
取安全系数为1.5,则许用屈服极限为
σ =
240
=160MPa 综合螺栓受力,得到与下底和法兰连接的螺栓所受的力为
p∙Sσ==102.9MPa
上底与中底采用螺钉连接,校核如下
σ=
p∙S
=14.8MPa
9. 整体结构分析
为便于设计加工,结构上多采用焊接进行连接。该发动机是单件生产,因此大部分零件均采用机械加工,零件尽量采用标准件。发动机燃烧室燃气温度较高,对材料耐高温性能均有很高的要求。
9.1头部结构
该发动机采用三底两腔的经典结构。上底和中底通过螺钉连接。该发动机仅用于地面,故可以采用简单且可重复使用的电点火方式。
气氧和酒精分别通过不同的管路进入两个腔体,并从三个直流撞击式喷嘴喷入燃烧室燃烧。
9.2喷注器
本发动机是气氧酒精双组元发动机,采用直流撞击式双组元喷注器。作为燃料的酒精通过自击式喷嘴喷注进入燃烧室,这样液滴细度高好而且混合距离短。作为氧化剂的气氧,从结构简单原则敖律,让其通过中间的直流喷嘴喷出,与酒精混合即可。
喷注盘面较小,故喷嘴数量不宜过多,可选择三个喷注器并均匀分布在点火器出口的圆周上。
9.3点火器
该发动机用于地面试验,故可以使用电点火,结构简单,可靠性高,可重复使用,点火组元仍然采用气氧和酒精,使其燃烧产物同推力室一致,预燃室通过螺纹与上底连接。点火器的燃气通过其喉部进入燃烧室点燃主推进剂。
9.4推力室
燃烧室为圆筒形,圆筒段采用的材料为1Cr18Ni9Ti 。
圆筒段采用的冷却方式为外冷却,冷却通道为铣槽结构,有56个冷却通道。 该发动机为试验用,其收缩段采用双圆弧型面,扩张段也采用双圆弧型面。与燃烧室为一体,冷却方式为外冷却,喉部有18个冷却通道,。
喷管冷却套采用两个对称的版套筒,用焊接的方式连接而成。
9.6密封结构
该发动机中需要考虑的密封较多,为简单起见,在没有高温燃气处的密封均采用国家标准液压气动用o 型橡胶密封圈,在接触高温燃气处采用矩形紫铜密封圈。
10. 感悟
经过此次课程设计,我对于火箭发动机设计方面的理论知识有了更深的理解,同时也对于机械设计的流程有了进一步的掌握。在本次设计中,张老师曾多次向我提出某零件如何安装,如何定位,如何确保密封的提问,让我意识到自己在多方面考虑火箭发动机设计问题方面还有着很多的不足。而此次的课程设计也正是为我提供了这样一个实践的机会,来了解设计人员在设计火箭发动机时所要考虑的多方面问题。不仅如此,此次火箭发动机课程设计的重头戏——手工绘制设计图也同样让我感受到自己在工程制图方面的不足,同时也为我提供了锻炼的机会。总之,在此次课程设计中,我的收获很大,从理论计算到实际制图的能力都有了不小的提升。感谢学院安排的此次课程,感谢张黎辉老师的悉心教导。
参考文献
[1] 邢继发,刘国球,等. 世界导弹与航天发动机大全. 北京:军事科学出版社,1999
[2] 蔡国飙,李家文,等. 液体火箭发动机设计. 北京:北京航空航天大学出版社,2011
[3] 王之栎,王大康,等. 机械设计综合课程设计. 北京:北京航空航天大学出版社,2008
[4]D.K.休泽尔等. 液体火箭发动机现代工程设计. 北京:宇航出版社,2004 [5] 杨世铭,陶文铨,等. 传热学. 北京:高等教育出版社,1998
课程设计说明书
院(系)名 称: 宇 航 学 院 学 生 姓 名: 李 东 来 学 号:12151075
专 业 名 称: 飞行器动力工程(航天) 指 导 教 师:张黎辉
2016.1.22
课程设计任务书
一、课程设计题目:
设计实验用液体火箭发动机推力室
二、课程设计题目的原始数据及设计技术要求
推力:500N
燃料:气氧+75%酒精 余氧系数:α=0.8 燃烧室压力:2MPa 出口压力:0.1MPa
三、课程设计任务:
1进行热力计算、推力室结构参数计算:确定圆柱形燃烧室直径、 长度,喉部直径,喷管收敛段、扩张段长度,喷管出口直径。 2进行喷嘴设计、推力室水冷却计算。 3 详细设计并绘制推力室部件总图。 4 零件设计: 5 撰写设计说明书。
四、课程设计日期:自2015年12月14日至2016年1月22日 学生:李东来指导教师:张黎辉 班级:121516教研室主任:
目录
1. 设计参数 ............................................................................................................ 1 2. 推力室参数计算结果 ........................................................................................ 1 3. 推力室结构参数计算 ........................................................................................ 1 4. 推力室头部设计 ................................................................................................ 3 4.1 燃料喷嘴设计 ................................................................................................ 3 4.2 氧化剂喷嘴: ................................................................................................ 3 5. 推力室身部设计 ................................................................................................ 4
5.1 推力室圆筒段冷却计算 ........................................................................ 4
5.1.1 燃气的气动参数 ......................................................................... 4 5.1.2 计算燃气与内壁面的对流换热密度 ......................................... 4 5.1.3 计算燃气与内壁面的辐射热流密度 ......................................... 5 5.1.4 计算总热流密度、总热流量及冷却剂流量 ............................. 6 5.1.5 确定冷却通道参数 ..................................................................... 6 5.1.6 计算内壁面和外壁面温度 ......................................................... 6 5.2 推力室喉部冷却计算 ............................................................................ 7
5.2.1 燃气的气动参数 ......................................................................... 7 5.2.2 计算燃气与内壁面的对流换热密度 ......................................... 7 5.2.3 计算燃气与内壁面的辐射热流密度 ......................................... 8 5.2.4 计算总热流密度、总热流量及冷却剂流量 ............................. 8 5.2.5 确定冷却通道参数 ..................................................................... 9 5.2.6 计算内壁面和外壁面温度 ......................................................... 9
6. 推力室强度校核 .............................................................................................. 10
6.1推力室圆筒段强度校核 ....................................................................... 10 6.2喷管强度校核 ....................................................................................... 10 7. 点火器设计 ...................................................................................................... 11 8. 螺栓强度校核 .................................................................................................. 12 9. 整体结构分析 .................................................................................................. 12
9.1头部结构 ............................................................................................... 12 9.2喷注器 ................................................................................................... 13 9.3点火器 ................................................................................................... 13 9.4推力室 ................................................................................................... 13 9.6密封结构 ............................................................................................... 13 10. 感悟 ................................................................................................................ 14 参考文献 ............................................................................................................. 15
1. 设计参数
推力:Ftc=500N 推进剂:气氧+75%酒精 余氧系数:α=0.8 燃烧室压力:pc=2MPa 出口压力:pe=1atm
2. 推力室参数计算结果
化学当量混合比rm0=1.485,实际混合比rmc=1.188,地面理论比冲Istct =2317.87m/s,特征速度c∗=1649.35m/s,扩张比εe=3.635。
3. 推力室结构参数计算
取燃烧室效率ηC=0.96;喷管效率ηn=0.96。 推力室总质量流量为:qmc=从而得出:
推力室氧化剂质量流量:qmoc=qmc×r
rmc+1
mc
FtcIstct ηcηn
=0.234kg/s
=0.127kg/s
推力室燃料质量流量:qmfc=qmc−qmoc=0.107kg/s 喷管喉部面积:At=
c∗⋅qmcpc
=1.930×10−4m2
喉部直径Dt=
4Atπ
=15.68mm,圆整取Dt=16mm ,则At=2.01×10−4m2
取燃烧室的特征长度L=2.4m
−43
燃烧室容积Vc=L×At=4.8255×10m
利用燃烧室收缩比求燃烧室直径
取燃烧室收缩比为εc=16则燃烧室横截面积为:Ac=εcAt=3.217 ×10−3m2燃烧室直径为:Dc= c⋅Dt=64mm 设计推力室喷管双圆弧收敛段型面
选择R1=1.5Rt=11.76mm,圆整取R1=12mm 选择R 2,取ρ=2.5,则R2=ρRc=ρ cRt=80mm 则收敛段长度为:Lc2
=Rt k+ρ c − ρ−1 c+k+1
2
2
式中,k=1.5,ρ=2.5,εc=16,则计算得Lc2=61.96mm。圆整取Lc2=62mm
以R1和R2所作圆弧切点的位置为 =c2=8.09mm
H=Lc2− =53.91mm
2
y=kRt+Rt− k2Rt− 2=11.14mm
收敛段容积为:Vc2=1.17982385×10−4m3 燃烧室圆柱段的长度为:Lc1=
Vc−Vc2Ac
=113.33mm,圆整取114mm
喷管扩张段直径De= e⋅Dt=30.505mm,圆整取32mm 喷管扩张段与喉部截面之间用半径R3=Dt=16mm的圆弧过渡
n=1.7930,求大圆弧相对半径 喷管出口角取2βe=150查得喷管相对长度L
0=R
De
2Ln+ 1.5−2D−1
t
2
nsinβe− 1.5− cosβe 2 1−L
2D
t
e
=4.56
0Dt=72.97mm,圆整取73mm 则:R0=R
nDt=28.688mm,圆整取29mm Ln=L
X0=Ln+R0sinβe=38.528mm,圆整取39mm
Y0=R0cosβe−βm=sin −1
De2
=56.375mm,圆整取57mm =22.46°。
Ln+R0sinβe
R0+Dt
4. 推力室头部设计
采用带切向孔的直流-离心式喷嘴,燃料采用切向式离心喷嘴,氧化剂采用直流式喷嘴;排布方式:中间1个喷嘴,外圈均布3个喷嘴,燃料和氧化剂喷嘴数量为nf=n0=4。
4.1燃料喷嘴设计
已知:75%酒精密度:ρf=877.3kg∕m3 酒精喷嘴压降∆p0=0.2pc=0.4MPa
根据经验数据确定流量系数。取l d=3,得到μ=0.8 则有
Anf==5.049×10−6m2
df= n
nf
fπ
=1.035mm,圆整取1.1mm
lf=3df=3.3mm
该喷嘴为自击式,故αs始终为0,取喷嘴偏转角度αf=45°。
4.2氧化剂喷嘴:
氧化剂喷嘴质量流量
qmo =qmoc nf=
0.127
kg/s=0.042kg/s 选取喷嘴压降Δpo=0.4MPa,取流量系数μo=0.8。 气氧压力
pin=pc+Δp=2.4MPa
标准状况下,氧气多变指数k=1.4,ρst=1.43kg∕m3。由状态方程
ppst
= ρ得
st
ρ
k
pin1∕k
ρo= ρst=12.17kg∕m3
st
气氧的喷出速度为
2kpc
k−1
2×1.42
1.4−1
wo=
RTin 1− in
=
×259.8×273.15× 1− =52.94m/s
由气体直流喷嘴的质量流量方程得
A=qno
mo o=0.032×106mm2=81.49mm2oo 喷孔直径d4A
o= non0
π=5.88mm,圆整取do=6mm,lo=18mm
由几何关系可知壁厚b=
2
lf+d =3.1mm,圆整取b=4mm
5. 推力室身部设计
5.1推力室圆筒段冷却计算
5.1.1燃气的气动参数
圆筒段燃气温度T st =3006K 燃气多变指数k=1.15
燃气定压比热容cp=2.296kJ kg∙K 燃气粘度μ=1.0×10−4Pa⋅s 燃气普朗特数Pr=0.5881
5.1.2计算燃气与内壁面的对流换热密度
圆筒段横截面积
12A=πDc=3.217×10−3m2
喷管喉部过渡平均半径
R =0.5 R 1+R 3 =14mm
假设内壁温度:T wg =500K 。
利用巴兹法计算燃气与内壁面的对流换热系数: 根据T
T wg
st
=0.166,查表得到考虑附面层内燃气性能变化的修正系数σ=1.52。
燃气与内壁面的对流换热系数
0.026μ0.2cppc0.8Dt0.1At0.9
g= ς=1067W/(m2⋅K)
Dt
ns燃气与内壁面的对流换热密度qk= g Tst−Twg =2.675×106W∕m2
5.1.3计算燃气与内壁面的辐射热流密度
c1
根据D=1.78,查得气体对整个壁面辐射的平均射线长
c
L
l=0.85Dc=54.4mm
水蒸气分压pH2O=pc⋅
nHOnn
=1.12MPa =0.34MPa
二氧化碳分压pCO2=pc⋅
nCO计算得到pH2Ol=0.061MPa⋅mpCO2l=0.018MPa⋅m
查图得水蒸气发射率ε0H 2O =0.014,指数关系n=1+kH2OpH2O=1.6,则水蒸气的实际发射率为
εH2O=1− 1−ε0H2O =1− 1−0.014 1.6=0.0223
查图得二氧化碳发射率εC O2=0,则总的发射率为
εg=εH2O+εCO2−εH2OεCO2=0.0223
壁面发射率一般取为ε0w =0.8,则实际有效壁面发射率为
εW=ε0W 1+ 1−ε0W 1−εg =0.956
由于壁面温度较低,故壁面对燃气的辐射可以忽略,因此燃气辐射热流密度
n
为
Tg4Twg4Tg4
qr=5.67⋅εw⋅ εg −αg =5.67⋅εw⋅εg =9.870×104W∕m2
5.1.4计算总热流密度、总热流量及冷却剂流量
总热流密度q=qk+qr=2.774×106W∕m2
取推力室圆筒段壁厚δ=1mm,则外壁面温度 TWf=TWg−总热流量Φ=q⋅2πRclc=6.36×104W
若要求冷却水通过冷却通道时的温升为40K ,则冷却水流量为
Φ
q co ==0.381kg/s
pco co 冷却水的温度可以取为
T co =
20+602
qδλw
=492K
+273 K =313K
5.1.5确定冷却通道参数
推力室内壁面及肋条材料为1Cr18Ni9Ti ,导热系数λW =23.26W ∕ m ⋅K ,内壁厚δ=1mm ,冷却通道高度h=3mm,肋条厚b=1.5mm,冷却通道宽a=2mm,则当量直径ⅆe=a+b=2.4mm,冷却通道数n=57。
2a
5.1.6计算内壁面和外壁面温度
冷却通道面积 冷却剂流速
e=
A=na =3.42×10−4m2 Vt=
μco
qcoA⋅ρco
=1.114m/s
1000×1.25×2.4×10−3
0.658×10冷却剂雷诺数 冷却剂普朗特数 冷却剂努塞尔数
ρcoVtde
==4063
Pr=
Cpcoμcoλco
=4.31
0.4
Nu f =0.023R 0.8=31.8 e Pr
冷却剂和外壁面的对流换热系数
f=
Nufλcode
=8321W∕m2
液体壁面温度 气体壁面温度
TWf=η
q
p f
+Tf=480K
Twg=λ+Twf=488K
w
qδ
由计算结果可知,推力室圆筒段气体壁面Twg=484K小于假定的温度500K, 相差2.4%小于5%,符合冷却要求。
5.2推力室喉部冷却计算
5.2.1燃气的气动参数
喉部燃气温度 燃气多变指数 燃气定压比热容 燃气粘度
T st =2845K 。 k=1.1478。 C p =2.28kJ ∕kg ⋅K μ=0.0001Pa ⋅s Pr=0.5967 pt=1.1552MPa
燃气普朗特数 喉部压力
5.2.2计算燃气与内壁面的对流换热密度
喉部横截面积
A=4πDt2=2.01×10−4m2
1
喷管喉部过渡半径 R =0.5 R 1+R 3 =14mm 假设内壁温度:
T wg =740K 。
利用巴兹法计算燃气与内壁面的对流换热系数 根据T
T wg
st
=0.26,查表得到考虑附面层内燃气性能变化的修正系数σ=1.35。
燃气与内壁面的对流换热系数
0.026μ0.2Cppc0.8Dt0.1At0.9
g= ς=11320W/(m2⋅K)
t
ns燃气与内壁面的对流换热密度
qk= g Tst−Twg =2.38×107W∕m2
5.2.3计算燃气与内壁面的辐射热流密度
c1
根据D=1.78,查得气体对整个壁面辐射的平均射线长
c
L
l=0.85Dc=54.4mm
水蒸气分压
pH2O=pt⋅pCO2=pt⋅
nHOnnCOn
=0.670MPa =0.208MPa
pCO2l=0.011MPa⋅m
二氧化碳分压 计算得到
pH2Ol=0.036MPa⋅m
查图得水蒸气发射率ε0H 2O =0.01,指数关系n=1+kH2OpH2O=1.4,则水蒸气的实际发射率为
εH2O=1− 1−ε0H2O =0.014
查图得二氧化碳发射率εC O2=0,则总的发射率为
εg=εH2O+εCO2−εH2OεCO2=0.014
壁面发射率一般取为ε0w =0.8,则实际有效壁面发射率为
εW=ε0W 1+ 1−ε0W 1−εg =0.958
由于壁面温度较低,故壁面对燃气的辐射可以忽略,因此燃气辐射热流密度为
Tg4Twg4Tg4qr=5.67⋅εw⋅ εg −αg =5.67⋅εw⋅εg =49820W∕m2
n
5.2.4计算总热流密度、总热流量及冷却剂流量
总热流密度
q=qk+qr=2.39×107W∕m2
取推力室圆喉部壁厚δ=1mm,则外壁面温度
TWf=TWg
qδ2.39×107×1×10−3−= 740−K=672K w取喉部截面长度10mm 总热流量
Φ=q⋅S=1.15×104W
Φ1.15×104ΔT ===6.5K
co co 取定性温度313K
冷却水流量为q co =0.421kg/s ,则喉部内冷却水温升:
5.2.5确定冷却通道参数
推力室内壁面及肋条材料为1Cr18Ni9Ti ,导热系数λW =23.26W ∕ m ⋅K ,内壁厚δ=1mm ,冷却通道高度h=3mm,肋条厚b=1.5mm,冷却通道宽a=2mm,则当量直径ⅆe=a+b=2.4mm,冷却通道数n=14。
2a
5.2.6计算内壁面和外壁面温度
冷却通道面积 冷却剂流速
A=14a =8.4×10−5m2
co
Vt=A⋅ρ=5.01m/s
co
q
冷却剂雷诺数 冷却剂普朗特数 冷却剂努塞尔数
Re=Pr=
ρcoVtdeμcoCpcoμcoλco
=18357 =4.31
0.4Nu f =0.023R 0.8=106 e Pr
冷却剂和外壁面的对流换热系数
Nufλco
f==28223W∕m2
e
液体壁面温度 气体壁面温度
TWf=η
q
p f
+Tf=683K
Twg=λ+Twf=727K
w
qδ
由计算结果可知,推力室喉部气体壁面Twg=727K大于假定的温度740K, 相差1.7%小于5%,,且当Twg=727K时1Cr18Ni9Ti 不会失效,故符合冷却要求。
6. 推力室强度校核
6.1推力室圆筒段强度校核
圆筒段所选材料为1Cr18Ni9Ti ,查机械手册得:抗拉强度为ςb=550MPa,屈服极限为ςs=200MPa
将推力室圆筒段视为承压薄壁圆筒,其周向应力为
pcDc
ς1=ςt==32MPa
1轴向应力为
ς2=ςx=
径向应力为
ς3=ςr=pc=2 MPa
推力室圆筒段内壁温度约400K,温度较高,材料的屈服极限有所下降,故采用较大的安全系数,取n=2.5
许用应力为
ςs =
由第三强度理论有
ςr3=ς1−ς3=30MPa
由第四强度理论有
ςr4
(ς1−ς2)+(ς2−ς3)+(ς3−ς1)= =30.3MPa
2
2
2
F
=1.22MPa ςs
=80MPa 由此推力室圆柱段满足强度条件
6.2喷管强度校核
喷管所选材料为1Cr18Ni9Ti ,查机械手册得:抗拉强度为ςb=550MPa,屈服极限为ςs=200MPa
喷管喉部周向应力
ptDt
ς1=ςt==4.58MPa
1
喷管喉部轴向应力为
ς2=ςx=
喷管喉部径向应力为
ς3=ςr=pt=1.207MPa
F514==4.72MPa 推力室圆筒段内壁温度约640K,温度更高,材料的屈服极限下降更多,故采用更大的安全系数,取n=5
许用应力为
ςs =
由第三强度理论有
ςr3=ς1−ς3=3.483MPa
由第四强度理论有
(ς1−ς2)+(ς2−ς3)+(ς3−ς1)
ςr4= =3.41MPa
2
2
2
ςs
=40MPa 由此推力室喷管段满足强度条件
综上,推力室各段在保证发动机性能和冷却要求的前提下均满足强度条件。
7. 点火器设计
点火器采用棒料加工而成,在棒料的一侧边钻孔,用来安装酒精管嘴和电点火塞,气氧管嘴安装在上部平面。为保证点火腔有一定的压力,应设计燃气通道喉部,具体计算如下:
At(点火)=
C∗∙q点火点火
=2.831×10−5m2
dt(点火)=6mm
8. 螺栓强度校核
燃烧室压力为2MPa,酒精腔压力为2.4MPa,气氧腔压力为2.5MPa,喷管出口压强为1atm 。本地面试验用气氧酒精发动机采用强度等级为4.6级受拉螺栓,螺栓材料为钢,其强度极限为
ςb=400MPa
则其屈服极限为
ςs=0.6×400=240MPa
取安全系数为1.5,则许用屈服极限为
σ =
240
=160MPa 综合螺栓受力,得到与下底和法兰连接的螺栓所受的力为
p∙Sσ==102.9MPa
上底与中底采用螺钉连接,校核如下
σ=
p∙S
=14.8MPa
9. 整体结构分析
为便于设计加工,结构上多采用焊接进行连接。该发动机是单件生产,因此大部分零件均采用机械加工,零件尽量采用标准件。发动机燃烧室燃气温度较高,对材料耐高温性能均有很高的要求。
9.1头部结构
该发动机采用三底两腔的经典结构。上底和中底通过螺钉连接。该发动机仅用于地面,故可以采用简单且可重复使用的电点火方式。
气氧和酒精分别通过不同的管路进入两个腔体,并从三个直流撞击式喷嘴喷入燃烧室燃烧。
9.2喷注器
本发动机是气氧酒精双组元发动机,采用直流撞击式双组元喷注器。作为燃料的酒精通过自击式喷嘴喷注进入燃烧室,这样液滴细度高好而且混合距离短。作为氧化剂的气氧,从结构简单原则敖律,让其通过中间的直流喷嘴喷出,与酒精混合即可。
喷注盘面较小,故喷嘴数量不宜过多,可选择三个喷注器并均匀分布在点火器出口的圆周上。
9.3点火器
该发动机用于地面试验,故可以使用电点火,结构简单,可靠性高,可重复使用,点火组元仍然采用气氧和酒精,使其燃烧产物同推力室一致,预燃室通过螺纹与上底连接。点火器的燃气通过其喉部进入燃烧室点燃主推进剂。
9.4推力室
燃烧室为圆筒形,圆筒段采用的材料为1Cr18Ni9Ti 。
圆筒段采用的冷却方式为外冷却,冷却通道为铣槽结构,有56个冷却通道。 该发动机为试验用,其收缩段采用双圆弧型面,扩张段也采用双圆弧型面。与燃烧室为一体,冷却方式为外冷却,喉部有18个冷却通道,。
喷管冷却套采用两个对称的版套筒,用焊接的方式连接而成。
9.6密封结构
该发动机中需要考虑的密封较多,为简单起见,在没有高温燃气处的密封均采用国家标准液压气动用o 型橡胶密封圈,在接触高温燃气处采用矩形紫铜密封圈。
10. 感悟
经过此次课程设计,我对于火箭发动机设计方面的理论知识有了更深的理解,同时也对于机械设计的流程有了进一步的掌握。在本次设计中,张老师曾多次向我提出某零件如何安装,如何定位,如何确保密封的提问,让我意识到自己在多方面考虑火箭发动机设计问题方面还有着很多的不足。而此次的课程设计也正是为我提供了这样一个实践的机会,来了解设计人员在设计火箭发动机时所要考虑的多方面问题。不仅如此,此次火箭发动机课程设计的重头戏——手工绘制设计图也同样让我感受到自己在工程制图方面的不足,同时也为我提供了锻炼的机会。总之,在此次课程设计中,我的收获很大,从理论计算到实际制图的能力都有了不小的提升。感谢学院安排的此次课程,感谢张黎辉老师的悉心教导。
参考文献
[1] 邢继发,刘国球,等. 世界导弹与航天发动机大全. 北京:军事科学出版社,1999
[2] 蔡国飙,李家文,等. 液体火箭发动机设计. 北京:北京航空航天大学出版社,2011
[3] 王之栎,王大康,等. 机械设计综合课程设计. 北京:北京航空航天大学出版社,2008
[4]D.K.休泽尔等. 液体火箭发动机现代工程设计. 北京:宇航出版社,2004 [5] 杨世铭,陶文铨,等. 传热学. 北京:高等教育出版社,1998