线化欧拉方程的高阶间断有限元数值解法研究

采用高阶间断有限元法于非结构网格上针对复杂外形数值求解声学控制方程------线化欧拉方程. 背景流场采用有限体积法于结构网格求得, 一种高精度数据传递方法将基于有限体积法的背景流场数据传递到声场计算所采用的较为稀疏的非结构网格上, 保证了背景流场信息的完整和精确. 为提高计算效率, 采用了一种更为直接的Quadrature-FreeImplementation技术以及网格分区并行技术. 数值结果表明采用高阶的情况下即使在稀疏的网格上也可以捕捉到细微的声场结构.      

机动弹头高超声速粘性绕流数值模拟

采用空间二阶精度的交替方向隐式分解的NND格式求解完全气体假定下的非定常薄层近似NaVierStokes方程,并采用抛物化的椭圆型方程生成复杂带翼弹头的空间网格.最后给出了机动弹头在M∞=7.3,α=20°下的流场计算结果.计算结果表明本文所采用的计算网格及计算方法可以适用于高超声速复杂流场的数值模拟.     

CRACK FACE CONTACT PROBLEM ANALYSIS USING THE SCALED BOUNDARY FINITE ELEMENT METHOD

比例边界有限元侧面上有任意荷载时,将侧面载荷分解成关于径向方向局部坐标的多项式函数的和,推导给出了考虑侧面载荷存在的新型形函数,并基于该形函数推导了刚度矩阵和等效节点载荷列阵.首次对比例边界有限元法求解裂纹面接触问题进行了研究,运用Lagrange乘子引入接触界面约束条件,推导给出了比例边界有限元求解裂纹面接触问题的控制方程.将裂纹面单元分为非裂尖单元和含有侧面的裂尖单元.在非裂尖单元中的裂纹面,裂纹面作为多边形单元的边界,边界上的接触力可等效到节点上,通过在节点上构造Lagrange乘子,采用点对点接触约束进行处理.对于含有侧面的裂尖单元,在整个侧面上构造Lagrange乘子的插值场,采用边对边接触约束进行处理.对三个不同的接触约束状态下的算例进行了数值计算,通过与解析解及有限元软件ABAQUS计算结果的对比,验证了本文提出的比例边界有限元点对点和边对边接触求解裂纹面接触问题的精确性与有效性.     

比例边界有限元法求解裂纹面接触问题

比例边界有限元侧面上有任意荷载时,将侧面载荷分解成关于径向方向局部坐标的多项式函数的和,推导给出了考虑侧面载荷存在的新型形函数,并基于该形函数推导了刚度矩阵和等效节点载荷列阵.首次对比例边界有限元法求解裂纹面接触问题进行了研究,运用Lagrange乘子引入接触界面约束条件,推导给出了比例边界有限元求解裂纹面接触问题的控制方程.将裂纹面单元分为非裂尖单元和含有侧面的裂尖单元.在非裂尖单元中的裂纹面,裂纹面作为多边形单元的边界,边界上的接触力可等效到节点上,通过在节点上构造Lagrange乘子,采用点对点接触约束进行处理.对于含有侧面的裂尖单元,在整个侧面上构造Lagrange乘子的插值场,采用边对边接触约束进行处理.对三个不同的接触约束状态下的算例进行了数值计算,通过与解析解及有限元软件ABAQUS计算结果的对比,验证了本文提出的比例边界有限元点对点和边对边接触求解裂纹面接触问题的精确性与有效性.     

非结构动网格在多介质流体数值模拟中的应用

采用非结构动网格方法对含多介质的流场进行数值模拟.采用改进的弹簧方法来处理由于边界运动而产生的网格变形.采用基于格心的有限体积方法求解守恒型的ALE(Arbitrary Lagrangiall-Eulerian)方程,控制面通量的计算采用HLLC(Hartem,Lax,van Leer,Contact)方法,采用几何构造的方法使空间达到二阶精度,时间离散采用四阶Runge-Kutta方法.物质界面的处理采用虚拟流体方法.本文对含动边界的激波管、水下爆炸等流场进行数值模拟,取得较好的结果,不同时刻界面的位置和整个扩张过程被准确模拟.     

三维复杂流动的间断有限元方法模拟

本文基于三维可压缩Euler方程,采用基于Runge-Kutta时间离散的间断有限元方法(RKDG方法),对三维前台阶、三维Riemann问题和球Riemann等问题进行了模拟。结果表明,本文的RKDG方法能够在很少的网格内清晰地捕捉到三维复杂流场中的激波和接触间断;同时,将球Riemann问题中z=0.4平面压强沿到对称轴距离的分布与文献中的近似精确解相比,吻合较好,这也验证了本文的RKDG方法不仅能够进行三维复杂流场的定性描述,也能够应用于三维复杂流场的定量计算。     

A NEW SPATIAL THIN-WALLED BEAM ELEMENT INCLUDING TRANSVERSE AND TORSIONAL SHEAR DEFORMATION

基于Timoshenko梁及Benscoter薄壁杆件理论,建立了考虑剪切变形、弯扭耦合以及翘曲剪应力影响的空间任意开闭口薄壁截面梁单元. 通过引入单元内部结点,对弯曲转角和翘曲角采用三节点Lagrange独立插值的方法,考虑了剪切变形和翘曲剪应力的影响并避免了横向剪切锁死问题;借助载荷作用下薄壁梁的截面运动分析,在位移和应变方程中考虑了弯扭耦合的影响. 通过数值算例将该单元的计算结果与理论解以及商用有限元软件和其他文献中的数值解进行对比和验证,结果对比表明该薄壁梁单元具有良好的精度和收敛性.     

预应力混凝土空间组合式滑移模型

考虑预应力钢筋和混凝土间的粘结滑移效应,结合组合式模型和分离式模型的优点,提出了一种预应力混凝土空间组合式滑移模型。通过预应力钢筋和混凝土单元交点处的虚拟结点,将模拟预应力钢筋和混凝土交界面的无厚度环形粘结单元嵌入组合单元中。粘结单元忽略预应力钢筋-混凝土间的径向相对位移,通过内、外表面的切向位移差值表征预应力钢筋的相对滑移量。基于位移有限元框架,根据预应力钢筋、混凝土、粘结单元各自的本构模型,由虚功原理推导出三者对组合单元刚度矩阵的贡献,建立了该空间组合式滑移模型的有限元平衡方程,并给出了有限元计算流程。在该模型中,预应力钢筋能以任意形式穿过混凝土单元而不需考虑其走向,网格划分灵活;对单个粘结单元而言,在平衡方程中缩减了2个虚拟结点的6个自由度,自由度数减少了8%。     

基于CEL有限元法的大管径钢管桩贯入海基过程的承载力分析

鉴于耦合欧拉-拉格朗日(Coupled Eulerian-Lagrangian, CEL)有限元方法在计算分析接触大变形问题时的优势,建立了大管径钢管桩贯入海基的CEL有限元模型,即在被贯入的海基层采用欧拉有限元法建模,其它层采用拉格朗日有限元法建模,两者之间通过罚函数的接触算法进行耦合．在欧拉有限元模型中,采用算子分裂技术．为兼顾计算模型的精度和效率,本文先对模型的网格密度和贯入速度进行稳定性分析,以此确定合适的网格模型和适宜的贯入速度,同时验证CEL方法在贯入模拟中的可行性．随后针对不同土体类型和不同土体参数变化,探究土体贯入阻力随贯入深度的变化规律．计算结果表明砂土的贯入阻力最大,粉砂土次之,粘土的贯入阻力最小;且贯入阻力对土体内摩擦角的变化最为敏感,对土体粘聚力的变化比较敏感,而对土体模量的变化不敏感．     

大展弦比无人机的静气弹问题计算及分析

普遍采用大展弦比机翼的无人机,其气动弹件的问题显得尤为突出.本文主要研究基于非结构弹性网格体系的欧拉方程CFD计算及其在大展弦比无人机静气动弹性问题中的应用.针对三维非结构运动网格技术进行了研究和开发,在此基础上,利用计算流体力学的方法,发展了一套具有一定通用性的、适用于非结构网格的Euler方程流场求解器,并综合上面的技术,进一步通过耦合结构力学方程,对大展弦比无人机的静气动弹性问题进行了,计算和分析.     

Simulation of liquid sloshing in curved‐wall containers with arbitrary Lagrangian–Eulerian method

There are many challenges in the numerical simulation of liquid sloshing in horizontal cylinders and spherical containers using the finite element method of arbitrary Lagrangian–Eulerian (ALE) formulation: tracking the motion of the free surface with the contact points, defining the mesh velocity on the curved wall boundary and updating the computational mesh. In order to keep the contact points slipping along the curved side wall, the shape vector in each time advancement is defined to modify the kinematical boundary conditions on the free surface. A special function is introduced to automatically smooth the nodal velocities on the curved wall boundary based on the liquid nodal velocities. The elliptic partial differential equation with Dirichlet boundary conditions can directly rezone the inner nodal velocities in more than a single freedom. The incremental fractional step method is introduced to solve the finite element liquid equations. The numerical results that stemmed from the algorithm show good agreement with experimental phenomena, which demonstrates that the ALE method provides an efficient computing scheme in moving curved wall boundaries. This method can be extended to 3D cases by improving the technique to compute the shape vector. Copyright © 2007 John Wiley & Sons, Ltd.     

An extended finite element method for the simulation of particulate viscoelastic flows

We present an extended finite element method (XFEM) for the direct numerical simulation of the flow of viscoelastic fluids with suspended particles. For moving particle problems, we devise a temporary arbitrary Lagrangian–Eulerian (ALE) scheme which defines the mapping of field variables at previous time levels onto the computational mesh at the current time level. In this method, a regular mesh is used for the whole computational domain including both fluid and particles. A temporary ALE mesh is constructed separately and the computational mesh is kept unchanged throughout the whole computations. Particles are moving on a fixed Eulerian mesh without any need of re-meshing. For mesh refinements around the interface, we combine XFEM with the grid deformation method, in which nodal points are redistributed close to the interface while preserving the mesh topology. Our method is verified by comparing with the results of boundary fitted mesh problems combined with the conventional ALE scheme. The proposed method shows similar accuracy compared with boundary fitted mesh problems and superior accuracy compared with the fictitious domain method. If the grid deformation method is combined with XFEM, the required computational time is reduced significantly compared to uniform mesh refinements, while providing mesh convergent solutions. We apply the proposed method to the particle migration in rotating Couette flow of a Giesekus fluid. We investigate the effect of initial particle positions, the Weissenberg number, the mobility parameter of the Giesekus model and the particle size on the particle migration. We also show two-particle interactions in confined shear flow of a viscoelastic fluid. We find three different regimes of particle motions according to initial separations of particles.     

Adaptive large eddy simulation with moving grids

The paper presents adaptive mesh moving methods for large eddy simulation (LES) of turbulent flows. With this approach, a given number of grid points is redistributed with respect to an appropriately selected criterion. The Arbitrary Lagrangian–Eulerian formulation is applied to solve the governing equation on moving grids employing a collocated finite volume formulation. A dynamic moving mesh partial differential equation based on a variational principle is solved for the corner points of the grid by means of a dedicated solver. Adaptation is performed in a statistical sense so that statistical quantities of interest are employed. Various LES-specific design criteria and combination of them are proposed, such as the time-averaged gradient of streamwise velocity, turbulent kinetic energy and production rate. These are investigated in the framework of elementary and balanced monitor functions. These are tested for the three-dimensional flow in a channel with periodic constrictions. The numerical results are compared to a highly resolved LES reference solution. The independence of the moving mesh method from the initial LES is shown, and its potential to improve the efficient resolution of turbulent flow features is demonstrated.     

黏弹性人工边界等效荷载计算的改进方法

黏弹性人工边界在场地地震反应和结构-地基动力相互作用等问题的计算中已得到了广泛的应用.地震波在黏弹性人工边界中的输入是通过将地震波转化为作用于人工边界处的等效载荷来实现的.计算等效节点载荷的常规方法默认边界节点对应区域内的应力为均布力,但实际上该节点对应区域内的应力分布通常是不均匀的.本文在有限元方法结合黏弹性局部人工边界的显式时域波动方法的基础上,建立了无限域散射问题地震波等效载荷计算的一种改进方法.该方法采用细化网格与应力积分相结合的方法计算人工边界等效节点力,有效地降低了人工边界上等效节点力的计算误差.以不同角度入射地震波的二维算例为例,算例给出的波场位移云图和节点位移时程曲线验证了本文方法的有效性,其计算精度与网格尺寸和地震波入射角度密切相关,且网格越小、入射角度越小,计算精度越高.对于相同的网格尺寸,本文采用方法的计算精度明显高于常规方法,尤其是对于斜入射问题优势更为明显.     

A marker‐and‐cell approach to free surface 2‐D multiphase flows

This work describes a methodology to simulate free surface incompressible multiphase flows. This novel methodology allows the simulation of multiphase flows with an arbitrary number of phases, each of them having different densities and viscosities. Surface and interfacial tension effects are also included. The numerical technique is based on the GENSMAC front‐tracking method. The velocity field is computed using a finite‐difference discretization of a modification of the Navier–Stokes equations. These equations together with the continuity equation are solved for the two‐dimensional multiphase flows, with different densities and viscosities in the different phases. The governing equations are solved on a regular Eulerian grid, and a Lagrangian mesh is employed to track free surfaces and interfaces. The method is validated by comparing numerical with analytic results for a number of simple problems; it was also employed to simulate complex problems for which no analytic solutions are available. The method presented in this paper has been shown to be robust and computationally efficient. Copyright © 2012 John Wiley & Sons, Ltd.     

Efficient modeling of cable-pulley system with friction based on arbitrary-Lagrangian-Eulerian approach

In conventional modeling of a cable-pulley system, the cable must be finely meshed with Lagrangian elements for valid contact detections with pulleys, leading to extremely low efficiency. The sliding joint method based on the arbitrary-LagrangianEulerian(ALE) formulation still lacks an efficient cable element, and in particular, modeling of friction between a sliding joint and the cable has not been studied. This paper presents efficient multi-body modeling of a cable-pulley system with friction. A variablelength cable element with a node movable along the cable, which is described with ALE,is developed to mesh the cable. A transitional cable element is then proposed to model the contact part of the cable by fixing its two nodes to the two corresponding locations of the pulley. Friction of the cable-pulley is derived as a simple law of tension decay and embedded in the multi-body system modeling. It is simplified as a generalized friction force acting only on the arc-length coordinate. This approach can use a rough mesh on the cable, and is free of contact detections, thus significantly saving computation time.Several examples are presented to validate the proposed method, and show its effectiveness in real engineering applications.     

Validation of a fluid–structure interaction model for a bileaflet mechanical heart valve

A fluid-structure interaction (FSI) model for heart valve simulation is presented. In a partitioned framework, separate fluid and structure solvers are weakly coupled, which in combination with the use of artificial compressibility in the fluid solver, leads to a stable and efficient approach. An Arbitrary Lagrangian Eulerian formulation is employed in the fluid solver to permit the accurate calculation of shear stresses next to the valve boundary. The mesh quality is maintained through a combination of smoothing and local remeshing in 3D. The FSI algorithm is validated on experiments of an idealised quasi-2D mechanical heart valve, and the efficiency of the remeshing approach is demonstrated on a realistic 3D heart-valve geometry.     

Prediction of extrudate swell in polymer melt extrusion using an Arbitrary Lagrangian Eulerian (ALE) based finite element method

Accurate prediction of extrudate (die) swell in polymer melt extrusion is important as this helps in appropriate die design for profile extrusion applications. Extrudate swell prediction has shown significant difficulties due to two key reasons. The first is the appropriate representation of the constitutive behavior of the polymer melt. The second is regarding the simulation of the free surface, which requires special techniques in the traditionally used Eulerian framework. In this paper we propose a method for simulation of extrudate swell using an Arbitrary Lagrangian Eulerian (ALE) technique based finite element formulation. The ALE technique provides advantages of both Lagrangian and Eulerian frameworks by allowing the computational mesh to move in an arbitrary manner, independent of the material motion. In the present method, a fractional-step ALE technique is employed in which the Lagrangian phase of material motion and convection arising out of mesh motion are decoupled. In the first step, the relevant flow and constitutive equations are solved in Lagrangian framework. The simpler representation of polymer constitutive equations in a Lagrangian framework avoids the difficulties associated with convective terms thereby resulting in a robust numerical formulation besides allowing for natural evolution of the free surface with the flow. In the second step, mesh is moved in ALE mode and the associated convection of the variables due to relative motion of the mesh is performed using a Godunov type scheme. While the mesh is fixed in space in the die region, the nodal points of the mesh on the extrudate free surface are allowed to move normal to flow direction with special rules to facilitate the simulation of swell. A differential exponential Phan Thien Tanner (PTT) model is used to represent the constitutive behavior of the melt. Using this method we simulate extrudate swell in planar and axisymmetric extrusion with abrupt contraction ahead of the die exit. This geometry allows the extrudate to have significant memory for shorter die lengths and acts as a good test for swell predictions. We demonstrate that our predictions of extrudate swell match well with reported experimental and numerical simulations.     

An analytically-numerical approach for the analysis of an interface crack with a contact zone in a piezoelectric bimaterial compound

An interface crack with an artificial contact zone at the right-hand side crack tip between two dissimilar finite-sized piezoelectric materials is considered under remote mixed-mode loading. To find the singular electromechanical field at the crack tip, an asymptotic solution is derived in connection with the conventional finite element method. For mechanical loads, the stress intensity factors at the singular points are obtained. As a particular case of this solution, the contact zone model (in Comninou’s sense) is derived. A simple transcendental equation and an asymptotic formula for the determination of the real contact zone length are derived. The dependencies of the contact zone lengths on external load coefficients are illustrated in graphical form. For a particular case of a short crack with respect to the dimensions of the bimaterial compound, the numerical results are compared to the exact analytical solutions, obtained for a piezoelectric bimaterial plane with an interface crack.     

Finite element analysis of flow and heat transfer with moving free surface using fixed grid system

A numerical study was performed on flow and heat transfer involving moving free surfaces that occurs in mold filling processes such as casting and injection molding. In these problems, the calculation domain changes continuously and the numerical treatment of the moving interface tends to cause artificial diffusion. Among the solution algorithms based on the Eulerian method, the volume-of-fluid (VOF) method was used because the method is simple and efficient in handling the complex flow patterns inside the cavity. To solve the transport equation of free surface without artificial smearing of the interface the baby-cell method was employed in the geometric reconstruction of the free surface. Furthermore, a predictor–corrector method was adopted in the time integration of volume-of-fluid (VOF) transport equation to increase the accuracy. The proposed scheme was verified through several benchmark problems. In order to show the capability of the proposed method, several three-dimensional mold filling processes were solved. The current algorithm was applied to the floating body problem. Three-dimensional floating body problems were tested.