8. THE CONTACT DYNAMICS METHOD

Size: px
Start display at page:

Download "8. THE CONTACT DYNAMICS METHOD"

Transcription

1 8 THE CONTACT DYNAMICS METHOD 8 Introducton The Contact Dynamcs ( CD ) method was ntroduced at the begnnng of the 990es by M Jean and JJ Moreau (Jean és Moreau, 992; Jean, 999) The method turned out to be extremely fast when smulatng granular flows, rad avalanches, segregaton, vbraton roblems of granular materals etc Most dscrete element methods take nto consderaton the deformablty of the elements eg by concentratng t nto the contacts, lke n the case of BALL-tye models, or by alyng an nternal fnte element mesh lke n UDEC The calculaton of the contact forces between the elements are based on the stffness characterstcs The logc of CD s dfferent The elements are consdered to be erfectly rgd, and the contact forces are not related to any stffness data The contact forces should ensure the eulbrum of the elements, and n addton, they must not volate reurements lke eg Coulomb lmt for frcton, or no tenson n cohesonless contacts, but ther calculaton does not aly any consttutve relatons For statcally hghly ndetermnate systems lke eg granular assembles n general, there exst several vald force systems satsfyng the eulbrum condtons; CD roduces randomly one of them, and several eually vald solutons can be receved f the roblem s calculated agan and agan wth the Contact Dynamcs Method No wonder that ts alcaton s not very wdesread n the engneerng ractce; on the other hand, snce the method s comutatonally very effcent for the smulaton of dynamc roblems, t s rather oular n the granular hyscs lterature The orgnal aers on CD were rather dffcult to decher The aer of Unger and Kertész (2003) was very dfferent from them n ths resect: the authors gave a clear, codewrter-orented ntroducton to the man lne of thought of the method, gvng a sgnfcant hel ths way for those who wanted to wrte ther own code and also for those who just wanted to understand how the method worked whch they were alyng n ther researches The ntroducton below s also based on ths ublcaton 82 The elements and the contacts Unlke n other methods where the basc unt of the analyss s the element or a art of the element around the node, n CD the basc unt s the ar of neghbourng elements The am of the calculatons s to teratvely fnd the sutable contact forces n all ars n the system The elements n the orgnal verson of the CD method were crcular and erfectly rgd Later several ublcatons aeared n the lterature on CD models wth non-crcular and even on deformable elements; but for smlcty, the ntroducton below wll focus only on erfectly rgd, erfectly shercal three-dmensonal elements These elements have sx degrees-of-freedom (3 translatons and 3 rotatons) A reference ont s defned for each element concdng wth ts centre of gravty Focus now on a selected ar of two elements close to (erhas n contact wth) each other We shall collect the state varables and the mechancal characterstcs of ths ar 77

2 Let g denote the dstance between the elements and If ths dstance euals to zero, the two elements get nto a ont-lke contact (denote t by c), and a concentrated contact force can be transmtted between the elements Unlke n most dscrete element technues, now the contact and the elements do not deform, so contactng elements must move n such a way that there would be no overla between them: g cannot become negatve g Fgure 4 The dstance g between the elements and The force actng on by s f c, and the force exressed on by s f c : f c F F F c x c y c z () t () t () t and f c F F F c x c y c z () t () t () t, and these can be reduced to the reference onts of and wth the hel of the B c (t) and B c (t) transton matrces, resectvely: B c c ( t) c c ; B ( t) c c 0 rz ry 0 rz ry c c c c rz 0 r x rz 0 r x c c c c ry rx 0 ry rx 0 The vectors r c and r c ont from the reference onts to the contact The reduced forces are then: c c c c c c f ( t) B ( t) f ( t) ; f ( t) B ( t) f ( t) red Summarze these two forces nto the 2-scalar reduced contact force vector of the ar:: red f f () t f c red c red () t () t Ths vector wll take art n the euatons of moton of the ar 78

3 n c Fgure 5 The normal vector of the contact Assgn an (n, t, w) local coordnate frame to the contact: n s the common unt normal of the contact lane, ontng towards element (see Fgure 5); t and w are an arbtrary ar of orthogonal unt vectors n the contact lane The vector v denotes the velocty vector of the (, ) ar: v v () t v () t () t, whch can be exressed as the tme dervatves of the dslacements: v dux ( t) dux( t) duy( t) duy( t) duz ( t) duz( t) ( t) ; v ( t) dx ( t) dx( t) dy ( t) dy( t) dz ( t) dz ( t) Consder now those two materal onts, c and c, whch form the contact Ther veloctes, v c és v c, can be exressed agan wth the hel of the B c (t) and B c (t) transton matrces: c ct v ( t) B ( t) v ( t) c ct v ( t) B ( t) v ( t) The relatve velocty of the contact s the dfference of the veloctes of the two materal onts c and c: 79

4 whch yelds x () t c c μ ( t) y ( t) : v ( t) v ( t) z () t c c ct ct μ ( t) v ( t) v ( t) B ( t) v ( t) B ( t) v ( t) Ths vector shows the relatve translaton of ont c wth resect to ont c, so f the two onts are already n contact, the n-comonent of ths vector must be ostve to exclude overlang Contacts n the CD method are modelled n the followng way: () If the two elements are n contact, e f g = 0, then a contact force can be transmtted If g > 0, then there s no contact, and no contact force exst n the ar The case g < 0 s not ossble (2) The normal comonent of the contact force, N, can only be negatve, e ct ct comressonal, but otherwse ts magntude ( N f ( n) vagy N f n) s arbtrary The (N, g ) relaton s shown n Fgure 6: N g Fgure 6 Normal force versus ga wh n CD (3) The magntude of the tangental T (whch has a t and a w comonent) s lmted by the Coulomb frcton law: T N where s the frcton coeffcent Fgure 7 llustrates ths lmtaton: the vector of the tangental force must ont from the orgn ether to nsde the cone (non-sldng contact) or just to the surface of the cone (sldng contact) 80

5 N arctan T t T w Fgure 7 Columb lmt for the tangental force As long as T s below the frcton lmt, the tangental comonent of the relatve translaton must be zero (elastc deformaton not ossble) When reachng the frcton lmt, the ontact starts to slde: the μ () t relatve velocty has to be just ooste to the drecton of the T force (the tangental relatve translaton must not have a comonent erendcularly to the frctonal force reachng the lmt), and the magntude of T s eual to N The magntude of the relatve translaton s not lmted n the contact model, ths can be determned from knematcal consderatons 83 The euatons of moton CD s a tme-steng method Its fundamental unknowns are the tme-deendent ostons and veloctes of the elements; however, they are not comled nto hyervectors lke n other dscrete element methods, because they are analysed searately ar by ar Assume that at t the state of the system s known: as far as numercally ossble, the ostons and veloctes of the elements are gven: u ( t ) u ; v ( t ) v, and the external, ext forces actng on the elements ( f ) and contact forces for all c ( f c ) are also known The external forces are reduced to the reference onts; the contact forces act n the ont-lke contacts The tme-deendence of the external forces s also known (eg that the gravtatonal force s constant), so the external forces are gven also n t t t (, ext f ) From these data the state of the system at t + (contact forces, the ostons and the veloctes of the elements) s searched for CD ales the mlct verson of the Euler method for ths The basc ste for the and ar can be wrtten as: v v ( M ) f : t v v ( M ) f u u v : t u u v ; 8

6 Here f and f denote the resultants of the external and all contact forces actng on and resectvely, beng reduced to the reference onts (Note that accordng to the mlct scheme, the veloctes and acceleratons belongng to the end of the tme nterval are consdered to be vald along the whole nterval) These reduced forces are, külső k k f B f f ( k ) :, külső k k f f B f ( k ) Summaton over ndex k runs along all contacts of element, ncludng the just analysed contact c as well Smlarly, ndex k runs along all contacts of The transton matrces are assumed to be constant durng the (t, t + ) tme nterval, and eual to ther values at t Indeed, f the dslacement ncrements are small durng the tmeste, the modfcaton of the vectors ontng from the reference onts to the contacts s neglgble Collect the mass and rotatonal nerta of the elements nto the matrces M and M, whch have the followng form n the case of shercal elements: m m m M I I I (Note that because of the shercal symmetry, these matrces are constant n tme) In order to determne the velocty of the ar ( v ), the resultants f and f should be known So, n addton to the external forces actng at t +, the contact forces should also be known at the end of the tmeste The Contact Dynamcs models search for these contact forces wth the hel of an teratve solver (whch has to be erformed over and over agan at every tmeste, as the contact forces change wth tme) 84 The teratve solver The solver swees along all ars of neghbourng or nearly contactng elements When consderng a gven ar, an aroxmaton s gven (based on the euatons of moton of the ar), so that the condtons assumed on the mechancal behavour would be satsfed: no overla; Coulomb-frcton etc Then a next ar s consdered After all ars were swet over, the solver starts the ars from the begnnng, and t s reeated over and over agan, untl the next aroxmatons are already suffcently close to the revous ones It means that the contact forces belongng to t + have been found, and the next tme ste can follow 82

7 The aroxmaton of the contact force n the ar (, ) s based on the euatons of moton of that ar Before turnng onto the detals, a few notatons have to be ntroduced: For element, reduce to the reference ont all those forces (external and contact forces) actng at t +, excet from the force exressed by element through contact c:, no _ c, ext k fred, : f f red, k c The f k contact force s only an actual aroxmaton of the force ndeed actng n contact k at t + ; t receves new and new values durng the teratons (At the begnnng of the analyss of the tme ste the contact forces are aroxmated to be the same as ther fnal, just determned values at the end of the revous tme ste, whch s the same as the begnnng of the just analysed tmeste) Smlarly, reduce all the forces actng at t + on excet from that force actng n c to the reference ont of :, no _ c, ext k fred, : f f red,, k c and collect the two vectors nto a hyervector:, no _ c, no _ c f red, fred, :, no _ c fred, Summarze the two transton matrces belongng to c nto a hyermatrx: c B B : c B and the matrces of nerta of and nto a block-dagonal matrx, whose nverse s: M 0 M : 0 M Later the followng two matrces wll also be necessary: and : T M B M B, M : B M B T And now the euatons of moton of the ar (, ) can be comled Frst, the euatons belongng to the end of the tme nterval can searately be wrtten as:, no _ c c v v f red, f red, M, no _ c c t v v fred, fred, Multly both sdes by B T from the left: T, no _ c B M fred, t T c B M B f c c (t was taken nto consderaton that f f ) After some rearrangements: 83

8 , _ t t T no c c B M fred, M f It s easy to notce that on the left sde the vector n the arentheses means that relatve velocty whch would occur n the contact at t + f f c s zero, e f there s no force n the contact Ths vector wll have a secal mortance n the forthcomng dervaton, so a secal notaton s gven to t:, no _ c, _ : T no c t B M fred, The euatons of moton can now be wrtten as: t M f, no _ c c where and f c are the unknowns So the euatons of moton gve the relaton between the unknown contact force and the unknown relatve velocty belongng to the contact Ths wll be the startng ont of the forthcomng calculatons Fnally the normal and tangental comonents of the relatve velocty vector of the contact wll be needed: T n ; n (Remember that n tw n denoted the velocty of the materal ont c relatve to the materal ont c Hence a ostve n means ncreasng ga between the two materal onts) Snce, no _ c the vector belongng to the tme nstant t + can drectly be calculated from the already, no _ c exstng aroxmatons of all other contact forces excet from c, the comonents and, no _ c tw, can also be determned, whle the comonents of the vector n, are unknowns The unknown and f c vectors are determned n three stes: Ste Frst decde whether the two elements wll be n contact at t + : calculate how large wll the ga be between them, assumng zero contact force:, no _ c, no _ c g g t n, A ostve result means that there wll be no contact at t +, and the analyss of another ar can mmedately follow A negatve result, on the other hand, means that wthout a contact force the elements and would overla, so an f c contact force s needed to avod the overla In ths case Ste 2 follows Ste 2 The contact force should modfy the veloctes of the two elements n such a way that nstead of overlang, they would exactly touch each other at the end of the tme ste In Ste 2 the am s to determne f c that satsfes the followng two condtons: () at t + the gawh between and s exactly zero: g t n, 0 () the contact does not slde, so the tangental comonent of the relatve translaton s zero: tw, 0 To satsfy these two condtons, the relatve velocty of the contact should be: 84

9 g n t (the ngatve sgn means that f the gawh was larger than zero, then should get closer to to touch t) The f c has to be such a force that f contnuously actng between and durng (t, t + ), at t + the relatve velocty would be just eual to μ From the euatons of moton, ths force turns out to be eual to: c, no _ c f M g n t t Now the ueston s whether ths force volates the consttutve condtons There were two condtons on the comonents of the contact forces The frst one reured the normal force a comresson Ths s automatcally satsfed because of Ste The second one was the Coulomb-condton: T c N c If ths holds for the calculated f c, then the analyss of the (, ) ar s ready, and a next ar can follow However, f the tangental comonent exceeds the frcton lmt, then the calculated contact force cannot be transmtted n the contact: the contact sldes, whch means that the tangental comonent of s not zero, and the calculaton based on zero tangental comonent should be corrected Ths correcton s done n Ste 3 Ste 3 In a sldng contact the tangental force comonent has to satsfy the followng to condtons, and as the thrd condton the euatons of moton: () The contact s sldng, so the magntude of the tangental force comonent s eual to the Coulomb-lmt: T c N c () The drecton of the tangental relatve velocty s just ooste to the drecton of the tangental force comonent: c T tw, c T tw, () the euatons of moton: c, no _ c f M g n tw, t t From these condtons the unknowns f c és can be calculated, and the analyss of the (, ) ar s ready The next ar can follow These calculatons ntroduced above gve an aroxmaton for the contact force n a ar, assumng that all other contact forces are unchanged and kee ther values last aroxmated When turnng to the next ar, the latest aroxmatons n other ars are aled Proceedng from ar to ar ths way, an aroxmaton s receved for the whole system of contact forces By sweeng through the comlete set of contacts over and over agan, the results (at least, hoefully) get closer and closer to what should really exst at t + The modfcatons caused by the consecutve teraton cycles cause smaller and smaller modfcatons n the contact forces; and the teraton can be termnated as the modfcatons 85

10 decrease under a rescrbed threshold Now the state belongng to t + has been found, and a new tme ste can be analysed The convergence of the method s an oen ueston n the lterature The order accordng to whch the ars are consdered wthn an teraton ste s random; the only reurement s that every ar should be consdered once In the next teraton stes the orderng s dfferent, rescrbed also by a random number generator If the same roblem s analysed twce, by startng the random number generator from two dfferent laces, the two resultng contact force systems wll be dfferent Ths s evdent for an engneer, artcularly f eulbrum s searched for: for a statcally hghly ntedermnate system several eulbrated force systems can be found, and wthout flexblty data the correct one cannot be selected There are several ublcatons dealng wth the ndetermnacy of the results Exerences show that though the order of the ars greatly affect the ndvdual contact forces and even the toology of the system, the overall, macro characterstcs lke average stress tensor or freuency dagram of contact force magntudes reman the same, aart from slght statstcal devatons 85 Alcatons The Contact Dynamcs Method s rather oular among hyscsts studyng granular dynamcs roblems (eg Daudon et al, 997; Radja et al, 998; Unger et al, 2004), but not aled yet n the engneerng ractce, artly because of the doubts on ts relablty and the random nature of the results Questons 8 Wrte the euatons of moton of ar n CD, and exlan the meanng of the uanttes n t! 82 Descrbe the mechancs of the contacts n CD! 83 Exlan the analyss of a tme ste n CD! 86

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Week3, 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 information

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure 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 information

Week 9 Chapter 10 Section 1-5

Week 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 information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

Translational Equations of Motion for A Body Translational equations of motion (centroidal) for a body are m r = f.

Translational 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 information

CHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics)

CHAPTER 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 information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 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 information

One Dimensional Axial Deformations

One 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 information

Spin-rotation coupling of the angularly accelerated rigid body

Spin-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 information

Week 11: Chapter 11. The Vector Product. The Vector Product Defined. The Vector Product and Torque. More About the Vector Product

Week 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 information

Finite Element Modelling of truss/cable structures

Finite 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 information

Indeterminate pin-jointed frames (trusses)

Indeterminate 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

χ 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 information

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

ELASTIC 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 information

Difference Equations

Difference 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 information

The classical spin-rotation coupling

The 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 information

PHYS 705: Classical Mechanics. Newtonian Mechanics

PHYS 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 information

Physics 53. Rotational Motion 3. Sir, I have found you an argument, but I am not obliged to find you an understanding.

Physics 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 information

Lecture 12: Discrete Laplacian

Lecture 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 information

Moments 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 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 information

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

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 information

Physics 207 Lecture 6

Physics 207 Lecture 6 Physcs 207 Lecture 6 Agenda: Physcs 207, Lecture 6, Sept. 25 Chapter 4 Frames of reference Chapter 5 ewton s Law Mass Inerta s (contact and non-contact) Frcton (a external force that opposes moton) Free

More information

Advanced Topics in Optimization. Piecewise Linear Approximation of a Nonlinear Function

Advanced Topics in Optimization. Piecewise Linear Approximation of a Nonlinear Function Advanced Tocs n Otmzaton Pecewse Lnear Aroxmaton of a Nonlnear Functon Otmzaton Methods: M8L Introducton and Objectves Introducton There exsts no general algorthm for nonlnear rogrammng due to ts rregular

More information

Modeling of Dynamic Systems

Modeling 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 information

In this section is given an overview of the common elasticity models.

In 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 information

Chapter 8. Potential Energy and Conservation of Energy

Chapter 8. Potential Energy and Conservation of Energy Chapter 8 Potental Energy and Conservaton of Energy In ths chapter we wll ntroduce the followng concepts: Potental Energy Conservatve and non-conservatve forces Mechancal Energy Conservaton of Mechancal

More information

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Physics 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 information

ENGN 40 Dynamics and Vibrations Homework # 7 Due: Friday, April 15

ENGN 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 information

Iterative General Dynamic Model for Serial-Link Manipulators

Iterative 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 information

First 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.

First 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 information

MMA and GCMMA two methods for nonlinear optimization

MMA and GCMMA two methods for nonlinear optimization MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons

More information

( ) 2 ( ) ( ) Problem Set 4 Suggested Solutions. Problem 1

( ) 2 ( ) ( ) Problem Set 4 Suggested Solutions. Problem 1 Problem Set 4 Suggested Solutons Problem (A) The market demand functon s the soluton to the followng utlty-maxmzaton roblem (UMP): The Lagrangean: ( x, x, x ) = + max U x, x, x x x x st.. x + x + x y x,

More information

Lecture Notes on Linear Regression

Lecture Notes on Linear Regression Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume

More information

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

THE 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 information

coordinates. Then, the position vectors are described by

coordinates. 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 information

6. Hamilton s Equations

6. Hamilton s Equations 6. Hamlton s Equatons Mchael Fowler A Dynamcal System s Path n Confguraton Sace and n State Sace The story so far: For a mechancal system wth n degrees of freedom, the satal confguraton at some nstant

More information

Physics 181. Particle Systems

Physics 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

PY2101 Classical Mechanics Dr. Síle Nic Chormaic, Room 215 D Kane Bldg

PY2101 Classical Mechanics Dr. Síle Nic Chormaic, Room 215 D Kane Bldg PY2101 Classcal Mechancs Dr. Síle Nc Chormac, Room 215 D Kane Bldg s.ncchormac@ucc.e Lectures stll some ssues to resolve. Slots shared between PY2101 and PY2104. Hope to have t fnalsed by tomorrow. Mondays

More information

Chapter - 2. Distribution System Power Flow Analysis

Chapter - 2. Distribution System Power Flow Analysis Chapter - 2 Dstrbuton System Power Flow Analyss CHAPTER - 2 Radal Dstrbuton System Load Flow 2.1 Introducton Load flow s an mportant tool [66] for analyzng electrcal power system network performance. Load

More information

The Feynman path integral

The 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 information

NMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582

NMT 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 information

EN40: Dynamics and Vibrations. Homework 7: Rigid Body Kinematics

EN40: 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 information

12. The Hamilton-Jacobi Equation Michael Fowler

12. 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 information

Mathematical Preparations

Mathematical 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 information

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

Irregular 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 information

A particle in a state of uniform motion remain in that state of motion unless acted upon by external force.

A 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 information

SCALARS AND VECTORS All physical quantities in engineering mechanics are measured using either scalars or vectors.

SCALARS 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 information

Lagrange Multipliers. A Somewhat Silly Example. Monday, 25 September 2013

Lagrange 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 information

Physics 111: Mechanics Lecture 11

Physics 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 information

So far: simple (planar) geometries

So 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

PHYS 705: Classical Mechanics. Calculus of Variations II

PHYS 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 information

CHAPTER 14 GENERAL PERTURBATION THEORY

CHAPTER 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 information

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

COMPOSITE 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 information

Managing Capacity Through Reward Programs. on-line companion page. Byung-Do Kim Seoul National University College of Business Administration

Managing Capacity Through Reward Programs. on-line companion page. Byung-Do Kim Seoul National University College of Business Administration Managng Caacty Through eward Programs on-lne comanon age Byung-Do Km Seoul Natonal Unversty College of Busness Admnstraton Mengze Sh Unversty of Toronto otman School of Management Toronto ON M5S E6 Canada

More information

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

The equation of motion of a dynamical system is given by a set of differential equations. That is (1) Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence

More information

Rigid body simulation

Rigid 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 information

A PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.

A PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY. Proceedngs of the th Brazlan Congress of Thermal Scences and Engneerng -- ENCIT 006 Braz. Soc. of Mechancal Scences and Engneerng -- ABCM, Curtba, Brazl,- Dec. 5-8, 006 A PROCEDURE FOR SIMULATING THE NONLINEAR

More information

Module 3: Element Properties Lecture 1: Natural Coordinates

Module 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 information

C/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/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 information

DUE: WEDS FEB 21ST 2018

DUE: 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 information

Chapter 11 Angular Momentum

Chapter 11 Angular Momentum Chapter 11 Angular Momentum Analyss Model: Nonsolated System (Angular Momentum) Angular Momentum of a Rotatng Rgd Object Analyss Model: Isolated System (Angular Momentum) Angular Momentum of a Partcle

More information

11. Dynamics in Rotating Frames of Reference

11. 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 information

Chapter Eight. Review and Summary. Two methods in solid mechanics ---- vectorial methods and energy methods or variational methods

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 information

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

Salmon: 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 information

Please review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam.

Please review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam. Please revew the followng statement: I certfy that I have not gven unauthorzed ad nor have I receved ad n the completon of ths exam. Sgnature: Instructor s Name and Secton: (Crcle Your Secton) Sectons:

More information

Lecture Note 3. Eshelby s Inclusion II

Lecture 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 information

Part C Dynamics and Statics of Rigid Body. Chapter 5 Rotation of a Rigid Body About a Fixed Axis

Part 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 information

Confidence intervals for weighted polynomial calibrations

Confidence intervals for weighted polynomial calibrations Confdence ntervals for weghted olynomal calbratons Sergey Maltsev, Amersand Ltd., Moscow, Russa; ur Kalambet, Amersand Internatonal, Inc., Beachwood, OH e-mal: kalambet@amersand-ntl.com htt://www.chromandsec.com

More information

1 Matrix representations of canonical matrices

1 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

Work is the change in energy of a system (neglecting heat transfer). To examine what could

Work is the change in energy of a system (neglecting heat transfer). To examine what could Work Work s the change n energy o a system (neglectng heat transer). To eamne what could cause work, let s look at the dmensons o energy: L ML E M L F L so T T dmensonally energy s equal to a orce tmes

More information

Physics 207: Lecture 20. Today s Agenda Homework for Monday

Physics 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 information

763622S ADVANCED QUANTUM MECHANICS Solution Set 1 Spring c n a n. c n 2 = 1.

763622S ADVANCED QUANTUM MECHANICS Solution Set 1 Spring c n a n. c n 2 = 1. 7636S ADVANCED QUANTUM MECHANICS Soluton Set 1 Sprng 013 1 Warm-up Show that the egenvalues of a Hermtan operator  are real and that the egenkets correspondng to dfferent egenvalues are orthogonal (b)

More information

A Quadratic Cumulative Production Model for the Material Balance of Abnormally-Pressured Gas Reservoirs F.E. Gonzalez M.S.

A Quadratic Cumulative Production Model for the Material Balance of Abnormally-Pressured Gas Reservoirs F.E. Gonzalez M.S. Natural as Engneerng A Quadratc Cumulatve Producton Model for the Materal Balance of Abnormally-Pressured as Reservors F.E. onale M.S. Thess (2003) T.A. Blasngame, Texas A&M U. Deartment of Petroleum Engneerng

More information

Time-Varying Systems and Computations Lecture 6

Time-Varying Systems and Computations Lecture 6 Tme-Varyng Systems and Computatons Lecture 6 Klaus Depold 14. Januar 2014 The Kalman Flter The Kalman estmaton flter attempts to estmate the actual state of an unknown dscrete dynamcal system, gven nosy

More information

Non-Ideality Through Fugacity and Activity

Non-Ideality Through Fugacity and Activity Non-Idealty Through Fugacty and Actvty S. Patel Deartment of Chemstry and Bochemstry, Unversty of Delaware, Newark, Delaware 19716, USA Corresondng author. E-mal: saatel@udel.edu 1 I. FUGACITY In ths dscusson,

More information

From Biot-Savart Law to Divergence of B (1)

From Biot-Savart Law to Divergence of B (1) From Bot-Savart Law to Dvergence of B (1) Let s prove that Bot-Savart gves us B (r ) = 0 for an arbtrary current densty. Frst take the dvergence of both sdes of Bot-Savart. The dervatve s wth respect to

More information

Tensor Analysis. For orthogonal curvilinear coordinates, ˆ ˆ (98) Expanding the derivative, we have, ˆ. h q. . h q h q

Tensor 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 information

Chapter Newton s Method

Chapter Newton s Method Chapter 9. Newton s Method After readng ths chapter, you should be able to:. Understand how Newton s method s dfferent from the Golden Secton Search method. Understand how Newton s method works 3. Solve

More information

829. An adaptive method for inertia force identification in cantilever under moving mass

829. 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 information

Foundations of Arithmetic

Foundations of Arithmetic Foundatons of Arthmetc Notaton We shall denote the sum and product of numbers n the usual notaton as a 2 + a 2 + a 3 + + a = a, a 1 a 2 a 3 a = a The notaton a b means a dvdes b,.e. ac = b where c s an

More information

LAB # 4 - Torque. d (1)

LAB # 4 - Torque. d (1) LAB # 4 - Torque. Introducton Through the use of Newton's three laws of moton, t s possble (n prncple, f not n fact) to predct the moton of any set of partcles. That s, n order to descrbe the moton of

More information

A Quadratic Cumulative Production Model for the Material Balance of Abnormally-Pressured Gas Reservoirs F.E. Gonzalez M.S.

A Quadratic Cumulative Production Model for the Material Balance of Abnormally-Pressured Gas Reservoirs F.E. Gonzalez M.S. Formaton Evaluaton and the Analyss of Reservor Performance A Quadratc Cumulatve Producton Model for the Materal Balance of Abnormally-Pressured as Reservors F.E. onale M.S. Thess (2003) T.A. Blasngame,

More information

Please review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam.

Please review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam. ME 270 Summer 2014 Fnal Exam NAME (Last, Frst): Please revew the followng statement: I certfy that I have not gven unauthorzed ad nor have I receved ad n the completon of ths exam. Sgnature: INSTRUCTIONS

More information

Uncertainty and auto-correlation in. Measurement

Uncertainty and auto-correlation in. Measurement Uncertanty and auto-correlaton n arxv:1707.03276v2 [physcs.data-an] 30 Dec 2017 Measurement Markus Schebl Federal Offce of Metrology and Surveyng (BEV), 1160 Venna, Austra E-mal: markus.schebl@bev.gv.at

More information

Implicit Integration Henyey Method

Implicit Integration Henyey Method Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure

More information

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Lectures - 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 information

Affine transformations and convexity

Affine transformations and convexity Affne transformatons and convexty The purpose of ths document s to prove some basc propertes of affne transformatons nvolvng convex sets. Here are a few onlne references for background nformaton: http://math.ucr.edu/

More information

Lecture 5 Decoding Binary BCH Codes

Lecture 5 Decoding Binary BCH Codes Lecture 5 Decodng Bnary BCH Codes In ths class, we wll ntroduce dfferent methods for decodng BCH codes 51 Decodng the [15, 7, 5] 2 -BCH Code Consder the [15, 7, 5] 2 -code C we ntroduced n the last lecture

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

MEV442 Introduction to Robotics Module 2. Dr. Santhakumar Mohan Assistant Professor Mechanical Engineering National Institute of Technology Calicut

MEV442 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 information

Problem Set 9 Solutions

Problem Set 9 Solutions Desgn and Analyss of Algorthms May 4, 2015 Massachusetts Insttute of Technology 6.046J/18.410J Profs. Erk Demane, Srn Devadas, and Nancy Lynch Problem Set 9 Solutons Problem Set 9 Solutons Ths problem

More information

An Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors

An 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 information

Mechanics Physics 151

Mechanics Physics 151 Mechancs hyscs 151 Lecture Canoncal Transformatons (Chater 9) What We Dd Last Tme Drect Condtons Q j Q j = = j, Q, j, Q, Necessary and suffcent j j for Canoncal Transf. = = j Q, Q, j Q, Q, Infntesmal CT

More information

Second Order Analysis

Second 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 information

Chapter 3. r r. Position, Velocity, and Acceleration Revisited

Chapter 3. r r. Position, Velocity, and Acceleration Revisited Chapter 3 Poston, Velocty, and Acceleraton Revsted The poston vector of a partcle s a vector drawn from the orgn to the locaton of the partcle. In two dmensons: r = x ˆ+ yj ˆ (1) The dsplacement vector

More information

2. PROBLEM STATEMENT AND SOLUTION STRATEGIES. L q. Suppose that we have a structure with known geometry (b, h, and L) and material properties (EA).

2. PROBLEM STATEMENT AND SOLUTION STRATEGIES. L q. Suppose that we have a structure with known geometry (b, h, and L) and material properties (EA). . PROBEM STATEMENT AND SOUTION STRATEGIES Problem statement P, Q h ρ ρ o EA, N b b Suppose that we have a structure wth known geometry (b, h, and ) and materal propertes (EA). Gven load (P), determne the

More information

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Transfer 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 information

SCHOOL OF COMPUTING, ENGINEERING AND MATHEMATICS SEMESTER 2 EXAMINATIONS 2011/2012 DYNAMICS ME247 DR. N.D.D. MICHÉ

SCHOOL OF COMPUTING, ENGINEERING AND MATHEMATICS SEMESTER 2 EXAMINATIONS 2011/2012 DYNAMICS ME247 DR. N.D.D. MICHÉ s SCHOOL OF COMPUTING, ENGINEERING ND MTHEMTICS SEMESTER EXMINTIONS 011/01 DYNMICS ME47 DR. N.D.D. MICHÉ Tme allowed: THREE hours nswer: ny FOUR from SIX questons Each queston carres 5 marks Ths s a CLOSED-BOOK

More information

Inductance Calculation for Conductors of Arbitrary Shape

Inductance 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 information

More metrics on cartesian products

More metrics on cartesian products More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of

More information