A Method to Approximate the Steady-State Creep Response of Three-Dimensional Pipe Bend Finite Element Models Under Internal Pressure Loading Using Two-Dimensional Axisymmetric Models

2013 ◽  
Vol 136 (1) ◽  
Author(s):  
J. P. Rouse ◽  
W. Sun ◽  
T. H. Hyde ◽  
A. Morris ◽  
W. Montgomery

Pipe bends are regions of geometric discontinuities in the pipe systems used in power plants and most industry recorded failures have been located around similar regions. Understanding these potential locations of weakness is therefore of great interest for the safe and economic operation of piping components. Increased predictive accuracy would assist in component design, condition monitoring, and retirement strategy decisions. Modeling of piping components for finite element analysis (FEA) is complicated by the variation of the cross section dimensions (changes in wall thicknesses or cross section ovality) around the pipe bend due to the manufacturing procedure implemented. Quantities such as peak rupture stress and creep rupture life can be greatly affected by these geometric variation (Rouse, J. P., Leom, M. Z., Sun, W., Hyde, T. H., Morris, A., “Steady-state Creep Peak Rupture Stresses in 90 Pipe Bends with Manufacture Induced Cross Section Dimension Variations”International Journal of Pressure Vessels and Piping, Volumes 105–106, May–June 2013, pp. 1–11). Three dimensional (3D) models can be used to approximate to the realistic level of detail found in pipe bends. These simulations may however be computationally expensive and could take a considerable amount of time to complete. Two dimensional (2D) axisymmetric models are relatively straight forward to produce and quick to run, but of course cannot represent the full geometric complexity around the pipe bend. A method is proposed that utilises multiple 2D axisymmetric pipe bend models to approximate the result of a 3D analysis through interpolation, thus exploiting the greatly reduced computing time observed for the 2D models. The prediction of peak rupture stress (both magnitude and location) is assessed using a simple power law material model. Comments are made on the applicability of the proposed procedure to a range of bends angles (90 deg, 60 deg, and 30 deg), as well as the effect of the stress exponent (n) and tri-axial (α) material constants. Provided that peak stresses do not occur at the bend/straight interface, the magnitude and location of the peak rupture stress can be predicted by the 2D axisymmetric interpolation method with a typical percentage difference of less than 1%.

Author(s):  
Laura Galuppi ◽  
Gianni Royer-Carfagni

Prandtl's membrane analogy for the torsion problem of prismatic homogeneous bars is extended to multi-material cross sections. The linear elastic problem is governed by the same equations describing the deformation of an inflated membrane, differently tensioned in regions that correspond to the domains hosting different materials in the bar cross section, in a way proportional to the inverse of the material shear modulus. Multi-connected cross sections correspond to materials with vanishing stiffness inside the holes, implying infinite tension in the corresponding portions of the membrane. To define the interface constrains that allow to apply such a state of prestress to the membrane, a physical apparatus is proposed, which can be numerically modelled with a two-dimensional mesh implementable in commercial finite-element model codes. This approach presents noteworthy advantages with respect to the three-dimensional modelling of the twisted bar.


Author(s):  
Francis H. Ku ◽  
Trevor G. Hicks ◽  
William R. Mabe ◽  
Jason R. Miller

Two-dimensional (2D) and three-dimensional (3D) weld-induced residual stress finite element analyses have been performed for 2-inch Schedule 80 Type-304 stainless steel pipe sections joined by a multi-layer segmented-bead pipe weld. The analyses investigate the similarities and differences between the two modeling approaches in terms of residual stresses and axial shrinkage induced by the pipe weld. The 2D analyses are of axisymmetric behavior and evaluate two different pipe end constraints, namely fixed-fixed and fixed-free, while the 3D analysis approximates the non-axisymmetric segmented welding expected in production, with fixed-free pipe end constraints. Based on the results presented, the following conclusions can be drawn. The welding temperature contour results between the 2D and 3D analyses are very similar. Only the 3D analysis is capable of simulating the non-axisymmetric behavior of the segmented welding technique. The 2D analyses yield similar hoop residual stresses to the 3D analysis, and closely capture the maximum and minimum ID surface hoop residual stresses from the 3D analysis. The primary difference in ID surface residual stresses between the 2D fixed-fixed and 2D fixed-free constraints cases is the higher tensile axial stresses in the pipe outside of the weld region. The 2D analyses under-predict the maximum axial residual stress compared to the 3D analysis. The 2D ID surface residual stress results tend to bound the averaged 3D results. 2D axisymmetric modeling tends to significantly under-predict weld shrinkage. Axial weld shrinkage from 3D modeling is of the same magnitude as values measured in the laboratory on a prototypic mockup.


1984 ◽  
Vol 142 ◽  
pp. 451-466 ◽  
Author(s):  
A. C. Robinson ◽  
P. G. Saffman

The three-dimensional linear stability of a rectilinear vortex of elliptical cross-section existing as a steady state in an irrotational straining field is studied numerically in the case of finite strain. It is shown that the instability predicted analytically for weak strain persists for finite strain and that the weak-strain results continue to be quantitatively valid for finite strain. The dependence of the growth rates of the unstable modes on the strain and the axial-disturbance wavelength is discussed. It is also shown that a three-dimensional instability is always more unstable than a two-dimensional instability in the range of parameters of most interest.


2004 ◽  
Vol 127 (4) ◽  
pp. 782-787 ◽  
Author(s):  
B. Zettl ◽  
W. Szyszkowski ◽  
W. J. Zhang

This paper discusses the finite element method (FEM) based modeling of the behavior of typical right circular flexure hinges used in planar compliant mechanisms. Such hinges have traditionally been approximated either by simple beams in the analytical approach or very often by two-dimensional (2D) plane stress elements when using the FEM. The three-dimensional (3D) analysis presented examines these approximations, focusing on systematic errors due to 2D modeling. It is shown that the 2D models provide only the lower (assuming the plane stress state) or the upper (assuming the plane strain state) limits of the hinge’s stiffness. The error of modeling a particular hinge by 2D elements (with either the plane stress or the plane strain assumptions) depends mainly on its depth-to-height ratio and may reach up to about 12%. However, this error becomes negligible for hinges with sufficiently high or sufficiently low depth-to-height ratios, in which either the plane strain or stress states dominate respectively. It is also shown that the computationally intensive 3D elements can be replaced, without sacrificing accuracy, by numerically efficient 2D elements if the material properties are appropriately manipulated.


1994 ◽  
Vol 361 ◽  
Author(s):  
V.A. Alyoshin ◽  
E.V. Sviridov ◽  
V.I.M. Hukhortov ◽  
I.H. Zakharchenko ◽  
V.P. Dudkevich

ABSTRACTSurface and cross-section relief evolution of ferroelectric epitaxial (Ba,Sr)TiO3 films rf-sputtered on (001) HgO crystal cle-avage surface versus the oxygen worKing gas pressure P and subst-rate temperature T were studied. Specific features of both three-dimensional and two-dimensional epitaxy mechanisms corresponding to various deposition conditions were revealed. Difference between low and high P-T-value 3D epitaxy was established. The deposition of films with mirror-smooth surfaces and perfect interfaces is shown to be possible.


1992 ◽  
Vol 241 ◽  
pp. 587-614 ◽  
Author(s):  
T. Dracos ◽  
M. Giger ◽  
G. H. Jirka

An experimental investigation of plane turbulent jets in bounded fluid layers is presented. The development of the jet is regular up to a distance from the orifice of approximately twice the depth of the fluid layer. From there on to a distance of about ten times the depth, the flow is dominated by secondary currents. The velocity distribution over a cross-section of the jet becomes three-dimensional and the jet undergoes a constriction in the midplane and a widening near the bounding surfaces. Beyond a distance of approximately ten times the depth of the bounded fluid layer the secondary currents disappear and the jet starts to meander around its centreplane. Large vortical structures develop with axes perpendicular to the bounding surfaces of the fluid layer. With increasing distance the size of these structures increases by pairing. These features of the jet are associated with the development of quasi two-dimensional turbulence. It is shown that the secondary currents and the meandering do not significantly affect the spreading of the jet. The quasi-two-dimensional turbulence, however, developing in the meandering jet, significantly influences the mixing of entrained fluid.


2021 ◽  
Vol 5 (4) ◽  
pp. 53-60
Author(s):  
Daniel Gurgul ◽  
Andriy Burbelko ◽  
Tomasz Wiktor

This paper presents a new proposition on how to derive mathematical formulas that describe an unknown Probability Density Function (PDF3) of the spherical radii (r3) of particles randomly placed in non-transparent materials. We have presented two attempts here, both of which are based on data collected from a random planar cross-section passed through space containing three-dimensional nodules. The first attempt uses a Probability Density Function (PDF2) the form of which is experimentally obtained on the basis of a set containing two-dimensional radii (r2). These radii are produced by an intersection of the space by a random plane. In turn, the second solution also uses an experimentally obtained Probability Density Function (PDF1). But the form of PDF1 has been created on the basis of a set containing chord lengths collected from a cross-section.The most important finding presented in this paper is the conclusion that if the PDF1 has proportional scopes, the PDF3 must have a constant value in these scopes. This fact allows stating that there are no nodules in the sample space that have particular radii belonging to the proportional ranges the PDF1.


Author(s):  
S. Khajehpour ◽  
R. G. Sauve´ ◽  
N. Badie

A method has been developed to incorporate the local three-dimensional shell behavior of two concentric tubes in the two-dimensional beam modeling of the problem. The two dimensional modeling of fuel channels in CANDU pressurized heavy water nuclear reactors is used in lieu of a more accurate three dimensional finite element approach in order to reduce the on-line simulation time which greatly affects the SLAR (Spacer Location And Repositioning) maintenance operation cost during outage. However, effort must be made to include the three-dimensional shell behavior of these channels into the two-dimensional modeling. In recent studies a nonlinear force-dependent model for contact stiffness between the calandria tube and pressure tube has been developed. However, local deformation of calandria the tube at spacer locations due to in-reactor creep leads to settling of the spacer into the calandria tube that consequently reduces the gap between the two tubes. In this work, the effect of local deformation (elastic and creep) of calandria tubes on modeling of contact at spacer locations is assessed using a three dimensional finite element code. The result is incorporated into a two-dimensional beam model of the problem as a reduction in size of the spacers that separate the two tubes. It is shown that the proposed method increases the accuracy of prediction of contact time and the spacer. In general, the method described in this paper suggests a way to incorporate local shell deformation into beam models of slender shell structure.


1944 ◽  
Vol 11 (3) ◽  
pp. A149-A161
Author(s):  
Gabriel Kron

Abstract This paper presents equivalent circuits representing the partial differential equations of the theory of elasticity for bodies of arbitrary shapes. Transient, steady-state, or sinusoidally oscillating elastic-field phenomena may now be studied, within any desired degree of accuracy, either by a “network analyzer,” or by numerical- and analytical-circuit methods. Such problems are the propagation of elastic waves, determination of the natural frequencies of vibration of elastic bodies, or of stresses and strains in steady-stressed states. The elastic body may be non-homogeneous, may have arbitrary shape and arbitrary boundary conditions, it may rotate at a uniform angular velocity and may, for representation, be divided into blocks of uneven length in different directions. The circuits are developed to handle both two- and three-dimensional phenomena. They are expressed in all types of orthogonal curvilinear reference frames in order to simplify the boundary relations and to allow the solution of three-dimensional problems with axial and other symmetry by the use of only a two-dimensional network. Detailed circuits are given for the important cases of axial symmetry, cylindrical co-ordinates (two-dimensional) and rectangular co-ordinates (two- and three-dimensional). Nonlinear stress-strain relations in the plastic range may be handled by a step-by-step variation of the circuit constants. Nonisotropic bodies and nonorthogonal reference frames, however, require an extension of the circuits given. The circuits for steady-state stress and small oscillation phenomena require only inductances and capacitors, while the circuits for transients require also standard (not ideal) transformers. A companion paper deals in detail with numerical and experimental methods to solve the equivalent circuits.


Sign in / Sign up

Export Citation Format

Share Document