基于特征分裂有限元准隐格式的共轭传热整体耦合数值模拟方法

共轭传热现象在科学和工程领域中大量存在. 随着计算能力的发展, 对共轭传热现象进行准确有效的数值模拟, 成为科学研究和工程设计上的重要挑战.共轭传热数值模拟的方法可以分为两大类: 分区耦合和整体耦合.本文采用有限元法对共轭传热问题进行整体耦合模拟. 固体传热求解采用标准的伽辽金有限元方法.流动求解采用基于特征分裂的有限元方法. 该方法是一种重要的求解流动问题的有限元方法, 可以使用等阶有限元. 该方法的准隐格式与其他格式相比, 具有时间步长大的特点. 将稳定项中的时间步长与全局时间步长分开, 改进了准隐格式的稳定性. 基于改进的特征分裂有限元方法的准隐格式, 发展了一种层流共轭传热数值模拟的整体耦合方法. 采用这种方法可以将流体计算域和固体计算域作为一个整体划分有限元网格, 并且所有变量都可以采用相同的插值函数, 从而有利于程序的实现. 通过对典型问题的模拟, 验证了这种方法的准确性. 本工作还研究了固体区域时间步长对定常共轭传热问题数值模拟收敛性的影响.      

薄膜结构流固耦合的CFD数值模拟研究

基于弱耦合分区求解策略,在CompaqVisualFortran6．5环境下搭建了薄膜结构三维流固耦合效应的CFD数值模拟平台。程序采用模块化编程思想,主要包含几何建模、流体分析、结构分析和数据交换四个模块。其中几何建模模块采用自行编制的膜结构找形分析程序,流体分析模块采用经过二次开发的计算流体力学软件FLUENT6．0,结构分析模块采用自行编制的膜结构动力分析程序MDLFX;在数据交换模块中,编制了基于薄板样条法的插值计算程序,以实现流固交界面上不同区域网格间的数据传递问题,编制了基于代数法和迭代法的动网格变形程序,以实现流固耦合运算中的动网格更新。基于该软件平台,对单向柔性屋盖和鞍形膜结构屋盖进行了流固耦合数值模拟,验证了方法的有效性。     

航行体水下发射流固耦合效应分析

对于水下发射过程来说,掌握水动力载荷形成机理与结构响应特征是一个亟待解决的问题.研究该问题需要考虑含相变的复杂多相流动,变约束的结构运动以及这二者之间的耦合效应.本文采用松耦合的方法,以流体求解器为主体,将自编的固体结构程序接入流体求解器中,在每个时间步长内分别对流体动力学方程和固体结构动力学方程进行求解,通过流固界面之间的数据交换实现耦合计算.其中,流体求解器基于雷诺平均纳维斯托克斯方程,采用单流体模型处理多相流问题,引入空化模型描述空化相变,采用修正的湍流模型模拟混合物的湍流效应,并采用动网格技术处理移动边界问题.航行体的刚体运动和结构振动分开求解.结构求解器采用等效梁模型描述结构的振动,通过坐标变换给出了随体坐标系下的结构振动方程,求解方法采用时域积分法.所建立的流固耦合方法不仅能够捕捉到自然空化的演化情况,还可获得航行体所受水动力、结构振动响应以及截面的弯矩,获得了实验的验证.基于该方法研究了结构刚度、发射速度对空泡溃灭与结构振动耦合效应的影响规律.结果表明,同步溃灭是影响结构载荷的主要因素,包括溃灭压力幅值,溃灭压力作用位置,以及溃灭压力与结构振动的相位关系.     

风与柔性结构流固耦合作用的预处理方法研究

针对强耦合方法求解风与柔性结构流固耦合作用时,大量计算资源都耗费在对强耦合方程求解中这一弊端,本文研究了强耦合方程的预处理求解方法。在风与柔性结构流固耦合作用的强耦合整体方程的基础上,将时空离散和线性化后的类似结构方程看成是鞍点问题,首先推导得到了类似结构方程的预处理矩阵;再基于此推导出了强耦合整体方程的预处理矩阵。首先采用预处理方法对经典二维流固耦合问题进行了计算,验证了提出的预处理矩阵的正确性;然后对风与三维膜结构的流固耦合作用进行了分析,评估了所提出预处理方法的相关计算参数。计算结果表明,所提出的预处理方法可使强耦合整体方程的求解在计算精度和计算效率上都得到较大提升,证明本文提出的预处理方法适用于风与柔性结构的流固耦合分析。     

无限长圆柱体中广义电磁热弹耦合问题研究

基于L-S广义热弹性理论,研究了处于磁场中无限长理想圆柱导体在边界受热冲击作用时的电磁热弹耦合问题.建立了广义电磁热弹耦合的有限元方程,为避免积分变换方法求解带来的精度丟失.采用将有限元方程直接在时间域求解的方法,得到了圆柱体中的温度、位移、应力、感应磁场和感应电场的分布规律,反映了热的波动性及电磁热弹的耦合效应.结果表明,将有限元方程直接在时间域求解,可以获得各物理量的准确分布.得到温度在热波波前处的阶跃,准确地反应热波的波动效应.     

求解三维流固耦合问题的一种全隐全耦合区域分解并行算法

三维流固耦合问题的非结构网格数值算法在很多工程领域都有重要应用,目前现有的数值方法主要基于分区算法,即流体和固体区域分别进行求解,因此存在收敛速度较慢以及附加质量导致的稳定性问题,此外,该类算法的并行可扩展性不高,在大规模应用计算方面也受到一定限制．本文针对三维非定常流固耦合问题,提出一种基于区域分解的全隐全耦合可扩展并行算法．首先基于任意拉格朗日-欧拉框架建立流固耦合控制方程,然后时间方向采用二阶向后差分隐式格式、空间方向采用非结构稳定化有限元方法进行离散．对于大规模非线性离散系统,构造一种结合非精确Newton法、Krylov子空间迭代法与区域分解Schwarz预条件子的Newton-Krylov-Schwarz (NKS)并行求解算法,实现流体、固体和动网格方程的一次性整体求解．采用弹性障碍物绕流的标准测试算例对数值方法的准确性进行了验证,数值性能测试结果显示本文构造的全隐全耦合算法具有良好的稳定性,在不同的物理参数下具有良好的鲁棒性,在“天河二号”超级计算机上,当并行规模从192增加到3072个处理器核时获得了91%的并行效率．性能测试结果表明本文构造的NKS算法有望应用于复杂...     

流体-结构耦合问题的有限元并行计算研究

流体-结构耦合问题广泛存在于各种工程领域,本文采用ALE显式有限元法求解该类问题,并对该方法的并行性进行讨论。同时根据流体-结构耦合问题与ALE显式有限元的计算特点,在坐标递归分区方法的基础上设计并程序实现了基于流体-结构耦合均衡的分区算法。通过与坐标递归分区方法的计算结果相比较,对于流体-结构耦合问题的求解,耦合均衡并行分区方法具有更好的加速比和并行效率。     

精密离心机流-热-固多场耦合计算方法

离心机是惯性导航系统加速度计的标定设备,转盘变形将严重影响标定精度。基于耦合理论研究了离心机流-热-固多场耦合问题的数值计算方法,分析了此类问题的全耦合数学模型,研究了求解数学方程的强、弱耦合算法及其应用。根据耦合方式不同定义了三种分析方法:直接法、顺序迭代法和顺序耦合法,建立了它们与强、弱耦合算法的联系。验证方法正确性后,对比离心机温度场强、弱耦合算法的计算结果发现,强耦合算法计算时间长,但模拟细节好,更接近真实情况,弱耦合算法计算时间短,但精度较低,故采用强耦合算法计算离心机流场和温度场。三种分析方法计算热-结构耦合的结果相近,故选择顺序耦合法计算离心机结构变形。     

飞行器流场与结构温度场耦合数值分析

利用高精度的有限体积迎风格式TVD来求解超声速流场，利用有限元方法求解固壁的温度场，流场与结构温度场互为边界条件交换数据，实现流场解算与温度场解算的耦合数值分析，耦合方法研究表明气动加热问题属于强耦合，耦合迭代直接影响温度平衡时间，因此采用小时问步的耦合计算可以准确预测高速飞行器的气动加热。     

MPS-FEM数值分析带自由面的流固耦合问题

随着计算科学的发展,研究人员为探索流固耦合问题的物理机理而提出了众多的数值方法。其中,耦合的移动粒子半隐式方法 MPS(Moving Particle Semi-Implicit method)和有限单元法FEM(Finite Element method)为流固耦合问题的数值仿真工作提供了新的途径。本文所有流场的数值模拟工作均采用课题组自主开发的无网格法求解器MLParticle-SJTU来完成。该求解器在原始的MPS法基础上,对核函数、压力梯度模型、压力泊松方程的求解和自由面判断方式等方面进行了改进。此外,在该求解器框架内,基于FEM法拓展了针对结构场进行求解的功能。首先,对MPS和FEM方法的理论模型及其耦合策略进行了介绍。然后,采用该自研MPS-FEM耦合求解器,数值模拟了溃坝流动对弹性结构的冲击及其相互作用的标准问题。通过将结构变形及自由面波型变化等结果与已发表结果进行对比,验证了该求解器在处理带自由面剧烈变化的粘性流体和柔性变形结构的耦合作用问题上的可行性。     

Aeroelastic stability analysis of a flexible over-expanded rocket nozzle using numerical coupling by the method of transpiration

The aim of this paper is to present and compare two different approaches for aeroelastic stability analysis of a flexible over-expanded rocket nozzle. The first approach is based on the aeroelastic stability models developed in a previous work, while the second uses the numerical fluid–structure coupling via the transpiration method. The aeroelastic frequencies of the nozzle obtained by various stability models are compared with those extracted from the numerical coupling by the method of transpiration. Both set of results show an overall good agreement.     

基于流形切空间插值的折叠翼参数化气动弹性建模

变体飞行器的气动弹性力学建模是当前先进飞行器设计的研究热点和难点. 然而传统的气动弹性动力学建模方法对于具有结构参变特性的变体飞行器气动弹性力学研究存在建模效率低、计算复杂等问题. 本研究提出了一种基于流形切空间插值的可折叠式变体机翼参数化气动弹性建模方法. 首先, 该方法建立若干个典型折叠角下的折叠翼结构有限元模型, 通过流形切空间插值方法建立折叠翼参数化结构动力学模型. 其次, 采用偶极子网格法得到参数化非定常气动力模型, 进而建立气动和结构相互耦合的折叠翼参数化气动弹性模型. 为了验证该参数化建模方法在折叠翼气动弹性分析中的准确性, 本文以一小展弦比折叠翼为研究对象, 从折叠翼自由振动时的参变模态特性、颤振边界预测两方面进行了算例验证, 并与直接计算方法进行了对比, 进一步验证了参数化气动弹性模型的有效性. 研究结果表明, 该参数化气动弹性模型对上述两类问题的计算精度与直接计算方法一致, 并且有着计算效率更高的优势.      

Aeroelastic Wind Forces on Flexible Bridge Girders

Non-stationary aeroelastic wind force processes on a streamlined bridge deck are modelled for use within computations in the range of statistically stationary girder oscillation and for the mathematical investigation of the limit of aeroelastic stability. Time Domain (TD) methods, as the direct integration on a finite-element (FE) model, are applied in structural analyses for taking into consideration both geometric and physical nonlinearities. The FE method requires a continuous form of the force model, which permits an arbitrary discretization in time. The present model is based on continuous step response functions of Wagner-type, taking advantage of modern advances of that theory. Special emphasis is put on the continuity and the existence of equivalent representations in time and frequency domains. A realistic wind load model must include measured force functions and effects from the wind turbulence. For that reason, the force functions are calibrated in wind tunnel tests on section models in grid turbulence and boundary layer flow. The aeroelastic in-phase and out-of-phase forces are identified from highly accurate calibrations experiments on section models up to high frequencies which permits to decrease the time step increment. Turbulence effects on the aeroelastic coefficients are included in mean terms.     

基于POD和代理模型的静气动弹性分析方法

静气动弹性问题考虑弹性结构与定常气动力间的相互耦合作用,对飞行器的性能和安全具有显著的影响.在现代飞行器设计阶段,计算流体力学(CFD)/计算结构力学(CSD)直接耦合方法是精确考察静气动弹性影响的重要手段.然而,基于CFD技术的气动力仿真手段在耦合过程中计算量大且耗时长,难以满足设计阶段的需求.因此,为了兼顾计算精度与效率,文章采用本征正交分解(POD)和Kriging代理模型相结合的模型降阶方法,替代CFD求解过程并耦合有限元分析(FEA)方法,建立了高效、准确的静气动弹性分析框架.相较于传统的以模态法为主的静气动弹性分析方法,该方法能够解决更为复杂的静气动弹性问题以及提供静气动弹性变形过程中的气动分布载荷.针对典型三维跨声速HIRENASD机翼模型开展的马赫数、迎角变化的算例验证表明:由建立的静气动弹性分析方法与CFD/CSD直接耦合方法计算得到机翼翼梢处的静变形量间的相对误差在5%以内;同时该方法预测静平衡位置处的气动分布载荷的误差在5%以内,静气动弹性分析的计算效率至少提升了6倍.     

Explicit/implicit fluid/structure staggered procedures with a structural predictor and fluid subcycling for 2D inviscid aeroelastic simulations

Field time integrators with second-order-accurate numerical schemes for both the fluid and the structure are considered for unsteady Euler aeroelastic computations. We show that if these schemes are simply coupled and used straightforwardly with subcycling, then accuracy and stability properties may be lost. We present new coupling staggered procedures where momentum conservation is enforced at the interface. This is done by using a structural predictor. Continuity of structural and fluid grid displacements is not satisfied at the fluid/structure interface. However, we show on a two-degree-of-freedom aerofoil that this new type of method has many advantages, e.g. accuracy of conservation at the interface and extended stability. The supersonic flutter of a flat panel is simulated in order to numerically prove that the algorithm gives accurate results with arbitrary subcycling for the fluid in the satisfying limit of 30 time steps per period of coupled oscillation. © 1997 John Wiley & Sons, Ltd.     

Nonlinear aeroelastic stability analysis of wind turbine blade with bending–bending–twist coupling

In this study, the nonlinear aeroelastic stability of wind turbine blade with bending–bending–twist coupling has been investigated for composite thin-walled structure with pretwist angle. The aerodynamic model used here is the differential dynamic stall nonlinear ONERA model. The nonlinear aeroelastic equations are reduced to ordinary equations by Galerkin method, with the aerodynamic force decomposition by strip theory. The nonlinear resulting equations are solved by a time-marching approach, and are linearized by small perturbation about the equilibrium point. The nonlinear aeroelastic stability characteristics are investigated through eigenvalue analysis, nonlinear time domain response, and linearized time domain response.     

An aeroelastic model for composite rotor blades with straight and swept tips. Part I: Aeroelastic stability in hover

An analytical model for predicting the aeroelastic behavior of composite rotor blades with straight and swept tips is presented. The blade is modeled by beam type finite elements along the elastic axis. A single finite element is used to model the swept tip. The non-linear equations of motion for the finite element model are derived using Hamilton's principle and based on a moderate deflection theory and accounts for: arbitrary cross-sectional shape, pretwist, generally anisotropic material behavior, transverse shears and out-of-plane warping. Numerical results illustrating the effects of tip sweep, anhedral and composite ply orientation on blade aeroelastic behavior are presented. It is shown that composite ply orientation has a substantial effect on blade stability. At low thrust conditions, certain ply orientations can cause instability in the lag mode. The flap-torsion coupling associated with tip sweep can also induce aeroelastic instability in the blade. This instability can be removed by appropriate ply orientation in the composite construction.     

A grid deformation technique for unsteady flow computations

A grid deformation technique is presented here based on a transfinite interpolation algorithm applied to the grid displacements. The method, tested using a two‐dimensional flow solver that uses an implicit dual‐time method for the solution of the unsteady Euler equations on deforming grids, is applicable to problems with time varying geometries arising from aeroelasticity and free surface marine problems. The present work is placed into a multi‐block framework and fits into the development of a generally applicable parallel multi‐block flow solver. The effect of grid deformation is examined and comparison with rigidly rotated grids is made for a series of pitching aerofoil test cases selected from the AGARD aeroelastic configurations for the NACA0012 aerofoil. The effect of using a geometric conservation law is also examined. Finally, a demonstration test case for the Williams aerofoil with an oscillating flap is presented, showing the capability of the grid deformation technique. Copyright © 2000 John Wiley & Sons, Ltd.     

Parametric aeroelastic modeling based on component modal synthesis and stability analysis for horizontally folding wing with hinge joints

To investigate the aeroelastic stability of a folding wing effectively, a parametric aeroelastic analysis approach is proposed. First, the fixed interface component modal synthesis is used to derive the structural dynamic equation for a folding wing, in which the elastic connection is considered. The unsteady aerodynamic model is established by the doublet lattice method (DLM), and the aeroelastic model is achieved from integration of the DLM with the component modal analysis. The generalized aerodynamic influence coefficient matrix is established by modes kept and constraint modes of each component. The aeroelastic stability of a folding wing is investigated based on the Gram matrix in control theory. The effectiveness of the proposed method is verified via comparison with traditional flutter eigenvalue analysis for both extended and folded configurations. The proposed method identifies coupled modes and improves computational efficiency when compared to classical aeroelastic stability analysis methods, such as the p–k method.     

Aeroelastic analysis of helicopter rotor blade in hover using an efficient reduced-order aerodynamic model

This paper presents a coupled flap–lag–torsion aeroelastic stability analysis and response of a hingeless helicopter blade in the hovering flight condition. The boundary element method based on the wake eigenvalues is used for the prediction of unsteady airloads of the rotor blade. The aeroelastic equations of motion of the rotor blade are derived by Galerkin's method. To obtain the aeroelastic stability and response, the governing nonlinear equations of motion are linearized about the nonlinear steady equilibrium positions using small perturbation theory. The equilibrium deflections are calculated through the iterative Newton–Raphson method. Numerical results comprising steady equilibrium state deflections, aeroelastic eigenvalues and time history response about these states for a two-bladed rotor are presented, and some of them are compared with those obtained from a two-dimensional quasi-steady strip aerodynamic theory. Also, the effect of the number of aerodynamic eigenmodes is investigated. The results show that the three-dimensional aerodynamic formulation has considerable impact on the determination of both the equilibrium condition and lead-lag instability.