膨胀效应对激波/湍流边界层干扰的影响

采用直接数值模拟方法对来流马赫数2.9,30°激波角的入射激波与膨胀角湍流边界层干扰问题进行了数值研究。入射激波在壁面上的名义入射点固定在膨胀角角点,膨胀角角度分别取为0°、2°、5°和10°。通过改变膨胀角角度,考察了膨胀效应对干扰区内复杂流动现象的影响规律和作用机制,如分离泡、物面压力脉动特性、膨胀区湍流边界层和物面剪切应力脉动场等。研究发现,膨胀角角度的增大使得分离区流向长度和法向高度急剧降低,尤其是在强膨胀效应下分离泡形态呈现整体往下游偏移的双峰结构。物面压力脉动功率谱结果表明,膨胀角为2°和5°时,分离激波的非定常运动仍表征为大尺度低频振荡,而膨胀角为10°,强膨胀效应极大地抑制了分离激波的低频振荡,加速了下游再附边界层物面压力脉动的恢复过程。膨胀区湍流边界层雷诺剪切应力各象限事件贡献和出现概率呈现逐步恢复到上游湍流边界层的趋势,Görtler-like流向涡结构展向和法向尺度变化剧烈,同时在近壁区将诱导生成大量小尺度流向涡。此外,物面剪切应力脉动场的本征正交分解分析指出,膨胀效应的影响体现在低阶模态能量的急剧降低从而使得高阶模态的总体贡献相对升高。     

带非线性支撑的转子有限元模型求解方法

用数值方法研究了非线性支撑的柔性转子系统的动学行为,提出了一种将有限元与非线性支撑结合的模型和求解方法。利用有限元法(FEM)构建转轴和转盘部分的模型,通过矩阵进行组合;利用离散元方法对包含滚动轴承和挤压油膜阻尼器(SFD)的支撑部分进行建模,此部分包含4个单元,分别为轴承内圈、外圈、SFD内圈和支撑鼠笼。有限元部分和离散元部分通过轴端节点相连,仿真过程中轴端位移传递给非线性支撑部分,支撑部分通过位移计算得到的非线性力反过来作用于有限元转子轴端部分。为了耦合求解有限元转子和非线性支撑组成的数学模型,提出了一种综合的迭代求解方法,克服传统的有限元求解方法对轴端隐性非线性支撑的求解局限性。由于转轴部分采用了Timoshenko梁单元建模,对比与简单转子模型,可以考虑陀螺力矩和轴的柔性特征,更能体现非线性支撑对振动真实影响。在建立的20个轴单元的有限元转子模型中,非线性响应更多体现在靠近非线性支撑的节点1和节点21处,响应频谱中靠近轴端的节点能体现出滚动轴承的2倍和3倍变柔振动频率。     

潜入式喷管喉衬界面间隙优化设计

为了解决固体火箭发动机潜入式喷管喉衬在热载荷与内压联合作用下界面间隙优化问题,数值模拟方法求解了喷管热结构响应.采用流体仿真软件,计算了潜入式喷管流动过程的稳态流场,获取了高温燃气流动参数与固体壁面对流换热系数.并采用三维有限元结构计算平台,编写了非均布壁面压力载荷与非均布对流换热系数子程序,求解了喉衬前后搭接界面间隙...     

MTPS蜂窝夹芯结构传热性能及热应力分析

对于金属防热结构的蜂窝夹芯结构的稳态情况,基于材料的全灰体假设,同时考虑热传导和热辐射两种传热形式对温度场的耦合作用,利用热流量守恒建立了蜂窝芯层温度场的非线性积分方程。离散化后利用数值方法得到方程组的数值解。对于美国兰利研究中心的实验结果,与本文方法的对比计算结果基本吻合。进一步,利用计算结果讨论了给定面板温度边界情况下,下面板、柱体层的灰度、蜂窝结构长径比对夹芯温度场的影响。并根据温度场和近似的应力分析模型,用半解析结果讨论了稳态情况下蜂窝芯层上的热应力。     

模态综合的直接变换法

基于固定界面模态综合的Craig-Bampton方法,提出了模态综合的直接变换方法,该法在形式上更简洁,适用于大规模动力学系统的模态分析和动力学模型减缩。对采用直接变换法模态综合过程的计算误差做了分析,给出了计算误差界,阐明了模态综合后总体结构的固有频率与子结构保留主模态个数之间的关系。针对子结构界面自由度太大的情况,提出了直接变换法的分部算法,不需组集系统总刚度、总质量矩阵,仅在子结构的基础上完成系统的减缩。计算实例表明该算法的有效性及计算的稳定性。     

高精度反射面板制造中边值问题的求解方法(英文)

采用三层铝板两层蜂窝芯的夹层结构制造高精度反射面板是一种先进、复杂的技术。其中,单层铝板的受力分析和夹层板的回弹计算是成形精度的根本保障。利用Fourier级数表示单层板应力函数、利用幂级数表示夹层板挠度函数,分别对高精度反射面板制造中单层板压贴和夹层板回弹两类边值问题进行了求解。与有限元模拟结果对比表明,本文提出的单层板压贴膜应力计算Fourier级数法具有较高精度;与试验结果对比表明,本文提出的夹层板回弹量幂级数法具有满意精度。     

有、无障碍物爆震管起爆过程的CE/SE模拟

用时空守恒元和解元(CE/SE)法模拟有无障碍物情况下爆震管二维爆震波的形成和传播过程,分析爆震波的不同机理,研究在较低能量下,点火区温度、爆震管管径以及障碍物对爆燃-爆震转捩过程(DDT)的影响.控制方程采取无量纲形式,用拓展的隐式梯形法实现刚性源项积分.研究表明,点火能量与点火压力相关,在较低的点火能量下,可以通过提高点火温度、减小爆震管径及加入障碍物加强爆震过程,减小爆燃-爆震转捩过程DDT时间和距离.为研究爆震波点火和DDT过程提供了理论依据.      

结冰翼型表面明冰的压力分布和应力计算

在不同来流速度和攻角下,比较分析了冰型外部压力分布的变化规律,计算分析了外部压力分布对表面冰型的力学特性影响规律,研究了冰型和蒙皮界面间最大剪切应力的变化规律,并结合剪切强度随界面温度变化的实验曲线,分析了最大剪切应力对冰层脱落的贡献.计算结果表明:最大剪切应力是剪切强度的5.48%,因此气动力对冰层的脱落贡献不大;当开启热除冰装置后,由于蒙皮表面温度的不断上升,冰型和蒙皮界面间剪切强度的下降,当表面温度为-2℃时,最大剪切应力超过剪切强度的10%,此位置可能率先产生破裂,从而加速冰层的脱落.      

实测环境温度下固体火箭发动机二维温度场及应力场分析

监测了某舰载导弹发射贮运箱内温度,用有限元法求解了实测环境温度作用下发动机内二维温度分布和应力分布情况。结果表明:发动机内温度与贮箱内交变环境温度相比,温度峰值有所衰减,相位有所滞后;发动机星尖处应力最大,是比较容易产生裂纹的地方。     

高温合金松弛过程回归分析方法

基于线性过程回归分析理论,建立了高温合金应力松弛过程回归分析方法,给出了松弛过程的回归方程和高置信水平、高可靠度的单侧置信下限曲线.高温合金松弛过程回归分析方法可以将应力松弛过程作为一个整体进行回归分析,充分利用了不同时刻剩余应力试验数据间的纵向信息,具有信息量大、精度高的特点.而且当试样数量较少时,传统方法得到的剩余应力置信下限曲线可能会发生畸变,而该方法则可以避免这一现象的发生.工程应用表明,在试样数相同的情况下,该方法比传统方法具有更高的精度,而在精度相同的情况下,则可减少大量试样.