基于DSPN的航天测控系统任务可靠性仿真建模

提出了基于确定与随机Petri网（deterministic and stochastic Petri nets,DSPN）的航天测控系统（tracking,telemetry and command,TT&C）任务可靠性定量分析方法,旨在对相关航天测控方案进行可靠性预计.通过对TT&C系统任务剖面进行时序弧段划分,考虑实际系统中测控单元阶段依赖、单元故障可修以及各单元参与任务起止时间不同等其他建模方法难以处理的复杂因素,建立了“单元层-系统逻辑层-阶段层”3层相互关联的TT&C系统任务可靠性DSPN模型.通过对模型仿真运行,实现了对给定测控方案下TT&C系统任务可靠性定量化评估.分析表明:仿真结果随着仿真次数增加逐渐收敛,与Markov解析方法求得的精确值对比误差控制在1%以内.      

弹塑性梯度材料热机耦合分析的有限元方法

陶瓷/金属功能梯度材料(FGMs)的组份及性能梯度分布,在航空航天高温热结构中有重要应用。材料性能在空间上的连续变化给这类结构强度设计的有限元分析造成了一定困难,针对上述问题,提出一种热弹塑性梯度有限元方法。该方法通过编写有限元软件ABAQUS的材料用户子程序(UMAT)来实现,在UMAT中将FGMs的弹塑性本构模型参数定义为空间坐标的函数;用热弹塑性梯度有限元方法对陶瓷-FGMs-金属典型的三明治板结构进行热机械分析。结果表明:FGMs组分布函数对其热应力分布有很大的影响,特别是在三明治结构的界面上,影响更为突出,这和现有文献的实验结果相符。     

基于库伦摩擦效应的箔片轴承极限承载力理论与试验

基于有限单元法建立了考虑库伦摩擦的波箔型径向气体箔片轴承的箔片结构模型,采用有限差分法和有限单元法耦合求解Reynolds方程和气膜厚度方程,通过求解轴颈达到极限偏心率时的轴承极限承载力,研究了箔片结构库伦摩擦效应对轴承极限承载力的影响规律,并搭建了轴承极限承载力测试试验台,利用温度法测量了两个具有不同轴承壳内表面粗糙度的波箔型径向气体箔片轴承的轴承极限承载力.通过对比分析仿真结果与试验结果表明:轴承壳圆柱孔内表面粗糙度为0.4μm的轴承在10000r/min和20000r/min下,轴承极限承载力分别为15.5N和42.3N;而表面粗糙度为1.6μm的轴承极限承载力为10.9N和29.6N,这是由于波纹箔片和轴承壳体之间的库伦摩擦力增大了波纹箔片的刚度,因此增大箔片结构摩擦因数使得轴承极限承载力降低,并且仿真结果变化趋势和试验结果变化趋势吻合.      

基于概率的航空发动机飞行换算率改进算法

针对斯贝MK202发动机应力标准对零部件工作寿命评定中未包含各种随机因素的影响,考虑零部件疲劳性能参数的分散性和发动机整机载荷谱的随机性,将实测数据统计分析、零部件应力有限元分析以及蒙特卡罗模拟分析有机结合,提出了基于概率的航空发动机飞行换算率改进算法,以解释航空发动机零部件实际工作寿命存在的差异.经范例验证表明:确定性寿命评定方法的飞行换算率计算结果与改进算法下的平均值相当;对于单次飞行换算率,改进算法的单次飞行换算率平均值比确定性方法的计算结果小2%;对于综合飞行换算率,改进算法的综合飞行换算率平均值比确定性方法的结果小1%;且在改进算法下单次飞行换算率和综合飞行换算率服从正态分布.     

基于CFD/CSD方法的蜻蜓柔性翼气动特性分析

给出了一种基于计算流体力学/计算结构力学（CFD/CSD）的双向流固耦合方法.通过交替数字二叉树（ADT）搜索技术识别流固网格之间的宿主-受体关系.采用局部插值算法完成两套网格系统之间的数据交换,并使用Delaunay图映射方法来完成气动网格的移动.将自编的非线性结构有限元程序、接口程序与南京航空航天大学（NUAA）微型飞行器中心的流体计算程序3D2MUFS相连接,应用于蜻蜓柔性翼拍动飞行的气动计算中.计算结果表明:柔性变形使得蜻蜓翼的时均举力系数从0.31提高到0.53,时均推力系数从0.07提高到0.13,证实了柔性变形能改善扑翼的气动性能.      

含分层损伤复合材料层合板振动特性

针对复合材料层合板分层损伤区域上、下子板的畸变模态,采用自定义矩阵单元模拟其损伤区的接触刚度,建立了一种合理的层合板分层损伤有限元振动分析模型;在此基础上研究了分层深度和分层大小对复合材料层合板振动特性的影响.数值模拟结果与实验结果的对比表明:采用的自定义矩阵单元可以有效地模拟层合板的分层损伤,模态计算值与实验值的最大误差为10.67%,最小误差为0.34%;分层深度和分层大小对复合材料层合板振动特性有较大影响,随分层深度变化,固有频率最多下降50%;随分层大小变化,前4阶固有频率最多下降12%.      

飞机橡皮囊成形零件毛坯精确预示方法研究

 橡皮囊成形是飞机钣金零件成形的主要方法之一,其初始毛坯形状的预示是橡皮囊成形中的一个难点问题。针对橡皮囊成形的工艺特点,提出了一种基于一步逆成形有限元法的飞机橡皮囊成形钣金零件毛坯快速精确预示方法。首先将设计好的飞机钣金零件作为最终构形,然后采用一步逆成形有限元法进行初始毛坯的快速反向模拟,最后将通过该方法获得的毛坯外形应用到实际飞机钣金橡皮囊成形中,实验结果表明该方法可以快速精确地预示出飞机橡皮囊成形钣金零件的初始毛坯形状。     

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

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

2/3层板结构空心风扇叶片设计与对比分析

针对空心结构风扇叶片设计技术在航空发动机上的应用,分析了2/3层板空心叶片结构特点,提出了1种基于气动数据源的空心叶片结构设计方法,并以某风扇叶片为例开展了2种结构空心叶片设计,有限元对比分析结果表明:2种结构空心叶片在静强度、振动特性方面各有优势,但2层板结构空心叶片可设计性更优,应用前景更好.     

基于SPN的跑道利用率模型分析与研究

采用随机Petri网(SPN)建立了跑道系统运行约束方程;通过构造同构马尔科夫链建立了系统处于各种状态时的稳态概率表达式;改变变迁的触发速率,分析跑道利用率的变化情况。证明了模型的可达性、有界性,且不存在死点。