A variational multiscale Newton–Schur approach for the incompressible Navier–Stokes equations

In the following paper, we present a consistent Newton–Schur (NS) solution approach for variational multiscale formulations of the time‐dependent Navier–Stokes equations in three dimensions. The main contributions of this work are a systematic study of the variational multiscale method for three‐dimensional problems and an implementation of a consistent formulation suitable for large problems with high nonlinearity, unstructured meshes, and non‐symmetric matrices. In addition to the quadratic convergence characteristics of a Newton–Raphson‐based scheme, the NS approach increases computational efficiency and parallel scalability by implementing the tangent stiffness matrix in Schur complement form. As a result, more computations are performed at the element level. Using a variational multiscale framework, we construct a two‐level approach to stabilizing the incompressible Navier–Stokes equations based on a coarse and fine‐scale subproblem. We then derive the Schur complement form of the consistent tangent matrix. We demonstrate the performance of the method for a number of three‐dimensional problems for Reynolds number up to 1000 including steady and time‐dependent flows. Copyright © 2009 John Wiley & Sons, Ltd.     

An implicit consistent algorithm for the integration of thermoviscoplastic constitutive equations in adiabatic conditions and finite deformations

The so-called viscoplastic consistency model, proposed by Wang, Sluys and de Borst, is extended here to the integration of a thermoviscoplastic constitutive equation for J₂ plasticity and adiabatic conditions. The consistency condition in this case includes not only strain rate but also the effect of temperature on the yield function. Using the backward Euler integration scheme to integrate the constitutive equations, an implicit algorithm is proposed, leading to generalized expressions of the classical return mapping algorithm for J₂ plasticity, both for the iterative calculation of the plastic multiplier increment and for the consistent tangent operator when strain rate and temperature are considered also as state variables of the hardening equation. The model was implemented in a commercial finite element code and its performance is demonstrated with the numerical simulation of four Taylor impact tests.     

一种非统一黏塑性本构的有限元实现

在与率相关黏塑性本构模型的背应力分量中加入静态回复项通常用来更精确地描述材料应力松弛效应．本文基于包含静态回复项的非统一黏塑性本构模型,提出了对应的隐式应力积分算法．该模型可简化为两个关于偏背应力和等效应力张量的方程组,通过牛顿迭代法进行求解．在应力积分算法中采用了Return-Mapping求解策略．此外还推导出了材料的一致切线刚度矩阵来判断有限元分析收敛性．在此基础上开发出了适用于ABAQUS工程软件的UMAT用户材料子程序．文中选取了高温下镍合金材料在伴随应力松弛的循环加载路径下的试验结果来验证算法的有效性．从比较结果中可以发现：用户子程序即使在大时间增量下也具有良好的收敛性,静态回复项的引入使得模型对应力松弛效应的模拟更加准确．     

Homogenization of two-phase elasto-plastic composite materials and structures: Study of tangent operators,cyclic plasticity and numerical algorithms

We develop homogenization schemes and numerical algorithms for two-phase elasto-plastic composite materials and structures. A Hill-type incremental formulation enables the simulation of unloading and cyclic loadings. It also allows to handle any rate-independent model for each phase. We study the crucial issue of tangent operators: elasto-plastic (or “continuum”) versus algorithmic (or “consistent”), and anisotropic versus isotropic. We apply two methods of extraction of isotropic tangent moduli. We compare mathematically the stiffnesses of various tangent operators. All rate equations are discretized in time using implicit integration. We implemented two homogenization schemes: Mori–Tanaka and a double inclusion model, and two plasticity models: classical J₂ plasticity and Chaboche’s model with non-linear kinematic and isotropic hardenings. We consider composites with different properties and present several discriminating numerical simulations. In many cases, the results are validated against finite element (FE) or experimental data. We integrated our homogenization code into the FE program ABAQUS using a user material interface UMAT. A two-scale procedure allows to compute realistic structures made of non-linear composite materials within reasonable CPU time and memory usage; examples are shown.     

牛顿迭代一致性算法及其在板弹塑性有限元分析中的应用

本文简略讨论了有限载荷增量弹塑性有限元分析中传统切线刚度法丧失精度和牛顿迭代平方收敛速度的原因,并提出保持牛顿迭代平方收敛速度、保证一阶精度和无条件稳定性的一致性算法.一致性算法具备以下两个特征:1)采用路径无关计算格式;2)采用一致弹塑性切线模量。根据一致性算法构造出以弯矩和曲率为基本变量的弹塑性板弯曲有限元NIDKQ元。数值结果表明NIDKQ元具有令人满意的精度,同时验证了有限载荷增量下牛顿迭代一致性算法的平方收敛率特性,而传统切线刚度法随着塑性区的扩展将大大降低收敛速度。     

Mathematical and numerical modeling of the non-associated plasticity of soils—Part 2: Finite element analysis

The method of the implicit standard material has allowed the formulation of a consistent mathematical model of the boundary value problem for the non-associated plasticity of soil. The mean accomplished steps are the achievement of the bipotential function, the recovering of the stress–strain relationship under a normality rule, introduction of the bifunctional and the proof of the solution existence. Here the mathematical model is discretized by the finite element method. First, the stress update scheme was formulated, the tangent matrix is explicitly derived and then the non-linear system is solved by the Newton–Raphson method where a new algorithm using a symmetrical tangent matrix is improved. This is in opposition to conventional non-associated plasticity, which uses a non-symmetric tangent matrix. Through the numerical examples we show the feasibility and the efficiency of the algorithm. It is also seen that we perform some studies of the numerical solutions, particularly the comparison between associated and non-associated limit load.     

Coupled thermoviscoplasticity of glassy polymers in the logarithmic strain space based on the free volume theory

The paper outlines a constitutive model for finite thermo-visco-plastic behavior of amorphous glassy polymers and considers details of its numerical implementation. In contrast to existing kinematical approaches to finite plasticity of glassy polymers, the formulation applies a plastic metric theory based on an additive split of Lagrangian Hencky-type strains into elastic and plastic parts. The analogy between the proposed formulation in the logarithmic strain space and the geometrically linear theory of plasticity, makes this constitutive framework very transparent and attractive with regard to its numerical formulation. The characteristic strain hardening of the model is derived from a polymer network model. We consider the particularly simple eight chain model, but also comment on the recently developed microsphere model. The viscoplastic flow rule in the logarithmic strain space uses structures of the free volume flow theory, which provides a highly predictive modeling capacity at the onset of viscoplastic flow. The integration of this micromechanically motivated approach into a three-dimensional computational model is a key concern of this work. We outline details of the numerical implementation of this model, including elements such as geometric pre- and post-transformations to/from the logarithmic strain space, a thermomechanical operator split algorithm consisting of an isothermal mechanical predictor followed by a heat conduction corrector and finally, the consistent linearization of the local update algorithm for the dissipative variables as well as its relationship to the global tangent operator. The performance of the proposed formulation is demonstrated by means of a spectrum of numerical examples, which we compare with our experimental findings.     

牛顿迭代一致性算法及其在板弹塑性有限元分析中的应用

本文简略讨论了有限载荷增量弹塑性有限元分析中传统切线刚度法丧失精度和牛顿迭代平方收敛速度的原因,并提出保持牛顿迭代平方收敛速度、保证一阶精度和无条件稳定性的一致性算法.一致性算法具备以下两个特征:1)采用路径无关计算格式;2)采用一致弹塑性切线模量。根据一致性算法构造出以弯矩和曲率为基本变量的弹塑性板弯曲有限元NIDKQ元。数值结果表明NIDKQ元具有令人满意的精度,同时验证了有限载荷增量下牛顿迭代一致性算法的平方收敛率特性,而传统切线刚度法随着塑性区的扩展将大大降低收敛速度。     

An eigenanalysis-based bifurcation indicator proposed in the framework of a reduced-order modeling technique for non-linear structural analysis

The Koiter–Newton method is a reduced order modeling technique which allows us to trace efficiently the entire equilibrium path of a non-linear structural analysis. In the framework of buckling the method is capable to handle snap-back and snap-through phenomena but may fail to predict reliably bifurcation branches along the equilibrium path. In this contribution we extend the original Koiter–Newton approach with a reliable and accurate bifurcation indicator which is based on an eigenanalysis of the reduced order tangent stiffness matrix. The proposed indicator has a negligible numerical effort since all computations refer to the reduced order model which is typically of very small dimension. The extension allows the identification of bifurcation points and a tracing of corresponding bifurcation branches in each sector of the equilibrium path. The performance of the method in terms of reliability, accuracy and computational effort is demonstrated with several examples.     

Accurate iteration-free mixed-stabilised formulation for laminar incompressible Navier–Stokes: Applications to fluid–structure interaction

Stabilised mixed velocity–pressure formulations are one of the widely-used finite element schemes for computing the numerical solutions of laminar incompressible Navier–Stokes. In these formulations, the Newton–Raphson scheme is employed to solve the nonlinearity in the convection term. One fundamental issue with this approach is the computational cost incurred in the Newton–Raphson iterations at every load/time step. In this paper, we present an iteration-free mixed finite element formulation for incompressible Navier–Stokes that preserves second-order temporal accuracy of the generalised-alpha and related schemes for both velocity and pressure fields. First, we demonstrate the second-order temporal accuracy using numerical convergence studies for an example with a manufactured solution. Later, we assess the accuracy and the computational benefits of the proposed scheme by studying the benchmark example of flow past a fixed circular cylinder. Towards showcasing the applicability of the proposed technique in a wider context, the inf–sup stable P2–P1 pair for the formulation without stabilisation is also considered. Finally, the resulting benefits of using the proposed scheme for fluid–structure interaction problems are illustrated using two benchmark examples in fluid-flexible structure interaction.     

Multiscale modeling of irradiated polycrystalline FCC metals

We propose a set of models for the post-irradiation deformation response of polycrystalline FCC metals. First, a defect- and dislocation-density based evolution model is developed to capture the features of irradiation-induced hardening as well as intra-granular softening. The proposed hardening model is incorporated within a rate-independent single crystal plasticity model. The result is a non-homogeneous deformation model that accounts for defect absorption on the active slip planes during plastic loading. The macroscopic non-linear constitutive response of the polycrystalline aggregate of the single crystal grains is then obtained using a micro–macro transition scheme, which is realized within a Jacobian-free multiscale method (JFMM). The Jacobian-free approach circumvents explicit computation of the tangent matrix at the macroscale by using a Newton–Krylov process. This has a major advantage in terms of storage requirements and computational cost over existing approaches based on homogenized material coefficients in which explicit Jacobian computation is required at every Newton step. The mechanical response of neutron-irradiated single and polycrystalline OFHC copper is studied and it is shown to capture experimentally observed grain-level phenomena.     

Finite viscoplasticity of amorphous glassy polymers in the logarithmic strain space

The paper outlines a new constitutive model and experimental results of rate-dependent finite elastic–plastic behavior of amorphous glassy polymers. In contrast to existing kinematical approaches to finite viscoplasticity of glassy polymers, the formulation proposed is constructed in the logarithmic strain space and related to a six-dimensional plastic metric. Therefore, it a priori avoids difficulties concerning with the uniqueness of a plastic rotation. The constitutive framework consists of three major steps: (i) A geometric pre-processing defines a total and a plastic logarithmic strain measures determined from the current and plastic metrics, respectively. (ii) The constitutive model describes the stresses and the consistent moduli work-conjugate to the logarithmic strain measures in an analogous structure to the geometrically linear theory. (iii) A geometric post-processing maps the stresses and the algorithmic tangent moduli computed in the logarithmic strain space to their nominal, material or spatial counterparts in the finite deformation space. The analogy between the formulation of finite plasticity in the logarithmic strain space and the geometrically linear theory of plasticity makes this framework very attractive, in particular regarding the algorithmic implementation. The flow rule for viscoplastic strains in the logarithmic strain space is adopted from the celebrated double-kink theory. The post-yield kinematic hardening is modeled by different network models. Here, we compare the response of the eight chain model with the newly proposed non-affine micro-sphere model. Apart from the constitutive model, experimental results obtained from both the homogeneous compression and inhomogeneous tension tests on polycarbonate are presented. Besides the load–displacement data acquired from inhomogeneous experiments, quantitative three-dimensional optical measurements of the surface strain fields are carried out. With regard to these experimental data, the excellent predictive quality of the theory proposed is demonstrated by means of representative numerical simulations.     

Stress resultants plasticity with general closest point projection

In this paper, general closest point projection algorithm is derived for the elastoplastic behavior of a cross-section of a beam finite element. For given section deformations, the section forces (stress resultants) and the section tangent stiffness matrix are obtained as the response for the cross-section. Backward Euler time integration rule is used for the solution of the nonlinear evolution equations. The solution yields the general closest projection algorithm for stress resultants plasticity model. Algorithmic consistent tangent stiffness matrix for the section is derived. Numerical verification of the algorithms in a mixed formulation beam finite element proves the accuracy and robustness of the approach in simulating nonlinear behavior.     

Determination of sensor positions for predictive maintenance of revolving machines

The monitoring by measurement and analysis of vibration is largely used to detect the defects in revolving machines. The determination of the best sensor positions is one of the main research goals in the field of predictive maintenance. This paper proposes a numerical methodology based on a finite element model and a spectral analysis in order to find optimum sensor positions. The bearing is a key component for the vibration propagation from the moving parts to static ones. An analytical bearing model and its numerical implementation in a finite element code are presented. The tangent stiffness matrix of the bearing element is obtained by the Newton–Raphson method and then used for the modal and spectral analyses. Several techniques are used to find the most sensitive zones to common defects. The proposed numerical approach correlate well with the experimental results. The numerical modeling of a grinder shows the interests in industrial applications.     

Effect of strain hardening in elastic–plastic transition behavior in a hemisphere in contact with a rigid flat

An axisymmetrical hemispherical asperity in contact with a rigid flat is modeled for an elastic–plastic material on the lines of the Kogut–Etsion Model (KE Model) and the Jackson–Green Model (JG Model). The present work extends the previous KE and JG works, accounting for the effect of realistic material behavior in terms of the varying yield strengths and the isotropic strain hardening behavior. The predicted results show that the transition behavior of the materials from the elastic–plastic to the fully plastic case is influenced by the yield strength and the tangent modulus (E_t) and such transition do not take place at specific values of interference ratios as suggested by the KE model. New empirical relations are proposed to determine the contact load and the contact area based on the analysis. Numerical results from the finite element modeling are also validated with an experimental ball on flat configuration approach.     

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

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

A modified upper bound approach to limit analysis for plane strain deeply cracked specimens

The classical upper bound approach of limit analysis is based on assumption of rigid blocks of deformation that move between lines of tangential displacement discontinuity. This assumption leads to considerable simplification but often at cost of higher estimate of the actual load. Moreover, in many cases, it does not give a correct shape of the plastic field. In order to overcome these limitations a modified upper bound approach is proposed in this article. The proposed approach is basically an energetic approach but unlike the classical upper bound approach it is capable of including presence of statically governed stress field. As an application, of proposed approach, theoretical plane strain solutions are presented for deeply cracked fracture mechanics specimens (single edge cracked specimen in pure bending – SE (PB), single edge cracked specimen in three-point bending – SE (B), and compact tension – C (T) specimens). Plane strain plasticity problem in rigid elastic–plastic mono-material (homogeneous) was solved to evaluate useful parameters like limit load, plastic eta function (η_p) and plastic rotation factor (r_p) and in bi-material (mismatch welds) to evaluate mismatch limit load, for deeply cracked specimens. New kinematically admissible velocity fields are proposed for SE (B) and C (T) specimens. Proposed theoretical solutions were confirmed by classical slip-line field solutions, wherever available, and by detailed elastic–plastic finite element analysis with Von-Mises yield criterion. Good agreement was found between proposed solutions and results obtained from the classical slip-line field theory and finite element analysis.     

A C2 continuous approximation to the Mohr–Coulomb yield surface

In spite of the development of more sophisticated constitutive models for soil, the Mohr–Coulomb yield criterion remains a popular choice for geotechnical analysis due to its simplicity and ease of use by practising engineers. The implementation of the criterion in finite element programs, however, presents some numerical difficulties due to the gradient discontinuities which occur at both the edges and the tip of the hexagonal yield surface pyramid. Furthermore, some implicit techniques utilising consistent tangent stiffness formulations are unable to achieve full quadratic convergence as the yield criteria is not C2 continuous. This paper extends the previous work of Abbo and Sloan (1995) through the introduction of C2 continuous rounding of the Mohr–Coulomb yield surface in the octahedral plane. This approximation, when combined with the hyperbolic approximation in the meridional plane (Abbo and Sloan, 1995), describes a yield surface that is C2 continuous at all stress states. The new smooth yield surface can be made to approximate the Mohr–Coulomb yield function as closely as required by adjusting only two parameters, and is suitable for consistent tangent stiffness formulations.