曲线物面边界条件在DGM中的应用研究

基于二维非结构网格,采用高精度间断Galerkin有限元方法数值求解定常Euler方程。为了有效克服传统的反射物面边界条件在物面处易于生成伪熵降低计算精度的缺陷,提出一种曲线物面边界条件,从而实现对物理模型而不是数值模型进行数值模拟。对经典圆柱亚音速无粘绕流进行数值模拟,结果表明:采用曲线物面边界条件之后,流场解的精度得到很好的提高;此外,在非常稀疏的网格上,通过提高基函数的阶次仍可以得到高精度的数值解。     

高阶精度间断Galerkin有限元方法研究

在二维非结构网格上,对高阶精度间断Galerkin有限元方法求解跨音速欧拉方程进行研究。运用间断有限元理论,采用施密特正交化多项式基函数对流场解进行近似描述,使用HLLC近似黎曼解方法计算网格单元边界处的数值通量,积分项通过高斯积分求解,时间推进采用经典四步Runge-Kutta方法,并引入斜率限制器,抑制数值振荡。对NACA0012翼型跨音速无粘流动进行数值模拟,数值结果表明：间断Galerkin有限元方法具有较高的精度,较小的数值耗散和较强的激波捕捉能力。     

含孔复合材料圆柱壳冲击破坏试验与有限元分析

 为了揭示轴向压缩载荷与径向冲击载荷共同作用下复合材料壳体开孔处裂纹的产生机理,开展了含圆孔复合材料圆柱壳冲击试验,并对冲击试验进行了有限元仿真分析。提出复杂冲击载荷作用下的动态响应分析方法,运用LS-DYNA对冲击载荷作用下含圆孔复合材料圆柱壳动态响应过程进行了模拟,采用含刚度退化的Chang-Chang失效准则预测复合材料圆柱壳破坏过程,得到的冲击加速度响应曲线及破坏区域与试验结果一致,验证了本文方法的正确性。对有限元模型进行动力学及静力学破坏分析,结果表明,径向冲击引起的环向拉应力是圆孔边缘破坏区域90°铺层纤维断裂与基体开裂的主要原因,而拉应力只引起0°铺层基体开裂。由破坏起始分析可知,将复合材料圆柱壳90°铺层含量由20%提高至50%,可使结构承载能力增加56%。     

含运动边界二维多体干扰流场的数值模拟

本文采用二阶精度Godonov扩展方法求解非定常：Euler方程，以基于格心的自适应非结构网格的显式有限体积法为基础，空间和时间都是二阶精度，使用HLLC近似黎曼解的方法计算网格单元边界处的守恒量通量，对作俯仰振动的NACA0012翼型绕流及激波管中运动、静止圆柱的流场进行了数值模拟，并对数值模拟的结果进行了分析，为进一步对含多体分离流场进行数值模拟打下了基础。     

平面温度场的延拓有限元分析

 平面温度场的延拓有限元分析 大多数复杂温度场问题的解析解往往是难以获得的,必须采用一些近似的数值计算方法。延拓有限元技术合理利用了单元邻近外点信息,以构造单元内高精度广义插值近似函数,从而在不增加网格剖分单元及结点数目的条件下,有效地改善温度场的分析精度。     

无网格方法关于导数计算的程序验证

为了提高无网格方法对空间导数的数值计算能力，提出并论证了新的适用于直角坐标系、柱坐标系和球坐标系的空间导数无网格算法。基于数值解与理论解的对比分析，对数值误差和收敛性进行了后验评估。评估结果表明:该算法对函数和导数的估计均为2阶精度，支撑域尺度对数值误差的大小有一定影响，但不影响数值计算的精度等级。在离散尺度h为0.01的条件下，对所选函数及其导数的数值计算相对误差不大于0.65%。      

靶板击穿极限速度的数值模拟分析与试验验证

利用工程有限元软件MSC.Dytran,对钢球冲击LY12铝合金靶板时的击穿极限速度进行了数值模拟研究.分析了接触力比例因子、模型的有限元网格密度及靶板材料参数估计误差对模拟结果的影响.对照试验结果,发现数值模拟得出的弹体冲击靶板的现象及击穿极限速度与试验结果基本吻合.说明用MSC.Dytran进行数值模拟,能够为结构的抗弹体冲击强度分析提供指导性意见.     

网格对高超声速钝头体表面热流数值模拟结果的影响

高超声速钝头体热流的精确预估，对高超声速飞行器热防护的设计具有重要意义，而网格划分尺度对预估飞行器表面热流精度具有重要的影响。应用理论分析及数值模拟方法，提出了一种基于壁面网格雷诺数及基于钝体特征长度的来流雷诺数网格划分方法，并给出了壁面网格雷诺数的取值范围。应用所提出的网格划分方法，对不同来流条件的高超声速半圆柱及球头表面热流进行了数值模拟，结果表明，应用所提出的网格划分方法，及给出的可供参考的壁面网格雷诺数的取值范围划分网格，在满足高超声速钝体表面热流精度要求的同时，有利于网格合理分布并提高数值计算效率。     

基于网格缩小的非结构网格梯度重构及制造解精度测试与验证

网格收敛性研究是验证与确认中的重要手段,本文针对一般非结构网格难以自相似加密的问题,实现了基于网格缩小的精度测试方法。在典型的各向同性和各向异性拉伸网格上,利用网格缩小精度测试方法分别考察了梯度重构精度以及制造解流动模拟精度,将网格缩小精度测试结果与传统的网格加密精度测试结果以及理论分析结果进行对比,验证了网格缩小精度测试方法与网格加密精度测试方法及理论分析的一致性。最后将网格缩小精度测试方法应用到各向异性拉伸网格粘性制造解精度测试中,得到了预期中的结果,并初步考察了梯度重构方法、网格类型对数值模拟精度的影响,显示出网格缩小精度测试方法在一般非结构网格精度测试上的优势,具有一定应用前景。     

基于IPDG方法的嵌套网格技术

在内罚间断伽辽金（IPDG）方法框架内应用嵌套网格技术求解了二维可压缩Navier-Stokes方程。通过把黏性通量作为辅助变量实现了Navier-Stokes方程的降阶,继而用间断有限元方法进行离散得到空间半离散方程。时间推进采用Newton-Krylov隐式方法。利用库埃特流动精确解验证了该方法的精度。在此基础上,对包括有黏NACA0012翼型绕流、定常和非定常的圆柱绕流在内的若干典型测试算例开展数值模拟,进一步验证了该方法的鲁棒性和可靠性。在嵌套网格中应用了h型网格自适应技术用于提高激波分辨率,对NACA0012翼型无黏跨声速绕流开展数值模拟,结果表明：自适应之后的单元数相比无自适应时只增加了8.4%,但是激波分辨率得到了显著提高,激波附近的计算结果也与实验值符合得更好,从而验证了该方法的有效性。     

机翼和翼身融合体跨音速绕流有限元算法的研究

从全速位方程出发研究了机翼和翼身融合体跨音速绕流的有限元算法。给出适用于有限元法和复杂外形的块结构网格生成技术和使单元系数矩阵积分数目减少2/3的方法;并发展了一种内存量少的总体系数矩阵存贮方式。为了提高有限元离散方程组解的精度和加速收敛,提出α—GMRES算法代替传统全速位方程有限元解法中所习用收敛慢的线松弛迭代解法。通过机翼和翼身融合体实例考核,效果是好的,为飞行器先进气动布局的数值计算提供了一种新的有效算法。     

Navier-Stokes方程两层稳定有限元算法分析

分析了定常Navier-Stokes方程的两种两层稳定有限元算法。它们将局部Gauss积分稳定化技术和两层算法的思想充分结合,采用低次等阶有限元P1-P1或Q1-Q1对N-S方程进行数值求解。误差分析和数值算例结果表明,当粗、细网格尺度H=O（h1/2）时,它们与在细网格上的单层有限元算法具有相同的收敛速度,而两层算法却节省了大量的计算时间。相比之下,Simple算法具有更高的计算效率。而且进一步发现Oseen算法能够对小粘性系数N-S方程进行有效求解。     

飞机风挡透明件的鸟撞分析

 本文研究了风挡透明件鸟撞动力响应的分析模型和方法。考虑了撞击载荷与风挡透明件动力响应的耦合作用和非线性效应。为了能方便地使用现有的非线性有限元分析程序并避免数值解的振荡现象,在研究中不采用常用的拉格朗日法和罚函数法,而是采用动量平衡法。同时,为了避免不断调整撞击接触点和修改网格,本文采用超弹性的橡皮材料来模拟鸟体,而不是象国外那样采用流体柱。研究和比较了衍生动载荷响应解,并提出相应的看法。     

平面热应力问题的组合杂交有限元方法

利用组合变分原理可以增强杂交元方法解的稳定性,建立了平面热应力问题的基于区域分解的组合变分原理.与弹性力学问题相比,右端项改变;并且提出组合杂交元来离散新的变分原理、形成刚度矩阵;进而引入能量协调条件,不仅简化了变分原理及其相应的刚度矩阵,而且减少有限元解的误差.数值结果表明,能量协调的4节点组合杂交元数值性能最佳:实现粗网格高精度,对网格畸变不敏感且可以克服闭锁现象.     

基于IPDG方法的超声速混合层流动数值模拟

为了满足超声速混合层高精度模拟需求,实现了基于内罚方法的间断伽辽金(IPDG)方法数值模拟。通过将黏性通量作为辅助变量使得Navier-Stokes方程降阶,并利用间断伽辽金方法进行空间离散,最后采用Newton-Krylov隐式方法对空间半离散方程进行时间推进。相对于有限体积法数值精度提高到了3阶。将该方法应用于对流马赫数为0.2的二维平面超声速混合层的数值模拟中,通过与实验数据的对比验证了方法的可靠性。数值结果显示了混合层中层流到湍流的发展过程。在此基础上,将基于数值解误差分布的自适应网格技术与IPDG方法结合起来,对比发现自适应网格数量减少了9倍,计算时间减少了8倍,大幅提高了方法的计算效率。      

基于打靶试验的风扇机匣包容能力评估方法

为保障飞行安全,航空发动机机匣需具有足够的抗冲击能力以包容高速旋转状态下丢失的叶片。针对某型涡扇发动机对开式风扇机匣包容性评估需求,提出了1种结合真实机匣打靶试验和有限元分析评估机匣包容能力的方法。通过使用真实机匣和真实叶片进行打靶试验获得风扇机匣的冲击损伤情况,并基于ANSYS/LS-DYNA进行了瞬态动力学有限元分析。结果表明:采用Johnson-Cook模型预测的机匣伤形状、尺寸以及叶片残余速度均与试验结果接近,验证了数值分析方法的准确性。采用验证过的数值分析方法开展旋转状态下机匣的包容性评估,发现由于撞击姿态差异和失效模式转变,风扇机匣可以包容以100%工作转速飞出的叶片,但机匣出现长裂纹,接近临界包容状态。所提出的方法可以在不具备部件包容试验条件的情况下,以较方便的形式对机匣包容能力可靠评估。     

3-D NUMERICAL SIMULATION OF UNDERWATER EXPLOSION BUBBLE

３－ＤＮＵＭＥＲＩＣＡＬＳＩＭＵＬＡＴＩＯＮＯＦＵＮＤＥＲＷＡＴＥＲＥＸＰＬＯＳＩＯＮＢＵＢＢＬＥＬｕＣｈｕａｎｊｉｎｇ（ＳｈａｎｇｈａｉＪｉａｏＴｏｎｇＵｎｉｖｅｒｓｉｔｙ，Ｓｈａｎｇｈａｉ，Ｃｈｉｎａ，２０００３０）Ａｂｓｔｒａｃｔ：Ｔｈｅｐｒｏ...     

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.     

Convergence Analysis of the Numerical Solution for Cathode Design of Aero-engine Blades in Electrochemical Machining

As a main difficult problem encountered in electrochemical machining （ECM）, the cathode design is tackled, at present, with various numerical analysis methods such as finite difference, finite element and boundary element methods. Among them, the finite element method presents more flexibility to deal with the irregularly shaped workpieces. However, it is very difficult to ensure the convergence of finite element numerical approach. This paper proposes an accurate model and a finite element numerical approach of cathode design based on the potential distribution in inter-electrode gap. In order to ensure the convergence of finite element numerical approach and increase the accuracy in cathode design, the cathode shape should be iterated to eliminate the design errors in computational process. Several experiments are conducted to verify the machining accuracy of the designed cathode. The experimental results have proven perfect convergence and good computing accuracy of the proposed finite element numerical approach by the high surface quality and dimensional accuracy of the machined blades.