Suspension Dynamics Modeling for (TAMA and) LIGO and KAGRA. Mark Barton July GW Seminar 7/24/15
|
|
- Darren Dickerson
- 5 years ago
- Views:
Transcription
1 Suspension Dynamics Modeling for (TAMA and) LIGO and KAGRA Mark Barton July GW Seminar 7/24/15
2 Prehistory Started modeling the X- pendulum I developed for TAMA300 at ICRR Originally MathemaNca v2.2 (now v10.1)!
3 Glasgow Matlab Models Calum Torrie modeled the GEO triple as part of his PhD (LIGO- P v1). Lots of hand- coded Matlab hard to check and hard to modify. Also lots of simplificanons to the physics. Ken Strain extended Calum's code to model first aligo QUAD prototype. These were being used to design aligo suspensions. Needed something to compare these models against and allow for future improvements.
4 The MathemaNca Toolkit, PendUNl.nb The X-pendulum model was developed into a general toolkit. Implemented as a Mathematica package, PendUtil.nb, for specifying different configurations (e.g., quad, triple etc) in a (relatively) user-friendly way Supported features: 6-DOF rigid bodies for masses (no internal modes) Springs described by an elasticity tensor and a vector of pre-load forces Massless wires (i.e., no violin modes) but detailed elasticity model from beam equation Arbitrary frequency-dependent damping on all sources of elasticity Symbolic up to the point of minimizing the potential to find the equilibrium position Calculates elasticity and mass matrices semi-numerically (symbolic partial derivatives of functions with mostly numeric coefficients) Eigenfrequencies and eigenmodes calculated numerically Arbitrary frequency dependent damping on each different elastic element Transfer functions Thermal noise plots Export of state-space matrices to Matlab and E2E
5 Based on standard method of mass/snffness matrices Express the potennal energy of the system in terms of the coordinates: E P = E P ( x 1,...x n ) = E P ( x) Express the kinenc energy of the system in terms of the coordinates and velocines: E K = E K ( x 1,...x n,!x 1,...!x n ) Minimize the potennal energy to find the equilibrium values of the coordinates. ( ( ) ) T x eq = x 1( eq),...x n eq Compute the matrix of second derivanves of potennal, a.k.a., the potennal energy matrix or the snffness matrix. E P = E P ( x eq ) x x eq ( ) T K( x x eq ) K : K ij = E P x i x j x=x eq Compute the matrix of second derivanves of kinenc energy, a.k.a., the kinenc energy matrix or the mass matrix: E K = 1 2!xT M!x M : M ij = E K!x i!x j!x=0 x=x eq Write matrix equanons of monon and solve for parncular numerical frequencies using!x = iωx = 2πifx!!x = (2π f ) 2 x M!!x = Kx + Cx external + f external x = ( K ( 2π f ) 2 M) 1 ( Cx external + f external ) Or, assume a sinusoidal solunon with no external forces Ke i = ω i 2 Me i Do a simultaneous diagonalizanon to obtain the eigenfrequencies f i = ω i 2π and eigenmodes e i : x i ( t) = x eq + e i e ω it See Goldstein, "Classical Mechanics", 3 rd Ed.
6 Rigid- body mass coordinates Six coordinates for center of mass (COM): x, y, z, yaw (Y), pitch (P), roll (R) Three extra "body" coordinates for non- COM points such as wire agachments = x space y space z space x COM y COM z COM = + x COM y COM z COM + cosy siny 0 siny cosy cos P 0 sin P sin P 0 cos P The model definer needs to provide a list of the COM coordinates of all the masses cos R sin R 0 sin R cos R cos P cosy cosy sin Psin R cos RsinY cos RcosY sin P + sin RsinY cos PsinY cos RcosY + sin Psin RsinY cos Rsin PsinY cosy sin R sin P cos Psin R cos P cos R x body y body z body x body y body z body
7 KineNc energy There is typically a linear term in terms of mass m plus angular term in terms of moment of inerna tensor I: ( ) E K = 1 2 m(!x 2 +!y 2 +!z 2 ) ω x ω y ω z I xx I xy I zx I xy I yy I yz I xz I yz I zz ω x ω y ω z I is in body coordinates so the angular velocity vector needs to be transformed: ω x ω y ω z = sin P 0 1 cos Psin R cos R 0 cos P cos R sin R 0 Y!!P!R The model definer needs to supply a MathemaNca expression for the kinenc energy for all masses, with terms similar to the above.
8 GravitaNonal energy E G = mgz The model definer has to supply a MathemaNca expression giving the total gravitanonal potennal for all masses.
9 Wire energy Wire energy is broken into four terms: Longitudinal stretch based on straight- line distance Bending near the endpoints Extra longitudinal stretch due to bending TJ Torsion E P(wire_torsional ) = 1 2 A + GJ 2 e Δθ T Model definer needs to supply a list of parameters for each wire, including Young's modulus, length, moment of area, agachment points, etc, etc. E P( wire_bending) = YM 1 2 E P(wire_longitudinal ) = 1 2 k w E P wire_ extra _ stretch ( ) = T L d 2 y dl 2 ( l( t) l 0 ) 2 dy dl L 2 dl 2 dl + YM L d 2 z dl 2 2 dl
10 Flexure correcnon Wire bending terms use some complicated 3D geometry and are slow to calculate. It turns out there's a shortcut: the wire behaves as if agached with a hinge at a distance a from actual break- off point, where a = EI T E=Young's modulus, I = second moment of area, T = tension This trick was not used in the MathemaNca toolkit but has been used in the exported Matlab code.
11 Springs Simple zero- length spring model connecnng a point on one mass to a point on another. 6x6 matrix of elasnc constants plus a 6x1 vector of pre- load forces: E P(spring) = 1 ( 2 Δx Δy Δz ΔY ΔP ΔR ) K xx K xy K xz K xy K xp K xr K xy K yy K yz K yy K yp K yr K xz K yz K zz K zy K zp K zr K xy K yy K zy K YY K YP K YR K xp K yp K zp K YP K PP K PR K xr K yr K zr K YR K PR K RR Model definer has to provide a list of springs with masses, agachment points and elasnc constant and pre- load informanon. Δx Δy Δz ΔY ΔP ΔR + f x f y f z f Y f P f R ( ) Δx Δy Δz ΔY ΔP ΔR
12 As well as the coordinates used in the normal mode analysis ("vars"), two other sorts are important: "params" the support and other objects which are stanonary during normal modes but move for transfer funcnons "floats" things like wire- spring juncnons with no associated mass Extra coordinates
13 SNffness matrix for extra coordinates It is convenient to calculate a master snffness matrix with parnal derivanves between all types of coordinates: "vars" "floats" K C XQ C XS K master = C QX Q C QS C SX C SQ S K, Q and S give the coupling among vars, floats and params within their respecnve groups. C XQ, C QS and C SX give the coupling between groups. Because the "floats" are dependent on the others, they need to be eliminated: "params" "vars" "floats" "params" K effective = K C XQ Q 1 C QX C SX(effective) = C SX C SQ Q 1 C QX S effective = S C SQ Q 1 C QS
14 Damping Lossiness in elasnc components can be represented by a complex elasnc constant: k k 0 ε ω ( ( ) + i ε ( ω )) The real and imaginary parts should sansfy the Krämers Krönig relanonship: ( ) 1 = 2 π PV ε ( x) ε ω dx ε ω x ω If losses are small, the real part can be assumed to be constant: k k 0 ( 1+ iφ ( f )) The model definer can specify a different damping funcnon for each elasnc component. The toolkit keeps track of which damping funcnon applies to which coefficients in the master snffness matrix. ( ) = 2 π PV ε ( x) 1 x ω dx
15 Wire/Fibre damping Wire φ usually has a frequency- independent "structural" term and a frequency dependent "thermoelasnc " term f f HHzL ThermoelasNc damping has a characterisnc peak at the Nmescale of heat flow across the wire. HLTS HSTS
16 Thermal Noise Suspension thermal noise is a potennal liminng factor in GW detectors. Noise is given in terms of damping by FluctuaNon DissipaNon Theorem: x 2 ( ( )) ( ω ) = 4k BT Re Y ω Y is the admigance (v/f) at a test point where the thermal noise is to be calculated. ω 2
17 DissipaNon DiluNon (i) In a simple mass- spring system, the quality factor of the oscillanon depends purely on the spring material: 1 Q = φ material However systems such as pendulums and wires are used in GW detectors because with certain geometries, the damping (dissipanon of energy) is "diluted": Q = D φ material 1 φ material
18 DissipaNon DiluNon (ii) The dissipanon dilunon factor D depends on the fracnon of the energy involved in first- order stress changes of the material. Pulling a pendulum or violin string sideways creates only a second- order or smaller stretch - > dilunon. Quick test: calculate the contribunon to the potennal matrix for each potennal term twice, once with and without all tensions zeroed.
19 CalculaNon Procedure (i) The model calculanon notebook is run. It loads the model defininon notebook, which loads the toolkit, the model defininon, and a default set of numerical values. The calculanon notebook can then selecnvely override some of the numerical values. The wire and spring lists are processed to create a list of potennal terms, each with a damping funcnon. The total potennal is computed, with numerical values for all quannnes except the "var" coordinates. (Usually the wire bending potennal terms are omiged for speed.) The total potennal is minimized to find the equilibrium posinon.
20 CalculaNon Procedure (ii) Numerical values for everything but the coordinates are subsntuted into the potennal terms. The potennal terms are differennated to find the snffness matrix elements. Each term is processed separately for two reasons: for speed (each term depends on only a few coordinates, so most derivanves are zero and don't need to be computed) to keep contribunons with different damping funcnons separate. The process is repeated with the tension switched off, to allow dissipanon dilunon to be calculated.
21 CalculaNon Procedure (iii) Versions of the snffness matrix with and without damping are calculated. The snffness matrix without damping (a totally numerical matrix) is used to calculate the eigenfrequencies and eigenmodes. The snffness matrix with damping (a mostly numerical matrix that is a MathemaNca funcnon of frequency, f) is used for all other purposes. MathemaNca funcnons are provided to allow transfer funcnons from/to selected inputs/outputs, thermal noise plots at selected test points, eigenmode shape plots etc. If the small but Nme- consuming wire bending/torsion potennal terms were omiged at the beginning, they are computed and added in.
22 Models for LIGO Models defined for all the LIGO suspensions (no KAGRA yet): QUAD: single chain of AdvLIGO quad pendulum, with 4 masses, 6 blade springs and 14 wires. BSFM, HSTS, HLTS: 3 masses, 6 blade springs and 10 wires. OFMC, TMTS: 2 masses, 2 or 4 blades, 6 wires HAUX, HTTS, OFIS: 1 mass, assorted blade/wire combinations Many toy models
23 2 blade springs 2 wires top mass 2 blade springs 4 wires upper intermediate mass 2 blade springs 4 wires intermediate mass 4 fibres (or wires) optic (or reaction mass) Example model: LIGO quad pendulum
24 Sample Output (i) Table of Mode Freqs/Shapes Model " TMproducNonTM" The aligo quad suspension main chain, with monolithic final stage (i.e., fused silica fibres suppornng the test mass)
25 Sample Output (ii) Individual Mode Shape
26 Sample Output (iii) Transfer Function
27 Sample Output (iv) Thermal Noise
28 State Space Terminology For simulanons, it is convenient to write the equanons of monon in state- space format. In tradinonal notanon:!x y = A B C D x: vector of state coordinates u: vector of inputs y: vector of outputs A: state matrix (maps state to rate of change) B: input matrix (maps u to rate of change) C: output matrix (maps state to output) D: feedthrough matrix (maps inputs directly to output) x u
29 SS State/Inputs/Outputs for Suspension Models In mechanical simulanons, the EOMs are second- order, so the full state needs to be the coordinates ("vars") plus the velocines: x x!x A good set of inputs is the support coordinates ("params") plus input forces on the main coordinates: u = s!s f x A good set of outputs is all the main coordinates plus the reacnon forces on the support: y = x f s!x!!x = A x!x + Bu
30 Full State Space for Suspension Models!x!!x 2 N x 1 = 0 N x N x I N x N x M 1 K N x N x M 1 E A N x N x 2 N x 2 N x x!x 2 N x N x N s 0 N x N s 0 N x N x M 1 C N x N s M 1 E B N x N s M 1 N x N x s!s f x ( 2 N s +N x ) 1 x f s ( N x +N s ) 1 = I N x N x C N s N x 0 N x N x 0 N s N x ( N x +N s ) 2 N x x!x 2 N x N x N s 0 N x N s 0 N x N x S N s N s E D N x N s 0 N s N x ( N x +N s ) ( N s +N x ) s!s f x ( 2 N s +N x ) 1 K is the snffness matrix; M is the mass matrix C is the coupling snffness matrix; S is the support snffness The E matrices are arrays of velocity damping coefficients. Sizes of matrices and vectors are given as <Nrows x Ncolumns>.
31 Matlab Export All the components of the SS A/B/C/D matrices can easily be exported from MathemaNca as numbers or Matlab code. For symbolic export, it's beger to export M and calculate M - 1 in Matlab.
32 Usage at LIGO MathemaNca models were used to study thermal noise and asymmetrical suspensions. The GEO Matlab models were retained but had the hand- wrigen state- space matrices replaced with improved code generated by the corresponding MathemaNca model. Matlab models were used for analysis of experimental data and for control design. All MathemaNca models are hosted on the SUS SVN (Subversion version control system) at hgps://redoubt.ligo- wa.caltech.edu/websvn/lisnng.php? repname=sus&path=%2ftrunk%2fcommon %2FMathemaNcaModels %2F#path_trunk_Common_MathemaNcaModels_ (needs LIGO.org credennals). Matlab models are in the same SVN but scagered about.
33 References LIGO- T Models of the Advanced LIGO Suspensions in MathemaNca LIGO- T Models of the Advanced LIGO Suspensions in MATLAB hgps://awiki.ligo- wa.caltech.edu/aligo/ Suspensions/OpsManual (needs LIGO credennals or see LIGO- E ). LIGO- T DissipaNon dilunon LIGO- T Flexure correcnons
34 KAGRA Sekiguchi- san has used some of the ideas and wrigen his own package. Much more automated and easy to use. Probably no need (certainly no plan) to redo KAGRA models with LIGO toolkit. Not too hard to do if there is a need for say detailed thermal noise plots.
Pendulum Modeling in Mathematica and Matlab
Pendulum Modeling in Mathematica and Matlab IGR Thermal Noise Group Meeting 4 May 2007 4 May 2007 The X-Pendulum Developed as a low frequency vibration isolator for TAMA 1D version 2D version 4 May 2007
More information3D Models of Vibration Isolation Systems in Mathematica
3D Models of Vibration Isolation Systems in Mathematica April 12 th,2011 Takanori Sekiguchi Concept Three-dimensional rigid-body models in Mathematica Inspired by Mark Barton s suspension models http://www.ligo.caltech.edu/~e2e/susmodels/
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 informationResonances, quality factors in LIGO pendulums: a compendium of formulas
Resonances, quality factors in LIGO pendulums: a compendium of formulas Gabriela González December 11, 2001 1 Introduction and notes There are simple formulas for the different pendulum and violin modes,
More informationDissipation Dilution
LASER INERFEROMEER GRAVIAIONAL WAVE OBSERVAORY LIGO Laboratory / LIGO Scientific Collaboration LIGO-711-draft75 ADVANCED LIGO May 7 Dissipation Dilution Mark Barton Distribution of this document: LIGO
More informationA DARK GREY P O N T, with a Switch Tail, and a small Star on the Forehead. Any
Y Y Y X X «/ YY Y Y ««Y x ) & \ & & } # Y \#$& / Y Y X» \\ / X X X x & Y Y X «q «z \x» = q Y # % \ & [ & Z \ & { + % ) / / «q zy» / & / / / & x x X / % % ) Y x X Y $ Z % Y Y x x } / % «] «] # z» & Y X»
More informationClosed-Form Solution Of Absolute Orientation Using Unit Quaternions
Closed-Form Solution Of Absolute Orientation Using Unit Berthold K. P. Horn Department of Computer and Information Sciences November 11, 2004 Outline 1 Introduction 2 3 The Problem Given: two sets of corresponding
More informationBASELINE SUSPENSION DESIGN FOR LIGO II - UPDATE
BASELINE SUSPENSION DESIGN FOR LIGO II - UPDATE Norna A Robertson University of Glasgow for the GEO suspension team LSC, Hanford 15th August 2000 LIGO-G000295-00-D GEO suspension team for LIGO II G Cagnoli,
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 informationAB-267 DYNAMICS & CONTROL OF FLEXIBLE AIRCRAFT
FLÁIO SILESTRE DYNAMICS & CONTROL OF FLEXIBLE AIRCRAFT LECTURE NOTES LAGRANGIAN MECHANICS APPLIED TO RIGID-BODY DYNAMICS IMAGE CREDITS: BOEING FLÁIO SILESTRE Introduction Lagrangian Mechanics shall be
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 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 informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Overview Milestones in continuum mechanics Concepts of modulus and stiffness. Stress-strain relations Elasticity Surface and body
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 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 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 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 informationLES of Turbulent Flows: Lecture 11 (ME EN )
1 LES of Turbulent Flows: Lecture 11 (ME EN 7960-003) Prof. Rob Stoll Department of Mechanical Engineering University of Utah Fall 2014 2 1- EquaNon Eddy viscosity Models EvoluNon of τ ij : Deardorff (ASME
More information12. Stresses and Strains
12. Stresses and Strains Finite Element Method Differential Equation Weak Formulation Approximating Functions Weighted Residuals FEM - Formulation Classification of Problems Scalar Vector 1-D T(x) u(x)
More informationAircraft Structural Analysis Summary
Aircraft Structural Analysis Summary 1. Basic quations 1.1 Definitions, conventions and basic relations Before we can start throwing around equations, we have to define some variables and state some conventions.
More informationA. H. Hall, 33, 35 &37, Lendoi
7 X x > - z Z - ----»»x - % x x» [> Q - ) < % - - 7»- -Q 9 Q # 5 - z -> Q x > z»- ~» - x " < z Q q»» > X»? Q ~ - - % % < - < - - 7 - x -X - -- 6 97 9
More informationApplications of Eigenvalues & Eigenvectors
Applications of Eigenvalues & Eigenvectors Louie L. Yaw Walla Walla University Engineering Department For Linear Algebra Class November 17, 214 Outline 1 The eigenvalue/eigenvector problem 2 Principal
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 informationM5 Simple Beam Theory (continued)
M5 Simple Beam Theory (continued) Reading: Crandall, Dahl and Lardner 7.-7.6 In the previous lecture we had reached the point of obtaining 5 equations, 5 unknowns by application of equations of elasticity
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 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 informationStructural Dynamics Lecture 4. Outline of Lecture 4. Multi-Degree-of-Freedom Systems. Formulation of Equations of Motions. Undamped Eigenvibrations.
Outline of Multi-Degree-of-Freedom Systems Formulation of Equations of Motions. Newton s 2 nd Law Applied to Free Masses. D Alembert s Principle. Basic Equations of Motion for Forced Vibrations of Linear
More information7.3 Design of members subjected to combined forces
7.3 Design of members subjected to combined forces 7.3.1 General In the previous chapters of Draft IS: 800 LSM version, we have stipulated the codal provisions for determining the stress distribution in
More informationEnergy Considerations
Physics 42200 Waves & Oscillations Lecture 4 French, Chapter 3 Spring 2016 Semester Matthew Jones Energy Considerations The force in Hooke s law is = Potential energy can be used to describe conservative
More informationIntroduction, Basic Mechanics 2
Computational Biomechanics 18 Lecture : Introduction, Basic Mechanics Ulli Simon, Lucas Engelhardt, Martin Pietsch Scientific Computing Centre Ulm, UZWR Ulm University Contents Mechanical Basics Moment
More informationSec. 14.3: Partial Derivatives. All of the following are ways of representing the derivative. y dx
Math 2204 Multivariable Calc Chapter 14: Partial Derivatives I. Review from math 1225 A. First Derivative Sec. 14.3: Partial Derivatives 1. Def n : The derivative of the function f with respect to the
More informationII&Ij <Md Tmlaiiiiiit, aad once in Ihe y a w Teataa m i, the vmb thatalmta oot Uiaapirit world. into as abode or wotld by them- CooTBOtioa
382 4 7 q X
More informationChapter 5 Elastic Strain, Deflection, and Stability 1. Elastic Stress-Strain Relationship
Chapter 5 Elastic Strain, Deflection, and Stability Elastic Stress-Strain Relationship A stress in the x-direction causes a strain in the x-direction by σ x also causes a strain in the y-direction & z-direction
More informationLOWELL WEEKLY JOURNAL
Y -» $ 5 Y 7 Y Y -Y- Q x Q» 75»»/ q } # ]»\ - - $ { Q» / X x»»- 3 q $ 9 ) Y q - 5 5 3 3 3 7 Q q - - Q _»»/Q Y - 9 - - - )- [ X 7» -» - )»? / /? Q Y»» # X Q» - -?» Q ) Q \ Q - - - 3? 7» -? #»»» 7 - / Q
More informationRotational Motion. Chapter 4. P. J. Grandinetti. Sep. 1, Chem P. J. Grandinetti (Chem. 4300) Rotational Motion Sep.
Rotational Motion Chapter 4 P. J. Grandinetti Chem. 4300 Sep. 1, 2017 P. J. Grandinetti (Chem. 4300) Rotational Motion Sep. 1, 2017 1 / 76 Angular Momentum The angular momentum of a particle with respect
More informationLecture 8 Analyzing the diffusion weighted signal. Room CSB 272 this week! Please install AFNI
Lecture 8 Analyzing the diffusion weighted signal Room CSB 272 this week! Please install AFNI http://afni.nimh.nih.gov/afni/ Next lecture, DTI For this lecture, think in terms of a single voxel We re still
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 information2.003 Engineering Dynamics Problem Set 6 with solution
.00 Engineering Dynamics Problem Set 6 with solution Problem : A slender uniform rod of mass m is attached to a cart of mass m at a frictionless pivot located at point A. The cart is connected to a fixed
More information1. What would be the value of F1 to balance the system if F2=20N? 20cm T =? 20kg
1. What would be the value of F1 to balance the system if F2=20N? F2 5cm 20cm F1 (a) 3 N (b) 5 N (c) 4N (d) None of the above 2. The stress in a wire of diameter 2 mm, if a load of 100 gram is applied
More informationWave and Elasticity Equations
1 Wave and lasticity quations Now let us consider the vibrating string problem which is modeled by the one-dimensional wave equation. Suppose that a taut string is suspended by its extremes at the points
More informationEQUIVALENT SINGLE-DEGREE-OF-FREEDOM SYSTEM AND FREE VIBRATION
1 EQUIVALENT SINGLE-DEGREE-OF-FREEDOM SYSTEM AND FREE VIBRATION The course on Mechanical Vibration is an important part of the Mechanical Engineering undergraduate curriculum. It is necessary for the development
More informationNon-Linear Response of Test Mass to External Forces and Arbitrary Motion of Suspension Point
LASER INTERFEROMETER GRAVITATIONAL WAVE OBSERVATORY -LIGO- CALIFORNIA INSTITUTE OF TECHNOLOGY MASSACHUSETTS INSTITUTE OF TECHNOLOGY Technical Note LIGO-T980005-01- D 10/28/97 Non-Linear Response of Test
More informationTwo Posts to Fill On School Board
Y Y 9 86 4 4 qz 86 x : ( ) z 7 854 Y x 4 z z x x 4 87 88 Y 5 x q x 8 Y 8 x x : 6 ; : 5 x ; 4 ( z ; ( ) ) x ; z 94 ; x 3 3 3 5 94 ; ; ; ; 3 x : 5 89 q ; ; x ; x ; ; x : ; ; ; ; ; ; 87 47% : () : / : 83
More informationLecture 3 The Concept of Stress, Generalized Stresses and Equilibrium
Lecture 3 The Concept of Stress, Generalized Stresses and Equilibrium Problem 3-1: Cauchy s Stress Theorem Cauchy s stress theorem states that in a stress tensor field there is a traction vector t that
More informationHIGHER-ORDER THEORIES
HIGHER-ORDER THEORIES THIRD-ORDER SHEAR DEFORMATION PLATE THEORY LAYERWISE LAMINATE THEORY J.N. Reddy 1 Third-Order Shear Deformation Plate Theory Assumed Displacement Field µ u(x y z t) u 0 (x y t) +
More informationStress, Strain, Mohr s Circle
Stress, Strain, Mohr s Circle The fundamental quantities in solid mechanics are stresses and strains. In accordance with the continuum mechanics assumption, the molecular structure of materials is neglected
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
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 information' Liberty and Umou Ono and Inseparablo "
3 5? #< q 8 2 / / ) 9 ) 2 ) > < _ / ] > ) 2 ) ) 5 > x > [ < > < ) > _ ] ]? <
More informationConsider a particle in 1D at position x(t), subject to a force F (x), so that mẍ = F (x). Define the kinetic energy to be.
Chapter 4 Energy and Stability 4.1 Energy in 1D Consider a particle in 1D at position x(t), subject to a force F (x), so that mẍ = F (x). Define the kinetic energy to be T = 1 2 mẋ2 and the potential energy
More informationContinuum mechanism: Stress and strain
Continuum mechanics deals with the relation between forces (stress, σ) and deformation (strain, ε), or deformation rate (strain rate, ε). Solid materials, rigid, usually deform elastically, that is the
More informationMechanics of Solids notes
Mechanics of Solids notes 1 UNIT II Pure Bending Loading restrictions: As we are aware of the fact internal reactions developed on any cross-section of a beam may consists of a resultant normal force,
More informationGG612 Lecture 3. Outline
GG61 Lecture 3 Strain and Stress Should complete infinitesimal strain by adding rota>on. Outline Matrix Opera+ons Strain 1 General concepts Homogeneous strain 3 Matrix representa>ons 4 Squares of line
More informationMass on a spring: including a driving force, and friction proportional to velocity (eg, from a dashpot).
Physics 06a, Caltech 8 October, 208 Lecture 6: Equilibria and Oscillations We study equilibria and linear oscillations. Motion of dynamical systems can be classified as secular (long-term, non-periodic)
More informationLecture 20: Modes in a crystal and continuum
Physics 16a, Caltech 11 December, 218 The material in this lecture (which we had no time for) will NOT be on the Final exam. Lecture 2: Modes in a crystal and continuum The vibrational modes of the periodic
More informationMemoir on the analytical expression of the elasticity and stiffness of curves of double curvature
Mémoire sur l expression analytique de élasticité et de la raideur des courbes à double courbure J Ecole poly Cahier XVII 10 (1815) 418-456 Memoir on the analytical expression of the elasticity and stiffness
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 informationPart 8: Rigid Body Dynamics
Document that contains homework problems. Comment out the solutions when printing off for students. Part 8: Rigid Body Dynamics Problem 1. Inertia review Find the moment of inertia for a thin uniform rod
More informationSTATICALLY INDETERMINATE STRUCTURES
STATICALLY INDETERMINATE STRUCTURES INTRODUCTION Generally the trusses are supported on (i) a hinged support and (ii) a roller support. The reaction components of a hinged support are two (in horizontal
More informationChapter 1 Fluid Characteristics
Chapter 1 Fluid Characteristics 1.1 Introduction 1.1.1 Phases Solid increasing increasing spacing and intermolecular liquid latitude of cohesive Fluid gas (vapor) molecular force plasma motion 1.1.2 Fluidity
More informationEDEXCEL NATIONAL CERTIFICATE/DIPLOMA MECHANICAL PRINCIPLES AND APPLICATIONS NQF LEVEL 3 OUTCOME 1 - LOADING SYSTEMS TUTORIAL 3 LOADED COMPONENTS
EDEXCEL NATIONAL CERTIICATE/DIPLOMA MECHANICAL PRINCIPLES AND APPLICATIONS NQ LEVEL 3 OUTCOME 1 - LOADING SYSTEMS TUTORIAL 3 LOADED COMPONENTS 1. Be able to determine the effects of loading in static engineering
More informationENG2000 Chapter 7 Beams. ENG2000: R.I. Hornsey Beam: 1
ENG2000 Chapter 7 Beams ENG2000: R.I. Hornsey Beam: 1 Overview In this chapter, we consider the stresses and moments present in loaded beams shear stress and bending moment diagrams We will also look at
More informationPURE BENDING. If a simply supported beam carries two point loads of 10 kn as shown in the following figure, pure bending occurs at segment BC.
BENDING STRESS The effect of a bending moment applied to a cross-section of a beam is to induce a state of stress across that section. These stresses are known as bending stresses and they act normally
More informationName :. Roll No. :... Invigilator s Signature :.. CS/B.TECH (CE-NEW)/SEM-3/CE-301/ SOLID MECHANICS
Name :. Roll No. :..... Invigilator s Signature :.. 2011 SOLID MECHANICS Time Allotted : 3 Hours Full Marks : 70 The figures in the margin indicate full marks. Candidates are required to give their answers
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 information3. The linear 3-D elasticity mathematical model
3. The linear 3-D elasticity mathematical model In Chapter we examined some fundamental conditions that should be satisfied in the modeling of all deformable solids and structures. The study of truss structures
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 informationkx m x B N 1 C L, M Mg θ
.004 MODELING DYNAMICS AND CONTROL II Spring 00 Solutions to Problem Set No. 7 Problem 1. Pendulum mounted on elastic support. This problem is an execise in the application of momentum principles. Two
More informationCoupled Oscillators. 1 Introduction. 2 Theory. PHY 300 Lab 2 Fall 2012
Coupled Oscillators 1 Introduction In this experiment you are going to observe the normal modes of oscillation of several different mechanical systems, first on the air tracks and then using some coupled
More informationChapter 3. Load and Stress Analysis
Chapter 3 Load and Stress Analysis 2 Shear Force and Bending Moments in Beams Internal shear force V & bending moment M must ensure equilibrium Fig. 3 2 Sign Conventions for Bending and Shear Fig. 3 3
More informationFundamentals Physics. Chapter 15 Oscillations
Fundamentals Physics Tenth Edition Halliday Chapter 15 Oscillations 15-1 Simple Harmonic Motion (1 of 20) Learning Objectives 15.01 Distinguish simple harmonic motion from other types of periodic motion.
More informationPhysics 351, Spring 2017, Homework #4. Due at start of class, Friday, February 10, 2017
Physics 351, Spring 2017, Homework #4. Due at start of class, Friday, February 10, 2017 Course info is at positron.hep.upenn.edu/p351 When you finish this homework, remember to visit the feedback page
More informationLANMARK UNIVERSITY OMU-ARAN, KWARA STATE DEPARTMENT OF MECHANICAL ENGINEERING COURSE: MECHANICS OF MACHINE (MCE 322). LECTURER: ENGR.
LANMARK UNIVERSITY OMU-ARAN, KWARA STATE DEPARTMENT OF MECHANICAL ENGINEERING COURSE: MECHANICS OF MACHINE (MCE 322). LECTURER: ENGR. IBIKUNLE ROTIMI ADEDAYO SIMPLE HARMONIC MOTION. Introduction Consider
More information2 Introduction to mechanics
21 Motivation Thermodynamic bodies are being characterized by two competing opposite phenomena, energy and entropy which some researchers in thermodynamics would classify as cause and chance or determinancy
More informationECEN 420 LINEAR CONTROL SYSTEMS. Lecture 6 Mathematical Representation of Physical Systems II 1/67
1/67 ECEN 420 LINEAR CONTROL SYSTEMS Lecture 6 Mathematical Representation of Physical Systems II State Variable Models for Dynamic Systems u 1 u 2 u ṙ. Internal Variables x 1, x 2 x n y 1 y 2. y m Figure
More informationOSEM Coil/Magnet/Flag Calculation
KAGRA 10 September 2016 OSEM Coil/Magnet/Flag Calculation Mark Barton Distribution of this document: JGW-DCC This is an internal working note of the KAGRA collaboration. http://gwcenter.icrr.u-tokyo.ac.jp
More informationMECH 5312 Solid Mechanics II. Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso
MECH 5312 Solid Mechanics II Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso Table of Contents Preliminary Math Concept of Stress Stress Components Equilibrium
More informationDevelopment of ground based laser interferometers for the detection of gravitational waves
Development of ground based laser interferometers for the detection of gravitational waves Rahul Kumar ICRR, The University of Tokyo, 7 th March 2014 1 Outline 1. Gravitational waves, nature & their sources
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 informationShafts: Torsion of Circular Shafts Reading: Crandall, Dahl and Lardner 6.2, 6.3
M9 Shafts: Torsion of Circular Shafts Reading: Crandall, Dahl and Lardner 6., 6.3 A shaft is a structural member which is long and slender and subject to a torque (moment) acting about its long axis. We
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 informationRock Rheology GEOL 5700 Physics and Chemistry of the Solid Earth
Rock Rheology GEOL 5700 Physics and Chemistry of the Solid Earth References: Turcotte and Schubert, Geodynamics, Sections 2.1,-2.4, 2.7, 3.1-3.8, 6.1, 6.2, 6.8, 7.1-7.4. Jaeger and Cook, Fundamentals of
More information6. Bending CHAPTER OBJECTIVES
CHAPTER OBJECTIVES Determine stress in members caused by bending Discuss how to establish shear and moment diagrams for a beam or shaft Determine largest shear and moment in a member, and specify where
More informationUnit - 7 Vibration of Continuous System
Unit - 7 Vibration of Continuous System Dr. T. Jagadish. Professor for Post Graduation, Department of Mechanical Engineering, Bangalore Institute of Technology, Bangalore Continuous systems are tore which
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 informationInstitute of Structural Engineering Page 1. Method of Finite Elements I. Chapter 2. The Direct Stiffness Method. Method of Finite Elements I
Institute of Structural Engineering Page 1 Chapter 2 The Direct Stiffness Method Institute of Structural Engineering Page 2 Direct Stiffness Method (DSM) Computational method for structural analysis Matrix
More informationOWELL WEEKLY JOURNAL
Y \»< - } Y Y Y & #»»» q ] q»»»>) & - - - } ) x ( - { Y» & ( x - (» & )< - Y X - & Q Q» 3 - x Q Y 6 \Y > Y Y X 3 3-9 33 x - - / - -»- --
More informationMAE 323: Lecture 1. Review
This review is divided into two parts. The first part is a mini-review of statics and solid mechanics. The second part is a review of matrix/vector fundamentals. The first part is given as an refresher
More informationRobot Control Basics CS 685
Robot Control Basics CS 685 Control basics Use some concepts from control theory to understand and learn how to control robots Control Theory general field studies control and understanding of behavior
More informationChapter 12. Static Equilibrium and Elasticity
Chapter 12 Static Equilibrium and Elasticity Static Equilibrium Equilibrium implies that the object moves with both constant velocity and constant angular velocity relative to an observer in an inertial
More informationChapter 2 Governing Equations
Chapter Governing Equations Abstract In this chapter fundamental governing equations for propagation of a harmonic disturbance on the surface of an elastic half-space is presented. The elastic media is
More informationResponse Spectrum Analysis Shock and Seismic. FEMAP & NX Nastran
Response Spectrum Analysis Shock and Seismic FEMAP & NX Nastran Table of Contents 1. INTRODUCTION... 3 2. THE ACCELEROGRAM... 4 3. CREATING A RESPONSE SPECTRUM... 5 4. NX NASTRAN METHOD... 8 5. RESPONSE
More informationHIGHER-ORDER THEORIES
HIGHER-ORDER THEORIES Third-order Shear Deformation Plate Theory Displacement and strain fields Equations of motion Navier s solution for bending Layerwise Laminate Theory Interlaminar stress and strain
More informationDROP-WEIGHT SYSTEM FOR DYNAMIC PRESSURE CALIBRATION
DROP-WIGHT SYSTM FOR DYNAMIC PRSSUR CALIBRATION Antti Lakka and Sari Saxholm MIKS, Centre for Metrology and Accreditation P.O. Box 9 (Tekniikantie 1) FI-2151 spoo Finland telephone: +358-29-554456, e-mail:
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 informationBasic Energy Principles in Stiffness Analysis
Basic Energy Principles in Stiffness Analysis Stress-Strain Relations The application of any theory requires knowledge of the physical properties of the material(s) comprising the structure. We are limiting
More informationToward a novel approach for damage identification and health monitoring of bridge structures
Toward a novel approach for damage identification and health monitoring of bridge structures Paolo Martino Calvi 1, Paolo Venini 1 1 Department of Structural Mechanics, University of Pavia, Italy E-mail:
More informationProfessor Terje Haukaas University of British Columbia, Vancouver Notation
Notation This document establishes the notation that is employed throughout these notes. It is intended as a look-up source during the study of other documents and software on this website. As a general
More informationPhysics 106a, Caltech 16 October, Lecture 5: Hamilton s Principle with Constraints. Examples
Physics 106a, Caltech 16 October, 2018 Lecture 5: Hamilton s Principle with Constraints We have been avoiding forces of constraint, because in many cases they are uninteresting, and the constraints can
More informationName: Fall 2014 CLOSED BOOK
Name: Fall 2014 1. Rod AB with weight W = 40 lb is pinned at A to a vertical axle which rotates with constant angular velocity ω =15 rad/s. The rod position is maintained by a horizontal wire BC. Determine
More information