Axisymmetric stress analysis
|
|
- Leo Lang
- 5 years ago
- Views:
Transcription
1 5 Axisymmetric stress analysis 5.1 Introduction The problem of stress distribution in bodies of revolution (axisymmetric solids) under axisymmetric loading is of considerable practical interest. The mathematical problems presented are very similar to those of plane stress and plane strain as, once again, the situation is two dimensional. 1,2 By symmetry, the two components of displacements in any plane section of the body along its axis of symmetry define completely the state of strain and, therefore, the state of stress. Such a cross-section is shown in Fig If r and z denote respectively the radial and axial coordinates of a point, with u and u being the corresponding displacements, it can readily be seen that precisely the same displacement functions as those used in Chapter 4 can be used to define the displacements within the triangular element i, j, m shown. The volume of material associated with an element is now that of a body of revolution indicated in Fig. 5.1, and all integrations have to be referred to this. The triangular element is again used mainly for illustrative purposes, the principles developed being completely general. In plane stress or strain problems it was shown that internal work was associated with three strain components in the coordinate plane, the stress component normal to this plane not being involved due to zero values of either the stress or the strain. In the axisymmetrical situation any radial displacement automatically induces a strain in the circumferential direction, and as the stresses in this direction are certainly non-zero, this fourth component of strain and of the associated stress has to be considered. Here lies the essential difference in the treatment of the axisymmetric situation. The reader will find the algebra involved in this chapter somewhat more tedious than that in the previous one but, essentially, identical operations are once again involved, following the general formulation of Chapter Element characteristics Displacement function Using the triangular shape of element (Fig. 5.1) with the nodes i, j, m numbered in the
2 Element characteristics 1 13 Fig. 5.1 Element of an axisymmetric solid. anticlockwise sense, we define the nodal displacement by its two components as ai={;}, and the element displacements by the vector { ai ae = aj am (5.1) Obviously, as in Sec , a linear polynomial can be used to define uniquely the displacements within the element. As the algebra involved is identical to that of Chapter 4 it will not be repeated here. The displacement field is now given again by Eq. (4.7): u = { :} = [IN,, INj, INm]ae (5.3) with ai + bir + c,z N. =, etc. 2A and I a two-by-two identity matrix. In the above ai = rjzm - r z. mj bi = zj - Z, ci = r, - rj etc., in cyclic order. Once again A is the area of the element triangle. (5.4)
3 114 Axisyrnmetric stress analysis Strain (total) As already mentioned, four components of strain have now to be considered. These are, in fact, all the non-zero strain components possible in an axisymmetric deformation. Figure 5.2 illustrates and defines these strains and the associated stresses. The strain vector defined below lists the strain components involved and defines them in terms of the displacements of a point. The expressions involved are almost self-evident and will not be derived here. The interested reader can consult a standard elasticity textbook3 for the full derivation. We thus have = su Using the displacement functions defined by Eqs (5.3) and (5.4) we have E = Ba' = [Bi, B,, Bm]ae in which B. = Fig. 5.2 Strains and stresses involved in the analysis of axisymmetric solids.
4 Element characteristics 1 15 With the B matrix now involving the coordinates r and z, the strains are no longer constant within an element as in the plane stress or strain case. This strain variation is due to the EO term. If the imposed nodal displacements are such that u is proportional to r then indeed the strains will all be constant. In addition, constant E, and -yrr strains may be deduced from a linear displacement. This is the only state of displacement coincident with a constant strain condition and it is clear that the displacement function satisfies the basic criterion of Chapter Initial strain (thermal strain) In general, four independent components of the initial strain vector can be envisaged: Eo = { Trio (5.7) Although this can, in general, be variable within the element, it will be convenient to take the initial strain as constant there. The most frequently encountered case of initial strain will be that due to thermal expansion. For an isotropic material we shall have then where 8 is the average temperature rise in an element and a is the coefficient of thermal expansion. The general case of anisotropy need not be considered since axial symmetry would be impossible to achieve under such circumstances. A case of some interest in practice is that of a stratified material, similar to the one discussed in Chapter 4, in which the plane of isotropy is normal to the axis of symmetry (Fig. 5.3). Here, two different expansion coefficients are possible: one in the axial direction a, and another in the plane normal to it, ar. Now the initial thermal strain becomes Practical cases of such stratified anisotropy often arise in laminated or fibreglass construction of machine components.
5 116 Axisymmetric stress analysis Fig. 5.3 Axisymmetrically stratified material Elasticity matrix The elasticity matrix D linking the strains E and the stresses (r in the standard form [Eq. (2.511, (r = { ".) = D(E - E ~ + ) (ro Trz needs now to be derived. The anisotropic 'stratified' material will be considered first, as the isotropic case can be simply presented as a special form. Anisotropic, stratified, material (Fig. 5.3) With the z-axis representing the normal to the planes of stratification we can rewrite Eqs (4.19) (again ignoring the initial strains and stresses for convenience) as & =----- Or "202 v100 &= v2ar az u20r r El E2 El E2 E2 E2 E#=---- vlur " Ti-, Ti-z = - El E2 El G2 (5.10)
6 Element characteristics 1 17 Writing again we have on solving for the stresses, that I 2 n(1 - nvq), ny(1 + VI), n(vl +nv$), VI, D,!!? ny(1 +v), d n(1 -nv;,, 0 sym. Isotropic material For an isotropic material we can obtain the D matrix by taking and El=E2=E or n=l VI = v2 = v and using the well-known relationship between isotropic elastic constants Substituting in Eq. (5.1 1) we now have _ G2 -_ =m=- 1 E2 E 2( 1 + v) 1 - v, v, v, E v, 1 - v, v, D= (l+v)(l-2v)[ 6: v, 1-v, 0, 0 0, (1-2v)/2 (5.11) (5.12) The stiffness matrix The stiffness matrix of the element ijm can now be computed according to the general relationship (2.13). Remembering that the volume integral has to be taken over the whole ring of material we have f KG = 27r BTDBjr dr dz (5.13) with B given by Eq. (5.6) and D by either Eq. (5.11) or Eq. (5.12), depending on the material. The integration cannot now be performed as simply as was the case in the plane stress problem because the B matrix depends on the coordinates. Two possibilities exist: the first is that of numerical integration and the second of an explicit multiplication and term-by-term integration. The simplest numerical integration procedure is to evaluate all quantities for a centroidal point r= ri +rj +r, zi + zj + z, 3 and?= 3
7 118 Axisymmetric stress analysis In this case we have simply as a first approximation K; = 27rBTDB,TA (5.14) with A being the triangle area and B the value of the strain-displacement matrix at the centroidal point. More elaborate numerical integration schemes could be used by evaluating the integrand at several points of the triangle. Such methods will be discussed in detail in Chapter 9. However, it can be shown that if the numerical integration is of such an order that the volume of the element is exactly determined by it, then in the limit of subdivision, the solution will converge to the exact an~wer.~ The one point integration suggested here is of this type, as it is well known that the volume of a body of revolution is given exactly by the product of the area and the path swept around by its centroid. With the simple triangular element used here a fairly fine subdivision is in any case needed for accuracy and most practical programs use the simple approximation which, surprisingly perhaps, is in fact usually superior to exact integration (see Chapter lo). One reason for this is the occurrence of logarithmic terms in the exact formulation. These involve ratios of the type ri/rm and, when the element is at a large distance from the axis, such terms tend to unity and evaluation of the logarithm is inaccurate External nodal forces In the case of the two-dimensional problems of the previous chapter the question of assigning of the external loads was so obvious as not to need further comment. In the present case, however, it is important to realize that the nodal forces represent a combined effect of the force acting along the whole circumference of the circle forming the element node. This point was already brought out in the integration of the expressions for the stiffness of an element, such integrations being conducted over the whole ring. Thus, if R represents the radial component of force per unit length of the circumference of a node at a radius r, the external force which will have to be introduced in the computation is 2rrR In the axial direction we shall, similarly, have 2rrZ to represent the combined effect of axial forces Nodal forces due to initial strain Again, by Eq. (2.13), fe = -27r BTDEor dr dz (5.15) J
8 Element characteristics 1 19 or noting that is constant, fp = -2II( [BTrdrdz)D&o (5.16) The integration should be performed in a similar manner to that used in the determination of the stiffness. It will readily be seen that, again, an approximate expression using a centroidal value is Initial stress forces are treated in an identical manner. fp = -27rBTD&oFA (5.17) Distributed body forces Distributed body forces, such as those due to gravity (if acting along the z-axis), centrifugal force in rotating machine parts, or pore pressure, often occur in axisymmetric problems. Let such forces be denoted by b= { i:} (5.18) per unit volume of material in the directions of r and z respectively. By the general equation (2.13) we have fp = -2n[IN,( ::}rdrdr (5.19) Using a coordinate shift similar to that of Sec it is easy to show that the first approximation, if the body forces are constant, results in (5.20) Although this is not exact the error term will be found to decrease with reduction of element size and, as it is also self-balancing, it will not introduce inaccuracies. Indeed, as will be shown in Chapter 10, the convergence rate is maintained. If the body forces are given by a potential similar to that defined in Sec , Le., (5.21) and if this potential is defined linearly by its nodal values, an expression equivalent to Eq. (4.42) can again be determined. In many problems the body forces vary proportionately to r. For example in rotating machinery we have centrifugal forces b, = w 2 pr where w is the angular velocity and p the density of the material. (5.22)
9 120 Axisymmetric stress analysis Fig. 5.4 Stresses in a sphere subject to an internal pressure (Poisson s ratio v = 0.3: (a) triangular mesh - centroidal values; (b) triangular mesh - nodal averages; (c) quadrilateral mesh obtained by averaging adjacent triangles.
10 Some illustrative examples Evaluation of stresses The stresses now vary throughout the element, as will be appreciated from Eqs (4.5) and (4.6). It is convenient now to evaluate the average stress at the centroid of the element. The stress matrix resulting from Eqs (5.6) and (2.3) gives there, as usual, 8 = DBae - DE^ + oo (5.23) It will be found that a certain amount of oscillation of stress values between elements occurs and better approximation can be achieved by averaging nodal stresses or recovery procedures of Chapter Some illustrative examples Test problems such as those of a cylinder under constant axial or radial stress give, as indeed would be expected, solutions which correspond to exact ones. This is again an obvious corollary of the ability of the displacement function to reproduce constant strain conditions. A problem for which an exact solution is available and in which almost linear stress gradients occur is that of a sphere subject to internal pressure. Figure 5.4(a) shows the centroidal stresses obtained using rather a coarse mesh, and the stress oscillation around the exact values should be noted. (This oscillation becomes even more pronounced at larger values of Poisson s ratio although the exact solution is independent of it.) In Fig. 5.4(b) the very much better approximation obtained by averaging the stresses at nodal points is shown, and in Fig. 5.4(c) a further improvement is given by element averaging. The close agreement with the exact solution even for the very coarse subdivision used here shows the accuracy achievable. The displacements at nodes compared with the exact solution are given in Fig Fig. 5.5 Displacements of internal and external surfaces of sphere under loading of Fig. 5.4.
11 122 Axisymmetric stress analysis Fig. 5.6 Sphere subject to steady-state heat flow (100 "C internal temperature, 0 C external temperature): (a) temperature and stress variation on radial section; (b) 'quadrilateral' averages. Fig. 5.7 A reactor pressure vessel. (a) 'Quadrilateral' mesh used in analysis; this was generated automatically by a computer. (b) Stresses due to a uniform internal pressure (automatic computer plot). Solution based on quadrilateral averages. (Poisson's ratio v = 0.1 5).
12 Early practical applications 123 In Fig. 5.6 thermal stresses in the same sphere are computed for the steady-state temperature variation shown. Again, excellent accuracy is demonstrated by comparison with the exact solution. 5.4 Early practical applications Two examples of practical applications of the programs available for axisymmetrical stress distribution are given here A prestressed concrete reactor pressure vessel Figure 5.7 shows the stress distribution in a relatively simple prototype pressure vessel. Due to symmetry only one-half of the vessel is analysed, the results given here referring to the components of stress due to internal pressure. In Fig. 5.8 contours of equal major principal stresses caused by temperature are shown. The thermal state is due to steady-state heat conduction and was itself found by the finite element method in a way described in Chapter 7. Fig. 5.8 A reactor pressure vessel. Thermal stresses due to steady-state heat conduction. Contours of major principal stress in pounds per square inch. (interior temperature 400 "C, exterior temperature 0 "C, a=5x 10-6/"C.E=2.58x 1061b/in2, v=o.15). -
13 124 Axisymmetric stress analysis Foundation pile Figure 5.9 shows the stress distribution around a foundation pile penetrating two different strata. This non-homogeneous problem presents no difficulties and is treated by the standard approach given in this chapter in which the quadrilateral elements shown are assemblies of two triangles and the results are averaged. 5.5 Non-symmetrical loading The method described in the present chapter can be extended to deal with nonsymmetrical loading. If the circumferential loading variation is expressed in circular harmonics then it is still possible to focus attention on one axial section although the nodal degrees of freedom are now increased to three. Details of this process are described in references 5 and 6 and in Chapter 9 of Volume Axisymmetry - plane strain and plane stress In the previous chapter we noted that plane stress and strain analysis was done in terms of three stress and strain components and, indeed, both cases could be generally incorporated in a single program with an indicator changing appropriate constants in the matrix D. Doing this loses track of the a, component in the plane strain case which has to be separately evaluated. Further, special expressions [viz. Eq had to be used to introduce initial strains. This is inconvenient (especially when non-linear constitutive laws are used), and an alternative of writing the plane strain case in terms of four stress-strain components as a special case of axisymmetric analysis is highly recommended. If the axisymmetric strain definition of Eq. (5.5) is examined, we note that r = 03 gives E~ 0 and plane strain conditions are obtained. Thus, if we ignore the terms in B associated with E ~ replace, the coordinates randz by xandy and further change the volume of integration 27rr to 1 the plane strain formulation becomes available from the axisymmetric plane strain directly. Plane stress conditions can similarly be incorporated, requiring in addition substitution of the axisymmetric D matrix by Eqs (4.13) or (4.19) augmented by an appropriate zero row and column. Thus, at the cost of additional storage of the fourth stress and strain component, all the cases discussed can be incorporated in a single format.
14 Axisymmetry - plane strain and plane stress 125 Fig. 5.9 (a) A pile in stratified soil. Irregular mesh and data for the problem. (b) A pile in stratified soil. Plot of vertical stresses on horizontal sections. Solution also plotted for Boussinesq problem obtained by making E, = E2 = EP,l,, and this is compared with exact values.
15 126 Axisymmetric stress analysis References 1. R.W. Clough. Chapter 7 of Stress Analysis (eds O.C. Zienkiewicz and G.S. Holister), Wiley, R.W. Clough and Y.R. Rashid. Finite element analysis of axi-symmetric solids. Proc. ASCE, 91, EM.l, 71, S. Timoshenko and J.N. Goodier. Theory of Elasticity. 2nd ed., McGraw-Hill, B.M. Irons. Comment on Stiffness matrices for section element by I.R. Raju and A.K. Rao. JAZAA, 7, 156-7, E.L. Wilson. Structural analysis of axisymmetric solids. JAZAA, 3, , O.C. Zienkiewicz. The Finite Element Method. 3rd ed., McGraw-Hill, 1977.
CRITERIA FOR SELECTION OF FEM MODELS.
CRITERIA FOR SELECTION OF FEM MODELS. Prof. P. C.Vasani,Applied Mechanics Department, L. D. College of Engineering,Ahmedabad- 380015 Ph.(079) 7486320 [R] E-mail:pcv-im@eth.net 1. Criteria for Convergence.
More informationApplications in Solid Mechanics
Companies, 4 CHAPTER 9 Applications in Solid 9. INTRODUCTION The bar and beam elements discussed in Chapters 4 are line elements, as only a single coordinate axis is required to define the element reference
More informationInterpolation Functions for General Element Formulation
CHPTER 6 Interpolation Functions 6.1 INTRODUCTION The structural elements introduced in the previous chapters were formulated on the basis of known principles from elementary strength of materials theory.
More informationENGN 2340 Final Project Report. Optimization of Mechanical Isotropy of Soft Network Material
ENGN 2340 Final Project Report Optimization of Mechanical Isotropy of Soft Network Material Enrui Zhang 12/15/2017 1. Introduction of the Problem This project deals with the stress-strain response of a
More informationNumerical integration in the axisymmetric finite element formulation
Advances in Engineering Software 31 (2000) 137 141 Short Communication Numerical integration in the axisymmetric finite element formulation J.D. Clayton a, J.J. Rencis b, * a The George W. Woodruff School
More informationROTATING RING. Volume of small element = Rdθbt if weight density of ring = ρ weight of small element = ρrbtdθ. Figure 1 Rotating ring
ROTATIONAL STRESSES INTRODUCTION High centrifugal forces are developed in machine components rotating at a high angular speed of the order of 100 to 500 revolutions per second (rps). High centrifugal force
More informationARTICLE A-8000 STRESSES IN PERFORATED FLAT PLATES
ARTICLE A-8000 STRESSES IN PERFORATED FLAT PLATES Delete endnote 18, which says "Express metric values in exponential form" A-8100 INTRODUCTION A-8110 SCOPE (a) This Article contains a method of analysis
More informationStresses Analysis of Petroleum Pipe Finite Element under Internal Pressure
ISSN : 48-96, Vol. 6, Issue 8, ( Part -4 August 06, pp.3-38 RESEARCH ARTICLE Stresses Analysis of Petroleum Pipe Finite Element under Internal Pressure Dr.Ragbe.M.Abdusslam Eng. Khaled.S.Bagar ABSTRACT
More informationChapter 3. Load and Stress Analysis. Lecture Slides
Lecture Slides Chapter 3 Load and Stress Analysis 2015 by McGraw Hill Education. This is proprietary material solely for authorized instructor use. Not authorized for sale or distribution in any manner.
More informationMathematics FINITE ELEMENT ANALYSIS AS COMPUTATION. What the textbooks don't teach you about finite element analysis. Chapter 3
Mathematics FINITE ELEMENT ANALYSIS AS COMPUTATION What the textbooks don't teach you about finite element analysis Chapter 3 Completeness and continuity: How to choose shape functions? Gangan Prathap
More informationFig. 1. Circular fiber and interphase between the fiber and the matrix.
Finite element unit cell model based on ABAQUS for fiber reinforced composites Tian Tang Composites Manufacturing & Simulation Center, Purdue University West Lafayette, IN 47906 1. Problem Statement In
More informationAnalytical Strip Method for Thin Isotropic Cylindrical Shells
IOSR Journal of Mechanical and Civil Engineering (IOSR-JMCE) e-issn: 2278-1684,p-ISSN: 2320-334X, Volume 14, Issue 4 Ver. III (Jul. Aug. 2017), PP 24-38 www.iosrjournals.org Analytical Strip Method for
More informationCIVL 8/7117 Chapter 12 - Structural Dynamics 1/75. To discuss the dynamics of a single-degree-of freedom springmass
CIV 8/77 Chapter - /75 Introduction To discuss the dynamics of a single-degree-of freedom springmass system. To derive the finite element equations for the time-dependent stress analysis of the one-dimensional
More informationTheoretical Manual Theoretical background to the Strand7 finite element analysis system
Theoretical Manual Theoretical background to the Strand7 finite element analysis system Edition 1 January 2005 Strand7 Release 2.3 2004-2005 Strand7 Pty Limited All rights reserved Contents Preface Chapter
More informationAdvanced Structural Analysis EGF Cylinders Under Pressure
Advanced Structural Analysis EGF316 4. Cylinders Under Pressure 4.1 Introduction When a cylinder is subjected to pressure, three mutually perpendicular principal stresses will be set up within the walls
More informationDAMPING OF GENERALIZED THERMO ELASTIC WAVES IN A HOMOGENEOUS ISOTROPIC PLATE
Materials Physics and Mechanics 4 () 64-73 Received: April 9 DAMPING OF GENERALIZED THERMO ELASTIC WAVES IN A HOMOGENEOUS ISOTROPIC PLATE R. Selvamani * P. Ponnusamy Department of Mathematics Karunya University
More informationExample-3. Title. Description. Cylindrical Hole in an Infinite Mohr-Coulomb Medium
Example-3 Title Cylindrical Hole in an Infinite Mohr-Coulomb Medium Description The problem concerns the determination of stresses and displacements for the case of a cylindrical hole in an infinite elasto-plastic
More informationApplication of pseudo-symmetric technique in dynamic analysis of concrete gravity dams
Application of pseudo-symmetric technique in dynamic analysis of concrete gravity dams V. Lotfi Department of Civil and Environmental Engineering, Amirkabir University, Iran Abstract A new approach is
More informationStress analysis of a stepped bar
Stress analysis of a stepped bar Problem Find the stresses induced in the axially loaded stepped bar shown in Figure. The bar has cross-sectional areas of A ) and A ) over the lengths l ) and l ), respectively.
More informationContents as of 12/8/2017. Preface. 1. Overview...1
Contents as of 12/8/2017 Preface 1. Overview...1 1.1 Introduction...1 1.2 Finite element data...1 1.3 Matrix notation...3 1.4 Matrix partitions...8 1.5 Special finite element matrix notations...9 1.6 Finite
More informationME FINITE ELEMENT ANALYSIS FORMULAS
ME 2353 - FINITE ELEMENT ANALYSIS FORMULAS UNIT I FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PROBLEMS 01. Global Equation for Force Vector, {F} = [K] {u} {F} = Global Force Vector [K] = Global Stiffness
More informationBending of Simply Supported Isotropic and Composite Laminate Plates
Bending of Simply Supported Isotropic and Composite Laminate Plates Ernesto Gutierrez-Miravete 1 Isotropic Plates Consider simply a supported rectangular plate of isotropic material (length a, width b,
More informationTIME-DEPENDENT BEHAVIOR OF PILE UNDER LATERAL LOAD USING THE BOUNDING SURFACE MODEL
TIME-DEPENDENT BEHAVIOR OF PILE UNDER LATERAL LOAD USING THE BOUNDING SURFACE MODEL Qassun S. Mohammed Shafiqu and Maarib M. Ahmed Al-Sammaraey Department of Civil Engineering, Nahrain University, Iraq
More informationInfluence of the filament winding process variables on the mechanical behavior of a composite pressure vessel
Influence of the filament winding process variables on the mechanical behavior of a composite pressure vessel G. Vargas 1 & A. Miravete 2 1 Universidad Pontificia Bolivariana, Facultad de Ingeniería Mecánica,
More informationAxisymmetric Stress Analysis of Internally Pressurized Rotating Cylinder using Finite Element Method
Axisymmetric Stress Analysis of Internally Pressurized Rotating Cylinder using Finite Element Method A Project Report Submitted in partial fulfillment for the award of the degree Of BACHELOR OF TECHNOLOGY
More informationPlane and axisymmetric models in Mentat & MARC. Tutorial with some Background
Plane and axisymmetric models in Mentat & MARC Tutorial with some Background Eindhoven University of Technology Department of Mechanical Engineering Piet J.G. Schreurs Lambèrt C.A. van Breemen March 6,
More informationUsing MATLAB and. Abaqus. Finite Element Analysis. Introduction to. Amar Khennane. Taylor & Francis Croup. Taylor & Francis Croup,
Introduction to Finite Element Analysis Using MATLAB and Abaqus Amar Khennane Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Croup, an informa business
More informationHeat Transfer Analysis of Machine Tool Main Spindle
Technical Paper Heat Transfer Analysis of Machine Tool Main Spindle oshimitsu HIRASAWA Yukimitsu YAMAMOTO CAE analysis is very useful for shortening development time and reducing the need for development
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationEFFECTS OF THERMAL STRESSES AND BOUNDARY CONDITIONS ON THE RESPONSE OF A RECTANGULAR ELASTIC BODY MADE OF FGM
Proceedings of the International Conference on Mechanical Engineering 2007 (ICME2007) 29-31 December 2007, Dhaka, Bangladesh ICME2007-AM-76 EFFECTS OF THERMAL STRESSES AND BOUNDARY CONDITIONS ON THE RESPONSE
More informationLevel 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method
9210-203 Level 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method You should have the following for this examination one answer book No additional data is attached
More informationLecture 8. Stress Strain in Multi-dimension
Lecture 8. Stress Strain in Multi-dimension Module. General Field Equations General Field Equations [] Equilibrium Equations in Elastic bodies xx x y z yx zx f x 0, etc [2] Kinematics xx u x x,etc. [3]
More informationAppendix A: Details of Deformation Analyses
Appendix A: Details of Deformation Analyses This appendix presents details of theory and procedures for some aspects of the analyses presented in Chapter 4. Specifically this appendix consists of the following:
More informationDESIGN OF LAMINATES FOR IN-PLANE LOADING
DESIGN OF LAMINATES FOR IN-PLANOADING G. VERCHERY ISMANS 44 avenue F.A. Bartholdi, 72000 Le Mans, France Georges.Verchery@m4x.org SUMMARY This work relates to the design of laminated structures primarily
More informationParametric Study Of The Material On Mechanical Behavior Of Pressure Vessel
Parametric Study Of The Material On Mechanical Behavior Of Pressure Vessel Dr. Mohammad Tariq Assistant Professor Mechanical Engineering Department SSET, SHIATS-Deemed University, Naini, Allahabad, U.P.,
More informationAdvantages, Limitations and Error Estimation of Mixed Solid Axisymmetric Modeling
Advantages, Limitations and Error Estimation of Mixed Solid Axisymmetric Modeling Sudeep Bosu TATA Consultancy Services Baskaran Sundaram TATA Consultancy Services, 185 Lloyds Road,Chennai-600086,INDIA
More informationCHAPTER 8: Thermal Analysis
CHAPER 8: hermal Analysis hermal Analysis: calculation of temperatures in a solid body. Magnitude and direction of heat flow can also be calculated from temperature gradients in the body. Modes of heat
More informationEffect of fibre shape on transverse thermal conductivity of unidirectional composites
Sādhanā Vol. 4, Part 2, April 25, pp. 53 53. c Indian Academy of Sciences Effect of fibre shape on transverse thermal conductivity of unidirectional composites B RAGHAVA RAO,, V RAMACHANDRA RAJU 2 and
More informationThermo Mechanical Analysis of AV1 Diesel Engine Piston using FEM
Journal of Advanced Engineering Research ISSN: 2393-8447 Volume 2, Issue 1, 2015, pp.23-28 Thermo Mechanical Analysis of AV1 Diesel Engine Piston using FEM Subodh Kumar Sharma 1, *, P. K. Saini 2, N. K.
More informationPRESSURE VESSELS & PRESSURE CABINS FOR BLENDED WING BODIES
PRESSURE VESSELS & PRESSURE CABINS FOR BLENDED WING BODIES F.J.J.M.M. Geuskens, O.K. Bergsma 2, S. Koussios 2 & A. Beukers 3 PhD Researcher, 2 Associate professor, 3 Professor / DPCS, TU Delft Kluyverweg,
More informationON THE PREDICTION OF EXPERIMENTAL RESULTS FROM TWO PILE TESTS UNDER FORCED VIBRATIONS
Transactions, SMiRT-24 ON THE PREDICTION OF EXPERIMENTAL RESULTS FROM TWO PILE TESTS UNDER FORCED VIBRATIONS 1 Principal Engineer, MTR & Associates, USA INTRODUCTION Mansour Tabatabaie 1 Dynamic response
More informationPlane Strain Test for Metal Sheet Characterization
Plane Strain Test for Metal Sheet Characterization Paulo Flores 1, Felix Bonnet 2 and Anne-Marie Habraken 3 1 DIM, University of Concepción, Edmundo Larenas 270, Concepción, Chile 2 ENS - Cachan, Avenue
More informationMAE 323: Chapter 6. Structural Models
Common element types for structural analyis: oplane stress/strain, Axisymmetric obeam, truss,spring oplate/shell elements o3d solid ospecial: Usually used for contact or other constraints What you need
More informationExercise: concepts from chapter 8
Reading: Fundamentals of Structural Geology, Ch 8 1) The following exercises explore elementary concepts associated with a linear elastic material that is isotropic and homogeneous with respect to elastic
More informationTuesday, February 11, Chapter 3. Load and Stress Analysis. Dr. Mohammad Suliman Abuhaiba, PE
1 Chapter 3 Load and Stress Analysis 2 Chapter Outline Equilibrium & Free-Body Diagrams Shear Force and Bending Moments in Beams Singularity Functions Stress Cartesian Stress Components Mohr s Circle for
More informationPhysics 342 Lecture 23. Radial Separation. Lecture 23. Physics 342 Quantum Mechanics I
Physics 342 Lecture 23 Radial Separation Lecture 23 Physics 342 Quantum Mechanics I Friday, March 26th, 2010 We begin our spherical solutions with the simplest possible case zero potential. Aside from
More informationTwo Dimensional State of Stress and Strain: examples
Lecture 1-5: Two Dimensional State of Stress and Strain: examples Principal stress. Stresses on oblique plane: Till now we have dealt with either pure normal direct stress or pure shear stress. In many
More informationTruss Structures: The Direct Stiffness Method
. Truss Structures: The Companies, CHAPTER Truss Structures: The Direct Stiffness Method. INTRODUCTION The simple line elements discussed in Chapter introduced the concepts of nodes, nodal displacements,
More information4 The matrix approach
4 The matrix approach 4. Introduction It is not a coincidence that the finite element method blossomed at the time when computers were being developed, because it calls for a vast amount of computation
More informationA Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials
Dublin, October 2010 A Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials FracMan Technology Group Dr Mark Cottrell Presentation Outline Some Physical
More informationAnalyses of Soil-Structure Systems
Analyses of Soil-Structure Systems Julio García Paul C. Rizzo Associates, Inc., 105 Mall Blvd. Suite 270E, Monroeville, PA 15146, USA Phone: +1 (412)856-9700, Fax: +1 (412)856-9749, Email: julio.garcia@rizzoassoc.com
More informationSolution Methods. Steady State Diffusion Equation. Lecture 04
Solution Methods Steady State Diffusion Equation Lecture 04 1 Solution methods Focus on finite volume method. Background of finite volume method. Discretization example. General solution method. Convergence.
More informationAn Analytical Solution for Hoop Tension in Liquid Storage Cylindrical Tanks
International Journal of Engineering and Applied Sciences (IJEAS) ISSN: 2394-3661, Volume-4, Issue-7, July 2017 An Analytical Solution for Hoop Tension in Liquid Storage Cylindrical Tanks Anand Daftardar,
More informationMIXED RECTANGULAR FINITE ELEMENTS FOR PLATE BENDING
144 MIXED RECTANGULAR FINITE ELEMENTS FOR PLATE BENDING J. N. Reddy* and Chen-Shyh-Tsay School of Aerospace, Mechanical and Nuclear Engineering, University of Oklahoma, Norman, Oklahoma The paper describes
More informationThe Rotating Inhomogeneous Elastic Cylinders of. Variable-Thickness and Density
Applied Mathematics & Information Sciences 23 2008, 237 257 An International Journal c 2008 Dixie W Publishing Corporation, U. S. A. The Rotating Inhomogeneous Elastic Cylinders of Variable-Thickness and
More informationUniversity of Illinois at Urbana-Champaign College of Engineering
University of Illinois at Urbana-Champaign College of Engineering CEE 570 Finite Element Methods (in Solid and Structural Mechanics) Spring Semester 03 Quiz # April 8, 03 Name: SOUTION ID#: PS.: A the
More informationACCURATE FREE VIBRATION ANALYSIS OF POINT SUPPORTED MINDLIN PLATES BY THE SUPERPOSITION METHOD
Journal of Sound and Vibration (1999) 219(2), 265 277 Article No. jsvi.1998.1874, available online at http://www.idealibrary.com.on ACCURATE FREE VIBRATION ANALYSIS OF POINT SUPPORTED MINDLIN PLATES BY
More informationMITOCW MITRES2_002S10nonlinear_lec15_300k-mp4
MITOCW MITRES2_002S10nonlinear_lec15_300k-mp4 The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources
More informationAnalysis of asymmetric radial deformation in pipe with local wall thinning under internal pressure using strain energy method
Analysis of asymmetric radial deformation in pipe with local wall thinning under internal pressure using strain energy method V.M.F. Nascimento Departameto de ngenharia Mecânica TM, UFF, Rio de Janeiro
More information4. BEAMS: CURVED, COMPOSITE, UNSYMMETRICAL
4. BEMS: CURVED, COMPOSITE, UNSYMMETRICL Discussions of beams in bending are usually limited to beams with at least one longitudinal plane of symmetry with the load applied in the plane of symmetry or
More informationLecture 21: Isoparametric Formulation of Plane Elements.
6.6. Rectangular Plane Stress/Strain Element. he CS element uses a triangular shape. he 3 nodes of the CS element allow us to employ linear multivariate approximations to the u and v displacements. he
More informationStress and Displacement Analysis of a Rectangular Plate with Central Elliptical Hole
Stress and Displacement Analysis of a Rectangular Plate with Central Elliptical Hole Dheeraj Gunwant, J. P. Singh mailto.dheerajgunwant@gmail.com, jitenderpal2007@gmail.com, AIT, Rampur Abstract- A static
More informationTriangular Plate Displacement Elements
Triangular Plate Displacement Elements Chapter : TRIANGULAR PLATE DISPLACEMENT ELEMENTS TABLE OF CONTENTS Page. Introduction...................... Triangular Element Properties................ Triangle
More informationThe Finite Element Method for Solid and Structural Mechanics
The Finite Element Method for Solid and Structural Mechanics Sixth edition O.C. Zienkiewicz, CBE, FRS UNESCO Professor of Numerical Methods in Engineering International Centre for Numerical Methods in
More informationJEPPIAAR ENGINEERING COLLEGE
JEPPIAAR ENGINEERING COLLEGE Jeppiaar Nagar, Rajiv Gandhi Salai 600 119 DEPARTMENT OFMECHANICAL ENGINEERING QUESTION BANK VI SEMESTER ME6603 FINITE ELEMENT ANALYSIS Regulation 013 SUBJECT YEAR /SEM: III
More information7.6 Stress in symmetrical elastic beam transmitting both shear force and bending moment
7.6 Stress in symmetrical elastic beam transmitting both shear force and bending moment à It is more difficult to obtain an exact solution to this problem since the presence of the shear force means that
More informationFEA A Guide to Good Practice. What to expect when you re expecting FEA A guide to good practice
FEA A Guide to Good Practice What to expect when you re expecting FEA A guide to good practice 1. Background Finite Element Analysis (FEA) has transformed design procedures for engineers. Allowing more
More informationModule-6: Laminated Composites-II. Learning Unit-1: M6.1. M 6.1 Structural Mechanics of Laminates
Module-6: Laminated Composites-II Learning Unit-1: M6.1 M 6.1 Structural Mechanics of Laminates Classical Lamination Theory: Laminate Stiffness Matrix To this point in the development of classical lamination
More informationIn all of the following equations, is the coefficient of permeability in the x direction, and is the hydraulic head.
Groundwater Seepage 1 Groundwater Seepage Simplified Steady State Fluid Flow The finite element method can be used to model both steady state and transient groundwater flow, and it has been used to incorporate
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 13
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras (Refer Slide Time: 00:25) Module - 01 Lecture - 13 In the last class, we have seen how
More informationLecture 12: Finite Elements
Materials Science & Metallurgy Part III Course M6 Computation of Phase Diagrams H. K. D. H. Bhadeshia Lecture 2: Finite Elements In finite element analysis, functions of continuous quantities such as temperature
More informationPractice Final Examination. Please initial the statement below to show that you have read it
EN175: Advanced Mechanics of Solids Practice Final Examination School of Engineering Brown University NAME: General Instructions No collaboration of any kind is permitted on this examination. You may use
More informationComputational Stiffness Method
Computational Stiffness Method Hand calculations are central in the classical stiffness method. In that approach, the stiffness matrix is established column-by-column by setting the degrees of freedom
More informationHyperbolic Soil Bearing Capacity
1 Introduction Hyperbolic Soil Bearing Capacity This example involves analyzing the bearing capacity of a round footing. The example is useful for illustrating several SIGMA/W features and procedures,
More informationDHANALAKSHMI COLLEGE OF ENGINEERING, CHENNAI DEPARTMENT OF MECHANICAL ENGINEERING ME 6603 FINITE ELEMENT ANALYSIS PART A (2 MARKS)
DHANALAKSHMI COLLEGE OF ENGINEERING, CHENNAI DEPARTMENT OF MECHANICAL ENGINEERING ME 6603 FINITE ELEMENT ANALYSIS UNIT I : FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PART A (2 MARKS) 1. Write the types
More informationEFFECT OF ELLIPTIC OR CIRCULAR HOLES ON THE STRESS DISTRIBUTION IN PLATES
EFFECT OF ELLIPTIC OR CIRCULAR HOLES ON THE STRESS DISTRIBUTION IN PLATES OF WOOD OR PLYWOOD CONSIDERED AS ORTHOTROPIC MATERIALS Information Revied and Reaffirmed March 1956 No. 1510 EFFECT OF ELLIPTIC
More informationA new closed-form model for isotropic elastic sphere including new solutions for the free vibrations problem
A new closed-form model for isotropic elastic sphere including new solutions for the free vibrations problem E Hanukah Faculty of Mechanical Engineering, Technion Israel Institute of Technology, Haifa
More informationEffect of embedment depth and stress anisotropy on expansion and contraction of cylindrical cavities
Effect of embedment depth and stress anisotropy on expansion and contraction of cylindrical cavities Hany El Naggar, Ph.D., P. Eng. and M. Hesham El Naggar, Ph.D., P. Eng. Department of Civil Engineering
More informationGeneral method for simulating laboratory tests with constitutive models for geomechanics
General method for simulating laboratory tests with constitutive models for geomechanics Tomáš Janda 1 and David Mašín 2 1 Czech Technical University in Prague, Faculty of Civil Engineering, Czech Republic
More informationThermoelastic Stresses in a Rod Subjected to Periodic Boundary Condition: An Analytical Treatment
Thermoelastic Stresses in a Rod Subjected to Periodic Boundary Condition: An Analytical Treatment Ahmet N. Eraslan Department of Engineering Sciences Middle East Technical University Ankara 06531, Turkey
More informationFrequently Asked Questions
Frequently Asked Questions Why do we have to make the assumption that plane sections plane? How about bars with non-axis symmetric cross section? The formulae derived look very similar to beam and axial
More informationSTRESSES WITHIN CURVED LAMINATED BEAMS OF DOUGLAS-FIR
UNITED STATES DEPARTMENT OF AGRICULTURE. FOREST SERVICE - FOREST PRODUCTS LABORATORY - MADISON, WIS. STRESSES WITHIN CURVED LAMINATED BEAMS OF DOUGLAS-FIR NOVEMBER 1963 FPL-020 STRESSES WITHIN CURVED LAMINATED
More informationGeneralized Finite Element Methods for Three Dimensional Structural Mechanics Problems. C. A. Duarte. I. Babuška and J. T. Oden
Generalized Finite Element Methods for Three Dimensional Structural Mechanics Problems C. A. Duarte COMCO, Inc., 7800 Shoal Creek Blvd. Suite 290E Austin, Texas, 78757, USA I. Babuška and J. T. Oden TICAM,
More informationBenefits of Collaboration between Centrifuge Modeling and Numerical Modeling. Xiangwu Zeng Case Western Reserve University, Cleveland, Ohio
Benefits of Collaboration between Centrifuge Modeling and Numerical Modeling Xiangwu Zeng Case Western Reserve University, Cleveland, Ohio ABSTRACT There is little doubt that collaboration between centrifuge
More informationMECH 576 Geometry in Mechanics September 16, 2009 Using Line Geometry
MECH 576 Geometry in Mechanics September 16, 2009 Using Line Geometry 1 Deriving Equations in Line Coordinates Four exercises in deriving fundamental geometric equations with line coordinates will be conducted.
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 11
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Module - 01 Lecture - 11 Last class, what we did is, we looked at a method called superposition
More informationVIBRATION AND DAMPING ANALYSIS OF FIBER REINFORCED COMPOSITE MATERIAL CONICAL SHELLS
VIBRATION AND DAMPING ANALYSIS OF FIBER REINFORCED COMPOSITE MATERIAL CONICAL SHELLS Mechanical Engineering Department, Indian Institute of Technology, New Delhi 110 016, India (Received 22 January 1992,
More informationFinite Element Method-Part II Isoparametric FE Formulation and some numerical examples Lecture 29 Smart and Micro Systems
Finite Element Method-Part II Isoparametric FE Formulation and some numerical examples Lecture 29 Smart and Micro Systems Introduction Till now we dealt only with finite elements having straight edges.
More informationMethods of Analysis. Force or Flexibility Method
INTRODUCTION: The structural analysis is a mathematical process by which the response of a structure to specified loads is determined. This response is measured by determining the internal forces or stresses
More informationStiffness Matrices, Spring and Bar Elements
CHAPTER Stiffness Matrices, Spring and Bar Elements. INTRODUCTION The primary characteristics of a finite element are embodied in the element stiffness matrix. For a structural finite element, the stiffness
More informationLecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity
Lecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity by Borja Erice and Dirk Mohr ETH Zurich, Department of Mechanical and Process Engineering, Chair of Computational Modeling
More informationMATERIAL PROPERTIES. Material Properties Must Be Evaluated By Laboratory or Field Tests 1.1 INTRODUCTION 1.2 ANISOTROPIC MATERIALS
. MARIAL PROPRIS Material Properties Must Be valuated By Laboratory or Field ests. INRODUCION he fundamental equations of structural mechanics can be placed in three categories[]. First, the stress-strain
More informationCurvature of a Cantilever Beam Subjected to an Equi-Biaxial Bending Moment. P. Krulevitch G. C. Johnson
UCRL-JC-30440 PREPRINT Curvature of a Cantilever Beam Subjected to an Equi-Biaxial Bending Moment P. Krulevitch G. C. Johnson This paper was prepared for submittal to the Materials Research Society Spring
More informationDevelopment of Truss Equations
CIVL 7/87 Chapter 3 - Truss Equations - Part /53 Chapter 3a Development of Truss Equations Learning Objectives To derive the stiffness matri for a bar element. To illustrate how to solve a bar assemblage
More informationStretching of a Prismatic Bar by its Own Weight
1 APES documentation (revision date: 12.03.10) Stretching of a Prismatic Bar by its Own Weight. This sample analysis is provided in order to illustrate the correct specification of the gravitational acceleration
More informationPost Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method
9210-220 Post Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method You should have the following for this examination one answer book scientific calculator No
More informationChapter 5 Structural Elements: The truss & beam elements
Institute of Structural Engineering Page 1 Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 2 Chapter Goals Learn how to formulate the Finite Element Equations
More informationCHAPTER 5. Beam Theory
CHPTER 5. Beam Theory SangJoon Shin School of Mechanical and erospace Engineering Seoul National University ctive eroelasticity and Rotorcraft Lab. 5. The Euler-Bernoulli assumptions One of its dimensions
More informationDynamic Analysis Contents - 1
Dynamic Analysis Contents - 1 TABLE OF CONTENTS 1 DYNAMIC ANALYSIS 1.1 Overview... 1-1 1.2 Relation to Equivalent-Linear Methods... 1-2 1.2.1 Characteristics of the Equivalent-Linear Method... 1-2 1.2.2
More information1 Exercise: Linear, incompressible Stokes flow with FE
Figure 1: Pressure and velocity solution for a sinking, fluid slab impinging on viscosity contrast problem. 1 Exercise: Linear, incompressible Stokes flow with FE Reading Hughes (2000), sec. 4.2-4.4 Dabrowski
More information