航空发动机主轴承接触应力精确仿真计算方法

采用稳态热分析与动态接触力学分析相结合的方法以提高计算精度,考虑轴承摩擦生热及瞬态冲击效应的影响,得到经热应力修正的轴承应力分布及变化特征。研究表明:最大动态接触应力位于钢球表面,外圈滚道次之,内圈滚道最小;轴承接触应力具有明显动态冲击响应特性,分布及大小具有随机变化性;接触应力由主应力与剪切应力共同构成,随转速及轴向载荷增大而增大。      

流线曲率法在多级跨声速轴流压气机特性预测中的应用

为研究多级跨声速压气机的分析问题,以通流理论为基础,采用了一系列适用于跨声速压气机的攻角、落后角和损失等经验模型,发展了一套基于流线曲率法的通流计算程序来预测跨声速压气机流场及其工作特性。为提高经验模型的预测精度,考虑到真实压气机中复杂的三维流动效应,针对部分早期模型进行了合理改进,包括改进了落后角模型使其适用于更大弯度范围叶型,以及采用一种更为合理的可变结构激波损失预测模型。针对两台跨声速压气机算例进行了计算校验,并将校验结果与实验值和三维数值计算进行对比。对比表明,设计工况下总压比最大计算误差为4.1%,效率误差为1.1%,在非设计工况特性预测和展向流场参数计算中也能得到和实验值相符的变化趋势,该通流计算方法可为现代跨声速轴流多级压气机特性分析提供具有参考价值的预测结果。     

利用渗透边界模型分析三维内转式进气道启动性能

进气道抽吸区域一般包含大量抽吸孔,这些抽吸孔的网格前处理异常繁复、离散求解及CFD仿真困难。为避免这些问题,利用集成了渗透边界模型的数值仿真软件AHL3D模拟小孔抽吸,获得了小孔抽吸对三维内转式进气道Ma4~6内的启动性能的影响。结果表明:同等条件下,渗透边界与抽吸孔仿真的机体侧壁面压力曲线基本重合,且进气道喉部参数最大差别小于1.5%,说明利用渗透边界模型研究抽吸对进气道启动性能的影响具有可行性;边界层抽吸位于分离泡最高压力点附近时,可实现进气道宽马赫数范围(Ma4~5.5)的启动;Ma5条件下,开孔率在0.1左右,进气道实现启动,且启动后流量抽吸率低于1%;抽吸背压为6.5倍来流静压时,进气道实现启动,启动后流量损失几乎为0,压力分布规律与远场初始化得到的启动流场完全一致。     

预混气爆轰马赫反射实验研究

为研究斜爆轰马赫反射的实验现象以及不稳定性对马赫反射的影响,在矩形爆轰管道中,对3种预混气的马赫反射现象进行了实验研究,取得烟膜记录。烟膜结果显示,三波点高度在楔角附近与无反应三激波理论类似,随后逐渐偏向反应三激波理论,三波点实际高度由反应实际放热量决定,三种气体偏转始末距离与结束偏转距离之比平均值分别为0.42,0.49,0.31,与三者活化能大小关系相符合。胞格结构随初始压力增大而渐趋规律,胞格尺寸也逐渐减小,由CJ区域到过驱区域胞格结构存在明显的转变。横波是影响斜爆轰稳定的重要因素,压力增大可抑制分子的无规则运动,使横波出现的位置更加集中,胞格更为规律。预混气爆轰不稳定度大的气体横波出现的时空位置更不易确定,胞格结构更不规律,胞格尺寸较小。     

一种TBCC 可调喷管流动特性的实验研究

针对一种喉道面积与出口面积独立可调的串联布局TBCC喷管,进行了内外部流动特性的实验研究,并将所得到的结果与仿真结果进行了对比。结果表明,该可调喷管能够在宽广的落压比和通流流量范围内正常工作,在起飞状态、过渡状态、巡航状态下,其流动结构建立正常,推力系数较高,分别为0.90,0.96和0.97。并且在各典型状态下,实验、仿真获得的流动结构、沿程静压分布曲线均具有较好的一致性。     

基体特性对HTPB推进剂力学性能的影响

按工程生产工艺要求,制得hydroxyl terminated polybutadiene (HTPB)固化胶片,通过控制试验环境温度（-50、-35、-20、0、20、35℃）,得到了基于Prony级数形式的具有不同松弛特性的复合基体材料的松弛行为描述,为细观模型的建立提供了必要参数。随后基于分子动力学方法,结合具有内聚本构的黏接单元,建立了推进剂的细观计算模型。为验证所建数值模型的准确性,对固体推进剂在常温条件下进行宏观松弛力学试验,将仿真结果与试验数值的对比,试验结果与仿真结果相差在20%范围之内。最后对具有不同随机分布及不同基体松弛特性的推进剂细观模型进行有限元计算,结果表明:颗粒随机性不影响推进剂的宏观力学行为,而基体松弛特性显著影响固体推进剂的宏观力学性能,基体松弛特性获取环境温度与推进剂的宏观初始模量、延伸率以及强度均呈指数型关系。      

复合材料机身C型柱准静态压溃仿真及失效模式

为了研究复合材料机身薄壁C型柱结构轴向压溃吸能特性、失效模式及C型柱多层壳单元建模方法,建立多层壳单元有限元模型,基于准静态轴向压溃实验结果进行对比验证。结果表明:C型柱多层壳单元模型能够在一定程度上模拟层间分层失效及压溃过程中的局部弯曲变形和层束弯曲失效模式;仿真与实验的载荷-位移曲线吻合性较好,压溃初始峰值载荷,比吸能以及压溃均值偏差较小;但C型柱结构压溃初始峰值载荷较大,载荷效率较低,需通过优化设计进一步降低其初始载荷峰值。     

AZ31B镁合金双道次热变形过程流变应力特征模拟研究

建立了一类镁合金双道次热变形过程的二维元胞自动机模型。考虑到镁合金具有六重对称的密排六方晶体结构,模型采用了六边形元胞网格。另外,通过计算位错密度演变和设定晶粒形核长大规则,模型将热变形过程中动态再结晶、静态回复、静态再结晶、亚动态再结晶和晶粒长大等单个物理过程耦合到一起。为了验证模型,将其应用到AZ31B镁合金单道次和双道次热压缩过程。探讨了不同变形温度、应变速率、道次间隔和预应变对应力-应变曲线的影响。将计算结果和实验结果进行了对比,吻合较好。     

基于结构动力学特性的损伤识别方法

基于等效模态应变/动能理论,提出了一种利用实际结构的测试数据识别结构中损伤位置的方法。在此基础上,研究了利用模型修正技术识别结构中损伤强度的方法。分别以一个单损伤平板结构和多损伤平板结构为例,通过仿真分析了以上方法的有效性。结果表明,以上方法可以有效识别结构中的损伤位置和损伤强度。     

高雷诺数槽道湍流的壁面模化大涡模拟研究

选取基准壁湍流的槽道流动,研究了多种模型的壁面模化大涡模拟.模型包括经典的大涡模拟、Spalart-Allmaras、分离涡模拟和一种动态混合模型.基于摩擦速度的雷诺数范围为395~12000,采用3组粗糙网格,流向和展向维数分别同取37,49和65,法向维数保证y+(1)~1.主要研究平均速度、雷诺切应力分布、详细分析了各模型的特性差异并展示了相应的湍流结构.研究表明:在高雷诺数粗糙网格下,大涡模拟失去求解精度,分离涡模拟出现对数律不匹配,动态混合模型的计算接近直接数值模拟,其对数率区可解应力约占雷诺切应力的93%,边界层外层可解应力约占99%.这说明合适的混合模型可以在经济成本下保证计算精度,具有解决实际问题的潜力.