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

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

基于表面粘结损伤的复合固体推进剂细观损伤数值模拟

为更准确地模拟固体推进剂的细观损伤,研究复合固体推进剂的非线性力学性能,采用分子动力学方法对高填充比复合固体推进剂细观模型进行建模,验证了Surface-based cohesive方法的可行性,并替代经典粘结单元模拟了复合固体推进剂细观颗粒与基体之间界面损伤的产生和聚合.通过对比固体推进剂细观模型的计算结果,研究了复合固体推进剂颗粒粒径、位置的随机分布对固体推进剂细观损伤及宏观力学性能的影响.结果表明,采用Surface-based cohesive方法对固体推进剂细观模型较易进行有限元前处理,易收敛,计算结果比较合理.     

复合固体推进剂颗粒脱湿损伤参数反演

为确定复合固体推进剂颗粒脱湿损伤参数,采用分子动力学方法建立了复合固体推进剂颗粒夹杂模型,根据Surface-based cohesive方法,在高氯酸铵(AP)颗粒与基体之间的界面处设置接触损伤,对复合固体推进剂的细观损伤过程进行数值仿真计算。将数值仿真得到的应力-应变曲线与试验曲线进行对比,建立颗粒脱湿损伤参数的优化目标函数,通过分步迭代计算,对颗粒脱湿损伤参数进行反演识别。结果表明:根据反演参数计算复合固体推进剂细观损伤过程,颗粒夹杂模型表现出的宏观应力-应变关系与试验结果吻合良好,说明双折线损伤模型可以近似表征复合固体推进剂的脱湿损伤。     

固体火箭发动机药柱的有限元网格划分研究

讨论了混合网格划分功能在对固体火箭发动机药柱进行有限元分析中的应用,对同一推进剂药柱进行了两种不同形式的网格划分,并对计算时间进行了比较,提出了推进剂药柱网格划分应注意的一些细节方面。     

基于细观力学的复合固体推进剂模量预估方法

应用分子动力学方法得到了复合固体推进剂颗粒在基体内随机分布的填充模型。针对推进剂所含颗粒大小不同,采用细观力学有限元法和Mori-Tanaka分析法相结合的方法计算了推进剂的模量,计算所得模量与不同颗粒体积分数推进剂的单向拉伸试验结果比较吻合。此外,该方法可以获得推进剂细观结构应力,从而为推进剂的细观损伤分析奠定基础。     

基于原位X射线成像的推进剂损伤演化表征

为实现对硝酸酯增塑聚醚（NEPE）固体推进剂的细观损伤及其演化行为的可视化表征，基于自主研制的原位力学试验系统与第3代高分辨同步辐射X射线三维成像技术，对拉伸速率为0.1 mm/s下的NEPE固体推进剂单调拉伸过程进行三维成像原位观测和表征，获取初始状态固体推进剂的细微观形貌及其随拉伸载荷的演化特征，提取并分析其内部典型损伤的体积与球度随载荷的演化规律。结果表明：同步辐射原位成像技术能够准确获取NEPE固体推进剂的细观结构特征，可以基于灰度差异对固体推进剂的AP颗粒、Al颗粒、基体以及缺陷等实现特征的准确识别。研究发现，NEPE固体推进剂内部缺陷主要有2种，一种为颗粒内的孔洞，一种为AP颗粒/基体界面的初始脱湿。推进剂细观损伤首先表现为初始的界面脱湿形成的孔隙：在拉伸载荷较小时，推进剂的损伤形式主要表现为较大的AP颗粒脱湿形成的孔隙；在拉伸载荷较大时，除AP颗粒脱湿外可以观察到Al颗粒的脱湿，大量的AP颗粒脱湿后形成的孔隙相互融合，最终导致固体推进剂宏观断裂。对固体推进剂内部孔隙定量化表征的结果表明，孔隙体积随着拉伸载荷增大而增加，而球度变化与初始缺陷相关，有初始缺陷时，球度呈现单调减...     

复合固体推进剂过渡相(中间相)力学模型

在总结归纳复合固体推进剂力学性能研究的基础上,吸取键合剂作用模型的优点,结合新近提出的网络与形态结构理论,并引入颗粒充填复合材料中的成熟理论,提出了一个描述复合固体推进剂力学性能的过渡相（中间相）模型。通过力学,数学推导建立了复合固体推进剂的结构因素与力学性能表征参数（σ,ε）的定性关系,这与验证主前人的结果基本一致,其中提出的关于过渡相模量与基体模量匹配的调节方法,为进一步提高复合固体推进剂的力     

含铝固体复合推进剂颗粒团聚现象

含Al固体复合推进剂在燃烧过程中，燃烧表面存在Al颗粒团聚现象，其对火箭发动机性能产生重要影响。通过分子动力学算法建立固体复合推进剂三维颗粒微观模型，实现了模拟三维空间内AP颗粒和Al颗粒的随机分布特征。并建立了颗粒的团聚模型，分析研究Al颗粒在随机填充模型中的团聚现象，得到了不同的临界分离距离下Al颗粒粒径分布规律，并将之与实验数据进行对比，最终总结出了合适的临界分离距离与Al颗粒直径的经验公式。然后总结出了Al颗粒在气固交界面发生团聚后粒径的分布规律，设计了Al颗粒在气固交界面的Rosin-Rammler概率分布函数，分析得到其均匀性系数与特征系数分别为1.445 3与86.49，可用于数值模拟计算固体燃料燃烧面退移过程中表面Al颗粒喷射的初始粒径。     

数字全息在固体推进剂铝燃烧三维测量中的应用研究

为获得固体推进剂铝颗粒动态燃烧的精细化过程,将数字全息技术用于固体推进剂铝燃烧三维动态过程测量中,成功地解决了传统光学显微成像法景深过小的问题。搭建了固体推进剂铝燃烧数字全息测量系统,分别在0.1MPa和1.0MPa压强条件下开展了推进剂中铝颗粒燃烧过程的三维测量工作,获得了铝颗粒动态燃烧过程的全息图像。研究结果表明:数字全息技术能够获得不同截面处的颗粒信息,并可跟踪颗粒的动态燃烧过程,准确地获得颗粒粒径信息,真正实现三维动态测量;全息法能够清晰地分辨出动态燃烧过程中十几微米至几百微米量级的铝颗粒,测量误差小于8%;跟踪单一铝颗粒的动态燃烧过程,可获得颗粒粒径、空间分布、粒径分布、燃烧火焰区域、运动速度以及氧化帽的动态生成过程等变化规律。     

基于原位拉伸的推进剂/衬层界面力学性能研究

为了解释某种HTPB复合固体推进剂/衬层粘接试件载荷-位移曲线的"双峰"特征,设计了一种基于推进剂/衬层微型试件的原位拉伸试验方法,根据界面细观破坏形态与载荷变化过程,提出了界面处颗粒脱湿及基体断裂过程分别对应两个载荷峰的假设。采用改进的并联Maxwell元件模型对界面断裂行为进行了模拟计算,重现了断裂过程的载荷-位移曲线的"双峰"特征,验证了假设的合理性。通过模拟计算实测曲线,给出了推进剂基体与颗粒之间的近似粘接强度等参数,为推进剂/衬层粘接系统细观材料参数计算方法提供了一种参考。     

二维任意域内基于节点的局部网格生成算法

凸域内基于节点的局部网格生成算法,克服了基于节点的有限元方法的网格生成可能产生的不一致性。将该基于节点的局部网格生成算法的适用范围拓展到二维任意域。另外,提出了通过使用约束Delaunay路径来划分任意域的区域划分算法,该算法使得在并行实现网格生成的过程中各处理器之间无需通信,从而大大提高了节点给定情形下有限元方法网格生成的并行效率。     

快速实现二次曲面片三角网格剖分的新方法

三维表面的三角网格剖分,对图形学和数值分析来说,都是极其重要的一项预处理技术.利用二维动态、带约束的Delaunay三角网格剖分算法,可以有效地实现二维空间中带边界约束的点集的三角网格剖分.提出了对二次曲面进行快速三角网格剖分的一种新方法.     

3D-Computaional Mesh Generation around a Propeller by Elliptic Differential Equation System

In this paper, an analytic method is applied for the generation of 3D mesh system in order for Navier-Stokes equations around ATP (Propfan). One of the advantages of this method is that mesh lines have strong differentiability. The differential equation used is Poisson type and the right hand side is called control function, and this is able to control the degree of mesh line clustering. Here, the form of control function was contrived to cluster near the solid surfaces. By this method, several mesh lines are laid in the boundary layer above the blade surfaces.     

二维有限元网格全自动生成方法

主要研究二维有限元网格全自动生成方法。提出一种新的节点生成方法 ,通过实例验证该方法可靠。采用三角形形状因子评定三角形网格质量 ,可减少三角化的时间 ,优化三角网格质量。     

大型筒段构件自动铺带成型技术

在研究自动铺带成型技术的基础上,提出了一种用于复合材料筒段构件自动铺带成型轨迹规划的设计方案。根据自动铺带成型系统的独有特点,提出专用于自动铺带成型系统的数控(NC)代码约定。在研究自动铺带成型CAD/CAM技术的基础上,开发针对大型筒段复合材料构件的专用CAD/CAM软件,并在自主开发的自动铺带成型设备上进行铺带成型实验验证。经实验验证,此自动铺带机专用CAD/CAM软件的开发相当成功,满足了大型复合材料筒段构件的实际生产需求。     

飞机简化模型水上迫降入水的数值研究

在飞机水上迫降过程中,其尾部吸力对姿态角的变化具有重要影响.基于Fluent软件,利用动网格技术、用户自定义函数(UDF)、六自由度(6DOF)模型、流体体积(VOF)模型,首先对二维圆柱、三维椭圆抛物体垂直入水以及三维平板水平划水进行数值模拟,并与试验结果对比验证模拟方法的可靠性;然后对NACA TN 2929简化飞机模型进行水上迫降数值模拟.研究发现:尾部吸力使飞机俯仰姿态角发生剧烈变化,入水初期尾部吸力产生的抬头力矩使飞机姿态角快速增加,当姿态角达到最大时,尾部吸力产生的抬头力矩很小,在重力作用下飞机姿态角快速减小.本文所采用的三维数值模拟方法能够较好地模拟飞机水上迫降过程中的吸力和姿态角变化.     

基于网格参数化变形的单级涡轮多学科可靠性设计优化

基于自由变形技术发展了一种可以协调实现流场、结构分析网格参数化变形的方法,并实现单级涡轮的多学科可靠性设计优化。针对流场分析网格展向积叠的特征,将三维网格变形转化为不同叶高二维截面的网格变形;根据二维流场分析网格的拓扑结构,设计二维控制面,基于叶型设计参数建立叶型控制点、其他控制点伴随叶型控制点移动,实现了网格的参数化变形,避免了网格畸变、提高了网格变形的质量;不同叶高的二维控制面组成控制体,实现三维流场分析网格的参数化变形。采用相同的控制体进行结构分析网格的变形,实现了耦合界面网格协调变形。基于Kriging代理模型进行单级涡轮的多学科可靠性优化,在满足可靠性约束的情况下效率提高4.97%、寿命提高40.86%,证明了方法的有效性。      

针刺C/C圆筒构件内/外压试验与仿真研究

根据某扩张段采用的针刺C/C复合材料制备圆筒构件试验件,设计外压和内压试验工装,对其开展内外压试验研究,总结结构应力分布规律,并与有限元仿真结果进行对比,仿真结果与试验结果一致性较好,获得了准确高精度3D圆筒构件水压仿真计算模型。试验及仿真结果对于改进固体火箭发动机针刺C/C复合材料喷管扩张段结构设计,修正材料本构模型参数,预估结构安全裕度有着重要的参考价值。     

一类二维非结构网格DSMC方法的实现策略及其应用

研究了一类二维非结构网格DSMC方法的实现策略。在数据结构方面，设计了局部化的数据组织方式，节约了内存与计算时间。发展了一种跟踪模拟分子迁移的算法，该算法仅需少量的逻辑运算与代数运算，不仅可以快速跟踪模拟分子在网格之间的迁移，而且可以准确判别分子与物面是否相互作用，搜索过程中的附带信息给出了分子与物面碰撞的精确时间与位置，避免了重新计算。引入碰撞距离的思想，既减少存储又保证正确的模拟结果。在程序编制过程中，我们充分展现了Fortran90高级语言的主要特性，引入动态数组、指针、链表以及派生类型数据，编制了计算程序。最后对过渡流域高超声绕流进行了数值试验。     

基于工程应用的有限元网格划分研究

首先介绍和分析了划分有限元网格的基本原则,然后基于MSC.Patran有限元软件,结合3个不同的工程实际问题进行网格划分分析,总结了划分有限元网格的一般过程,最后对划分有限元网格进行了小结。