Finite element analysis of elasto-plastic plate bending problems using transition rectangular plate elements

In this work, the finite element analysis of the elasto-plastic plate bending problems is carried out using transition rectangular plate elements. The shape functions of the transition plate elements are derived based on a practical rule. The transition plate elements are all quadrilateral and can be used to obtain efficient finite element models using minimum number of elements. The mesh convergence rates of the models including the transition elements are compared with the regular element models. To verify the developed elements, simple tests are demonstrated and various elasto-plastic problems are solved. Their results are compared with ANSYS results.

NEAR CRACK LINE ELASTIC-PLASTIC ANALYSIS FOR A FINITE PLATE LOADED BY TWO PAIRS OF ANTIPLANE POINT FORCES

In this paper, the improved near crack line analysis method proposed in Refs. [1] and [2] is used to investigate a center crack loaded by two pairs of antiplane point forces in a finite plate in an elastic-perfectly plastic solid. And the analytical solutions are obtained, that is elastic-plastic stress fields near the crack line and the law that the length of the plastic zone along the crack line is varied with an external loads tr,ld the bearing capacity of a finite plate with a center crack. The results of this paper are sufficiently precise near the crack line, because the assumptions of the small scale yielding theory have not been used and no other assumptions have been taken.

Finite element analysis of stress at implant-bone interface of dental implants with different structures

The effect of structure,elastic modulus and thickness of lower modulus layer in porous titanium implants on the stress distribution at the implant-bone interface was investigated.Three-dimensional finite element models of different titanium implants were constructed.The structures of the implants included the whole lower modulus style （No.1）,bio-mimetic style （No.2）,the whole lower modulus style in cancellous bone （No.3） and the whole dense style No.4.The stress distributions at bone-implant interface under static loading were analyzed using Ansys Workbench 10.0 software.The results indicated that the distribution of interface stress is strongly depended on the structure of the implants.The maximum stresses in cancellous bone and root region of implant No.2 are lower than those in the other three implants.A decrease in the modulus of the low modulus layer facilitates the interface stress transferring.Increasing the thickness of the low modulus layer can reduce the stress and induce a more uniform stress distribution at the interface.Among the four implants,biomimetic style implant No.2 is superior in transferring implant-bone interface stress to surrounding bones.

SUPER-ELASTIC CONSTITUTIVE MODEL CONSIDERING PLASTICITY AND ITS FINITE ELEMENT IMPLEMENTATION

Based on the experimental results of super-elastic NiTi alloy, a three-dimensional super-elastic constitutive model including both of stress-induced martensite transformation and plasticity is constructed in a framework of general inelasticity. In the proposed model, transformation hardening, reverse transformation of stress-induced martensite, elastic mismatch between the austenite and martensite phases, and temperature-dependence of transformation stress and elastic modulus of each phase are considered. The plastic yielding of martensite occurred under high stress is addressed by a bilinear isotropic hardening rule. Drucker-Prager-typed transformation surfaces are employed to describe the asymmetric behavior of NiTi alloy in tension and compression. The prediction capability of the proposed model is verified by comparing the simulated results with the correspondent experimental ones. Based on backward Euler＇s integration, a new expression of consistent tangent modulus is derived. The proposed model is then implemented into a finite element package ABAQUS by user-subroutine UMAT. Finally, the validity of such implementation was verified by some numerical samples.

Finite Element Implementation of the Exponential Drucker-Prager Plasticity Model for Adhesive Joints

This paper deals with the numerical implementation of the exponential Drucker-Parger plasticitymodel in the commercial finite element software,ABAQUS,via user subroutine UMAT for adhesive joint simulations.The influence of hydrostatic pressure on adhesive strength was investigated by a modified Arcan fixture designed particularly to induce a different state of hydrostatic pressure within an adhesive layer.The developed user subroutine UMAT,which utilizes an associated plastic flow during a plastic deformation,can provide a good agreement between the simulations and the experimental data.Better numerical stability at highly positive hydrostatic pressure loads for a very high order of exponential function can also be achieved compared to when a non-associated flow is used.

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.

Finite Element Analysis of Elastic Deformation in Blade Machining

In the process of milling the half-open impeller, the milling force ACTS on the half-open impeller will cause the machining deformation error of the half-open impeller. The impeller material studied in this paper is 7050-t7451 aluminum alloy, whose elastic modulus is 70.3 GPa, with low stiffness. Therefore, during milling, it is easy to cause elastic deformation of the cutter, thus causing machining errors and reducing machining surface quality. In addition, due to the low stiffness of blade surface, the deformation of blade surface is much larger than that of hub surface. Therefore, this paper only studies the deformation law of blade surface in the process of milling.

Bending of small-scale Timoshenko beams based on the integral/differential nonlocal-micropolar elasticity theory: a finite element approach

A novel size-dependent model is developed herein to study the bending behavior of beam-type micro/nano-structures considering combined effects of nonlocality and micro-rotational degrees of freedom. To accomplish this aim, the micropolar theory is combined with the nonlocal elasticity. To consider the nonlocality, both integral (original) and differential formulations of Eringen’s nonlocal theory are considered. The beams are considered to be Timoshenko-type, and the governing equations are derived in the variational form through Hamilton’s principle. The relations are written in an appropriate matrix-vector representation that can be readily utilized in numerical approaches. A finite element (FE) approach is also proposed for the solution procedure. Parametric studies are conducted to show the simultaneous nonlocal and micropolar effects on the bending response of small-scale beams under different boundary conditions.

Finite Element Analysis of the Steel Reinforced Plastic Pipe

The steel reinforced plastic pipe is a kind of green environmental protection pipelines with double-sides corrosionresisting and better withstanding to medium working pressure. The structure and technical process of this pipe are described briefly in this paper, and the finite element analysis has been done for the sake of understanding the distributions of stress and displacement inside this pipe under hydrostatic pressure. The analysis results are very important for safety application of the steel reinforced plastic pipe.

Finite element analysis of slope stability by expanding the mobilized principal stress Mohr's circles-Development, encoding and validation

In recent years,finite element analysis is increasingly being proposed in slope stability problems as a competitive method to traditional limit equilibrium methods(LEMs)which are known for their inherent deficiencies.However,the application of finite element method(FEM)to slope stability as a strength reduction method(SRM)or as finite element limit analysis(FELA)is not always a success for the drawbacks that characterize both methods.To increase the performance of finite element analysis in this problem,a new approach is proposed in this paper.It consists in gradually expanding the mobilized stress Mohr’s circles until the soil failure occurs according to a prescribed non-convergence criterion.The present approach called stress deviator increasing method(SDIM)is considered rigorous for three main reasons.Firstly,it preserves the definition of the factor of safety(FOS)as the ratio of soil shear strength to the mobilized shear stress.Secondly,it maintains the progressive development of shear stress resulting from the increase in the principal stress deviator on the same plane,on which the shear strength takes place.Thirdly,by introducing the concept of equivalent stress loading,the resulting trial stresses are checked against the violation of the actual yield criterion formed with the real strength parameters rather than those reduced by a trial factor.The new numerical procedure was encoded in a Fortran computer code called S^(4)DINA and verified by several examples.Comparisons with other numerical methods such as the SRM,gravity increasing method(GIM)or even FELA by assessing both the FOS and contours of equivalent plastic strains showed promising results.

A MULTI-COUPLED FINITE ELEMENT ANALYSIS OF RESISTANCE SPOT WELDING PROCESS

A two-dimensional axisymmetric finite element model is developed to analyze the transient thermal and mechanical behaviors of the Resistance Spot Welding （RSW） process using commercial software ANSYS. Firstly a direct-coupled electrical-thermal Finite Element Analysis （FEA） is performed to analyze the transient thermal characteristics of the RSW process. Then based on the thermal results a sequential coupled thermo-elastic-plastic analysis is conducted to determine the mechanical features of the RSW process. The thermal history of the whole process and the temperature distribution of the weldment are obtained through the analysis. The mechanical features, including the distributions of the contact pressure at both the faying surface and the electrode-workpiece interface, the stress and strain distributions in the weldment and their changes during the RSW process, the deformation of the weldment and the electrode displacement are also calculated.

Elastic-plastic analytical solution for centric crack loaded by two pairs of point shear forces in finite plate

The near crack line analysis method was used to investigate a centric crack loaded by two pairs of point shear forces in a finite plate, and the analytical solution was obtained. The solution includes the unit normal vector of the elastic-plastic boundary near the crack line, the elastic-plastic stress fields near the crack line, the variations of the length of the plastic zone along the crack line with an external load, and the bearing capacity of a finite plate with a centric crack loaded by two pairs of point shear forces. The results are sufficiently precise near the crack line because the assumptions of small scale yielding theory have not been made and no other assumptions are taken.

Analysis of the elastic-plastic indentation properties of materials with varying ratio of hardness to Young’s modulus

The elastic-plastic indentation properties of materials with varying ratio of hardness to Young’s modulus(H/E) were analyzed with the finite element method. And the indentation stress fields of materials with varying ratio H/E on the surface were studied by the experiment. The results show that the penetration depth, contact radius, plastic pile-up and the degree of elastic recovery depend strongly on the ratio H/E. Moreover, graphs were established to describe the relationship between the elastic-plastic indentation parameters and H/E. The established graphs can be used to predict the H/E of materials when compared with experimental data.

Seismic elastic-plastic time history analysis and reliability study of quayside container crane

Quayside container crane is a kind of huge dimension steel structure,which is the major equipment used for handling container at modern ports.With the aim to validate the safety and reliability of the crane under seismic loads,besides conventional analysis,elastic-plastic time history analysis under rare seismic intensity is carried out.An ideal finite element（FEM） elastic-plastic mechanical model of the quayside container crane is presented by using ANSYS codes.Furthermore,according to elastic-plastic time history analysis theory,deformation,stress and damage pattern of the structure under rare seismic intensity are investigated.Based on the above analysis,the established reliability model according to the reliability theory,together with seismic reliability analysis based on Monte-Carlo simulation is applied to practical analysis.The results show that the overall structure of the quayside container crane is generally unstable under rare seismic intensity,and the structure needs to be reinforced.

An elastic-plastic analysis of short-leg shear wall structures during earthquakes

Short-leg shear wall structures are a new form of building structure that combine the merits of both frame and shear wall structures. Its architectural features, structure bearing and engineering cost are reasonable. To analyze the elastic-plastic response of a short-leg shear wall structure during an earthquake, this study modified the multiple-vertical-rod element model of the shear wall, considered the shear lag effect and proposed a multiple-vertical-rod element coupling beam model with a new local stiffness domain. Based on the principle of minimum potential energy and the variational principle, the stiffness matrixes of a short-leg shear wall and a coupling beam are derived in this study. Furthermore, the bending shear correlation for the analysis of different parameters to describe the structure, such as the beam height to span ratio, short-leg shear wall height to thickness ratio, and steel ratio are introduced. The results show that the height to span ratio directly affects the structural integrity; and the short-leg shear wall height to thickness ratio should be limited to a range of approximately 6.0 to 7.0. The design of short-leg shear walls should be in accordance with the ＂strong wall and weak beam＂ principle.

PLASTIC LIMIT ANALYSIS OF INCOMPATIBLE FINITE ELEMENT METHOD

This paper describes an incompatible finite element model satisfying the consistency condition of energy to solve the numerical precision problem of finite element solution in perfectly plastic analysis. In this paper the reason and criterion of the application of the model to plastic limit analysis are discussed, and an algorithm of computing plastic limit load is given.

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 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.

Finite element analyses of slope stability problems using non-associated plasticity

In recent years, finite element analyses have increasingly been utilized for slope stability problems. In comparison to limit equilibrium methods, numerical analyses do not require any definition of the failure mechanism a priori and enable the determination of the safety level more accurately. The paper compares the performances of strength reduction finite element analysis(SRFEA) with finite element limit analysis(FELA), whereby the focus is related to non-associated plasticity. Displacement-based finite element analyses using a strength reduction technique suffer from numerical instabilities when using non-associated plasticity, especially when dealing with high friction angles but moderate dilatancy angles. The FELA on the other hand provides rigorous upper and lower bounds of the factor of safety(FoS) but is restricted to associated flow rules. Suggestions to overcome this problem, proposed by Davis(1968), lead to conservative FoSs; therefore, an enhanced procedure has been investigated. When using the modified approach, both the SRFEA and the FELA provide very similar results. Further studies highlight the advantages of using an adaptive mesh refinement to determine FoSs. Additionally, it is shown that the initial stress field does not affect the FoS when using a Mohr-Coulomb failure criterion.

RIGID FINITE ELEMENT AND LIMIT ANALYSIS

According to the lower-bound theorem of limit analysis the Rigid Finite Element Meth-od(RFEM)is applied to structural limit analysis and the linear programmings for limit analysis are deducedin this paper.Moreover,the Thermo-Parameter Method(TPM)and Parametric Variational principles(PVP)are used to reduce the computational effort while maintaining the accuracy of solutions.A better solution isalso obtained in this paper.