Fatigue PDF
Fatigue PDF
CONTENTS
1 INTRODUCTION
2.2 FEM MODEL & BOUNDARY CONDITION FOR THERMAL ANALYSIS AROUND TOP DISHED
END
2.4 FEM MODEL AND BOUNDARY CONDITION FOR FATIGUE STRESS ANALYSIS AROUND TOP
DISHED END
2.5 RESULTS FROM FATIGUE STRESS ANALYSIS AROUND TOP DISHED END AND FATIGUE
LIFE CALCULATION
2.6 FEM MODEL & BOUNDARY CONDITION FOR THERMAL ANALYSIS AROUND BOTTOM
DISHED END
2.8 FEM MODEL AND BOUNDARY CONDITION FOR FATIGUE STRESS ANALYSIS AROUND
BOTTOM DISHED END
2.9 RESULTS FROM FATIGUE STRESS ANALYSIS AROUND BOTTOM DISHED END AND
FATIGUE LIFE CALCULATION
3.2 FEM MODEL AND BOUNDARY CONDITIONS FOR STRESS ANALYSIS AROUND NOZZLE
A1/N1
3.4 FEM MODEL AND BOUNDARY CONDITIONS FOR STRESS ANALYSIS AROUND NOZZLE N2
3.8 FATIGUE ANALYSIS FOR INTERNAL DIS TRIBUTOR PIPE FOR NOZZLE N1 & N2
4 CONCLUSIONS
5 REFERENCES
LIST OF TABLES
LIST OF FIGURES
FIGURE 3.DISPLACEMENT PROFILE FOR TOP DISHED END (THERMAL AND PRESSURE LOAD )
FIGURE 4.MAXIMUM SHEAR STRESS CONTOUR FOR TOP DISHED END (THERMAL AND PRESSURE
LOAD)
FIGURE 7.DISPLACEMENT PROFILE FOR BOTTOM DISHED END (THERMAL AND PRESSURE LOAD
)
FIGURE 8.MAXIMUM SHEAR STRESS CONTOUR FOR BOTTOM DISHED END (THERMAL AND
PRESSURE LOAD)
FIGURE 9.DISPLACEMENT PROFILE FOR TOP DISHED END (THERMAL, PRESSURE AND NOZZLE
LOAD)
FIGURE 10.MAXIMUM SHEAR STRESS CONTOUR FOR TOP DISHED END (THERMAL , PRESSURE
AND NOZZLE LOAD)
FIGURE 11.DISPLACEMENT PROFILE FOR BOTTOM DISHED END (THERMAL , PRESSURE AND
NOZZLE LOAD)
-4-
NOVEMBER 2001
FIGURE 12.MAXIMUM SHEAR STRESS CONTOUR FOR BOTTOM DISHED END (THERMAL ,
PRESSURE AND NOZZLE LOAD)
FIGURE 13.MAXIMUM SHEAR STRESS CONTOUR FOR SKIRT TO BOTTOM DISHED END JOINT
-5-
NOVEMBER 2001
1 INTRODUCTION
Three number of Adsorbers (4C-1301 A/B/C) are manufactured for WOODSIDE ENERGY
LIMITED, AUSTRALIA as per AS1210 -1997. Each of the equipment is vertical with inside
diameter of 4000 mm, shell thickness of 109 mm and top / bottom dished end thickness of 58 mm.
Each Adsorber undergoes 1350 minute operating cycle consisting of temperature ramp-up/ramp-
down in the range of 21°C to 320°C, and corresponding internal pressure variation from 5780 kPaa
to 5370 kPaa .The temperature fluctuation take place in not less than 20 minutes time. Adsorbers
Mechanical design of vessel, nozzles is required to consider the effects of above mentioned
• Fatigue life calculations for vessel shell, dished ends, nozzles N1/A1 (36”)&N2 (24”) for above
thermal, pressure cycle using FEA module of NISA version 9.0 software.
• Stress check for nozzles N1/A1 & N2 for thermal, internal pressure and nozzle local loads
• Transient thermal analysis to arrive at time step corresponding to severest temperature gradient
• Peak stress and fatigue life calculation for above time step.
-6-
NOVEMBER 2001
The maximum temperature ramp up at the regeneration gas heater outlet is from 21°c to 320°c
in not less than 0.33hr and vice versa for maximum temperature ramp down. The internal
Nozzle local loads are varying as per the following. These loads are assumed to include the
thermal movement of piping attached & as a result no additional constraint are required.
-7-
NOVEMBER 2001
2.2 FEM MODEL & BOUNDARY CONDITION FOR THERMAL ANALYSIS AROUND TOP DISHED
END
Top section of the vessel including nozzle N1/A1, dished end and shell portion is modeled as shown
in Fig.1. All dimensions are taken in corroded condition. The model is having 1000 second order,
Heating case: inside gas temperature of 320°C and initial vessel temperature of 21°C.
Cooling case: inside gas temperature of 21°C and initial vessel temperature of 320°C.
• Convective film heat transfer coefficients on inside surfaces is assumed as follows. These are
320°C
21°C
Fig.2 Shows temperature distribution around top dished end at time step of 1300 seconds from
vessel.
2.4 FEM MODEL AND BOUNDARY CONDITION FOR FATIGUE STRESS ANALYSIS AROUND TOP
DISHED END
Same FEA model as that used for thermal analysis is used for stress analysis. The boundary
• Separate runs are taken with internal pressure of 5780 kPaa and 5370 kPaa is applied at 1300
2.5 RESULTS FROM FATIGUE STRESS ANALYSIS AROUND TOP DISHED END AND FATIGUE
LIFE CALCULATION
Fig.3 shows typical deformation profile of the model under given loading at 1300 seconds. Fig. 4
shows typical maximum shear stress contours for the model at 1300 seconds.
Since the principle stress direction are varying with time, as per [3], clause SC2.2.2, following table
lists stress ranges at various time steps and at various locations. The stress ranges are calculated
-9-
NOVEMBER 2001
based on stress output at 1300 seconds as base case, as this gives highest combination for stress
ranges.
As shown in above table, the maximum peak stress range in the model due to loading cycle is 25.07
Kg/mm² (245.9 MPa) around inner corner of nozzle A1/N1 & 3.88 Kg/mm² (38.09 MPa) at outer
corner to top dished end junction. To these stresses the fluctuating stress range due to noz. Local
loads (calculated separately) are added. The max. shear stresses due to nozzle local loads are
Alt. Stress range due to nozzle Srij2 Kg/mm²(MPa) 0.47 (4.61) 2.66 (26.09)
loads
Alternating stress amplitude Salt=0.5x Srij Kg/mm²(MPa) 12.77 (125.27) 3.27 (32.09)
- 10 -
NOVEMBER 2001
curve
Corrected alternating stress Salt’ = Salt x Egraph Kg/mm²(MPa) 14.76 (144.85) 3.78 (37.09)
amplitude /E
2.6 FEM MODEL & BOUNDARY CONDITION FOR THERMAL ANALYSIS AROUND BOTTOM
DISHED END
Bottom section of the vessel including nozzle N2, dished end, shell portion and skirt is modeled as
shown in Fig.5. All dimensions are taken in corroded condition .The model is having 2124 second
Heating case: inside gas temperature of 320°C and initial vessel temperature of 21°C.
Cooling case: inside gas temperature of 21°C and initial vessel temperature of 320°C.
- 11 -
NOVEMBER 2001
• Convective film heat transfer coefficients on inside surfaces is assumed as follows. These are
Fig.6 Shows temperature distribution around bottom dished end at time step of 1220 seconds from
vessel.
2.8 FEM MODEL AND BOUNDARY CONDITION FOR FATIGUE STRESS ANALYSIS AROUND
Same FEA model as that used for thermal analysis is used for stress analysis. The boundary
• Separate runs are taken with internal pressure of 5780 kPaa and 5370 kPaa is applied at 1220
- 12 -
NOVEMBER 2001
2.9 RESULTS FROM FATIGUE STRESS ANALYSIS AROUND BOTTOM DISHED END AND
Fig.7 shows typical deformation profile of the model under given loading at 1220 seconds. Fig.8
shows typical maximum shear stress contours for the model at 1220 seconds.
Since the principle stress direction are varying with time, as per [3], clause SC2.2.2, following table
lists stress ranges at various time steps and at various locations. The stress ranges are calculated
based on stress output at 1220 seconds as base case, as this gives highest combination for stress
ranges.
* since peak stress range is exceeding 3 Sm, membrane + bending stress range excluding peak
As shown in above table the maximum peak stress range due to loading cycle is 39.08 Kg/mm²
(383.37 MPa) around inner corner of nozzle N2 & 14.52 Kg/mm² (142.44 MPa) about outer
- 13 -
NOVEMBER 2001
corner and bottom dished end junction. To these stresses the fluctuating stress range due to nozzle
local loads (calculated separately) are added. The max. shear stress due to nozzle local loads are
Alternating stress range due to Srij1 Kg/mm²(MPa) 39.08 (383.37) 14.52 (142.44)
Alt. Stress range due to nozzle Srij2 Kg/mm²(MPa) 1.68 (16.48) 13.6 (133.41)
loads
Alternating stress amplitude Salt=0.5x Srij Kg/mm²(MPa) 20.38 (199.92) 14.06 (137.92)
curve
amplitude h/E
- 14 -
NOVEMBER 2001
As seen from section 2.5 and 2.9 the fatigue life of vessel subjected to thermal cycle as per section
2.1 is observed to be 10,300 cycles, which is more than required 7,787 cycles.
The stresses around nozzles N1/A1 and N2 are checked for thermal, internal pressure and nozzle
local load using FEA. Special axisymmetric elements are used for this analysis, which allow for
application of non-axisymmetric loading. This is used to simulate axial force and circumferential
/longitudinal moments nozzle local loads. The stresses due to shear loads (i.e. longitudinal
/circumferential forces and torsional moment) are calculated separately using classical method.
These two stresses are superimposed algebraically to arrive at maximum stress intensity, which is
- 15 -
NOVEMBER 2001
3.2 FEM MODEL AND BOUNDARY CONDITIONS FOR STRESS ANALYSIS AROUND NOZZLE
A1/N1
Same FEA model as that used for fatigue life calculations is used. The element type used allows
for application of non-axisymmetric loads as described in 3.1 earlier. The boundary conditions used
are as follows-
• Temperature distribution at time step of 1300 seconds as obtained from thermal analysis is
applied on nodes.
• Loads due to FA, ML & MC are applied as equivalent pressure load on nozzle free face.
Fig.9 shows deformation profile of model under given loading.Fig.10 shows maximum shear stress
contour for model. The maximum stress is located at nozzle to dished end junction.
- 16 -
NOVEMBER 2001
Calculations for stresses due to nozzle local loads FL, FC and MT on nozzle A1/N1 are as follows-
3.4 FEM MODEL AND BOUNDARY CONDITIONS FOR STRESS ANALYSIS AROUND NOZZLE N2
Same FEA model as that used for fatigue life calculations is used. The element type used allows
for application of non-axisymmetric loads as described in 3.1 earlier. The boundary conditions used
are as follows-
- 17 -
NOVEMBER 2001
• Temperature distribution at time step of 1200 seconds from the beginning of thermal ramp-
• Loads due to FA, ML & MC are applied as equivalent pressure load on nozzle free face.
Fig.11 shows deformation profile of model under given loading.Fig.12 shows maximum shear stress
contour for model. The maximum stress is located at nozzle to dished end junction.
Calculations for stresses due to nozzle local loads FL, FC and MT on nozzle N2 are as follows-
- 18 -
NOVEMBER 2001
- 19 -
NOVEMBER 2001
The other nozzle stresses are calculated by stress index method given in (article 4-6 of ASME sect.
** Form stress in opening for fatigue evaluation of ASME sec VIII Div. 2 article 4-6
fatigue life calculation for flange and flange bolting due to internal pressure cycle is as given in table
below.
- 20 -
NOVEMBER 2001
& N2 :
The construction of this internal distributor pipe (N1&N2) isolates it from the stresses arising from
the fatigue life calculations for welded attachments are calculated based on cyclic stress range
obtained from FEA results for plain shell. The calculations are given as follows
As given in above table, calculated fatigue life is more than required fatigue life for attachments.
The bed support components are subjected to reversing loading depending on direction of
regeneration cycle. The calculation of stresses for these components is as per method described in
Induced stress range
*Max. stress Stress near non-integral Sa'
Sa Fatigue Required
range concentration (Sa x 207E3 /
4P-1105-
Sr. no. BedMDR Part4. Following table shows fatigue life calculation
component (S induced for these components.
attachments (S
(S range) factor (assumed) induced=S range x
179.02E3)
life life
/2) (Mpa) (Cycles) (Cycles)
(Mpa) (SCF) (Mpa)
SCF) (Mpa)
1 Bed support cleat 115.17 3 345.51 172.75 199.75 12,000 7,787
2 Bed support beam 168.71 1 168.71 84.36 97.54 410,000 7,787
- 21 -
NOVEMBER 2001
Stresses around butt welds in vessel are calculated based on FEA result for plain shell, the stress
range obtained from FEA in plain shell is multiplied by assumed stress concentration factor of 3
(the actual stress concentration factor will be less than 3, however 3 is a conservative assumption)
to arrive at peak stress range around welds. The calculations are shown as below.
As shown in above table these welds can be used in “as welded" condition, without smooth
dressing.
- 22 -
NOVEMBER 2001
4 CONCLUSIONS
In the present analysis, vessel shell, dished ends and nozzle A1/N1 and N2 are checked for cyclic
loading using FEA, for calculating fatigue life. Also the two nozzles A1/N1,N2 and skirt are
checked for combined thermal, pressure and nozzle local loads. The stresses from FEA results are
used to derive fatigue life for other nozzles and attachments. The flange and bolting are also
• The maximum peak stress range is found to be around inner corner of nozzle N2 to bottom
• Fatigue life of the vessel is calculated to be 10,300 cycles, which is more than required 7,787
cycles.
• The stresses around nozzle A1/N1,N2 and skirt for combined pressure, thermal and nozzle local
• The fatigue life for other nozzles, attachments, flanges and bolting is found to be satisfactory.
5 REFERENCES
- 23 -
NOVEMBER 2001
- 24 -
NOVEMBER 2001
- 25 -
NOVEMBER 2001
- 26 -
NOVEMBER 2001
FIGURE 3. DISPLACEMENT PROFILE FOR TOP DISHED END (THERMAL AND PRESSURE
LOAD )
FIGURE 4.
- 27 -
NOVEMBER 2001
FIGURE 4 MAXIMUM SHEAR STRESS CONTOUR FOR TOP DISHED END (THERMAL AND
PRESSURE LOAD)
- 28 -
NOVEMBER 2001
- 29 -
NOVEMBER 2001
- 30 -
NOVEMBER 2001
- 31 -
NOVEMBER 2001
FIGURE 8. MAXIMUM SHEAR STRESS CONTOUR FOR BOTTOM DISHED END (THERMAL
- 32 -
NOVEMBER 2001
FIGURE 9. DISPLACEMENT PROFILE FOR TOP DISHED END (THERMAL, PRESSURE AND
NOZZLE LOAD)
- 33 -
NOVEMBER 2001
FIGURE 10. MAXIMUM SHEAR STRESS CONTOUR FOR TOP DISHED END (THERMAL,
- 34 -
NOVEMBER 2001
FIGURE 11. DISPLACEMENT PROFILE FOR BOTTOM DISHED END (THERMAL, PRESSURE AND
NOZZLE LOAD)
- 35 -
NOVEMBER 2001
FIGURE 12. MAXIMUM SHEAR STRESS CONTOUR FOR BOTTOM DISHED END (THERMAL,
- 36 -
NOVEMBER 2001
FIGURE 13. MAXIMUM SHEAR STRESS CONTOUR FOR SKIRT TO BOTTOM DISHED END
- 37 -
NOVEMBER 2001
- 38 -
NOVEMBER 2001
- 39 -