径厚比对锥柱形金属膜片屈曲和后屈曲状态的影响

通过有限元方法研究了半径及径厚比对锥柱形金属膜片翻转过程的影响,并选取典型的半径及径厚比样本点对分析结论进行了数值验证.研究结果表明:膜片半径太小或太大均表现出一定的尺寸效应,小半径膜片的尺寸效应更显著;半径为203~403mm的膜片其屈曲和后屈曲状态表现出明显的一致性;对于确定半径的膜片,径厚比太大,膜片承载能力太小而导致翻转效率较低,径厚比太小,膜片屈曲载荷呈指数级显著增加而不易翻转,径厚比为160~400的膜片其翻转效率较好,并且后屈曲状态表现出一致性;选取典型的样本点对结论进行了数值验证,一致性较好.      

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

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

螺旋管内超临界航空煤油流动阻力特性

研究了系统压力、质量流速及螺旋管结构形式对超临界压力下航空煤油RP-3在螺旋管内的流动阻力特性.实验结果表明:螺旋管局部阻力系数变化曲线由螺旋管内流体临界雷诺数、拟临界温度分为明显的3个部分:工质温度低于拟临界温度且管内雷诺数小于螺旋管临界雷诺数时,局部损失系 数与雷诺数和螺旋直径与管径比D/d_in相关;流体温度低于拟临界温度且雷诺数大于临界雷诺数时,局部阻力系数只与D/d_in相关; 流体温度大于拟临界温度时,局部阻力系数除与雷诺数和D/d_in相关外,同时还与密度、黏性的大幅变化相关.另外,基于实验结果,提出了一种用压力、温度和螺旋直径与管径比进行修正的局部阻力系数关联式,拟合结果与实验结果相比,具有较好的一致性.      

声载荷作用下高温薄壁结构响应特性分析

先进的航空航天器表面结构暴露在严酷的工作载荷环境中,包括复杂的机械力载荷、压力载荷、声载荷和热载荷等,航空航天器表面结构可简化为薄壁结构,在复合载荷作用下结构以非线性方式响应,呈现出复杂的响应特性。首先以热弹性力学、板壳理论及结构稳定性原理为基础,建立热及噪声载荷联合作用下薄壁板运动模态方程,讨论了薄壁结构跳变响应的机理,进而运用等价线性化方法求解模态方程,在此基础上分析了热及噪声载荷对薄壁结构屈曲的影响,进而探讨了热及噪声载荷作用下薄壁结构的非线性响应统计特性。计算了四边简支高温钛合金薄板在声载荷下四个关键点处的均方应变,为进一步开展薄壁结构声疲劳寿命估算和强度设计奠定了基础。     

局部特征尺度分解改进算法

针对局部特征尺度分解估计包络曲线时未考虑其凹凸性的问题,提出了一种改进的算法。首先,利用3次样条插值方法进行均值点估计,结合前向以及后向插值2种方法,给出了双向形式;其次,参照标准算法,给出了具体的改进算法实现步骤;最后,通过仿真实例验证了有效性。     

热声载荷作用下薄壁结构非线性响应分析和试验验证

针对热声载荷作用下薄壁结构大挠度非线性响应问题,开展了固支约束金属薄壁板结构热声激励试验及数值模拟分析。通过计算结果与试验结果对比,表明两者结果存在一致性,进而验证了薄壁板在热声载荷作用下动态响应计算方法和数值模型的有效性。在此基础上,针对加筋板结构完成了多种热声载荷组合作用下的动力学响应计算,获得了时域位移响应。重点对该结构在后屈曲状态下的3种典型振动形式进行分析,总结出热载荷与声载荷之间的相对强度决定了板的不同跳变形式。采用统计分析方法建立了位移响应的概率谱密度函数（PDF）并绘图,清楚地显示了后屈曲板的PDF表现出双峰现象。使用功率谱密度（PSD）函数分析了响应频率和峰值随着温度升高的变化,并确定了板的软化和硬化区域。总结了屈曲前/后结构特定区域拉应力和压应力的变化规律,并阐述了造成这种变化的原因。本文工作可对热声载荷作用下薄壁结构响应分析和动强度设计提供参考依据。     

有局部稀薄气体效应的高超声速流动数值模拟

近空间高超声速飞行器当飞行高度和速度足够高时,其流场计算可能要考虑稀薄气体效应,传统的计算流体力学(CFD)方法预测的阻力和升阻比将不够准确。而现有的模拟稀薄气体流动的计算方法由于其计算量巨大,难以在工程实际中应用。因此需要发展能用于近空间高超声速飞行器流场的可行、可靠的计算方法。陈杰和赵磊在文献[1]中针对边界层中既有强剪切而气体分子自由程又相对较大的情况进行分析,提出了刻画此类局部稀薄效应的无量纲参数Zh,并提出了在传统CFD中通过采用依赖于Zh参数的等效黏性系数考虑局部稀薄效应对阻力计算影响的研究思路。因此,本文尝试将此等效黏性系数纳入CFD模型中,以在70km高空,以马赫数15飞行的小迎角钝平板为例,来检验计算方法是否合理可行。结果表明:和传统的CFD方法所得结果相比,新模型计算的阻力减小,升阻比增加,其改进的方向与现有飞行试验结果定性相符,且所增加的计算时间非常有限,可方便地应用于现有的计算空气动力学中。     

基于ABAQUS二次开发变角度层合板屈曲特性分析

基于ABAQUS二次开发,探究周期函数曲线纤维路径变角度复合材料层合板的屈曲特性。首先,以正弦曲线为基本参考路径为例,得到纤维角度的变化规律;然后,利用Python编写变角度复合材料层合板的有限元分析前处理程序,开发出ABAQUS交互式界面;最后,利用本文开发的GUI插件,对经典定角度和变角度层合板进行屈曲分析。结果表明,变角度层合板的屈曲载荷有很大提高,并且随着幅值参数A和周期参数T逐渐增大,一阶线性屈曲均是呈先增大后减小的趋势,因此基于铺层稳定性,参数应控制在相应的范围。采用Python对ABAQUS进行二次开发,从而实现变角度层合板的自动建模和计算分析,为实际工程研究提供了研究思路和流程,有一定的实践意义。     

局部规则涂层对气冷叶片测温影响的数值研究

为研究微细热电偶测温安装涂层对其测量精度的影响机理,基于气固耦合传热(CHT)方法,建立了接近真实气膜冷却叶片的三维模型,选择SST(shear stress transport) γ -θ湍流模型,对未喷涂氧化铝涂层与在叶片吸力面前缘、压力面前缘、压力面尾缘处喷涂氧化铝涂层后的气冷叶片表面温度进行数值研究。无涂层工况的数值计算结果与实验结果误差控制在5%以下,验证了数值方法的有效性。研究结果表明:涂层对近壁面燃气的流动特性产生显著影响;三种有涂层工况的数值计算结果表明,压力面前缘处比吸力面前缘处温差波幅小,且比压力面尾缘处测温误差降低35.5%;在压力面及吸力面前缘处,最佳测温区域为涂层中部至下端之间的部位,而在压力面尾缘处,最佳测温区域尽量选择靠近涂层上端部位。      

低雷诺数翼型局部振动非定常气动特性

针对低雷诺数翼型特殊的气动特性,采用基于动网格的非定常数值模拟方法,研究翼型表面不同弦向位置的局部蒙皮以不同频率及振幅振动时对低雷诺数翼型气动特性及流场结构的影响,揭示蒙皮振动增升减阻的机理。研究表明,在低雷诺数条件下局部蒙皮振动可有效提高翼型气动特性,与刚性翼型相比蒙皮局部振动可使翼型升力系数提高,阻力系数降低,升阻比提高。振动位置对翼型气动特性及流场结构有显著的影响,振动表面位于翼型前缘附近或位于层流分离泡中心时可有效控制翼型层流分离,从而提高翼型气动特性。振动频率对翼型表面层流分离及转捩位置均有显著的影响,随着振动频率增加,翼型气动特性出现最优值。与刚性翼型相比,表面振动使翼型转捩位置略向上游移动,摩擦阻力增加,但振动使等效翼型相对厚度减小,压差阻力明显减小。在小幅振动范围内,随着振幅增加,流场非定常特性更加显著,翼型升阻比增加。