高超声速流场与结构温度场一体化计算方法

在对国内外流场与结构温度场一体化计算方法的不足进行细致讨论的基础上,提出了一种高超声速流场与结构温度场一体化计算方法.采用统一的积分方程组作为气动加热和结构传热物理过程的控制方程,对整个物理场进行统一的迎风格式有限体积方法离散,给出了流场与结构交界面上温度、温度梯度及导热系数等参数的计算方法.在时间推进方面,定常状态采用多步龙格库塔迭代格式,非定常状态则采用双时间步长方法.采用发展的一体化算法对二维圆管模型的气动加热和结构传热问题进行了数值模拟.结果表明:2s时驻点物面温度为390K,与其他文献的误差在3.1K范围内;稳定时驻点物面温度为647K.      

全复合材料主承力梁结构研制与分析验证

根据结构特点和受力特性,研制了一种全复合材料主承力梁结构。给出了主梁的铺层设计,详述了其成形工艺。基于MSC.Patran/Nastran建立了主梁有限元分析模型,采用线性求解器对多级载荷作用下的主梁进行了静强度分析,给出了结构应变与位移结果。同时设计了强度试验方案,对主梁结构进行了静强度试验验证,采集了应变与位移信息,确定了结构的响应。结果表明,全复合材料主承力梁结构满足强度设计要求,为复合材料应用于主承力结构提供了依据。     

复合材料机身壁板屈曲优化设计

针对受压加筋曲壳的后屈曲特点,提出了一个后屈曲载荷的快速计算方法。用一部分有效蒙皮代替蒙皮失稳后的承载能力,利用欧拉失稳模型估算后屈曲失稳临界载荷。以外形参数给定的机身加筋壳结构为例,基于蒙皮局部失稳的简化模型,利用二次序列法优化求解结构的最佳布局。然后按照之前提出的后屈曲载荷的快速计算方法,设计筋条的截面刚度,最终得到最轻的机身结构质量。     

基于直通型CT试样测试方法的断裂韧度研究

依据ASTM用于CT试样获取断裂韧度载荷分离理论规则化方法提出量纲一载荷分离法,从而完善用于断裂韧度测试的载荷分离理论,且根据弹塑性有限元精细计算结果以及直通型CT变形几何关系,提出裂纹嘴张开位移与加载线位移的转换公式。采用不同初始裂纹长度CT钝裂纹试样对公式的有效性和精确性进行试验验证,进而实现对1Cr12Mo、30Cr2Ni4Mo V以及P91管材10Cr9Mo1VNb N三种合金直通型CT试样的JR阻力曲线测定,结果表明,试验值与公式计算值的最大相对误差均在5%以内。     

基于Gambit前处理的气膜冷却火焰筒壁温分析

为更好实现航空发动机燃烧室气膜冷却结构火焰筒的壁温分析和冷却结构优化,针对原有的火焰筒2维壁温计算程序开展了2次开发工作,形成了基于GAMBIT前处理的火焰筒壁温分析程序。新的计算过程采用Gambit软件对几何模型进行前处理,生成三角形网格、指定边界条件分组,通过编制前处理模块代码,对导出的网格文件进行解析,进一步将各种信息导入已有的有限元壁温计算程序,完成壁温计算分析。给出了采用以上方法进行火焰筒壁温计算和优化的实例,结果表明:该方法较为有效地克服了原方法的各项缺点,其交互性强,大幅度提高了工作效率。     

基于大变形蠕变分析的持久寿命预测方法

提出一种基于真应力真应变弹塑性蠕变本构模型和大变形有限元分析的高温构件持久寿命预测方法.该方法利用以真应力-真应变表示的材料高温拉伸应力-应变曲线建立材料的弹塑性模型,基于蠕变曲线建立蠕变本构模型,并采用大变形有限元方法计算高温构件在给定载荷下的变形响应曲线,根据其响应曲线的变化趋势来确定构件持久寿命.通过TC11钛合金缺口试件500℃下的持久试验对上述方法进行验证,并与三种基于小变形分析的持久寿命预测方法进行对比.结果表明:本工作提出的方法可以较准确地预测TC11缺口试件的高温蠕变响应和持久寿命,其预测精度优于基于关键点断裂应变、缺口净截面平均有效应力以及骨点应力的小变形有限元分析的寿命预测方法.     

ZL205A铝合金大型复杂筒体构件淬火过程模拟与变形预报

基于ABAQUS软件采用有限元方法建立ZL205A铝合金大型复杂构件淬火过程的温度场与应力场模型,分析构件在淬火过程中不同时刻的温度场分布.根据构件结构特征,选择4个特征点获得淬火过程冷却曲线,并对构件厚壁、加强筋与薄壁部分的冷却曲线进行研究,使用有限元模拟方法预报构件变形和残余应力.ZL205A铝合金大型复杂构件淬火变形主要为径向变形,通过比较淬火后残余应力的6个分量可知正应力远大于剪应力.模拟厚壁在不同换热系数下的淬火效果,发现增加换热系数可以提高厚壁淬火冷却速度,但是作用比较有限.     

异形型腔电解加工阴极的数字化设计

建立了异形型腔电解加工阴极设计和加工过程仿真的数学物理模型,并在MATLAB平台上利用PDETool工具箱编制了有限元求解和动态加工仿真程序,针对某型发动机扩压器叶间流道异形型腔设计制造了阴极和工装夹具,选择并优化了加工参数,试验加工出的型腔试件能均匀留出0.15～0.30 mm的加工余量,满足了实际加工的工序要求.     

高阶辐射换热有限元算法研究

随着航空航天技术的发展,飞行器热结构所需承受的温度越来越高,辐射换热变得非常重要;当前大部分商业软件对于辐射问题的有限元计算方法还是基于一致表面温度和辐射热流假设,使得计算精度和网格密度的矛盾越来越严重.进行了选用高阶单元、采用高斯积分精确计算单元表面变辐射热流方法的研究,从而摆脱了一致表面温度和辐射热流的假设,使得在相同网格密度的情况下计算精度大大提高;同时,从包含辐射换热问题的有限元计算方程出发,采用与有限元数值计算时相同的积分方案,只在独立的积分点处计算辐射热流,克服了积分方法计算效率低的缺点.经与ANSYS的计算结果对比,应用辐射热流积分方法于高阶单元能大大提高计算精度;并且在相同计算精度条件下,此方法的计算效率更高,具有一定的实用价值.     

栅格冀导弹流场混合网格N-S方程数值计算

采用结构和非结构混合网格技术，对栅格翼导弹黏性流场进行数值模拟，预测其气动特性。计算的马赫数为0．7～2．5。用有限体积法和LU-SGS算法求解N-S方程。计算与实验数据以及栅格翼附近流场结构进行了比较。