Open access

# 3D Nonlinear Finite Element Plastic Analysis of Cylindrical Vessels Under In-Plane Moment

Written By

B.H. Wu

Submitted: 02 December 2011 Published: 10 October 2012

DOI: 10.5772/46166

From the Edited Volume

## Finite Element Analysis - Applications in Mechanical Engineering

Chapter metrics overview

View Full Metrics

## 1. Introduction

Cylindrical vessels with nozzles are common structural components in many industries, such as power engineering, petrochemical, etc. Under the applied pressure and piping loads, a high stress concentration is caused by the geometric discontinuity. The connection region of the vessel and nozzle will become the weakest location of the entire structure. Therefore, it is necessary to have an accurate design method for the structure. The plastic limit design method is one possibility. Indesign, a gross plastic deformation is prevented by restricting the allowable load relative to the plastic limit load of the vessels. Therefore, the key step in limit design is to determine the plastic limit load of the vessels under different loads.

The objective of this paper is to determine the plastic limit moment by both experiment and finite element analysis for cylindrical vessels under in-plane moment loading on the nozzle. Based on these results, a parametric analysis is carried out and an empirical formula is proposed.

## 2. Experimental study

### 2.1. Model vessels

Three model vessels with different d/D ratios were designed and fabricated for the experimental study. Every model vessel consisted of a cylinder, nozzle, and flanges for fixing and loading bars. Figure 1 shows a typical configuration, and dimensions for the model vessels are listed in Table 1.

### Table 1.

Dimensions of model vessels

The materials of the cylinder and the nozzle are Q235-A (low carbon steel, similar to A36-77) and 20# (low carbon steel, similar to A106-80 GrA), respectively. Detailed chemical composition and mechanical properties of materials are given in Reference [14]. Figure 2shows the engineering stress-strain curve of materials for Q235-A and 20# Steel.

### 2.2. Experimental setup

The local strains and nozzle displacements under in-plane moment on nozzle were measured to obtain the elastic stresses and deformation characteristics of the model vessels.The in-plane moment is applied as a force at the end of the nozzle, and the moment arm is the distance to the surface of the vessel.

Strain values at the typical measurement points of the cylinders and nozzles were measured using two element strain rosettes. Strain gages were mounted in the axial direction on the outside and inside surfaces for the cylinders, and on the outside surface only for the nozzles. Figure 3 illustrates the detailed locations of the strain gages for test vessel No. L2.

Displacement values under in-plane moment (longitudinal loading direction) for the selected measurement points on the nozzle were measured using mechanical displacement sensors to obtain the load-displacement relationship in the elastic and plastic stages of the nozzles. A total of three displacement sensors were installed at locations along the nozzle length direction. The locations of the displacement sensors are indicated in Figure1.

Figure 4is a photograph of test vessel No. L1during the test.The figure shows an obvious longitudinal deformation of the nozzle under the in-plane moment on the nozzle.

### 2.3. Experimentalresults

#### 2.3.1. Plastic deformation behavior

Figure 5 indicates the elastic-plastic load deformation response at the three measurement points on the nozzle for model vesselNo. L2.

The intersection area of the nozzle and cylinder also produced an evident plastic deformation: the root of the nozzle moved into the vessel wall in the radial direction; as the cylinder yielded on the compression side of the cylinder and nozzle. Figure 6 shows an actual deformation state for test vessels No.L1and No.L2 under in-plane moment of 4.2kN.m and 7.5kN.m.

#### 2.3.2. Plastic limit moment

Experimental plastic limit moments were obtained by the use of load versus displacement plots of the measurement points on the nozzle and load against strain curves of the key gauges located near the junction of the cylinder and nozzle. The plastic limit load is defined by applying the twice elastic slope criterion provided by the ASME Boiler and Pressure Vessel Code[1].

Figure 7 illustrates some typical load-displacement curves and corresponding limit moment for the sensor No.3.

Figure 8 shows some typical load-strain curves and relevant limit moment of thestrain gage No.6 as an example for the experiment vessels.

Table 2 is a summary and comparison of the experimental plastic limit moment for the three model vessels. From the table, it is seen that the limit moments obtained by two measurement methods (displacement measuring and strain measuring) are in agreement. The maximum difference between the limit moments obtained from the two methods is 9.76%. For the same model vessel, the limit moments obtained using different measuring points also are in good agreement regardless ofwhether the measuring method is displacement measuring or strain measuring.

### Table 2.

Experimental results of plastic limit moment (kN.m)

## 3. Finite element analysis

### 3.1. Models and mesh

A static nonlinear finite element analysis of the experimental model vessels was carried out using the ANSYS code [15].Three-dimensional isoperimetric solid elements defined by eight nodes were used to generate the FEA mesh of the cylinder, nozzle and weld seam between the cylinder and nozzle. Two elements were used across the thickness of the shell.Due to the symmetry of the structure and loading, only one-half of the test vessel was modeled. Figure 9 illustrates the finite element mesh of the model vessel No.L2.

The boundary conditions for the numerical simulation are the same as those used in the experimental portion of this study. One end of the cylinder is fixed, while the other end is free. All nodes on the symmetry section (longitudinal plane of the vessel) are constrained against displacement in the direction perpendicular to the symmetry plane. The FEA models were loaded a force at the end of the nozzle along the longitudinal direction of the vessels like in the experiment. For the purpose of comparing the results of the element analysis with those of the experiments, the same materials and loading increments as those of the test vessels were used.As a nonliner geometry analysis, the multilinear isotropic type of hardening behavior is used.A multi-linear elastic-plastic material model defined by nine points, as given in Table 3, was employed.

Material model

### 3.3. Analysis results

#### 3.3.1. Plastic deformation

For the same loading condition as the experiment, the nozzle produced an obvious bending deformation and a depression of the cylinder on thecompression side (in the longitudinal section)occurs. Figure 10 shows the local deformation characteristics which are consistent with those of the experiment (see Figure 4) for model vessel No.L2.

#### 3.3.2. Plastic limit moment

Figure 11 indicates the load-displacement plots of simulation vessels from the sensor No.3 and corresponding limit moments on the nozzles.

At the same location as that of the experiment, loadstrain curves were plated in Figure 12, and corresponding limit moments for simulation vessels of strain gage No.6 were also shown in Figure 12.

The plastic limit moments for the three model vessels by FEA are listed in Table 4.

### Table 4.

FEA Results for the Plastic Limit Moment (kN.m)

## 4. Summary and comparison

From Tables 2 and 4, it is seen that the plastic limit moments from the displacement measurements are consistent with those determined from the strain measurements, whether by experimental method or finite element analysis.

Using an average of the values from different measuring points and displacement sensors, a summary and comparison of the twice-elastic-slope plastic limit load by experiment and finite element analysis for the three model vessels are shown in Table 5.

### Table 5.

Summary and comparison of plastic limit moment

From Table 5, it can be seen that the results for the plastic limit moment from experiment and finite element analysis including displacement and strain measuring techniques are in good agreement. The difference between experiment and finite element analysis is within 10%.

## 5. Parametric analysis and correlation equation

### 5.1. Parametric finite element modeling

On the basis of the previous studies, parametric modeling was performed. The sets of parametersof the analysis models used in this part of the study are given in Table 6. A total of 64 configurations were analyzed to investigate the relationship between the various geometric parameters and plastic limit moments of the vessel-nozzle structures. The data in Table 7 were obtained from the finite element analysisof 64 models.

### 5.2. Correlation equation

A correlation equation for plastic limit load of cylindrical vessels under in-plane moment on the nozzle can be obtained by using the data from Table 7 and the software package Statistica(non-linear regression).The resulting correlation equation for the plastic limit moment on the nozzle is as follows:

MiL=[0.0148+0.704(dD)+1.845(dD)2](dD)1.761(tT)0.555(DT)0.476MbE1

Where, Mb is the limit load of a straight nozzle,

Mb=d2tσs[8](kN.m)E2

This equation holds for the following range of parameters:

0.1d/D0.40.25t/T150D/T125500mmD2000mmE3

### Table 6.

Primary parameters of the FEmodels

### Table 7.

Limit loads of the FE models(kN.m)

In order to verify the accuracy of the correlation equation, the plastic limit moments of three models were calculated and compared with those from the experimental and FEA studies as seen in Table 8.

### Table 8.

Comparison of plastic limit moment(kN.m)

## 6. Conclusions

Experiments and comparative 3D-nonlinear finite element analysesof cylindrical vessel-nozzle connections subjected to an in-plane moment on the nozzle were carried out. In addition, an extensive geometric parameter study of such joints was performed with FEA.The following general conclusions can be drawn from these studies:

1. The experimental and 3D-nonlinear finite element analyses including load-strain and load-displacement of the nozzle responses provide effective and reliable determination methods for the pressure vessel-nozzle (or main-branch pipe) connections subjected to external loads on the nozzle.

2. The plastic limit moments determined from experimental and finite element analysis studies are in good agreement. The maximum error of the average experimental versus FE calculated limit moments is about 10%.

3. The limit momentsfrom load-displacement and load-strain curves are in reasonable agreement for the analyzed geometry.

4. The plastic limit moments from experimental and finite element analysis studies are greatly increased with an increase of d/D, t/T ratios of the model vessels.

5. The correlation equation provided in this paper can be used to determine,within reasonable limits, the plastic limit moment of cylinder-nozzle connections subjected to an in-plane load on the nozzle.

## References

1. 1. ASME Boiler and Pressure Vessel Code,2004edition, section VIII, Division 2, American Society of Mechanical Engineering, New York.
2. 2. Ellyin,F.,1976, “Experimental Investigation of Limit Loads of Nozzles in Cylindrical Vessels,”Welding Research Council Bulletin,No.219, WRC, New York.
3. 3. Ellyin,F.,1977, “An Experimental Study of Elastic-plastic Response of Branch Pipe Tee Connection Subjected to Internal Pressure, External Couples and Combined Loadings,” Welding Research Council Bulletin,No.230, WRC, New York.
4. 4. Schroeder,J.,1985, “Experimental Limit Couples for Branch Moment Loads on 4 -in.ANSI B16.9 Tees,”Welding Research Council Bulletin,304WRC, New York.
5. 5. Junker,A.T.,1985, “Finite Element Determination of the Limit Load of a Nozzle in a Cylindrical Vessel Due to In-plane and Out-of-plane Moments,” Swanson Engineering Associates Corporation (SEAC) Report No.SEAC-TR-312,Rev.2, to the Pressure Vessel Research Committee, WRC, New York.
6. 6. Moffat,D.G.,1985Experimental Stress Analysis of Four Fabricated Equal Diameter Branch Pipe Intersections Subjected to Moment Loadings and Implications on Branch Junction Design,”ProcInstn. Mech.Engrs., 4261
7. 7. Moffat,D.G.,1991, Mwenifumbo, J.A.M.,Xu,S.H. and Mistry,T., “Effective Stress Factors for Piping Branch Junctions Due to Internal Pressure and External Moment Loads,” Journal of Strain Analysis, Vol.26,No.2,PP.85~101.
8. 8. Rodabaugh,E.C.,1988, “A Review of Area Replacement Rules for Pipe Connections in Pressure Vessels and Piping,”Welding Research Council Bulletin,No.335, WRC, New York.
9. 9. Rodabaugh, E.C., Interpretive report on limit analysis and plastic behavior of piping products, Welding Research Council Bullition No.254,1976
10. 10. Kalnins,A.,2001, “Guidelines for Sizing of Vessels by Limit Analysis,”Welding Research Council Bulletin,No.464, WRC, New York.
11. 11. ChapuliotS.MoulinD.PlancgPlancg,D.,2002Mechanical Behavior of a Branch Pipe Subjected to Out-of-plane Bending LoadASME Journal of Pressure Vessel Technology1241
12. 12. MouradHashem. M.MaherY.Y.A.,2002Limit-Load analysis of Pipe Bend under Out-of-plane Moment and Internal Pressure,” ASME Journal of Pressure Vessel Technology, 1247
13. 13. Widera,G.E.O.,Wei,Z.,1998, “Parametric Finite Element of Large Diameter Shell Intersection,” Part 2, External Loading, PVRC Project No.96-20AS, WRC, New York.
14. 14. Sang,Z.F.,Wang,Z.L.,Xue,L.P.,Widera,G..E.O.,2005, “Plastic Limit Loads of Nozzles in Cylindrical Vessels under Out-of-plane Moment Loading,” International Journal of Pressure Vessels and Piping, 82,PP.638-648.
15. 15. Swanson Analysis System Inc., “ANSYS Engineering Analysis Systems User’s Manual”.

Written By

B.H. Wu

Submitted: 02 December 2011 Published: 10 October 2012