A New Classification Method of Stress Modes in Hybrid Stress Finite Element

The paper presents a new method for classifying the stress modes in hybrid stress finite element in terms of natural stress modes in finite element and the rank analysis of matrix G in forming element It reveals the relation among the different assumed stress field, and gives the general method in forming stress field Comparing with the method of eigenvalue analysis, the new method is more efficient

THE STRESS SUBSPACE OF HYBRID STRESS ELEMENT AND THE DIAGONALIZATION METHOD FOR FLEXIBILITY MATRIX H

The following is proved: 1) The linear independence of assumed stress modes is the necessary and sufficient condition for the nonsingular flexibility matrix; 2) The equivalent assumed stress modes lead to the identical hybrid element. The Hilbert stress subspace of the assumed stress modes is established. So, it is easy to derive the equivalent orthogonal normal stress modes by Schmidt's method. Because of the resulting diagonal flexibility matrix, the identical hybrid element is free from the complex matrix inversion so that the hybrid efficiency, is improved greatly. The numerical examples show that the method is effective.

A quadrilateral membrane hybrid stress element with drilling degrees of freedom

A new kind of quadrilateral assumed stress hy- brid membrane element with drilling degrees of freedom and a traction-free inclined side has been developed based on an extended Hellinger-Reissner principle which is established by expanding the essential terms of the assumed stress field as polynomials in the natural coordinates of the element. The homogeneous equilibrium equations are imposed in a variational sense through the internal displacements which are also expanded in the natural coordinates, while the tractionfree conditions along the inclined side are satisfied exactly. The use of such special element in the finite element solution is shown to be highly accurate when only a very coarse element mesh is used for plates with V-shaped rounded notches or inclined sides.

Determination of stress intensity factor with direct stress approach using finite element analysis

In this article, a direct stress approach based on finite element analysis to determine the stress intensity factor is improved. Firstly, by comparing the rigorous solution against the asymptotic solution for a problem of an infinite plate embedded a central crack, we found that the stresses in a restrictive interval near the crack tip given by the rigorous solution can be used to determine the stress intensity factor, which is nearly equal to the stress intensity factor given by the asymptotic solution. Secondly, the crack problem is solved numerically by the finite element method. Depending on the modeling capability of the software, we designed an adaptive mesh model to simulate the stress singularity. Thus, the stress result in an appropriate interval near the crack tip is fairly approximated to the rigorous solution of the corresponding crack problem. Therefore, the stress intensity factor may be calculated from the stress distribution in the appropriate interval, with a high accuracy.

COMPUTATION OF SUPER-CONVERGENT NODAL STRESSES OF TIMOSHENKO BEAM ELEMENTS BY EEP METHOD

The newly proposed element energy projection(EEP) method has been applied to the computation of super_convergent nodal stresses of Timoshenko beam elements.General formulas based on element projection theorem were derived and illustrative numerical examples using two typical elements were given.Both the analysis and examples show that EEP method also works very well for the problems with vector function solutions.The EEP method gives super_convergent nodal stresses,which are well comparable to the nodal displacements in terms of both convergence rate and error magnitude.And in addition,it can overcome the “shear locking” difficulty for stresses even when the displacements are badly affected.This research paves the way for application of the EEP method to general one_dimensional systems of ordinary differential equations.

A hybrid-stress element based on Hamilton principle

A novel hybrid-stress finite element method is proposed for constructing simple 4-node quadrilateral plane elements, and the new element is denoted as HH4-3fl here. Firstly, the theoretical basis of the traditional hybrid-stress elements, i.e., the Hellinger-Reissner variational principle, is replaced by the Hamilton variational principle, in which the number of the stress variables is reduced from 3 to 2. Secondly, three stress parameters and corresponding trial functions are introduced into the system equations. Thirdly, the displacement fields of the conventional bilinear isoparametric element are employed in the new models. Finally, from the stationary condition, the stress parameters can be expressed in terms of the displacement parameters, and thus the new element stiffness matrices can be obtained. Since the required number of stress variables in the Hamilton variational principle is less than that in the Hellinger-Reissner variational principle, and no additional incompatible displacement modes are considered, the new hybrid-stress element is simpler than the traditional ones. Furthermore, in order to improve the accuracy of the stress solutions, two enhanced post-processing schemes are also proposed for element HH4-3β. Numerical examples show that the proposed model exhibits great improvements in both displacement and stress solutions, implying that the proposed technique is an effective way for developing simple finite element models with high performance.

3-D FINITE ELEMENT ANALYSIS OF THE EFFECT OF WELDING RESIDUAL STRESS ON HYDROGEN DIFFUSION IN HYDROGEN CONTAINED ENVIRONMENT

The hydrogen distribution of 16MnR steel weldment in hydrogen contained environment was calculated using the finite element method （ FEM）. The effect of welding residual stress on hydrogen diffusion has been discussed using a 3-D sequential coupling finite element analysis procedure complied by Abaqus code. The hydrogen diffusion coefficient in weld metal, the heat affected zone （HAZ）, and the base metal of the 16MnR steel weldment were measured using the electrochemical permeation technique. The hydrogen diffusion without the effect of stress was also calculated and compared. Owing to the existence of welding residual stress, the hydrogen concentration was obviously increased and the hydrogen wouM diffuse and accumulate in the higher stress region.

3D Finite Element Numerical Simulation of Residual Stresses on Electron Beam Welded BT20 Plates

A three-dimensional finite-element model (FEM) used for calculating electron beam (EB) welding temperature and stresses fields of thin plates of BT20 titanium has been developed in which the nonlinear thermophysical and thermo-mechanical properties of the material has been considered. The welding temperature field, the distributions of residual stresses in as-welded (AW) and electron beam local post-weld heat treatment (EBLPWHT) conditions have been successfully simulated. The results show that: (1) In the weld center, the maximum magnitude of residual tensile stresses of BT20 thin plates of Ti alloy is equal to 60%- 70% of its yield strength σs. (2) The residual tensile stresses in weld center can be even decreased after EBLPWHT and the longitudinal tensile stresses are decreased about 50% compared to joints in AW conditions. (3) The numerical calculating results of residual stresses by using FEM are basically in agreement with the experimental results. Combined with numerical calculating results, the effects of electron beam welding and EBLPWHT on the distribution of welding residual stresses in thin plates of BT20 have been analyzed in detail.

Orthogonal basic deformation mode method for zero-energy mode suppression of hybrid stress elements

A set of basic deformation modes for hybrid stress finite elements are directly derived from the element displacement field. Subsequently, by employing the so-called united orthogonal conditions, a new orthogonalization method is proposed. The result- ing orthogonal basic deformation modes exhibit simple and clear physical meanings. In addition, they do not involve any material parameters, and thus can be efficiently used to examine the element performance and serve as a unified tool to assess different hybrid elements. Thereafter, a convenient approach for the identification of spurious zero-energy modes is presented using the positive definiteness property of a flexibility matrix. More- over, based on the orthogonality relationship between the given initial stress modes and the orthogonal basic deformation modes, an alternative method of assumed stress modes to formulate a hybrid element free of spurious modes is discussed. It is found that the orthogonality of the basic deformation modes is the sufficient and necessary condition for the suppression of spurious zero-energy modes. Numerical examples of 2D 4-node quadrilateral elements and 3D 8-node hexahedral elements are illustrated in detail to demonstrate the efficiency of the proposed orthogonal basic deformation mode method.

DEFORMATION RIGIDITY OF ASSUMED STRESS MODES IN HYBRID ELEMENTS

The new methods to determine the zero-energy deformation modes in the hybrid elements and the zero-energy stress modes in their assumed stress fields are presented by the natural deformation modes of the elements. And the formula of the additional element deformation rigidity due to additional mode into the assumed stress field is derived. Based on, it is concluded in theory that the zero-energy stress mode cannot suppress the zero-energy deformation modes but increase the extra rigidity to the nonzero-energy deformation modes of the element instead. So they should not be employed to assume the stress field. In addition, the parasitic stress modes will produce the spurious parasitic energy and result the element behaving over rigidity. Thus, they should not be used into the assumed stress field even though they can suppress the zero-energy deformation modes of the element. The numerical examples show the performance of the elements including the zero-energy stress modes or the parasitic stress modes.

Finite-element 3D modeling of stress patterns around a dipping fault

Stresses in a block around a dipping fracture simulating a damage zone of a fault are reconstructed by finite-element modeling. A fracture corresponding to a fault of different lengths, with its plane dipping at different angles, is assumed to follow a lithological interface and to experience either compression or shear. The stress associated with the destruction shows an asymmetrical pattern with different distances from the highest stress sites to the fault plane in the hanging and foot walls. As the dip angle decreases,the high-stress zone becomes wider in the hanging wall but its width changes negligibly in the foot wall.The length of the simulated fault and the deformation type affect only the magnitude of maximum stress,which remains asymmetrical relative to the fault plane. The Lh/Lfratio, where Lhand Lfare the widths of high-stress zones in the hanging and foot walls of the fault, respectively, is inversely proportional to the fault plane dip. The arithmetic mean of this ratio over different fault lengths in fractures subject to compression changes from 0.29 at a dip of 80°to 1.67 at 30°. In the case of shift displacement, ratios are increasing to 1.2 and 2.94, respectively.Usually they consider vertical fault planes and symmetry in a damage zone of faults. Following that assumption may cause errors in reconstructions of stress and fault patterns in areas of complex structural setting. According geological data, we know the structures are different and asymmetric in hanging and foot walls of fault. Thus, it is important to quantify zones of that asymmetry. The modeling results have to be taken into account in studies of natural faults, especially for practical applications in seismic risk mapping, engineering geology, hydrogeology, and tectonics.

3-D finite element modeling for evolution of stress field and interaction among strong earthquakes in Sichuan-Yunnan region

Based on the latest achievement about activities of geological structure, a 3-D finite-element model containing four layers of upper crust, lower crust (two layers) and upper mantle is established in the paper. By repeated tests and revisions, the boundary conditions of the model are determined. And then the background stress field, the stress field caused by fault creep and the stress field triggered by strong earthquake in Sichuan-Yunnan region, as well as their dynamic variations are calculated. The results indicate that the latter earthquake often occurs in the area with positive Coulomb rupture stress change associated with the former one, the former earthquake has a triggering effect on the latter one to a certain extent, and strong earthquake often occur in groups under the background of high stress, which is of great significance for distinguishing seismic anomalies, as well as for improving the level of earthquake prediction.

A Finite Element Solution for Barrel Dynamic Stress

With the APDL language of ANSYS finite element analysis software, the solution program for barrel dynamic stress is developed. The paper describes the pivotal problems of dynamic strength design and provides a foundation for realizing the engineering and programming of barrel dynamic strength design.

CONTACT FRICTION ANALYSIS AND STRESS OSCILLATION SUPPRESSION WITH A SIMPLE INTERFACE ELEMENT

A simple interface element for analyzing contact friction problems is developed. Taking nodal displacements and contact stresses as unknowns, this element can simulate frictional slippage, decoupling and re-bonding of two bodies initially mating or having gaps at a common interface. The method is based on the Finite Element Method and load incremental theory. The geometric and static constraint conditions on contact surfaces are treated as additional conditions and are included in stiffness equations. This simple element has the advantages of easy implementation into standard finite element programs and fast speed for convergence as well as high accuracy for stress distribution in interface. Undesirable stress oscillations are also investigated whenever large stress gradients exist over the contact surfaces. Exact integration or the conventional Gauss integration scheme used to evaluate the interpolation function matrix of the interface element is found to be the source of the oscillations. Eigenmode analysis demonstrates that the stress behavior of an interface element can be improved by using the Newton-Cotes integration scheme. Finally, the test example of a strip footing problem is presented.

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.

Finite Element Analysis of Thermal Stresses in Ceramic/Metal Gradient Thermal Barrier Coatings

This paper studied the thermal stresses of ceramicl metal gradient thermal barrier coating which combines the conceptions of ceramic thermal barrier coating （TBC） and functionally gradient material （FGM）. Thermal stresses and residual thermal stresses were calculated by an ANSYS finite element analysis software. Negative thermal expansion coefficient method was proposed and element birth and death method was applied to analyze the residual thermal stresses which have non-uniform initial temperature field. The numerical results show a good agreement with the analytical results and the experimental results.

A 3D finite element modeling on the relationship between crustal stress and earthquake activity in Liaoning and its vicinity

Based on the latest result in research on 3D seismic wave velocity structure of crust and uppermost mantle and taking geological setting and fracture zones into consideration, a 3D geological model for the studied region is built up. The boundary constraint and force loading boundary condition for the model are determined according to the characteristics of crustal stress field deduced from earthquake focal mechanism and in-situ stress measurement data. Using linear elastic material model a 3D finite element modeling is conducted to study the characteristics of crustal stress field. A comparison analysis between the simulated stress field and earthquake locations reveals that the moderate and strong earthquakes generally occurred in the zones with high shear stress gradient. Furthermore, the paper notices a few potential earthquake-prone regions.

Autologous nerve graft repair of different degrees of sciatic nerve defect:stress and displacement at the anastomosis in a three-dimensional finite element simulation model

In the repair of peripheral nerve injury using autologous or synthetic nerve grafting, the mag- nitude of tensile forces at the anastomosis affects its response to physiological stress and the ultimate success of the treatment. One-dimensional stretching is commonly used to measure changes in tensile stress and strain; however, the accuracy of this simple method is limited. There- fore, in the present study, we established three-dimensional finite element models of sciatic nerve defects repaired by autologous nerve grafts. Using PRO E 5.0 finite element simulation software, we calculated the maximum stress and displacement of an anastomosis under a 5 N load in 10-, 20-, 30-, 40-mm long autologous nerve grafts. We found that maximum displacement increased with graft length, consistent with specimen force. These findings indicate that three-dimensional finite element simulation is a feasible method for analyzing stress and displacement at the anas- tomosis after autologous nerve grafting.

Finite element analysis of stresses and interface crack in TBC systems

Thermal barrier coatings have been used on high temperature components. Due to high stresses leading to unpredictable failure, a transient thermal-structural finite element solution was employed to analyze the stress distribution and J-integral at the interface between the bond coating and thermally growing oxide(TGO) in the EB-PVD thermal barrier coatings subjected to thermal loadings. The effects of some environmental and material parameters were studied, such as thermal convection coefficient, ceramic elastic modulus and TGO thickness. The results show that the stresses and J-integral values are impacted by these parameters.

SINGULAR INTEGRAL EQUATIONS AND BOUNDARY ELEMENT METHOD OF CRACKS IN THERMALLY STRESSED PLANAR SOLIDS

Using the method of the boundary integral equation, a set of singular integral equations of the hear transfer problems and the thermo-elastic problems of a crack embedded in a two-dimensional finite body is derived, and then,its numerical method is proposed by the numerical method of the singular integral equations combined with boundary element method. Moreover, the singular nature of temperature gradient field near the crack front is proved by the main-part analysis method of the singular integral equation, and the singular temperature gradients are exactly obtained. Finally, several typical examples calculated.