QJ 1489-1988 固体火箭发动机内弹道计算方法和计算机程序.pdf
《QJ 1489-1988 固体火箭发动机内弹道计算方法和计算机程序.pdf》由会员分享,可在线阅读,更多相关《QJ 1489-1988 固体火箭发动机内弹道计算方法和计算机程序.pdf(26页珍藏版)》请在麦多课文档分享上搜索。
1、。J中华人民共和国航天工业部部标准QJ 1489-88 固体火箭发动机内弹道计算方法和计算机程序1988-04-25发布1988-12-01实施中华人民共和国院去工业部发布El支1 主题内容与适用范理.h警们(1) 2 引用标准.u.( 1 ) 3 零维内弹道计算.a. ( 1 ) 4 计算旅程序山.川. 川.,. . .川. 川. . . 川. 川. .川.川. . . . . 黝录A二术仨i远击(补充件).川. . .川. 山. ,叠. .川. ,. . 川.川.川. .,山鑫川.川川. 川.川.川. 川川川.( 9 ) 董宵E录Bt刮计i-3窍?宇-饥程序SMIBC.FOR(补充件 .
2、oH. ( 10) 黝录Ctj去ti(参考件.HO. . .川.( 17 ) 中华人民共和昌航天工业部部标准程J1489-88理体火箭发动辄内弹道计算芳:圭和计算扭程序1 主噩内容与适用范本标准规定了固体火箭发动机内弹道零维tr算方法和计算机程序,本标准适用于固体火箭发动机的内弹道计算a2 引黯标准QJ 1041 段体发动机应强一时i湾、推力-1l1间数据主主理娩范QJ 912 复合国体推选弗j药条燃速的水下声发射测定方法QJ 915 复合国体推选剂药条燃速、然速压力兹数和然速温度敏惑系数的测定方法QJ 13吉3固体火箭发动机热力学计算方法剥计算机程序事维内部遭计算3.1 已知.戴辛DliI童
3、在关系3.1.1 已知参戴已知参数J且表1 竞主1 符号i计算程序捋号!名称单位铭Hi强?室主生Y. 革.m燃烧室燃气灼比唤出 念au要靠速温度敏惑系数9Io/k 40 602 To R古约燃速系辈主.01 cm/ .1.2 已知画.羡慕a. 肉厚Wj一燃商4白的对应关系把药柱沟厚分成j(吉、1,2.bf)份,并绘出与路气相对王室的燃商.1&Job. 根据模拟发动视或条件辍近、喉持材料相同皇宫有关试验绘苦哥哥主部燃烧量!J.tj-tj始采系曲线如愿1所示.其中!J.dtk.也是已知多数,既可以庭主量也可以外插.A d.t I _,/,/,/ (C时I_-I I IAdtK G 乞且tl哥t(S
4、) 因16dJ-il关系2 QJ1489-88 C. 温Jf下的然速关系由QJ9t5或QJ912提供最差连关系s, r=向012 式卒,a,约单位,按麦2换雾。主$。J95 QJ 912 给出岛手口a , 表2 本标准交用.的单位a , 4 m.羽4.085-,-, (MPa)国位1mC口2o 4 2 剖4一一一-一一一一-一一( kgf/cm .) . s 。1614( kgf!豆m.) S 3.2 注lMP.=l/乱。98.0665KgfJc扭a计篝公式计算药控初远为Tgj温度下的内罢手道计算公式如下z是去速系数 ( c血/( (MPa n) S) ) a. b. C. 矗,=aoe(T,
5、gi-Ta主瞬时压强P., (MPa) ,且L._ 1_0_0 pcj= C a.Pgr.可巳*c.o_二二一一-Ai (10) j=吉、1、2、.bf4, pcj= l a -一-一且一一一号、自在由.665招邻两点平均应缉F叶(主Pa)-4EL3 /j垃p,严EFc汁pcti_l) j=l, 2、bf暧管喉径串串时t重重!J.dti ( cm ) ( 1 ) 1 - n /10. ( 2 ) l-n .川(2a). ( 3 fldtj=Kt fj_1 u. ( 4 ) j=2、3、bI式中z,id KL- n,aLk . n. . (4a) t. 3 QJ1489-88 领管瞬时喉径dj(
6、cm) dt_-i =dto +.dtj j=2、3、bf.(5 ) 仁堪言管探对眼部部积4 (cm) Atj=士3fdlj=2. 3 .也f.( 6 ) 也相邻两点时间间隔Mr(s) Wi-W, .1t=ti- tj_l =一.J , _l-! j = 1、2、.b . ., ( 7 ) h. 然烧时间Ij (s) bf Ij= E M , j=。. 辛苦始喉逐比J , ap叫四j=o、1、2,bf. . ( 8 ) A F-L n-HUM-. . (川4p, 3.3 计算步骤a. j=O对,把,hhAt代入公式主得p吨,又因fj_l=0,所以fe=f;JI/J=号,从觅得Peo。b. 当j
7、=lR苦,把4b、4t.代入公式2得Pel,再出公式3得pclo然后由公式7和E得t1 , 由公式4革和B得.dtl、d程,:1t 1 Q C. 当j=2对,令d.,=dh,4,=4口,由公式Z幸在3得Pe和PCJ.。再出公式7和B得1l丽可解应当j二2对的一系列参数ed.当j=bf对,Pcj=Pbh fj=tbo e. 当j捉对,j=bf于1、才2,af-l,af时,pj按下式近似计算pcj=-一1-AEbL一一一一_.(10 ) 二号之(tj-tb) e 0: Ve+l 式中sr=v7-二-.)万石,-1). ( 11 ) Vc -r 1 求婚时首先在d.;-tj曲线上太子t.之后找一些小
8、步长仇。20)约1;得.,d:h11; , 然后代入公式10得.Pcjo当j=afR苦,Paf=O.3号MPa.如果P叶=pcf a;_l O.30MPa, 而Pcj=Pe(lln O f .,. (1吟10 . F - -. 卢丁可CFCF Pcj4ti., .主,.,. 09. . ( 13 ) 10 ri III. 沪、!cj和/tj分裂采用s.m!s、kgf!cm和cm2时,则s。80.66只.1 Fj=rl飞一一一7.:!_ PC Ar; .H. ( 13a) l c-式中z-7 当采用实际比冲1.和计算推力系数c,对,推力系数园子可以等于lOo当由QJl393提供的理论比1中f川和
9、理论特征速度由计算得的CF由为理论握主力系主立时,罢d公式12、13中击口自pdyo 自2IDPC口u自m口kuu aadyc口128) sbp口DJMMVA吨uvgdmv刷M恻rrnHUeaeo mott1) aJGCS口raanno 91、rreoaa-qmt rhhL orccdg -4 l23456 7. 8 9 15 oen(6,+Ie细胞,tn ) ,;r 1奄.(巾,15) forma专(8x. 驾.李榈2、生舰想嚣移响. 10 = teoad (唱(a) c t1 11r te (6, 155) 12: rie ,/n)?CnJ) 171 rllad (咿,(a)oV 18 I
10、 f (1)P. e. y) then 19. oen(7, i le- datal , status蝠,0 l d ) 20 = read (7. (9f1口.5,5千10,1.7110.3)pn. muc,受c.a02, CT. 2:0嘀变k,211 1 d奄0,st日,cl.c1 1. c12,bf. P口,企口,.币2,9 C, 9f, 22 2 m.l.奄.2231 etse 47 9 c oo l A宁ac p 口变ew 日专do p ka t 4 0b 进HVWdmlc du、,Ll Ol-c 翩.,8日lulc 贷币,a$2 食98冒也.D-419 2。币进牛、7ac dA崎U
11、m 嗣nuelo ;、4Z向旧鸥Jwdy,) i: hen olf i te (6. (4x.511口.5)n , MuC. l:c,a02,。千ilse w r ii:e (巴(4)E fpefl !loto 60 oPn (4, file-da专a2 s仓a:-!:us-old) 1 ead (4. (8口矛8.2) 。伶en(白,f i L e帽da-l:a3,status old) 俨ead(白, (8口币1口.2)ab :I,) l eo . d 飞衔(3) ) 1,1 lrite (6, 48) 手orma专(9x,同.8:.s=O.3 .0百.02.千黯.手2abp:.口.。do
12、 62 k1蜀.:;2,2巾1J,I+1.j,. ()J/2口JWD W口国/口OP-0篇b口细、Jya篇lml11) ) 126: caLL ea%(l), iSP , plis) 1271 call cz(l),c树,oc128币(1)辑9f*is糊980.665*ai:(1)+C1)/cx1291 endf 13日Ido 1口okl.;2.2湖11,V+1131. jf (t(kl-1l而.t口)8日,8日,81132 80 dt(kl)-dt口1331 30 啦。82134. 81 d,(kl)-dt日+变(kl-1)+tk 135 82 at (kl)捕。25+3.14159:26事
13、dt(k!)华伽21361 i f (VP. Ilq. y) then 1371 (kl)(9c+al*zo事cl事.b(kl)/.,(kll响(1日.0)嘟嘟(2G咿PO)巾*1:391 1 (1.口/(1.口-po)/1口口00日日139 f(kll-.币咿C手巾P(kD .进(k1)/1口,口1401 else 141 .(kl)皿 CX) 1441 1,) then 162 1 eps踊0.31631 else 164: ps-3. 06 1651 endi千1661 I千(2中1,0)+1).L. S) -i: heoli 1671 i:b-i(2伽!+1)/仓(2*i,+1)169
14、, f奄b=i手t(2丰1#+1)/1:(2华1#+1)169i 90吧。154170: I?ndi千12 QJ1489-88 171 do 102 kll, 15 172= 1口2Wf i 1: 02 (6, 155) 1731 write (丰,110) 174. 110手ormai:(4x.、-.入Idtl口) 175, r e ad (. (f 1 O. 5) ) d t l 0 176= 9am-5rt(muc)潮(2.0/(muc+1. 0)咿嗡(muc+1.0) /飞2.。中(muc-1.0) 177: kl=2*jw+l 1781 d奄L=dHO179= tl (kl)=O.
15、0 18日.115 klmkl+l 191: i千(pe手302.p (kl-l). and.庐e千.L e. P (kl-2) then 182: te手=(t(kl-1)-t(kl-2)/(p(kl-l)-p(kl-2)巾(pef-p(kl-2)+毛(kl-2)1831 dtef-d专(kl-2)+咆手巾奄k1841 ate币=0.25中3.1415926丰d变e手中.2185: if (P. eq. y) then 196 fef=9千巾C手划.e手巾a奄o2f/1口口187= else 188: call caz(pe手.i sPdl i sp) 189: caLL caz(pe手.C
16、X. pcx) 190: fef-9学申isp*98日665+pe+*奄e-f/cx191= endi手192= ppef-(p(kl-2)+pe千)/2 193. fpe手-(f(kl-2)+才fef)/2 1941 IPe币冒jp(kl-2)+ppef*dtl195: i手ef-if奄(kl-2)+fpe千响dtl196: i手(pP.eq. y) then 1971 3I l Lef-sll (kl-2)+ppeftoate手中dtl巾1日00口/01198. lse 199. call caz(pe币,cx , pcx) 2口口:9 l l e手9lt(kl-2)+ppe千*atET巾
17、d变l*980.665/cx2口1:end i于2口2:endi学203. 125奄1(kl)=tl(kl-1)+dtl204. t(kl)=奄(2+i、以+1)+t1(kl) 2口5d奄(kl)-dtO+专(kl)丰奄k2日6.仓(kl)0.25巾3.1415926*d奄(k1) +中22071 p(kl)-p(2刷iw+1)/(exp(go3m申期2+c1+a奄(kl)中t1(kl)/b币、208 i手(PVP.eq. ,)奄hen209 千(kl)=9手0手巾p(k1) +03奄(k1) /1日-口21日.else 211= c o3 ll c o3 z(p(kl). isp, pl i
18、sp) 212 caLl caz(p(kU , cx , pcx) 213. 手(k1)=g手巾5P丰98口.665p(k1)*at(kl)/cx214 endi手215= pp飞kl)=(p(kl-1)+p(k1)/2216 千P飞k1) (学(k1-1、+千飞kU)/2 217: i p ckl丁=ip(k1-1)+pp(kl)事dtl218. i千变(kl)=if仓(k1-1)+fp(k1)+d专L219= i千(pP.e o. y) then 220: gLl (kl)=gll (k1-1)+P(kl)巾at(kl) +dH响10000/c1221 eLse 2221 call caz
19、 (P (k1) , cx , pcx) 223 gL l (kl)=3 l L (k1-1) +pp (kl) +at (kl) *d奄L巾98日665/cx224. endi币225 = w俨i变e(+, (2千10.2)奄(kl),p(kl) 226-i手,.飞kl)-eps)120,12日,115227 120 d变la-dt:l牟eps/p(k1-1) 229;专1(k 1)司专L飞kl)-d奄l.229 奄飞kl).奄1(k1)+t (2* iw+1) 23日,d专1、dtO+奄(k1)巾tk13 QJ1469-88 ,E -需-帚帚踵pp b)b八阿臼,AM、川、J,川呻呻叫M川
20、M叫川川x内kkku队2dde1461ec/C1品晰、t111棚t/m/2k3)2辙mcm曲111,、,boC/C*日takk9X)K3/4T/)at牛、叫他aauamdf94fLC品AMA吃牛飞+、,PPD包k9kl)LdaaE问问咧-XURKK/14飞庸曲曲,、954)走k22dylthvc1lpcdl)/串llLl+b2b3*kpp)LkkkiMW/a34?3356钱341136$曲1Vhb-e-a2细b-e?埠。2nt$11ukpnpP曲llenaL变?呼叫/啦?bg日9471t74?唯79eatv6kk41ehvxp咿2Kleeff手4TI?k4JdTea4477学45h山pc6P+
21、h+C牛?Cht5罩tB3-t)P-PtTPtmw七P/舍tTb4l毡1k4b于a41p咱血+elllt/)112)l牛Laeatmb嗣lzl变ee4Tt鹏lKSV内白、1)1FJ-x一/、Jlk/-atpbp4pak-pkttpbhvehva.,C91111-eL)-KCtetbj阳tpt3yhv由mm皿wkkklA守千tl,X)YTm)P)PPK严TmplP审华SSP111,(Spt+川wor4Tt/+口舍5545巴9口1$2311c2345JTppeekKKt-PLUwttoe/WPO55(55-5FbLb669白gMAgs2qceel千qLLt*/电lfail-qlxm-xlx-x-x
22、lx-xmlv411447lxf筝申中P4TP-e。事白,$24T/1bRe/sezA峙C14A峰A峙A饰A啡4c49490hv4T22474飞tgEZZe2efk2+巳566Lb邑巴6FW6(6(6(k6k6LK辅aep俗aa9帽且鹏P忽thve食1C内V变(JK)剿幽CC事)LVLclpdIR?PJTPttt变+、进ttt舍监鼠t1)P4V-1-bkpkk4?+(aeeaeepea-ee-Gaeaeaeaealeasem。e息klleit-bRKKeL姐姐3翩翩事ma幽mttt变tm-tm全mtmtmtmetm-tmztr2m2LKUKSLLKdttgLLdbb千dTaaE而bkrfror
- 1.请仔细阅读文档,确保文档完整性,对于不预览、不比对内容而直接下载带来的问题本站不予受理。
- 2.下载的文档,不会出现我们的网址水印。
- 3、该文档所得收入(下载+内容+预览)归上传者、原创作者;如果您是本文档原作者,请点此认领!既往收益都归您。
下载文档到电脑,查找使用更方便
5000 积分 0人已下载
下载 | 加入VIP,交流精品资源 |
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- QJ 1489 1988 固体 火箭发动机 弹道 计算方法 计算机 程序
