Applied. Attalysis. FitBute Element. JLarry «/. Segerlind. JOHN Wn.EY & SONS New York Chichester Brisbane Toronto Singapore
|
|
- Mitchell Harrison
- 5 years ago
- Views:
Transcription
1 Applied FitBute Element Attalysis JLarry «/. Segerlind Associate Professor Department of Agricultural Engineering Michigan State University JOHN Wn.EY & SONS New York Chichester Brisbane Toronto Singapore
2 6 Tarsian af JVancircular Sections The preceding four chapters cover the discretizing of the body, the interpolating polynomial for a single element, the interpolation polynomials for a discretized body, and the derivation of the finite element equations. Each of these chapters contains basic information that is related to the finite element method. This chapter is a transition from the theory to the implementation. Its objec tive is to illustrate the steps involved in imple menting the finite element method. This objective is accomplished by obtaining a numerical solu tion for a problem that involves the torsion of a noncircular section. There are two reasons for selecting the torsion of noncircular sections to illustrate the im plementation of the finite element method. The element equations are relatively straightforward to derive. The [K] matrix is easy to calculate and the boundary integrals do not have to be evaluated because the boundary values are pre scribed. These properties prevent the evaluation of the element equations from detracting from our main objective. The other reason is that the torsion of a noncircular section involves concepts important to both mechanics-oriented and field- 87
3 88 Torsion of Nonclrcular Sections problem-oriented individuals. The theory will be familiar to an individual with a mechanics background while the governing differential equation is similar to those governing heat transfer and groundwater flow. 6.1 GENERAL THEORY, TORSION OF NONCIRCULAR SECTIONS One theory* for calculating the shear stresses in a noncircular shaft subjected to a twisting moment T about the z axis (Fig. 6.1) states that the shear stress components at any point can be calculated using 9<J) dy d(j> dx (6.1) where (p is a stress function. The governing differential equation is 1 G dx^ 1 9 G (6.2) with </) = (6.3) on the boundary. The physical parameters in (6.2) are the shear modulus of the Fig. 6.1 The shear stress components in a noncircular section subjected to a torque loading. *There are two theories for calculating the shear stresses that result when a noncircular shaft is subjected to torsion. One of them was developed by St. Venant; the other was proposed by Prantl. Both theories are discussed by Fung (1965). The minimization formulation of Prandtl's theory, which is used in this chapter, is also discussed by Temoshenko and Goodier (197).
4 General Theory, Torsion of Nonclrcular Sections 89 material, G, (N/cm^), and the angle of twist per unit length, 6 (rad/cm). This formulation does not have the applied torque, T (N m) in the governing equa tion. Instead, T is calculated once <p is known, using ^=2f <j>da ''Area (6.4) The stress function represents a surface covering the cross section of the shaft (Fig. 6.2). The twisting moment is proportional to the volume under this surface while the shear stresses are related to the gradients in the a and y coordinate directions. Fig. 6.2 The <f) surface and a related shear stress component. Equation 6.2 is usually written d^cl>^ = dx 9y (6.5) when the shaft is composed of a single material. The integral formulation of (6.5) is dv (6.6) which according to our discussion in Chapter 5 can be rewritten as = j\\{gv[d]{g]-(2ge)<^ dv (6.7)
5 9 Torsion of Nonclrcuiar Sections where 9</) dx d<j> dy and [/)] = K.. K The column vector (g) is related to the shear stress components in this application, while [D] becomes the identity matrix since K^^ = K^^ = 1. The minimization of x with respect to (O) results in the set of linear equations E E Sf S f {'^G9)dV (6.8) e = 1 I/(^) / l/(^) e=\ where is defined by (4.1) and the gradient matrix, is defined by (5.39). The evaluation of (6.8) commences as soon as we select a shape and divide it into elements. 6.2 ASSEMBLAGE OF THE ELEMENT MATRICES The square shaft (Fig. 6.3a) will be used to illustrate the evaluation and the assemblage of the element matrices into a set of linear equations. This shaft has four axes of symmetry, therefore, only one eighth of the total cross section needs to be analyzed. This fractional portion of the cross section is divided into four elements as shown in Fig. 6.3b. Four elements are not sufficient to obtain an accurate answer, but they are enough to illustrate the assemblage technique that is one of our objectives. Following the methods of Chapter 4, the element interpolating polynomials are <^( = + ^) C >(2) =, + 7V2^2^ yvi2^5 + ^ (6.9) cj>(3) = GO J -h A^2^^> A^i^^4 + A^5^^^5 + ^ =, A^i'^^4 4 A^i^^5 + yv6^^^6 The general equation for the element stiffness matrix is J y since it has already been noted that [/)] = [/] for the torsion of a noncircular
6 Assemblage of the Element Matrices 91 Axes of symmetry Shaft is 1 cm2 G =.8 X 1^ N/cm2 = 1 deg in ICQ cm (4).25 cm f * i) (3) +.25 cm X (1). (2) A V -.25 cm cm - Fig. 6.3 Element subdivi sion for the torsion of a square shaft. section. Evaluation of involves differentiating with respect to x and y. Confining our attention to element one dx dx am" 9x 1 2^( [M" o] (6.1a, b) d(j> (1) dy dy dy dy 2A (1)[cS" 4'^ o]
7 92 Torsion of Noncircuiar Sections The gradient matrix is [5<"] = (1).(1) ^2 cv^ (6.11) The area for this element is 32 and 2/4"> = 16 The b and c coefficients are b\'^=y2-y,= -.25 c\'^=x^-x2 = b['^=y,-y,=.25 c(" = A',-X4 = M'>= 7,-72= 4" = ^2-'^i = Substituting these values into (6.11) yields [5<"] = (6.12) and the product is IS [fi<'>]^[5">]= or (6.13) The element stiffness matrix is the integral of (6.13). Since the matrix product consists of constant values, it may be removed from under the
8 Assemblage off the Element Matrices 93 integral yielding (1) assuming a unit thickness. Recalling that (6.13) yields 1/32 and combining this with (6.14) The evaluation of the volume integral (/'')= [ 2G<"l9 ^ ^ JyO) yvi'* dv is straightforward if we employ the area coordinate system discussed in Chapter 3. Defining the area coordinates as L, = Af)'\ T2=A')" and = (6.15) the volume integral becomes ^ L, ^ ' JyO) dv (6.16) Assuming a unit thickness and using the area integral for area coordinates, (3.43), yields
9 94 Torsion of Nonclrcular Sections (/")- (6.17) Substitution of the values of, and yields* in (6.18) The system of equations for element one is = (/'>} or " $ <I' (6.19a ) <^5 ^>6 A similar set of equations can be calculated for each of the other elements using an identical procedure. The resulting matrices are given below. = (/«'} ' (6.19b ) -1 1 % 29.7 ^6 *The units are N/cm^ for G, cm^ for twists 1 deg in 1 cm. and rad/cm for. ^ = 7r/18x 1/1, since the shaft
10 Assemblage of the Element Matrices 95 [ A: ]{f } = {/"*} " ' 1-1 <^ % (6.19c ) a. ' 1 $ $ (6.19d ) -1 :! -1 % ^ The final system of equations is the algebraic sum of the element equations, which is 1-1 P $, b j C.O (6.2 ) The nodal valu^ $3, $5, and are zero,.sin h these nodes are on the external boundary. Modifibaftbtrxif (672) and solution yields <I), = <f>3 = -\ o - >1 "' i- ^>2=16 <I>5 = (6.21) $4= $6 = The modification of (6.2) is discussed in the next chapter where we consider the computer techniques that related to the implementation of the finite element method. The (j) surface for this set of nodal values is presented in Fig. 6.4.
11 96 Torsion of Nonclrcular Sections 4), = = Fig. 6.4 The nodal values for the four-element torsion problem. The determination of the nodal values is a major step in the solution of the problem. In most problems, however, a set of element resultants must be calculated. The shear stress values in each element and the twisting torque T which produces the angle of twist 6 are of interest in the present example. The calculation of the element resultants is discussed in the next section. 6.3 CONVENTIONAL ELEMENT RESULTANTS The gradients of the nodal parameter, (f>, are the important quantities in the torsion problem because the shear stresses are related to these gradients by d<f) dy and dxj) a;c The shear stress values are readily calculated because the gradient matrix for
12 Conventional Element Resultants 97 each element has already been evaluated. The gradient matrix for element one is (6.11). d(j> dx d<j) dy 2A^'^ 6^" c{" c^" 4" Combining (6.12) and (6.21) { = therefore, 4^)= = N/cm^ dy ' 4=-.^ dx N/cm^ The stress components for the other elements are Element 2 = ON/cm^. T = N/cm^ Element 3 = N/cm^, t = 494. N/cm^ Element 4 = ON/cm^. T- =494. N/cm^ These values are shown schematically in Fig. 6.5.
13 98 Torsion of Nonclrcular Sections Positive shear stress components Fig. 6.5 The element shear stress components for the to u r- element torsion problem. All values are in Newtons per square centimeter. The shear stress values calculated for each element are constant over the region of the element because the element interpolation polynomial was linear in j\: andy. The failure to obtain gradients that vary with x andy within an element is a definite disadvantage of the simplex elements. There are three ways to improve the stress values obtained for this example. First, a larger number of elements can be used. As the size of the element decreases, the existence of a constant value within the element becomes more realistic. An alternative approach is to use a triangular element with more nodes and an interpolating polynomial with quadratic and cubic terms. Differentiation will then yield gradients that are a function of the coordinate directions. A third approach is to utilize conjugate approximation theory. This theory makes possible the determination of the stress values at the nodal points and as a function of x and y within the element. The implementation of this theory is discussed in the next section. Another resultant of interest is the magnitude of the twisting torque T, which is given in (6.4) as This integral is equivalent to ^=2f (i>da ''Area =2/ Ȧ rea da (6.22)
14 Conventional Element Resultants 99 where the equations for are given in (6.9). Starting with element one if <l>^'^da=2f fyvl" /VjO A'i" Ol J^O) $2 <1.3 da (6.23) or if <p('>da=l[^vf [N^'^VdA Ja(^) (6.24) which is identical to the integral in (6.16). We can immediately conclude that if <#,('>^^ = ^[<1)] ^(1) 2A (1) -(, + (^2 + ^4) (6.25) Substituting the nodal values yields 2 I = ^4 ( )= ^4 3 3 (51.79) Likewise, for the other elements, 2 2A^^^ = ^($2 + <E'3 + ^5)= ^(16) if < > ^(^>2 + <f5+ <i>4) = ^(283.63) Ja O) J r (A\ 2A^^^ 2A^^^ 2jy^dA = ^(4.4 +$5 + $,)= ^(123.63) Summing these relationships and noting that the area of each element is the
15 1 Torsion of Nonclrcuiar Sections same produces E T=y f = ^( ) r= ^(169.5)=^^ We must also recall that the region we have subdivided into elements is only one eighth of the total cross section; therefore, the total twisting moment is (169) Moment = 87=8 = N-cm 48 It takes a torque of 178 N-cm to produce a twist of 1 deg in a 1-cm square steel shaft that is 1 cm long. The accuracy of this result is highly questionable, however, because of the coarseness of the grid. In fact, the theoretical value* is N-cm. Our answer is 9.5 percent below this value. 6.4 CONSISTENT ELEMENT RESULTANTS The failure to obtain gradients that are a function of ;c and y is a definite disadvantage of the linear interpolation functions. The gradient and any quant ity related to the gradient has a constant value within an element. Various averaging techniques have been used to obtain a value of the desired quantity at a node. One procedure uses the average of the gradient values in each element that surrounds the node. The correct procedure for obtaining nodal values of the element resultants is to employ conjugate approximation theory (Oden and Brauchli, 1971). This theory yields element resultant values consistent with the approximating polynomial(s) for a scalar or vector quantity.. The derivation of the conjugate approximation theory is beyond the scope of this book. The implementation of this theory, however, is straightforward and will be illustrated relative to the four-element torsion problem that we have been analyzing. The nodal values of the element resultants are obtained by solving the system of equations [C]{a} = {7?} (6.26) *The relationship between the applied torque and the angle of twist for a square of dimension 2a is given by 7=.146 (la)"^(timoshenko and Goodier, 197, equation 17, page 313). For our example, 2a = 1 and 7=.146 Gff= N-cm.
16 Consistent Element Resultants 11 where [C] and (/?) are the sum of element matrices defined by [cw]= r J y (6.27) and [/^n= f Jy (6.28) where a is the conventional element resultant. The evaluation of presents no problems because it is similar to (6.8), since a is constant within an element. The evaluation of is performed most easily using area coordinates, Confining our attention to element one, yields [,yv<'>] = [L, L, L, ] The product becomes L.Lj [;V<')]^[7V<'>] = L,L, Ll ^2^3 L1L3 L2L3 Ll Integration over the area using (3.43) yields ['"'1 = 1? (6.29) assuming a unit thickness. The column vector for element one is {,(.)j (1)1 = r r '77.67" = ^( (6.3)
17 12 Torsion of Nonclrcular Sections The shear stress component was used in the previous calculations because it has the largest numerical value within each element. The shear stress com ponent could also be considered. The only quantity that changes is the column vector (/?), which would have to be reevaluated using the numerical values for The element matrices for each element are summarized in (6.29) to (6.33). [c ] = 4(2) w ' ,. (6.31) [C(3)].(3)1 = ^ 12 " ' = ^(3) (6.32) 4(4) IT ' ' (6.33) The final system of equations is obtained by summing the element equations which yields J_ ' < < < (6.34) The element areas cancel since all elements have the same area.
18 Consistent Element Resultants 13 The nodal resultant values are (a) ^ = [ 7.9,436.5,724.1,353.6,671.4,475.5 ] It is also possible to obtain equations for the element resultants that give the variation of these values over the element area. There is a very limited use for this information and the details will not be presented here. The use of conjugate approximation theory requires the solution of another system of algebraic equations that is identical in size to those used to obtain the nodal values. This is a definite disadvantage when we are solving problems involving a large number of elements. An alternative to solving the complete system of equations is discussed in Chapter 7 when we reconsider the torsion of a square shaft with an increased number of elements. The final steps of the finite element method have been illustrated in this chapter while we solved a specific problem. We have shown how the element matrices are evaluated and how the element resultants are determined once the nodal values are known. We shall discuss some of the computer techniques used when we implement the finite element method in the next chapter. We shall consider specific application areas in Chapters 8 to 12. PROBLEMS 41. Verify the element stiffness and force matrix in equation (a) 6.19b, (b) 6.19c, and (c) 6.19d. 42. Verify the shear stress values for element (a) two, (b) three, and (c) four of the example discussed in Section Determine the numerical values in {R} for the consistent element re sultants when the values of are desired at the nodes. References Fung, Y. C., Foundations of Solid Mechanics, Prentice-Hall, Englewood Cliffs, N.J., 1965, pp Oden, J. T., and H. J. Brauchli, "On the Calculation of Consistent Stress Distributions in Finite Element Approximations," International Journal for Numerical Methods in Engineering, 1971, Vol. 3, pp Timoshenko, S. P., and J. N. Goodier, Theory of Elasticity, McGraw-Hill, New York, 197, pp
19 14 Torsion of Nonclrcuiar Sections Uncited Reference Owen, D. R. J., and O. C. Zienkiewicz, "Torsion of Axi-Symmetric Solids of Variable Diameter Including Acceleration Effects," International Journal for Numerical Methods in Engineering, 1974, Vol. 8, pp Yamada, Y., T. Kawai, and N. Yoshemura, "Analysis of the Elastic-Plastic- Problems by the Matrix Displacement Method," Proceedings of the Second Conference on Matrix Methods in Structural Mechanics (AFFDL-TR-68-15), Wright-Patterson Air Force Base, Dayton, Ohio, 1968.
Stress 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 informationMEI solutions to exercise 4 1
MEI-55 solutions to exercise Problem Solve the stationary two-dimensional heat transfer problem shown in the figure below by using linear elements Use symmetry to reduce the problem size The material is
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 informationONE - DIMENSIONAL INTEGRATION
Chapter 4 ONE - DIMENSIONA INTEGRATION 4. Introduction Since the finite element method is based on integral relations it is logical to expect that one should strive to carry out the integrations as efficiently
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 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 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 informationMechanical Design in Optical Engineering
Torsion Torsion: Torsion refers to the twisting of a structural member that is loaded by couples (torque) that produce rotation about the member s longitudinal axis. In other words, the member is loaded
More informationTORSIONAL RIGIDITY OF WOOD COMPOSITE I-JOISTS Daniel Hindman. Harvey B. Manbeck. John J. Janowiak
TORSIONAL RIGIDITY OF WOOD COMPOSITE I-JOISTS Daniel Hindman Assistant Professor Department of Wood Science and Forest Products Virginia Polytechnic Institute and State University Brooks Forest Products
More informationFinite Element Analysis of Saint-Venant Torsion Problem with Exact Integration of the Elastic-Plastic Constitutive
Finite Element Analysis of Saint-Venant Torsion Problem with Exact Integration of the Elastic-Plastic Constitutive Equations W. Wagner Institut für Baustatik Universität Karlsruhe (TH) Kaiserstraße 12
More informationLinearized theory of elasticity
Linearized theory of elasticity Arie Verhoeven averhoev@win.tue.nl CASA Seminar, May 24, 2006 Seminar: Continuum mechanics 1 Stress and stress principles Bart Nowak March 8 2 Strain and deformation Mark
More informationLecture 7: The Beam Element Equations.
4.1 Beam Stiffness. A Beam: A long slender structural component generally subjected to transverse loading that produces significant bending effects as opposed to twisting or axial effects. MECH 40: Finite
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 informationQuintic beam closed form matrices (revised 2/21, 2/23/12) General elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationFlexural-Torsional Buckling of General Cold-Formed Steel Columns with Unequal Unbraced Lengths
Proceedings of the Annual Stability Conference Structural Stability Research Council San Antonio, Texas, March 21-24, 2017 Flexural-Torsional Buckling of General Cold-Formed Steel Columns with Unequal
More information4. SHAFTS. A shaft is an element used to transmit power and torque, and it can support
4. SHAFTS A shaft is an element used to transmit power and torque, and it can support reverse bending (fatigue). Most shafts have circular cross sections, either solid or tubular. The difference between
More informationMECHANICS OF MATERIALS
2009 The McGraw-Hill Companies, Inc. All rights reserved. Fifth SI Edition CHAPTER 3 MECHANICS OF MATERIALS Ferdinand P. Beer E. Russell Johnston, Jr. John T. DeWolf David F. Mazurek Torsion Lecture Notes:
More informationChapter 5 Torsion STRUCTURAL MECHANICS: CE203. Notes are based on Mechanics of Materials: by R. C. Hibbeler, 7th Edition, Pearson
STRUCTURAL MECHANICS: CE203 Chapter 5 Torsion Notes are based on Mechanics of Materials: by R. C. Hibbeler, 7th Edition, Pearson Dr B. Achour & Dr Eng. K. El-kashif Civil Engineering Department, University
More informationGeneral elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
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 informationFurther Linear Elasticity
Torsion of cylindrical bodies Further Linear Elasticity Problem Sheet # 1. Consider a cylindrical body of length L, the ends of which are subjected to distributions of tractions that are statically equivalent
More informationLecture 2: Finite Elements
Materials Science & Metallurgy Master of Philosophy, Materials Modelling, Course MP7, Finite Element Analysis, H. K. D. H. Bhadeshia Lecture 2: Finite Elements In finite element analysis, functions of
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 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 informationStress Analysis Lecture 3 ME 276 Spring Dr./ Ahmed Mohamed Nagib Elmekawy
Stress Analysis Lecture 3 ME 276 Spring 2017-2018 Dr./ Ahmed Mohamed Nagib Elmekawy Axial Stress 2 Beam under the action of two tensile forces 3 Beam under the action of two tensile forces 4 Shear Stress
More informationMMJ1153 COMPUTATIONAL METHOD IN SOLID MECHANICS PRELIMINARIES TO FEM
B Course Content: A INTRODUCTION AND OVERVIEW Numerical method and Computer-Aided Engineering; Phsical problems; Mathematical models; Finite element method;. B Elements and nodes, natural coordinates,
More informationVIBRATION PROBLEMS IN ENGINEERING
VIBRATION PROBLEMS IN ENGINEERING FIFTH EDITION W. WEAVER, JR. Professor Emeritus of Structural Engineering The Late S. P. TIMOSHENKO Professor Emeritus of Engineering Mechanics The Late D. H. YOUNG Professor
More informationDiscontinuous Distributions in Mechanics of Materials
Discontinuous Distributions in Mechanics of Materials J.E. Akin, Rice University 1. Introduction The study of the mechanics of materials continues to change slowly. The student needs to learn about software
More informationChapter Objectives. Copyright 2011 Pearson Education South Asia Pte Ltd
Chapter Objectives To determine the torsional deformation of a perfectly elastic circular shaft. To determine the support reactions when these reactions cannot be determined solely from the moment equilibrium
More informationChapter 1 General Introduction Instructor: Dr. Mürüde Çelikağ Office : CE Building Room CE230 and GE241
CIVL222 STRENGTH OF MATERIALS Chapter 1 General Introduction Instructor: Dr. Mürüde Çelikağ Office : CE Building Room CE230 and GE241 E-mail : murude.celikag@emu.edu.tr 1. INTRODUCTION There are three
More informationMECHANICS OF MATERIALS
GE SI CHAPTER 3 MECHANICS OF MATERIALS Ferdinand P. Beer E. Russell Johnston, Jr. John T. DeWolf David F. Mazurek Torsion Lecture Notes: J. Walt Oler Texas Tech University Torsional Loads on Circular Shafts
More informationExact and Numerical Solution of Pure Torsional Shaft
Australian Journal of Basic and Applied Sciences, 4(8): 3043-3052, 2010 ISSN 1991-8178 Exact and Numerical Solution of Pure Torsional Shaft 1 Irsyadi Yani, 2 M.A Hannan, 1 Hassan Basri, and 2 E. Scavino
More informationLecture 8: Assembly of beam elements.
ecture 8: Assembly of beam elements. 4. Example of Assemblage of Beam Stiffness Matrices. Place nodes at the load application points. Assembling the two sets of element equations (note the common elemental
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 informationComb Resonator Design (2)
Lecture 6: Comb Resonator Design () -Intro. to Mechanics of Materials Sh School of felectrical ti lengineering i and dcomputer Science, Si Seoul National University Nano/Micro Systems & Controls Laboratory
More informationRetrospectives on My Studies of Solid Mechanics (II)
Retrospectives on My Studies of Solid Mechanics (II) - a new energy method based on the stationary principle of total energy By Tadahiko Kawai + and Etsu Kazama ++ ABSTRACT A new energy method is proposed
More informationCIVL 7/8117 Chapter 4 - Development of Beam Equations - Part 2 1/34. Chapter 4b Development of Beam Equations. Learning Objectives
CIV 7/87 Chapter 4 - Development of Beam Equations - Part /4 Chapter 4b Development of Beam Equations earning Objectives To introduce the work-equivalence method for replacing distributed loading by a
More informationOn the torsion of functionally graded anisotropic linearly elastic bars
IMA Journal of Applied Mathematics (2007) 72, 556 562 doi:10.1093/imamat/hxm027 Advance Access publication on September 25, 2007 edicated with admiration to Robin Knops On the torsion of functionally graded
More informationApplications in Fluid Mechanics
CHAPTER 8 Applications in Fluid 8.1 INTRODUCTION The general topic of fluid mechanics encompasses a wide range of problems of interest in engineering applications. The most basic definition of a fluid
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 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 informationComb resonator design (2)
Lecture 6: Comb resonator design () -Intro Intro. to Mechanics of Materials School of Electrical l Engineering i and Computer Science, Seoul National University Nano/Micro Systems & Controls Laboratory
More informationSteps in the Finite Element Method. Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen
Steps in the Finite Element Method Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen General Idea Engineers are interested in evaluating effects such as deformations, stresses,
More informationSection 6: PRISMATIC BEAMS. Beam Theory
Beam Theory There are two types of beam theory aailable to craft beam element formulations from. They are Bernoulli-Euler beam theory Timoshenko beam theory One learns the details of Bernoulli-Euler beam
More informationTorsion/Axial Illustration: 1 (3/30/00)
Torsion/Axial Illustration: 1 (3/30/00) Table of Contents Intro / General Strategy Axial: Different Materia The Displacement Method 1 2 Calculate the Stresses General Strategy The same structure is loaded
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 06
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 06 In the last lecture, we have seen a boundary value problem, using the formal
More informationHaus, Hermann A., and James R. Melcher. Electromagnetic Fields and Energy. Englewood Cliffs, NJ: Prentice-Hall, ISBN:
MIT OpenCourseWare http://ocw.mit.edu Haus, Hermann A., and James R. Melcher. Electromagnetic Fields and Energy. Englewood Cliffs, NJ: Prentice-Hall, 989. ISBN: 978032490207. Please use the following citation
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 informationLab Exercise #3: Torsion
Lab Exercise #3: Pre-lab assignment: Yes No Goals: 1. To evaluate the equations of angular displacement, shear stress, and shear strain for a shaft undergoing torsional stress. Principles: testing of round
More informationUNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES
UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES A Thesis by WOORAM KIM Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the
More informationLaminated Beams of Isotropic or Orthotropic Materials Subjected to Temperature Change
United States Department of Agriculture Forest Service Forest Products Laboratory Research Paper FPL 375 June 1980 Laminated Beams of Isotropic or Orthotropic Materials Subjected to Temperature Change
More information5. What is the moment of inertia about the x - x axis of the rectangular beam shown?
1 of 5 Continuing Education Course #274 What Every Engineer Should Know About Structures Part D - Bending Strength Of Materials NOTE: The following question was revised on 15 August 2018 1. The moment
More informationUnit 15 Shearing and Torsion (and Bending) of Shell Beams
Unit 15 Shearing and Torsion (and Bending) of Shell Beams Readings: Rivello Ch. 9, section 8.7 (again), section 7.6 T & G 126, 127 Paul A. Lagace, Ph.D. Professor of Aeronautics & Astronautics and Engineering
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationDispersion relation for transverse waves in a linear chain of particles
Dispersion relation for transverse waves in a linear chain of particles V. I. Repchenkov* It is difficult to overestimate the importance that have for the development of science the simplest physical and
More informationA BRIEF HISTORY OF THE BEGINNING OF THE FINITE ELEMENT METHOD
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, VOL. 39, 3761-3774 (1996) A BRIEF HISTORY OF THE BEGINNING OF THE FINITE ELEMENT METHOD K. K. GUPTA' Integrated System, NASA, Dryden Flight Research
More information3. Numerical integration
3. Numerical integration... 3. One-dimensional quadratures... 3. Two- and three-dimensional quadratures... 3.3 Exact Integrals for Straight Sided Triangles... 5 3.4 Reduced and Selected Integration...
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 informationFIRST INTERNATIONAL SEMINAR DEEP AND HIGH STRESS MINING 6-8 NOVEMBER 2002 PERTH, AUSTRALIA. Potential. T. Wiles Mine Modelling Pty Ltd, Australia
FIRST INTERNATIONAL SEMINAR ON DEEP AND HIGH STRESS MINING 6-8 NOVEMBER 22 PERTH, AUSTRALIA Loading System Stiffness A Parameter to Evaluate Rockburst Potential T. Wiles Mine Modelling Pty Ltd, Australia
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 of Thin Beams by PIM and RPIM methods. *B. Kanber¹, and O. M. Tufik 1
APCOM & ISCM -4 th December, 23, Singapore Vibration of Thin Beams by PIM and RPIM methods *B. Kanber¹, and O. M. Tufik Mechanical Engineering Department, University of Gaziantep, Turkey. *Corresponding
More informationSupplementary figures
Supplementary figures Supplementary figure 1: The fang model (longitudinal section) used for the stiffness analysis. The base of the fang is considered fixed (i.e. zero displacement constraint) and the
More informationJeff Brown Hope College, Department of Engineering, 27 Graves Pl., Holland, Michigan, USA UNESCO EOLSS
MECHANICS OF MATERIALS Jeff Brown Hope College, Department of Engineering, 27 Graves Pl., Holland, Michigan, USA Keywords: Solid mechanics, stress, strain, yield strength Contents 1. Introduction 2. Stress
More informationWorkshop 8. Lateral Buckling
Workshop 8 Lateral Buckling cross section A transversely loaded member that is bent about its major axis may buckle sideways if its compression flange is not laterally supported. The reason buckling occurs
More informationChapter 5: Torsion. 1. Torsional Deformation of a Circular Shaft 2. The Torsion Formula 3. Power Transmission 4. Angle of Twist CHAPTER OBJECTIVES
CHAPTER OBJECTIVES Chapter 5: Torsion Discuss effects of applying torsional loading to a long straight member (shaft or tube) Determine stress distribution within the member under torsional load Determine
More informationEMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 3 Torsion
EMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 3 Torsion Introduction Stress and strain in components subjected to torque T Circular Cross-section shape Material Shaft design Non-circular
More informationCE 221: MECHANICS OF SOLIDS I CHAPTER 1: STRESS. Dr. Krisada Chaiyasarn Department of Civil Engineering, Faculty of Engineering Thammasat university
CE 221: MECHANICS OF SOLIDS I CHAPTER 1: STRESS By Dr. Krisada Chaiyasarn Department of Civil Engineering, Faculty of Engineering Thammasat university Agenda Introduction to your lecturer Introduction
More informationReview of Strain Energy Methods and Introduction to Stiffness Matrix Methods of Structural Analysis
uke University epartment of Civil and Environmental Engineering CEE 42L. Matrix Structural Analysis Henri P. Gavin Fall, 22 Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods
More informationUsing the finite element method of structural analysis, determine displacements at nodes 1 and 2.
Question 1 A pin-jointed plane frame, shown in Figure Q1, is fixed to rigid supports at nodes and 4 to prevent their nodal displacements. The frame is loaded at nodes 1 and by a horizontal and a vertical
More information[7] Torsion. [7.1] Torsion. [7.2] Statically Indeterminate Torsion. [7] Torsion Page 1 of 21
[7] Torsion Page 1 of 21 [7] Torsion [7.1] Torsion [7.2] Statically Indeterminate Torsion [7] Torsion Page 2 of 21 [7.1] Torsion SHEAR STRAIN DUE TO TORSION 1) A shaft with a circular cross section is
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 informationChapter 10 INTEGRATION METHODS Introduction Unit Coordinate Integration
Chapter 0 INTEGRATION METHODS 0. Introduction The finite element analysis techniques are always based on an integral formulation. At the very minimum it will always be necessary to integrate at least an
More informationSolving the torsion problem for isotropic matrial with a rectangular cross section using the FEM and FVM methods with triangular elements
Solving the torsion problem for isotropic matrial with a rectangular cross section using the FEM and FVM methods with triangular elements Nasser M. Abbasi. June 0, 04 Contents Introduction. Problem setup...................................
More informationWORCESTER POLYTECHNIC INSTITUTE
WORCESTER POLYTECHNIC INSTITUTE MECHANICAL ENGINEERING DEPARTMENT STRESS ANALYSIS ES-2502, C 2012 Lecture 17: 10 February 2012 General information Instructor: Cosme Furlong HL-151 (508) 831-5126 cfurlong@wpi.edu
More informationTORSION TEST. Figure 1 Schematic view of torsion test
TORSION TEST 1. THE PURPOSE OF THE TEST The torsion test is performed for determining the properties of materials like shear modulus (G) and shear yield stress ( A ). 2. IDENTIFICATIONS: Shear modulus:
More informationAN INTEGRATED KIRCHHOFF PLATE ELEMENT BY GALERKIN METHOD FOR THE ANALYSIS OF PLATES ON ELASTIC FOUNDATION
AN INTEGRATED IRCHHOFF PLATE ELEMENT BY GALERIN METHOD FOR THE ANALYSIS OF PLATES ON ELASTIC FOUNDATION Ragesh.P.P, V.Mustafa 2, T.P.Somasundaran 3 (Research Scholar, National Institute of Technology Calicut,
More informationCIVL222 STRENGTH OF MATERIALS. Chapter 6. Torsion
CIVL222 STRENGTH OF MATERIALS Chapter 6 Torsion Definition Torque is a moment that tends to twist a member about its longitudinal axis. Slender members subjected to a twisting load are said to be in torsion.
More informationFLEXIBILITY METHOD FOR INDETERMINATE FRAMES
UNIT - I FLEXIBILITY METHOD FOR INDETERMINATE FRAMES 1. What is meant by indeterminate structures? Structures that do not satisfy the conditions of equilibrium are called indeterminate structure. These
More informationFINITE ELEMENT METHOD FOR TWO-DIMENSIONAL ELASTICITY PROBLEM HANIMAH OTHMAN
FINITE ELEMENT METHOD FOR TWO-DIMENSIONAL ELASTICITY PROBLEM HANIMAH OTHMAN A dissertation submitted in partial fulfilment of the requirements for the award of the degree of Master of Science (Mathematics)
More informationDynamic Analysis of Coupling Vehicle-Bridge System Using Finite Prism Method
Dynamic Analysis of Coupling Vehicle-Bridge System Using Finite Prism Method A. T. Saeed and Zhongfu Xiang Abstract To investigate the transient responses of bridges under moving vehicles, Finite Prism
More informationEngineering Sciences 241 Advanced Elasticity, Spring Distributed Thursday 8 February.
Engineering Sciences 241 Advanced Elasticity, Spring 2001 J. R. Rice Homework Problems / Class Notes Mechanics of finite deformation (list of references at end) Distributed Thursday 8 February. Problems
More informationTwo hours. To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER. 29 May :45 11:45
Two hours MATH20602 To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER NUMERICAL ANALYSIS 1 29 May 2015 9:45 11:45 Answer THREE of the FOUR questions. If more
More informationChapter 12 Plate Bending Elements. Chapter 12 Plate Bending Elements
CIVL 7/8117 Chapter 12 - Plate Bending Elements 1/34 Chapter 12 Plate Bending Elements Learning Objectives To introduce basic concepts of plate bending. To derive a common plate bending element stiffness
More informationCOMPUTATIONAL ELASTICITY
COMPUTATIONAL ELASTICITY Theory of Elasticity and Finite and Boundary Element Methods Mohammed Ameen Alpha Science International Ltd. Harrow, U.K. Contents Preface Notation vii xi PART A: THEORETICAL ELASTICITY
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 informationFIXED BEAMS IN BENDING
FIXED BEAMS IN BENDING INTRODUCTION Fixed or built-in beams are commonly used in building construction because they possess high rigidity in comparison to simply supported beams. When a simply supported
More informationVORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS
The 4 th World Conference on Earthquake Engineering October -7, 008, Beijing, China VORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS K.
More information3D Stress Tensors. 3D Stress Tensors, Eigenvalues and Rotations
3D Stress Tensors 3D Stress Tensors, Eigenvalues and Rotations Recall that we can think of an n x n matrix Mij as a transformation matrix that transforms a vector xi to give a new vector yj (first index
More informationGyroscopic matrixes of the straight beams and the discs
Titre : Matrice gyroscopique des poutres droites et des di[...] Date : 29/05/2013 Page : 1/12 Gyroscopic matrixes of the straight beams and the discs Summarized: This document presents the formulation
More information4.5 The framework element stiffness matrix
45 The framework element stiffness matri Consider a 1 degree-of-freedom element that is straight prismatic and symmetric about both principal cross-sectional aes For such a section the shear center coincides
More informationChapter 6: Momentum Analysis
6-1 Introduction 6-2Newton s Law and Conservation of Momentum 6-3 Choosing a Control Volume 6-4 Forces Acting on a Control Volume 6-5Linear Momentum Equation 6-6 Angular Momentum 6-7 The Second Law of
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 informationEarthquake Loads According to IBC IBC Safety Concept
Earthquake Loads According to IBC 2003 The process of determining earthquake loads according to IBC 2003 Spectral Design Method can be broken down into the following basic steps: Determination of the maimum
More informationInternational Journal of Advanced Engineering Technology E-ISSN
Research Article INTEGRATED FORCE METHOD FOR FIBER REINFORCED COMPOSITE PLATE BENDING PROBLEMS Doiphode G. S., Patodi S. C.* Address for Correspondence Assistant Professor, Applied Mechanics Department,
More informationA short review of continuum mechanics
A short review of continuum mechanics Professor Anette M. Karlsson, Department of Mechanical ngineering, UD September, 006 This is a short and arbitrary review of continuum mechanics. Most of this material
More information(48) CHAPTER 3: TORSION
(48) CHAPTER 3: TORSION Introduction: In this chapter structural members and machine parts that are in torsion will be considered. More specifically, you will analyze the stresses and strains in members
More informationTorsion. Click here to check the Animation. Introduction. Basic Assumptions. Assumptions. Torsion Formula. Stress Formula. Stresses on Inclined Planes
Torsion ntroduction Basic Assumtions Torsion Formula Stresses on nclined Planes Angle of Twist in Torsion Torsion of Circular Elastic Bars: Formulae Click here to check the Animation Assumtions Stress
More information6.1 Formulation of the basic equations of torsion of prismatic bars (St. Venant) Figure 6.1: Torsion of a prismatic bar
Module 6 Torsion Learning Objectives 6.1 Formulation of the basic equations of torsion of prismatic bars (St. Venant) Readings: Sadd 9.3, Timoshenko Chapter 11 e e 1 e 3 Figure 6.1: Torsion of a prismatic
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 informationUNIT-I Introduction & Plane Stress and Plane Strain Analysis
SIDDHARTH INSTITUTE OF ENGINEERING & TECHNOLOGY:: PUTTUR (AUTONOMOUS) Siddharth Nagar, Narayanavanam Road 517583 QUESTION BANK (DESCRIPTIVE) Subject with Code : Advanced Solid Mechanics (18CE1002) Year
More information