Nonconforming finite elements for the equation of planar elasticity

Two new locking-free nonconforming finite elements for the pure displacement planar elasticity problem are presented. Convergence rates of the elements are uniformly optimal with respect to A. The energy norm and L2 norm errors are proved to be O（h2） and O（h3）, respectively. Numerical tests confirm the theoretical analysis.

NEW SECOND ORDER NONCONFORMING TRIANGULAR ELEMENT FOR PLANAR ELASTICITY PROBLEMS

In the use of finite element methods to the planar elasticity problems,one diffculty is to overcome locking when elasticity constant λ→∞.In the case of traction boundary condition,another diffculty is to make the discrete Korn＇s second inequality valid.In this paper,a triangular element is presented.We prove that this element is locking-free,the discrete Korn＇s second inequality holds and the convergence order is two.

Multiscale Finite Element Method for Coupling Analysis of Heterogeneous Magneto-Electro-Elastic Structures in Thermal Environment

Magneto-electro-elastic (MEE) materials, a new type of composite intelligent materials, exhibit excellent multifield coupling effects. Due to the heterogeneity of the materials, it is challenging to use the traditional finite element method (FEM) for mechanical analysis. Additionally, the MEE materials are often in a complex service environment, especially under the influence of the thermal field with thermoelectric and thermomagnetic effects, which affect its mechanical properties. Therefore, this paper proposes the efficient multiscale computational method for the multifield coupling problem of heterogeneous MEE structures under the thermal environment. The method constructs a multi-physics field with numerical base functions (the displacement, electric potential, and magnetic potential multiscale base functions). It equates a single cell of heterogeneous MEE materials to a macroscopic unit and supplements the macroscopic model with a microscopic model. This allows the problem to be solved directly on a macroscopic scale. Finally, the numerical simulation results demonstrate that compared with the traditional FEM, the multiscale finite element method (MsFEM) can achieve the purpose of ensuring accuracy and reducing the degree of freedom, and significantly improving the calculation efficiency.

A Locking-free Nonconforming Finite Element for Planar Linear Elasticity

We propose a locking-free nonconforming finite element method to solve for the displacement variation in the pure displacement boundary value problem of planar linear elasticity. The method proposed in this paper is robust and optimal, in the sense that the convergence estimate in the energy is independent of the Lame parameter λ.

A LOCKING-FREE ANISOTROPIC NONCONFORMING FINITE ELEMENT FOR PLANAR LINEAR ELASTICITY PROBLEM

The main aim of this article is to study the approximation of a locking-free anisotropic nonconforming finite element for the pure displacement boundary value problem of planar linear elasticity. The optimal error estimates are obtained by using some novel approaches and techniques. The method proposed in this article is robust in the sense that the convergence estimates in the energy and L^2-norms are independent-of the Lame parameter λ.

Immersed Interface Finite Element Methods for Elasticity Interface Problems with Non-Homogeneous Jump Conditions

In this paper,a class of new immersed interface finite element methods (IIFEM) is developed to solve elasticity interface problems with homogeneous and non-homogeneous jump conditions in two dimensions.Simple non-body-fitted meshes are used.For homogeneous jump conditions,both non-conforming and conforming basis functions are constructed in such a way that they satisfy the natural jump conditions. For non-homogeneous jump conditions,a pair of functions that satisfy the same non-homogeneous jump conditions are constructed using a level-set representation of the interface.With such a pair of functions,the discontinuities across the interface in the solution and flux are removed;and an equivalent elasticity interface problem with homogeneous jump conditions is formulated.Numerical examples are presented to demonstrate that such methods have second order convergence.

Implementing the Node Based Smoothed Finite Element Method as User Element in Abaqus for Linear and Nonlinear Elasticity

In this paper,the node based smoothed-strain Abaqus user element(UEL)in the framework of finite element method is introduced.The basic idea behind of the node based smoothed finite element(NSFEM)is that finite element cells are divided into subcells and subcells construct the smoothing domain associated with each node of a finite element cell[Liu,Dai and Nguyen-Thoi(2007)].Therefore,the numerical integration is globally performed over smoothing domains.It is demonstrated that the proposed UEL retains all the advantages of the NSFEM,i.e.,upper bound solution,overly soft stiffness and free from locking in compressible and nearly-incompressible media.In this work,the constant strain triangular(CST)elements are used to construct node based smoothing domains,since any complex two dimensional domains can be discretized using CST elements.This additional challenge is successfully addressed in this paper.The efficacy and robustness of the proposed work is obtained by several benchmark problems in both linear and nonlinear elasticity.The developed UEL and the associated files can be downloaded from https://github.com/nsundar/NSFEM.

Quasi-optimal complexity of adaptive finite element method for linear elasticity problems in two dimensions

This paper introduces an adaptive finite element method （AFEM） using the newest vertex bisection and marking exclusively according to the error estimator without special treatment of oscillation. By the combination of the global lower bound and the localized upper bound of the posteriori error estimator, perturbation of oscillation, and cardinality of the marked element set, it is proved that the AFEM is quasi-optimal for linear elasticity problems in two dimensions, and this conclusion is verified by the numerical examples.

An explicit finite element－finite difference method for analyzing the effect of visco－elastic local topography on the earthquake motion

An explicit finite element-finite difference method for analyzing the effects of two-dimensional visco-elastic localtopography on earthquake ground motion is prOPosed in this paper. In the method, at first, the finite elementdiscrete model is formed by using the artificial boundary and finite element method, and the dynamic equationsof local nodes in the discrete model are obtained according to the theory of the special finite element method similar to the finite difference method, and then the explicit step-by-step integration formulas are presented by usingthe explicit difference method for solving the visco-elastic dynamic equation and Generalized Multi-transmittingBoundary. The method has the advantages of saving computing time and computer memory space, and it is suitable for any case of topography and has high computing accuracy and good computing stability.

ANALYSES ON NONLINEAR COUPLING OF MAGNETO-THERMO-ELASTICITY OF FERROMAGNETIC THIN SHELL—II:FINITE ELEMENT MODELING AND APPLICATION

Based on the generalized vaxiational principle of magneto-thermo-elasticity of a ferromagnetic thin shell established （see, Analyses on nonlinear coupling of magneto-thermo- elasticity of ferromagnetic thin shell--I）, the present paper developed a finite element modeling for the mechanical-magneto-thermal multi-field coupling of a ferromagnetic thin shell. The numerical modeling composes of finite element equations for three sub-systems of magnetic, thermal and deformation fields, as well as iterative methods for nonlinearities of the geometrical large-deflection and the multi-field coupling of the ferromagnetic shell. As examples, the numerical simulations on magneto-elastic behaviors of a ferromagnetic cylindrical shell in an applied magnetic field, and magneto-thermo-elastic behaviors of the shell in applied magnetic and thermal fields are carried out. The results are in good agreement with the experimental ones.

Three-Dimensional Thermo-Elastic-Plastic Finite Element Method Modeling for Predicting Weld-Induced Residual Stresses and Distortions in Steel Stiffened-Plate Structures

The objective of the present paper is to develop nonlinear finite element method models for predicting the weld-induced initial deflection and residual stress of plating in steel stiffened-plate structures. For this purpose, three-dimensional thermo-elastic-plastic finite element method computations are performed with varying plate thickness and weld bead length (leg length) in welded plate panels, the latter being associated with weld heat input. The finite element models are verified by a comparison with experimental database which was obtained by the authors in separate studies with full scale measurements. It is concluded that the nonlinear finite element method models developed in the present paper are very accurate in terms of predicting the weld-induced initial imperfections of steel stiffened plate structures. Details of the numerical computations together with test database are documented.

RATIONAL FINITE ELEMENT METHOD FOR ELASTIC BENDING OF REISSNER PLATES

In this paper; some deformation patterns defined by differential equations of the elastic system are introduced into the revised functional for the incompatible elements. And therefore the rational FEM, which is perfect combination of the analytic methods and numeric methods, has been presented. This new approach satisfies not only the mechanical requirement of the elements but also the geometric requirement of the complex structures. What's more the quadrilateral element obtained accordingly for the elastic bending of thick plates demonstrates such advantages as high precision for computation and availability of accurate integration for stiffness matrices.

RELIABILITY OF ELASTO-PLASTIC STRUCTURE USING FINITE ELEMENT METHOD

A solution of probabilistic FEM for elastic-plastic materials is presented based on the incremental theory of plasticity and a modified initial stress method. The formulations are deduced through a direct differentiation scheme. Partial differentiation of displacement, stress and the performance function can be iteratively performed with the computation of the mean values of displacement and stress. The presented method enjoys the efficiency of both the perturbation method and the finite difference method, but avoids the approximation during the partial differentiation calculation. In order to improve the efficiency, the adjoint vector method is introduced to calculate the differentiation of stress and displacement with respect to random variables. In addition, a time-saving computational method for reliability index of elastic-plastic materials is suggested based upon the advanced First Order Second Moment (FOSM) and by the usage of Taylor expansion for displacement. The suggested method is also applicable to 3-D cases.

AN ANISOTROPIC LOCKING-FREE NONCONFORMING TRIANGULAR FINITE ELEMENT METHOD FOR PLANAR LINEAR ELASTICITY PROBLEM

The main aim of this paper is to study the nonconforming linear triangular Crouzeix- Raviart type finite element approximation of planar linear elasticity problem with the pure displacement boundary value on anisotropic general triangular meshes satisfying the maximal angle condition and coordinate system condition. The optimal order error estimates of energy norm and L2-norm are obtained, which are independent of lame parameter λ. Numerical results are given to demonstrate the validity of our theoretical analysis.Mathematics subject classification： 65N30, 65N15.

Transient analysis of 1D inhomogeneous media by dynamic inhomogeneous finite element method

The dynamic inhomogeneous finite element method is studied for use in the transient analysis of one dimensional inhomogeneous media. The general formula of the inhomogeneous consistent mass matrix is established based on the shape function. In order to research the advantages of this method, it is compared with the general finite element method. A linear bar element is chosen for the discretization tests of material parameters with two fictitious distributions. And, a numerical example is solved to observe the differences in the results between these two methods. Some characteristics of the dynamic inhomogeneous finite element method that demonstrate its advantages are obtained through comparison with the general finite element method. It is found that the method can be used to solve elastic wave motion problems with a large element scale and a large number of iteration steps.

Application of scaled boundary finite element method in static and dynamic fracture problems

The scaled boundary finite element method （SBFEM） is a recently developed numerical method combining advantages of both finite element methods （FEM） and boundary element methods （BEM） and with its own special features as well. One of the most prominent advantages is its capability of calculating stress intensity factors （SIFs） directly from the stress solutions whose singularities at crack tips are analytically represented. This advantage is taken in this study to model static and dynamic fracture problems. For static problems, a remeshing algorithm as simple as used in the BEM is developed while retaining the generality and flexibility of the FEM. Fully-automatic modelling of the mixed-mode crack propagation is then realised by combining the remeshing algorithm with a propagation criterion. For dynamic fracture problems, a newly developed series-increasing solution to the SBFEM governing equations in the frequency domain is applied to calculate dynamic SIFs. Three plane problems are modelled. The numerical results show that the SBFEM can accurately predict static and dynamic SIFs, cracking paths and load-displacement curves, using only a fraction of degrees of freedom generally needed by the traditional finite element methods.

Numerical simulation of coal-bed methane transfer with finite element method

The mathematical model and the numerical simulation for the transfer of coal bed methane were established based on the combination of the porous flow theory and elastic plastic mechanics theory and the numerical solution was given, together with the consideration of the fluid solid interaction between the coal bed gas and coal framework. Then the dispersion for the equation of gas porous flow and coal seam distortion was carried out and the functional analysis equation was obtained. Finally, the coupling solution was educed and calculated by finite element method(FEM) on a model example.

Surrounding rock deformation analysis of underground caverns with multi-body finite element method

Discontinuous deformation problems are common in rock engineering. Numerical analysis methods based on system models of the discrete body can better solve these problems. One of the most effective solutions is discontinuous deformation analysis （DDA） method, but the DDA method brings about rock embedding problems when it uses the strain assumption in elastic deformation and adopts virtual springs to simulate the contact problems. The multi-body finite element method （FEM） proposed in this paper can solve the problems of contact and deformation of blocks very well because it integrates the FEM and multi-body system dynamics theory. It is therefore a complete method for solving discontinuous deformation problems through balance equations of the contact surface and for simulating the displacement of whole blocks. In this study, this method was successfully used for deformation analysis of underground caverns in stratified rock. The simulation results indicate that the multi-body FEM can show contact forces and the stress states on contact surfaces better than DDA, and that the results calculated with the multi-body FEM are more consistent with engineering practice than those calculated with DDA method.

Finite Element Based Analysis of Magnetic Forces between Planar Spiral Coils

This paper elaborates on the magnetic forces between current carrying planar spiral coils. Direct and concentric rings methods are employed in order to calculate the magnetic force between these coils. The results obtained by two calculation methods show the efficiency of the replaced rings method in both simplicity and calculation time. Simula-tions using the Finite Element Method (FEM) are carried out to analyze the distribution of the magnetic flux density around the coils. Also, coils with precise size have been constructed and tested. The experimental results as well as the results obtained by FEM are used to validate the accuracy of the calculations.

An explicit finite element method for dynamic analysis in three-medium coupling system and its application

In this paper, an explicit finite element method to analyze the dynamic responses of three-medium coupled systems with any terrain is developed on the basis of the numerical simulation of the continuous conditions on the bounda-ries among fluid saturated porous medium, elastic single-phase medium and ideal fluid medium. This method is a very effective one with the characteristic of high calculating speed and small memory needed because the formulae for this explicit finite element method have the characteristic of decoupling, and which does not need to solve sys-tem of linear equations. The method is applied to analyze the dynamic response of a reservoir with considering the dynamic interactions among water, dam, sediment and basement rock. The vertical displacement at the top point of the dam is calculated and some conclusions are given.