A 3-D,symmetric, finite element formulation of the Biot equations with application to acoustic wave propagation through an elastic porous medium

A weak solution of the coupled, acoustic-elastic, wave propagation problem for a flexible porous material is proposed for a 3-D continuum. Symmetry in the matrix equations; with respect to both volume, i.e. ‘porous frame’–‘pore fluid’, and surface, i.e. ‘porous frame/pore fluid’–‘non-porous media’, fluid–structure interaction; is ensured with only five unknowns per node; fluid pore pressure, fluid-displacement potential and three Cartesian components of the porous frame displacement field. Taking Biot's general theory as starting point, the discretized form of the equations is derived from a weighted residual statement, using a standard Galerkin approximation and iso-parametric interpolation of the dependent variables. The coupling integrals appearing along the boundary of the porous medium are derived for a number of different surface conditions. The primary application of the proposed symmetric 3-D finite element formulation is modelling of noise transmission in typical transportation vehicles, such as aircraft, cars, etc., where porous materials are used for both temperature and noise insulation purposes. As an example of an application of the implemented finite elements, the noise transmission through a double panel with porous filling and different boundary conditions at the two panel boundaries are analysed. © 1998 John Wiley & Sons, Ltd.     

A stochastic finite element method for real eigenvalue problems

A stochastic finite element method (SFEM) based on local averages of a random vector field is developed for both distinct and repeated eigenvalues. Formulae for the variances and covariances of the eigenvalues and eigenvectors are derived. It is shown in a numerical example that, as the number of elements increases, solutions obtained from the present SFEM formulation converge much faster than those obtained from the SFEM formulation based on mid-point discretization.     

A general high‐order finite element formulation for shells at large strains and finite rotations

For hyperelastic shells with finite rotations and large strains a p‐finite element formulation is presented accommodating general kinematic assumptions, interpolation polynomials and particularly general three‐dimensional hyperelastic constitutive laws. This goal is achieved by hierarchical, high‐order shell models. The tangent stiffness matrices for the hierarchical shell models are derived by computer algebra. Both non‐hierarchical, nodal as well as hierarchical element shape functions are admissible. Numerical experiments show the high‐order formulation to be less prone to locking effects. Copyright © 2003 John Wiley & Sons, Ltd.     

A scaled boundary finite element formulation for dynamic elastoplastic analysis

This study presents the development of the scaled boundary finite element method (SBFEM) to simulate elastoplastic stress wave propagation problems subjected to transient dynamic loadings. Material nonlinearity is considered by first reformulating the SBFEM to obtain an explicit form of shape functions for polygons with an arbitrary number of sides. The material constitutive matrix and the residual stress fields are then determined as analytical polynomial functions in the scaled boundary coordinates through a local least squares fit to evaluate the elastoplastic stiffness matrix and the residual load vector semianalytically. The treatment of the inertial force within the solution of the nonlinear system of equations is also presented within the SBFEM framework. The nonlinear equation system is solved using the unconditionally stable Newmark time integration algorithm. The proposed formulation is validated using several benchmark numerical examples.     

A new finite‐element formulation for electromechanical boundary value problems

In this paper, a new finite‐element formulation for the solution of electromechanical boundary value problems is presented. As opposed to the standard formulation that uses scalar electric potential as nodal variables, this new formulation implements a vector potential from which components of electric displacement are derived. For linear piezoelectric materials with positive definite material moduli, the resulting finite‐element stiffness matrix from the vector potential formulation is also positive definite. If the material is non‐linear in a fashion characteristic of ferroelectric materials, it is demonstrated that a straightforward iterative solution procedure is unstable for the standard scalar potential formulation, but stable for the new vector potential formulation. Finally, the method is used to compute fields around a crack tip in an idealized non‐linear ferroelectric material, and results are compared to an analytical solution. Copyright © 2002 John Wiley & Sons, Ltd.     

A scaled boundary finite element formulation for poroelasticity

This paper develops the scaled boundary finite element formulation for applications in coupled field problems, in particular, to poroelasticity. The salient feature of this formulation is that it can be applied over arbitrary polygons and/or quadtree decomposition, which is widely employed to traverse between small and large scales. Moreover, the formulation can treat singularities of any order. Within this framework, 2 sets of semianalytical, scaled boundary shape functions are used to interpolate the displacement and the pore fluid pressure. These shape functions are obtained from the solution of vector and scalar Laplacian, respectively, which are then used to discretise the unknown field variables similar to that of the finite element method. The resulting system of equations are similar in form as that obtained using standard procedures such as the finite element method and, hence, solved using the standard procedures. The formulation is validated using several numerical benchmarks to demonstrate its accuracy and convergence properties.     

Eulerian formulation using stabilized finite element method for large deformation solid dynamics

This paper describes an Eulerian formulation for large deformation solid dynamics. In the present Eulerian approach, an advective equation is solved using the Stream‐Upwind/Petrov–Galerkin finite element method. The Eulerian finite element method is applied to path‐dependent solid analyses such as impact bar and ductile necking problems. These computational results using the Eulerian finite element method are compared with the results obtained from using the Lagrangian finite element method and an Eulerian formulation based on a finite difference method. Copyright © 2007 John Wiley & Sons, Ltd.     

A finite element method for stability problems in finite elasticity

In this paper a finite element method is developed to treat stability problems in finite elasticity. For this purpose the constitutive equations are formulated in principal stretches which allows a general representation of the derivatives of the strain energy function with respect to the principal stretches. These results can then be used to derive an efficient numerical scheme for the computation of singular points.     

Dispersion analysis of the gradient weighted finite element method for acoustic problems in one,two, and three dimensions

This paper reports a detailed analysis on the numerical dispersion error in solving one-, two-, and three-dimensional acoustic problems governed by the Helmholtz equation using the gradient weighted finite element method (GW-FEM) in comparison with the standard FEM and the modified methods presented in the literatures. The discretized system equations derived based on the gradient weighted operation corresponding to the considered method are first briefed. The discrete dispersion relationships relating the exact and numerical wave numbers defined in different dimensions are then formulated, which will be further used to investigate the dispersion effect mainly caused by the approximation of field variables. The influence of nondimensional wave number and wave propagation angle on the dispersion error is detailedly studied. Comparisons are made with the classical FEM and high-performance algorithms. Results of both theoretical and numerical experiments show that the present method can effectively reduce the pollution effect in computational acoustics owning to its crucial effectiveness in handing the dispersion error in the discrete numerical model.     

A finite element formulation for piezoelectric shell structures considering geometrical and material non‐linearities

In this paper, we present a non‐linear finite element formulation for piezoelectric shell structures. Based on a mixed multi‐field variational formulation, an electro‐mechanical coupled shell element is developed considering geometrically and materially non‐linear behavior of ferroelectric ceramics. The mixed formulation includes the independent fields of displacements, electric potential, strains, electric field, stresses, and dielectric displacements. Besides the mechanical degrees of freedom, the shell counts only one electrical degree of freedom. This is the difference in the electric potential in the thickness direction of the shell. Incorporating non‐linear kinematic assumptions, structures with large deformations and stability problems can be analyzed. According to a Reissner–Mindlin theory, the shell element accounts for constant transversal shear strains. The formulation incorporates a three‐dimensional transversal isotropic material law, thus the kinematic in the thickness direction of the shell is considered. The normal zero stress condition and the normal zero dielectric displacement condition of shells are enforced by the independent resultant stress and the resultant dielectric displacement fields. Accounting for material non‐linearities, the ferroelectric hysteresis phenomena are considered using the Preisach model. As a special aspect, the formulation includes temperature‐dependent effects and thus the change of the piezoelectric material parameters due to the temperature. This enables the element to describe temperature‐dependent hysteresis curves. Copyright © 2011 John Wiley & Sons, Ltd.     

A control volume mixed finite element formulation for tetrahedral elements

The control‐volume mixed finite element method is formulated for and applied to a computational domain consisting of a tetrahedral partitioning to solve the steady groundwater flow equations. Test functions consistent with piecewise constant and piecewise linear pressure distributions are used in the formulation. Comparisons are made with a standard mixed finite element formulation using lowest‐order Raviart Thomas basis functions. Results suggest that the control‐volume based formulation is a viable alternative to the standard formulation. Copyright © 2015 John Wiley & Sons, Ltd.     

A novel finite element formulation for frictionless contact problems

This article advocates a new methodology for the finite element solution of contact problems involving bodies that may undergo finite motions and deformations. The analysis is based on a decomposition of the two-body contact problem into two simultaneous sub-problems, and results naturally in geometrically unbiased discretization of the contacting surfaces. A proposed two-dimensional contact element is specifically designed to unconditionally allow for exact transmission of constant normal traction through interacting surfaces.     

Energy conserving and dissipative time finite element schemes for N‐body stiff problems

Petrov–Galerkin finite element method is adopted to develop a family of temporal integrators, which preserves the feature of energy conservation or numerical dissipation for non‐linear N‐body dynamical systems. This leads to an enhancement of numerical stability and the integrators may therefore offer some advantage for the numerical solution of stiff systems in long‐term simulations. Dynamically tuneable numerical integration is exploited to improve the accuracy of the time‐stepping schemes. Representative simulations for simple non‐linear systems show the performance of the schemes in controlling over or damping out unresolved high frequencies. Copyright © 2004 John Wiley & Sons, Ltd.     

ODDLS: A new unstructured mesh finite element method for the analysis of free surface flow problems

This paper introduces a new stabilized finite element method based on the finite calculus (Comput. Methods Appl. Mech. Eng. 1998; 151 :233–267) and arbitrary Lagrangian–Eulerian techniques (Comput. Methods Appl. Mech. Eng. 1998; 155 :235–249) for the solution to free surface problems. The main innovation of this method is the application of an overlapping domain decomposition concept in the statement of the problem. The aim is to increase the accuracy in the capture of the free surface as well as in the resolution of the governing equations in the interface between the two fluids. Free surface capturing is based on the solution to a level set equation. The Navier–Stokes equations are solved using an iterative monolithic predictor–corrector algorithm (Encyclopedia of Computational Mechanics. Wiley: New York, 2004), where the correction step is based on imposing the divergence‐free condition in the velocity field by means of the solution to a scalar equation for the pressure. Examples of application of the ODDLS formulation (for overlapping domain decomposition level set) to the analysis of different free surface flow problems are presented. Copyright © 2008 John Wiley & Sons, Ltd.     

A new fast finite element method for dislocations based on interior discontinuities

A new technique for the modelling of multiple dislocations based on introducing interior discontinuities is presented. In contrast to existing methods, the superposition of infinite domain solutions is avoided; interior discontinuities are specified on the dislocation slip surfaces and the resulting boundary value problem is solved by a finite element method. The accuracy of the proposed method is verified and its efficiency for multi‐dislocation problems is illustrated. Bounded core energies are incorporated into the method through regularization of the discontinuities at their edges. Though the method is applied to edge dislocations here, its extension to other types of dislocations is straightforward. Copyright © 2006 John Wiley & Sons, Ltd.     

A four-node solid shell element formulation with assumed strain

A set of four-node shell element models based on the assumed strain formulation is considered here. The formulation allows for changes in the shell thickness. As a result, the kinematics of deformation are described by purely vectorial variables, without using rotational angles. The present study investigates the use of bubble function displacements and the assumed strain field. Careful selection of the assumed strain terms generates an element whose order of numerical integration does not increase even when the bubble function displacements are added. Results for the four-node element without any bubble function terms show sensitivity to element distortion. Use of the bubble functions with a carefully chosen assumed strain field greatly improves element performance. © 1998 John Wiley & Sons, Ltd.     

A finite element solution of acoustic propagation in rigid porous media

This paper deals with the acoustical behaviour of a rigid porous material. A finite element method to compute both the response to an harmonic excitation and the free vibrations of a three‐dimensional finite multilayer system consisting of a free fluid and a rigid porous material is considered. The finite element used is the lowest order face element introduced by Raviart and Thomas, that eliminates the spurious or circulation modes with no physical meaning. For the porous medium a Darcy's like model and the Allard–Champoux model are taken into account. The numerical results show that the finite element method allows us to compute the response curve for the coupled system and the complex eigenfrequencies. Some of them have a small imaginary part but there are also overdamped modes. Copyright © 2004 John Wiley & Sons, Ltd.     

Stabilized finite element method for viscoplastic flow: formulation with state variable evolution

A stabilized, mixed finite element formulation for modelling viscoplastic flow, which can be used to model approximately steady‐state metal‐forming processes, is presented. The mixed formulation is expressed in terms of the velocity, pressure and state variable fields, where the state variable is used to describe the evolution of the material's resistance to plastic flow. The resulting system of equations has two sources of well‐known instabilities, one due to the incompressibility constraint and one due to the convection‐type state variable equation. Both of these instabilities are handled by adding mesh‐dependent stabilization terms, which are functions of the Euler–Lagrange equations, to the usual Galerkin method. Linearization of the weak form is derived to enable a Newton–Raphson implementation into an object‐oriented finite element framework. A progressive solution strategy is used for improving convergence for highly non‐linear material behaviour, typical for metals. Numerical experiments using the stabilization method with hierarchic shape functions for the velocity, pressure and state variable fields in viscoplastic flow and metal‐forming problems show that the stabilized finite element method is effective and efficient for non‐linear steady forming problems. Finally, the results are discussed and conclusions are inferred. Copyright © 2002 John Wiley & Sons, Ltd.