固体火箭发动机燃烧室与内弹道计算
龙格库塔法计算固体火箭发动机内弹道
!计?算?星?型³装痢?药?的?几?何?尺?寸?!参?数簓符?号?说μ明¶!-----------------------------------------------------------------------------------!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 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**2write(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)!调獭?用?子哩?程²序´计?算?装痢?药?内¸弹獭?道台?曲¸线?callinternal_ballistics0(d,len,e1,n,thet,epsilon,r,r1,l,a,pn,Ctz,At,k,density_p,mpeff,Isp) !输?出?星?型³装痢?药?的?几?何?参?数簓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!该?子哩?程²序´用?于 ?计?算?零?维?变?截?面?燃?烧?装痢?药?的?内¸弹獭?道台?!利?用?龙ⅷ?格?-库a塔t法ぁ?计?算?内¸弹獭?道台?曲¸线?!可°用?于 ?计?算?侵?蚀骸?燃?烧?效§应畖下?的?内¸弹獭?道台?曲¸线?!--------------------------------------------------------------------------!d外猘径?,len长¤度¯,e1平?均·肉╝厚?,n星?角?数簓,thet星?边?夹D角?,epsilon角?度¯系μ数簓!r过y度¯圆2弧?半?径?,r1星?角?圆2弧?半?径?和³l特?征¶尺?寸?,key表括?示?是?否?考?虑?侵?蚀骸?燃?烧?!density_p推?进?剂®的?密¹度¯,k比括?热¯?比括?ga系μ数簓,C特?征¶速·度¯,mpeff药?柱·质±量?!a速·度¯系μ数簓,pn压1强?指?数簓,At喉³部?面?积y,Pc燃?烧?室酣?的?设Θ?计?压1力 ,Isp表括?示?理え?论?比括?冲?!dt时骸?间?步?长¤,e燃?层?厚?度¯,time时骸?间?,Sa同?一?时骸?刻²轴®向´的?各¶节¸点?的?燃?面?,Apa通 ?气?面?积y!Poc各¶节¸点?的?压1力 ,Pe喷?管¹出?口¸压1力 ,P_I压1力 冲?量?,ep侵?蚀骸?比括?!P_av平?均·压1力 ,F_av平?均·推?力 ,I0总哩?冲?,Im重?量?比括?冲?Iv体?积y比括?冲?,Cf推?力 系μ数簓!--------------------------------------------------------------------------subroutineinternal_ballistics0(d,len,e1,n,thet,epsilon,r,r1,l,a,pn,Ctz,At,k,density_p,mpeff,Isp)implicit nonereal,parameter::Pi=3.14integer :: n,keyreal(kind=8),intent(in) :: d,len,e1,thet,epsilon,r,r1,lreal(kind=8),intent(in) :: density_p,k,Ctz,mpeff,Ispreal(kind=8) :: a,pn,At,dt,e,time,Sa,Apa,Poc,Pe,P_I,F real(kind=8) :: ep,P_av,F_av,I0,Im,Iv,Cf,gareal(kind=8) :: f1,f2,f3,f4Pe=101325.0P_I=0.0!药?柱·的?能¹量?特?性?参?量?读®入?open(3,file="inter_parameter.dat")read(3,*) keyread(3,*) dtread(3,*) Pocclose(3)!计?算?点?火e压1强?时骸?的?推?力 值μPoc=Poc*101325.0ga=sqrt(k)*(2/(k+1))**((k+1)/(k-1)/2)Cf=ga*sqrt(2*k*(1-(Pe/Poc)**((k-1)/k))/(k-1))F=Cf*Poc*Ate=0.0time=0.0open(10,file="export_Poc.dat")open(20,file="export_F.dat")write(10,"(f8.5,2X,f10.5)") time,Poc/101325.0write(20,"(f8.5,2X,f15.5)") time,Fdo while(.true.)!计?算?燃?面?和³通 ?气?面?积yif(e<=r1) thenSa=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-e)*Pi/n/l)*l*len Apa=(n*((1-epsilon)*Pi/n+sin(epsilon*Pi/n)*(cos(epsilon*Pi/n)-&sin(epsilon*Pi/n)*cotan(thet/2)))+2*n*(r+e)*(sin(epsilon*Pi/n)/&sin(thet/2)+(1-epsilon)*Pi/n)/l+n*(r+e)**2*(Pi/n+Pi/2-thet/2-&cotan(thet/2))/l**2+n*(r1-e)**2*(thet/2+cotan(thet/2)-Pi/2)/l**2)*l**2else if(e>r1.and.e<=(l*sin(epsilon*Pi/n)/cos(thet/2)-r)) thenSa=2*n*(sin(epsilon*Pi/n)/sin(thet/2)+(1-epsilon)*Pi/n+&(e+r)*(Pi/n+Pi/2-thet/2-cotan(thet/2))/l)*l*lenApa=(n*((1-epsilon)*Pi/n+sin(epsilon*Pi/n)*(cos(epsilon*Pi/n)-&sin(epsilon*Pi/n)*cotan(thet/2)))+2*n*(r+e)*(sin(epsilon*Pi/n)/&sin(thet/2)+(1-epsilon)*Pi/n)/l+n*(r+e)**2*(Pi/n+Pi/2-thet/2-cotan(thet/2))/l**2)*l**2 else if(e>(l*sin(epsilon*Pi/n)/cos(thet/2)-r).and.e<=e1) thenSa=2*n*((1-epsilon)*Pi/n+(r+e)*(Pi/n+asin(l*sin(epsilon*Pi/n)/(e+r)))/l)*l*len Apa=n*((1-epsilon)*Pi*(1+(r+e)/l)**2/n+sin(epsilon*Pi/n)*(sqrt((r+e)**2/l**2-& sin(epsilon*Pi/n)**2)+cos(epsilon*Pi/n))+(r+e)**2*&(epsilon*Pi/n+asin(l*sin(epsilon*Pi/n)/(r+e)))/l**2)*l**2elseSa=0.0end if!计?算?侵?蚀骸?比括?if(key==1) thenif(Sa/Apa<=72.9) thenep=1.0elseep=1.3128-1.3249e-2*Sa/Apa+1.5527e-4*(Sa/Apa)**2-4.3868e-7*(Sa/Apa)**3end ifelseep=1.0end if!采°用?4阶¬的?龙ⅷ?格?库a塔t法ぁ?计?算?f1=ga**2*Ctz**2*(density_p*Sa*ep*a*Poc**pn-Poc*At/Ctz)/(Apa*len)f2=ga**2*Ctz**2*(density_p*Sa*ep*a*(Poc+dt*f1/2)**pn-(Poc+dt*f1/2)*At/Ctz)/(Apa*len )f3=ga**2*Ctz**2*(density_p*Sa*ep*a*(Poc+dt*f2/2)**pn-(Poc+dt*f2/2)*At/Ctz)/(Apa*len )f4=ga**2*Ctz**2*(density_p*Sa*ep*a*(Poc+dt*f3)**pn-(Poc+dt*f3)*At/Ctz)/(Apa*len) Poc=Poc+dt*(f1+2.0*f2+2.0*f3+f4)/6.0Cf=ga*sqrt(2*k*(1-(Pe/Poc)**((k-1)/k))/(k-1))F=Cf*Poc*Attime=time+dte=e+dt*ep*a*Poc**pnif(e>e1.and.Poc<101325.0) exitwrite(10,"(f8.5,2X,f10.5)") time,Poc/101325.0write(20,"(f8.5,2X,f15.5)") time,FP_I=P_I+Poc*dtend doclose(10)close(20)!循-环«结®束?P_av=P_I/timeCf=ga*sqrt(2*k*(1-(Pe/P_av)**((k-1)/k))/(k-1))F_av=Cf*P_av*AtI0=F_av*timeIm=I0/mpeff/9.81Iv=I0/(Apa*len)!输?出?固²体?火e箭y发ぁ?动ˉ机¸的?工¤作痢?参?数簓值μopen(30,file="export_propellant.dat")write(30,"('平?均·压1力 p_av=',f9.5)") p_av/101325.0write(30,"('平?均·推?力 F_av=',f10.5,1X,'KN')") F_av/1000.0 write(30,"('总哩?冲?I0=',f10.5,1X,'KN.s')") I0/1000.0write(30,"('重?量?比括?冲?Im=',f10.5,1X,'s')") Imwrite(30,"('体?积y比括?冲?Iv=',f10.5,1X,'KN.s/m3')") Iv/1000.0 write(30,"('工¤作痢?时骸?间?time=',f8.5,1X,'s')") timeclose(30)returnend subroutine。
三维药柱固体火箭发动机内弹道数值模拟_程雪玲
n
∑ Aeb =
A
i b
+Aeb ⊥ (e
=1
,2,
… , k)
(4 )
i =1
式中 Aeb⊥为垂直于 Z 轴的燃面面积 。
∑ Aeb ⊥
=
1 2
m
(ΔY ⊥, j
j =1
+ΔY
⊥, j +1)×ΔXj
(j =1 , 2 , … , m ) (5 )
式中 ΔY ⊥, j 和 ΔY ⊥, j +1为燃层厚度内的 Y 向线
1
段
p0 =
aρg rc *
Ab At
1 -n
其它时刻的头部压强 pmc 为
rρg r A b
p
m c
=
p
m c
-1
+ t
1-
p
m c
-1
A
c*
t
(j 12 )
×
RT V
b
Δt
(13 )
由于上升段的压强对时间的斜率较大 , 在计算
头部压强时采用下式 :
p
m c
=Hale Waihona Puke peqRT Vb=
1
s-的燃1 -
段长度 , 如图 2 所示 。
图 2 翼柱型药柱横截面图
F ig .2 Cross section diagram of fin-cylindric g rain in fin-slot section
2 .2 分段解析法计算一维内弹道 采用分段解析法计算内弹道 , 视燃烧产物的流
动在微元时间内为定常流 , 而在两段微元时间之间 为非定常流 , 即视燃烧产物为准定常流的分段解析
固体火箭发动机内弹道学 方丁酉
固体火箭发动机内弹道学方丁酉固体火箭发动机内弹道学方丁酉近年来,随着航天事业的蓬勃发展,固体火箭发动机内弹道学成为越来越受关注的领域。
方丁酉,中国工程院院士、火箭技术专家,曾在这一领域做出了卓越的贡献。
本文将从简单到复杂、由浅入深地探讨固体火箭发动机内弹道学,以帮助读者全面、深入地理解这一话题。
一、固体火箭发动机内弹道学的概念固体火箭发动机内弹道学是研究固体火箭发动机内燃烧过程、燃料燃烧特性、燃烧产物排放、工作环境等方面的学科。
方丁酉指出,固体火箭发动机内弹道学是固体火箭技术中的重要基础理论,对于提高固体火箭发动机的性能、可靠性和运载能力具有重要意义。
二、固体火箭发动机内弹道学的理论基础在固体火箭发动机内弹道学的研究中,燃烧动力学是重要的理论基础之一。
燃烧动力学研究了燃料在燃烧过程中的变化规律,以及燃烧反应对火箭发动机内部气体流动和压力变化的影响。
燃烧产物排放和燃烧室内部气体动力学也是固体火箭发动机内弹道学的重要内容之一。
研究这些理论基础可以帮助我们更好地理解固体火箭发动机内部的工作原理和特性。
三、固体火箭发动机内弹道学的关键技术在固体火箭发动机内弹道学研究中,燃烧稳定性和效率是两个关键技术。
燃烧稳定性是指在固体火箭发动机工作过程中保持燃烧的稳定性,避免出现燃烧不均匀或燃烧失稳等问题。
而燃烧效率则是指在燃烧过程中尽可能提高燃料的利用率,减少燃料的浪费。
方丁酉在固体火箭发动机内弹道学的研究中,提出了一系列有效的技术方案,使固体火箭发动机在燃烧稳定性和效率方面取得了显著的进展。
四、固体火箭发动机内弹道学的应用前景固体火箭发动机内弹道学的研究成果已经在我国的航天事业中得到了广泛的应用。
在长征系列火箭、嫦娥探月工程等多个航天工程中,固体火箭发动机内弹道学的研究成果为提高火箭的性能和可靠性做出了重要贡献。
未来,随着我国航天事业的不断发展,固体火箭发动机内弹道学将继续发挥重要作用,推动我国航天事业迈向新的高度。
论业余模型火箭发动机设计方法5.31(1)
CH4
2.91E-07
2.77E-07
4.03E-07
CO
1.91E-01
1.85E-01
1.62E-01
CO2
1.54E-01
1.57E-01
1.75E-01
H
4.34E-06
2.38E-06
燃料名称 氧化剂 粘结剂
理论最大比冲
KNSU KNO3
蔗糖
137
KNSB KNO3 山梨醇
164
KNDX KNO3 葡萄糖
160
燃速数据如下
表2
图3
对于上述的燃料,可通过 Cprepop 软件进行相关参数的求解,下图为 Cprepop 软件界面
4
中青为科技小组
论业余模型火箭发动机设计方法
图4
下面以 KNSU 为例,在软件界面中输入相应的配比,这里硝酸钾对蔗糖之比为取 13:7。 可以计算化学平衡流与化学冻结流下的燃料产物组分与相关参数,如燃烧温度 Tf ,燃
固体火箭发动机的特点是通过消耗自身携带的燃料,向后排出高速运动的工质推动 载荷进行运动。一般而言,模型火箭要求其发动机制造价格低廉,结构简单,可快速重 新装填入模型火箭中进行发射作业。故现对一般的固体火箭发动机结构给予适当简化。 模型火箭发动机典型工作时间通常在 1.0s 左右,所以对于火箭发动机内壁的热防护可以 基本不必考虑,隔热层并没有安装的必要。对于市面上大多数模型固体火箭发动机而言, 通常都采用简单的圆孔尾喷管代替一般固体火箭发动机的拉法尔喷管;采用简单的深/ 浅内孔燃烧的管状装药以方便批量生产。可以说,模型火箭牺牲了一部分性能满足了上 述的价格低廉,结构简单的特点。图 1 所示为一般模型火箭发动机结构:
2 模型固体火箭发动机一般设计流程
固体火箭发动机燃烧室三维流动数值计算
固体火箭发动机燃烧室三维流动数值计算近年来,固体火箭发动机燃烧室三维流动数值计算已成为航天领域
的重要研究方向之一。
在这个背景下,本文通过对一个典型的固体火
箭发动机燃烧室进行了三维流动数值计算,并分析了其内部气体流场
的特点和参数变化规律。
具体而言,本文采用了基于有限元方法和RANS模型相结合的求解策略,在考虑边界条件、壁面摩擦等影响因
素后,得到了该火箭发动机燃烧室内部压力、温度、速度等关键物理
量随时间推移的变化情况,并对其进行了详尽分析和讨论。
根据结果
显示,该固体火箭发动机燃烧室内存在着复杂多变且高强度的湍流现象,在不同位置处形成了不同类型及大小尺度上占优势地位的涡旋结构。
同时,在整个过程中气体温度呈现出明显非均匀性分布,并随着
时间逐渐上升;而压力则表现出先增大后减小再稳定下降趋势。
此外,在不同位置处还存在着明显差异且互相作用影响较大的气体速度场。
总之,本文通过对固体火箭发动机燃烧室三维流动数值计算与分析,
深入揭示并探究了其内部气态流场特性及参数变化规律。
这些结果对
进一步提高航天技术水平以及设计更加可靠安全有效率的固体火箭发
动机具有重要意义和参考价值。
固体推进剂能量计算方法
固体推进剂能量计算方法一 固体推进剂能量计算原理 1,基本假设在火箭发动机工作时,固体推进剂的化学潜能转换为燃气的动能,经历了推进剂燃烧和燃烧产物膨胀两个过程。
发动机的实际工作过程是非常复杂的。
其复杂性在于:由于存在热损失,难以保证燃烧过程是等压绝热的;燃烧产物在燃烧室内分布是不均匀的;对于含铝、含镁、含硼推进剂或含有某些金属化合物的性能添加剂的推进剂,存在凝聚相产物,这些凝相产物在喷管膨胀过程中导致两相流损失;喷管流动难以保证等熵条件等等。
为了反映固体推进剂能量转换过程的本质,抓住主要矛盾,在进行其理论性能预估时,进行了一些基本假设。
(1) 在燃烧室中,推进剂的燃烧反应达到化学平衡,且燃烧过程为等压绝热过程,即热力学中的等焓过程;而且燃烧产物的分布是均匀的。
(2) 燃气为理想气体,凝相产物的体积忽略不计。
(3) 喷管中燃气的流动过程为绝热可逆过程,即为等熵过程;燃气在喷管中的流动为一维定常流,即在喷管的任一截面上,燃气的组成及各性能参数的分布是均匀的。
(4) 不考虑凝聚相燃烧产物的两相流损失。
2, 基本方程 (1) 质量守恒方程常见的固体推进剂是由C 、H 、O 、N 、Cl 、Al 等元素构成的某些化学物质的混合物。
对于这样一个复杂的系统,假设固体推进剂的燃烧产物共有n 种,而固体推进剂所含有的元素共l 种。
对j 元素的质量守恒方程可表达成:()11,2,,nij ij i a xb j l ===⋅⋅⋅∑ (1)式中,ij a 为混合物系中第i 种产物含j 种元素的原子摩尔数,它由i 燃烧产物的分子式得到;i x 为单位质量燃烧产物中第i 种产物的摩尔数;j b 为单位质量推进剂中含第j 种元素的原子摩尔数,它由推进剂的假想化学式得到。
(2) 能量守恒方程根据假设(1),燃烧室内燃烧为等焓过程,则有p c H H = (2)式中,p H 为单位质量推进剂在初温0T 时的总焓(通常取0298T K =); c H 为单位质量推进剂燃烧产物在平衡火焰温度c T 下的总焓。
基于Pro/E软件进行固体火箭发动机内弹道计算的方法初探
基于Pro/E软件进行固体火箭发动机内弹道计算的方法初探作者:王磊代义李宏来源:《中国新技术新产品》2009年第10期摘要:固体火箭发动机内弹道计算给出的推力、质量流量随时间变化的曲线一直是导弹或卫星外弹道设计的计算依据。
为适应需求,固体火箭发动机装药的几何形状需要反复调整,复杂的几何形状使内弹道计算成为一项繁琐的工作。
本文利用Pro/E软件强大的三维建模能力,通过参数化设计模拟发动机燃烧的全过程,编制简单程序实现发动机零维内弹道计算,为繁琐的计算提供了一条便捷的途径。
关键词:固体火箭发动机;内弹道计算;Pro/E软件1 引言固体火箭发动机是一种采用固体推进剂的化学火箭动力装置,在航天技术中有着广泛的应用。
内弹道计算是固体火箭发动机设计的一项重要内容,它提供的推力、流量等性能参数是导弹总体进行射程计算、散布分析及导弹姿态控制的重要依据[1]。
在满足内弹道特性和保证药柱结构完整性的前提下,设计出具有较高装填分数的药柱对于提高发动机性能具有重要的意义。
药柱几何形状的演变过程,是从简单管状药柱演变到贴壁浇注内孔(星形、车轮形等)燃烧药柱,再到具有较高装填特性的翼柱、锥柱等三维药柱。
内弹道计算的关键在于得出药柱燃面随时间的变化曲线。
有着简单几何形状的药柱,计算燃面都有较为成熟的方法。
为了探索三维药柱燃面的计算方法,国内外有关学者开展了许多工作。
目前,国外内弹道计算最为著名的方法是美国NASA的SPP法,国内较常用的是通用坐标法、有限元素法和边界坐标法。
国内的这三种方法均为数值法,由于使用原函数对复杂药柱进行精确定义很困难,通常难以精确表示燃烧过程中十分复杂的图形,以及连接各部分药柱圆角的复杂变化,所以内弹道计算的精度较差。
同时,它们的计算结果对药柱网格划分的依赖性很强,网格划分不合适,所得结果也不正确[2]。
这些方法计算得到的燃面波动性也较大,药柱初始设计的修改很复杂,数据前后处理也抽象复杂。
这些因素导致内弹道的计算成为一项繁琐复杂的工作。
固体推进剂能量计算方法
固体推进剂能量计算方法一 固体推进剂能量计算原理 1,基本假设在火箭发动机工作时,固体推进剂的化学潜能转换为燃气的动能,经历了推进剂燃烧和燃烧产物膨胀两个过程。
发动机的实际工作过程是非常复杂的。
其复杂性在于:由于存在热损失,难以保证燃烧过程是等压绝热的;燃烧产物在燃烧室内分布是不均匀的;对于含铝、含镁、含硼推进剂或含有某些金属化合物的性能添加剂的推进剂,存在凝聚相产物,这些凝相产物在喷管膨胀过程中导致两相流损失;喷管流动难以保证等熵条件等等。
为了反映固体推进剂能量转换过程的本质,抓住主要矛盾,在进行其理论性能预估时,进行了一些基本假设。
(1) 在燃烧室中,推进剂的燃烧反应达到化学平衡,且燃烧过程为等压绝热过程,即热力学中的等焓过程;而且燃烧产物的分布是均匀的。
(2) 燃气为理想气体,凝相产物的体积忽略不计。
(3) 喷管中燃气的流动过程为绝热可逆过程,即为等熵过程;燃气在喷管中的流动为一维定常流,即在喷管的任一截面上,燃气的组成及各性能参数的分布是均匀的。
(4) 不考虑凝聚相燃烧产物的两相流损失。
2, 基本方程 (1) 质量守恒方程常见的固体推进剂是由C 、H 、O 、N 、Cl 、Al 等元素构成的某些化学物质的混合物。
对于这样一个复杂的系统,假设固体推进剂的燃烧产物共有n 种,而固体推进剂所含有的元素共l 种。
对j 元素的质量守恒方程可表达成:()11,2,,nij ij i a xb j l ===⋅⋅⋅∑ (1)式中,ij a 为混合物系中第i 种产物含j 种元素的原子摩尔数,它由i 燃烧产物的分子式得到;i x 为单位质量燃烧产物中第i 种产物的摩尔数;j b 为单位质量推进剂中含第j 种元素的原子摩尔数,它由推进剂的假想化学式得到。
(2) 能量守恒方程根据假设(1),燃烧室内燃烧为等焓过程,则有p c H H = (2)式中,p H 为单位质量推进剂在初温0T 时的总焓(通常取0298T K =); c H 为单位质量推进剂燃烧产物在平衡火焰温度c T 下的总焓。
固体火箭发动机零维两相内弹道研究
固体火箭发动机零维两相内弹道研究陈军【摘要】为方便应用两相内弹道流动模型对火箭发动机进行性能预示以及提高性能预示精度,利用火箭喷管内的两相流动性能计算公式,建立了零维两相内弹道模型,包括零维两相内弹道微分方程和平衡压强公式,给出了模型中涉及到的两相特性参数的计算方法.相比于一维两相内弹道模型,该零维内弹道模型简单且满足必要的精度,适于工程应用的快速估算.利用该模型对某远程火箭发动机进行了内弹道计算与分析,计算结果与实验数据吻合良好,表明该两相内弹道模型可以有效地降低纯气相模型引起的理论与实际之间的模型偏差,有利于快速计算固体推进剂火箭发动机的两相内弹道参数以及提高预示精度.%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随着高能推进剂在火箭发动机中的普遍应用,两相流动对发动机性能的影响越来越受到重视。
《火箭发动机》 7 内弹道
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
高速旋转固体火箭发动机内流场数值计算
高速旋转固体火箭发动机内流场数值计算郝雯;封锋;罗盟;谢爱元【摘要】为计算高速旋转固体火箭发动机内流场特性,采用RNG k-ε湍流模型,对某端面和内孔同时燃烧管状装药固体火箭发动机燃烧室-喷管统一内流场进行了三维流场数值计算.采用用户自定义函数UDF编程给出了质量入口边界.计算结果表明,高速旋转固体火箭发动机内流场切向速度分布较复杂,在燃烧室前封头呈现出了准Rankie组合涡分布,在喷管内呈现明显的准强迫涡特性.燃气切向速度在燃烧室和喷管喉部分别达到峰值.燃烧室前封头处燃气轴向速度变化剧烈,出现了2个峰值,并在近轴线处出现了轴向速度的负值.高速旋转对前封头与喷喉结构强度与热防护带来严重影响.【期刊名称】《固体火箭技术》【年(卷),期】2013(036)004【总页数】5页(P481-485)【关键词】高速旋转;固体火箭发动机;组合涡;切向速度;轴向速度【作者】郝雯;封锋;罗盟;谢爱元【作者单位】南京理工大学机械工程学院,南京210094;南京理工大学机械工程学院,南京210094;南京理工大学机械工程学院,南京210094;南京理工大学机械工程学院,南京210094【正文语种】中文【中图分类】V4350 引言火箭增程及底排-火箭复合增程是当今炮射弹药增程的两种主要方式。
旋转稳定型弹药飞行时的转速可达到20 000 r/min或以上。
这种高速旋转的工作条件显著改变发动机的装药燃烧规律[1]及发动机内弹道参数[2],主要表现在缩短固体火箭发动机的燃烧时间、增加推力与压强、降低比冲[3-5]。
此外,高速旋转所致的强旋流动现象对发动机结构的热防护带来负面效应[6]。
文献[7]中,采用某固体火箭发动机进行高速旋转实验研究,最后拟合出了动态燃速随装药燃面相对半径和转速变化的经验式。
本文参照此文献中拟合的燃速公式,在Fluent软件基础上,使用用户自定义函数UDF,给出质量入口边界,对端面和内孔同时燃烧的管状装药固体火箭发动机进行了三维内流场数值计算。
固体火箭发动机内弹道学 方丁酉
一、引言固体火箭发动机内弹道学方丁酉,是一门深奥而又具有重要意义的科学。
固体火箭发动机是航天领域中不可或缺的一部分,其内弹道学方丁酉的研究对于提高火箭发动机的性能和可靠性具有重要意义。
在本文中,将深入探讨固体火箭发动机内弹道学方丁酉的重要性、原理和应用。
二、固体火箭发动机内弹道学理论的重要性固体火箭发动机内弹道学方丁酉,是对固体火箭发动机内部燃烧过程的运动状态、燃烧特性、瞬态过程等进行研究的学科。
通过对固体火箭发动机内弹道学方丁酉的研究,可以深入了解燃烧室内各种参数的变化规律,确定燃烧室内流场的性质,以及改善发动机的设计和性能等方面提供重要参考。
固体火箭发动机内弹道学方丁酉的研究对提高火箭发动机的推进效率、降低燃料消耗、提高火箭发射成功率等具有重要意义。
通过对固体火箭发动机内弹道学方丁酉的研究,可以进一步提高火箭的响应速度和控制精度,提高火箭的战场生存能力和作战效能。
固体火箭发动机内弹道学方丁酉的研究内容主要包括燃烧室内燃气运动的多相流动、燃烧过程的热力学性质、燃烧产物的组成及其对流场和燃烧室壁面的影响等。
通过对这些内容的研究,可以从微观和宏观两个层面全面了解固体火箭发动机内弹道学方丁酉的运行机理,进而进行针对性的改进和应用。
固体火箭发动机内弹道学方丁酉的原理在航天领域中有着广泛的应用。
在燃烧室内燃气运动的多相流动方面,可以通过优化燃烧室的结构和燃烧室内部的燃烧过程,提高火箭的推进效率;在燃烧产物的组成及其对流场和燃烧室壁面的影响方面,可以减少火箭发动机的燃料消耗,提高火箭的响应速度和控制精度等。
四、固体火箭发动机内弹道学方丁酉的个人观点对于固体火箭发动机内弹道学方丁酉这一学科,我个人认为其研究的深入和广泛应用,对于航天领域的发展和火箭技术的提升有着极其重要的作用。
只有通过对固体火箭发动机内弹道学方丁酉的深入研究,才能够更好地理解火箭发动机内部燃烧过程的运动状态、燃烧特性、瞬态过程等,从而更好地改进火箭的设计和性能,提高火箭的推进效率、降低燃料消耗、提高火箭发射成功率等。
固体火箭发动机设计
word文档下载后可任意复制编辑第1章绪论1.1设计背景固体火箭发动机与液体火箭发动机和其他化学能火箭发动机相比,具有很多的优点,因而它被广泛的用作各类小型、近程的军用火箭和战术导弹的动力装置。
近几十年来,由于高能推进剂的出现,先进的装药设计和大型药柱浇注工艺的采用,优异的壳体材料和耐烧蚀材料的问世,以及高效而可靠的推力矢量控制装置的研制成功,已在很大程度上克服了固体火箭发动机的缺点,更由于其结构简单,使它在竞争中显示更加优势的地位。
目前,固体火箭发动机除了用于军事用途外,也用于其他的很多方向。
研制和使用新型的高能推进剂,进一步提高推进剂的综合性能,发展无烟推进剂是火箭推进技术主要的研究和发展方向。
总之,随着固体推进技术在航天领域和导弹技术中应用不断发展,会有更多的新课题出现,许多技术问题有待开发。
所以,对固体火箭发动机的研究有十分重要的意思。
1.2固体火箭发动机简介1.2.1 固体火箭发动机基本结构固体火箭发动机主要由固体推进剂、燃烧室、喷管和点火装置等四大部分组成。
图1.1为固体火箭发动机示意图。
1、推进剂装药装药是装入燃烧室中的具有一定形状和尺寸的推进剂药柱的总称,它是固体火箭发动机的能源。
由于装药的燃烧,化学能转化为动能,并且向外做工功,从而推动发动机的运动。
常用的固体推进剂有三类:双基推进剂、复合推进剂word文档下载后可任意复制编辑和改性双基推进剂。
固体推进剂包含有燃烧剂和氧化剂,它自身能够形成封闭的化学反应系统。
2、燃烧室燃烧室里面装载了固体推进剂,是发生化学反应的场所。
它主要由起支承作用的燃烧室壳体和起热防护作用的内绝热层组成,而燃烧室壳体一般由筒体和前后封头组成。
大部分燃烧室都制作成圆柱形,他是主要的受力场所。
燃烧室材料大多采用强度很高的材料,也有采用玻璃纤维缠绕加树脂成型的玻璃钢结构,以大幅度减轻燃烧室壳体的重量。
1——药柱;2——燃烧室;3——喷管;4——点火装置。
图1.1 固体火箭发动机示意图3、喷管在喷管里气流的势能转化为动能,从而使气流加速流动,并保持一定的燃烧室压力,它主要由壳体和热防护层组成。
固体火箭发动机燃烧室及内弹道计算
(A A V dx )(V dx ) 2 dt x x A ( pA)dxdt p dxdt x x
整理后有:
A ( AV ) ( pA AV 2 ) p t x x
燃气的流速远小于当地音速
p 2 tb
推进剂密度远大于燃气密度 流动参数发生显著变化的时间要远 小于装药燃烧时间
以上条件在SRM中是完全可以满足的。
(2) 准定常假设 假设燃气在燃烧室内 流动时,燃气流动的时间
划分为多个时间段,则在
某一时间段内各气流参数 不随时间的变化而变化
燃烧面积 . m . dm
dm p r Ab x
A V (V dx)dt p r b dxdt x x
. . m+dm
整理后有:
A ( A) ( AV ) p r b t x x
dx
(2)动量方程
根据动量守恒:在dt时间内,微元体中动量的变化量=燃气 流动迁移在微元体内产生的动量增量+作用于两截面上的总压力 冲量+作用于微元体侧表面上压力冲量在x轴上的分量,即
2
H0
H
0
2H 0
V
静焓、总焓和流速的变化关系
,总焓和流速的关系如图。
熵:物理学上指热能除以温度所得的商,标志热量转化为功的程度,是表征物质 内部状态的物理量。
滞止状态 V =0 p=p0 T=T0 ρ= ρ0
1
2
V p T ρ
1
静参数与滞止参数 一一对应
2
内侧面燃烧装药发动机示意图
H 01 H 02 H 0 c pT0 根据能量方程可知:
《火箭发动机》 7 内弹道 共21页PPT资料
在上面的分析中,认为燃烧室是一个充满高压燃烧气体的容器, 不考虑燃气的流动和燃烧室内的压强分布,室内各点的压强都相等。 这样,整个燃烧室压强同时随时间变化,与该点的位置坐标x无关,这 就是所谓“零维”的压强变化。对于燃气流速很小的燃烧室来说,压 强计算可以看作是一个“零维”问题来处理。但是,对装填密度较大 的侧面燃烧装药,燃气在通道中的流动沿轴向产生很大的速度,因此, 压强沿轴向有显著的变化。这种情况下,必须考虑压强在燃烧室中的 分布,应作为“一维”问题来进行压强计算。
由发动机实验所测得的 燃烧室压强一时间曲线可见, 燃烧室压强的变化有三个阶 段,如右图所示:
1.发动机起动阶段(上升段) 这包括点火和压强建立过程。首先 依靠点火装置中点火药点燃并燃烧生成的高温气体充满燃烧室,一方 面使燃烧室压强上升到点火压强;另一方面加热推进剂表面,点燃主 装药,这就是点火过程。当主装药全面点燃后,燃气质量生成量迅速 增大,并在瞬时超过喷管的质量流量,使燃烧室的压强迅速增加,同 时又促使喷管流量的增加,不断地与燃气生成量趋于相对平衡。最后, 燃烧室压强达到其相对稳定值,这个相对稳定值的压强称为工作压强。 这个压强建立的过程即称为发动机启动阶段。对一般发动机来说,这 个过程在几十毫秒内完成。
第七章 固体火箭发动机的内弹道计算
一、内弹道计算的任务 二、燃烧室压强的变化 三、零维内弹道计算的微分方程 四、平衡压强及其影响因素 五、燃烧室压强—时间曲线的简化计算
_固体火箭发动机结构
Δc
De ——外径的下偏差值; Di ——内径的上偏差值;
ri
min——强度计算的最小壁厚;
re
re Δc
1 De Di c 2
δ c1
δ c2
ri
D i/2
可能不满足强度要求
D i/2
燃烧室图纸尺寸 ( De图 Di图 ) / 2 c min
根据燃烧室壁厚作强度校核
燃烧室载荷分析 燃气压力 旋转时离心惯性力
运输时振动冲击力
弹道上运动的惯性力
2016/10/9
尾翼式火箭弹燃烧室壳体壁厚计算 计算假设: 忽略外部大气压强 忽略切向惯性力、摆动惯性力以及空气动力和力矩
忽略燃烧室壳体两端轴向力的差异,认为两端拉力相等
壳体为内壁受均布压力的密封容器
2016/10/9
对钢材u=0.3
R 2 0.34Pm
R2 0.34Pm
2
边缘处:
考虑周边并非固支,燃烧室有变形及受热
2
t r
3u R Pm 4 2 3 R r Pm 4
2016/10/9
尾翼式火箭弹燃烧室壳体壁厚计算
(a)按厚壁圆筒
2 ri2 pm r 1 e t 2 应力分布: re ri2 r2 2 ri2 pm r r 2 2 1 e2 re ri r
σx σr σt σr
固体火箭发动机简图(浇注)
1—顶盖;2—点火装置;3—燃烧室壳体;4—药柱; 5—底盖;6—喷管;7—石墨衬套;8—堵盖。
2016/10/9
5.1 燃烧室设计
固体火箭发动机设计
工作时间: tk ≤1.8 秒 工作温度:TH = -40℃~50℃ 燃烧室外径: DH ≤0.426 米,发动机总长: L ≤2.52 米 附加条件:选用一种推进剂,装药采用多根管状药柱,为确
保同温度下的推力在允许范围内,可采用更换喷管。
三、课程设计任务:
1 分析原始条件,确定设计方案,进行发动机结构设计
2.1 发动机结构形式 ............................................................................................ 5 2.2 推进剂及药型选择 ........................................................................................ 5 2.3 壳体材料 ........................................................................................................ 6 2.4 发动机主要设计参数选择 ............................................................................ 7 三、 发动机的装药设计: ..................................................................................... 9 3.1 药柱基本参数 ................................................................................................ 9 3.2 药柱尺寸的确定 .......................................................................................... 10 四、 发动机燃烧室设计 ....................................................................................... 12 4.1 燃烧室结构,材料以及尺寸 ...................................................................... 12 4.2 燃烧室应力分析以及强度校核 .................................................................. 15 五、 发动机喷管设计 ........................................................................................... 15 5.1 喷管的气动设计 .......................................................................................... 15 5.2 喷管壁厚 ...................................................................................................... 17 5.3 喷管的热防护 .............................................................................................. 17 5.4 喷管堵盖 ...................................................................................................... 17 六、 喉部设计及校验计算 ................................................................................... 18 6.1 设计喉部尺寸 dt .......................................................................................... 18 6.2 温度区间—喉衬尺寸 如下表示: ............................................................ 25 七、 内弹道计算 ................................................................................................... 25 7.1 内弹道计算基本方程 .................................................................................. 25
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
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. 流动特点
燃烧室中的燃气流动具有以下特点: 燃烧室中推进剂装药不断燃烧,产生新的燃气加入主流,因此燃 气在燃烧室中的流动是一个有质量加入的流动过程; 在推进剂装药燃烧过程中,燃气的通道面积不断增大,故燃气的 流动参数是位置与时间的函数,因而是非定常流动; 燃气流过形状复杂的通道或流线变化剧烈时,将产生一定的能量 损失; 燃气的流动特性与推进剂的燃烧特性密切相关。例如燃气的压强 和流速影响推进剂的燃速,而燃速又反过来影响燃气压强和流速; 由此可见,在发动机工作过程中,燃气在燃烧室中的流动中极其 复杂的。为了研究方便,特作如下假设: 推进剂燃烧在装药燃烧表面上瞬时完成,燃气的化学组分和热力 性质均匀一致,并沿燃烧表面的外法线方向注入通道。 燃气为完全气体,服从完全气体状态方程,而且比热不变。 燃气在装药通道中作一维流动。 燃气与外界无热和功的交换。
因而,燃气在燃烧室内的流动可视为一维定常绝热加质流动。
2. 控制方程
燃气在燃烧室内流动时,应遵守自然界的普遍规律,下面将采 用微元体法,导出燃气运动的基本方程。
在侧面燃烧装药的燃气通
道中,取长度为dx的微元体,
它由通道左右两截面A(x,t)和
A(x+dx,t)和装药的微元燃烧 A
表面 Ab / x dx 所围成。在
根据能量守恒:在dt时间内,微元体中能量的变化量=燃气 流动迁移在微元体内产生的能量增量+推进剂燃烧加入的能量, 即
t
A(
Ec
V2 2
)
dxdt
AV
(Ec
V2 2
)dt
(
x
dx)( A
A x
dx)(V
V2
)
dx x
T T dx x
V V dx x
(1)连续方程/质量方程
根据质量守恒定律:在dt时间内,微元体中质量的变化量= 燃气流动迁移在微元体内产生的质量增量+装药燃烧进入微元体 的质量,即:
( Adx)dt AVdt ( dx)(A A dx)
I p : 单位质量推进剂所含的总能量(总焓)
Ic Ec P / P / 压力势能
(4)几何燃烧方程
在dt时间内,微元体体积的变化量=装药燃烧掉的体积量,
即:
A dtdx rdt Ab dx
x
x
整理后有:
A r Ab x x
(5)状态方程
根据燃气为完全气体的假设,燃气的状态方程为
燃气在燃烧室中的流动十分复杂,其流动界面是装药的燃 烧表面和燃烧室的内表面。在装药燃烧表面上不断生成燃气加 入燃气通道,即装药的燃烧表面是燃气的生成源。同时装药燃 烧过程中,装药燃烧表面以一定的燃烧速度沿其内法线方向不 断推移,使装药通道截面积不断扩大。由此可见,燃烧室中燃 气流动过程与装药的燃烧过程密切相关,并相互影响。
p RT
以上5个方程为有质量加入的一维非定常流动的基本方程 组,它可用来描述SRM燃烧室中装药通道内的燃气流动过程, 用数值法可求解五个未知参数 p, ,T,V , A 的变化规律。
3. 准定常流动假设
在一定条件下,可以忽略气体非定常流动所引起的不均匀
性,将上述5个方程简化为定常流动方程组,从而得到方程的
( AV )dxdt AV 2dt ( dx)
t
x
( A A dx)(V V dx)2 dt
x
x
( pA)dxdt p A dxdt
x
x
整理后有:
( AV ) ( pA AV 2 ) p A
t
x
x
(3) 能量方程
解析解。 (1) 定常假设 严格上,燃气流动通道在随时间变化,属于非定常流动,但
只要满足 p / t 0,V / t 0,A/ t 0, 即可处理为定常流动。
若A / t 0, 欲使p / t 0,V / t 0,则 V a
燃气的流速远小于当地音速
固体火箭发动机燃烧室及 内弹道计算
燃气在燃烧室中的流动
一、燃烧室中燃气流动的特点及控制方程 二、燃烧室中燃气流动参数的计算
一、燃烧室中燃气流动的特点及控制方程
在固体火箭发动机中,推进剂燃烧生成的燃气通过燃烧室 中的装药通道从喷管高速排出。由于燃烧室和喷管各有不同的 结构特点,故燃气在燃烧室和喷管中的流动也具有不同的特点。 燃烧室中有固体推进剂装药、点火装置和挡药板等零件。固体 推进剂装药在燃烧室中燃烧成为具有一定压力和温度的燃气。 该燃气以一定的流速通过燃烧室内的燃气通道进入喷管。本章 讨论燃气在燃烧室中流动的基本规律。
p
两个截面上,燃气压强、密度、 T
温度、流速分别为 p, ,T,V 和 V
p p / x dx, / x dx,
T T / x dx,V V / x dx,
P和V的正方向与x正向一致。
Ab dx x
dx
A A dx x
p p dx x
若p / t 0, V / t 0, 欲使A / t 0,则
p
推进剂密度远大于燃气密度
2 tb
流动参数发生显著变化的时间要远 小于装药燃烧时间