An edge‐based strain smoothing particle finite element method for large deformation problems in geotechnical engineering

2020 ◽  
Vol 44 (7) ◽  
pp. 923-941 ◽  
Author(s):  
Yin‐Fu Jin ◽  
Wei‐Hai Yuan ◽  
Zhen‐Yu Yin ◽  
Yung‐Ming Cheng
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Large Deformation ◽  
Geotechnical Engineering ◽  
Particle Finite Element Method ◽  
Strain Smoothing ◽  
Edge Based ◽  
Large Deformation Problems ◽  
Element Method

Application of the particle finite element method for large deformation consolidation analysis

2019 ◽  
Vol 36 (9) ◽  
pp. 3138-3163 ◽  
Author(s):  
Wei-Hai Yuan ◽  
Wei Zhang ◽  
Beibing Dai ◽  
Yuan Wang
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Large Deformation ◽  
Excess Pore Pressure ◽  
Triangular Element ◽  
Particle Finite Element Method ◽  
Content Type ◽  
Modified Cam Clay ◽  
Large Deformation Problems ◽  
Element Method

Purpose Large deformation problems are frequently encountered in various fields of geotechnical engineering. The particle finite element method (PFEM) has been proven to be a promising method to solve large deformation problems. This study aims to develop a computational framework for modelling the hydro-mechanical coupled porous media at large deformation based on the PFEM. Design/methodology/approach The PFEM is extended by adopting the linear and quadratic triangular elements for pore water pressure and displacements. A six-node triangular element is used for modelling two-dimensional problems instead of the low-order three-node triangular element. Thus, the numerical instability induced by volumetric locking is avoided. The Modified Cam Clay (MCC) model is used to describe the elasto-plastic soil behaviour. Findings The proposed approach is used for analysing several consolidation problems. The numerical results have demonstrated that large deformation consolidation problems with the proposed approach can be accomplished without numerical difficulties and loss of accuracy. The coupled PFEM provides a stable and robust numerical tool in solving large deformation consolidation problems. It is demonstrated that the proposed approach is intrinsically stable. Originality/value The PFEM is extended to consider large deformation-coupled hydro-mechanical problem. PFEM is enhanced by using a six-node quadratic triangular element for displacement and this is coupled with a four-node quadrilateral element for modelling excess pore pressure.

A Concept of Cell-Based Smoothed Finite Element Method Using 10-Node Tetrahedral Elements (CS-FEM-T10) for Large Deformation Problems of Nearly Incompressible Solids

2019 ◽  
Vol 17 (02) ◽  
pp. 1845009
Author(s):  
Yuki Onishi
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Large Deformation ◽  
Hydrostatic Stress ◽  
Reaction Force ◽  
Strain Smoothing ◽  
Smoothed Finite Element Method ◽  
Incompressible Solids ◽  
Large Deformation Problems ◽  
Element Method

A new concept of smoothed finite element method (S-FEM) using 10-node tetrahedral (T10) elements, CS-FEM-T10, is proposed. CS-FEM-T10 is a kind of cell-based S-FEM (CS-FEM) and thus it smooths the strain only within each T10 element. Unlike the other types of S-FEMs [node-based S-FEM (NS-FEM), edge-based S-FEM (ES-FEM), and face-based S-FEM (FS-FEM)], CS-FEM can be implemented in general finite element (FE) codes as a piece of the element library. Therefore, CS-FEM-T10 is also compatible with general FE codes as a T10 element. A concrete example of CS-FEM-T10 named SelectiveCS-FEM-T10 is introduced for large deformation problems of nearly incompressible solids. SelectiveCS-FEM-T10 subdivides each T10 element into 12 four-node tetrahedral (T4) subelements with an additional dummy node at the element center. Two types of strain smoothing are conducted for the deviatoric and hydrostatic stress evaluations and the selective reduced integration (SRI) technique is utilized for the stress integration. As a result, SelectiveCS-FEM-T10 avoids the shear/volumetric locking, pressure checkerboarding, and reaction force oscillation in nearly incompressible solids. In addition, SelectiveCS-FEM-T10 has a relatively long-lasting property in large deformation problems. A few examples of large deformation analyses of a hyperelastic material confirm the good accuracy and robustness of SelectiveCS-FEM-T10. Moreover, an implementation of SelectiveCS-FEM-T10 in the FE code ABAQUS as a user-defined element (UEL) is conducted and its capability is discussed.

scholarly journals GPU-accelerated smoothed particle finite element method for large deformation analysis in geomechanics

2021 ◽  
Author(s):  
Wei Zhang ◽  
Zhi-hao Zhong ◽  
Chong ◽  
Wei-hai Yuan ◽  
Wei Wu
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Parallel Computing ◽  
Large Deformation ◽  
Computational Efficiency ◽  
Particle Finite Element Method ◽  
Memory Space ◽  
Strain Smoothing ◽  
Smoothed Particle ◽  
Element Method

Particle finite element method (PFEM) is an effective numerical tool for solving large-deformation problems in geomechanics. By incorporating the node integration technique with strain smoothing into the PFEM, we proposed the smoothed particle ?nite element method (SPFEM). This paper extends the SPFEM to three-dimensional cases and presents a SPFEM executed on graphics processing units (GPUs) to boost the computational efficiency. The detailed parallel computing strategy on GPU is introduced. New computation formulations related to the strain smoothing technique are proposed to save memory space in the GPU parallel computing. Several benchmark problems are solved to validate the proposed approach and to evaluate the GPU acceleration performance. Numerical examples show that with the new formulations not only the memory space can be saved but also the computational efficiency is improved. The computational cost is reduced by 70% for the double-precision GPU parallel computing with the new formulations.

Accurate Viscoelastic Large Deformation Analysis Using F-Bar Aided Edge-Based Smoothed Finite Element Method for 4-Node Tetrahedral Meshes (F-BarES-FEM-T4)

2019 ◽  
Vol 17 (02) ◽  
pp. 1845003 ◽  
Author(s):  
Yuki Onishi ◽  
Ryoya Iida ◽  
Kenji Amaya
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Large Deformation ◽  
Pressure Oscillation ◽  
Deformation Analysis ◽  
Smoothed Finite Element Method ◽  
Edge Based ◽  
Smoothing Procedure ◽  
Element Method

A state-of-the-art tetrahedral smoothed finite element method, F-barES-FEM-T4, is demonstrated on viscoelastic large deformation problems. The stress relaxation of viscoelastic materials brings near incompressibility when the long-term Poisson’s ratio is close to 0.5. The conventional hybrid 4-node tetrahedral (T4) elements cannot avoid the shear locking and pressure checkerboarding issues, meanwhile F-barES-FEM-T4 can suppress these issues successfully by adopting the edge-based smoothed finite element method (ES-FEM) with the aid of the F-bar method and the cyclic smoothing procedure. A few examples of analyses verify that F-barES-FEM-T4 is locking-free and pressure oscillation-free in viscoelastic analyses as well as in nearly incompressible hyperelastic or elastoplastic analyses.

F-Bar Aided Edge-Based Smoothed Finite Element Method with 4-Node Tetrahedral Elements for Static Large Deformation Elastoplastic Problems

2019 ◽  
Vol 16 (05) ◽  
pp. 1840010 ◽  
Author(s):  
Yuki Onishi
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Large Deformation ◽  
Element Formulation ◽  
Tetrahedral Elements ◽  
Smoothed Finite Element Method ◽  
New Type ◽  
Displacement Based ◽  
Edge Based ◽  
Element Method

A new type of smoothed finite element method (S-FEM), F-barES-FEM-T4, is demonstrated in static large deformation elastoplastic cases. F-barES-FEM-T4 combines the edge-based S-FEM (ES-FEM) and the node-based S-FEM (NS-FEM) for 4-node tetrahedral (T4) elements with the aid of the F-bar method in order to resolve the major issues of Selective ES/NS-FEM-T4. As well as most of the other S-FEMs, F-barES-FEM-T4 inherits pure displacement-based formulation and thus has no increase in DOF. Moreover, the cyclic smoothing procedure introduced in F-barES-FEM-T4 is effective to adjust the smoothing level so that pressure checkerboarding (oscillation) is suppressed reasonably. Some examples of static large deformation analyses for elastoplastic materials proof the excellent performance of F-barES-FEM-T4 in contrast to the conventional hybrid T4 element formulation.

scholarly journals Numerical Simulation of Hypervelocity Impact FEM-SPH Algorithm Based on Large Deformation of Material

2016 ◽  
Author(s):  
Aimin Yang ◽  
Jinze Li ◽  
Hengheng Qu ◽  
Yuhang Pan ◽  
Yanhong Kang ◽  
...  
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Large Deformation ◽  
Particle Method ◽  
The Finite Element Method ◽  
Smoothed Particle ◽  
Sph Algorithm ◽  
Coupling Algorithm ◽  
Large Deformation Problems ◽  
Element Method

In this paper, we first discuss the research status and application progress of the finite element method and the smoothed particle method. By analyzing the advantages of the smoothed particle method and the finite element method, a new coupling algorithm, namely FEM-SPH algorithm, is proposed. By the method of comparison, it shows that finite element method and SPH method in the simulation of large deformation problems each have advantages and disadvantages, the finite element method smoothed particle coupling algorithm is effective to achieve the performance of high computational efficiency and can naturally simulate large deformation problems across. In the process of calculation, the large deformation unit can be freely into an algorithm to facilitate the calculation accuracy and efficiency of three methods of numerical simulation. Through the study found, FEM-SPH algorithm not only overcome the defect of smooth particle tensile instability, but also overcomes the problem of low efficiency of finite element computation. To further test the FEM-SPH algorithm has advantages in the practical engineering, we have carried out the actual test to the example of the super high speed collision, concluded that, since the target of most of the computational domain is always finite element, smoothed particle focused only in contact with the projectile and target of local area, particle number is not much, the whole calculation process just ten minutes, computational efficiency has been greatly improved, at the same time in the simulation of large deformation, the advantage is very obvious .This provides a criterion for the actual project, depending on the specific material deformation mode and choose a more appropriate conversion algorithm.

Sign in / Sign up

Export Citation Format

Share Document