基于损伤的HTPB推进剂／衬层界面内聚法则构建

为提高丁羟（ HTPB）推进剂／衬层界面数值仿真结果的准确性,首先改进了单搭接试件,通过剪切实验获取界面断裂参数,而后分别采用双线型模型和指数型模型对试件进行数值研究,进行影响分析,并在此基础上建立了一种基于损伤变量的自定义内聚力模型,使模型具有了明确的物理含义。结果表明,改进的单搭接试件能够测定柔韧粘接件的II型断裂参数;II型界面断裂时,内聚力曲线形式与实验曲线形式一致;该基于损伤的内聚力模型能够比双线性模型和指数型模型更准确地反映HTPB推进剂／衬层界面的II型断裂性质。     

基于扩展有限元的改性双基推进剂的开裂过程模拟研究

为研究和模拟改性双基推进剂(CMDB)在准静态条件下的断裂特性,利用扩展有限元XFEM技术,建立了中心直裂缝含加载平台的圆盘形三维模型。基于线粘弹性理论,通过单轴拉伸实验直接获取改性双基推进剂的最大拉应力8.21MPa,作为仿真模型中裂纹粘聚区的断裂强度,模拟了固体推进剂在静态加载条件下裂纹实时扩展过程。结果表明实验和仿真过程中裂纹均首先呈I型扩展,最后以复合断裂形式扩展失效。讨论了推进剂在压缩断裂过程中裂纹扩展区的应力变化对裂纹扩展的影响,发现实验预期结果与仿真结果符合性良好,说明扩展有限元法能够为改性双基推进剂断裂过程的数值模拟提供新的方法。     

冲压增程弹进气道型面气动优化研究

以某型冲压增程弹为研究对象,为了使进气道总压恢复系数最大,利用自动优化平台,集成了参数化建模和进气道流场数值模拟,以总压恢复系数最大为优化目标,运用组合优化策略对双锥进气道的锥角进行气动优化,优化后的进气道总压恢复系数提高了11.96%.数值模拟结果与试验数据吻合较好,证明优化结果可信.结果表明,通过对进气道锥角的优化可快速找出最佳锥角组合.     

基于Hopkinson技术的CMDB推进剂动态起裂韧性实验

为了研究改性双基（CMDB）推进剂的动态起裂韧性,利用Hopkinson实验技术对CMDB推进剂的含切缝半圆形（NSCB）试件进行了冲击实验.运用脉冲整形技术,获得了CMDB推进剂的载荷-位移曲线.采用改进型柔度变化率法得到NSCB推进剂试件的起裂点.通过ABAQUS有限元软件对NSCB试件的无量纲结构因子进行了标定,并获得了NSCB的动态起裂韧性.利用电镜扫描设备分析了NSCB试件断裂的微观机理,研究结果表明:在加载率为5.41×105~8.94×105MPa·m1/2·s-1范围内,CMDB推进剂的动态起裂韧性具有显著的加载率敏感性.结合电镜扫描图分析,随着加载率的增大,断面内部颗粒微观结构的破坏形貌逐渐恶化,消耗的能量相应增大.      

长尾喷管传热及烧蚀特性研究

针对长尾喷管的热防护问题,本文采用实验和数值仿真结合的方法,研究了采用SQ-2推进剂固体火箭发动机复合结构长尾喷管的传热以及烧蚀特性。仿真程序采用基于格心的有限体积法,对流体域求解Navier-Stokes方程,对固体域求解热传导方程,湍流模型采用k-ω SST模型,通过保证热流密度大小相等、温度连续的方法实现耦合传热计算;在以上仿真计算方法的基础上建立了基于热解动力学的变热物性模型描述碳/酚醛材料的体积烧蚀,并采用简化的异相反应模型对喷管热化学烧蚀量进行计算。结果表明在收敛段和等直段区域压强变化较小,炭化层厚度和烧蚀量与温度分布相似, 并在距离等直段入口约1.6mm位置达到极小值。碳/酚醛材料炭化速率随时间减小,烧蚀速率随时间增长趋于稳定,而喉衬区域C/C复合材料烧蚀速率随时间增加。     

NEPE推进剂燃烧表面铝团聚物动态行为研究

针对NEPE推进剂燃烧表面铝团聚物的动态行为，在1,2,3MPa的N2环境中，对NEPE推进剂中铝团聚物在燃烧表面和脱离燃面后的动态燃烧过程进行了研究，提出了燃面与铝颗粒联合运动的简化理论模型。结果表明，铝颗粒的团聚分为堆积、聚集和团聚三个过程，铝团聚物在燃烧表面以及脱离燃面后均可能发生二次团聚，铝团聚物的二次团聚过程通常会形成大尺寸的铝团聚物。通过燃面与铝颗粒协同运动的简化理论模型，认为铝团聚的团聚时间受到压强、铝颗粒的体积分数和半径的影响。同尺寸铝团聚物的团聚时间随着压强的增大而较小，且在高压环境中的影响程度降低，与理论模型一致。同时压强会影响铝团聚物脱离燃面后的随流运动速率，同一压强条件下，直径越小的初始团聚物随流运动速率越快，直径越大的初始团聚物随流运动速率越慢。     

应用TTM网格研究冲压增程弹丸进气道内外流场

基于TTM贴体网格生成法,研究了某冲压增程弹丸超音速进气道内外流场数值模拟中单区域网格生成问题,流场求解采用了二阶隐式TVD格式。NACA0012翼型跨音速流场算例结果表明,所采用的算法和编制的程序可靠,能够适用于具有亚、跨、超音速三种流态共存的复杂流场数值计算。对超音速进气道数值模拟结果显示,得到的流场复杂波系结构是合理的,来流马赫数和攻角变化显著影响着进气道性能。     

EPDM薄膜橡胶包覆材料的粘-超弹本构模型研究

为了描述一种新型三元乙丙(EPDM)橡胶薄膜包覆材料的力学性能,进行了不同拉伸速率(0.1~5mm/s)下的单轴拉伸试验和多步松弛试验。研究了材料的应力-应变响应,建立了有限变形下的粘-超弹本构模型。模型由超弹与粘弹两部分组成:超弹部分基于多项式应变能函数推导获得,粘弹部分则采用Maxwell模型来描述。通过单轴拉伸试验和多步松弛试验分别获得粘弹参数和超弹参数。将所建模型预测结果与试验结果相对比。结果表明:模型能很好地描述材料在20%工程应变范围内的力学特性,预测结果与试验结果最大偏差仅为5.6%,验证了模型的可靠性。     

弧形翼-身组合体空气动力特性的数值模拟研究

利用高精度ＴＶＤ格式成功地完成了标准ＴＴＣＰ弧形翼－身组合体在马赫数为０．８～２．３６的湍流流场计算,通过对翼表面压力的积分求得自诱导转动力矩,与相应的风洞实验数据基本一致。流场的计算结果显示：自诱导转动力矩的变化主要是由于翼前缘和翼后半部分产生的法向力的差异造成的。最后进行了不同翼形参数的弧形翼－身组合体湍流流场计算,发现翼前缘削尖角角度增大和翼弦长增大,使总法向力增大。     

基于多物理场耦合的双脉冲发动机点火过程数值模拟

为研究双脉冲固体火箭发动机Ⅱ脉冲点火瞬态过程,发展一套多物理场耦合求解器。流体控制方程基于有限体积法求解,时间推进采用双时间步LU-SGS(Lower Upper Symmetric Guass Seidel)方法;固体推进剂表面温度基于耦合传热方法计算;结构动力学运动方程基于有限元方法离散,采用经典的Newmark格式进行时间推进,流固耦合采用松耦合算法,并通过算例验证求解器的可靠性。计算结果表明:该求解器能够数值模拟Ⅱ脉冲启动过程中的点火药气体冲击、燃气非定常流动及金属膜片机械响应过程,获得金属膜片的破裂时间和压强;且随着点火质量流率增加,推进剂装药首次点燃时间和金属膜片破裂时间变短,膜片破裂压强降低;金属膜片破裂时间和压强不仅与作用在其表面的压强载荷大小相关,而且与压强载荷加载的过程相关;金属膜片厚度越薄,膜片破裂时间越短,膜片轴向位移越大,膜片破裂压强越低。