Technical Report TR05
|
|
- Tamsyn Hall
- 5 years ago
- Views:
Transcription
1 Techncal Report TR05 An Introducton to the Floatng Frame of Reference Formulaton for Small Deformaton n Flexble Multbody Dynamcs Antono Recuero and Dan Negrut May 11, 2016
2 Abstract Ths techncal report descrbes the bass of the nonlnear fnte element floatng frame of reference formulaton. Keywords: Floatng Frame of Reference Formulaton, Nonlnear Fnte Elements, Flexble Multbody Systems 1
3 Contents 1 Introducton 3 2 Floatng Frame of Reference Knematcs of the FFR Equatons of Moton Fnte Element FFR Formulaton Local Deformaton and Reference Frames Descrpton of Local Deformaton Reference Condtons Knematcs of an FFR Fnte Element Model Order Reducton
4 1 Introducton The dea of a floatng frame of reference (FFR) may have ts orgn n the analyss of astrodynamc systems. Such systems experence strong couplng of nerta forces, due to hgh acceleratons n ts reference (rgd body) moton, and flexblty due to the lght nature of the structures nvolved (see Ref. [2]). The dea of a body-wse co-rotated frame or floatng frame can be summarzed as the splttng of the overall moton of bodes that experence small deformaton nto a frame that captures the dynamcs of mean rgd body moton and a superposed flexble moton that defnes the deformaton state. Ths method requres one flexble body to have one sngle FFR to descrbe ts dynamcs, n contrast wth the co-rotatonal formulaton. In addton, the fact that the body s (structure s) elastcty can be defned va a constant stffness matrx opens up the possblty of usng general model order reducton technques to reduce the order of a lnear fnte element model. In what follows, the bass of the FFR formulaton s presented, whereas many dervatons are left to the reader as an exercse or gven as a reference. 2 Floatng Frame of Reference 2.1 Knematcs of the FFR From now on, we wll assume that Euler parameters wll be our choce of rotaton coordnates. Even though other choces can be made, t s, n any case, mportant to guarantee that sngulartes wll be avoded. The knematc state of the frame of reference of a body may be wrtten as [ qr = R T θ TT], (1) where q r s a set of Cartesan coordnates that locates the orgn of the FFR of body, O, wth respect to the global frame, θ s a set of Euler parameters that descrbe the orentaton of such a frame n the global coordnate system. A pont P n body can be descrbed, usng FFR knematcs, as r P = R + A ( θ ) ū P, (2) where A denotes the orentaton matrx of the FFR and ū P s the vector that defnes the poston of pont P wth respect to the FFR (superscrpt s dropped for smplcty). Ths vector can be splt nto ts reference poston and ts elastc dsplacement as follows ū P = ū P o + ū P f. (3) ū P o represents the undeformed poston of P, whereas ts elastc dsplacement, measured n the FFR, s denoted by the vector ū P f (see Fg. 1). The elastc dsplacement vector ūp f may be decomposed nto the product of a space-dependent matrx and a vector of tme-dependent flexble coordnates, as follows ū P f = S ( ) ū P o q f, (4) 3
5 where S (ū P o ) s a shape functon matrx referred to the reference (undeformed) confguraton and qf s the vector of flexble coordnates of body. By spellng out the knematc equatons n (3) and (4), the expresson of the deformed poston of the pont P n body may be wrtten as rp = R + A ( ) ū P o + S qf (5) Velocty expressons of a deformed pont P wthn the FFR formulaton are gven below: ū P f = S ( ū P o ) q f, (6) ṙ P = Ṙ + A ( ū P o + S q f) + A S q f. (7) Note that Eqs. (6)-(7) take nto account that ū P o and S do not depend on tme. Fgure 1: Basc knematcs of the FFR formulaton. The floatng frame axes are denoted by X Y Z, the ntal, unstressed confguraton s depcted wth sold lnes, whereas the current confguraton translated, rotated, and deformed s shown n dashed lnes. 2.2 Equatons of Moton Here, we present the dervaton equatons of moton of the system n abbrevated manner, omttng some dervatons. The nterested reader can consult the book by Shabana [8] to fnd out more on some detaled dervatons of the FFR formulaton equatons. Before movng on wth dervatons, we arrange, for convenence, the velocty vector n the followng form: Ṙ ṙp = [ I B A S ] θ = L q, (8) where matrx B s a lnear operator that acts on the rotaton coordnates, as follows A ū P = B θ, and B = A ū P Ḡ. The tlde denotes the skew symmetrc matrx 4 q f
6 operaton on a vector, and G s a lnear operator that transforms tme dervatves of rotaton parameters to the angular velocty vector of the reference moton of body. The form of G depends upon the selecton of rotatonal parameters. The equatons of moton of the FFR formulatons are to be derved usng Lagrange s equatons. For a flexble body descrbed usng a floatng frame, these equatons take the form d dt ( ) T T ( T q q ) T + C T q λ = Q e + Q a, (9) where T s the knetc energy, q and q are the generalzed coordnates and veloctes of body, respectvely, C q s the Jacoban of the constrant equatons, λ s the vector of Lagrange multplers, and Q e and Q a are the vectors of elastc and appled forces, respectvely. Note that, n general, the vectors of Lagrange multplers and appled forces may lnk the flexble body wth other (flexble) bodes and/or other nteractng systems. Mass Matrx The total knetc energy of body may be expressed as a volume ntegral: T = 1 ρ ṙ T 2 P ṙp dv (10) whch, n terms of the generalzed coordnates, may be wrtten as V T = 1 2 qt M q. (11) The mass matrx M s nonlnear n the rotaton coordnates and takes the followng form M = V ρ I B A S m RR m Rθ m Rf B T B B T A S dv = m θθ m θf, (12) sym. S T S sym. where submatrces have been named for convenence. Cross terms n the mass matrx, e.g. m Rf and m θf, ndcate that couplng between reference and flexble coordnates s captured. The Quadratc Velocty Vector The FFR formulaton ntroduces complex nerta terms whch represent Corols and centrfugal forces. Dervng the nerta terms of Lagrange equatons n Eq. (9), one obtans m ff ( ) d T T ( ) T T [ ( = M q + Ṁ q q T M q ) dt q q q } {{} Q v(quadratc velocty vector) ] T (13) 5
7 [ The fnal expresson for the quadratc velocty vector Q v = (Q v ) T R (Q v ) T θ terms of generalzed coordnates, reads ] T, (Q v ) T f n ( Q v ) ( Q v ) R = A ρ θ = ḠT ρ ( ) Q = v f ρ V V V S T [ ) ] ( ω 2ū + 2 ω S q f dv (14) [ ū T ) ] ( ω 2ū + 2 ū T ω S q f dv (15) [ ) ] ( ω 2ū + 2 ω S q f dv (16) and may be obtaned by dervng equaton (13) or va the Prncple of Vrtual Work. The Euler parameter dentty Ḡ θ = 0 s used above to smplfy the expressons. Checkng the correctness of (14) s left to the reader as an exercse. As an advanced topc, the ready may learn on the accurate calculaton of FFR velocty-dependent nerta terms n [9]. Appled Forces To nclude forces and moments n the system t s necessary to obtan ther generalzed counterparts by usng FFR knematcs. As an example, here we descrbe how to compute the generalzed force assocated wth a pont force. The vrtual work of a force descrbed n the global frame, F a, appled at P may be expressed as (17) δw a = F aδr P = Q aδq. (18) As a frst step, we wrte the varaton of the poston vector of a pont P n an FFR body as δrp = [I A ū ] δr Ḡ A S δθ. (19) Equaton (19) can be plugged nto (18) to yeld the followng generalzed force vector ( Q a ) R = F P a, ( Q a ) δq f θ = ḠT ū T A T F P a, ( Q a ) f = ST A T F P a, (20) where space-dependent vector and matrx, ū and S, respectvely, must be partcularzed at the poston of pont P. Note that the FFR formulaton does not mpose lmtatons as to the locaton or nature of the load. 6
8 Equatons of Moton In general form, the FFR equatons of moton of a body can be wrtten as follows M q + D q + K q + C T q λ = Q a + Q v, (21) whch, n turn, may be expressed n terms of submatrces as m RR m Rθ m Rf R Ṙ R m θθ m θf θ θ θ + sym. m ff q f 0 0 D ff q f 0 0 K ff q f C T R (Q a ) R (Q v ) R C θ λ = (Q a ) θ + (Q v ) θ. C T q f (Q a ) q f (Q v ) q f The new terms K ff and D ff refer to the stffness and dampng matrces of the structure. These two matrces may be obtaned by modelng the structure as a lnear system usng, for nstance, the fnte element method. 3 Fnte Element FFR Formulaton As shown n the prevous secton, the FFR formulaton leads to a local lnear problem that enables the use of classcal fnte element technology: Flexble coordnates can be consdered as fnte element degrees of freedom. Wthn the context of FFR, fnte elements contanng nfntesmal angles stll accurately descrbe rgd body moton. Ths fact, together wth the use of ntermedate element coordnate systems, allows to obtan exact modelng of rgd body nerta for flexble bodes. The ntroducton of ntermedate element frames, as descrbed by Shabana n [8], elmnates the ssue of ensurng that actual small strans are seen as small strans when descrbed n the body s floatng frame of reference. Complex geometrc shapes can be modeled by usng one sngle, co-rotated floatng frame per body due to the ntroducton of ntermedate frames. The co-rotatonal formulaton, whch uses element-wse coordnates to descrbe the state of fnte elements, also appeared due to need for modelng complex geometres [4]. The floatng frame of reference must be attach to the flexble body by prescrbng a set of condtons that wll ensure that no rgd body s degrees of freedom reman n the structure. That s, the FFR s degrees of freedom must descrbe reference moton, whereas the fnte elements coordnates must only descrbe dsplacements that lead to deformaton. Ths s accomplshed by mposng reference condtons, whch lnk the fnte element model to the FFR. The success of the FFR n practcal multbody applcatons owes much to the fact the t enables the use of lnear model order reducton technques, that s, numercal methods amng 7
9 at reducng the number of coordnates of a lnear system whle preservng statc/dynamc nformaton of the orgnal, large-scale model. Ths dea wll be elaborated on n a later subsecton. 3.1 Local Deformaton and Reference Frames We wll ntroduce a more complete notaton n ths subsecton to refer to fnte elements. The poston of a fnte element j n a body s expressed as follows: where e j 0 and e j f e j = e j 0 + e j f, (22) denote the rgd body and flexble coordnates of a fnte element j n a body. Several coordnate systems needed for a general problem are ntroduced: 1) X Y Z s a body coordnate system that represents the overall moton of the body but need not be rgdly attached to the flexble body. 2) An element-wse reference system X j Y j Z j s rgdly attached to the fnte element j. 3) Fnally, at the locaton of the body coordnate system, an ntermedate element coordnate system wth axes ntally parallel to X j Y j Z j s placed, whch s denoted by X j. Ths ntermedate frame s necessary to guarantee exact modelng of rgd body knematcs [8]. A way to nterpret the need for ths system s the followng: Snce the FFR formulaton assumes one sngle reference frame per body, there s a need to ensure that we account for the reference orentaton at the unstressed confguraton. Note that the ntermedate element coordnate system does not ntroduce addtonal coordnates n the formulaton snce t s defned by the geometry of the flexble body at the ntal confguraton. Y j Z j Fgure 2: The ntermedate frame of reference X j Y j Z j s shown n the fgure as a rotated coordnated system whch s parallel to the ntal fnte element frame, X j Y j Z j. The reference frame X j Y j Z j allows the transformaton of the fnte element coordnates n the element system, e j, nto the body reference frame, qj n. 8
10 3.2 Descrpton of Local Deformaton Connectvty condtons n the flexble body are establshed, as n any fnte element model, by dentfyng the nodes shared by elements. From now on, the connectvty problem wll be dsregarded and the dervatons wll be presented element-wse, wthout loss of generalty. To descrbe the knematcs of the deformaton of a fnte element j n a body, we need frst to defne the orentaton matrx descrbng the orentaton of the ntermedate element coordnate system wth respect to the body coordnate system or floatng frame, C j. The assumed dsplacement feld w j of an element can be expressed n the body frame as ū j = C j w j = C j S j e j, (23) where S j s the shape functon matrx of an element j n a body and e j element j s coordnates. Note that w j s the vector of s defned wth respect to the ntermedate element may be defned n the coordnate system (see Fg. 2). The vector of nodal coordnates e j body coordnate system as = C jt qn j, (24) n whch qn j s the vector of element j s nodal coordnates expressed n the body coordnate system. Makng use of Eqs. (23)-(24), one gets e j ū j = C j w j = C j S j C jt q j n. (25) Equaton (25) s key snce t expresses the poston coordnates of an arbtrary pont on the fnte element wth respect to the orgn of the body coordnate system [7]. 3.3 Reference Condtons Cartesan and rotatonal coordnates descrbe the locaton and atttude of the floatng frame of reference. Flexble coordnates, on the other hand, descrbe small deformaton of the body. However, so far, we have skpped descrbng how to elmnate the rgd body modes assocated wth the fnte element coordnates. In other words, the model descrbng the deformaton of the flexble body must be attached to the floatng frame by usng a set of reference condtons that must fully elmnate relatve rgd body moton. Reference condtons can mposed n dfferent ways; some of them are summarzed below: Free body. If the flexble body s not connected to other bodes, a sensble soluton to the problem of assgnng reference condtons to the flexble body s to keep the orgn of the coordnate system at the nstantaneous mass center and, smultaneously, ensure that ts axes are prncpal, that s, they defne the three products of nerta reman zero as the body deforms [1]. The reference condtons for mean axes may be mposed by selectng a proper set of algebrac constrants [6]. Clamped condtons. A smpler way of defnng reference condtons s by fully constranng the body coordnate system to be located and orented as one node at a fnte element. For nstance, the body frame may be fully constraned to the node lnkng a rotor blade to the rotor tself. 9
11 Smply supported. For beams, t s possble to attach the body coordnate frame to the flexble body n such a way that ts locaton s placed at one end of the beam and one of ts axes always ponts to the other end. Note that ths s the case of many robotc flexble manpulators. The choce of reference condtons must be consstent wth the problem at hand snce, n a practcal multbody system problem, the way flexble bodes are connected to other bodes may determne the way they deform. In other words, reference condtons must be selected such that they respect both the boundary condtons of the flexble body wthn the system and the separate model used to descrbe the lnear deformaton problem (e.g. the fnte element model). 3.4 Knematcs of an FFR Fnte Element For convenence, we rewrte Equaton (25) n the followng manner ū j = N j q j n, (26) where N j s a space-dependent matrx: N j = C j w j = C j S j C jt. Usng the new knematc descrpton that accounts for the use of fnte element technques, the poston vector r j of an arbtrary pont on element j of a body may be wrtten as r j = R + A ū j = R + A N j qn j = R + A N ( ) j qo + Brq j f, (27) where qo s the vector of nodal coordnates n the undeformed state,.e. they locate materal ponts wth respect to the body frame n the unformed state, qf s the vector of nodal deformatons, and Br s a matrx that may apply a set of reference condtons on the deformaton coordnates, e.g., n order to attach a node to the locaton of the reference frame. Equaton (27) may be dfferentated wth respect to tme to obtan ṙ j = Ṙ + A ( ω ū ) + A N j B r q j f. (28) Note that n Eq. (28) only matrx N j has spatal dependency. The knetc and elastc terms of the equatons of moton of the fnte element FFR formulaton can therefore be obtaned usng Eqs. (27)-(28). A set of nerta shape ntegrals nvolvng the fnte element shape functons and ther relatve orentaton wth respect to the body frame of reference may be calculated at a preprocessng stage as descrbed n Ref. [8]. Usng the knematc equatons presented n ths subsecton, Eqs. (12), (14), (20), and (21) can be rederved for the fnte element FFR formulaton. 3.5 Model Order Reducton The equatons of moton of the fnte element FFR formulaton may contan an excessve number of coordnates owng to the detaled fnte element models used to descrbe a flexble 10
12 body s dynamcs. The sngle frame of reference ntroduced to capture flexble body dynamcs allows for a separaton of rgd body and flexble coordnates. Snce flexble coordnates n the FFR formulaton descrbe lnear deformaton, any model order reducton technque may be appled. Tradtonally, methods such as Component Mode Synthess, Guyan condensaton, and purely egenmode analyss have been profusely employed to reduce the number of coordnates of the flexble body. In the past decade, there has been an nterest for the use of modern model order reducton technques n multbody system dynamcs, see [3, 5]. Proper use of model order reducton allows to analyze complex system at a hgh level of accuracy and reduced computatonal cost. References [1] OP Agrawal and AA Shabana. Applcaton of deformable-body mean axs to flexble multbody system dynamcs. Computer Methods n Appled Mechancs and Engneerng, 56(2): , [2] B Fraejs De Veubeke. The dynamcs of flexble bodes. Internatonal Journal of Engneerng Scence, 14(10): , [3] Jörg Fehr and Peter Eberhard. Smulaton process of flexble multbody systems wth non-modal model order reducton technques. Multbody System Dynamcs, 25(3): , [4] CA Felppa and B Haugen. A unfed formulaton of small-stran corotatonal fnte elements: I. theory. Computer Methods n Appled Mechancs and Engneerng, 194(21): , [5] P Koutsovasls and M Betelschmdt. Comparson of model reducton technques for large mechancal systems. Multbody System Dynamcs, 20(2): , [6] Parvz E Nkravesh and Y-shh Ln. Use of prncpal axes as the floatng reference frame for a movng deformable body. Multbody System Dynamcs, 13(2): , [7] AA Shabana. Dynamcs of large scale flexble mechancal systems. PhD thess, Ph. D. dssertaton, Unversty of Iowa, Iowa Cty, [8] Ahmed A Shabana. Dynamcs of multbody systems. Cambrdge unversty press, [9] Karm Sherf and Karn Nachbagauer. A detaled dervaton of the velocty-dependent nerta forces n the floatng frame of reference formulaton. Journal of Computatonal and Nonlnear Dynamcs, 9(4):044501,
Physics 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 informationPhysics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1
P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the
More informationIterative General Dynamic Model for Serial-Link Manipulators
EEL6667: Knematcs, Dynamcs and Control of Robot Manpulators 1. Introducton Iteratve General Dynamc Model for Seral-Lnk Manpulators In ths set of notes, we are gong to develop a method for computng a general
More informationcoordinates. Then, the position vectors are described by
Revewng, what we have dscussed so far: Generalzed coordnates Any number of varables (say, n) suffcent to specfy the confguraton of the system at each nstant to tme (need not be the mnmum number). In general,
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 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 informationSpin-rotation coupling of the angularly accelerated rigid body
Spn-rotaton couplng of the angularly accelerated rgd body Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 E-mal: louaelzen@gmal.com November 1, 2017 All Rghts Reserved. Abstract Ths paper s
More informationAn Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors
An Algorthm to Solve the Inverse Knematcs Problem of a Robotc Manpulator Based on Rotaton Vectors Mohamad Z. Al-az*, Mazn Z. Othman**, and Baker B. Al-Bahr* *AL-Nahran Unversty, Computer Eng. Dep., Baghdad,
More informationCanonical transformations
Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,
More informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
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 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 informationProfessor Terje Haukaas University of British Columbia, Vancouver The Q4 Element
Professor Terje Haukaas Unversty of Brtsh Columba, ancouver www.nrsk.ubc.ca The Q Element Ths document consders fnte elements that carry load only n ther plane. These elements are sometmes referred to
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 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 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 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 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 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 informationPHYS 705: Classical Mechanics. Calculus of Variations II
1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary
More informationNovember 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project
SE 8 Fnal Project Story Shear Frame u m Gven: u m L L m L L EI ω ω Solve for m Story Bendng Beam u u m L m L Gven: m L L EI ω ω Solve for m 3 3 Story Shear Frame u 3 m 3 Gven: L 3 m m L L L 3 EI ω ω ω
More informationMEV442 Introduction to Robotics Module 2. Dr. Santhakumar Mohan Assistant Professor Mechanical Engineering National Institute of Technology Calicut
MEV442 Introducton to Robotcs Module 2 Dr. Santhakumar Mohan Assstant Professor Mechancal Engneerng Natonal Insttute of Technology Calcut Jacobans: Veloctes and statc forces Introducton Notaton for tme-varyng
More informationCHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics)
CHAPTER 6 LAGRANGE S EQUATIONS (Analytcal Mechancs) 1 Ex. 1: Consder a partcle movng on a fxed horzontal surface. r P Let, be the poston and F be the total force on the partcle. The FBD s: -mgk F 1 x O
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 informationRigid body simulation
Rgd bod smulaton Rgd bod smulaton Once we consder an object wth spacal etent, partcle sstem smulaton s no longer suffcent Problems Problems Unconstraned sstem rotatonal moton torques and angular momentum
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 informationThe classical spin-rotation coupling
LOUAI H. ELZEIN 2018 All Rghts Reserved The classcal spn-rotaton couplng Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 louaelzen@gmal.com Abstract Ths paper s prepared to show that a rgd
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 informationLecture 8 Modal Analysis
Lecture 8 Modal Analyss 16.0 Release Introducton to ANSYS Mechancal 1 2015 ANSYS, Inc. February 27, 2015 Chapter Overvew In ths chapter free vbraton as well as pre-stressed vbraton analyses n Mechancal
More informationAppendix B. The Finite Difference Scheme
140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton
More informationClassical Mechanics Virtual Work & d Alembert s Principle
Classcal Mechancs Vrtual Work & d Alembert s Prncple Dpan Kumar Ghosh UM-DAE Centre for Excellence n Basc Scences Kalna, Mumba 400098 August 15, 2016 1 Constrants Moton of a system of partcles s often
More informationKinematics of Fluids. Lecture 16. (Refer the text book CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlines) 17/02/2017
17/0/017 Lecture 16 (Refer the text boo CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlnes) Knematcs of Fluds Last class, we started dscussng about the nematcs of fluds. Recall the Lagrangan and Euleran
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationNUMERICAL RESULTS QUALITY IN DEPENDENCE ON ABAQUS PLANE STRESS ELEMENTS TYPE IN BIG DISPLACEMENTS COMPRESSION TEST
Appled Computer Scence, vol. 13, no. 4, pp. 56 64 do: 10.23743/acs-2017-29 Submtted: 2017-10-30 Revsed: 2017-11-15 Accepted: 2017-12-06 Abaqus Fnte Elements, Plane Stress, Orthotropc Materal Bartosz KAWECKI
More informationMathematical Preparations
1 Introducton Mathematcal Preparatons The theory of relatvty was developed to explan experments whch studed the propagaton of electromagnetc radaton n movng coordnate systems. Wthn expermental error the
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 informationCHAIN DYNAMIC FORMULATIONS FOR MULTIBODY SYSTEM TRACKED VEHICLES
2012 NDIA GROUND VEHICLE SYSEMS ENGINEERING AND ECHNOLOGY SYMOSIUM MODELING &SIMULAION, ESING AND VALIDAION (MSV) MINI-SYMOSIUM AUGUS 14-16, MICHIGAN CHAIN DYNAMIC FORMULAIONS FOR MULIBODY SYSEM RACKED
More informationMechanics Physics 151
Mechancs Physcs 151 Lecture 3 Lagrange s Equatons (Goldsten Chapter 1) Hamlton s Prncple (Chapter 2) What We Dd Last Tme! Dscussed mult-partcle systems! Internal and external forces! Laws of acton and
More informationPhysics 106a, Caltech 11 October, Lecture 4: Constraints, Virtual Work, etc. Constraints
Physcs 106a, Caltech 11 October, 2018 Lecture 4: Constrants, Vrtual Work, etc. Many, f not all, dynamcal problems we want to solve are constraned: not all of the possble 3 coordnates for M partcles (or
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 informationLagrange Multipliers. A Somewhat Silly Example. Monday, 25 September 2013
Lagrange Multplers Monday, 5 September 013 Sometmes t s convenent to use redundant coordnates, and to effect the varaton of the acton consstent wth the constrants va the method of Lagrange undetermned
More information829. An adaptive method for inertia force identification in cantilever under moving mass
89. An adaptve method for nerta force dentfcaton n cantlever under movng mass Qang Chen 1, Mnzhuo Wang, Hao Yan 3, Haonan Ye 4, Guola Yang 5 1,, 3, 4 Department of Control and System Engneerng, Nanng Unversty,
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 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 informationPoisson brackets and canonical transformations
rof O B Wrght Mechancs Notes osson brackets and canoncal transformatons osson Brackets Consder an arbtrary functon f f ( qp t) df f f f q p q p t But q p p where ( qp ) pq q df f f f p q q p t In order
More informationA particle in a state of uniform motion remain in that state of motion unless acted upon by external force.
The fundamental prncples of classcal mechancs were lad down by Galleo and Newton n the 16th and 17th centures. In 1686, Newton wrote the Prncpa where he gave us three laws of moton, one law of gravty,
More informationUIC University of Illinois at Chicago
DSL Dynamc Smulaton Laboratory UIC Unversty o Illnos at Chcago FINITE ELEMENT/ MULTIBODY SYSTEM ALGORITHMS FOR RAILROAD VEHICLE SYSTEM DYNAMICS Ahmed A. Shabana Department o Mechancal and Industral Engneerng
More informationMoments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.
Moments of Inerta Suppose a body s movng on a crcular path wth constant speed Let s consder two quanttes: the body s angular momentum L about the center of the crcle, and ts knetc energy T How are these
More informationCHAPTER 14 GENERAL PERTURBATION THEORY
CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves
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 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 informationChapter 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 informationTHE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD
Journal of Appled Mathematcs and Computatonal Mechancs 7, 6(3), 7- www.amcm.pcz.pl p-issn 99-9965 DOI:.75/jamcm.7.3. e-issn 353-588 THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS
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 informationENGI9496 Lecture Notes Multiport Models in Mechanics
ENGI9496 Moellng an Smulaton of Dynamc Systems Mechancs an Mechansms ENGI9496 Lecture Notes Multport Moels n Mechancs (New text Secton 4..3; Secton 9.1 generalzes to 3D moton) Defntons Generalze coornates
More informationInductance Calculation for Conductors of Arbitrary Shape
CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors
More informationDifference Equations
Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1
More informationThe Feynman path integral
The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space
More informationClassical Field Theory
Classcal Feld Theory Before we embark on quantzng an nteractng theory, we wll take a dverson nto classcal feld theory and classcal perturbaton theory and see how far we can get. The reader s expected to
More informationA NUMERICAL COMPARISON OF LANGRANGE AND KANE S METHODS OF AN ARM SEGMENT
Internatonal Conference Mathematcal and Computatonal ology 0 Internatonal Journal of Modern Physcs: Conference Seres Vol. 9 0 68 75 World Scentfc Publshng Company DOI: 0.4/S009450059 A NUMERICAL COMPARISON
More informationIrregular vibrations in multi-mass discrete-continuous systems torsionally deformed
(2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected
More informationPhysics 111: Mechanics Lecture 11
Physcs 111: Mechancs Lecture 11 Bn Chen NJIT Physcs Department Textbook Chapter 10: Dynamcs of Rotatonal Moton q 10.1 Torque q 10. Torque and Angular Acceleraton for a Rgd Body q 10.3 Rgd-Body Rotaton
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 informationTHE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions
THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructons by George Hardgrove Chemstry Department St. Olaf College Northfeld, MN 55057 hardgrov@lars.acc.stolaf.edu Copyrght George
More informationCOMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
More informationLie Group Formulation of Articulated Rigid Body Dynamics
Le Group Formulaton of Artculated Rgd Body Dynamcs Junggon Km 11/9/2012, Ver 1.01 Abstract It has been usual n most old-style text books for dynamcs to treat the formulas descrbng lnearor translatonal
More informationVIBRATION FATIGUE DESIGN METHODOLOGY OF A LARGE SCALE HEAVY DUTY ROBOT
ICSV14 Carns Australa 9-12 July, 2007 VIBRATION FATIGUE DESIGN METHODOLOGY OF A LARGE SCALE HEAVY DUTY ROBOT Jong Hw Seo 1, Jae Chul Hwang 1, Yong Won Cho 1, Dong Il Km 1, Hong Jae Ym 2 1 Robotcs Technology
More informationDimensional Synthesis of Wristed Binary Hands
Dmensonal Synthess of Wrsted Bnary Hands Neda Hassanzadeh Department of Mechancal Engneerng, Idaho State Unversty 921 S. 8th. Ave, Pocatello, ID 83209, USA e-mal:hassneda@su.edu Alba Perez-Graca Assocate
More informationON MECHANICS WITH VARIABLE NONCOMMUTATIVITY
ON MECHANICS WITH VARIABLE NONCOMMUTATIVITY CIPRIAN ACATRINEI Natonal Insttute of Nuclear Physcs and Engneerng P.O. Box MG-6, 07725-Bucharest, Romana E-mal: acatrne@theory.npne.ro. Receved March 6, 2008
More informationTensor Smooth Length for SPH Modelling of High Speed Impact
Tensor Smooth Length for SPH Modellng of Hgh Speed Impact Roman Cherepanov and Alexander Gerasmov Insttute of Appled mathematcs and mechancs, Tomsk State Unversty 634050, Lenna av. 36, Tomsk, Russa RCherepanov82@gmal.com,Ger@npmm.tsu.ru
More informationHidden Markov Models & The Multivariate Gaussian (10/26/04)
CS281A/Stat241A: Statstcal Learnng Theory Hdden Markov Models & The Multvarate Gaussan (10/26/04) Lecturer: Mchael I. Jordan Scrbes: Jonathan W. Hu 1 Hdden Markov Models As a bref revew, hdden Markov models
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 informationCSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography
CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve
More information12. The Hamilton-Jacobi Equation Michael Fowler
1. The Hamlton-Jacob Equaton Mchael Fowler Back to Confguraton Space We ve establshed that the acton, regarded as a functon of ts coordnate endponts and tme, satsfes ( ) ( ) S q, t / t+ H qpt,, = 0, and
More informationThe Analysis of Coriolis Effect on a Robot Manipulator
Internatonal Journal of Innovatons n Engneerng and echnology (IJIE) he Analyss of Corols Effect on a Robot Manpulator Pratap P homas Assstant Professor Department of Mechancal Engneerng K G Reddy college
More informationElshaboury SM et al.; Sch. J. Phys. Math. Stat., 2015; Vol-2; Issue-2B (Mar-May); pp
Elshabour SM et al.; Sch. J. Phs. Math. Stat. 5; Vol-; Issue-B (Mar-Ma); pp-69-75 Scholars Journal of Phscs Mathematcs Statstcs Sch. J. Phs. Math. Stat. 5; (B):69-75 Scholars Academc Scentfc Publshers
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 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 informationCOMPLEX NUMBERS AND QUADRATIC EQUATIONS
COMPLEX NUMBERS AND QUADRATIC EQUATIONS INTRODUCTION We know that x 0 for all x R e the square of a real number (whether postve, negatve or ero) s non-negatve Hence the equatons x, x, x + 7 0 etc are not
More informationEVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES
EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES Manuel J. C. Mnhoto Polytechnc Insttute of Bragança, Bragança, Portugal E-mal: mnhoto@pb.pt Paulo A. A. Perera and Jorge
More informationOn a direct solver for linear least squares problems
ISSN 2066-6594 Ann. Acad. Rom. Sc. Ser. Math. Appl. Vol. 8, No. 2/2016 On a drect solver for lnear least squares problems Constantn Popa Abstract The Null Space (NS) algorthm s a drect solver for lnear
More informationHandout # 6 (MEEN 617) Numerical Integration to Find Time Response of SDOF mechanical system. and write EOM (1) as two first-order Eqs.
Handout # 6 (MEEN 67) Numercal Integraton to Fnd Tme Response of SDOF mechancal system State Space Method The EOM for a lnear system s M X + DX + K X = F() t () t = X = X X = X = V wth ntal condtons, at
More informationSection 8.3 Polar Form of Complex Numbers
80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the
More informationLecture 14: Forces and Stresses
The Nuts and Bolts of Frst-Prncples Smulaton Lecture 14: Forces and Stresses Durham, 6th-13th December 2001 CASTEP Developers Group wth support from the ESF ψ k Network Overvew of Lecture Why bother? Theoretcal
More information9 Characteristic classes
THEODORE VORONOV DIFFERENTIAL GEOMETRY. Sprng 2009 [under constructon] 9 Characterstc classes 9.1 The frst Chern class of a lne bundle Consder a complex vector bundle E B of rank p. We shall construct
More informationConvexity preserving interpolation by splines of arbitrary degree
Computer Scence Journal of Moldova, vol.18, no.1(52), 2010 Convexty preservng nterpolaton by splnes of arbtrary degree Igor Verlan Abstract In the present paper an algorthm of C 2 nterpolaton of dscrete
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 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 informationEEE 241: Linear Systems
EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they
More informationPhysics 207: Lecture 20. Today s Agenda Homework for Monday
Physcs 207: Lecture 20 Today s Agenda Homework for Monday Recap: Systems of Partcles Center of mass Velocty and acceleraton of the center of mass Dynamcs of the center of mass Lnear Momentum Example problems
More informationReport on Image warping
Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.
More informationLagrangian Field Theory
Lagrangan Feld Theory Adam Lott PHY 391 Aprl 6, 017 1 Introducton Ths paper s a summary of Chapter of Mandl and Shaw s Quantum Feld Theory [1]. The frst thng to do s to fx the notaton. For the most part,
More informationAP Physics 1 & 2 Summer Assignment
AP Physcs 1 & 2 Summer Assgnment AP Physcs 1 requres an exceptonal profcency n algebra, trgonometry, and geometry. It was desgned by a select group of college professors and hgh school scence teachers
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 informationKernel Methods and SVMs Extension
Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general
More informationSalmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2
Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to
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 informationSupplemental document
Electronc Supplementary Materal (ESI) for Physcal Chemstry Chemcal Physcs. Ths journal s the Owner Socetes 01 Supplemental document Behnam Nkoobakht School of Chemstry, The Unversty of Sydney, Sydney,
More informationSo far: simple (planar) geometries
Physcs 06 ecture 5 Torque and Angular Momentum as Vectors SJ 7thEd.: Chap. to 3 Rotatonal quanttes as vectors Cross product Torque epressed as a vector Angular momentum defned Angular momentum as a vector
More information