某风洞主体结构的有限元分析

由于某风洞主体结构为典型的薄壁结构，所以在设计阶段必须考虑其静强度、静刚度、动态特性以及支座形式和布置问题。笔者阐述了风洞主体结构的网格划分方法，根据模型简化原则，通过MSC／Patran建立了有限元计算模型；同时，通过CAE软件MSC／Nastran，根据先进的模态计算方法-Lanczos法，对该风洞主体结构进行了有限元分析，获得了静力学以及模态分析结果。结果表明：该风洞主体结构满足静强度与静刚度要求，总体刚度分布合理，模态频率分布良好，但应对部分部件进行进一步优化。     

三维复杂外形非结构网格的生成及应用

以 Delaunay方法为基础生成了三维复杂外形的四面体非结构网格。这种方法的特点是先布结点 ,然后按 Delaunay准则将其连接起来。与其他方法相比 ,其好处是能较好地根据计算需要控制网格点的疏密 ,提高计算效率 ,其运行过程也较稳定。在布点过程中采用了背景网格 ,使点的生成自动而合理。文末给出了 M6机翼和简化航天飞机的数值算例。     

三维非结构混合网格区域分解及数据处理

讨论了三维非结构混合网格区域分解技术及子分区显示技术。在三维复杂外形的流场计算中,为了得到较高的计算效率及保证数值模拟的准确性,经常要采用混合网格及并行计算技术。本文针对三维非结构混合网格的自动分区及子网格的数据处理进行了细致的讨论,编写了相应的软件。对三维实际外形非结构混合网格的处理结果表明,本文发展的方法及软件系统具有通用性强、人为干涉参数少、自动化程度高等特点,可作为通用软件推广使用。     

有限元频响分析中网格划分技术研究

阐述了利用有限元进行频响分析过程中存在的一个突出问题；简单介绍了波动理论的一些基本公式，并提出有限元频响分析时划分网格的参考标准；首先通过对简单结构(包括梁、矩形板)的有限元计算以及与解析解的对比初步验证了该方法的有效性，然后给出了两个工程算例的应用情况。     

三维非结构混合网格高超声速流场并行计算

在非结构混合网格上对三维高超声速化学非平衡粘性绕流进行了基于PC-Cluster的分布式并行数值模拟.本采用区域分裂思想,研究了三维非结构混合网格区域自动分解技术,并以此为基础对高超声速化学非平衡绕流进行了并行数值计算.控制方程为多组分N-S方程,空间离散采用有限体积格心格式,时间推进为显式Runge-Kutta格式.化学非平衡动力学模型为七组元带电离反应模型,对化学反应源项进行了点隐式处理,温度场的计算采用牛顿迭代法.在PC-Cluster上对三维双椭球模型的高超声速绕流流场进行了基于区域分解技术的并行数值模拟,所得数值结果与参考献中的结果作了对比验证。     

复合材料蜂窝夹层结构的总体稳定性研究

采用两种有限元模型和工程方法分别对复合材料蜂窝夹层结构在压缩、剪切载荷作用下的总体稳定性进行了计算,根据各计算结果与试验结果进行了对比分析。结果表明:对于承受压缩载荷的结构,采用工程计算法能够较好的预估结构的屈曲临界载荷,而对于承受剪切载荷的结构,采用三维有限元法能较好的计算结构的屈曲载荷。     

HLLC方法在三维非结构网格上的应用

以三维非结构网格的显式有限体积法为基础，采用格心方法以及基于近似黎曼解的二阶精度Godonov扩展方法求解三维Euler方程，使用HLLC(Harten，Lax，van Leer,Contact)近似黎曼解的方法计算网格单元边界处的守恒量通量。为了验证方法的可行性，利用三维解算器对NACA0012翼型绕流进行了数值计算，结果令人满意。在此基础上对某型弹的流场进行了数值模拟，取得了比较好的结果。     

三维网格动态细化技术在切削仿真中的应用

为解决切削加工有限元仿真模型求解速度与精度难以平衡的问题，提出了一种适用于三维有限元模型的网格动态细化算法。该算法的主要功能包括网格细化区域判断、网格细化以及新旧网格物理场传递。采用Python语言对Abaqus软件进行二次开发，将该算法应用于Ti2AlNb钛合金车削加工仿真之中，并最终通过实验验证了仿真模型的准确性。与采用局部网格细化的仿真模型计算结果对比，采用网格动态细化技术的仿真模型求解切削力误差增加了7.3%，求解最大应力误差增加0.2%，求解速度提升174.2%。实现了在保证仿真计算精度的基础上，有效提高运算速度。     

用加权余量法建立三维涡流场的有限元方程

本文应用加权余量法推导了三维涡流场的有限元离散化方程,使得三维涡流场的有限元方程表达方便,概念清楚,计算较为简捷。     

某发动机推力室内壁鼓起原因分析研究

针对某发动机推力室内壁鼓起烧蚀断裂问题,截取了推力室喉部区域进行三维的有限元建模分析,分别用静力分析以及动态响应分析计算了结构在压力及温度载荷作用下的塑性变形情况;同时计算了结构在高温阶段的稳态蠕变以及加速蠕变的蠕变变形。计算发现结构在喉部上游发生了塑性变形,引起内壁鼓起,但未超过其延伸率,而在加速蠕变阶段会产生很大的蠕变变形,导致了结构的断裂发生,计算预测的位置及形状和实际情况基本一致。     

反复荷载下钢筋混凝土框架节点非线性分析

论述了反复荷载下分析、计算普通钢筋混凝土及轻骨料钢筋混凝土框架节点的非线性有限元分析方法 ,采用等效单轴应变的增量正交异性混凝土本构关系模型和反复荷载下应力应变关系的“焦点模型” ,编制二维非线性有限元分析程序 ,基于四个混凝土边节点的试验结果 ,对钢筋混凝土和轻骨料钢筋混凝土框架节点进行非线性有限元对比分析 ,计算的结果与试验吻合较好。     

高速碰撞中Lagrange有限元方法及其应用

根据连续介质力学守恒方程和有限元离散方法,针对三维高速碰撞现象的特性,应用瞬时最小闰能原理,完整地推导了高速碰撞有限元基本列式和离散方程,并给出了高速碰撞问题的中心差分求解算法,讨论了非线性Ｌａｇｒａｎｇｅ动力有限元程序数值计算中应重视的问题,本文给出的典型算例说明本方法对高速碰撞问题分析的有效性。     

大挠度复合材料柔性梁的静力学分析与试验

提供了一种处理具有结构耦合的大挠度复合材料柔性梁的分析方法。在将三维弹性问题分解为二维线性剖面特性分析和一维非线性分析的基础上，采用有限元法计算梁的二维剖面刚度特性，编制了具有多个单元库的程序软件，它能有效地处理复杂剖面形状，对参考剖面点的翘曲位移未加任何限制；一维非线性分析时用Ｅｕｌｅｒ角描述梁的空间变形位置，建立的非线性微分方程组适用于任意大变形情况。文中给出了详细的分析与试验结果，并着重分析和试验了根部安装角分别为０°，４５°的两种情况，揭示了结构耦合对变形的影响，试验与分析具有良好的一致性，表明本文方法及处理简便、可行。     

Numerical Simulation of ATPS Parachute Transient Dynamics Using Fluid-Structure Interaction Method

In order to simulate and analyze the dynamic characteristics of the parachute from advanced tactical parachute system(ATPS),a nonlinear finite element algorithm and a preconditioning finite volume method are employed and developed to construct three dimensional parachute fluid-structure interaction(FSI)model.Parachute fabric material is represented by membrane-cable elements,and geometrical nonlinear algorithm is employed with wrinkling technique embedded to simulate the large deformations of parachute structure by applying the NewtonRaphson iteration method.On the other hand,the time-dependent flow surrounding parachute canopy is simulated using preconditioned lower-upper symmetric Gauss-Seidel(LU-SGS)method.The pseudo solid dynamic mesh algorithm is employed to update the flow-field mesh based on the complex and arbitrary motion of parachute canopy.Due to the large amount of computation during the FSI simulation,massage passing interface(MPI)parallel computation technique is used for all those three modules to improve the performance of the FSI code.The FSI method is tested to simulate one kind of ATPS parachutes to predict the parachute configuration and anticipate the parachute descent speeds.The comparison of results between the proposed method and those in literatures demonstrates the method to be a useful tool for parachute designers.     

嵌套网格粘性流动数值模拟用于风洞洞壁干扰研究

运用嵌套网格技术和Navier Stokes数值模拟对机翼半模和翼身组合体试验时风洞的四壁干扰进行综合模拟、评估和修正。计算格式在空间上采用中心有限体积离散 ,在时间上采用多步Runge Kutta时间步长格式进行积分。计算结果证明了该方法的可行性和优越性。     

基于响应面方法的多目标有限元模型修正技术研究

研究基于响应面方法的多目标模型修正技术,并以一个实际结构为例进行修正。论述了该技术的主要步骤及其中参数选择、试验设计、响应面拟合、多目标优化和确认准则等关键技术的理论基础。设计结构并进行相关模态试验,用仿真软件MSC.Nastran建立该结构的有限元模型,然后基于该技术对模型进行修正。比较并分析试验结果与修正后仿真结果,通过两者前四阶频率与振型的比较,证明了该技术的有效性和合理性。     

最小二乘和奇异值分解法在边界元法中的应用

采用最小二乘与奇异值分解结合的方法，给出求解系数矩阵不满秩的线性代数方程组的数值方法，进而将此方法应用于边界无法中，处理给定外力的第一边值问题，特别地用于处理给定外力的三维裂纹问题。此外，本文还给出求解三维有限体裂纹问题的超奇异积分方程组，并使用有限部积分与边界元法为其建立了数值法。最后计算了若于典型例子的应力强度因子，数值结果与现有文献的相比，符合很好。     

基于多尺度方法的飞行器结构强度分析

建立了一套飞行器结构多尺度分析方法,能够较为高效、准确地分析飞行器结构的力学行为,确定危险区域以及研究损伤模式。采用层次多尺度法,分别建立了飞行器整体结构、局部舱段结构、单钉连接模型三种有限元模型,对飞行器结构进行分析。建立了宏观变量与细观响应之间的信息传递和反馈。对三向正交碳/碳机织复合材料的宏观刚度和强度性能进行了预测,建立了宏观损伤起始包络线。采用协同多尺度法对单钉连接模型进行渐进损伤分析。研究表明该方法具有良好适用性,能够较为准确地分析结构的损伤模式,为飞行器结构设计提供参考。     

分析二维导体电磁散射问题的FEM/POPTD方法——TE模(英文)

提出了一种分析计算带有裂缝的二维电大尺寸复杂导体目标电磁散射问题的混合方法—— FEM/ PO-PTD方法。该方法采用基于棱边的有限元法 ( Edge-basedFEM)为低频方法 ,物理光学法与物理绕射理论 ( PO-PTD)为高频方法 ,通过耦合技术将两者结合在一起。本文将该方法应用于带有缝隙的二维导电柱 TE模的电磁散射特性分析 ,计算结果与有关文献的数据一致性很好 ,从而验证了该方法的准确性。文中还给出了另外几种截面的导电柱体雷达截面的计算曲线。理论分析与计算结果表明 ,本文提出的混合方法与其他计算同类问题的方法相比 ,能节省计算机存储单元、提高计算速度     

凸轮弹性连杆复合机构的动力学修改

本文以胶印机递纸机构为研究对象，该机构为一两自由度凸轮连杆复合机构。在建立该机构运动弹性动力学（ＫＥＤ）有限元分析模型和现场袂测固有频率的基础上，提出了一种利用实 固有频率数据修正复合机构有限元动力模型的方法。其基本思想是利用机构各个运动位置上的固有频率实测值与有限元模型的计算值之间的加权误差建立目标函数，然后利用优化方法使目标函数极小化，从而求出有限元模型参数的修正值。计算结果表明该方法是简便有