BAR & TRUSS FINITE ELEMENT. Direct Stiffness Method
|
|
- Caroline Austin
- 5 years ago
- Views:
Transcription
1 BAR & TRUSS FINITE ELEMENT Drect Stness Method FINITE ELEMENT ANALYSIS AND APPLICATIONS INTRODUCTION TO FINITE ELEMENT METHOD What s the nte element method (FEM)? A technqe or obtanng approxmate soltons o derental eqatons. Partton o the doman nto a set o smple shapes (element) Approxmate the solton sng pecewse polynomals wthn the element F xx xy bx 0 x y xy yy Strctre by 0 x y Element Pecewse-Lnear Approxmaton x
2 INTRODUCTION TO FEM cont. How to dscretze the doman? Usng smple shapes (element) All elements are connected sng nodes Nodes Elements Solton at Element s descrbed sng the vales at Nodes,, 6, and 5 (Interpolaton). Elements and share the solton at Nodes and 6. INTRODUCTION TO FEM cont. Methods Drect method: Easy to nderstand, lmted to D problems Varatonal method Weghted resdal method Obectves Determne dsplacements, orces, and spportng reactons Wll consder only statc problem
3 -D SYSTEM OF SPRINGS F 5 F 6 5 F Bodes move only n horzontal drecton External orces, F, F, and F, are appled No need to dscretze the system (t s already dscretzed!) Rgd body (ncldng walls) NODE Sprng ELEMENT 5 SPRING ELEMENT Element e Consst o Nodes and, Sprng constant k (e) Force appled to the nodes: e e, () e e, () e Dsplacement and Elongaton: e e e e Force n the sprng: P k k Relaton b/w sprng orce and nodal orces: Eqlbrm: ( e ) e e e e 0 or e P e 6
4 Sprng Element e SPRING ELEMENT cont. Relaton between nodal orces and dsplacements e e e e k k Matrx notaton: k: stness matrx q: vector o DOFs : vector o element orces e e k e e k k k k [ k ] q k q 7 Stness matrx SPRING ELEMENT cont. It s sqare as t relates to the same nmber o orces as the dsplacements. It s symmetrc. It s snglar,.e., determnant s eqal to zero and t cannot be nverted. It s postve sem-dente Observaton For gven nodal dsplacements, nodal orces can be calclated by [ k ] q ( e ) ( e ) ( e ) For gven nodal orces, nodal dsplacements cannot be determned nqely e k k e k k e e 8
5 Element eqaton and assembly SYSTEM OF SPRINGS cont. F F 5 F 6 5 () k k () k k () k k () k k () k k () k k () k k () () k 0 k 0 k k () 0 k 0 k SYSTEM OF SPRINGS cont. () k k () () () () k k 0 k k k k k k () k k () 0 k k 0 0 ( ) 0 k 0 k () () kk k k 0 0 () () () () k k 0 k k k k k k () k k () () k k kk 0 0 () 0 k 0 k (5) k5 k5 (5) k5 k 5 () () kk k k 0 0 () () () k k 0 k k k k () () (5) k k k k k5 k5 0 () (5) 0 k k5 k k
6 (6) k6 k6 (6) k6 k SYSTEM OF SPRINGS cont. () () kk k k 0 0 () () () k k k k k k 0 () () (5) k k k k5 k k5 0 () (5) (6) 0 k k5 k k5 k6 k 6 (6) k6 k6 5 5 F 5 F 6 5 F Relaton b/w element orces and external orce Force eqlbrm F e e e e e e SYSTEM OF SPRINGS cont. 0 F,,... ND At each node, the smmaton o element orces s eqal to the appled, external orce F () () F 5 F 6 5 (5) F R () () () () () F () () (5) F () (5) (6) F (6) R 5 5 6
7 SYSTEM OF SPRINGS cont. Assembled System o Matrx Eqaton: k k k k 0 0 R k k k k k k 0 F k k k k5 k k5 0 F 0 k k5 k k5 k6 k 6 F k6 k6 5 R5 [ ]{ Q } { F } s s s [ s ] s sqare, symmetrc, snglar and postve sem-dente. When dsplacement s known, orce s nknown 5 0 R and R 5 are nknown reacton orces SYSTEM OF SPRINGS cont. Imposng Bondary Condtons Ignore the eqatons or whch the RHS orces are nknown and strke ot the correspondng rows n [ s ]. Elmnate the colmns n [ s ] that mltply nto zero vales o dsplacements o the bondary nodes. k k k k 0 0 R k k k k k k 0 F k k k k5 k k5 0 F 0 k k k k k k F k6 k6 5 R5 7
8 Global Matrx Eqaton SYSTEM OF SPRINGS cont. k k k k k F k k k k5 k 5 F k k5 k k5 k 6 F [ ]{ Q} { F} Global Stness Matrx [] sqare, symmetrc and postve dente and hence non-snglar Solton { Q} [ ] { F} Once nodal dsplacements are obtaned, sprng orces can be calclated rom e e e e P k k 5 UNIAXIAL BAR For general naxal bar, we need to dvde the bar nto a set o elements and nodes Elements are connected by sharng a node Forces are appled at the nodes (dstrbted load mst be converted to the eqvalent nodal orces) Assemble all elements n the same way wth the system o sprngs Solve the matrx eqaton or nodal dsplacements p(x) F Calclate stress and stran sng nodal dsplacements Statcally ndetermnate p(x) Statcally determnate F x 6 8
9 Two-orce member Only constant cross-secton Element orce s proportonal to relatve dspl Frst node: second code: Force-dsplacement relaton D BAR ELEMENT () e Node x L =EA/L Node A () e AE ( ) L AE ( ) L Smlar to the sprng element 7 Matrx notaton D BAR ELEMENT cont. AE L Ether orce or dsplacement (not both) mst be gven at each node. Example: = 0 and = 00 N. What happens when and are gven? Nodal eqlbrm Eqlbrm o orces actng on Node I () e Node =EA/L ( ) ( ) ( ) { e } [ k e ]{ q e } Node () e F 0 F In general F e F e Element e (e) Node (e+) Element e+ 8 9
10 Assembly D BAR ELEMENT cont. Smlar process as sprng elements Replace all nternal nodal orces wth External Appled Nodal Force Obtan system o eqatons [ ]{ Q } { F } s s s Property o [ s ] [ s ]: Strctral stness matrx {Q s } Vector o nodal DOFs {F s }: Vector o appled orces Sqare, symmetrc, postve sem-dente, snglar, non-negatve dagonal terms Applyng bondary condtons Remove rgd-body moton be xng DOFs Strkng-the-nodes and strkng-the-colmns (Reer to sprnt elements) [ ]{ Q} { F} []: Global stness matrx {Q} Vector o nknown nodal DOFs {F}: Vector o known appled orces 9 D BAR ELEMENT cont. Applyng bondary condtons cont. [] s sqare, symmetrc, postve dente, non-snglar, nvertble, and postve dagonal terms Can obtan nqe {Q} Element orces Ater solvng nodal dsplacements, the element orce can be calclated AE P L P A Element stress Reacton Forces P AE P L Note P = P Use [ s ]{Q s } = {F s }: the rows that have been deleted (strke-the-rows) Or, se e F e 0 0
11 EXAMPLE elements and nodes At node : F () () () Eqaton or each element: F Element Element x Element () () () () () () N F Element Element N N Element N F F EXAMPLE cont. How can we combne derent element eqatons? (Assembly) Frst, prepare global matrx eqaton: Dsplacement vector Stness matrx Appled orce vector Wrte the eqaton o element n the correspondng locaton 0 0 () ()
12 Wrte the eqaton o element : EXAMPLE cont. () 0 0 () Combne two eqatons o elements and 0 0 () () () 0 () Wrte the eqaton o element EXAMPLE cont () () 0 0 Combne wth other two elements () F 0 0 () () () F ( ) () F 0 0 () F 0 0 Strctral Stness Matrx
13 EXAMPLE cont. Sbsttte bondary condtons and solve or the nknown dsplacements. Let = 50 N/cm, = 0 N/cm, = 70 N/cm and = 0 N. F F 50 ( ) 0 70 F F nowns: F, F,, and Unknowns: F, F,, and ( ) 0 70 F F EXAMPLE cont. Remove zero-dsplacement colmns: and F 0 0 F 0 70 Remove nknown orce rows: F and F Now, the matrx shold not be snglar. Solve or and. Usng and, Solve or F and F.. cm 0. cm F 0 0 N F N 6
14 EXAMPLE cont. Recover element data () () Element orce () () () () N F = 0 N. cm 0. cm -8 N 7 EXAMPLE Statcally ndetermnate bars E = 00 GPa R L F = 0,000 N A = 0 m, A = 0 m Element stness matrces: () [ k ] A B C 0.5 m () [ k ] Assembly 0 F , F F 0. m R R 8
15 Applyng BC EXAMPLE cont ,000.0 m 7 Element orces or Element stresses AE P L Reacton orces () 7 P 0, N () 7 P 50 5,556N () RL P N () RR P 5,556N,, 9 PLANE TRUSS ELEMENT What s the derence between D and D nte elements? D element can move x- and y-drecton ( DOFs per node). However, the stness can be appled only axal drecton. Local Coordnate System D FE ormlaton can be sed a body-xed local coordnate system s constrcted along the length o the element Y The global coordnate system (X and Y axes) s chosen to represent the entre strctre X The local coordnate system (x and y axes) s selected to algn the x-axs along the length o the element x EA x L cm 8 cm 50 N 0 5
16 PLANE TRUSS ELEMENT cont. Element Eqaton (Local Coordnate System) Axal drecton s the local x-axs. D element eqaton x 0 0 y v EA x L 0 0 y v { } [ k]{ q} [ k] s sqare, symmetrc, postve sem-dente, and non-negatve dagonal components. How to connect to the neghborng elements? Cannot connect to other elements becase LCS s derent Use coordnate transormaton y x Local coordnates v y x v x Global coordnates x COORDINATE TRANSFORMATION Transorm to the global coord. and assemble cos sn v sn cos v cos sn v sn cos v v v Transormaton matrx cos sn 0 0 v sn cos 0 0 v 0 0 cos sn v 0 0 sn cos v local global { q} [ T]{ q} Transormaton matrx 6
17 COORDINATE TRANSFORMATION cont. The same transormaton or orce vector x cos sn 0 0 x y sn cos 0 0 y 0 0 cos sn x x y 0 0 sn cos y local global Property o transormaton matrx T T { } [ T]{ } { } [ T] T { } [ ] [ ] T { } [ T]{ } ELEMENT STIFFNESS IN GLOBAL COORD. Element x 0 0 y EA v L 0 0 x y v element stness matrx { } [ k]{ q} Transorm to the global coordnates [ T]{ } [ k][ T]{ q} y x v N { } [ T] [ k][ T] { q} global global v N x x [ k] [ T] [ k][ T] { } [ k]{ q} 7
18 ELEMENT STIFFNESS IN GLOBAL COORD. cont. Element stness matrx n global coordnates [ k] T [ k] T [ k] EA L T cos cos sn cos cos sn cos sn sn cos sn sn cos cos sn cos cos sn cos sn sn cos sn sn Depends on Yong s modls (E), cross-sectonal area (A), length (L), and angle o rotaton () Axal rgdty = EA Sqare, symmetrc, postve sem-dente, snglar, and non-negatve dagonal terms 5 EXAMPLE N 50 N Two-bar trss Dameter = 0.5 cm E = 00 6 N/cm Element 8 cm Element In local coordnate () () () { } [ k ]{ q } N cm Element N x 0 0 y EA v L 0 0 x y v y x v N v N x =.7 o E = 0 x 0 6 N/cm A = pr = 0.09 cm L =. cm x 6 8
19 Element cont. EXAMPLE cont. Element eqaton n the global coordnates () x () y v 050 () x () y v y () () () { } [ k ]{ q } x Element () x () 0 0 y v 85 () x () y 0 0 v = 90 o E = 0 x 0 6 N/cm A = pr = 0.09 cm L = 8 cm N v N v x x 7 EXAMPLE cont. Assembly Ater transormng to the global coordnates Element F x F y v F x F y v F x F y v Bondary Condtons Nodes and are xed. Node has known appled orces: F x = 50 N, F y = 0 N Element 8 9
20 EXAMPLE cont. Bondary condtons (strkng-the-colmns) F x F y v F x F y Strkng-the-rows v Solve the global matrx eqaton v 8.80 cm.80 cm 9 Spport reactons EXAMPLE cont. F x F y 8.80 N F x F y The reacton orce s parallel to the element length (two-orce member) Element orce and stress (Element ) Need to transorm to the element local coordnates v v v
21 EXAMPLE cont. Element orce and stress (Element ) cont. Element orce can only be calclated sng local element eqaton x EA N y x L y There s no orce components n the local y-drecton In x-drecton, two orces are eqal and opposte The orce n the second node s eqal to the element orce Normal stress = 60. / 0.09 = 8 N/cm. 60. N 60. N OTHER WAY OF ELEMENT FORCE CALCULATION Element orce or plane trss AE AE P L L Wrte n terms o global dsplacements AE P l mv l mv L l mv v AE L l cos m sn
22 EXAMPLE Drectly assemblng global matrx eqaton (applyng BC n the element level) Element property & drecton cosne table Elem AE/L -> l = cos m = sn > > > Snce and v wll be deleted ater assembly, t s not necessary to keep them v v l lm l lm () () EA lm m lm m v L l lm l lm lm m lm m v k F 5 () () EA l lm L lm m v k v SPACE TRUSS ELEMENT A smlar extenson sng coordnate transormaton DOF per node, v, and w x, y, and z Y v Y X w v N x Element stness matrx s 6x6 x w N Z X Z FE eqaton n the local coord. x AE x L { } [ k]{ q}
23 SPACE TRUSS ELEMENT cont. Relaton between local and global dsplacements Each node has DOFs (, v, w ) v l m n 0 0 0w { q} [ T] { q} l m n () ( 6) (6) v w Drecton cosnes x x y y z z l cos x, m cos y, n cosz L L L L x x y y z z 5 SPACE TRUSS ELEMENT cont. Relaton between local and global orce vectors Stness matrx x l 0 y m 0 z n 0 x x 0 l x y 0 m z 0 n { } [ T] T { } T T { } [ k]{ q } [ T] { } [ T] [ k][ T]{ q } { } [ k]{ q} l lm ln l lm ln m mn lm m mn v EA n ln mn n w [ k] L l lm ln sym m mn v n w [ k] [ T] T [ k][ T] 6
24 THERMAL STRESSES Temperatre change cases thermal stran L L L No stress, no stran No stress, thermal stran (a) at T = T re Thermal stress, no stran (b) at T = T re + T Constrants case thermal stresses Thermo-elastc stress-stran relatonshp = E T = + T E Thermal expanson coecent 7 THERMAL STRESSES cont. Force-dsplacement relaton L L P = AE T AE AET L L Fnte element eqaton Thermal orce vector ( ) ( ) ( ) ( ) { e } [ k e ]{ q e } { e T } ( ) 0 e v { T } AET For plane trss, transorm to the global coord. 0 v { } [ k]{ q} { } [ k]{ q} { } { } T T l m v { T} AET l mv [ ]{ Q } { F } { F } s s s Ts 8
Chapter Eight. Review and Summary. Two methods in solid mechanics ---- vectorial methods and energy methods or variational methods
Chapter Eght Energy Method 8. Introducton 8. Stran energy expressons 8.3 Prncpal of statonary potental energy; several degrees of freedom ------ Castglano s frst theorem ---- Examples 8.4 Prncpal of statonary
More informationAE/ME 339. K. M. Isaac. 8/31/2004 topic4: Implicit method, Stability, ADI method. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.
AE/ME 339 Comptatonal Fld Dynamcs (CFD) Comptatonal Fld Dynamcs (AE/ME 339) Implct form of dfference eqaton In the prevos explct method, the solton at tme level n,,n, depended only on the known vales of,
More informationOne Dimensional Axial Deformations
One Dmensonal al Deformatons In ths secton, a specfc smple geometr s consdered, that of a long and thn straght component loaded n such a wa that t deforms n the aal drecton onl. The -as s taken as the
More information1. Solve Problem 1.3-3(c) 2. Solve Problem 2.2-2(b)
. Sole Problem.-(c). Sole Problem.-(b). A two dimensional trss shown in the figre is made of alminm with Yong s modls E = 8 GPa and failre stress Y = 5 MPa. Determine the minimm cross-sectional area of
More informationFinite Element Modelling of truss/cable structures
Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures
More informationMEMBRANE ELEMENT WITH NORMAL ROTATIONS
9. MEMBRANE ELEMENT WITH NORMAL ROTATIONS Rotatons Mst Be Compatble Between Beam, Membrane and Shell Elements 9. INTRODUCTION { XE "Membrane Element" }The comple natre of most bldngs and other cvl engneerng
More informationIndeterminate pin-jointed frames (trusses)
Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all
More informationINDETERMINATE STRUCTURES METHOD OF CONSISTENT DEFORMATIONS (FORCE METHOD)
INETNTE STUTUES ETHO OF ONSISTENT EFOTIONS (FOE ETHO) If all the support reactons and nternal forces (, Q, and N) can not be determned by usng equlbrum equatons only, the structure wll be referred as STTIY
More informationC PLANE ELASTICITY PROBLEM FORMULATIONS
C M.. Tamn, CSMLab, UTM Corse Content: A ITRODUCTIO AD OVERVIEW mercal method and Compter-Aded Engneerng; Phscal problems; Mathematcal models; Fnte element method. B REVIEW OF -D FORMULATIOS Elements and
More informationSE Story Shear Frame. Final Project. 2 Story Bending Beam. m 2. u 2. m 1. u 1. m 3. u 3 L 3. Given: L 1 L 2. EI ω 1 ω 2 Solve for m 2.
SE 8 Fnal Project Story Sear Frame Gven: EI ω ω Solve for Story Bendng Beam Gven: EI ω ω 3 Story Sear Frame Gven: L 3 EI ω ω ω 3 3 m 3 L 3 Solve for Solve for m 3 3 4 3 Story Bendng Beam Part : Determnng
More informationApplication to Plane (rigid) frame structure
Advanced Computatonal echancs 18 Chapter 4 Applcaton to Plane rgd frame structure 1. Dscusson on degrees of freedom In case of truss structures, t was enough that the element force equaton provdes onl
More informationModule 11 Design of Joints for Special Loading. Version 2 ME, IIT Kharagpur
Module 11 Desgn o Jonts or Specal Loadng Verson ME, IIT Kharagpur Lesson 1 Desgn o Eccentrcally Loaded Bolted/Rveted Jonts Verson ME, IIT Kharagpur Instructonal Objectves: At the end o ths lesson, the
More informationCOMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD
COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,
More informationQuantum Mechanics I - Session 4
Quantum Mechancs I - Sesson 4 Aprl 3, 05 Contents Operators Change of Bass 4 3 Egenvectors and Egenvalues 5 3. Denton....................................... 5 3. Rotaton n D....................................
More informationMECHANICS OF MATERIALS
Fourth Edton CHTER MECHNICS OF MTERIS Ferdnand. Beer E. Russell Johnston, Jr. John T. DeWolf ecture Notes: J. Walt Oler Texas Tech Unversty Stress and Stran xal oadng Contents Stress & Stran: xal oadng
More informationPhysics 2A Chapters 6 - Work & Energy Fall 2017
Physcs A Chapters 6 - Work & Energy Fall 017 These notes are eght pages. A quck summary: The work-energy theorem s a combnaton o Chap and Chap 4 equatons. Work s dened as the product o the orce actng on
More informationDescription of the Force Method Procedure. Indeterminate Analysis Force Method 1. Force Method con t. Force Method con t
Indeternate Analyss Force Method The force (flexblty) ethod expresses the relatonshps between dsplaceents and forces that exst n a structure. Prary objectve of the force ethod s to deterne the chosen set
More informationResearch Note NONLINEAR ANALYSIS OF SEMI-RIGID FRAMES WITH RIGID END SECTIONS * S. SEREN AKAVCI **
Iranan Jornal of Scence & Technology, Transacton, Engneerng, Vol., No., pp 7-7 rnted n The Islamc Repblc of Iran, 7 Shraz Unversty Research Note NONLINER NLYSIS OF SEMI-RIGID FRMES WIT RIGID END SECTIONS
More informationPhysics 53. Rotational Motion 3. Sir, I have found you an argument, but I am not obliged to find you an understanding.
Physcs 53 Rotatonal Moton 3 Sr, I have found you an argument, but I am not oblged to fnd you an understandng. Samuel Johnson Angular momentum Wth respect to rotatonal moton of a body, moment of nerta plays
More informationAPPENDIX F A DISPLACEMENT-BASED BEAM ELEMENT WITH SHEAR DEFORMATIONS. Never use a Cubic Function Approximation for a Non-Prismatic Beam
APPENDIX F A DISPACEMENT-BASED BEAM EEMENT WITH SHEAR DEFORMATIONS Never use a Cubc Functon Approxmaton for a Non-Prsmatc Beam F. INTRODUCTION { XE "Shearng Deformatons" }In ths appendx a unque development
More informationSecond Order Analysis
Second Order Analyss In the prevous classes we looked at a method that determnes the load correspondng to a state of bfurcaton equlbrum of a perfect frame by egenvalye analyss The system was assumed to
More informationUNIVERSITY OF BOLTON RAK ACADEMIC CENTRE BENG(HONS) MECHANICAL ENGINEERING SEMESTER TWO EXAMINATION 2017/2018 FINITE ELEMENT AND DIFFERENCE SOLUTIONS
OCD0 UNIVERSITY OF BOLTON RAK ACADEMIC CENTRE BENG(HONS) MECHANICAL ENGINEERING SEMESTER TWO EXAMINATION 07/08 FINITE ELEMENT AND DIFFERENCE SOLUTIONS MODULE NO. AME6006 Date: Wednesda 0 Ma 08 Tme: 0:00
More informationNMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582
NMT EE 589 & UNM ME 48/58 ROBOT ENGINEERING Dr. Stephen Bruder NMT EE 589 & UNM ME 48/58 7. Robot Dynamcs 7.5 The Equatons of Moton Gven that we wsh to fnd the path q(t (n jont space) whch mnmzes the energy
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationPHYS 1443 Section 004 Lecture #12 Thursday, Oct. 2, 2014
PHYS 1443 Secton 004 Lecture #1 Thursday, Oct., 014 Work-Knetc Energy Theorem Work under rcton Potental Energy and the Conservatve Force Gravtatonal Potental Energy Elastc Potental Energy Conservaton o
More informationC PLANE ELASTICITY PROBLEM FORMULATIONS
C M.. Tamn, CSMLab, UTM Corse Content: A ITRODUCTIO AD OVERVIEW mercal method and Compter-Aded Engneerng; Phscal problems; Mathematcal models; Fnte element method. B REVIEW OF -D FORMULATIOS Elements and
More informationIn this section is given an overview of the common elasticity models.
Secton 4.1 4.1 Elastc Solds In ths secton s gven an overvew of the common elastcty models. 4.1.1 The Lnear Elastc Sold The classcal Lnear Elastc model, or Hooean model, has the followng lnear relatonshp
More informationModule 3: Element Properties Lecture 1: Natural Coordinates
Module 3: Element Propertes Lecture : Natural Coordnates Natural coordnate system s bascally a local coordnate system whch allows the specfcaton of a pont wthn the element by a set of dmensonless numbers
More informationCIVL 7/8111 Chapter 1 - Introduction to FEM 1/21
CIV 7/8 Chapter - to FEM / The nte element method has become a powerl tool or the nmercal solton o a wde range o engneerng problems. Applcatons range rom deormaton and stress analyss o atomotve, arcrat,
More informationMEEM 3700 Mechanical Vibrations
MEEM 700 Mechancal Vbratons Mohan D. Rao Chuck Van Karsen Mechancal Engneerng-Engneerng Mechancs Mchgan echnologcal Unversty Copyrght 00 Lecture & MEEM 700 Multple Degree of Freedom Systems (ext: S.S.
More information= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system.
Chapter Matlab Exercses Chapter Matlab Exercses. Consder the lnear system of Example n Secton.. x x x y z y y z (a) Use the MATLAB command rref to solve the system. (b) Let A be the coeffcent matrx and
More informationLinear Momentum. Center of Mass.
Lecture 6 Chapter 9 Physcs I 03.3.04 Lnear omentum. Center of ass. Course webste: http://faculty.uml.edu/ndry_danylov/teachng/physcsi Lecture Capture: http://echo360.uml.edu/danylov03/physcssprng.html
More informationPhysics 181. Particle Systems
Physcs 181 Partcle Systems Overvew In these notes we dscuss the varables approprate to the descrpton of systems of partcles, ther defntons, ther relatons, and ther conservatons laws. We consder a system
More informationˆ (0.10 m) E ( N m /C ) 36 ˆj ( j C m)
7.. = = 3 = 4 = 5. The electrc feld s constant everywhere between the plates. Ths s ndcated by the electrc feld vectors, whch are all the same length and n the same drecton. 7.5. Model: The dstances to
More informationFEA Solution Procedure
EA Soltion Procedre (demonstrated with a -D bar element problem) EA Procedre for Static Analysis. Prepare the E model a. discretize (mesh) the strctre b. prescribe loads c. prescribe spports. Perform calclations
More informationPart C Dynamics and Statics of Rigid Body. Chapter 5 Rotation of a Rigid Body About a Fixed Axis
Part C Dynamcs and Statcs of Rgd Body Chapter 5 Rotaton of a Rgd Body About a Fxed Axs 5.. Rotatonal Varables 5.. Rotaton wth Constant Angular Acceleraton 5.3. Knetc Energy of Rotaton, Rotatonal Inerta
More informationEN40: Dynamics and Vibrations. Homework 7: Rigid Body Kinematics
N40: ynamcs and Vbratons Homewor 7: Rgd Body Knematcs School of ngneerng Brown Unversty 1. In the fgure below, bar AB rotates counterclocwse at 4 rad/s. What are the angular veloctes of bars BC and C?.
More informationTranslational Equations of Motion for A Body Translational equations of motion (centroidal) for a body are m r = f.
Lesson 12: Equatons o Moton Newton s Laws Frst Law: A artcle remans at rest or contnues to move n a straght lne wth constant seed there s no orce actng on t Second Law: The acceleraton o a artcle s roortonal
More informationSCALARS AND VECTORS All physical quantities in engineering mechanics are measured using either scalars or vectors.
SCALARS AND ECTORS All phscal uanttes n engneerng mechancs are measured usng ether scalars or vectors. Scalar. A scalar s an postve or negatve phscal uantt that can be completel specfed b ts magntude.
More informationTensor Analysis. For orthogonal curvilinear coordinates, ˆ ˆ (98) Expanding the derivative, we have, ˆ. h q. . h q h q
For orthogonal curvlnear coordnates, eˆ grad a a= ( aˆ ˆ e). h q (98) Expandng the dervatve, we have, eˆ aˆ ˆ e a= ˆ ˆ a h e + q q 1 aˆ ˆ ˆ a e = ee ˆˆ ˆ + e. h q h q Now expandng eˆ / q (some of the detals
More informationECEN 5005 Crystals, Nanocrystals and Device Applications Class 19 Group Theory For Crystals
ECEN 5005 Crystals, Nanocrystals and Devce Applcatons Class 9 Group Theory For Crystals Dee Dagram Radatve Transton Probablty Wgner-Ecart Theorem Selecton Rule Dee Dagram Expermentally determned energy
More information2.3 Nilpotent endomorphisms
s a block dagonal matrx, wth A Mat dm U (C) In fact, we can assume that B = B 1 B k, wth B an ordered bass of U, and that A = [f U ] B, where f U : U U s the restrcton of f to U 40 23 Nlpotent endomorphsms
More informationDUE: WEDS FEB 21ST 2018
HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant
More informationENGN 40 Dynamics and Vibrations Homework # 7 Due: Friday, April 15
NGN 40 ynamcs and Vbratons Homework # 7 ue: Frday, Aprl 15 1. Consder a concal hostng drum used n the mnng ndustry to host a mass up/down. A cable of dameter d has the mass connected at one end and s wound/unwound
More information1 Matrix representations of canonical matrices
1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:
More informationχ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body
Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown
More informationHere is the rationale: If X and y have a strong positive relationship to one another, then ( x x) will tend to be positive when ( y y)
Secton 1.5 Correlaton In the prevous sectons, we looked at regresson and the value r was a measurement of how much of the varaton n y can be attrbuted to the lnear relatonshp between y and x. In ths secton,
More informationMomentum. Momentum. Impulse. Momentum and Collisions
Momentum Momentum and Collsons From Newton s laws: orce must be present to change an object s elocty (speed and/or drecton) Wsh to consder eects o collsons and correspondng change n elocty Gol ball ntally
More informationFrame element resists external loads or disturbances by developing internal axial forces, shear forces, and bending moments.
CE7 Structural Analyss II PAAR FRAE EEET y 5 x E, A, I, Each node can translate and rotate n plane. The fnal dsplaced shape has ndependent generalzed dsplacements (.e. translatons and rotatons) noled.
More informationAPPENDIX A Some Linear Algebra
APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,
More informationA Tale of Friction Basic Rollercoaster Physics. Fahrenheit Rollercoaster, Hershey, PA max height = 121 ft max speed = 58 mph
A Tale o Frcton Basc Rollercoaster Physcs Fahrenhet Rollercoaster, Hershey, PA max heght = 11 t max speed = 58 mph PLAY PLAY PLAY PLAY Rotatonal Movement Knematcs Smlar to how lnear velocty s dened, angular
More informationWeek 9 Chapter 10 Section 1-5
Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,
More informationSolutions to Problem Set 6
Solutons to Problem Set 6 Problem 6. (Resdue theory) a) Problem 4.7.7 Boas. n ths problem we wll solve ths ntegral: x sn x x + 4x + 5 dx: To solve ths usng the resdue theorem, we study ths complex ntegral:
More informationChat eld, C. and A.J.Collins, Introduction to multivariate analysis. Chapman & Hall, 1980
MT07: Multvarate Statstcal Methods Mke Tso: emal mke.tso@manchester.ac.uk Webpage for notes: http://www.maths.manchester.ac.uk/~mkt/new_teachng.htm. Introducton to multvarate data. Books Chat eld, C. and
More informationCME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 13
CME 30: NUMERICAL LINEAR ALGEBRA FALL 005/06 LECTURE 13 GENE H GOLUB 1 Iteratve Methods Very large problems (naturally sparse, from applcatons): teratve methods Structured matrces (even sometmes dense,
More informationArmy Ants Tunneling for Classical Simulations
Electronc Supplementary Materal (ESI) for Chemcal Scence. Ths journal s The Royal Socety of Chemstry 2014 electronc supplementary nformaton (ESI) for Chemcal Scence Army Ants Tunnelng for Classcal Smulatons
More information3/31/ = 0. φi φi. Use 10 Linear elements to solve the equation. dx x dx THE PETROV-GALERKIN METHOD
THE PETROV-GAERKIN METHO Consder the Galern solton sng near elements of the modfed convecton-dffson eqaton α h d φ d φ + + = α s a parameter between and. If α =, we wll have the dscrete Galern form of
More informationNEWTON S LAWS. These laws only apply when viewed from an inertial coordinate system (unaccelerated system).
EWTO S LAWS Consder two partcles. 1 1. If 1 0 then 0 wth p 1 m1v. 1 1 2. 1.. 3. 11 These laws only apply when vewed from an nertal coordnate system (unaccelerated system). consder a collecton of partcles
More informationPeriod & Frequency. Work and Energy. Methods of Energy Transfer: Energy. Work-KE Theorem 3/4/16. Ranking: Which has the greatest kinetic energy?
Perod & Frequency Perod (T): Tme to complete one ull rotaton Frequency (): Number o rotatons completed per second. = 1/T, T = 1/ v = πr/t Work and Energy Work: W = F!d (pcks out parallel components) F
More informationModeling of Dynamic Systems
Modelng of Dynamc Systems Ref: Control System Engneerng Norman Nse : Chapters & 3 Chapter objectves : Revew the Laplace transform Learn how to fnd a mathematcal model, called a transfer functon Learn how
More informationPHYS 705: Classical Mechanics. Newtonian Mechanics
1 PHYS 705: Classcal Mechancs Newtonan Mechancs Quck Revew of Newtonan Mechancs Basc Descrpton: -An dealzed pont partcle or a system of pont partcles n an nertal reference frame [Rgd bodes (ch. 5 later)]
More informationTEST 5 (phy 240) 2. Show that the volume coefficient of thermal expansion for an ideal gas at constant pressure is temperature dependent and given by
ES 5 (phy 40). a) Wrte the zeroth law o thermodynamcs. b) What s thermal conductvty? c) Identyng all es, draw schematcally a P dagram o the arnot cycle. d) What s the ecency o an engne and what s the coecent
More informationStructure and Drive Paul A. Jensen Copyright July 20, 2003
Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.
More informationChapter 3 and Chapter 4
Chapter 3 and Chapter 4 Chapter 3 Energy 3. Introducton:Work Work W s energy transerred to or rom an object by means o a orce actng on the object. Energy transerred to the object s postve work, and energy
More informationChapter 07: Kinetic Energy and Work
Chapter 07: Knetc Energy and Work Conservaton o Energy s one o Nature s undamental laws that s not volated. Energy can take on derent orms n a gven system. Ths chapter we wll dscuss work and knetc energy.
More informationInner Product. Euclidean Space. Orthonormal Basis. Orthogonal
Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationComplex Variables. Chapter 18 Integration in the Complex Plane. March 12, 2013 Lecturer: Shih-Yuan Chen
omplex Varables hapter 8 Integraton n the omplex Plane March, Lecturer: Shh-Yuan hen Except where otherwse noted, content s lcensed under a BY-N-SA. TW Lcense. ontents ontour ntegrals auchy-goursat theorem
More informationWeek 11: Chapter 11. The Vector Product. The Vector Product Defined. The Vector Product and Torque. More About the Vector Product
The Vector Product Week 11: Chapter 11 Angular Momentum There are nstances where the product of two vectors s another vector Earler we saw where the product of two vectors was a scalar Ths was called the
More informationLecture Note 3. Eshelby s Inclusion II
ME340B Elastcty of Mcroscopc Structures Stanford Unversty Wnter 004 Lecture Note 3. Eshelby s Incluson II Chrs Wenberger and We Ca c All rghts reserved January 6, 004 Contents 1 Incluson energy n an nfnte
More informationVirtual Work 3rd Year Structural Engineering
Vrtual Work 3rd Year Structural Engneerng 2/ Dr. Coln Capran Contents. Introducton... 4. General... 4.2 Background... 5 2. The Prncple of Vrtual Work... 4 2. Defnton... 4 2.2 Vrtual Dsplacements... 7 2.3
More informationProf. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model
EXACT OE-DIMESIOAL ISIG MODEL The one-dmensonal Isng model conssts of a chan of spns, each spn nteractng only wth ts two nearest neghbors. The smple Isng problem n one dmenson can be solved drectly n several
More informationLesson 5: Kinematics and Dynamics of Particles
Lesson 5: Knematcs and Dynamcs of Partcles hs set of notes descrbes the basc methodology for formulatng the knematc and knetc equatons for multbody dynamcs. In order to concentrate on the methodology and
More informationForce = F Piston area = A
CHAPTER III Ths chapter s an mportant transton between the propertes o pure substances and the most mportant chapter whch s: the rst law o thermodynamcs In ths chapter, we wll ntroduce the notons o heat,
More informationFirst day August 1, Problems and Solutions
FOURTH INTERNATIONAL COMPETITION FOR UNIVERSITY STUDENTS IN MATHEMATICS July 30 August 4, 997, Plovdv, BULGARIA Frst day August, 997 Problems and Solutons Problem. Let {ε n } n= be a sequence of postve
More informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
More informationPHYS 1441 Section 002 Lecture #16
PHYS 1441 Secton 00 Lecture #16 Monday, Mar. 4, 008 Potental Energy Conservatve and Non-conservatve Forces Conservaton o Mechancal Energy Power Today s homework s homework #8, due 9pm, Monday, Mar. 31!!
More informationPhysics 2A Chapter 3 HW Solutions
Phscs A Chapter 3 HW Solutons Chapter 3 Conceptual Queston: 4, 6, 8, Problems: 5,, 8, 7, 3, 44, 46, 69, 70, 73 Q3.4. Reason: (a) C = A+ B onl A and B are n the same drecton. Sze does not matter. (b) C
More informationSolutions to selected problems from homework 1.
Jan Hagemejer 1 Soltons to selected problems from homeork 1. Qeston 1 Let be a tlty fncton hch generates demand fncton xp, ) and ndrect tlty fncton vp, ). Let F : R R be a strctly ncreasng fncton. If the
More informationFirst Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force.
Secton 1. Dynamcs (Newton s Laws of Moton) Two approaches: 1) Gven all the forces actng on a body, predct the subsequent (changes n) moton. 2) Gven the (changes n) moton of a body, nfer what forces act
More informationSolutions to exam in SF1811 Optimization, Jan 14, 2015
Solutons to exam n SF8 Optmzaton, Jan 4, 25 3 3 O------O -4 \ / \ / The network: \/ where all lnks go from left to rght. /\ / \ / \ 6 O------O -5 2 4.(a) Let x = ( x 3, x 4, x 23, x 24 ) T, where the varable
More informationI have not received unauthorized aid in the completion of this exam.
ME 270 Sprng 2013 Fnal Examnaton Please read and respond to the followng statement, I have not receved unauthorzed ad n the completon of ths exam. Agree Dsagree Sgnature INSTRUCTIONS Begn each problem
More informationPrinciple of virtual work
Ths prncple s the most general prncple n mechancs 2.9.217 Prncple of vrtual work There s Equvalence between the Prncple of Vrtual Work and the Equlbrum Equaton You must know ths from statc course and dynamcs
More informationPhysics 5153 Classical Mechanics. Principle of Virtual Work-1
P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal
More informationFormulas for the Determinant
page 224 224 CHAPTER 3 Determnants e t te t e 2t 38 A = e t 2te t e 2t e t te t 2e 2t 39 If 123 A = 345, 456 compute the matrx product A adj(a) What can you conclude about det(a)? For Problems 40 43, use
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationSome Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)
Some Comments on Acceleratng Convergence of Iteratve Sequences Usng Drect Inverson of the Iteratve Subspace (DIIS) C. Davd Sherrll School of Chemstry and Bochemstry Georga Insttute of Technology May 1998
More informationC/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1
C/CS/Phy9 Problem Set 3 Solutons Out: Oct, 8 Suppose you have two qubts n some arbtrary entangled state ψ You apply the teleportaton protocol to each of the qubts separately What s the resultng state obtaned
More informationPHYS 1441 Section 002 Lecture #15
PHYS 1441 Secton 00 Lecture #15 Monday, March 18, 013 Work wth rcton Potental Energy Gravtatonal Potental Energy Elastc Potental Energy Mechancal Energy Conservaton Announcements Mdterm comprehensve exam
More informatione - c o m p a n i o n
OPERATIONS RESEARCH http://dxdoorg/0287/opre007ec e - c o m p a n o n ONLY AVAILABLE IN ELECTRONIC FORM 202 INFORMS Electronc Companon Generalzed Quantty Competton for Multple Products and Loss of Effcency
More informationNew Method for Solving Poisson Equation. on Irregular Domains
Appled Mathematcal Scences Vol. 6 01 no. 8 369 380 New Method for Solvng Posson Equaton on Irregular Domans J. Izadan and N. Karamooz Department of Mathematcs Facult of Scences Mashhad BranchIslamc Azad
More information11. Dynamics in Rotating Frames of Reference
Unversty of Rhode Island DgtalCommons@URI Classcal Dynamcs Physcs Course Materals 2015 11. Dynamcs n Rotatng Frames of Reference Gerhard Müller Unversty of Rhode Island, gmuller@ur.edu Creatve Commons
More informationIntroduction to elastic wave equation. Salam Alnabulsi University of Calgary Department of Mathematics and Statistics October 15,2012
Introdcton to elastc wave eqaton Salam Alnabls Unversty of Calgary Department of Mathematcs and Statstcs October 15,01 Otlne Motvaton Elastc wave eqaton Eqaton of moton, Defntons and The lnear Stress-
More informationPHYS 1443 Section 002 Lecture #20
PHYS 1443 Secton 002 Lecture #20 Dr. Jae Condtons for Equlbru & Mechancal Equlbru How to Solve Equlbru Probles? A ew Exaples of Mechancal Equlbru Elastc Propertes of Solds Densty and Specfc Gravty lud
More informationFormal solvers of the RT equation
Formal solvers of the RT equaton Formal RT solvers Runge- Kutta (reference solver) Pskunov N.: 979, Master Thess Long characterstcs (Feautrer scheme) Cannon C.J.: 970, ApJ 6, 55 Short characterstcs (Hermtan
More informationLinear Algebra and its Applications
Lnear Algebra and ts Applcatons 4 (00) 5 56 Contents lsts avalable at ScenceDrect Lnear Algebra and ts Applcatons journal homepage: wwwelsevercom/locate/laa Notes on Hlbert and Cauchy matrces Mroslav Fedler
More informationAdvanced topics in Finite Element Method 3D truss structures. Jerzy Podgórski
Advanced topics in Finite Element Method 3D trss strctres Jerzy Podgórski Introdction Althogh 3D trss strctres have been arond for a long time, they have been sed very rarely ntil now. They are difficlt
More informationGeneral displacement arch-cantilever element method for stress analysis of arch dam
Water Scence and Engneerng, 009, (): 3-4 do:0.388/j.ssn.674-370.009.0.004 http://kkb.hhu.edu.cn e-mal: wse@hhu.edu.cn General dsplacement arch-cantlever element method for stress analyss of arch dam Hao
More informationLectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix
Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could
More informationU.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017
U.C. Berkeley CS94: Beyond Worst-Case Analyss Handout 4s Luca Trevsan September 5, 07 Summary of Lecture 4 In whch we ntroduce semdefnte programmng and apply t to Max Cut. Semdefnte Programmng Recall that
More information