and test its strength. The tested results are consistent
with the finite element simulation results.
2MODELING OF THE EXTRUSION
CYLINDER
2.1 Cylinder Material Model
The stress-strain curve of the material is shown in
Figure 1. It clearly demonstrates the Bauschinger-
effect. Analysis procedure for autofrettage process
involves overstrain which is typically based on
Tresca or Von. Mises failure criteria. Here, using a
bilinear isotropic hardening model approximating
the real material behaviour is shown in Figure 1; has
been used in which E1 is the slope of the linear line
in the tensile elastic region (modulus of elasticity),
H1 is the slope of the linear line in the tensile plastic
region, E2 is the slope of the linear line in the
compression elastic region (modulus of elasticity),
and H2 is the slope of the linear line in the
compression plastic region. This material’s constants
are as follows: E1=206GPa; H1=2GPa; E2=207GPa;
H2=57GPa; σy1 =1103 MPa, σy2
=1100 MPa, where, v, σy1, and σy2 are the
Poisson's ratio , tensile yield stress, compression
yield stress, respectively.
300
600
900
1.20.6
1.8
2.4
3.0
-0.6
-300
-1200
Strain / %
-600
-900
uniaxial tension-compression
Fitting model
Calculating model
E
1
E
2
H
1
H
2
Figure 1: Material strain-stress curve.
2.2 Finite Element Model
The finite element model of TSEC has been
constructed in ANSYS APDL 15.0. The geometric
model is 1/4 of circumference of the cylinder, and
the element used here is 3D SOLID 186, and
mapped mesh method be used for easy simulation of
boring bore.
The cylinder body is mainly subjected to radial
force. During shrinking and autofrettage, the contact
surfaces are all one-way contact behaviors between
the surfaces, and the material of the three-layer
cylinder is also the same, so the contact can be
regarded as the “surface to surface” contact model of
the “Flexible-Flexible”. The target surface is the
inner surface of the middle layer and the outer layer;
The element type used is Target170. The contact
surface is the outer surface and middle surface, the
element type is Contact174. The geometric
dimensions of every cylinders was the actual
dimensions, including the initial mutual penetration,
and penetration tolerance was 0.0001 mm, and the
contact algorithm was chosen to be the Augmented
Lagrangian Method.
2.3 Boundary Conditions
Analytical model was the 1/4 of the extrusion
cylinder, so it is necessary to apply symmetry B. C.
on the corresponding surface. For easy solution,
assuming that the extrusion cylinder was in the plane
strain state, then axial displacement constraints has
been applied the both ends of cylinder.
3 SIMULATION AND RESULT
ANALYSIS
In the ANSYS environment, a number of
simulations have been performed on the shrink-fit,
autofrettage, boring and loading processes of
cylinder to calculate the residual stress distribution
in the shrink-fitted and multiple autofrettaged
compound cylinder and the maximum elastic
pressure capacity. Here the autofrettage pressure was
1800 MPa (pressurized system up to maximum
pressure).
First, having simulated the shrink-fit process of
the three-layer cylinder, residual hoop stress(RHS) is
shown in Figure 2(a). It reveals that the inner layer is
subjected to compressive stress with a large
value,and stress at the inner bore area (working area)
is 1150.94MPa, and the inner part of middle cylinder
is subjected to tension, and the outer part is
subjected to compression, and the stress value is
small, and the outer layer is subjected tension. Mises
stress is shown in Figure 2(b), and the maximum
stress is 1028.67MPa, and it is lower than
compressive yield limit of material, on this
condition, the entire cylinder is in elastic state.