竖直流道宽度对气泡运动行为影响的数值模拟

用数值方法模拟了竖直通道宽度对气泡在液体中的非定常运动、变形以及传热特性的影响。在这个模拟中,界面跟踪采用了VOF方法,并采用PL IC进行界面重构。主流场计算采用有限容积方法将控制方程离散,其中扩散项采用中心差分格式,对流项采用一阶迎风格式。用成熟的S IM PLE算法求解N-S方程的速度与压力的耦合问题。引入CSF模型处理运动界面的表面张力。利用所编制的程序计算了竖直流道中的单个气泡的形状、运动特性以及气泡内外流场与传热特性,并对竖直通道宽度在不同情况下,对气泡的形状、运动特性以及传热特性进行了进一步的研究。得到了一系列有价值的结果,并与实验结果比较。表明数值模拟结果与实验结果吻合的较好。     

粘性流中气泡的界面跟踪模拟及动力学分析

气泡在运动过程中的动力学特性与流场的压力分布密切相关。本文用界面跟踪法(Front-Tracking Method,FTM)对气泡界面及其形状进行数值计算,模拟了粘性流体中单气泡上升的界面形态,并与文献中的实验结果进行对比。结果表明,在同一数量级上的厄特沃什数OE、莫顿数OM形状变化规律基本一致,计算模型的正确性得到验证。对同轴两气泡上升过程进行计算,与VOF的模拟结果进行对比,FTM方法得出的气泡形态变化相同,界面更加清晰;进一步研究同轴两个气泡上升动力学特性,使用清晰的界面位置对应压力场进行分析,从压力的角度分析了其变形、速度变化。结果表明:上下两个气泡都有一个压力中心区域使得气泡运动速度大幅度变化,气泡因界面所处区域的压力梯度分布不均呈现出由球形-子弹状-椭圆形的形态变化;随着时间推移,两个相互作用的气泡间的相互作用力在减弱,逐渐融合。     

一种基于几何重构的Youngs-VOF耦合水平集追踪方法

针对界面追踪方法中拉格朗日方法和欧拉--拉格朗日方法计算效率低、不适用大变形、不能应用于三维数值计算模型等问题,研究了一种效率高、界面清晰、适用于三维模型的计算气液两相界面迁移特性的欧拉运动界面追踪方法,该方法将"米"状相邻单元Youngs方法用于运动界面重构,将Youngs-VOF和水平集通过几何方法耦合,提高运动界面精度,克服了VOF和水平集方法存在的缺陷,避免了利用高阶导数本身的稳定性去求解水平集对流方程和距离函数方程."米"状相邻单元Youngs方法避免了数值耗散、数值色散性以及非线性效应引起的捕捉界面模糊的情况.Youngs-VOF耦合水平集方法既保证了计算界面时的稳定性,与拉格朗日方法相比又提高了计算效率.利用Youngs-VOF耦合水平集方法与VOF方法对单个气泡在水中上升过程数值计算与实验对比并对经典剪切流场中圆形运动界面模型的数值计算,验证了Youngs-VOF耦合水平集方法的有效性并比VOF方法捕捉界面更清晰、锐利;通过对溃坝--自由表面流动过程数值计算并与实验进行对比,验证了Youngs-VOF耦合水平集方法的稳定性以及对三维数值模型的适用性.      

一种基于几何重构的Youngs-VOF耦合水平集追踪方法

针对界面追踪方法中拉格朗日方法和欧拉–拉格朗日方法计算效率低、不适用大变形、不能应用于三维数值计算模型等问题,研究了一种效率高、界面清晰、适用于三维模型的计算气液两相界面迁移特性的欧拉运动界面追踪方法,该方法将"米"状相邻单元Youngs方法用于运动界面重构,将Youngs-VOF和水平集通过几何方法耦合,提高运动界面精度,克服了VOF和水平集方法存在的缺陷,避免了利用高阶导数本身的稳定性去求解水平集对流方程和距离函数方程."米"状相邻单元Youngs方法避免了数值耗散、数值色散性以及非线性效应引起的捕捉界面模糊的情况. Youngs-VOF耦合水平集方法既保证了计算界面时的稳定性,与拉格朗日方法相比又提高了计算效率.利用Youngs-VOF耦合水平集方法与VOF方法对单个气泡在水中上升过程数值计算与实验对比并对经典剪切流场中圆形运动界面模型的数值计算,验证了Youngs-VOF耦合水平集方法的有效性并比VOF方法捕捉界面更清晰、锐利;通过对溃坝–自由表面流动过程数值计算并与实验进行对比,验证了Youngs-VOF耦合水平集方法的稳定性以及对三维数值模型的适用性.     

基于Level Set方法对油水和气水两相界面的数值模拟

为研究两相界面迁移特性,基于LevelSet方法,建立了求解非定常不可压缩两相界面流动的数值方法。计算中使用结构化网格采用LevelSet函数捕捉两相界面。通过对经典算例的模拟,验证了数值方法对界面捕捉的有效性和精确性。模拟了油滴在水中上升、变形,与油层融合的过程,研究了气泡在产生、发展、脱离阶段的变形机理,和表面张力系数对气泡形状的影响。计算结果发现,表面张力系数越大,气泡在发展阶段持续的时间就越长,膨胀的程度也越大,并在脱离时刻,气泡的体积也越大,为进一步研究两相界面迁移特性提供了新的途径。     

三维气泡与刚性壁面的相互作用研究

基于势流理论建立水下爆炸气泡运动三维模型,采用边界积分法求解拉普拉斯方程,得到气泡的变形及位置,并在计算过程中引入弹性网格技术,避免了因网格扭曲而导致的数值发散,进而模拟了刚性壁面附近三维气泡的动态特性。在数值模拟过程中,将本文计算值与实验数据进行对比分析,结果表明,计算值与实验数据吻合良好。在此基础上,分别模拟了弱浮力、强Bjerknes力,强浮力、弱Bjerknes力以及浮力与Bjerknes力相当时壁面附近气泡的运动特征,并将各种工况的计算结果与基于开尔文冲量理论(Kelvin Impulse)的Blake准则进行对比分析与讨论,得到了不同参数下气泡的运动特征。     

水下爆炸气泡射流载荷计算的新方法研究

由于存在强冲击、多相界面复杂运动等强非线性效应,目前对于水下爆炸气泡运动的计算方法无法给出较为可信的水射流运动特征及其载荷形式.本文基于多相可压缩流体的five-equation计算模型,引入界面函数和界面密度压缩技术,采用5阶WENO重构与HLLC近似Riemann求解器进行空间离散,时间离散采用3阶TVD Runge-Kutta法.通过上述方法来提高水下爆炸气泡溃灭过程中气-液界面的计算精度,捕捉水射流的产生、发展以及水锤冲击等典型演化过程.计算结果表明,新的计算方法能够对水下近壁面爆炸气泡的溃灭过程进行有效的模拟,初步揭示了水射流直接载荷的特征,为水射流的产生机理及其损伤效应的研究提供技术支撑.     

近边界三维水下爆炸气泡动态特性研究

模拟了近壁、近自由面的水下爆炸气泡的非线性动态特性,假定水下爆炸气泡脉动阶段的流场是无旋、不可压缩的,采用高阶曲面三角形单元离散三维气泡表面,用边界积分法求解气泡的运动,在计算奇异积分时通过重新构造双层位势的主值积分消除双层奇异积分的奇异性,得到更精确的结果,并通过合理的加权方法精确的求解边界面上各节点的真实速度,结合弹性网格技术（elastic mesh technique, EMT）得到优化速度,在整个模拟过程中不需要采用数值光顺。将本文的三维模型与轴对称模型进行的对比分析表明,两种模型计算结果吻合很好,并用三维模型模拟了气泡与自由表面及圆筒的相互作用,水下爆炸气泡在自由表面及圆筒的联合作用下呈现出强非线性。     

简单Green函数法模拟三维水下爆炸气泡运动

假定水下爆炸气泡脉动阶段的流场是无旋、不可压缩的，运用势流理论导出气泡边界面运动的控制方程，采用高阶曲面三角形单元离散了维气泡表面，用边界积分法求解气泡的运动．并将计算结果与Rayleigh-Plesset气泡模型和试验数据进行对比分析，分析结果表明高阶曲面单元能够高精度的模拟水下爆炸气泡运动，且比线性单元有多方面的优越性．分别模拟了有、无重力场和刚壁时对气泡运动的影响，并预测了气泡在流场中膨胀、坍塌、迁移、射流形成等苇要动力学行为，同时建立了水下爆炸气泡与圆柱简相互作用的三维模型，模拟了自由液面、圆柱筒附近三维气泡的动力学特性．     

复杂边界附近气泡的动态特性研究

本文假设气泡周围流场为无粘、无旋、不可压缩的理想流体,基于势流理论,运用边界元法模拟近边界水下爆炸气泡的动力学特征,建立气泡、壁面以及自由面三者之间复杂耦合动力学模型,在气泡运动模拟过程中引入数值光顺技术及弹性网格技术(EMT),避免因网格扭曲而导致的数值发散,开发相应的三维计算程序,并与自由表面附近气泡运动的实验数据进行对比分析,计算结果表明本文的计算方法及程序具有较高的精度.在此基础上,用本文开发的三维程序模拟了单个及多个气泡与自由表面及圆筒等复杂边界的相互作用,其中包括水面漂浮结构和水下结构,气泡在自由表面及结构的联合作用下呈现出强非线性.本主文旨在为相关复杂边界附近气泡动力学特性研究提供参考.     

A large deformation framework for compressible viscoelastic materials: Constitutive equations and finite element implementation

For modeling the constitutive properties of viscoelastic solids in the context of small deformations, the so-called three-parameter solid is often used. The differential equation governing the model response may be derived in a thermodynamically consistent way considering linear spring-dashpot elements. The main problem in generalizing constitutive models from small to finite deformations is to extend the theory in a thermodynamically consistent way, so that the second law of thermodynamics remains satisfied in every admissible process. This paper concerns with the formulation and constitutive equations of finite strain viscoelastic material using multiplicative decomposition in a thermodynamically consistent manner. Based on the proposed constitutive equations, a finite element (FE) procedure is developed and implemented in an FE code. Subsequently, the code is used to predict the response of elastomer bushings. The finite element analysis predicts displacements and rotations at the relaxed state reasonably well. The response to coupled radial and torsional deformations is also simulated.     

梯度塑性的有限元分析及应变局部化模拟

对梯度塑性连续体提出了一个有限元方法．内状态变量的Ｌａｐｌａｃｉａｎ的确定基于它在求积点邻域的最小二乘方多项式近似．具体地考虑了具有一点求积和Ｈｏｕｒｇｌａｓｓ控制特点的基于胡海昌－Ｗａｓｈｉｚｕ变分原理的混合应变元和单元平均意义下的ｖｏｎ－Ｍｉｓｅｓ屈服准则．解析地导出了梯度塑性下一致性单元切线刚度矩阵和速率本构方程的一致性积分算法．在所建议的非局部化途径中求积点的一致性条件在非局部化意义下逐点精确满足．数值例题表明所提出的非经典连续体的有限元方法求解应变局部化问题的有效性     

Macroscopical non-linear material model for ferroelectric materials inside a hybrid finite element formulation

A new approach for modeling hysteretic non-linear ferroelectric ceramics is presented, based on a fully ferroelectric/ferroelastic coupled macroscopic material model. The material behavior is described by a set of yield functions and the history dependence is stored in internal state variables representing the remanent polarization and the remanent strain. For the solution of the electromechanical coupled boundary value problem, a hybrid finite element formulation is used. Inside this formulation the electric displacement is available as nodal quantity (i.e. degree of freedom) which is used instead of the electric field to determine the evolution of remanent polarization. This involves naturally the electromechanical coupling. A highly efficient integration technique of the constitutive equations, defining a system of ordinary differential equations, is obtained by a customized return mapping algorithm. Due to some simplifications of the algorithm, an analytical solution can be calculated. The automatic differentiation technique is used to obtain the consistent tangent operator. Altogether this has been implemented into the finite element code FEAP via a user element. Extensive verification tests are performed in this work to evaluate the behavior of the material model under pure electrical and mechanical as well as coupled and multi-axial loading conditions.     

高温条件下混凝土墙体的非线性分析研究

本文建立了一个弹塑性-损伤耦合本构模型用于数值模拟高温下混凝土的真实破坏过程。导出了一个利用Newton-Raphson迭代的一般的直接应力返回映射算法。同时求解应力向量和塑性、损伤的内状态变量。并推导了用于化学-热-湿-力学耦合分析的全局守恒方程Newton-Raphson迭代过程的一致性切线模量矩阵。建议了一个用于弹塑性-损伤耦合分析的两级求解过程。给出的数值例题结果显示了所提出的数法和公式的正确性,表明了所发展的弹塑性-损伤耦合本构模型在模拟高温下混凝土墙体中复杂破坏过程的能力。     

Thermo-mechanically coupled model of diffusionless phase transformation in austenitic steel

A constitutive model of thermo-mechanically coupled finite strain plasticity considering martensitic phase transformation is presented. The model is formulated within a thermodynamic framework, giving a physically sound format where the thermodynamic mechanical and chemical forces that drive the phase transformation are conveniently identifiable. The phase fraction is treated through an internal variable approach and the first law of thermodynamics allows a consistent treatment of the internal heat generation due to dissipation of inelastic work. The model is calibrated against experimental data on a Ni–Cr steel of AISI304-type, allowing illustrative simulations to be performed. It becomes clear that the thermal effects considered in the present formulation have a significant impact on the material behavior. This is seen, not least, in the effects found on forming limit diagrams, also considered in the present paper.     

Lateral buckling of beams with shear deformations – A hyperelastic formulation

The equilibrium and buckling equations are derived for the lateral buckling of a prismatic straight beam. A consistent finite strain constitutive law is used, which is based on a hyperelastic model for an isotropic material. The kinematics of the cross-sectional deformations are based on a Timoshenko type beam displacement of the cross-sectional plane using Euler angles and two shear finite rotations coupled with warping taken normal to the displaced plane. Also derived are the second order approximations to the displacements, curvatures, twist and internal actions. The constitutive relationships for the internal actions reveal new coupling terms between the bending moments, torsion and bimoment, which are functions of the cross-sectional warping and shear deformations. New Wagner type nonlinear torsion terms are derived which are functions of the warping of the cross-sectional plane, and are coupled to the twisting and shear deformations of the cross-section. Solutions are determined for the lateral buckling of a prismatic monosymmetric beam under pure bending and the flexural–torsional buckling under axial compression. For the flexural–torsional buckling problem it is found that the Euler type column buckling formula is consistent with Haringx’s column buckling formula while the torsional buckling formula is different to conventional equations. The second variation of the total potential is also derived. The effects of shear deformations are explored by examining the non-dimensional lateral buckling equation for a simply supported beam.     

On internal dissipation inequalities and finite strain inelastic constitutive laws: Theoretical and numerical comparisons

Internal dissipation always occurs in irreversible inelastic deformation processes of materials. The internal dissipation inequalities (specific mathematical forms of the second law of thermodynamics) determine the evolution direction of inelastic processes. Based on different internal dissipation inequalities several finite strain inelastic constitutive laws have been formulated for instance by Simo [Simo, J.C., 1992. Algorithms for static and dynamic multiplicative plasticity that preserve the classical return mapping schemes of the infinitesimal theory. Computer Methods in Applied Mechanics and Engineering 99, 61–112]; Simo and Miehe [Simo, J.C., Miehe, C., 1992. Associative coupled thermoplasticity at finite strains: formulation, numerical analysis and implementation. Computer Methods in Applied Mechanics and Engineering 98, 41–104]; Lion [Lion, A., 1997. A physically based method to represent the thermo-mechanical behavior of elastomers. Acta Mechanica 123, 1–25]; Reese and Govindjee [Reese, S., Govindjee, S., 1998. A theory of finite viscoelasticity and numerical aspects. International Journal of Solids and Structures 35, 3455–3482]; Lin and Schomburg [Lin, R.C., Schomburg, U., 2003. A finite elastic–viscoelastic–elastoplastic material law with damage: theoretical and numerical aspects. Computer Methods in Applied Mechanics and Engineering 192, 1591–1627]; Lin and Brocks [Lin, R.C., Brocks, W., 2004. On a finite strain viscoplastic theory based on a new internal dissipation inequality. International Journal of Plasticity 20, 1281–1311]; and Lin and Brocks [Lin, R.C., Brocks, W., 2005. An extended Chaboche’s viscoplastic law at finite strains: theoretical and numerical aspects. Journal of Materials Science and Technology 21, 145–147]. These constitutive laws are consistent with the second law of thermodynamics. As the internal dissipation inequalities are described in different configurations or coordinate systems, the related constitutive laws are also formulated in the corresponding configurations or coordinate systems. Mathematically, these constitutive laws have very different formulations. Now, a question is whether the constitutive laws provide identical constitutive responses for the same inelastic constitutive problems. In the present work, four types of finite strain viscoelastic and endochronically plastic laws as well as three types of J₂-plasticity laws are formulated based on four types of dissipation inequalities. Then, they are numerically compared for several problems of homogeneous or complex finite deformations. It is demonstrated that for the same inelastic constitutive problem the stress responses are identical for deformation processes without rotations. In the simple shear deformation process with large rotation, the presented viscoelastic and endochronically plastic laws also show almost identical stress responses up to a shear strain of about 100%. The three laws of J₂-plasticity also produce the same shear stresses up to a shear strain of 100%, while different normal stresses are generated even at infinitesimal shear strains. The three J₂-plasticity laws are also compared at three complex finite deformation processes: billet upsetting, cylinder necking and channel forming. For the first two deformation processes similar constitutive responses are obtained, whereas for the third deformation process (with large global rotations) significant differences of constitutive responses can be observed.     

Simulation of damage evolution in ductile metals undergoing dynamic loading conditions

A set of constitutive equations for large rate-dependent elastic-plastic-damage materials at elevated temperatures is presented to be able to analyze adiabatic high strain rate deformation processes for a wide range of stress triaxialities. The model is based on the concepts of continuum damage mechanics. Since the material macroscopic thermo-mechanical response under large strain and high strain rate deformation loading is governed by different physical mechanisms, a multi-dissipative approach is proposed. It incorporates thermo-mechanical coupling effects as well as internal dissipative mechanisms through rate-dependent constitutive relations with a set of internal variables. In addition, the effect of stress triaxiality on the onset and evolution of plastic flow, damage and failure is discussed.Furthermore, the algorithm for numerical integration of the coupled constitutive rate equations is presented. It relies on operator split methodology resulting in an inelastic predictor-elastic corrector technique. The explicit finite element program LS-DYNA augmented by an user-defined material subroutine is used to approximate boundary-value problems under dynamic loading conditions. Numerical simulations of dynamic experiments with different specimens are performed and good correlation of numerical results and published experimental data is achieved. Based on numerical studies modified specimens geometries are proposed to be able to detect complex damage and failure mechanisms in Hopkinson-Bar experiments.     

Amplitude dependence of filler-reinforced rubber: Experiments,constitutive modelling and FEM – Implementation

The focus of the present paper is the experimental investigation, the constitutive representation and the numerical simulation of the amplitude dependence of filler-reinforced elastomers. A standard way to investigate the dynamic properties of viscoelastic materials is via the dynamic modulus which is obtained from stress signals due to harmonic strain excitations. Based on comprehensive experimental data, an amplitude-dependent constitutive model of finite viscoelasticity is developed. The model is based on a modified Maxwell chain with process-dependent viscosities which depend on additional internal state variables. The evaluation of this thermodynamically consistent model is possible in both the time domain, via stress-time signals, and in the frequency domain, via the dynamic modulus. This property is very profitable for the parameter identification process. The implementation of the constitutive model into the commercial finite element code ANSYS with the user-programmable feature (UPF) USERMAT for large deformations in updated Lagrange formulation is presented. This implementation allows simulating the time-dependent behaviour of rubber components under arbitrary transient loading histories. Due to physical and geometrical nonlinearities, these simulations are not possible in the frequency domain. But, transient FEM computations of large loading histories are sometimes not possible in an acceptable time. In the context of the parameter identification the fundamental ideas are presented, how this problem has been solved. Transient FEM simulations of real rubber components are also shown to visualize the properties of the model in the context of the transient material behaviour.     

A 3D finite strain phenomenological constitutive model for shape memory alloys considering martensite reorientation

Most devices based on shape memory alloys experience both finite deformations and non-proportional loading conditions in engineering applications. This motivates the development of constitutive models considering finite strain as well as martensite variant reorientation. To this end, in the present article, based on the principles of continuum thermodynamics with internal variables, a three-dimensional finite strain phenomenological constitutive model is proposed taking its basis from the recent model in the small strain regime proposed by Panico and Brinson (J Mech Phys Solids 55:2491–2511, 2007). In the finite strain constitutive model derivation, a multiplicative decomposition of the deformation gradient into elastic and inelastic parts, together with an additive decomposition of the inelastic strain rate tensor into transformation and reorientation parts is adopted. Moreover, it is shown that, when linearized, the proposed model reduces exactly to the original small strain model.