复合固体推进剂细观结构建模及脱黏过程数值模拟

为研究复合固体推进剂损伤演化规律,基于分子动力学颗粒填充算法构建了HTPB(hydroxyl terminated polybutadiene)推进剂细观结构模型,通过在AP(ammonium perchlorate)颗粒/HTPB基体界面处引入黏接接触替代传统的黏接单元,并基于Hooke Jeeves的参数优化算法反演得到颗粒/基体界面处内聚力模型参数,利用双线性和自定义指数型损伤内聚力模型模拟了AP颗粒和HTPB基体黏接界面处损伤的萌生、发展、聚合直至宏观裂纹破坏的过程。通过数值仿真与实验结果对比发现,指数型损伤内聚力模型比双线性模型能更准确描述推进剂单轴拉伸过程中颗粒与HTPB基体界面间脱黏过程。最后对比了多阶段加载实验结果与仿真结果曲线,发现两者变化趋势基本一致,最大偏差仅为10%,验证了所建细观模型的可靠性及反演所得界面参数的准确性。      

Experimental research on rotating detonation with liquid hypergolic propellants

This paper describes experimental research into the initiation and propagation of rotating detonation for liquid Nitrogen TetrOxide (NTO) and liquid MonoMethylHydrazine (MMH). An annular rocket-type combustor without nozzle was designed to investigate detonation combustion. The propellants were injected through unlike impingement injectors. The combustion flame fronts and pressure waves were detected using optical diagnostics and dynamic pressure sensors, respectively. The propagation of rotating detonation was established spontaneously by increasing the mass flow rate of propellants. The velocity of propagation of the flame fronts and pressure waves was nearly equal and reaches supersonic speed. Two different detonation combustion patterns are present, single wave mode and double waves mode. And in double waves mode, the two detonation waves are always counter-rotating. The possibility of rotating detonation initiation in a combustor with nozzle was also checked. Stable rotating detonation can be initialized and sustained at similar operating conditions.     

HTPB固体推进剂老化性能检测新方法

引入了光致正电子湮灭分析方法 （PIPA,Photon Induced Positron Analysis）,对HTPB固体推进剂老化性能进行检测研究。介绍了PIPA的原理、PIPA试验平台的搭建以及PIPA的数值处理方法,并用511KeV能谱的FWHM值ε（半高宽,Full Width Half Maximum）表征HTPB固体推进剂老化引起的微观变化,所得结论与固体推进剂的常规测试一致,证明了PIPA用于固体推进剂老化性能无损检测的可行性。     

幂律流体液膜破裂的线性稳定性分析

为对凝胶推进剂撞击式喷嘴产生的液膜的破裂行为进行研究,使用线性稳定性分析方法并结合幂定律本构方程对幂律流体液膜的稳定性进行了研究,并对流体物性参数、流变参数等对液膜流动稳定性的影响规律进行了分析.结果发现:稠度系数、流动指数和表面张力越大,液膜表面扰动波的增长率越小,波长越大.而较高的液膜速度、气液密度比和液膜厚度都使扰动波增长率变大,但波长随液膜速度和气液密度比的增大而减小,而液膜厚度的情况却相反.理论计算与实验值的对比表明,线性稳定性分析能较准确地预测液膜的破碎长度,但预测的表面波长却大于实验值.      

基于细观颗粒夹杂模型的复合固体推进剂松弛模量预测

为更准确地预测不同固体颗粒体积分数的复合固体推进剂的松弛模量,采用了分子动力学方法对不同体积分数的复合固体推进剂细观模型进行建模.根据有限元理论及细观力学均匀化方法,计算在定应变工况下复合固体推进剂细观模型的平均应力随时间的变化,从而有效地预测复合固体推进剂的松弛模量.该方法有效地体现了随填充颗粒体积分数的增大,复合固体推进剂瞬时模量逐渐增大的变化规律及颗粒随机分布对复合固体推进剂瞬时模量的影响.将其应用到复合固体推进剂的设计过程中,可有效降低设计成本,缩短设计周期.      

含损伤粘弹性本构及其在有限元分析中的实现

为研究固体发动机药柱的损伤粘弹性行为,针对一种耦合损伤的三维非线性粘弹性本构模型,采用积分算法对其进行了数值离散,导出了本构方程的增量有限元形式.针对该增量型本构方程对应的非线性方程组,建立了采用Newton-Raphson迭代法求解的应力更新方法,并基于Abaqus软件和Fortran语言开发了相应材料本构子程序,最后应用算例对所设计数值方法及程序进行了验证与分析.结果表明,本文有限元解与对应解析解吻合良好,保持了很高的精度,所设计方法及程序有效,可应用于实际型号发动机药柱损伤结构分析.     

复合固体推进剂颗粒的三维网格自动生成

为了分析圆柱颗粒对复合固体推进剂性能的影响,必须对复合固体推进剂中三维圆柱颗粒进行网格划分。根据"力平衡"网格划分原理设计了一种适合于三维圆柱颗粒的网格自动剖分方法。编程运行结果表明,方法能自动完成划分,网格质量较高,满足要求。     

近临界区液氮低速喷射特性研究

针对低温推进剂在超临界环境中的喷射特性,以液氮为模拟介质,基于SRK状态方程和LES模拟方法开展数值计算研究,获得了4MPa压力下液氮跨临界射流的形态特征,射流密度分布规律与试验结果吻合良好。计算结果表明：在临界点附近,液氮射流表面会形成高比热容屏障,抑制射流内部流体温度升高,从而维持射流核心区稳定;核心区射流表面涡的形成和发展由斜压效应主导,随着射流向下游发展,斜压效应、体积膨胀和粘性效应三者对涡量输运的贡献趋于同一水平;液氮射流破碎后形成并维持大的“高密度块”形态,随着温度升高,密度块逐渐扩散消失。     

助推器喷口火焰温度测量

火焰温度测量是一个复杂而困难的问题。接触式测温传感器在很多情况下并不适用。固体燃料助推器火焰中有大量固体粒子,这些固体粒子运动速度又很高,采用接触式测温传感器测量这种火焰温度问题较多。介绍用四支小型不调焦光电高温计从不同角度测量固体燃料助推器喷口火焰温度的试验装置和试验结果,并对小型不调焦光电高温计的基本原理、结构、性能作了简单介绍。     

干扰因素对补燃式液体火箭发动机性能的影响

给出了描述补燃循环液体火箭发动机静态特性的非线性数学模型,用拟牛顿法求解各性能参数,分析了干扰因素对液体火箭发动机参数的影响,所得结论可用于发动机试验结果分析、发动机可靠性分析、发动机故障分析,也可用于揭示发动机参数随各种干扰因素的变化规律.