FINITE ELEMENT ANALYSIS FOR CONSOLIDATION IN INTERACTION BETWEEN STRUCTURE AND SATURATED SOIL FOUNDATION

The consolidation analysis of interaction between structure and saturated soil foundation is discussed. With the use of substructure technique, the structure is condensed onto the interface of the soil, and then the consolidation governing equations to describe the interaction between soil and structure are derived, The solution with non-iterative algorithm is proposed in this paper. The pressure Master-Slave relation method is used to deal with the non-permeability conditions on soil boundaries. A numerical example is illustrated. Based on this paper, the interactive consolidation analysis between large structure and soil has been more practical.

NONLINEAR AND ELASTO-PLASTICITY CONSOLIDATION MODELS OF UNSATURATED SOIL AND APPLICATIONS

The non-linear constitutive model suggested by the authors and the Alonso's elasto-plasticity model of unsaturated soil modified by the authors are introduced into the consolidation theory of unsaturated soil proposed by CHEN Zheng-han, and the non-linear and the elasto-plasticity consolidation models of unsaturated soil are obtained. Programs related to the two consolidation models are designed, and a 2-D consolidation problem of unsaturated sail is solved using the programs, the consolidation process and the development of plastic;one under multi-grade bad are studied. The above research develops the consolidation theory of unsaturated soil to a new level.

A New Mixed Finite Element Method for Biot Consolidation Equations

In this paper,we study a new finite element method for poroelasticity problem with homogeneous boundary conditions.The finite element discretization method is based on a three-variable weak form with mixed finite element for the linear elasticity,i.e.,the stress tensor,displacement and pressure are unknown variables in the weak form.For the linear elasticity formula,we use a conforming finite element proposed in[11]for the mixed form of the linear elasticity and piecewise continuous finite element for the pressure of the fluid flow.We will show that the newly proposed finite element method maintains optimal convergence order.

NONLINEAR AND ELASTO_PLASTICITY CONSOLIDATION MODELS OF UNSATURATED SOIL AND APPLICATIONS

The non_linear constitutive model suggested by the authors and the Alonso’s elasto_plasticity model of unsaturated soil modified by the authors are introduced into the consolidation theory of unsaturated soil proposed by CHEN Zheng_han, and the non_linear and the elasto_plasticity consolidation models of unsaturated soil are obtained. Programs related to the two consolidation models are designed, and a 2_D consolidation problem of unsaturated soil is solved using the programs, the consolidation process and the development of plastic zone under multi_grade load are studied. The above research develops the consolidation theory of unsaturated soil to a new level.

3D finite element method （FEM） simulation of groundwater flow during backward erosion piping

Backward erosion piping is an important failure mechanism for cohesive water retaining structures which are founded on a sandy aquifer. At present, the prediction models for safety assessment are often based on 2D assumptions. In this work, a 3D numerical approach of the groundwater flow leading to the erosion mechanism of backward erosion piping is presented and discussed. Comparison of the 2D and 3D numerical results explicitly demonstrates the inherent 3D nature of the piping phenomenon. In addition, the influence of the seepage length is investigated and discussed for both piping initiation and piping progression. The results clearly indicate the superiority of the presented 3D numerical model compared to the established 2D approach. Moreover, the 3D numerical results enable a better understanding of the complex physical mechanism involved in backward erosion piping and thus can lead to a significant improvement in the safety assessment of water retaining structures.

Semi-analytical finite element method applied for characterizing micropolar fibrous composites

A semi-analytical finite element method(SAFEM),based on the two-scale asymptotic homogenization method(AHM)and the finite element method(FEM),is implemented to obtain the effective properties of two-phase fiber-reinforced composites(FRCs).The fibers are periodically distributed and unidirectionally aligned in a homogeneous matrix.This framework addresses the static linear elastic micropolar problem through partial differential equations,subject to boundary conditions and perfect interface contact conditions.The mathematical formulation of the local problems and the effective coefficients are presented by the AHM.The local problems obtained from the AHM are solved by the FEM,which is denoted as the SAFEM.The numerical results are provided,and the accuracy of the solutions is analyzed,indicating that the formulas and results obtained with the SAFEM may serve as the reference points for validating the outcomes of experimental and numerical computations.

P_1-nonconforming triangular finite element method for elliptic and parabolic interface problems

The lowest order Pl-nonconforming triangular finite element method （FEM） for elliptic and parabolic interface problems is investigated. Under some reasonable regularity assumptions on the exact solutions, the optimal order error estimates are obtained in the broken energy norm. Finally, some numerical results are provided to verify the theoretical analysis.

Reduced-order finite element method based on POD for fractional Tricomi-type equation

The reduced-order finite element method （FEM） based on a proper orthogo- nal decomposition （POD） theory is applied to the time fractional Tricomi-type equation. The present method is an improvement on the general FEM. It can significantly save mem- ory space and effectively relieve the computing load due to its reconstruction of POD basis functions. Furthermore, the reduced-order finite element （FE） scheme is shown to be un- conditionally stable, and error estimation is derived in detail. Two numerical examples are presented to show the feasibility and effectiveness of the method for time fractional differential equations （FDEs）.

Thermo-hydro-mechanical-air coupling finite element method and its application to multi-phase problems

In this paper, a finite element method （FEM）-based multi-phase problem based on a newly proposed thermal elastoplastic constitutive model for saturated/unsaturated geomaterial is discussed. A program of FEM named as SOFT, adopting unified field equations for thermo-hydro-mechanical-air （THMA） behavior of geomaterial and using finite element-finite difference （FE-FD） scheme for so/l-water-air three-phase coupling problem, is used in the numerical simulation. As an application of the newly proposed numerical method, two engineering problems, one for slope failure in unsaturated model ground and another for in situ heating test related to deep geological repository of high-level radioactive waste （HLRW）, are simulated. The model tests on slope failure in unsaturated Shirasu ground, carried out by Kitamura et al. （2007）, is simulated in the framework of soil-water-air three-phase coupling under the condition of constant temperature. While the in situ heating test reported by Munoz （2006） is simulated in the same framework under the conditions of variable temperature hut constant air pressure.

Anisotropic adaptive finite element method for magnetohydrodynamic flow at high Hartmann numbers

This paper presents an anisotropic adaptive finite element method （FEM） to solve the governing equations of steady magnetohydrodynamic （MHD） duct flow. A resid- ual error estimator is presented for the standard FEM, and two-sided bounds on the error independent of the aspect ratio of meshes are provided. Based on the Zienkiewicz-Zhu es- timates, a computable anisotropic error indicator and an implement anisotropic adaptive refinement for the MHD problem are derived at different values of the Hartmann number. The most distinguishing feature of the method is that the layer information from some directions is captured well such that the number of mesh vertices is dramatically reduced for a given level of accuracy. Thus, this approach is more suitable for approximating the layer problem at high Hartmann numbers. Numerical results show efficiency of the algorithm.

Self-adaptive one-dimensional nonlinear finite element method based on element energy projection method

The element energy projection （EEP） method for computation of super- convergent resulting in a one-dimensional finite element method （FEM） is successfully used to self-adaptive FEM analysis of various linear problems, based on which this paper presents a substantial extension of the whole set of technology to nonlinear problems. The main idea behind the technology transfer from linear analysis to nonlinear analysis is to use Newton＇s method to linearize nonlinear problems into a series of linear problems so that the EEP formulation and the corresponding adaptive strategy can be directly used without the need for specific super-convergence formulation for nonlinear FEM. As a re- sult, a unified and general self-adaptive algorithm for nonlinear FEM analysis is formed. The proposed algorithm is found to be able to produce satisfactory finite element results with accuracy satisfying the user-preset error tolerances by maximum norm anywhere on the mesh. Taking the nonlinear ordinary differential equation （ODE） of second-order as the model problem, this paper describes the related fundamental idea, the imple- mentation strategy, and the computational algorithm. Representative numerical exam- ples are given to show the efficiency, stability, versatility, and reliability of the proposed approach.

Stability assessment of landslide-prone road cut rock slopes in Himalayan terrain:A finite element method based approach

Large-scale slope destabilization could be aggravated due to swift urbanization and ever-rising demands of geoengineering projects such as dams,tunnels,bridges and widening roads.National Highway-58 connects Delhi to Badrinath in India,which passes through complex geomorphological and geological terrain and often encounters cut slopes susceptible to slope failures.In the present investigation,a detailed geotechnical appraisal is conducted along the road cut slopes from Rishikesh to Devprayag in the Himalayas.Twenty vulnerable road cut slopes were demarcated for detailed slope stability analysis using Phase2D finite element modeling simulator.Nonlinear generalized Hoek-Brown(GHB)criterion was adopted for stability analyses.Out of 20 slopes,five slopes(S6,S7,S18,S19 and S20)are unstable with factor of safety(FoS)less than or equal to 1,and thus needs immediate attention.The FoS values of four slopes(S2,S9,S13 and S17)lie between 1 and 1.3,i.e.marginally stable,and slopes S1,S3,S4,S5,S8,S10,Sll,S12,S14,S15 and S16 are stable.Mohr-Coulomb(MC)criterion was also adopted to compare the slope stability analysis with GHB criterion.The FoS calculated from GHB criterion is close to that using MC criterion for lower values of FoS whereas for higher values,the difference is marked.For the jointed rock in the Himalayan region,the nonlinear GHB criterion gives better results as compared to MC criterion and matches with the prevailing field conditions.Accordingly,some suggestions are proposed to strengthen the stability of cut slopes.

Mixed finite element and differential quadrature method for free and forced vibration and buckling analysis of rectangular plates

This paper presents a combined application of the finite element method （FEM） and the differential quadrature method （DQM） to vibration and buckling problems of rectangular plates. The proposed scheme combines the geometry flexibility of the FEM and the high accuracy and efficiency of the DQM. The accuracy of the present method is demonstrated by comparing the obtained results with those available in the literature. It is shown that highly accurate results can be obtained by using a small number of finite elements and DQM sample points. The proposed method is suitable for the problems considered due to its simplicity and potential for further development.

Self-adaptive strategy for one-dimensional finite element method based on EEP method with optimal super-convergence order

Based on the newly-developed element energy projection （EEP） method with optimal super-convergence order for computation of super-convergent results, an improved self-adaptive strategy for one-dimensional finite element method （FEM） is proposed. In the strategy, a posteriori errors are estimated by comparing FEM solutions to EEP super-convergent solutions with optimal order of super-convergence, meshes are refined by using the error-averaging method. Quasi-FEM solutions are used to replace the true FEM solutions in the adaptive process. This strategy has been found to be simple, clear, efficient and reliable. For most problems, only one adaptive step is needed to produce the required FEM solutions which pointwise satisfy the user specified error tolerances in the max-norm. Taking the elliptical ordinary differential equation of the second order as the model problem, this paper describes the fundamental idea, implementation strategy and computational algorithm and representative numerical examples are given to show the effectiveness and reliability of the proposed approach.

DEFORMATION ANALYSIS OF SHEET METAL SINGLE-POINT INCREMENTAL FORMING BY FINITE ELEMENT METHOD SIMULATION

Single-point incremental forming （SPIF） is an innovational sheet metal forming method without dedicated dies, which belongs to rapid prototyping technology. In generalizing the SPIF of sheet metal, the deformation analysis on forming process becomes an important and useful method for the planning of shell products, the choice of material, the design of the forming process and the planning of the forming tool. Using solid brick elements, the finite element method（FEM） model of truncated pyramid was established. Based on the theory of anisotropy and assumed strain formulation, the SPIF processes with different parameters were simulated. The resulted comparison between the simulations and the experiments shows that the FEM model is feasible and effective. Then, according to the simulated forming process, the deformation pattern of SPIF can be summarized as the combination of plane-stretching deformation and bending deformation. And the study about the process parameters＇ impact on deformation shows that the process parameter of interlayer spacing is a dominant factor on the deformation. Decreasing interlayer spacing, the strain of one step decreases and the formability of blank will be improved. With bigger interlayer spacing, the plastic deformation zone increases and the forming force will be bigger.

SELF-ADAPTIVE STRATEGY FOR ONE-DIMENSIONAL FINITE ELEMENT METHOD BASED ON ELEMENT ENERGY PROJECTION METHOD

Based on the newly-developed element energy projection （EEP） method for computation of super-convergent results in one-dimensional finite element method （FEM）, the task of self-adaptive FEM analysis was converted into the task of adaptive piecewise polynomial interpolation. As a result, a satisfactory FEM mesh can be obtained, and further FEM analysis on this mesh would immediately produce an FEM solution which usually satisfies the user specified error tolerance. Even though the error tolerance was not completely satisfied, one or two steps of further local refinements would be sufficient. This strategy was found to be very simple, rapid, cheap and efficient. Taking the elliptical ordinary differential equation of second order as the model problem, the fundamental idea, implementation strategy and detailed algorithm are described. Representative numerical examples are given to show the effectiveness and reliability of the proposed approach.

Experimental and Finite Element Method Studies of J-Lead Solder Joint Reliability

A comprehensive experimental and numerical study of solder joints for plastic leaded chip carrier (PLCC) 84-Pin, 1.27 mm pitch was carried out. The reliability of solder joints was assessed through accelerated thermal cycling at the temperature range of - 55℃-125℃. The samples were taken out to observe the evolution in microstructure, such as grain coarsening, initiation and propagation of cracks. It was found that the Pb-rich phases segregated gradually and formed a continuous layer adjacent to the intermetallic compound (IMC) layer with increasing the number of thermal cycles, resulting in cracks near the solder/lead interface. The response of stress and strain was studied using nonlinear finite element method (FEM), and the results agreed well with the experimental data.

Radial Force Analysis of Polydioxanone Coronary Stent by Finite Element Method

Coronary stent is used to treat stenosis artery by recovering the luminal diameter of artery and maintaining the normal blood flow. The geometry of coronary stent is an important factor for the radial force. In this study,the relation between the radial force of stent and crown angle was discussed. The result showed that there was no particular rule between the radial force of stent and the crown angle. The maximum radial force of stent was obtained when the crown angle was 50. 04° and the minimum radial force was got when the crown angle was 75°.

Modeling of contact surface morphology and dust particles by using finite element method

The effect of dust particles on electric contacts and a hazardous size range of hard dust particles using a rigid model were discussed before. As further research, elastic-plastic model of finite element analysis was established in this work, which is closer to real condition. In this work, the behavior of large size and small size particles, and the influence of particles hardness were investigated. The calculating result of small-size particles presents a general hazardous size coefficient for different contact surface morphology; for large-size particles, it presents a hazardous size coefficient for complicated composition of the dust. And the effect of the dust shape is also discussed.

PREDICTION FOR FORMING LIMIT OF AL2024T3 SHEET BASED ON DAMAGE THEORY USING FINITE ELEMENT METHOD

This paper presents the application of anisotropic damage theory to the study of forming limit diagram of A12024T3 aluminum alloy sheet. In the prediction of limiting strains of the aluminum sheet structure, a finite element cell model has been constructed. The cell model consists of two phases, the aluminum alloy matrix and the intermetallic cluster. The material behavior of the aluminum alloy matrix is described with a fully coupled elasto-plastic damage constitutive equation. The intermetallic cluster is assumed to be elastic and brittle. By varying the stretching ratio, the limiting strains of the sheet under biaxial stretching have been predicted by using the necking criterion proposed. The prediction is in good agreement with the experimental findings. Moreover, the finite element cell model can provide information for understanding the microscopic damage mechanism of the aluminum alloy. Over-estimation of the limit strains may result if the effect of material damage is ignored in the sheet metal forming study.