《火箭发动机》 7 内弹道

合集下载

内弹道基础概述内弹道学

内弹道基础概述内弹道学
内弹道基础概述内弹道学
➢火药的燃气生成函数
1)简单形状火药
Z/1 相对已燃厚度
S/S1 相对已燃表面
V/V1 相对已燃体积
内弹道基础概述内弹道学
Z1ZZ2
12Z3Z2
式中:
内弹道基础概述内弹道学
内弹道基础概述内弹道学
2)多孔火药
基于单孔的管状药接近定面燃烧的概念,为了使火 药具有增面燃烧性能,于是又产生了增加内孔的多孔 火药系列。但是多孔药与管状药不同,多孔火药在燃 完厚度的瞬间,火药却未全部燃尽,而是分裂成若干 碎粒。因此,多孔火药的燃烧存在两个阶段,即分裂 前的主体燃烧阶段和分裂后的碎粒燃烧阶段。多孔药 燃烧的增面性,只存在于主体燃烧阶段。在碎粒燃烧 阶段,则是强烈的减面性。
6 内弹道基础
内弹道基础概述内弹道学
6.1 概述
内弹道学是专门研究弹丸在膛内运动规律的科学。 研究的对象是膛内的射击现象,包括火药在膛内的 燃烧规律、弹丸运动的规律,以及膛内压力变化规律 等方面的内容。
内弹道基础概述内弹道学
6.2 膛内射击过程
➢枪炮发射系统
典型系统的包括身管、火药和弹丸。
1—炮闩;2--药室;3—坡膛;4—线膛
s sZ1sZ 1ZZb
内弹道基础概述内弹道学
➢弹丸运动方程
系统动量守恒:
m (v v 1 )v M 1 0 v
燃气和未燃药粒速度假设:
vv12 (vv1)2 vv1
可得: vv M2 v 1 Mm
md(vv1)S dt
pRx
M 2
1
M
m
(1
Rx
)
M2
mdvSP(1Rx )
内弹道基础概述内弹道学
内弹道基础概述内弹道学

固体火箭发动机内弹道性能的仿真研究_刘宝华

固体火箭发动机内弹道性能的仿真研究_刘宝华

2009年10月第10期电 子 测 试ELEC TRONIC T ES TOct.2009No.10固体火箭发动机内弹道性能的仿真研究刘宝华1,杨志菊2(1辽宁葫芦岛市海军飞行学院教研部 葫芦岛 125000,2辽宁葫芦岛92941部队 葫芦岛 125001) 摘 要:研究发动机燃烧室内压强随时间变化的规律,是固体火箭发动机工作过程分析的主要方面,一般采用实验法可以获得直观可靠的数据,但实验法耗资大、周期长,不易操作。

本文根据零维内弹道数学模型,运用龙格库塔法及M A T L A B语言对某型固体导弹发动机内弹道工作过程进行数值仿真,画出燃烧室内压强随时间变化的曲线,并进一步分析得出影响发动机内弹道性能的因素,仿真结果与发动机燃烧室内工作情况相符,影响因素分析与实验相一致,为固体火箭发动机内弹道性能研究提供了参考。

关键词:发动机;内弹道;数值仿真;燃烧室中图分类号:T N957 文献标识码:BSimulate research of solide rocket eng ine inside trajecto ry characteristicLiu Baohua1,Y ang Zhiju2(1T eaching and scientific resea rch ministry,N aval Flying A czdemy,Huludao125001,China;292941U nit,PL A,H uludao125001,China)A bstract:Researching the rule of the pressure with the time changing o f the eng ine's firebox is the m ost aspect of analy zing the process o f the solid missile engine.Co mmo nly,experimenta-tion can achieve the intuitionistic and credible data,but this method consume m uch tim e and money,and not easy to manipulate.According to zero-dimensional inside trajecto ry numerical simulate model,to o ur co untry some ty pe solid missile engine,applacating Lo nger-Kuta me thod and the lang uage of m atlab launches the inside trajecto ry numerical simulate of the w o rking pro cess,and draw ed the burner pressure-time curve,passing further analy sis o uted the influence facto rs of inside trajecto ry characteristic.The result of simulatio n is consistent to the fact in en-gine's firebox and analysis of the influence facto rs is consistent to experim entation.It offer s a new idea to researching the capability o f the ballistic trajecto ry.Keywords:engine;inside trajectory;numerical sim ulate;firebox0 引 言固体火箭发动机内弹道学的核心是研究发动机燃烧室内压强随时间变化的规律,是固体火箭发动机工作过程分析的主要方面。

固体火箭发动机内弹道学 方丁酉

固体火箭发动机内弹道学 方丁酉

固体火箭发动机内弹道学方丁酉固体火箭发动机内弹道学方丁酉近年来,随着航天事业的蓬勃发展,固体火箭发动机内弹道学成为越来越受关注的领域。

方丁酉,中国工程院院士、火箭技术专家,曾在这一领域做出了卓越的贡献。

本文将从简单到复杂、由浅入深地探讨固体火箭发动机内弹道学,以帮助读者全面、深入地理解这一话题。

一、固体火箭发动机内弹道学的概念固体火箭发动机内弹道学是研究固体火箭发动机内燃烧过程、燃料燃烧特性、燃烧产物排放、工作环境等方面的学科。

方丁酉指出,固体火箭发动机内弹道学是固体火箭技术中的重要基础理论,对于提高固体火箭发动机的性能、可靠性和运载能力具有重要意义。

二、固体火箭发动机内弹道学的理论基础在固体火箭发动机内弹道学的研究中,燃烧动力学是重要的理论基础之一。

燃烧动力学研究了燃料在燃烧过程中的变化规律,以及燃烧反应对火箭发动机内部气体流动和压力变化的影响。

燃烧产物排放和燃烧室内部气体动力学也是固体火箭发动机内弹道学的重要内容之一。

研究这些理论基础可以帮助我们更好地理解固体火箭发动机内部的工作原理和特性。

三、固体火箭发动机内弹道学的关键技术在固体火箭发动机内弹道学研究中,燃烧稳定性和效率是两个关键技术。

燃烧稳定性是指在固体火箭发动机工作过程中保持燃烧的稳定性,避免出现燃烧不均匀或燃烧失稳等问题。

而燃烧效率则是指在燃烧过程中尽可能提高燃料的利用率,减少燃料的浪费。

方丁酉在固体火箭发动机内弹道学的研究中,提出了一系列有效的技术方案,使固体火箭发动机在燃烧稳定性和效率方面取得了显著的进展。

四、固体火箭发动机内弹道学的应用前景固体火箭发动机内弹道学的研究成果已经在我国的航天事业中得到了广泛的应用。

在长征系列火箭、嫦娥探月工程等多个航天工程中,固体火箭发动机内弹道学的研究成果为提高火箭的性能和可靠性做出了重要贡献。

未来,随着我国航天事业的不断发展,固体火箭发动机内弹道学将继续发挥重要作用,推动我国航天事业迈向新的高度。

固体火箭发动机燃烧室与内弹道计算

固体火箭发动机燃烧室与内弹道计算

t
x
x
.
(V

V x
dx)dt

pr
Ab x
dxdt
m
整理后有:
t
(
A)

x
(
AV
)


pr
Ab x
燃烧面积
.. m+dm
. dm
dm


p
r
Ab x
dx
(2)动量方程
根据动量守恒:在dt时间内,微元体中动量的变化量=燃气 流动迁移在微元体内产生的动量增量+作用于两截面上的总压力 冲量+作用于微元体侧表面上压力冲量在x轴上的分量,即

x
(Ec

V2 2
)dx

dt

x
(
pAV
)dxdt


pr
Ab x
dxdt

Ip
整理后有:
t

A(Ec

V2 2
)

x

AV
(Ic

V2 2
)

prI
p
Ab x
Ec : 单位质量燃烧产物的内能(化学内能+热内能) Ic : 单位质量燃烧产物的总焓
1. 流动特点
燃烧室中的燃气流动具有以下特点: 燃烧室中推进剂装药不断燃烧,产生新的燃气加入主流,因此燃 气在燃烧室中的流动是一个有质量加入的流动过程; 在推进剂装药燃烧过程中,燃气的通道面积不断增大,故燃气的 流动参数是位置与时间的函数,因而是非定常流动; 燃气流过形状复杂的通道或流线变化剧烈时,将产生一定的能量 损失; 燃气的流动特性与推进剂的燃烧特性密切相关。例如燃气的压强 和流速影响推进剂的燃速,而燃速又反过来影响燃气压强和流速; 由此可见,在发动机工作过程中,燃气在燃烧室中的流动中极其 复杂的。为了研究方便,特作如下假设: 推进剂燃烧在装药燃烧表面上瞬时完成,燃气的化学组分和热力 性质均匀一致,并沿燃烧表面的外法线方向注入通道。 燃气为完全气体,服从完全气体状态方程,而且比热不变。 燃气在装药通道中作一维流动。 燃气与外界无热和功的交换。

瞬时平衡法计算固体火箭发动机内弹道曲线

瞬时平衡法计算固体火箭发动机内弹道曲线

!计?算?星?型æ装痢?药?的?几?何?尺?寸?!参?数簓符?号?说μ明é!-----------------------------------------------------------------------------------!n表括?示?星?角?数簓,num表括?示?将?推?进?剂á沿?肉╝厚?方?向ç分?为a几?等台?分?,m表括?示?选?择?压1力 值μ大洙?小?!d表括?示?外猘径?,len表括?示?长¤度â,thet表括?示?星?边?夹D角?,epsilon表括?示?角?度â系μ数簓!r表括?示?过y度â圆2弧?半?径?,r1表括?示?星?角?圆2弧?半?径?,R0表括?示?通 ?用?气?体?常£数簓!l表括?示?药?柱ê的?特?征é长¤度â,y0表括?示?初?始?特?征é参?数簓,y1表括?示?燃?尽?特?征é参?数簓!I0表括?示?总哩?冲?,F表括?示?推?力 ,Poc表括?示?燃?烧?室酣?的?工¤作痢?压1力!Isp表括?示?比括?冲?,density_p表括?示?密ì度â,k表括?示?比括?热â?比括?rspeed表括?示?燃?速ê,?pn表括?示?压1力 指?数簓!mpeff表括?示?有瓺效§装痢?药?量?,Cf表括?示?推?力 系μ数簓,Ctz表括?示?特?征é速ê度â,At表括?示?喉æ部?面?积y,a表括?示?燃?速ê系μ数簓!S表括?示?平?均ê燃?烧?面?积y,e1表括?示?平?均ê肉╝厚?,epsilon1表括?示?减?面?比括?epsilon2表括?示?增?面?比括?!foresmax表括?示?前©段?最?大洙?相à对?周ì边?长¤,backsmax表括?示?后è段?最?大洙?相à对?周ì边?长¤!smin表括?示?最?小?相à对?周ì边?长¤,thet1表括?示?周ì边?长¤取?得?最?小?值μ时骸?的?星?边?夹D角?!Ap表括?示?初?始?通 ?气?面?积y,J表括?示?初?始?通 ?气?参?量?,eta表括?示?装痢?填?系μ数簓,Af表括?示?剩骸?药?面?积y,etaf表括?示?剩骸?药?系μ数簓!de表括?示?每?一?份 肉╝厚?的?长¤度â,Sa表括?示?燃?烧?面?积y数簓组哩?Apa表括?示?通 ?气?面?积y数簓组哩?!------------------------------------------------------------------------------------program mainimplicit nonereal(kind=8),parameter :: Pi=3.14integer :: n,num,i,mreal(kind=8) :: d,len,thet,epsilon,r,r1real(kind=8) :: l,y0,y1real(kind=8) :: I0,F,Poc,Pe,R0real(kind=8) :: Isp,density_p,k,rspeed,Pnreal(kind=8) :: mpeff,Cf,Ctz,At,areal(kind=8) :: S,e1,epsilon1,epsilon2real(kind=8) :: foresmax,backsmax,smin,thet1real(kind=8) :: Ap,J,eta,Af,etafreal(kind=8) :: error,dereal(kind=8),allocatable :: Sa(:),Apa(:)!读á入?所ê需â要癮所ê用?参?数簓值μopen(3,file="design_parameter.dat")read(3,*) I0 !总哩?冲?read(3,*) F !推?力read(3,*) Poc !燃?烧?室酣?压1力read(3,*) Isp !比括?冲?read(3,*) density_p !推?进?剂á的?密ì度âread(3,*) k !比括?热â?比括?read(3,*) a !燃?速ê系μ数簓read(3,*) Ctz !读á取?特?征é速ê度âCtzread(3,*) Pn !压1力 指?数簓read(3,*) Pe !读á取?喷?管ì出?口ë的?压1力 Peread(3,*) R0 !读á取?通 ?用?气?体?常£数簓R0read(3,*) d,r !读á取?装痢?药?直ª径?d和æ过y度â圆2弧?半?径?rread(3,*) epsilon2 !读á取?增?面?比括╡psilon2read(3,*) r1 !读á取?星?角?圆2弧?半?径?r1read(3,*) thet1 !试?取?周ì边?长¤取?得?最?小?值μ时骸?的?星?边?夹D角?thet1read(3,*) thet !试?取?初?始?时骸?的?星?边?夹D角?thetread(3,*) n !读á取?星?角?数簓n(辍?,4,5,6,7,8)?read(3,*) epsilon !试?取?角?度â系μ数簓epsilonread(3,*) num !读á取?等台?分?肉╝厚?的?等台?分?数簓numread(3,*) m !m=1表括?示?最?小?压1力 给?定¨,m=2表括?示?最?大洙?压1力 给?定¨,m=0表括?示?平?均ê压1力 给?定¨close(3)allocate(Sa(0:num),Apa(0:num))!根ê据Y规?定¨的?总哩?冲?计?算?有瓺效§装痢?药?量?mpeffmpeff=1.02*I0/Isp!计?算?推?力 系μ数簓和æ喉æ部?面?积yCf=sqrt(k)*(2/(k+1))**((k+1)/(k-1)/2)*&sqrt(2*k*(1-(Pe/Poc)**((k-1)/k))/(k-1))At=F/Cf/(Poc*101325.0)!计?算?平?均ê肉╝厚?e1和æ计?算?特?征é长¤度â和æ平?均ê燃?烧?面?积yS!若?给?定¨的?是?最?大洙?压1力 则ç需â根ê据Y增?面?比括?计?算?最?小?燃?面?再ê计?算?平?均ê燃?面?值μ!y1一?般?取?值μ在ë1附?近ì,(0.8-1.2)if(m==0) then!平?均ê压1力 给?定¨S=At*(Poc*101325.0)**(1-Pn)/(Ctz*density_p*a)e1=mpeff/(density_p*S)else if(m==1) then!最?小?压1力 给?定¨S=At*(Poc*101325.0)**(1-Pn)/(Ctz*density_p*a)S=(epsilon2+1.0)*S/2.0e1=mpeff/(density_p*S)else!最?大洙?压1力 给?定¨S=At*(Poc*101325.0)**(1-Pn)/(Ctz*density_p*a)S=(1.0/epsilon2+1.0)*S/2.0e1=mpeff/(density_p*S)end ifl=d/2-e1-ry0=r/ly1=(e1+r)/l!根ê据Yepsilon2选?择?角?度â系μ数簓epsilon!计?算?取?得?最?小?周ì长¤时骸?的?thet1do while(.true.)error=thet1/2+cotan(thet1/2)-Pi/n-Pi/2if(error>=0.0001) thenthet1=thet1+0.00001else if(error<=-0.0001) thenthet1=thet1-0.00001elseexitend ifend do!计?算?角?度â系μ数簓epsilon的?准?确í?值μdo while(.true.)backsmax=2*n*((1-epsilon)*Pi/n+y1*(Pi/n+asin(sin(epsilon*Pi/n)/y1))) smin=2*n*(sin(epsilon*Pi/n)/sin(thet1/2)+(1-epsilon)*Pi/n) if(backSmax/Smin>epsilon2) thenepsilon=epsilon+0.001elseexitend ifend do!根ê据Y减?面?比括╡psilon1,?计?算?thetepsilon1=1/epsilon2foresmax=smin/epsilon1do while(.true.)error=2*n*(sin(epsilon*Pi/n)/sin(thet/2)+(1-epsilon)*Pi/n+&(r1+r)*(Pi/n+Pi/2-thet/2-cos(thet/2)/sin(thet/2))/l)-foresmax if(error>=0.0001) thenthet=thet+0.00001else if(error<=-0.0001) thenthet=thet-0.00001elseexitend ifend do!计?算?初?始?通 ?气?面?积y,?通 ?气?参?量?,?装痢?填?系μ数簓,药?柱ê长¤度âAp=(n*((1-epsilon)*Pi/n+sin(epsilon*Pi/n)*(cos(epsilon*Pi/n)-&sin(epsilon*Pi/n)*cotan(thet/2)))+2*n*r*(sin(epsilon*Pi/n)/sin(thet/2)+&(1-epsilon)*Pi/n)/l+n*r**2*(Pi/n+Pi/2-thet/2-cotan(thet/2))/l**2+&n*r1**2*(thet/2+cotan(thet/2)-Pi/2)/l**2)*l**2J=At/Apeta=4*(Pi*d**2/4-Ap)/(Pi*d**2)Af=(epsilon*Pi*(1+y1)**2-n*(sin(epsilon*Pi/n)*(sqrt(y1**2-sin(epsilon*Pi/n)**2)+& cos(epsilon*Pi/n)))-n*y1**2*(epsilon*Pi/n+asin(sin(epsilon*Pi/n)/y1)))*l**2etaf=4*Af/(pi*d**2)len=mpeff/density_p/(Pi*d**2/4-Ap-Af)!计?算?推?进?剂á的?燃?面?变?化ˉ规?律ã并¢输?出?结á果?de=e1/numopen(10,file="export_burnS.dat")open(20,file="export_Ap.dat")do i=0,num,1if((i*de)<=r1) thenSa(i)=2*n*(sin(epsilon*Pi/n)/sin(thet/2)+(1-epsilon)*Pi/n+&(r1+r)*(Pi/n+Pi/2-thet/2-cotan(thet/2))/l-(r1-i*de)*Pi/n/l)*l*len Apa(i)=(n*((1-epsilon)*Pi/n+sin(epsilon*Pi/n)*(cos(epsilon*Pi/n)-sin(epsilon*Pi/n)* cotan(thet/2)))+&2*n*(r+i*de)*(sin(epsilon*Pi/n)/sin(thet/2)+(1-epsilon)*Pi/n)/l+&n*(r+i*de)**2*(Pi/n+Pi/2-thet/2-cotan(thet/2))/l**2+n*(r1-i*de)**2*(thet/2+cotan(thet/2 )-Pi/2)/l**2)*l**2else if((i*de)>r1.and.(i*de)<=(l*sin(epsilon*Pi/n)/cos(thet/2)-r)) thenSa(i)=2*n*(sin(epsilon*Pi/n)/sin(thet/2)+(1-epsilon)*Pi/n+&(i*de+r)*(Pi/n+Pi/2-thet/2-cotan(thet/2))/l)*l*lenApa(i)=(n*((1-epsilon)*Pi/n+sin(epsilon*Pi/n)*(cos(epsilon*Pi/n)-sin(epsilon*Pi/n)* cotan(thet/2)))+&2*n*(r+i*de)*(sin(epsilon*Pi/n)/sin(thet/2)+(1-epsilon)*Pi/n)/l+&n*(r+i*de)**2*(Pi/n+Pi/2-thet/2-cotan(thet/2))/l**2)*l**2else if((i*de)>(l*sin(epsilon*Pi/n)/cos(thet/2)-r).and.(i*de)<=e1) thenSa(i)=2*n*((1-epsilon)*Pi/n+(r+i*de)*(Pi/n+asin(l*sin(epsilon*Pi/n)/(i*de+r)))/l)*l *lenApa(i)=n*((1-epsilon)*Pi*(1+(r+i*de)/l)**2/n+sin(epsilon*Pi/n)*(sqrt((r+i*de)**2/l* *2-&sin(epsilon*Pi/n)**2)+cos(epsilon*Pi/n))+(r+i*de)**2*(epsilon*Pi/n+&asin(l*sin(epsilon*Pi/n)/(r+i*de)))/l**2)*l**2end ifwrite(10,"(f10.5,2X,f15.5)") i*de,Sa(i)write(20,"(f10.5,2X,f15.5)") i*de,Apa(i)end doclose(10)close(20)!调獭?用?子哩?程å序ç计?算?装痢?药?内ë弹獭?道台?曲ë线?call internal_ballistics0(num,de,Sa,Apa,len,a,pn,Ctz,At,k,density_p,Isp,Pe)!输?出?星?型æ装痢?药?的?几?何?参?数簓open(30,file="export_star_geometry.dat")write(30,"(A6,f11.5,A2)") "mpeff=",mpeff,"Kg"write(30,"(A2,f8.5,A2)") "D=",d,"m"write(30,"(A2,f8.5,A2)") "L=",len,"m"write(30,"(A2,I2)") "n=",nwrite(30,"(A3,f8.5)") "θâ=",thetwrite(30,"(A3,f7.5)") "ε?=",epsilonwrite(30,"(A2,f7.5,A2)") "r=",r,"m"write(30,"(A3,f7.5,A2)") "r1=",r1,"m"write(30,"(A3,f8.5,A2)") "e1=",e1,"m"write(30,"(A3,f10.5,A2)") "Ap=",Ap,"㎡O"write(30,"(A2,f10.5)") "J=",Jwrite(30,"(A3,f7.5)") "η?=",etawrite(30,"(A3,f10.5,A2)") "Af=",Af,"㎡O"write(30,"(A4,f10.5)") "η?f=",etafwrite(30,"(A2,f10.5)") "l=",lwrite(30,"(A3,f10.5)") "y0=",y0write(30,"(A3,f10.5)") "y1=",y1close(30)stopend!该?子哩?程å序ç用?于 ?计?算?零?维?变?截?面?燃?烧?装痢?药?的?内ë弹獭?道台?!该?子哩?程å序ç利?用?瞬2时骸?平?衡a法ぁ?计?算?内ë弹獭?道台?!同?时骸?可ã以?计?算?侵?蚀骸?燃?烧?的?影 ?响å!--------------------------------------------------------------------------!Pc_eq表括?示?工¤作痢?段?的?平?衡a压1力 值μ,dPc_eq表括?示?压1力 修T正y值μ,rspeed表括?示?燃?速ê!time表括?示?工¤作痢?段?的?时骸?间?,Toc表括?示?燃?气?温?度â,Roc表括?示?燃?气?气?体?常£数簓,Pig表括?示?点?火e压1力!P1_eq表括?示?起e始?段?的?压1力 值μ,P2_eq表括?示?拖?尾2段?的?压1力 值μ,Cf表括?示?推?力 系μ数簓!time1表括?示?起e始?段?的?时骸?间?,time2表括?示?拖?尾2段?的?时骸?间?,Isp表括?示?理え?论?比括?冲?值μ!Feq表括?示?工¤作痢?段?的?推?力 值μ,F1表括?示?起e始?段?的?推?力 值μ,F2表括?示?拖?尾2段?的?推?力 值μ!Poc表括?示?燃?烧?室酣?压1力 值μ,Pe表括?示?喷?管ì出?口ë压1力!--------------------------------------------------------------------------subroutine internal_ballistics0(num,de,Sa,Apa,len,a,pn,Ctz,At,k,density_p,Isp,Pe) implicit noneinteger,parameter :: n=10integer :: num,i,keyreal(kind=8),intent(in) :: de,len,a,pn,Ctz,At,k,density_preal(kind=8),intent(in) :: Sa(0:num),Apa(0:num),Isp,Pereal(kind=8) :: Pc_eq(0:num),dPc_eq(0:num),rspeed(0:num)real(kind=8) :: time(0:num),Toc,Roc,Feq(0:num),Cfreal(kind=8),allocatable :: P1_eq(:),P2_eq(:),F1(:)real(kind=8),allocatable :: time1(:),time2(:),F2(:)real(kind=8) :: ga,Pigreal(kind=8),external :: funopen(3,file="inter_parameter.dat")read(3,*) key !读á入?key的?值μread(3,*) Pig !读á入?点?火e压1强?值μread(3,*) Toc !读á入?绝?热â?燃?烧?温?度â值μread(3,*) Roc !读á入?燃?气?的?气?体?常£数簓值μallocate(P1_eq(n),F1(n),time1(n))allocate(P2_eq(0:n-1),F2(0:n-1),time2(0:n-1))ga=sqrt(k)*(2/(k+1))**((k+1)/(k-1)/2)Pig=Pig*101325.0!-------------------------------------------------!计?算?工¤作痢?段?的?压1力 和æ时骸?间?值μ!-------------------------------------------------!计?算?各é燃?层?时骸?燃?烧?室酣?的?平?衡a压1力 Pc_eqdo i=0,num,1Pc_eq(i)=(density_p*Ctz*fun(Sa(i),Apa(i),key)*&a*Sa(i)/At)**(1/(1-pn))end do!计?算?各é燃?层?时骸?平?衡a压1力 的?修T正y值μdPc_eqdPc_eq=0.0do i=1,num,1dPc_eq(i)=Pc_eq(i)*(Sa(i)-Sa(i-1))/(Sa(i)*(1-pn))end do!计?算?各é燃?层?的?实害?际ä平?衡a压1力 值μdo i=0,num,1Pc_eq(i)=Pc_eq(i)-Apa(i)*len*fun(Sa(i),Apa(i),key)*&a*dPc_eq(i)/(de*Ctz*ga*ga*At)end do!计?算?各é燃?层?时骸?的?推?进?剂á燃?速ê值μrspeeddo i=0,num,1rspeed(i)=fun(Sa(i),Apa(i),key)*a*Pc_eq(i)**pnend do!计?算?各é燃?层?时骸?的?推?进?剂á的?推?力 值μFeqdo i=0,num,1Cf=0.0Cf=ga*sqrt(2*k*(1-(101325.0*Pe/Pc_eq(i))**((k-1)/k))/(k-1))Feq(i)=Cf*Pc_eq(i)*Atend do!计?算?燃?到?各é层?需â要癮的?时骸?间?值μtimetime=0.0do i=1,num,1time(i)=time(i-1)+de*2.0/(rspeed(i)+rspeed(i-1))end do!-----------------------------------------------------------!计?算?起e始?段?的?压1力 值μ、¢推?力 和æ时骸?间?!-----------------------------------------------------------do i=1,n,1P1_eq(i)=((Pc_eq(0)-Pig)/n)*(i-1)+Pigend dodo i=1,n,1Cf=0.0Cf=ga*sqrt(2*k*(1-(101325.0*Pe/P1_eq(i))**((k-1)/k))/(k-1))F1(i)=Cf*P1_eq(i)*Atend dodo i=1,n,1time1(i)=Apa(0)*len/(Ctz*ga*ga*At*(1-pn))*&log((Pc_eq(0)-Pig)/(Pc_eq(0)-P1_eq(i)))end do!-----------------------------------------------------------!计?算?加è上?起e始?段?时骸?工¤作痢?段?的?时骸?间?timedo i=0,num,1time(i)=time1(n)+time(i)end do!-----------------------------------------------------------!计?算?拖?尾2段?的?压1力 和æ时骸?间?值μ(按恪?照?绝?热â?膨ç胀æ过y程å计?算?,不?考?虑?有瓺剩骸?药?作痢?用?)!-----------------------------------------------------------do i=0,n-1,1P2_eq(i)=Pc_eq(num)-i*(Pc_eq(num)-101325.0)/(n-1) end dodo i=0,n-1,1Cf=0.0Cf=ga*sqrt(2*k*(1-(101325.0*Pe/P2_eq(i))**((k-1)/k))/(k-1)) F2(i)=Cf*P2_eq(i)*Atend dodo i=0,n-1,1time2(i)=time(num)+2.0*Apa(num)*len*&((P2_eq(i)/Pc_eq(num))**((1-k)/(2.0*k))-1.0)/&((k-1)*ga*At*sqrt(Roc*Toc))end do!将?平?衡a压1力 转羇化ˉ为a物?理え?大洙?气?压1的?形?式?do i=1,n,1P1_eq(i)=P1_eq(i)/101325.0end dodo i=0,num,1Pc_eq(i)=Pc_eq(i)/101325.0end dodo i=0,n-1,1P2_eq(i)=P2_eq(i)/101325.0end do!输?出?压1力 值μopen(10,file="export_inter_P.dat")do i=1,num+n+n-1,1if(i<=n-1) thenwrite(10,"(f15.5,2x,f15.5)") time1(i),P1_eq(i)else if(i>=n.and.i<num+n) thenwrite(10,"(f15.5,2x,f15.5)") time(i-n),Pc_eq(i-n) elsewrite(10,"(f15.5,2x,f15.5)") time2(i-num-n),P2_eq(i-num-n) end ifend doclose(10)!输?出?推?力 值μopen(20,file="export_inter_F.dat")do i=1,num+n+n-1,1if(i<=n-1) thenwrite(20,"(f15.5,2x,f15.5)") time1(i),F1(i)else if(i>=n.and.i<num+n) thenwrite(20,"(f15.5,2x,f15.5)") time(i-n),Feq(i-n)elsewrite(20,"(f15.5,2x,f15.5)") time2(i-num-n),F2(i-num-n) end ifend doclose(20)returnend subroutine!计?算?侵?蚀骸?比括?的?子哩?函ˉ数簓function fun(Sa,Apa,key)implicit noneinteger :: keyreal(kind=8) :: Sa,Apareal(kind=8) :: funif(key==1) thenif(Sa/Apa<=72.9) thenfun=1.0elsefun=1.3128-1.3249e-2*Sa/Apa+&1.5527e-4*(Sa/Apa)**2-4.3868e-7*(Sa/Apa)**3end ifelsefun=1.0end ifreturnend。

火箭发动机原理教学大纲

火箭发动机原理教学大纲

《火箭发动机原理》课程教学大纲课程代码:110132307课程英文名称:Solid Rocket Motor课程总学时:32 讲课:32 实验:0 上机:0适用专业:弹药工程与爆炸技术大纲编写(修订)时间:2017.10一、大纲使用说明(一)课程的地位及教学目标本门课程是弹药工程与爆炸技术专业的一门专业选修课。

固体火箭发动机是卫星、火箭、飞机、导弹等产品的动力装置,它在现代科学技术研究,国民经济的发展,人们日常生活的改善等方面有着很大的利用价值,在本专业中对于火箭、导弹或炮弹增程有着极其重要的作用。

通过本课程的学习,学生将达到以下要求:1.熟练掌握固体火箭发动机的基本结构、工作原理,燃气在喷管与燃烧室内的流动过程,掌握固体火箭发动机内弹道的计算方法。

2.掌握固体火箭发动机的总体结构设计方法。

3.要求学生能将所学知识灵活运用于产品的设计和生活实践当中。

(二)知识、能力及技能方面的基本要求要求学生理解并掌握《火箭发动机原理》这门课程,使学生对固体火箭发动机有一定的认识。

1.掌握固体火箭发动机原理的主要内容,包括固体火箭发动机的工作原理、固体火箭推进剂以及固体火箭推进剂在燃烧室中的燃烧过程、燃气在喷管中的流动过程、固体火箭发动机性能参数、固体火箭发动机的热力计算、固体火箭发动机的内弹道计算方法等方面的知识。

2.掌握固体火箭发动机设计的主要内容,包括固体火箭发动机的基本结构,主要设计参量的选择,发动机结构的初步设计等。

3.了解固体火箭发动机的应用及发展趋势,并能用所学知识指导在本领域的技术研究和产品的设计。

(三)实施说明1.教学方法:课堂讲授中重点对固体火箭发动机的基本概念,工作原理和设计方法进行讲解。

培养学生的思考能力和分析问题的能力。

在讲授中注意采用理论知识与实际应用相结合的方法,提高学生分析问题、解决问题的能力。

2.教学手段:在教学中主要采用电子教案、CAI 课件及多媒体教学系统等教学手段相结合。

火炮内弹道设计-毕业论文资料

火炮内弹道设计-毕业论文资料

1 绪论内弹道(internal ballistics)是弹道的一部分,内弹道研究弹丸从点火到离开发射器身管的行为。

内弹道学研究对各种身管武器都有重要意义。

击发方法:任何类型的身管武器第一步需要击发火药。

最早的枪支、大炮由一个一端密封的金属管组成。

1.1 内弹道学研究对象内弹道学是研究发射过程中枪炮膛内及火箭发动机内的火药燃烧、物质流动、能量转换、弹体运动和其它有关现象及其规律的弹道学分支学科。

燃烧的发射药产生具有很高压力的气体,使弹丸加速穿过炮膛,直到以预定初速离开炮口。

初速是具有一定质量和形状的弹丸最终要达到的整个射程的基础。

在设计火炮时必须进行计算以保证最正常、最有效地产生所需要的初速。

发射装药产生的能量用于完成好几种工作。

大部分能量用于赋予弹丸速度。

能量还消耗在做下述功上:使弹丸旋转,克服弹丸与膛壁之间的摩擦力,使发射药和发射药气体在膛内运动以及使火炮后坐部分后坐。

有些能量还以热能的形式损失在身管、炮尾、弹丸和药筒(如果使用药筒的话)上。

发射过程都是从点火开始,通过机械击发、电热或其他方式将点火药点燃,所产生的高温气体及灼热粒子再点燃火药装药,迅即扩展到整个装药表面,并同时沿着药粒厚度向内层燃烧。

燃烧进行在一个封闭的空间中,这个空间前由弹丸的弹带封闭,后有火炮所采用的紧塞装置封闭,紧塞装置用于防止火药气体从后面逸出。

在发射药气体的压力达到能使弹丸运动的程度之前,发射药的燃烧速度与膛压增加的速度是成正比例的。

所谓“弹丸启动压力”就是指使弹丸开始向前运动的压力。

当弹丸沿身管向前运动时,供发射药气体占用的空间增大,因此膛压的增加速度减小。

当空间增加所导致的压力的增加相等时,膛压达到最大值。

自此以后膛压开始下降,同时弹丸却在继续加速,甚至在发射药全部燃尽后弹丸仍在继续加速,只是加速度逐渐减小,弹丸一出炮口即变为减速。

下图说明膛内压力、弹丸膛内行程和弹丸速度间的关系。

内弹道学的研究对象,主要是有关点火药和火药的热化学性质,点火和火药燃烧的机理及规律;有关枪炮膛内火药燃气与固体药粒的混合流动现象,有关弹带嵌进膛线的受力变形现象,弹丸和枪炮身的运动现象;有关能量转换、传递的热力学现象和火药燃气与膛壁之间的热传导现象等。

固体火箭发动机零维两相内弹道研究

固体火箭发动机零维两相内弹道研究

固体火箭发动机零维两相内弹道研究陈军【摘要】为方便应用两相内弹道流动模型对火箭发动机进行性能预示以及提高性能预示精度,利用火箭喷管内的两相流动性能计算公式,建立了零维两相内弹道模型,包括零维两相内弹道微分方程和平衡压强公式,给出了模型中涉及到的两相特性参数的计算方法.相比于一维两相内弹道模型,该零维内弹道模型简单且满足必要的精度,适于工程应用的快速估算.利用该模型对某远程火箭发动机进行了内弹道计算与分析,计算结果与实验数据吻合良好,表明该两相内弹道模型可以有效地降低纯气相模型引起的理论与实际之间的模型偏差,有利于快速计算固体推进剂火箭发动机的两相内弹道参数以及提高预示精度.%To predict internal ballistic properties with two-phase flow for convenience and to improve prediction precision of internal ballistics in solid rocket motor(SRM),a model of zero-dimensional two-phase internal ballistics was built.The model includes internal ballistic differential equation and corresponding equilibrium pressure formula,and the computational methods of twophase property parameters in the model were pared with one-dimensional twophase internal ballistic model,the model is very simple and more suitable to quick computation in project applications,and more accurate.Based on the model,an internal ballistic computation and analysis for a long-range solid rocket engine were carried out.The computed results accord well with experimental data.The zero-dimensional two-phase internal ballistic model can availably reduce theoretical error caused by single-gas-phase model.By the model,two-phase internal ballisticproperties can be quickly computed,and the prediction precision of internal ballistics in SRM can be improved.【期刊名称】《弹道学报》【年(卷),期】2013(025)002【总页数】5页(P39-43)【关键词】内弹道;两相流;固体火箭发动机;固体推进剂【作者】陈军【作者单位】南京理工大学机械工程学院,南京210094【正文语种】中文【中图分类】V435随着高能推进剂在火箭发动机中的普遍应用,两相流动对发动机性能的影响越来越受到重视。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

1 1 n
M
1 1 n
1 1 n 1 1 M km 1 n p RT0
当Peq较小时,如Peq=9.8MPa时,可取ε=0,则此时有: Peq M
1 1 n
2. 影响平衡压强的因素: 1 * 从平衡压强的表达式 Peq 1- C p aKN ( æ ) / 1-n 可以看到:当 处于理想状态(即流量修正系数ψ=1、热损失修正系数χ=1与ε=0)时, 影响平衡压强的因素主要有两方面:一方面是推进剂的特性,如密度、 特征速度、燃速特性与压强指数;另一方面是发动机几何设计参数,如 面喉比等。 (1)推进剂特性对平衡压强的影响 推进剂特性对平衡压强有着决定性的影响。由于推进剂的种类很多, 设计中可根据要求选择推进剂,确定相应的平衡压强。例如,就常用的 推进剂而言,其密度ρp约为1600-1800kg/m3,其特征速度C*约为12001600m/s,其燃速则差别更大,速燃推进剂可高达每秒几十毫米,缓燃 推进剂只有每秒几毫米,甚至每秒1毫米。所有这些数据都直接影响发 动机的平衡压强。
在发动机实际工作过程中,燃气密度ρ远小于装药密度ρp, 因此,可忽略填充量,故微分方程又可简化为:
Vg dp mb mt RT0 dt
式中 其中
mb p Ab r p Ab apn ()
为沿装药全长的平均侵蚀比。
pAt pAt C* RT0
根据质量守恒原理,燃烧室内燃气生成率 mb 与燃气通过喷管 排出的质量流率 mt 之差应等于燃烧室内燃气质量变化率,即:
dmr m b mt dt
p : 推进剂装药密度
Ab : 装药燃烧面积
m b p rAb 其中: mr Vg
mr : 燃烧室内存留的燃气质量 Vg : 燃气占有的容积, 即燃烧室自由容积
根据对推进剂燃烧特性分析可知,常用燃速温度敏感系数σp来表 达燃速受初温的影响程度。从有关资料查得燃速温度敏感系数σp的值 和已知某初温Ti1下的燃速r1后,可利用如下公式求得任一初温Ti2的燃 速r2,即 r re p (Ti 2 Ti1 ) 2 1 并且,根据压强温度敏感系数的定义,可由初温Ti1下的平衡压强 估算初温Ti2下的平衡压强,近似公式为:
于相对稳定值。对于减面燃烧装药,燃气生成量不断减小,燃烧室压 强也逐渐减小。同样,由于喷管流量的不断下降使燃烧室压强与流量 处于相对稳定值。对于恒面燃烧装药,从燃气生成量与流量达到平衡 以后,由于燃烧表面积不变,燃气生成量与喷管流量的平衡可以一直 维持下去,压强也因而不变,直到整个装药燃烧结束。
3 拖尾阶段(结束段) 此时装药燃烧基本结束,燃气生成量近似 为零,只有燃气的排出。在此阶段,燃烧室内的燃气质量迅速减少, 因而压强迅速下降,直到与外界环境压强相等,排气停止,拖尾阶段 结束。这个压强下降过程又称为“后效”过程,是发动机工作的尾声。 以上三个阶段中,发动机工作阶段是火箭作为运载工具产生推进 动力的主要阶段。在大多数情况下,要求发动机性能相对稳定,燃面 变化尽可能小,尽量采用恒面燃烧,燃烧室压强随时间的变化比较小, 这样,可以作为定常或准定常问题来处理。与此相反,若在压强上升 段和拖尾段中,压强随时间的变化很大,则离准定常的条件更远,会 使问题的处理复杂化。
2.发动机工作阶段(工作段) 当燃烧室内已充满了高压的燃气, 燃气的生成量和喷管流量达到相对平衡,因而压强的变化比较平缓。 在这个阶段中,燃气生成量的变化主要决定于装药燃烧表面积的变化。 对于增面燃烧的装药,燃气生成量随燃面的增大而逐渐增加,燃烧室 压强也逐渐增加。与此同时,喷管流量的增大使燃烧室压强不断地处
机正常和稳定的工作,使推进剂的化学能充分转化为热能,要求燃 烧室压强必须高于推进剂完全燃烧的临界压强;从结构设计方面来 看,燃烧室是一个主要承受内压的部件,在进行各组件和药柱的强 度计算前,必须先确定燃烧室中可能出现的最大压强,其值的大 小,直接影响对燃烧室的强度要求和结构重量。 由此可见,在发动机设计过程中,首先确定推进剂成分,装药 几何尺寸和喷管喉径。计算出燃烧室压强随时间空变化的曲线;然 后求得发动机的推力随时间的变化规律和有关发动机的其它性能参 数以及进行发动机壳体结构设计和强度计算;最后,确定发动机设 计性能。有时,为达到总体设计要求,要反复多次地进行装药和喷 管几何尺寸的设计以及内弹道计算,以求得发动机的最佳设计。 总之,内弹道计算的任务是在确定推进剂成分、装药几何尺 寸、工作环境温度、喷管喉部直径等条件下,计算燃烧室压强随时 间的变化规律。
M 为装填参量
1 1-n 1 1n
km : 单位换算系数, p单位为Pa时,km 1; p单位为MPa时,km 106
则平衡压强为: Peq 1- C p aKN (æ ) / பைடு நூலகம்
*
[M (1 )]
由于ε<<1,则
Peq (1

1 n
)M
固体火箭发动机内弹道计算,即燃烧室压强时间曲线计算,是固 体火箭发动机设计中的一个重要环节。首先,由推力公式 F=CFPcAt 可见,燃烧室压强的变化规律可直接决定火箭弹的推力方案;其次, 对一定的装药来说,燃烧层的厚度是一定的,推进剂燃速受压强控制, 燃烧室压强越高,推进剂燃速越大,装药燃尽时间越短。因此,燃烧 室压强又是决定发动机工作时间的重要因素;此外,为保证发动
第七章 固体火箭发动机的内弹道计算
一、内弹道计算的任务 二、燃烧室压强的变化 三、零维内弹道计算的微分方程 四、平衡压强及其影响因素 五、燃烧室压强—时间曲线的简化计算
一、内弹道计算的任务
内弹道学是从枪炮技术中引来的一个术语。原意是研究发射过程 中弹丸在膛内的运动和膛内压强的变化。固体火箭发动机内弹道学的 核心是研究发动机燃烧室内压强随时间变化的规律。因此,内弹道计 算就是计算燃烧室内燃气压强—时间曲线,其最终目的是计算发动机 推力—时间曲线和质量流率—时间曲线,为火箭、导弹外弹道计算提 供依据。
在上面的分析中,认为燃烧室是一个充满高压燃烧气体的容器, 不考虑燃气的流动和燃烧室内的压强分布,室内各点的压强都相等。 这样,整个燃烧室压强同时随时间变化,与该点的位置坐标x无关,这 就是所谓“零维”的压强变化。对于燃气流速很小的燃烧室来说,压 强计算可以看作是一个“零维”问题来处理。但是,对装填密度较大 的侧面燃烧装药,燃气在通道中的流动沿轴向产生很大的速度,因此, 压强沿轴向有显著的变化。这种情况下,必须考虑压强在燃烧室中的 分布,应作为“一维”问题来进行压强计算。
Peq = 1 p p Ab a () RT0 At km
M

P k eq m p p RT0
p K N a ( æ ) RT0
1 1-n
km
其中: K N 为面喉比, K N Ab / At
引入流量修正系数ψ和热损失修正系数χ,则:
mt
其中C * : 特征速度,C *
RT0

于是又可得到计算p-t曲线的微分方程的又一形式:
pAt dp p Ab ap n () RT0 dt RT0
Vg
四、平衡压强及其影响因素
1.平衡压强的定义:由上面分析可知,当燃气每秒净增率等于质 量流率时,即: m mt时,dp/dt 0时,燃气处于平衡状态,对应的压 当 强称为平衡压强,记为Peq.
Vg
d dt
表示单位时间内改变燃气密度所需的燃气质量。
根据基本假设及完全气体的状态方程 p RT 可有:
Vg d Vg dp dt RT0 dt
整理后有:
Vg dp Ab r ( p ) mt RT0 dt
此即计算零维内弹道压强—时间(p-t)曲线的微分方程。
1.发动机起动阶段(上升段) 这包括点火和压强建立过程。首先 依靠点火装置中点火药点燃并燃烧生成的高温气体充满燃烧室,一方 面使燃烧室压强上升到点火压强;另一方面加热推进剂表面,点燃主 装药,这就是点火过程。当主装药全面点燃后,燃气质量生成量迅速 增大,并在瞬时超过喷管的质量流量,使燃烧室的压强迅速增加,同 时又促使喷管流量的增加,不断地与燃气生成量趋于相对平衡。最后, 燃烧室压强达到其相对稳定值,这个相对稳定值的压强称为工作压强。 这个压强建立的过程即称为发动机启动阶段。对一般发动机来说,这 个过程在几十毫秒内完成。
二、燃烧室压强的变化
在发动机工作过程中,一方面推进剂装药燃烧,不断生成燃气, 充填燃烧室自由容积。燃气生成的速度按每秒生成多少质量来计量, 并称之为燃气的质量生成量;另一方面,燃气经过喷管流出,以每秒 流出多少质量来计量,这称之为喷管的质量流量。如果燃气生成量超 过喷管的质量流量,燃烧室自由容积内的燃气质量不断积累而使压强 上升。反之,如果燃气生成量下降,低于喷管的质量流量,则燃烧室 压强就会下降。只有在一定的条件下,燃气生成量与喷管质量流量达 到相对平衡时,压强也达到相对稳定值。 由发动机实验所测得的 燃烧室压强一时间曲线可见, 燃烧室压强的变化有三个阶 段,如右图所示:
令 m Ab r ( p ) ,并称之为燃气每秒净增率,则燃烧室内
燃气压强变化情况取决于燃气每秒净增率与质量流率的大小。前者 是使压强增大的因素,后者是使压强减小的因素。
当 m mt时,dp/dt 0, 压强增大; 当 m mt时,dp/dt 0, 压强减小; 当 m mt时,dp/dt 0, 压强不变.
三、零维内弹道计算的微分方程
严格说来,内弹道计算与气流各参数(压强、温度、密度等)沿燃 烧室长度的分布随时间的变化规律有关,应使用一维非定常流的普遍 方程组来确定燃烧室内燃气的压强、温度及密度与坐标x及时间t的函 数关系,但这种方法在数学上会遇到一系列的困难。因此,在工程计 算中,通常简化为零维问题来处理。 零维内弹道计算的基本假设如下: • 燃气流动参数取其沿轴向的平均值; • 推进剂装药燃烧完全,燃烧产物组分不变,且燃烧温度等于推进 剂的等压燃烧温度; • • 燃气为完全气体,服从完全气体状态方程; 装药燃烧服从几何燃烧定律。
相关文档
最新文档