The thin plate theory assumes the following deformational kinematics:
|
|
- Clarence Lewis
- 6 years ago
- Views:
Transcription
1 MEG6007 (Fall, 2017) - Solutions of Homework # 8 (ue on Tuesay, 29 November 2017) 8.0 Variational Derivation of Thin Plate Theory The thin plate theory assumes the following eformational kinematics: u = u 0 (x, y) z x v = v 0 (x, y) z y w = w(x, y) (1) where u, v an w are the isplacements along the x, y an z coorinates; u 0 (x, y) an v 0 (x, y) represent the stretching the plate neutral surface; z is measure from the neutral surface of the plate, respectively. The strain-isplacement relations are obtaine as ɛ xx = u x = u 0 x z 2 w x 2 ɛ yy = v x = u 0 x z 2 w y 2 γ xy = v x + u y = v 0 x + u 0 y 2z x y γ xz = γ yz = 0 (2) For an isotropic material, the constitutive relations are moele by ɛ xx = 1 E (σ xx νσ yy ), ɛ yy = 1 E (σ yy νσ xx ), τ xy = Gγ xy (3) σ xx = E 1 ν (ɛ 2 xx + νɛ yy ), σ yy = E 1 ν (ɛ 2 yy + νɛ xx ), where σ xx, σ xx an τ xy are stress components; E an G are Young s moulus an shear moulus; an, ν is Poisson s ratio. Variational formulation of plate vibrations may be formulate employing the following expression: The bening strain energy of a plate can be expresse as δv b = V { E 1 ν 2 δɛ xx (ɛ xx +νɛ yy )+ E 1 ν 2 δɛ yy (ɛ yy +νɛ xx )+δγ xy G ɛ xy }V (4) 1
2 When the plate is sufficiently thin an the applie loas o no cause the neutral surface to stretch, we have u 0 = v 0 = 0, (5) 8.1 Carry out the necessary variation of the strain energy given in equation (4), using the expression D = Eh 3 12(1 ν 2 ) where D is calle plate bening rigiity, h is the plate thickness. an is the plate area. 8.2 Derive the kinetic energy an the external energy assuming that the applie force f(x, y, t) is istribute on the plate surface. 8.3 Obtain Hamilton s principle for a rectangular thin plate. n Ientify the bounary conitions from your variations. 8.4 Obtain the necessary characteristic equation for vibration frequencies assuming the following homogeneous solution form: w(x, y, t) = W (x, y)t (t), T (t) = e ωt W (x, y) = 1 sin αx sin γy + 2 sin αx cos γy + 3 cos αx sin γy + 4 cos αx cos γy (7) + 5 sinh α 1 x sinh γ 1 y + 6 sinh α 1 x cosh γ 1 y + 7 cosh α 1 x sinh γ 1 y + 8 cosh α 1 x cosh γ 1 y β 4 = ω2 m D, α 2 + γ 2 = β 2 α γ 2 1 = β 2 where m is the plate mass per unit area, m = ρh an obtain the characteristic equation for a simply supporte plate case. - Solution: (6) Substituting equations(1) -(3) into (4), carrying out the through-the-thickness 2
3 θy w y 4 b 3 θ x a 1 2 x Fig. 1. Plate geometry. Note that θ x = w y, θ y = w x spatial integration an imposing the thin plate assumption (5), one obtains: δv b = D { [ ( 2 w x + 2 w 2 y ) (1 ν) 2 w 2 y ] δ 2 w 2 x 2 + [ ( 2 w x + 2 w 2 y ) (1 ν) 2 w 2 x ] δ 2 w 2 y 2 + 2(1 ν) 2 w x y δ 2 w x y } D = Eh3 12(1 ν 2 ) (8) where spatial integration across the plate thickness z is carrie out, an where D is calle plate bening stiffness, h is the plate thickness, an is the plate area. 3
4 In carrying out the inicate variations, we make use of the ientities: G(w) δ x Γ = l G(w)δw Γ G(w) δ y Γ = m G(w)δw Γ G(w) x δw G(w) δw (9) y where l an m are the irectional cosines of the unit normal vector outwar from the surface contour Γ. Hence, δv b can be shown to be δv b = [l D( 2 w Γ x + ν 2 w ) + m (1 ν) D 2 w 2 y2 x y ] δ x Γ + [m D( 2 w Γ y + ν 2 w ) + l (1 ν) D 2 w 2 x2 x y ] δ y Γ [l D Γ x ( 2 w x + ν 2 w ) + l(1 ν) D 2 w ] δw Γ 2 y2 x y2 [m D Γ y ( 2 w y + ν 2 w ) + m(1 ν) D 2 2 x2 + D 4 w δw, 4 = 4 x x 2 y y 4 (11) w ] δw Γ (10) x 2 y In orer to ientify the natural an essential bounary conitions, we introuce the following efinitions: M xx = D( 2 w y 2 ) M yy = D( 2 w y 2 + ν 2 w x 2 ) M xy = D(1 ν) Q x = M xx x Q y = M yy y + M xy y + M xy x x y 2 (12) 4
5 so that (10) simplifies to δv b = Γ + + Γ { [l M xx + m M xy ] δ x + [m M yy + l M xy ] δ y } Γ { [l x M xx + l y M xy] + [m y M yy + m x M xy] } δw Γ D 4 w δw (13) In orer for the above expression to be vali for all possible bounary contour shapes incluing a circular plate, we introuce the following transformation (see Fig.(2)) x = l n m s y = m n + l, l = cos φ, m = sin φ (14) s Substituting the transformation (14) into (13) we obtain δv b = [ M n δ s n M ns δ s + Q n δw ] s + D 4 w δw M n =l 2 M xx + 2lm M xy + m 2 M yy M ns =lm (M yy M xx ) + (l 2 m 2 ) M xy Q n =l Q x + m Q y 2 = 2 n R n + 2 s, 1 2 R = φ s (15) Note that via (9) we have s M ns δ s s = M ns δw Γ s M ns δw s (16) s If the bounary contour Γ is smooth an close, the virtual work one by the ege twisting moment M ns shoul vanish. Otherwise, a gap along the contour woul exist. Remark: For rectangular plates, there exist iscontinuities at four corners, known as corner conitions, which states that M xy = 0. This conition must be satisfie either explicitly or implicitly. Thus, for a smooth close contour the variation of the strain energy becomes 5
6 Fig. 2. Tangential an normal irections at a plate bounary δv b = [ M n δ s n + (Q n + M ns s )δw ] s + D 4 w δw (17) The kinetic energy δt, the energy ue to the bounary unknown springs δv s, an the work performe by the applie forces δw can be expresse in a similar manner as for the case of beam in the form of t2 t 1 δt t = t2 m(x, y) ẅ(x, y, t) δw(x, t) t (18) t 1 where m(x, y) = ρh is the plate mass per unit plate area, an the kinetic energy ue to in-plane motions (u = zw xx, v = zw yy ) is neglecte as it is 6
7 associate with h 3 compare to the plate size. The variation of the potential energy of the unknown bounary forces an moments are given by δv s = s [ k w w δw + k θ n δ n ] s (19) The work performe by the applie force δw = f(x, y, t)δw(x, y, t) (20) Hamilton s principle for a continuum plate can be written as t2 t 1 + { s (D 4 w + mẅ f) δw [ (V n + k w w ) δw + (M n k θ n ) δ n ] s } t = 0 (21) where (V n, M n ) are given by V n = D n 2 w (1 ν) D M n = D 2 w + (1 ν) D( 1 R s ( 2 w n s 1 R s ) s + 2 w s ) 2 (22) The governing equation of motion for a plate can thus be obtaine as with the bounary conitions: D 4 w + mẅ = f (23) (V n + k w w ) = 0 (M n k θ n ) = 0 (24) Free Vibrations of a Rectangular Plate Free vibrations of a rectangular plate can be stuie by specializing the bounary conitions (24) accoring to the convention shown in Fig. (3). First, the solution of the homogeneous equation of motion for the plate with f = 0 in (23) may assume 7
8 y n = y, s = -x n = -x, s = -y b a n = x, s = y n = -y, s= x x Fig. 3. Bounary escriptions for a rectangular plate w(x, y, t) = W (x, y)e jωt ( 4 β 4 )W (x, y) = 0, β 4 = ω2 m D ( 2 + β 2 )( 2 β 2 )W (x, y) = 0 (25) 8
9 Hence, its solution consists of ( 2 + β 2 )W 1 = 0 W 1 = e jαx e jγy, α 2 + γ 2 = β 2 ( 2 β 2 )W 2 = 0 W 2 = e α1x e γ1y, α1 2 + γ1 2 = β 2 W = W 1 + W 2 W (x, y) = 1 sin αx sin γy + 2 sin αx cos γy + 3 cos αx sin γy + 4 cos αx cos γy + 5 sinh α 1 x sinh γ 1 y + 6 sinh α 1 x cosh γ 1 y + 7 cosh α 1 x sinh γ 1 y + 8 cosh α 1 x cosh γ 1 y (26) Specializing the bounary conition(24) to a rectangular case requires the conversion of (n, s) in terms of (x, y) along the four eges as shown in Fig. (3). The bounary conitions along the four eges are obtaine as follows.. long y = 0 we have (s = x, n = y): V y = D y 2 w + (1 ν) D 3 w [D y x 2 y 2 w + (1 ν) D 3 w y x + k ww ] 2 (x,0) = 0 M y = D( 2 w y 2 + ν 2 w x 2 ) [ D( 2 w y 2 + ν 2 w x 2 ) + k θ y ] (x,0) = 0 (27) B. long y = b we have (s = x, n = y): V y = [D y 2 w + (1 ν) D 3 w y x ] [D 2 y 2 w + (1 ν) D 3 w y x k ww ] 2 (x,b) = 0 M y = D( 2 w y 2 + ν 2 w x 2 ) [ D( 2 w y 2 + ν 2 w x 2 ) k θ y ] (x,b) = 0 (28) C. long x = 0 we have (s = y, n = x): V x = [D x 2 w + (1 ν) D 3 w x y ] [D 2 x 2 w + (1 ν) D M x = D( 2 w y 2 ) [ D( 2 w y 2 ) + k θ 3 w x y + k ww ] 2 (0,y) = 0 x ] (0,y) = 0 (29) 9
10 D. long x = a we have (s = y, n = x): V x = [D x 2 w + (1 ν) D 3 w x y 2 ] M x = D( 2 w y 2 ) [D x 2 w + (1 ν) D [ D( 2 w y 2 ) k θ 3 w x y k ww ] 2 (a,y) = 0 x ] (a,y) = 0 (30) where it is unerstoo that the istribute bounary springs (k w, k θ ) are efine along the four eges. Note that there are a total of eight bounary conitions for plate vibrations. Hence, the vibration moes W (x, y) must consist of eight unknown coefficients: Simply supporte along all four eges: The bounary conitions for this case is obtaine from (27) -(30) by setting k w an k θ = 0: long y = 0 : w (x,0) = 0 an M y = D( 2 w y 2 + ν 2 w x 2 ) (x,0) = 0 long y = b : w (x,b) = 0 an M y = D( 2 w y 2 + ν 2 w x 2 ) (x,b) = 0 long x = 0 : w (0,y) = 0 an M x = D( 2 w y 2 ) (0,y) = 0 long x = a : w (a,y) = 0 an M x = D( 2 w y 2 ) (a,y) = 0 (31) Note that along y = 0 we have w (x,0) = 0, which implies that w oes not vary along the plate ege (meaning with respect to x). Hence we have long y = 0 : w (x,0) = 0 implies x 2 (x,0) = 0 2 w y 2 (x,0) = 0 = 0 (32) pplying this observation to the three remaining eges, we arrive at the following bounary conitions for a simply supporte plate: 10
11 long y = 0 : w (x,0) = 0 an long y = b : w (x,b) = 0 an long x = 0 : w (0,y) = 0 an long x = a : w (a,y) = 0 an y 2 (x,0) = 0 y 2 (x,b) = 0 x 2 (0,y) = 0 x 2 (a,y) = 0 (33) Clearly, the above bounary conitions are satisfie by retaining 1 -term in (26): W (x, y) = 1 sin αx sin γy, α 2 + γ 2 = β 2, β 4 = ω2 m(x, y) D sin αa = 0, sin γb = 0 α m a = πm, m = 1, 2,... γ m a = πn, n = 1, 2,..., α 2 + γ 2 = β 2 β 2 mn = π 2 [( m a )2 + ( n b )2 ] (34) Therefore, the natural frequencies of a simply supporte plate are given by D ω mn = βmn 2 m(x, y) = π2 [( m a )2 + ( n D b )2 ] m(x, y) (35) Note that the frequencies of a simply supporte beam are obtaine from the above frequency expression by setting { b a 0, ν = 0} {D = Eh3 /12 = EI, m(x, y) = m(x)} ω n = ( mπ a )2 EI m(x, y) (36) Completely free along all four eges: 11
12 The bounary conitions for this case is obtaine from (27) -(30) by setting k w = 0 an k θ = 0: long y = (0, b) : y [ + (2 ν) y2 2 w x 2 ] = 0 y + ν 2 w 2 x 2 = 0 (37) long x = (0, a) : x [ + (2 ν) x 2 y 2 ] = 0 ( 2 w y 2 ) = 0 (38) It shoul be emphasize that, if an approximation is introuce for w(x, y), it must implicitly satisfy M ns = M xy = 0. Otherwise, the following must be explicitly enforce: M xy = D 2 w x y = 0 at the four corners (39) 12
Vibrations of Structures
Vibrations of Structures Module V: Vibrations of Plates Lesson 37: Dynamics of Plates: Variational Formulation Contents: 1. Introduction. Mathematical Modeling Keywords: Plate model, Variational formulation,
More informationStrength Analysis of CFRP Composite Material Considering Multiple Fracture Modes
5--XXXX Strength Analysis of CFRP Composite Material Consiering Multiple Fracture Moes Author, co-author (Do NOT enter this information. It will be pulle from participant tab in MyTechZone) Affiliation
More informationEXPONENTIAL FOURIER INTEGRAL TRANSFORM METHOD FOR STRESS ANALYSIS OF BOUNDARY LOAD ON SOIL
Tome XVI [18] Fascicule 3 [August] 1. Charles Chinwuba IKE EXPONENTIAL FOURIER INTEGRAL TRANSFORM METHOD FOR STRESS ANALYSIS OF BOUNDARY LOAD ON SOIL 1. Department of Civil Engineering, Enugu State University
More informationLecture 3: Development of the Truss Equations.
3.1 Derivation of the Stiffness Matrix for a Bar in Local Coorinates. In 3.1 we will perform Steps 1-4 of Logan s FEM. Derive the truss element equations. 1. Set the element type. 2. Select a isplacement
More informationFREE VIBRATION OF AXIALLY LOADED FUNCTIONALLY GRADED SANDWICH BEAMS USING REFINED SHEAR DEFORMATION THEORY
FREE VIBRATION OF AXIALLY LOADED FUNCTIONALLY GRADED SANDWICH BEAMS USING REFINED SHEAR DEFORMATION THEORY Thuc P. Vo 1, Adelaja Israel Osofero 1, Marco Corradi 1, Fawad Inam 1 1 Faculty of Engineering
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 informationMechanics Physics 151
Mechanics Physics 151 Lecture 3 Continuous Systems an Fiels (Chapter 13) Where Are We Now? We ve finishe all the essentials Final will cover Lectures 1 through Last two lectures: Classical Fiel Theory
More information05 The Continuum Limit and the Wave Equation
Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,
More informationGyroscopic matrices of the right beams and the discs
Titre : Matrice gyroscopique es poutres roites et es i[...] Date : 15/07/2014 Page : 1/16 Gyroscopic matrices of the right beams an the iscs Summary: This ocument presents the formulation of the matrices
More informationUNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich
UNIVERSITY OF SASKATCHEWAN ME 313.3 MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich A CLOSED BOOK EXAMINATION TIME: 3 HOURS For Marker s Use Only LAST NAME (printed): FIRST
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 information3D Elasticity Theory
3D lasticity Theory Many structural analysis problems are analysed using the theory of elasticity in which Hooke s law is used to enforce proportionality between stress and strain at any deformation level.
More informationVibration Analysis of Railway Tracks Forced by Distributed Moving Loads
IJR International Journal of Railway Vol. 6, No. 4 / December 13, pp. 155-159 The Korean Society for Railway Vibration Analysis of Railway Tracks Force by Distribute Moving Loas Sinyeob Lee*, Dongkyu Kim*,
More informationarxiv: v1 [physics.flu-dyn] 8 May 2014
Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa
More informationSEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by
SEMM Mechanics PhD Preliminary Exam Spring 2014 1. Consider a two-dimensional rigid motion, whose displacement field is given by u(x) = [cos(β)x 1 + sin(β)x 2 X 1 ]e 1 + [ sin(β)x 1 + cos(β)x 2 X 2 ]e
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
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 informationBasic Equations of Elasticity
A Basic Equations of Elasticity A.1 STRESS The state of stress at any point in a loaded bo is defined completely in terms of the nine components of stress: σ xx,σ yy,σ zz,σ xy,σ yx,σ yz,σ zy,σ zx,andσ
More informationMechanical Properties of Materials
Mechanical Properties of Materials Strains Material Model Stresses Learning objectives Understand the qualitative and quantitative description of mechanical properties of materials. Learn the logic of
More informationMath 2163, Practice Exam II, Solution
Math 63, Practice Exam II, Solution. (a) f =< f s, f t >=< s e t, s e t >, an v v = , so D v f(, ) =< ()e, e > =< 4, 4 > = 4. (b) f =< xy 3, 3x y 4y 3 > an v =< cos π, sin π >=, so
More informationSeparation of Variables
Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical
More informationA simple model for the small-strain behaviour of soils
A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:
More information3-D FEM Modeling of fiber/matrix interface debonding in UD composites including surface effects
IOP Conference Series: Materials Science an Engineering 3-D FEM Moeling of fiber/matrix interface eboning in UD composites incluing surface effects To cite this article: A Pupurs an J Varna 2012 IOP Conf.
More informationNon-Equilibrium Continuum Physics TA session #10 TA: Yohai Bar Sinai Dislocations
Non-Equilibrium Continuum Physics TA session #0 TA: Yohai Bar Sinai 0.06.206 Dislocations References There are countless books about islocations. The ones that I recommen are Theory of islocations, Hirth
More informationMore Examples Of Generalized Coordinates
Slides of ecture 8 Today s Class: Review Of Homework From ecture 7 Hamilton s Principle More Examples Of Generalized Coordinates Calculating Generalized Forces Via Virtual Work /3/98 /home/djsegal/unm/vibcourse/slides/ecture8.frm
More informationTopological Sensitivity Analysis for Three-dimensional Linear Elasticity Problem
Topological Sensitivity Analysis for Three-imensional Linear Elasticity Problem A.A. Novotny, R.A. Feijóo, E. Taroco Laboratório Nacional e Computação Científica LNCC/MCT, Av. Getúlio Vargas 333, 25651-075
More informationMath Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors
Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More informationDAMAGE DETECTIONS IN NONLINEAR VIBRATING THERMALLY LOADED STRUCTURES 1
11 th National Congress on Theoretical an Applie Mechanics, 2-5 Sept. 2009, Borovets, Bulgaria DAMAGE DETECTIONS IN NONLINEAR VIBRATING THERMALLY LOADED STRUCTURES 1 E. MANOACH Institute of Mechanics,
More informationShort Intro to Coordinate Transformation
Short Intro to Coorinate Transformation 1 A Vector A vector can basically be seen as an arrow in space pointing in a specific irection with a specific length. The following problem arises: How o we represent
More informationCHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES
CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES * Governing equations in beam and plate bending ** Solution by superposition 1.1 From Beam Bending to Plate Bending 1.2 Governing Equations For Symmetric
More informationChapter 4. Electrostatics of Macroscopic Media
Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1
More informationAssessment of the Buckling Behavior of Square Composite Plates with Circular Cutout Subjected to In-Plane Shear
Assessment of the Buckling Behavior of Square Composite Plates with Circular Cutout Sujecte to In-Plane Shear Husam Al Qalan 1)*, Hasan Katkhua 1) an Hazim Dwairi 1) 1) Assistant Professor, Civil Engineering
More informationAN INTRODUCTION TO AIRCRAFT WING FLUTTER Revision A
AN INTRODUCTION TO AIRCRAFT WIN FLUTTER Revision A By Tom Irvine Email: tomirvine@aol.com January 8, 000 Introuction Certain aircraft wings have experience violent oscillations uring high spee flight.
More informationABHELSINKI UNIVERSITY OF TECHNOLOGY
ABHELSINKI UNIVERSITY OF TECHNOLOGY TECHNISCHE UNIVERSITÄT HELSINKI UNIVERSITE DE TECHNOLOGIE D HELSINKI A posteriori error analysis for the Morley plate element Jarkko Niiranen Department of Structural
More informationShape functions in 1D
MAE 44 & CIV 44 Introuction to Finite Elements Reaing assignment: ecture notes, ogan.,. Summary: Prof. Suvranu De Shape functions in D inear shape functions in D Quaratic an higher orer shape functions
More informationCrack-tip stress evaluation of multi-scale Griffith crack subjected to
Crack-tip stress evaluation of multi-scale Griffith crack subjecte to tensile loaing by using periynamics Xiao-Wei Jiang, Hai Wang* School of Aeronautics an Astronautics, Shanghai Jiao Tong University,
More informationBioen MIDTERM. Covers: Weeks 1-4, FBD, stress/strain, stress analysis, rods and beams (not deflections).
Name Bioen 326 Midterm 2013 Page 1 Bioen 326 2013 MIDTERM Covers: Weeks 1-4, FBD, stress/strain, stress analysis, rods and beams (not deflections). Rules: Closed Book Exam: Please put away all notes and
More informationA note on the Mooney-Rivlin material model
A note on the Mooney-Rivlin material moel I-Shih Liu Instituto e Matemática Universiae Feeral o Rio e Janeiro 2945-97, Rio e Janeiro, Brasil Abstract In finite elasticity, the Mooney-Rivlin material moel
More informationConsider an elastic spring as shown in the Fig.2.4. When the spring is slowly
.3 Strain Energy Consider an elastic spring as shown in the Fig..4. When the spring is slowly pulled, it deflects by a small amount u 1. When the load is removed from the spring, it goes back to the original
More informationMODIFIED HYPERBOLIC SHEAR DEFORMATION THEORY FOR STATIC FLEXURE ANALYSIS OF THICK ISOTROPIC BEAM
MODIFIED HYPERBOLIC SHEAR DEFORMATION THEORY FOR STATIC FLEXURE ANALYSIS OF THICK ISOTROPIC BEAM S. Jasotharan * and I.R.A. Weerasekera University of Moratuwa, Moratuwa, Sri Lanka * E-Mail: jasos91@hotmail.com,
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 informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationExam 2 Review Solutions
Exam Review Solutions 1. True or False, an explain: (a) There exists a function f with continuous secon partial erivatives such that f x (x, y) = x + y f y = x y False. If the function has continuous secon
More informationApplications of the Plate Membrane Theory
Chapter 2 Applications of the Plate Membrane Theory In this chapter we will give solutions for plates, which are loaded only on their edges. This implies that no distributed forces p x and p y occur, and
More information6.642, Continuum Electromechanics Prof. Markus Zahn Lecture 1: Review of Maxwell s Equations
6.64, Continuum Electromechanics Prof. Markus Zahn Lecture 1: Review of Mawell s Equations I. Mawell s Equations in Integral Form in Free pace 1. Faraay s Law C E i s = - µ H a t i Circulation of E Magnetic
More informationLQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR
5 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES LQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR Tatsunori Kaneko* an Yasuto Asano* * Department of Mechanical Engineering,
More informationStructural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian
Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian ahmadian@iust.ac.ir Dynamic Response of MDOF Systems: Mode-Superposition Method Mode-Superposition Method:
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 informationPhysics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1
Physics 5153 Classical Mechanics The Virial Theorem an The Poisson Bracket 1 Introuction In this lecture we will consier two applications of the Hamiltonian. The first, the Virial Theorem, applies to systems
More informationCrack onset assessment near the sharp material inclusion tip by means of modified maximum tangential stress criterion
Focuse on Mechanical Fatigue of Metals Crack onset assessment near the sharp material inclusion tip by means of moifie maximum tangential stress criterion Onřej Krepl, Jan Klusák CEITEC IPM, Institute
More informationExample 3.7 Consider the undeformed configuration of a solid as shown in Figure 3.60.
162 3. The linear 3-D elasticity mathematical model The 3-D elasticity model is of great importance, since it is our highest order hierarchical model assuming linear elastic behavior. Therefore, it provides
More informationUnit 13 Review of Simple Beam Theory
MIT - 16.0 Fall, 00 Unit 13 Review of Simple Beam Theory Readings: Review Unified Engineering notes on Beam Theory BMP 3.8, 3.9, 3.10 T & G 10-15 Paul A. Lagace, Ph.D. Professor of Aeronautics & Astronautics
More informationG4003 Advanced Mechanics 1. We already saw that if q is a cyclic variable, the associated conjugate momentum is conserved, L = const.
G4003 Avance Mechanics 1 The Noether theorem We alreay saw that if q is a cyclic variable, the associate conjugate momentum is conserve, q = 0 p q = const. (1) This is the simplest incarnation of Noether
More informationChapter 2 Governing Equations
Chapter 2 Governing Equations In the present an the subsequent chapters, we shall, either irectly or inirectly, be concerne with the bounary-layer flow of an incompressible viscous flui without any involvement
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 information1 Static Plastic Behaviour of Beams
1 Static Plastic Behaviour of Beams 1.1 Introduction Many ductile materials which are used in engineering practice have a considerable reserve capacity beyond the initial yield condition. The uniaxial
More informationLecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations
Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:
More information1 Boas, p. 643, problem (b)
Physics 6C Solutions to Homework Set #6 Fall Boas, p. 643, problem 3.5-3b Fin the steay-state temperature istribution in a soli cyliner of height H an raius a if the top an curve surfaces are hel at an
More informationHomework 7 Due 18 November at 6:00 pm
Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine
More informationPDE Notes, Lecture #11
PDE Notes, Lecture # from Professor Jalal Shatah s Lectures Febuary 9th, 2009 Sobolev Spaces Recall that for u L loc we can efine the weak erivative Du by Du, φ := udφ φ C0 If v L loc such that Du, φ =
More informationHomework No. 1 MAE/CE 459/559 John A. Gilbert, Ph.D. Fall 2004
Homework No. 1 MAE/CE 459/559 John A. Gilbert, Ph.D. 1. A beam is loaded as shown. The dimensions of the cross section appear in the insert. the figure. Draw a complete free body diagram showing an equivalent
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 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 informationNDT&E Methods: UT. VJ Technologies CAVITY INSPECTION. Nondestructive Testing & Evaluation TPU Lecture Course 2015/16.
CAVITY INSPECTION NDT&E Methods: UT VJ Technologies NDT&E Methods: UT 6. NDT&E: Introduction to Methods 6.1. Ultrasonic Testing: Basics of Elasto-Dynamics 6.2. Principles of Measurement 6.3. The Pulse-Echo
More information'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21
Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting
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 information7.1 Support Vector Machine
67577 Intro. to Machine Learning Fall semester, 006/7 Lecture 7: Support Vector Machines an Kernel Functions II Lecturer: Amnon Shashua Scribe: Amnon Shashua 7. Support Vector Machine We return now to
More informationSymmetric Bending of Beams
Symmetric Bending of Beams beam is any long structural member on which loads act perpendicular to the longitudinal axis. Learning objectives Understand the theory, its limitations and its applications
More informationPartial Differential Equations
Chapter Partial Differential Equations. Introuction Have solve orinary ifferential equations, i.e. ones where there is one inepenent an one epenent variable. Only orinary ifferentiation is therefore involve.
More informationHydroelastic vibration of a rectangular perforated plate with a simply supported boundary condition
Fluid Structure Interaction and Moving Boundary Problems IV 63 Hydroelastic vibration of a rectangular perforated plate with a simply supported boundary condition K.-H. Jeong, G.-M. Lee, T.-W. Kim & J.-I.
More informationm (ft-lb/ft). Using the point-slope
ENGR 1990 Engineering athematics pplications of Derivatives E 560, E 570 Eample #1 Consier a long slener beam of length with a concentrate loa acting at istance a from the left en. Due to this loa, the
More informationG j dq i + G j. q i. = a jt. and
Lagrange Multipliers Wenesay, 8 September 011 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine
More informationNONLINEAR WAVE EQUATIONS ARISING IN MODELING OF SOME STRAIN-HARDENING STRUCTURES
NONLINEAR WAE EQUATIONS ARISING IN MODELING OF SOME STRAIN-HARDENING STRUCTURES DONGMING WEI Department of Mathematics, University of New Orleans, 2 Lakeshore Dr., New Orleans, LA 7148,USA E-mail: dwei@uno.edu
More information3.22 Mechanical Properties of Materials Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 3.22 Mechanical Properties of Materials Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Quiz #1 Example
More informationCalculus of variations - Lecture 11
Calculus of variations - Lecture 11 1 Introuction It is easiest to formulate the problem with a specific example. The classical problem of the brachistochrone (1696 Johann Bernoulli) is the search to fin
More informationVibro-acoustic response of FGM plates considering the thermal effects Tieliang Yang1, a, Qibai Huang1, *
3rd International Conference on Materials Engineering, Manufacturing Technology and Control (ICMEMTC 2016) Vibro-acoustic response of FGM plates considering the thermal effects Tieliang Yang1, a, Qibai
More informationLecture 1b. Differential operators and orthogonal coordinates. Partial derivatives. Divergence and divergence theorem. Gradient. A y. + A y y dy. 1b.
b. Partial erivatives Lecture b Differential operators an orthogonal coorinates Recall from our calculus courses that the erivative of a function can be efine as f ()=lim 0 or using the central ifference
More informationLecture 3: The Concept of Stress, Generalized Stresses and Equilibrium
Lecture 3: The Concept of Stress, Generalized Stresses and Equilibrium 3.1 Stress Tensor We start with the presentation of simple concepts in one and two dimensions before introducing a general concept
More information2. Rigid bar ABC supports a weight of W = 50 kn. Bar ABC is pinned at A and supported at B by rod (1). What is the axial force in rod (1)?
IDE 110 S08 Test 1 Name: 1. Determine the internal axial forces in segments (1), (2) and (3). (a) N 1 = kn (b) N 2 = kn (c) N 3 = kn 2. Rigid bar ABC supports a weight of W = 50 kn. Bar ABC is pinned at
More information1 Bending of beams Mindlin theory
1 BENDNG OF BEAMS MNDLN THEORY 1 1 Bending of beams Mindlin theory Cross-section kinematics assumptions Distributed load acts in the xz plane, which is also a plane of symmetry of a body Ω v(x = 0 m Vertical
More informationContinuum Mechanics Lecture 4 Fluid dynamics
Continuum Mechanics Lecture 4 Flui ynamics Prof. http://www.itp.uzh.ch/~teyssier Outline - Flui kinematics - Mass an momentum conservation laws - The energy equation - Real fluis - Ieal fluis - Incompressible
More informationPhysics 505 Electricity and Magnetism Fall 2003 Prof. G. Raithel. Problem Set 3. 2 (x x ) 2 + (y y ) 2 + (z + z ) 2
Physics 505 Electricity an Magnetism Fall 003 Prof. G. Raithel Problem Set 3 Problem.7 5 Points a): Green s function: Using cartesian coorinates x = (x, y, z), it is G(x, x ) = 1 (x x ) + (y y ) + (z z
More informationConservation Laws. Chapter Conservation of Energy
20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More informationProblem set 2: Solutions Math 207B, Winter 2016
Problem set : Solutions Math 07B, Winter 016 1. A particle of mass m with position x(t) at time t has potential energy V ( x) an kinetic energy T = 1 m x t. The action of the particle over times t t 1
More informationChain Rule. MATH 311, Calculus III. J. Robert Buchanan. Spring Department of Mathematics
3.33pt Chain Rule MATH 311, Calculus III J. Robert Buchanan Department of Mathematics Spring 2019 Single Variable Chain Rule Suppose y = g(x) and z = f (y) then dz dx = d (f (g(x))) dx = f (g(x))g (x)
More informationBasic Thermoelasticity
Basic hermoelasticity Biswajit Banerjee November 15, 2006 Contents 1 Governing Equations 1 1.1 Balance Laws.............................................. 2 1.2 he Clausius-Duhem Inequality....................................
More informationA Universal Model for Bingham Fluids with Two Characteristic Yield Stresses
A Universal Moel for Bingham Fluis with Two Characteristic Yiel Stresses NGDomostroeva 1 an NNTrunov DIMeneleyev Institute for Metrology Russia, StPetersburg 195 Moskovsky pr 19 February, 4, 9 Abstract:
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 informationMEC-E8001 FINITE ELEMENT ANALYSIS
MEC-E800 FINIE EEMEN ANAYSIS 07 - WHY FINIE EEMENS AND IS HEORY? Design of machines and structures: Solution to stress or displacement by analytical method is often impossible due to complex geometry,
More informationVUMAT for Fabric Reinforced Composites
VUMAT or Fabric Reinorce Composites. Introuction This ocument escribes a constitutive mo or abric reinorce composites that was introuce in Abaqus/Exicit 6.8. The mo has been imemente as a built-in VUMAT
More informationApplication of Laplace Iteration method to Study of Nonlinear Vibration of laminated composite plates
(3) 78 795 Application of Laplace Iteration method to Study of Nonlinear Vibration of laminated composite plates Abstract In this paper, free vibration characteristics of laminated composite plates are
More informationAnalysis of Infinitely Long Euler Bernoulli Beam on Two Parameter Elastic Foundation: Case of Point Load
Analysis of Infinitely Long Euler Bernoulli Beam on Two Parameter Elastic Foundation: Case of Point Load Dr. Mama, Benjamin Okwudili Department of Civil Engineering, University of Nigeria, Nsukka. benjamin.mama@unn.edu.ng
More informationThe Principle of Least Action
Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of
More informationTensors, Fields Pt. 1 and the Lie Bracket Pt. 1
Tensors, Fiels Pt. 1 an the Lie Bracket Pt. 1 PHYS 500 - Southern Illinois University September 8, 2016 PHYS 500 - Southern Illinois University Tensors, Fiels Pt. 1 an the Lie Bracket Pt. 1 September 8,
More informationDynamics of beams. Modes and waves. D. Clouteau. September 16, Department of Mechanical and Civil Engineering Ecole Centrale Paris, France
Dynamics of and waves Department of Mechanical and Civil Engineering Ecole Centrale Paris, France September 16, 2008 Outline 1 2 Pre-stressed Bernoulli-Euler in moving frames Kinematical assumption Resultant
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 informationName (Print) ME Mechanics of Materials Exam # 1 Date: October 5, 2016 Time: 8:00 10:00 PM
Name (Print) (Last) (First) Instructions: ME 323 - Mechanics of Materials Exam # 1 Date: October 5, 2016 Time: 8:00 10:00 PM Circle your lecturer s name and your class meeting time. Gonzalez Krousgrill
More informationPART I. Basic Concepts
PART I. Basic Concepts. Introduction. Basic Terminology of Structural Vibration.. Common Vibration Sources.. Forms of Vibration.3 Structural Vibration . Basic Terminology of Structural Vibration The term
More information