基于Euler／N—S方程的跨音速非线性静气动弹性问题研究

在C-H网格的基础上,采用Jameson的中心差分有限体积法求解Euler/N-S方程,采用结构影响系数法计算结构的弹性变形,用三角元面积加权法和常体积转换法(CVT)实现流固耦合,计算了弹性后掠机翼在跨音速状态的静气动弹性问题     

平移速度的变化对搅拌焊接材料流动的影响

基于有限元方法对不同搅拌头平移速度下搅拌摩擦焊接过程中焊接构件Al6061-T6材料的三维流动问题进行分析,研究了搅拌头平移速度的变化对等效塑性应变的分布和材料流动的影响。结果表明,搅拌头前方的材料向上涌起,被旋推到搅拌头后方并开始向下运动,这一过程是促使搅拌摩擦焊接顺利完成的主要原因。材料的等效塑性应变等值线与不同焊接区域的边界具有较好的对应关系。随着搅拌头平移速度的增加,垂直于焊缝方向上的截面上,等效塑性应变减小,但是材料的流动速度会随之增加。较低的搅拌头平移速度有利于焊缝中心线两侧材料的融合。     

求解Euler方程的无网格与有限体积混合算法研究

研究了无网格方法与有限体积法相结合的、求解Euler方程的混合算法.与单一的无网格方法相比较,由于在大部分计算区域采用了有限体积法,使得发展的算法在运算效率上能与有限体积法相当;同时在物体附近嵌入了无网格区,使得在处理几何外形上更具灵活性.有关算法涉及的无网格与网格区域搭接处理,通过在交界面处引入辅助卫星点和网格单元构成边界信息,实现了区域间的流动信息传递.Euler方程的空间导数分别在两区域用有限体积法和无网格法离散近似,时间方向都采用四步显式Runge-Kutta格式推进求解,数值模拟了喷管内流和绕翼型外流,并分别与整体有限体积法和整体无网格方法进行了比较.算例展示出,提出的混合算法能有效捕捉激波间断,且两区域等值线过渡光滑,算法效率如预期与有限体积法相当,表明该方法是可行的.     

钛合金焊接件疲劳寿命估算

基于有限单元法理论，采用ANSYS有限元分析软件，对航空发动机中常用钛合金焊接试件进行应力应变分析，得到了焊接接头的等效应力场和应变场，并确定出试件等效应力集中系数Kf，运用局部应力-应变法对光滑焊接试件估算疲劳寿命，并采用成组试验法对计算结果进行验证和分析。     

三维四向编织复合材料基本性能的有限元模拟

依据三维四向编织复合材料的结构特点,建立了能真实反映其空间构型的"双纽线"单元体胞模型,并以此为基础,采用有限元方法对C／SiC三维四向编织复合材料的基本性能参数进行了数值模拟,获得了该材料在制备冷却阶段的热变形行为以及残余应力和残余应变的分布状况,同时获得了制备后室温条件下该材料的等效弹性模量与等效热膨胀系数值。     

残余应力对加工变形影响的分析与模拟

为了研究材料去除过程中残余应力的重新分布及该过程所引起的工件变形,假设在该过程中没有引入另外的应力,采用理论计算和有限元模拟的方法对含有已知单向应力情况的矩形板进行了求解,计算了在剥层过程中的曲率、挠度和残余应力,并分析了剩余部分的残余应力变化趋势。结果表明理论解与有限元解获得了极好的吻合。采用有限元法进一步研究了材料去除过程中,单向和双向残余应力对加工变形的影响以及最终残余应力分布规律的异同。     

304不锈钢小冲孔蠕变试验与损伤数值模拟研究

对304不锈钢进行了恒定温度下多种载荷的小冲孔蠕变试验,在此基础上建立了304不锈钢小冲孔蠕变试样的有限元模型。应用改进的K-R蠕变损伤本构方程,分析了在650℃、恒载荷条件下的试样中心挠度、应变等随时间的变化规律以及试样蠕变损伤与位置的关系。结果表明,试验与数值模拟结果基本一致,试样中心蠕变曲线都具有明显的三个阶段,与单轴试验的应变曲线十分相似。整个试样在中心损伤比较严重,呈现明显局部化特征。试样的失效首先在下表面,距试样中心约1/10处,这与试验观察到的试样颈缩部位完全一致。     

弹塑性梯度材料热机耦合分析的有限元方法

陶瓷/金属功能梯度材料(FGMs)的组份及性能梯度分布,在航空航天高温热结构中有重要应用。材料性能在空间上的连续变化给这类结构强度设计的有限元分析造成了一定困难,针对上述问题,提出一种热弹塑性梯度有限元方法。该方法通过编写有限元软件ABAQUS的材料用户子程序(UMAT)来实现,在UMAT中将FGMs的弹塑性本构模型参数定义为空间坐标的函数;用热弹塑性梯度有限元方法对陶瓷-FGMs-金属典型的三明治板结构进行热机械分析。结果表明:FGMs组分布函数对其热应力分布有很大的影响,特别是在三明治结构的界面上,影响更为突出,这和现有文献的实验结果相符。     

基于库伦摩擦效应的箔片轴承极限承载力理论与试验

基于有限单元法建立了考虑库伦摩擦的波箔型径向气体箔片轴承的箔片结构模型,采用有限差分法和有限单元法耦合求解Reynolds方程和气膜厚度方程,通过求解轴颈达到极限偏心率时的轴承极限承载力,研究了箔片结构库伦摩擦效应对轴承极限承载力的影响规律,并搭建了轴承极限承载力测试试验台,利用温度法测量了两个具有不同轴承壳内表面粗糙度的波箔型径向气体箔片轴承的轴承极限承载力.通过对比分析仿真结果与试验结果表明:轴承壳圆柱孔内表面粗糙度为0.4μm的轴承在10000r/min和20000r/min下,轴承极限承载力分别为15.5N和42.3N;而表面粗糙度为1.6μm的轴承极限承载力为10.9N和29.6N,这是由于波纹箔片和轴承壳体之间的库伦摩擦力增大了波纹箔片的刚度,因此增大箔片结构摩擦因数使得轴承极限承载力降低,并且仿真结果变化趋势和试验结果变化趋势吻合.      

基于概率的航空发动机飞行换算率改进算法

针对斯贝MK202发动机应力标准对零部件工作寿命评定中未包含各种随机因素的影响,考虑零部件疲劳性能参数的分散性和发动机整机载荷谱的随机性,将实测数据统计分析、零部件应力有限元分析以及蒙特卡罗模拟分析有机结合,提出了基于概率的航空发动机飞行换算率改进算法,以解释航空发动机零部件实际工作寿命存在的差异.经范例验证表明:确定性寿命评定方法的飞行换算率计算结果与改进算法下的平均值相当;对于单次飞行换算率,改进算法的单次飞行换算率平均值比确定性方法的计算结果小2%;对于综合飞行换算率,改进算法的综合飞行换算率平均值比确定性方法的结果小1%;且在改进算法下单次飞行换算率和综合飞行换算率服从正态分布.