Bifurcation points and bifurcated branches by an asymptotic numerical method and Padé approximants

The aim of this work is to develop a reliable and fast algorithm to compute bifurcation points and bifurcated branches. It is based upon the asymptotic numerical method (ANM) and Padé approximants. The bifurcation point is detected by analysing the poles of Padé approximants or by evaluating, along the computed solution branch, a bifurcation indicator well adapted to ANM. Several examples are presented to assess the effectiveness of the proposed method, that emanate from buckling problems of thin elastic shells. Especially problems involving large rotations are discussed. Copyright © 2004 John Wiley & Sons, Ltd.     

An adaptive strategy for the bivariate solution of finite element problems using multivariate nested Padé approximants

Most engineering applications involving solutions by numerical methods are dependent on several parameters, whose impact on the solution may significantly vary from one to the other. At times an evaluation of these multivariate solutions may be required at the expense of a prohibitively high computational cost. In the present paper, an adaptive approach is proposed as a way to estimate the solution of such multivariate finite element problems. It is based upon the integration of so‐called nested Padé approximants within the finite element procedure. This procedure includes an effective control of the approximation error, which enables adaptive refinements of the converged intervals upon reconstruction of the solution. The main advantages lie in a potential reduction of the computational effort and the fact that the level of a priori knowledge required about the solution in order to have an accurate, efficient, and well‐sampled estimate of the solution is low. The approach is introduced for bivariate problems, for which it is validated on elasto‐poro‐acoustic problems of both academic and more industrial scale. It is argued that the methodology in general holds for more than two variables, and a discussion is opened about the truncation refinements required in order to generalize the results accordingly. Copyright © 2014 John Wiley & Sons, Ltd.     

Padé approximants and the modal connection: Towards increased robustness for fast parametric sweeps

To increase the robustness of a Padé‐based approximation of parametric solutions to finite element problems, an a priori estimate of the poles is proposed. The resulting original approach is shown to allow for a straightforward, efficient, subsequent Padé‐based expansion of the solution vector components, overcoming some of the current convergence and robustness limitations. In particular, this enables for the intervals of approximation to be chosen a priori in direct connection with a given choice of Padé approximants. The choice of these approximants, as shown in the present work, is theoretically supported by the Montessus de Ballore theorem, concerning the convergence of a series of approximants with fixed denominator degrees. Key features and originality of the proposed approach are (1) a component‐wise expansion which allows to specifically target subsets of the solution field and (2) the a priori, simultaneous choice of the Padé approximants and their associated interval of convergence for an effective and more robust approximation. An academic acoustic case study, a structural‐acoustic application, and a larger acoustic problem are presented to demonstrate the potential of the approach proposed.     

Studies of elastic and elastic–plastic J integral for mixed mode cracked plate under biaxial loading

Analyses of I–II mixed mode central cracked plate by finite element method are performed in this paper, and some different phenomena are found. First for I–II mixed mode crack, the distribution of J integral along crack tip thickness depends on biaxiality factors because of the existence of vertex (corner) singularity, which is unlike that for mode I or mode II crack. Then J integrals at middle layer keep constant for any cracked plates with different inclined angles β when the biaxiality ratio is equal to 1 or ?1, which implies that the inclined angle or the extent of I–II mixed mode has no effect on the J integral for positive or negative equal axial loading conditions. And the decreasing trend of J integral with the inclined angle β for biaxiality ratio λ being between?1 and 1 is just opposite with that for biaxiality ratio λ being larger than 1 and smaller than ?1. Finally, proposed h₁ (a/W, n, λ, β) of cracked plate with different inclined angles under different biaxial loading are calculated.     

Ductile damage analysis of elasto‐plastic shells at large inelastic strains

The objective of this contribution is to model ductile damage phenomena under consideration of large inelastic strains, to couple the corresponding constitutive law with a multi‐layer shell kinematics and to give finally an adequate finite element formulation. An elastic–plastic constitutive law is formulated by using a spatial hyperelasto‐plastic formulation based on the multiplicative decomposition of the deformation gradient. To include isotropic ductile damage the continuum damage model of Rousselier is modified so as to consider large strains and additionally extended by various void nucleation and macro‐crack criteria. In order to achieve numerical efficiency, elastic strains are supposed to be sufficiently small providing a numerical effective integration based on the backward Euler rule. Finite element formulation is enriched by means of the enhanced strain concept. Thus the well‐known deficiencies due to incompressible deformations and the inclusion of transverse strains are avoided. Several examples are given to demonstrate the performance of the algorithms developed concerning large inelastic strains of shells and ductile damage phenomena. Copyright © 2000 John Wiley & Sons, Ltd.     

Higher-order and higher floating-point precision numerical approximations of finite strain elasticity moduli

Two real-domain numerical approximation methods for accurate computation of finite strain elasticity moduli are developed and their accuracy and computational efficiency are investigated, with reference to hyperelastic constitutive models with known analytical solutions. The methods are higher-order and higher floating-point precision numerical approximation, the latter being novel in this context. A general formula for higher-order approximation finite difference schemes is derived and a new procedure is proposed to implement increased floating-point precision. The accuracy of the approximated elasticity moduli is investigated numerically using higher-order approximations in standard double precision and increased quadruple precision. It is found that, as the order of the approximation increases, the elasticity moduli tend toward the analytical solution. Using higher floating-point precision, the approximated elasticity moduli for all orders of approximation are found to be more accurate than the standard double precision evaluation of the analytical moduli. Application of the techniques to a finite element problem shows that the numerically approximated methods obtain convergence equivalent to the analytical method but require greater computational effort. It is concluded that numerical approximation of elasticity moduli is a powerful and effective means of implementing advanced constitutive models in the finite element method without prior derivation of difficult analytical solutions.     

A boundary condition in Padé series for frequency‐domain solution of wave propagation in unbounded domains

A boundary condition satisfying the radiation condition at infinity is frequently required in the numerical simulation of wave propagation in an unbounded domain. In a frequency domain analysis using finite elements, this boundary condition can be represented by the dynamic stiffness matrix of the unbounded domain defined on its boundary. A method for determining a Padé series of the dynamic stiffness matrix is proposed in this paper. This method starts from the scaled boundary finite‐element equation, which is a system of ordinary differential equations obtained by discretizing the boundary only. The coefficients of the Padé series are obtained directly from the ordinary differential equations, which are not actually solved for the dynamic stiffness matrix. The high rate of convergence of the Padé series with increasing order is demonstrated numerically. This technique is applicable to scalar waves and elastic vector waves propagating in anisotropic unbounded domains of irregular geometry. It can be combined seamlessly with standard finite elements. Copyright © 2006 John Wiley & Sons, Ltd.     

Tensor‐based Gauss–Jacobi numerical integration for high‐order mass and stiffness matrices

In this work, we choose the points and weights of the Gauss–Jacobi, Gauss–Radau–Jacobi and Gauss–Lobatto–Jacobi quadrature rules that optimize the number of operations for the mass and stiffness matrices of the high‐order finite element method. The procedure is particularly applied to the mass and stiffness matrices using the tensor‐based nodal and modal shape functions given in (Int. J. Numer. Meth. Engng 2007; 71 (5):529–563). For square and hexahedron elements, we show that it is possible to use tensor product of the 1D mass and stiffness matrices for the Poisson and elasticity problem. For the triangular and tetrahedron elements, an analogous analysis given in (Int. J. Numer. Meth. Engng 2005; 63 (2):1530–1558) was considered for the selection of the optimal points and weights for the stiffness matrix coefficients for triangles and mass and stiffness matrices for tetrahedra. Copyright © 2009 John Wiley & Sons, Ltd.     

Anisotropic Gurson‐Tvergaard‐Needleman plasticity and damage model for finite element analysis of elastic‐plastic problems

Implementation and analysis of the anisotropic version of the Gurson‐Tvergaard‐Needleman (GTN) isotropic damage criterion are performed on the basis of Hill's quadratic anisotropic yield theory with the definition of an effective anisotropic coefficient to represent the elastic‐plastic behavior of ductile metals. This study aims to analyze the extension of the GTN model suitable for anisotropic porous metals and to investigate the GTN model extension. An anisotropic damage model is implemented using the user material subroutine in ABAQUS/standard finite element code. The implementation is verified and applied to simulate a uniaxial tensile test on a commercially produced aluminum sheet material for three‐dimensional and plane stress test cases. Spherical and ellipsoidal micro voids are considered in the matrix material, and their effects on the uniaxial stress‐strain response of the material are analyzed. Hill's quadratic anisotropic yield theory predicts substantially large damage evolution and a low stress‐strain curve compared with those predicted by the isotropic model. An approximate model for anisotropic materials is proposed to avoid increased damage evolution. In this approximate model, Hill's anisotropic constants are replaced with an effective anisotropy coefficient. All model‐generated stress‐strain predictions are compared with the experimental stress‐strain curve of AA6016‐T4 alloy.     

The mixed finite element method for the quasi‐static and dynamic analysis of viscoelastic timoshenko beams

The quasi‐static and dynamic responses of a linear viscoelastic Timoshenko beam on Winkler foundation are studied numerically by using the hybrid Laplace–Carson and finite element method. In this analysis the field equation for viscoelastic material is used. In the transformed Laplace–Carson space two new functionals have been constructed for viscoelastic Timoshenko beams through a systematic procedure based on the Gâteaux differential. These functionals have six and two independent variables respectively. Two mixed finite element formulations are obtained; TB12 and TB4. For the inverse transform Schapery and Fourier methods are used. The numerical results for quasi‐static and dynamic responses of several visco‐elastic models are presented. Copyright © 1999 John Wiley & Sons, Ltd.     

Design sensitivity analysis in large deformation elasto‐plastic and elasto‐viscoplastic problems

Implicit time integration algorithm derived by Simo for his large‐deformation elasto‐plastic constitutive model is generalized, for the case of isotropy and associative flow rule, towards viscoplastic material behaviour and consistently differentiated with respect to its input parameters. Combining it with the general formulation of design sensitivity analysis (DSA) for non‐linear finite element transient equilibrium problem, we come at a numerically efficient, closed‐form finite element formulation of DSA for large deformation elasto‐plastic and elasto‐viscoplastic problems, with various types of design variables (material constants, shape parameters). The paper handles several specific issues, like the use of a non‐algorithmic coefficient matrix or sensitivity discontinuities at points of instantaneous structural stiffness change. Computational examples demonstrate abilities of the formulation and quality of results. Copyright © 2005 John Wiley & Sons, Ltd.     

Scale transition and enforcement of RVE boundary conditions in second‐order computational homogenization

Formulation of the scale transition equations coupling the microscopic and macroscopic variables in the second‐order computational homogenization of heterogeneous materials and the enforcement of generalized boundary conditions for the representative volume element (RVE) are considered. The proposed formulation builds on current approaches by allowing any type of RVE boundary conditions (e.g. displacement, traction, periodic) and arbitrary shapes of RVE to be applied in a unified manner. The formulation offers a useful geometric interpretation for the assumptions associated with the microstructural displacement fluctuation field within the RVE, which is here extended to second‐order computational homogenization. A unified approach to the enforcement of the boundary conditions has been undertaken using multiple constraint projection matrices. The results of an illustrative shear layer model problem indicate that the displacement and traction RVE boundary conditions provide the upper and lower bounds of the response determined via second‐order computational homogenization, and the solution associated with the periodic RVE boundary conditions lies between them. Copyright © 2007 John Wiley & Sons, Ltd.     

Validity of Lagrange (Bézier) and rational Bézier quads of degree 2

Finite elements of degree two or more are needed to solve various PDE problems. This paper discusses a method to validate such meshes for the case of quadrilateral elements of degree 2. The first section of this paper comes back to Bézier curve and Bézier quadrilateral patches of degree 2. The way in which a Bézier quad patch and a Q2 finite element quad are related is introduced. The two possible quads are discussed, the 9‐node (or complete) quad together with the 8‐node (or Serendipity) quad. A validity condition, the positivity of the Jacobian, is exhibited for these two elements. The discussion continues with a rational Bézier quad patch that can be used as a finite element. Extension to arbitrary degrees is given. Copyright © 2014 John Wiley & Sons, Ltd.     

Determination of high‐order fields for multianisotropic material antiplane V‐notches and inclusions by the asymptotic expansion technique and an overdeterministic method

In this paper, the singular behavior for anisotropic multimaterial V‐notched plates is investigated under antiplane shear loading condition. Firstly, the elastic governing equations are transformed into eigen ordinary differential equations through introducing the asymptotic expansions of displacements near the notch tip. The stress singularity exponents, including the higher‐order terms, and the corresponding eigen angular functions are then obtained by solving the established equations by using the interpolating matrix method. Thus, using the combination of the results from finite element analyses and the derived asymptotic expansion, an overdeterministic method is employed to calculate the amplitudes of the coefficients in the asymptotic expansions. Finally, the stress and displacement fields in the vicinity of the notch tip, consisting of both singular terms and higher‐order terms, are determined. The effects of material properties and geometry characteristic on the singular behaviour of the notch tip are discussed in detail.     

A mixed solid‐shell element for the analysis of laminated composites

This paper presents a versatile low order locking‐free mixed solid‐shell element that can be readily employed for a wide range of linear elastic structural analyses, that is, from thick isotropic structures to multilayer anisotropic composites. This solid‐shell element has eight nodes with only displacement degrees of freedom and few assumed stress parameters that provide very accurate interlaminar stress calculations through the element thickness. These elements can be stacked on top of each other to model multilayer structures, fulfilling the interlaminar stress continuity at the interlayer surfaces and zero traction conditions on the top and bottom surfaces of the laminate. The element formulation is based on the well‐known Fraeijs de Veubeke–Hu–Washizu mixed variational principle with enhanced assumed strains formulation and assumed natural strains formulation to alleviate the different types of locking phenomena in solid‐shell elements. The distinct feature of the present formulation is its ability to accurately calculate the interlaminar stress field in multilayer structures, which is achieved by the introduction of a constraint equation on the interlaminar stresses in the Fraeijs de Veubeke–Hu–Washizu principle‐based enhanced assumed strains formulation. The intelligent computer coding of the present formulation makes the present element appropriate for a wide range of structural analyses. To assess the present formulation's accuracy, a variety of popular numerical benchmark examples related to element convergence, mesh distortion, and shell and laminated composite analyses are investigated and the results are compared with those available in the literature. These benchmark examples reveal that the proposed formulation provides very good results for the structural analysis of shells and multilayer composites. Copyright © 2011 John Wiley & Sons, Ltd.     

Thermo‐elasto‐plastic finite element analysis of quasi‐state processes in Eulerian reference frames

An Eulerian finite element formulation for quasi‐state one way coupled thermo‐elasto‐plastic systems is presented. The formulation is suitable for modeling material processes such as welding and laser surfacing. In an Eulerian frame, the solution field of a quasi‐state process becomes steady state for the heat transfer problem and static for the stress problem. A mixed small deformation displacement elasto‐plastic formulation is proposed. The formulation accounts for temperature dependent material properties and exhibits a robust convergence. Streamline upwind Petrov–Galerkin (SUPG) is used to remove spurious oscillations. Smoothing functions are introduced to relax the non‐differentiable evolution equations and allow for the use of gradient (stiffness) solution scheme via the Newton–Raphson method. A 3‐dimensional simulation of a laser surfacing process is presented to exemplify the formulation. Results from the Eulerian formulation are in good agreement with results from the conventional Lagrangian formulation. However, the Eulerian formulation is approximately 15 times faster than the Lagrangian. Copyright © 2001 John Wiley & Sons, Ltd.     

A fictitious domain method for the simulation of thermoelastic deformations in NC‐milling processes

This paper presents a (higher‐order) finite element approach for the simulation of heat diffusion and thermoelastic deformations in NC‐milling processes. The inherent continuous material removal in the process of the simulation is taken into account via continuous removal‐dependent refinements of a paraxial hexahedron base‐mesh covering a given workpiece. These refinements rely on isotropic bisections of these hexahedrons along with subdivisions of the latter into tetrahedrons and pyramids in correspondence to a milling surface triangulation obtained from the application of the marching cubes algorithm. The resulting mesh is used for an element‐wise defined characteristic function for the milling‐dependent workpiece within that paraxial hexahedron base‐mesh. Using this characteristic function, a (higher‐order) fictitious domain method is used to compute the heat diffusion and thermoelastic deformations, where the corresponding ansatz spaces are defined for some hexahedron‐based refinement of the base‐mesh. Numerical experiments compared to real physical experiments exhibit the applicability of the proposed approach to predict deviations of the milled workpiece from its designed shape because of thermoelastic deformations in the process.     

Interface handling for three‐dimensional higher‐order XFEM‐computations in fluid–structure interaction

Three‐dimensional higher‐order eXtended finite element method (XFEM)‐computations still pose challenging computational geometry problems especially for moving interfaces. This paper provides a method for the localization of a higher‐order interface finite element (FE) mesh in an underlying three‐dimensional higher‐order FE mesh. Additionally, it demonstrates, how a subtetrahedralization of an intersected element can be obtained, which preserves the possibly curved interface and allows therefore exact numerical integration. The proposed interface algorithm collects initially a set of possibly intersecting elements by comparing their ‘eXtended axis‐aligned bounding boxes’. The intersection method is applied to a highly reduced number of intersection candidates. The resulting linearized interface is used as input for an elementwise constrained Delaunay tetrahedralization, which computes an appropriate subdivision for each intersected element. The curved interface is recovered from the linearized interface in the last step. The output comprises triangular integration cells representing the interface and tetrahedral integration cells for each intersected element. Application of the interface algorithm currently concentrates on fluid–structure interaction problems on low‐order and higher‐order FE meshes, which may be composed of any arbitrary element types such as hexahedra, tetrahedra, wedges, etc. Nevertheless, other XFEM‐problems with explicitly given interfaces or discontinuities may be tackled in addition. Multiple structures and interfaces per intersected element can be handled without any additional difficulties. Several parallelization strategies exist depending on the desired domain decomposition approach. Numerical test cases including various geometrical exceptions demonstrate the accuracy, robustness and efficiency of the interface handling. Copyright © 2009 John Wiley & Sons, Ltd.