基于四叉树网格下的浸入边界-格子Boltzmann方法

提出了一种基于非均匀四叉树网格下的浸入边界-格子Boltzmann方法.在不同层网格界面上,为了保证物理量的连续性,需要在时空方向上进行插值.由于四叉树同层网格在空间上步长相等、且相邻时间层上时间步长相等,所以在时空方向上采用平均值插值,这样做的优点是时空方向插值既可以达到二阶精度,又可以简化计算过程,节省资源,对任意边界加密下的网格,格子Boltzmann方法的实现比较容易.为了充分利用均匀笛卡尔网格的优势,物面边界的处理采用了速度修正法,与传统浸入边界-格子Boltzmann方法中的直接力法、动量交换法相比,无滑移边界条件得到了较好的保证.同时,把大涡模型加入到浸入边界-格子Boltzmann方法中,实现了在四叉树网格数据结构、边界处理技术、大涡模拟几种模型相结合下,绕障碍物的较大雷诺数流动的模拟.通过不可压缩粘性流中圆柱绕流算例验证,结果与其它方法结果吻合良好.     

机翼结构有限元快速建模技术研究

为提高飞机机翼结构初步设计阶段有限元建模的质量和效率,提出了机翼结构有限元参数化建模方法。采用模板技术将机翼结构几何建模过程中的方法、规则和专家经验等知识进行封装,在CATIA建模平台上通过Visual Basic开发了知识驱动的机翼结构几何参数化建模系统。研究了CAD与CAE系统之间的数据传递手段,读取上游模型信息,完成了几何模型的重构。基于PATRAN二次开发语言PCL进行了有限元网格自动剖分,快速生成了有限元模型,并实现了CAD模型向CAE模型转换的无缝集成。     

压力载荷下的结构拓扑-形状协同优化

 压力载荷作用下的结构轻量化设计是工程中的常见问题,由于压力加载面的可设计性,现有以固定载荷为基础的拓扑优化技术不能很好地处理这类问题。直接采用CAD参数化样条或B样条曲线描述压力加载面,通过拓扑和形状变量的联合优化满足了工程实际对结构轻量化与边界的功能性与光滑性设计要求。同时,为了避免结构边界形状变化时有限元网格刷新引起的定义拓扑伪密度变量的困难,用所提出的背景网格和密度点技术实现了每一步单元密度设计迭代结果的自动传递,并采用网格变形技术实现了形状设计变量灵敏度分析。采用4个数值算例验证了方法的有效性,其中发动机承力框架的设计结果充分说明该方法在航空结构设计中的重要应用价值。     

三维直角叉树切割网格Euler方程自适应算法

基于叉树数据结构，实现了一种用于三维直角叉树切割网格的自适应算法，包括对几何外,形以及对流场计算的自适应网格加密技术。在初始网格的生成过程中，根据相邻网格的物面法向向量间的差值，进行针对外形的自适应网格加密；在流场计算中，根据相邻网格间选定物理量梯度的变化，进行针对流场的自适应网格加密。详细地描述了三维直角叉树切割网格的生成过程，以及以任意网格的切割细分算法。在自适应过程中，分别采用了八叉树和全叉树的数据结构，八叉树是基本的数据结构，而全叉树的采用，使网格具有了各向异性的特征，从而大大的减少了自适应网格的数量。采用中心有限体积法，求解Euler方程，并运用上术方法，完成了对外形和流场的逢适应网格加密算法，获得了较好的数值计算结果，证明了自适应算法的正确性，体现了直角叉树切割网格自适应技术的有效性和实用性。     

面向设计的飞机翼面结构有限元自动建模技术

提出一种新的飞机翼面结构分析和优化的建模技术。引入“构件属性”的概念，利用“网状”数据结构来实施管理，并将图形用户界面（GUI）技术与CAD中的基本绘图功能相结合，从而形成具有领域特色的新一代快速和在一定条件下的自动建模系统。建模举例表明：此系统操作简单、使用方便、遵循结构分析与优化建模的专业做法。     

应用混合网格分析复杂积冰翼型气动性能

通过发展一种点-点对接的结构/非结构混合网格生成方法,避免了复杂积冰外形难以划分结构化网格的问题,并对复杂积冰翼型的气动性能进行了分析计算.在靠近积冰边界的内层采用了非结构网格以拟合复杂的边界,在非结构网格外采用结构化网格以节省存储空间和计算时间.为了检验网格对选用湍流模型的影响,整个流动区域分别采用SST和SA湍流模型,求解了雷诺平均N-S方程.数值计算结果表明SST模型更适于模拟复杂分离流动,积冰对翼型的气动性能造成了严重的影响.     

用动态非结构重叠网格法模拟三维多体相对运动绕流

 发展了动态非结构重叠网格技术。本文提出了基于阵面推进的相邻单元搜索算法和适合于高阶空间离散格式的网格间插值关系定义方法，前者提高了非结构重叠网格在模拟多体间具有大幅相对运动的非定常问题时网格处理的效率，而后者保证了网格间边界附近流场求解精度；本文非定常Euler方程的求解采用有限体积法，通量计算采用Osher格式，时间积分采用双时间步长推进。为验证本文提出的方法，分别对定常M6机翼绕流和机翼下外挂物投放过程的非定常绕流进行了模拟,结果表明本文提出的方法提高了非结构重叠网格在模拟多体间具有大幅相对运动问题时的效率，保证了网格边界处的精度，具有较高的工程应用价值。     

基于PCL的导弹吊挂强度分析的参数化方法

基于Patran/Nastran有限元分析软件的后台文件,采用PCL（Patran Comm and Language）对Patran/Nastran软件进行了二次开发,建立了用于导弹吊挂强度分析的参数化、自动化方法。提出的方法实现了吊挂结构强度分析过程中的参数化建模-自定义网格划分-材料边界选择-强度分析-结果查看等整个过程的自动化。方法的应用可以方便开展吊挂结构随意改变尺寸、材料参数、边界等的强度分析。并且方法对于任意结构在有限元分析中尺寸、网格、材料、边界等参数化建模及有限元分析均可适用。     

分析模型参数化建模在飞机多学科优化设计中的应用

 论述分析模型参数化建模在飞机多学科优化设计中的重要意义,提出一种机翼气动与结构分析模型参数化建模方法。优化设计变量与机翼气动计算网格的中间参数变量用于描述机翼气动网格变形情况,并建立结构有限元网格与优化设计变量的空间变化关系,实现机翼气动/结构多学科优化设计模型的参数化描述。气动性能计算采用基于N S方程的计算流体力学方法,利用试验设计方法（DOE）和响应曲面模型（RSM）技术构造气动响应曲面用于优化设计,结构计算调用MSC/NASTRAN完成。采用Pareto遗传算法对一大展弦比复合材料机翼模型进行了多学科优化设计分析,得到Pareto前沿面和可选方案以供决策者选择。     

强制偏流喷管的三维网络生成

运用泊松方程实强了制偏流喷管内的三维网络生成，方程中的源项采用了与边界形状和边界网格点分布无关的数学形式，源项内的自由参数由横向网格线与边界的交角及横向网格线在边界处的局部曲率两个约束条件来确定，实现了由边界上的网格分布直接控制内部的网格点，使生成的网格具有良好的正交性和合理性。     

Boundary condition modelling and identification for cantilever-like structures using natural frequencies

The actual boundary conditions of cantilever-like structures might be non-ideally clamped in engineering practice, and they can also vary with time due to damage or aging. Precise modelling of boundary conditions, in which both the boundary stiffness and the boundary mass should be modelled correctly, might be one of the most significant aspects in dynamic analysis and testing for such structures. However, only the boundary stiffness was considered in the most existing methods. In this paper, a boundary condition modelling and identification method for cantilever-like structures is proposed to precisely model both the boundary stiffness and the boundary mass using sensitivity analysis of natural frequencies. The boundary conditions of a cantilever-like structure can be parameterized by constant mass, constant rotational inertia,constant translational stiffness, and constant rotational stiffness. The relationship between natural frequencies and boundary parameters is deduced according to the vibration equation for the lateral vibration of a non-uniform beam. Then, an iterative identification formulation is established using the sensitivity analysis of natural frequencies with respect to the boundary parameters. The regularization technique is also used to solve the potential ill-posed problem in the identification procedure.Numerical simulations and experiments are performed to validate the feasibility and accuracy of the proposed method. Results show that the proposed method can be utilized to precisely model the boundary parameters of a cantilever-like structure.     

压气机转子叶片颤振边界的预测方法

通过发展的基于能量法的颤振数值预测方法得到了压气机转子叶片的颤振边界.将计算结构动力学(CSD)分析得到的叶片表面节点位移插值到耦合面的流体网格点上,并将设计的多层动网格技术应用于计算流体动力学(CFD)方法,实现叶栅振荡作用下的非定常分析,得到叶片表面的非定常气动功以及模态气动阻尼比.以某第一级压气机转子叶片为例,对...     

A 3D hybrid grid generation technique and a multigrid/parallel algorithm based on anisotropic agglomeration approach

A hybrid grid generation technique and a multigrid/parallel algorithm are presented in this paper for turbulence flow simulations over three-dimensional (3D) complex geometries. The hybrid grid generation technique is based on an agglomeration method of anisotropic tetrahedrons. Firstly, the complex computational domain is covered by pure tetrahedral grids, in which anisotropic tetrahedrons are adopted to discrete the boundary layer and isotropic tetrahedrons in the outer field. Then, the anisotropic tetrahedrons in the boundary layer are agglomerated to generate prismatic grids. The agglomeration method can improve the grid quality in boundary layer and reduce the grid quantity to enhance the numerical accuracy and efficiency. In order to accelerate the convergence history, a multigrid/parallel algorithm is developed also based on anisotropic agglomeration approach. The numerical results demonstrate the excellent accelerating capability of this multigrid method.     

基于结构张量特征的机体损伤区域动态划分

为提升飞机结构损伤维修分析的智能化水平,实现飞机结构损伤模型的高效重构,提出了基于结构张量特征的动态阈值损伤区域划分方法。首先,引入结构张量理论,特征表示图像局部颜色结构纹理,完备损伤图像特征表示;其次,通过构建共性分布的结构张量特征空间,从而使不同损伤图像具有一致性损伤区域划分流程;在此基础上,定义了动态阈值划分算子,并通过对动态阈值划分算子参数的计算,实现结构张量特征空间中的动态阈值划分;最后,选用不同飞机结构损伤实例图像对方法进行了对比验证。验证结果表明,相比于传统灰度方法、固定阈值划分和其他动态阈值划分算子,椭圆算子划分得到的损伤边界连贯完整,能够有效分割微小裂纹,噪点较少,质量较优。该方法运算高效、流程统一、特征表示完备,对不同类型结构损伤的适用性好。     

基于动态嵌套网格的飞行器外挂物投放的数值模拟

运用Euler方程和六自由度动力学方程数值模拟了外挂物投放的气动特性与运行姿态和轨迹。为了模拟外挂物从飞行器分离后的绕流特性和相对运动，采用了动态嵌套网格技术，并发展了一种洞边界建立和外边点的插值单元搜寻蛇行法，提高了计算效率。通过算例，证明了本文方法的有效性和合理性。     

Time-domain nonlinear aeroelastic analysis and wind tunnel test of a flexible wing using strain-based beam formulation

A theoretical formulation for time-domain nonlinear aeroelastic analysis of a flexible wing model is presented and validated by wind tunnel tests. A strain-based beam model for nonlinear structural analysis is combined with the Unsteady Vortex Lattice Method (UVLM) to form the complete framework for aeroelastic analysis. The nonlinear second-order differential equations are solved by an implicit time integration scheme that incorporates a Newton-Raphson sub-iteration technique. An advanced fiber optic sensing technique is firstly used in a wind tunnel for measuring large structural deformations. In the theoretical study, the nonlinear flutter boundary is determined by analyzing the transient response about the nonlinear static equilibrium with a series of flow velocities. The gust responses of the wing model at various gust frequencies are also studied. Comparisons of the theoretical and experimental results show that the proposed method is suitable for determining the nonlinear flutter boundary and simulating the gust response of flexible wings in the time domain.     

基于CFD/CSD耦合方法的旋翼气动弹性载荷计算分析

为提高直升机前飞状态下旋翼非定常气动弹性载荷的预估精度,在旋翼气动弹性综合分析方法中引入旋翼CFD模块,建立了一套基于CFD/CSD松耦合分析的计算方法和程序。为高效解决流固耦合方法中由于桨叶挥舞、扭转等弹性变形带来的旋翼贴体网格变形问题,采用基于代数变换方法的网格变形技术,桨叶运动变形量和旋翼气动力信息通过流固交接面传递。旋翼流场分析方法中,主控方程采用耦合S-A湍流模型的Navier-Stokes方程,围绕旋翼流场的网格采用结构嵌套网格方法生成,无黏通量计算采用Roe格式,时间推进采用双时间法。旋翼结构分析中,考虑旋翼配平,基于Hamilton变分原理和20自由度Timoshenko梁模型求解弹性旋翼非线性运动方程。分别对CSD和CFD方法进行验证,在此基础上,计算了SA349/2旋翼桨叶在前飞状态下的非定常气动力、挥舞弯矩、摆振弯矩和扭转力矩,并与飞行测试数据进行了对比。计算表明:CFD/CSD耦合方法可以显著提高旋翼非定常气动弹性载荷的分析精度,精确捕捉桨叶表面压强峰值、激波位置等,表明本文发展的旋翼CFD/CSD耦合方法可以有效地运用到旋翼气动弹性载荷的预测分析中。     

基于CAD模型的碳纤维复合材料天线反射面精度测量

针对天线反射面研制过程中的精度测量,提出一种基于CAD模型的自动测量技术和反射面数据处理技术的精度测量方法。在测点均匀布局规划和工件与模型对齐的基础上,自动获取反射面的完整测量数据,利用最优化算法使实测数据和理论模型达到最佳匹配,实现反射面面形精度的准确评价并获得变换矩阵。最终通过偏差矢量图表达天线反射面的面形偏差趋势,并利用获得的变换矩阵修正加工基准,有效保证了反射面的成型精度和加工精度,工程应用表明该方法较好地满足了反射面制造过程中的测量需求。