航发高温内腔温度场模型的建立与仿真分析

航空发动机内腔温度高、变化快，传统接触测试方法需要破坏结构或加载传感器于内腔获取温度数据。针对这种高温测量应用需求，从热传导和热弹性理论出发，理论分析了温度场作用下的钢质薄板变形特性。基于热弹性的基本方程和边界条件，研究了符合轴对称原则的金属薄板的应力、应变、位移、温度之间的数学关系，建立了多项式方程求解金属薄板温度位移模型。依据金属薄板材料参数和边界约束条件，有限元仿真分析了金属薄板的热变形状态，仿真结果与建立的模型计算数据比对验证基本一致，初步验证了温度计算模型的可靠性，为微变形内推温度测量提供可能。     

边界元法在薄板稳定问题中的应用

 <正> 1.基本公式 弹性薄板的稳定方程为式中D是板的弯曲刚度;λ是临界载荷因子;为双Laplace算子;N_1,N_2,N_3为板的中而力。 首先确定基本解。即:在无穷域上某一点作用一个单位力时,双调和微分算子的解,即板的挠度。与板的挠度相对应的基本解对应于方程     

集中力作用下悬臂梁几何中轴的弹性大挠度分析

以平截面假定为基础,并以几何中轴作为梁变形的参考轴,推导了大挠度情况下,弹性梁在轴力和弯矩联合作用下几何中轴的曲率——弯矩和轴力方程,并建立了求解几何中轴变形曲线的方程组,这样的方程组考虑了梁轴向变形的影响,克服了以中性轴作为梁变形参考轴时不能考虑轴向变形的不足。利用所建立的方程组,分析了受集中力作用的弹性悬臂梁大挠度问题,并同以中性轴作为梁变形参考轴的部分结果进行了对比。结果表明,当载荷较小时,2种理论计算出的弹性大挠度悬臂梁的变形曲线差别很小;而当载荷较大时,2种理论计算出的弹性大挠度悬臂梁的变形曲线存在比较明显的差异。     

航空薄壁板结构高温声疲劳分析

采用解析方法,分析了薄壁板结构在热载荷与噪声联合作用下的响应特性;依据薄板的大挠度非线性理论,采用以基频响应为主的单自由度模型,推导出了具有温度影响的非线性运动方程,并利用等价线性化方法将该方程线性化,在此基础上得出了应力响应的统计特性。     

弹性薄板弯曲的有限元分析法

弹性薄板弯曲是我们在制造过程中遇到的最常见的问题,为使问题简单化,便于理解,我们用有限元法分析薄析问题,用薄板的单元组成来取代原来的薄板,通常采用本角形和矩形薄板单元的组合代替壳体,由平面应力和弯曲应力状态加以组合得到薄板的应力状态。采用三角形单元,并将每一单元作为相应的平面应力单元和薄板弯曲单元的线性组合,介绍此类薄板弯曲问题的有限元法。     

Large Amplitude Flexural Vibration of the Orthotropic Composite Plate Embedded with Shape Memory Alloy Fibers

研究了形状记忆合金(SMA)纤维混杂复合材料大挠度层合板的非线性自由与受迫振动特性。基于描述SMA力学行为的Brinson理论以及层合板材料性能预测的混合率,建立了SMA纤维混杂复合材料大挠度层合板的本构方程,基于对称层合各向异性弹性板的非线性理论,建立了以横向挠度和应力函数表示的板的横向振动方程和相容方程。采用Galerkin近似解法将振动方程化为时间变量的含有3次非线性项的Duffing型常微分方程,采用谐波平衡法(HBM)获得系统的固有频率方程和强迫振动稳态频率响应方程。数值计算表明,非线性板自由振动频率比与激励温度的关系具有与线性板相同的特征,马氏体相向奥氏体相转变阶段温度对板的振动频响特性曲线的影响最显著,同时也讨论了SMA纤维含量、板的纵横比以及自由振动幅值对板的非线性频率比的影响。     

等离子电弧加热成形过程中薄板宽度的影响规律

当板厚、板长一定时采用不同宽度的1Cr18Ni9Ti薄板进行等离子电弧弯曲成形的试验,试验发现:随着薄板宽度的增长,薄板弯曲角度增大,当板宽的增大到一定值时,薄板弯曲角度趋于定值。随后在实验基础上,对等离子电弧加热过程中热影响区的弹塑性变形进行合理的简化,建立了薄板弯曲角度随薄板宽度变化的受力模型,该模型的计算结果与实验结果相符。     

含孔边裂纹板的弯曲断裂计算

 采用复变函数理论和边界配置方法 ,对含孔边裂纹板的弯曲断裂进行了分析计算。首先假设挠度的复变函数式 ,进而可以求板的内力。它们能满足一系列的基本方程和支配条件 ,仅板的边界条件需要考虑 ,并且可用边界配置法和最小二乘法近似满足。对孔双边裂纹问题进行了应力强度因子计算。数值算例表明 ,本文方法精度较高 ,计算量小 ,是一种有效的半解析、半数值计算方法。     

夹层板大挠度问题的力学表征与数值模拟

随着材料制备工艺的发展和成型加工技术的日趋完善,夹层板在现代航空领域结构设计中获得了广泛应用,以获得理想的重量及其它功能性指标,工程分析时往往利用薄板理论进行近似分析,这忽略了夹层板的良好的抗弯特性,给分析结果带来较大误差.在薄板大挠度问题及夹层板相关理论的基础上,提出分段叠加法分析大挠度夹层板的变形机理,并利用有限元进行了数值模拟,结合某型飞机蜂窝轮舱壁板静力试验验证了夹层板大挠度变形机理的合理性,证明了叠加法的有效性.     

复合材料斜交铺层矩形板弯曲的简化计算

复合材料斜交铺层层合板应用广泛,尤其是对称斜交铺层,制造工艺较简单。同时,若采用不同的铺层角,可得预期的两个垂直方向不同的弯曲刚度比。关于在小挠度下,薄板的弯曲问题已有多种解法,但一般说,计算较烦,计算量大。这里提出一种简化方法,计算简便,适于工程应用。试验表     

强制偏流喷管的数值仿真

采用薄层Ｎ－Ｓ方程成功地求解了强制流喷管内的三维复杂流场，对通量项和压力项用Ｖａｎ Ｌｅｅｒ发展的矢通量分裂法进行了离散，对离散方程采用线超松弛方法迭代求解，在亚音速与超音速气流共同存在的出口边界分情况进行了处理。计算结果验证了强制偏喷管具有自动补偿外界压力变化的能力，与实验结果基本吻合。     

面内弯曲载荷作用下两边简支两边固支复合材料层合板的屈曲

为扩充航空复合材料结构稳定性设计手段,基于Levy形式级数和有限差分法,提出了面内纯弯曲载荷作用下两边简支两边固支（SSCC）复合材料层合板临界屈曲载荷的数值解法。首先,应用有限差分法将挠度函数进行离散化处理;然后利用边界条件扩充方程组,最后将临界屈曲载荷的求解转化为广义特征值问题。通过数值算例分别在各向同性板和复合材料层合板上验证了本文方法,并与有限元方法（FEM）进行了比较。结果表明,本文方法具有很高的精度,为复合材料层合板的屈曲分析开辟了新的途径。     

梁/液体耦合振动的边界积分对称列式法

 利用格林积分导出的液体内部扰动压力函数作为测试函数,通过变分得到一组具有对称系数的代数方程,将该方程同梁的有限元振动方程联立后,即可求解梁在液体中的耦合振动问题。本方法具有未知数目少,数据存储和计算效率高的特点。     

N-S方程数值模拟方法研究及边界层中相干结构的非线性演化

建立了一套三维不可压流动的直接数值模拟方法。该算法在x及y向构造了基于非等间距网格的紧致有限差分格式和非线性项的迎风紧致型格式,在z向采用Fourier谱方法。并研究、提出了三维、非定常流体运动下游边界的无反射出流条件。另外,本文还构建了多涡结构的初始扰动理论模型,利用直接数值模拟的方法研究了边界层中多涡结构的非线性演化特点。结果表明多涡结构非线性演化中的许多特点,如雷诺应力分布、高剪切层的形成、喷射和扫掠等与壁湍流中相干结构的发展规律和现象十分相似。同时,通过对边界层中多涡结构的非线性演化问题的直接数值模拟,验证了该算法具有计算精度高,稳定性好,收敛速度快,出流边界影响小等优点。     

一种基于有限元法和弹性接触理论的齿轮啮合刚度改进算法

提出了一种将有限元法和弹性接触理论相结合的齿轮啮合刚度计算方法.该方法利用子结构法提取齿面原始柔度矩阵并分离出接触点弯曲变形,根据线弹接触变形解析公式计算接触变形,通过求解非线性变形协调方程得到齿轮时变啮合刚度和齿面载荷分布.以一对齿轮副为例,计算的啮合刚度与航空标准计算结果相差在6%以内.该方法发挥了有限元法在预测物体整体变形方面的优势,同时结合弹性接触理论能够准确计算局部接触变形的优点,与常规有限元法相比,能够有效地提升计算效率.由于接触变形问题的非线性,啮合刚度随总啮合力增加呈现非线性增大的趋势.      

直升机桨叶扬起下坠碰撞动响应计算

针对铰接式旋翼桨叶扬起下坠碰撞问题,采用中等变形梁模型,同时引入有限转角模拟桨叶绕铰的刚性运动。运用Hamilton原理建立桨叶扬起下坠碰撞动力学方程,采用隐式Newmark法计算桨叶在有初始挥舞角下的下坠碰撞响应。计算结果与试验数据的对比表明:如果采用定步长的方法计算响应,计算时间较长,计算结果会产生明显的延迟。为了解决上述问题,在碰撞前后设置了一碰撞缓冲区,采用了变步长的方法强制计算点在跳过此缓冲区时必须进入该区。此法在计算桨叶扬起下坠碰撞动响应时,能较好的解决计算延迟问题,计算时间大大减少,计算结果与试验数据吻合较好。     

Nonlinear deforming of laminated composite shells of revolution under finite deflections and normal’s rotation angles

The discretization of the boundary value problem for laminated composite shells is based on the finite difference approach using the regular mesh with the constant grid step and the difference operators of the second order of accuracy. The dynamic relaxation method is proposed for the solution of the nonlinear problem. The evolutionary equations of the dynamic relaxation are constructed, and the optimum parameters of the converging linear iterative process are estimated.     

非线性壁板颤振计算的子循环预测校正方法研究

为了解决流固耦合问题中弱耦合计算方法存在的滞后误差引起计算发散的问题, 将子循环方法与预测校正方法相结合, 使用有限元方法求解壁板动力响应, 使用Fluent求解受扰动的流场, 使用基于拉伸弹簧与扭转弹簧的动态网格方法更新流场网格, 分析了弹性壁板在超声速气流作用下的非线性动力响应, 结果表明此方法能准确预测壁板颤振的响应特性, 计算稳定性优于弱耦合方法, 可以保持长时间计算结果的稳定性.