A FINITE ELEMENT METHOD FOR STRESS ANALYSIS OR ELASTOPLASTIC BODY WITH POLYGONAL LINE STRAIN——HARDENING

In this paper, the stress-strain curve of material is fitted by polygonal line composed of three lines. According to the theory of proportional loading in elastoplasticity, we simplify the complete stress-strain relations, which are given by the increment theory of elastoplasticity. Thus, the finite element equation with the solution of displacement is derived. The assemblage elastoplastic stiffness matrix can be obtained by adding something to the elastic matrix, hence it will shorten the computing time. The determination of every loading increment follows the von Mises yield criteria. The iterative method is used in computation. It omits the redecomposition of the assemblage stiffness matrix and it will step further to shorten the computing time. Illustrations are given to the high-order element application departure from proportional loading, the computation of unloading fitting to the curve and the problem of load estimation.

THE APPLICATION OF COMPATIBLE STRESS ITERATIVE METHOD IN DYNAMIC FINITE ELEMENT ANALYSIS OF HIGH VELOCITY IMPACT

There is a common difficulty in elastic-plastic impact codes such as EPIC[2,3] NONSAP[4], etc.. Most of these codes use the simple linear functions usually taken from static problem to represent the displacement components. In such finite element formulation, the stress components are constant in each element and they are discontinuous in any two neighboring elements. Therefore, the bases of using the virtual work principle in such elements are unreliable. In this paper, we introduce a new method, namely, the compatible stress iterative method, to eliminate the above-said difficulty. The calculated examples show that the calculation using the new method in dynamic finite element analysis of high velocity impact is valid and stable, and the element stiffness can be somewhat reduced.

A novel twice-interpolation finite element method for solid mechanics problems

Formulation and numerical evaluation of a novel twice-interpolation finite element method （TFEM） is presented for solid mechanics problems. In this method, the trial function for Galerkin weak form is constructed through two stages of consecutive interpolation. The primary interpolation follows exactly the same procedure of standard FEM and is further reproduced according to both nodal values and averaged nodal gradients obtained from primary interpolation. The trial functions thus constructed have continuous nodal gradients and contain higher order polynomial without increasing total freedoms. Several benchmark examples and a real dam problem are used to examine the TFEM in terms of accuracy and convergence. Compared with standard FEM, TFEM can achieve significantly better accuracy and higher convergence rate, and the continuous nodal stress can be obtained without any smoothing operation. It is also found that TFEM is insensitive to the quality of the elemental mesh. In addition, the present TFEM can treat the incompressible material without any modification.

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.

Dynamic Crack Propagation Analysis Using Scaled Boundary Finite Element Method

The prediction of dynamic crack propagation in brittle materials is still an important issue in many engineering fields. The remeshing technique based on scaled boundary finite element method(SBFEM) is extended to predict the dynamic crack propagation in brittle materials. The structure is firstly divided into a number of superelements, only the boundaries of which need to be discretized with line elements. In the SBFEM formulation, the stiffness and mass matrices of the super-elements can be coupled seamlessly with standard finite elements, thus the advantages of versatility and flexibility of the FEM are well maintained. The transient response of the structure can be calculated directly in the time domain using a standard time-integration scheme. Then the dynamic stress intensity factor(DSIF) during crack propagation can be solved analytically due to the semi-analytical nature of SBFEM. Only the fine mesh discretization for the crack-tip super-element is needed to ensure the required accuracy for the determination of stress intensity factor(SIF). According to the predicted crack-tip position, a simple remeshing algorithm with the minimum mesh changes is suggested to simulate the dynamic crack propagation. Numerical examples indicate that the proposed method can be effectively used to deal with the dynamic crack propagation in a finite sized rectangular plate including a central crack. Comparison is made with the results available in the literature, which shows good agreement between each other.

SEMI-ELLIPTIC SURFACE CRACK IN AN ELASTIC SOLID WITH FINITE SIZE UNDER IMPACT LOADING

In this paper a semi-elliptic surface crack problem in an elastic solid of finite size under impact loading is investigated. An analysis is performed by means of fracture dynamics and the finite element method, and a three-dimensional finite element program is developed to compute the dynamic stress intensity factor. The results reveal that the effects of the solid＇s boundary surface, crack surface, material inertia and stress wave interactions play significant roles in dynamic fracture.

Numerical Analysis of Rolling-sliding Contact with the Frictional Heat in Rail

Thermal damage caused by frictional heat of rolling-sliding contact is one of the most important failure forms of wheel and rail. Many studies of wheel-rail frictional heating have been devoted to the temperature field, but few literatures focus on wheel-rail thermal stress caused by frictional heating. However, the wheel-rail creepage is one of important influencing factors of the thermal stress In this paper, a thermo-mechanical coupling model of wheel-rail rolling-sliding contact is developed using thermo-elasto-plastic finite element method. The effect of the wheel-rail elastic creepage on the distribution of heat flux is investigated using the numerical model in which the temperature-dependent material properties are taken into consideration. The moving wheel-rail contact force and the frictional heating are used to simulate the wheel rolling on the rail. The effect of the creepage on the temperature rise, thermal strain, residual stress and residual strain under wheel-rail sliding-rolling contact are investigated. The investigation results show that the thermally affected zone exists mainly in a very thin layer of material near the rail contact surface during the rolling-sliding contact. Both the temperature and thermal strain of rail increase with increasing creepage. The residual stresses induced by the frictional heat in the surface layer of rail appear to be tensile. When the creepage is large, the frictional heat has a significant influence on the residual stresses and residual strains of rail. This paper develops a thermo-meehanical coupling model of wheel-rail rolling-sliding contact, and the obtained results can help to understand the mechanism of wheel/rail frictional thermal fatigue.

NUMERICAL SIMULATION FOR THE SIMPLE SHEAR TEST IN AMORPHOUS GLASSY POLYMERS

In this paper , a driving stress finite element method ofelastic-plastic large deformation based on implicit time integratingalgorithm and an eight-chain molecular network model is used for thenumerical simulation of the simple shear test ofpolycarbonate(PC)materials. The simulated results are compared withexperimental ones. The strain localiztaion propagation for the shearband deforma- tion for simple shear deformation is investigatednumerically. The effects of microstructure parameters in the model onstrain softening and orientation hardening of the PC are discussed indetail.