下载文档原格式
固体火箭发动机内弹道学方丁酉固体火箭发动机内弹道学方丁酉近年来，随着航天事业的蓬勃发展，固体火箭发动机内弹道学成为越来越受关注的领域。
方丁酉，中国工程院院士、火箭技术专家，曾在这一领域做出了卓越的贡献。
本文将从简单到复杂、由浅入深地探讨固体火箭发动机内弹道学，以帮助读者全面、深入地理解这一话题。
一、固体火箭发动机内弹道学的概念固体火箭发动机内弹道学是研究固体火箭发动机内燃烧过程、燃料燃烧特性、燃烧产物排放、工作环境等方面的学科。
方丁酉指出，固体火箭发动机内弹道学是固体火箭技术中的重要基础理论，对于提高固体火箭发动机的性能、可靠性和运载能力具有重要意义。
二、固体火箭发动机内弹道学的理论基础在固体火箭发动机内弹道学的研究中，燃烧动力学是重要的理论基础之一。
燃烧动力学研究了燃料在燃烧过程中的变化规律，以及燃烧反应对火箭发动机内部气体流动和压力变化的影响。
燃烧产物排放和燃烧室内部气体动力学也是固体火箭发动机内弹道学的重要内容之一。
研究这些理论基础可以帮助我们更好地理解固体火箭发动机内部的工作原理和特性。
三、固体火箭发动机内弹道学的关键技术在固体火箭发动机内弹道学研究中，燃烧稳定性和效率是两个关键技术。
燃烧稳定性是指在固体火箭发动机工作过程中保持燃烧的稳定性，避免出现燃烧不均匀或燃烧失稳等问题。
而燃烧效率则是指在燃烧过程中尽可能提高燃料的利用率，减少燃料的浪费。
方丁酉在固体火箭发动机内弹道学的研究中，提出了一系列有效的技术方案，使固体火箭发动机在燃烧稳定性和效率方面取得了显著的进展。
四、固体火箭发动机内弹道学的应用前景固体火箭发动机内弹道学的研究成果已经在我国的航天事业中得到了广泛的应用。
在长征系列火箭、嫦娥探月工程等多个航天工程中，固体火箭发动机内弹道学的研究成果为提高火箭的性能和可靠性做出了重要贡献。
未来，随着我国航天事业的不断发展，固体火箭发动机内弹道学将继续发挥重要作用，推动我国航天事业迈向新的高度。
固体火箭发动机工作末期的内流场数值计算张俊; 高天宇; 高璞清; 王刚【期刊名称】《《弹箭与制导学报》》【年(卷),期】2018(038)003【总页数】5页(P58-62)【关键词】固体火箭发动机; 飞行过载; 颗粒轨道模型; 两相流; 燃气射流【作者】张俊; 高天宇; 高璞清; 王刚【作者单位】上海海洋大学工程学院上海201306; 中国航天科工集团第六研究院第41研究所呼和浩特010010【正文语种】中文【中图分类】V2310 引言固体火箭发动机的工作过程按照其燃烧室压强变化可大致分为点火建压期、稳定工作期和工作末期。
在发动机点火使其燃烧室达到一定工作压强后,进入持续时间相对较长的稳定工作期。
随着推进剂燃面的推移、肉厚不断减少,衬层和绝热层在高速高温燃气射流的烧蚀与冲刷作用下可能出现剥离、脱落、烧穿等现象。
尤其对于含金属的复合推进剂而言,产生的高速高温颗粒受到发动机的结构限制和飞行过载等条件的作用,将会使发动机后封头及喷管收敛段等部位的粒子沉积浓度升高明显,进一步加剧了这些部位的烧蚀。
因此,研究SRM的工作末期的多相内流场变化,确定流场中固相颗粒的分布规律、运动轨迹及沉积浓度等是一个必不可少的环节[1-3]。
目前,科研人员主要通过数值模拟研究这一问题,取得了一系列重要的研究成果[4-6]。
文献[7]对发动机的燃烧室进行冷流实验,对径向速度、脉动速度等物理量进行了实验测量;文献[8]将SST k-ω湍流模型及其改进形式用于燃烧室湍流流场的数值模拟,并将计算结果与Wilcox和S-A湍流模型进行了对比,得出SST k-ω计算的燃烧室径向速度分布与实验值一致性最好的结论;文献[9]对凝相粒子加速度进行了理论分析、数值计算及缩比发动机旋转模拟试验,分析了地面模拟过载试验中燃烧室粒子聚集区及旋转试验后的烧蚀部位;文献[10-11]对高含铝固体推进剂低压SRM尾流场复燃进行了数值模拟和实验研究,得到了铝粉引入、粒径和燃烧室压强对羽流温度的影响规律,并与地面试车结果进行了对比。
基于Ｐｒｏ／Ｅ软件进行固体火箭发动机内弹道计算的方法初探作者：王磊代义李宏来源：《中国新技术新产品》2009年第10期摘要:固体火箭发动机内弹道计算给出的推力、质量流量随时间变化的曲线一直是导弹或卫星外弹道设计的计算依据。
为适应需求,固体火箭发动机装药的几何形状需要反复调整,复杂的几何形状使内弹道计算成为一项繁琐的工作。
本文利用Pro/E软件强大的三维建模能力,通过参数化设计模拟发动机燃烧的全过程,编制简单程序实现发动机零维内弹道计算,为繁琐的计算提供了一条便捷的途径。
关键词:固体火箭发动机;内弹道计算;Pro/E软件1 引言固体火箭发动机是一种采用固体推进剂的化学火箭动力装置,在航天技术中有着广泛的应用。
内弹道计算是固体火箭发动机设计的一项重要内容,它提供的推力、流量等性能参数是导弹总体进行射程计算、散布分析及导弹姿态控制的重要依据[1]。
在满足内弹道特性和保证药柱结构完整性的前提下,设计出具有较高装填分数的药柱对于提高发动机性能具有重要的意义。
药柱几何形状的演变过程,是从简单管状药柱演变到贴壁浇注内孔(星形、车轮形等)燃烧药柱,再到具有较高装填特性的翼柱、锥柱等三维药柱。
内弹道计算的关键在于得出药柱燃面随时间的变化曲线。
有着简单几何形状的药柱,计算燃面都有较为成熟的方法。
为了探索三维药柱燃面的计算方法,国内外有关学者开展了许多工作。
目前,国外内弹道计算最为著名的方法是美国NASA的SPP法,国内较常用的是通用坐标法、有限元素法和边界坐标法。
国内的这三种方法均为数值法,由于使用原函数对复杂药柱进行精确定义很困难,通常难以精确表示燃烧过程中十分复杂的图形,以及连接各部分药柱圆角的复杂变化,所以内弹道计算的精度较差。
同时,它们的计算结果对药柱网格划分的依赖性很强,网格划分不合适,所得结果也不正确[2]。
这些方法计算得到的燃面波动性也较大,药柱初始设计的修改很复杂,数据前后处理也抽象复杂。
这些因素导致内弹道的计算成为一项繁琐复杂的工作。
固体火箭发动机瞬变参数辨识
方丁酉
【期刊名称】《航空动力学报》
【年(卷),期】1995(10)1
【摘要】叙述了固体火箭发动机参数辨识的方法，以反喷管打开后瞬变过程中推
进剂瞬变燃速、打开函数和喷喉收缩因子的辨识为例，说明了瞬变参数辨识的技术。
辨识采用的数学模型是零维不定常内弹道计算方程组，这是一个非线性系统的参数辨识问题。
优化方法采用阻尼最小二乘法。
给出了算例，经参数辨识准度分析，辨识结果是合理的。
【总页数】4页(P67-70)
【关键词】固体推进剂;火箭发动机;参数辨识
【作者】方丁酉
【作者单位】国防科技大学
【正文语种】中文
【中图分类】V448.153;V435.11
【相关文献】
1.基于遗传算法的固体火箭发动机参数辨识 [J], 樊超;张为华
2.固体火箭发动机喷管喉径瞬变值的计算和预示 [J], 余贞勇
3.固体火箭发动机参数辨识软件关键技术研究 [J], 魏瀚;鲍福廷;蔡强;刘旸
4.固体火箭发动机反喷管打开后的瞬变参数辩识研究 [J], 杨涛;张为华;代绪恒;方
丁酉
5.固体火箭发动机稳态燃速二维模型参数最优化辨识 [J], 李晓斌;张为华;王中伟因版权原因，仅展示原文概要，查看原文内容请购买。
固体火箭发动机装药临界弹性模量的计算文章题目：固体火箭发动机装药临界弹性模量的计算摘要：本文详细描述了固体火箭发动机装药临界弹性模量的计算方法。
首先，本文介绍了国内外固体火箭发动机技术发展现状，并简要分析了固体火箭发动机装药结构影响安全性和可靠性的关键因素。
其次，本文依据国际上公认的标准测试方法，以及无元素数值分析、有限单元分析软件模拟计算，测试计算固体火箭发动机装药临界弹性模量。
最后，本文利用TEST软件对所测试计算的结果进行验证，并给出对前人研究的改进建议。
关键词：固体火箭发动机；装药；临界弹性模量；数值分析；有限单元分析固体火箭发动机装药临界弹性模量的计算方法对固体火箭发动机安全性和可靠性有着重要的意义。
固体火箭发动机装药临界弹性模量的计算方法可以帮助研究者们准确地确定火箭发动机装药的各项质量指标，从而更好地控制固体火箭发动机装药的性能。
例如，通过使用精确的固体火箭发动机装药临界弹性模量，可以准确地计算推力、发动机总极限压力和器件的最大使用量。
此外，临界模量的计算方法还可以准确地预测包括热推力、阻力和持续时间等在内的复杂参数，以便在设计固体火箭发动机装药时尽量节省成本、提升性能和安全性。
由于固体火箭发动机装药临界弹性模量的计算方法是一个复杂的过程，因此研究者们可以采用数值分析、有限单元分析以及无元素分析等技术来提高计算的准确性。
例如，在数值分析中，研究者们可以精确计算弹性模量，同时也可以控制变形量，从而优化固体火箭发动机装药的可靠性。
在有限单元分析中，研究者们可以定义节点的质量、位移、温度和压力等参数，而在无元素分析中，研究者们可以控制风扇静压压力或流速以及装药以及固体火箭发动机腔体变形量等参数。
固体火箭发动机装药临界弹性模量的计算方法对确保固体火箭发动机的安全性和可靠性起着重要作用。
此外，其已经在多种国内外的火箭发动机上得到广泛的应用，也为固体火箭发动机的技术发展提供了重要的参考依据。
因此，借助固体火箭发动机装药临界弹性模量的计算方法，研究者们可以更好地控制火箭发动机的性能，从而提高安全性和可靠性。
弹道仿真中的数值计算方法研究
弹道仿真中的数值计算方法研究是一门比较复杂的学科，包括了实际的物理模型、数值计算方法和实验测试等。
它主要用于研究弹体在飞行过程中受到的力和偏移，以及其飞行轨迹。
在弹道仿真中，重力力、空气阻力以及其他相关因素都会影响弹体的运动，而这些因素也都必须通过数值计算来研究和模拟。
因此，数值计算方法是弹道仿真中不可或缺的部分，对于弹道仿真的精度有着重要的影响。
在弹道仿真中，常用的数值计算方法主要有两种，分别是有限差分法和有限元法。
有限差分法是将物理模型的微分方程表示为一系列的差分方程，并使用时间步长和空间步长来解决。
有限元法是将物理模型的微分方程表示为一系列的有限元方程，并使用有限元单元来解决。
在实际的弹道仿真中，数值计算方法的选择是非常重要的。
根据物理模型的不同，可以选择有限差分法或有限元法，以便尽可能准确地模拟弹体的运动。
此外，在弹道仿真中，还需要考虑时间步长和空间步长的选择，以及如何处理弹体碰撞、曲面反射等情况。
此外，在研究弹道仿真数值计算方法时，还需要考虑数值计算方法的可靠性和准确性。
例如，可以使用数值实
验或者实际试验来检验数值计算方法的准确性，也可以使用不同的数值计算方法来对比模拟结果。
总之，弹道仿真数值计算方法研究是一项复杂的学科，涉及到了物理模型、数值计算方法、实验测试以及可靠性和准确性等多方面的知识。
在实际的弹道仿真中，要选择合适的数值计算方法，并加以优化，以便尽可能准确地模拟弹体的运动。
固体火箭发动机的参数及其偏差预测摘要：本文弗研究固体火箭发动机的参数和其偏差预测技术。
在多层次实验结果的基础上，采用有限元方法建立参数及其偏差预测模型。
模型考虑了固体发动机的机械特性、材料性能、燃料/推力形成方式以及发动机的结构参数等因素，为火箭发动机的参数预测提供了更加准确的基础数据。
关键词：固体火箭发动机；参数偏差预测；有限元方法；发动机结构。
正文：1. 研究背景：火箭发动机是运载火箭飞行的关键设备之一，它的参数偏差对于其机械性能、热容量及稳定性有着重要影响。
因此，对固体火箭发动机参数偏差的准确预测有着重要意义。
2. 研究方法：本文采用有限元方法建立参数及其偏差预测模型，利用多层次实验获得的固体发动机的参数及其偏差输入，研究参数预测模型准确性。
3. 研究结果：研究结果表明，采用有限元方法进行参数及其偏差预测可以根据各种因素输入得到准确的预测结果，并且可以更好地反映火箭发动机的参数及其偏差的变化趋势，为发动机详细设计提供了可靠的数据支持。
4. 研究总结：本文研究表明，有限元方法在参数及其偏差预测方面的应用可以取得满意的精度，且可以更好地描述发动机参数及偏差的变化趋势，提供准确的计算结果，为发动机精细设计和维护提供了更加精确的数据支持。
固体火箭发动机参数偏差预测技术的实际应用一般要经历三个步骤：首先，进行规范的发动机计算，识别发动机的特性，包括材料耐久性、燃料系统、径向推重比等。
其次，使用有限元技术，对发动机参数进行大量仿真实验，例如：燃烧室压力、燃料组成、温度分布、噪声特性等，以收集发动机参数和其偏差的实验数据，并以此作为预测模型构建的基础数据。
最后，基于计算结果和实验数据，构建参数偏差预测模型，并通过模型验证精度，最终获得可靠的参数预测结果，以此来保证发动机的稳定性和效能。
此外，固体火箭发动机参数偏差预测技术也可以应用在多种其它发动机系统中，可以有效提升发动机的性能和稳定性，反映出技术的广泛应用价值。
高速旋转固体火箭发动机内流场数值计算郝雯;封锋;罗盟;谢爱元【摘要】为计算高速旋转固体火箭发动机内流场特性,采用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，给出质量入口边界，对端面和内孔同时燃烧的管状装药固体火箭发动机进行了三维内流场数值计算。
!计?算?星?型í装痢?药?的?几?何?尺?寸?!参?数簓符?号?说μ明÷!-----------------------------------------------------------------------------------!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)-foresmaxif(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!该?子哩?程ì序ò用?于 ?计?算?零?维?变?截?面?燃?烧?装痢?药?的?内ú弹獭?道台?!利?用?有瓺限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,Freal(kind=8)::ep,P_av,F_av,I0,Im,Iv,Cf,gaPe=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!采é用?有瓺限T差?分?法ぁ?计?算?压1强?Poc=Poc+(ga*ga*Ctz*Ctz*dt/(Apa*len))*&(density_p*Sa*ep*a*Poc**pn-Poc*At/Ctz)Cf=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。
