HTPB复合固体推进剂非线性本构研究

采用Ｓｃｈａｐｅｒｙ蠕变型非线性粘弹性本工用Ｍ－Ｃ方法,分析蠕变实验数据,研究了ＨＴＰＢ复合固体推进剂非线性本构关系。结果表明,该方法可以减少加卸载时数据波动的影响,降低参数拟合对实验条件的要求。Ｓｃｈａｐｅｒｙ非线性模型较好地反映了ＨＴＰＢ复合固体推进剂的力学性能,可以用于固体推进剂的力学分析。     

粉末合金(FGH95)蠕变/疲劳交互作用下寿命预测的损伤力学有限元分析

采用基于损伤力学的蠕变—疲劳交互作用下的寿命预测模型,应用损伤等效应力进行三维应力状态下的损伤计算,并考虑了压缩时闭合效应,利用ANSYS的二次开发工具APDL和UPFs开发程序,把基于损伤力学的寿命预测方法与ANSYS的结构分析结合起来,实现了对构件的损伤计算和寿命预测。针对粉末合金材料易含夹杂等初始缺陷的特点,提出了在寿命计算中通过单元初始损伤模拟初始缺陷对寿命影响的处理方法,探讨了考虑初始缺陷条件下的寿命预测,并利用试验对计算结果进行了对比验证。      

考虑多点保形的结构拓扑优化设计方法

为了抑制飞行器结构局部区域内多个控制点的翘曲变形,提出考虑结构保形的拓扑优化设计新方法。基于由局部区域内保形控制点生成的人工附加弱单元(AWE),定义AWE的变形能约束函数,从而引入了对结构局部区域的多点保形约束,定量衡量并抑制了该区域内结构的弹性翘曲变形程度。进一步讨论并分析了保形约束与整体刚度之间的消长关系以及优化过程中可能存在的传力路径畸变现象。拓扑优化数值算例的设计结果对比表明,保形拓扑优化方法相比现有拓扑优化方法,能够有效抑制多点间的相对位移和翘曲变形,实现保形设计效果。     

整体复合材料结构分层特性研究进展（一）

整体复合材料件越来越多地被用作现代航空航天器的主承力结构,而分层是其最易出现的破坏形式。本文作为整体复合材料结构分层特性研究综述的第一部分,重点介绍界面断裂力学和应力分析方法。论文首先对复合材料失效和分层研究进行了简要的回顾,然后详细介绍了界面断裂力学发展的历史、基本公式和在复合材料分层研究中的应用,并且介绍了几种界面裂纹尖端场的求解方法——J积分、修正的裂纹闭合积分（Modified Crack Closure Integral）和M积分法。其次,对应力分析在复合材料分层研究中的应用进行了简要的回顾。     

局部凸起对可压缩平板边界层稳定性的影响

在采用直接数值模拟(DNS)方法计算得到来流马赫数为4.5的带局部凸起的可压缩平板边界层基本流的基础上,基于线性稳定性理论(LST) 进行稳定性分析,通过求解扰动方程研究了不同频率的扰动波与不同高度的局部凸起的相互作用,定义了穿透系数来研究局部凸起对流动稳定性的影响.研究结果表明:当局部凸起高度小于0.2个边界层厚度时局部凸起对来流的影响是局部的;局部凸起使频率小于最不稳定波频率的扰动波更不稳定,对频率大于最不稳定波频率的扰动波有稳定作用,其中最不稳定波频率是指局部凸起中心位置处平板边界层中第2模态最不稳定波的频率;穿透系数定量刻画了局部凸起对稳定性的影响,高度为0.2个边界层厚度的局部凸起引起e-N方法中N值的修正量达到0.8,接近转捩位置处N值的10%,在转捩预测中必须考虑.      

Damage characteristics and constitutive modeling of the 2D C/SiC composite: Part Ⅱ–Material model and numerical implementation

In this work, a macroscopic non-linear constitutive model accounting for damage, inelastic strain and unilateral behavior is proposed for the 2D plain-woven C/Si C composite. A set of scalar damage variables and a new thermodynamic potential expression are introduced in the framework of continuum damage mechanics. In the deduced constitutive equations, the material's progressive damage deactivation behavior during the compression loading is described by a continuous function, and different deactivation rates under uniaxial and biaxial compression loadings are also considered. In damage evolution laws, the coupling effect among the damage modes and impediment effect of compression stress on the development of shear damage in different plane stress states are taken into account. Besides, the general plasticity theory is applied to describing the evolution of inelastic strain in tension and/or shear stress state. The Tsai–Wu failure criterion is adopted for strength analysis. Additionally, the material model is implemented as a user-defined material subroutine(UMAT) and linked to the ABAQUS finite element software, and its performance is demonstrated through several numerical examples.     

基于多孔介质局部非热平衡方法的刷式密封耦合传热特性

建立了刷式密封流动传热与摩擦接触的多物理场耦合分析的计算模型.推导了描述刷式密封在局部非热平衡条件下耦合传热过程的双能量方程,该方程既考虑了刷丝与气流的对流换热,同时也考虑了刷丝束内导热过程的各向异性.采用所建立的耦合计算模型数值预测了刷式密封在不同压比和转速下的温度分布.结果表明:刷式密封的摩擦热效应使刷丝与转子相接触的附近区域形成局部高温,刷丝的高温区域主要集中在围栏高度以下;泄漏气流与刷丝间强烈的对流换热使刷丝的温度分布和气流的温度分布非常相近;刷丝的最高温度随压比增大而升高,但围栏高度以下区域的平均温度随压比增大而降低;刷丝的最高温度和围栏高度以下区域的平均温度均随转速升高而升高;摩擦热效应引起的泄漏气流的温升效应使刷式密封泄漏量随转速升高而减小.      

针刺陶瓷基复合材料损伤本构模型及构件应力分析

基于连续介质损伤力学方法提出了适用于针刺陶瓷基复合材料的各向异性损伤本构模型.该模型考虑了拉伸与剪切损伤的影响.通过试验得到了材料的拉伸和剪切的应力-应变曲线,采用曲线拟合获得了本构模型的参数.应用用户子程序技术将本构模型写入商用有限元软件,计算所得拉伸和剪切的应力-应变曲线与试验曲线比较吻合,最大误差分别为5.62%和1.47%.使用该损伤本构模型及弹性模型计算了航空发动机尾喷管调节片在气动和温度载荷下的力学响应,两者计算所得应力分布相似,但损伤本构模型计算值明显小于弹性模型计算值.      

采用火焰面/反应进度变量方法模拟湍流燃烧

在自主开发的软件平台上添加了火焰面/反应进度变量方法,选取了两种反应进度变量分别对甲烷/空气的同轴射流燃烧室进行了数值模拟.将计算结果与稳态火焰面模型、文献参考值以及实验值进行了对比研究,结果表明:采用火焰面/反应进度变量方法能捕捉到稳态火焰面模型所不能预测到的火焰抬举等非稳态现象,得到的结果与文献参考值以及实验值出入很小,能准确反应真实的燃烧过程;另外不同反应进度变量的定义对计算结果有较大的影响,采用火焰面/反应进度变量方法对主要中间产物一氧化碳的预测比文献参考值更好,更接近于实验值.      

改进节点积分的无单元Galerkin法及其在流动问题中的应用

无单元Galerkin法需要在背景网格上积分,计算量大,且在求解对流占优问题时会出现非物理的数值伪振荡现象。为此,基于局部Taylor展开思想,采用节点处的局部Taylor展开计算积分,建立了局部Taylor展开积分无单元Galerkin法。该方法同时解决了标准的无单元Galerkin法计算量大和对流占优时会出现数值伪振荡的问题。一维定常对流扩散方程和二维Burgers方程的求解说明了该方法的有效性,且计算效率远高于无单元Galerkin法。