Lecture Note 3. Eshelby s Inclusion II
|
|
- Edwin Casey
- 6 years ago
- Views:
Transcription
1 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 sold Incluson energy by work method 5 3 Incluson energy n a fnte sold 7 4 Colonett s theorem 9 5 Fnte sold wth appled tractons 1 1
2 1 Incluson energy n an nfnte sold So far we have obtaned the expressons for the stress, stran and dsplacement feld both nsde and outsde the ncluson. An mportant queston s: what s the total elastc energy E of the sold contanng an ncluson? In ths and subsequent sectons, we derve the expressons for E, whch we refer to as the ncluson energy for brevty. However, we emphasze that E s the total elastc energy of the sold contanng an ncluson. E ncludes the elastc energy stored both nsde and outsde the ncluson. For example, f we obtan E as a functon of the ncluson sze, then the dervatve of E provdes the drvng force for the expanson (or shrnkage) of the ncluson. Notce that ths s the case only f E s the total elastc energy, not just the energy stored nsde the ncluson. There are two ways to obtan the expresson for the total energy E. Frst, we can ntegrate the elastc energy densty both nsde and outsde the ncluson, usng the feld expressons we have already obtaned. Second, we can obtan the elastc energy E by measurng the work done n a vrtual experment that transforms a sold system wth zero elastc energy to the sold contanng an ncluson. In ths secton, we take the frst approach. The work method s dscussed n the next secton, whch leads to dentcal results but may provde more physcal nsght. For clarty, let us ntroduce some symbols to descrbe the elastc felds nsde and outsde the ncluson. Let the elastc (stress, stran, dsplacement) felds nsde the ncluson be denoted by a superscrpt I, and the elastc felds outsde the ncluson (n the matrx) be denoted by a superscrpt M. Notce that whenever the superscrpt I or M s used, the felds only nclude the elastc component. For a homogeneous nfnte sold, the elastc felds n the matrx and the ncluson are, matrx ncluson e M = e c e I = e c e M = c I = c u M = u c u I = u c e x j Therefore, the total elastc energy s, E = 1 e I I d + 1 M e M d (1) Rewrtng E n terms of dsplacements, we have E = 1 4 (u I I,j + u I j,) d + 1 M (u M,j + u M 4 j,) d () and notng the symmetry of the stress tensor E = 1 u I I j, d + 1 M u M j, d (3) Now, the dervatve can be factored out usng the followng rule u,j = ( u j ),, u j (4)
3 E = 1 ( I u I j), I,u I j d + 1 ( M u M j ), M,u M j d (5) The body s assumed not to have any body forces actng on t, thus the dvergence of the stress tensor,,, s zero. Thus E = 1 ( u I I j), d + 1 ( M u M j ), d (6) We wsh to now use Gauss s theorem on ths equaton. We need to careful about the sgn of the unt normal vector that ponts outsde the ntegraton volume. Let the normal vector pontng out of the ncluson volume be n out. Let the unt normal vector pontng out of the outer surface of the matrx (at nfnty) be n. Applyng Gauss s theorem, E = 1 u I I jn out ds 1 M u M j n out ds + 1 M u M j n ds (7) S We expect that the surface ntegral over S should vansh as t approaches nfnty. To show ths, let S be a sphercal surface whose radus R approaches nfnty. Notce that D kl (x) = G,kl (x x )d (x ) (8) Because G kl (x x ) R 3 where R = x, for large R, then D kl (x) R 3. Therefore, ( ) 1 e M = O (9) R ( 3 ) 1 M = O (10) R 3 ds = O ( R ) (11) Thus S M u M j n out d 0 as R (1) Combnng the two ntegrals over, E = 1 ( ) I u I j M u M j n out ds (13) Although the stress across the ncluson nterface does not have to be contnuous, the tracton force across the nterface must be contnuous,.e., I n out whch leads to E = 1 = M n out I ( u I j u M j (14) ) n out ds (15) 3
4 From the defnton of (elastc dsplacement felds) u I j and u M j, we have Thus u I j u M j = (u c j e jkx k ) u c j = e jkx k (16) E = 1 I n out e jkx k ds (17) Therefore, we have expressed the total elastc energy E n terms of a surface ntegral over, the ncluson nterface. We can further smplfy ths expresson by transformng the ntegral back nto a volume ntegral (over the ncluson volume ). E = 1 ( ) I e jkx k d, = 1 ( ) e jk I, x k + x I k, d = 1 e jkkj I d = 1 e I d = 1 ( ) e c d (18) For an ellpsodal ncluson, the stress nsde s a constant, thus E = 1 ( ) c e = 1 I e (19) If the volume s not an ellpsod, we can stll wrte the energy n terms of the average stress n the ncluson E = 1 I e (0) where I 1 c (x) d (x) (1) Suppose that we wsh to account for how much of the energy s stored nsde the ncluson and how much s stored n the matrx. The energy store nsde the ncluson s E I = 1 e I I d For ellpsodal ncluson, the stress and stran are constant nsde, hence E I = 1 I e I = 1 I ( ) e c e 0 4
5 Snce the total elastc energy s E = 1 I e the elastc energy stored nsde the matrx must be, E M = E E I = 1 I e c Incluson energy by work method In ths secton, we re-derve the expressons n the prevous secton concernng the ncluson energy usng a dfferent approach. Rather than ntegratng the stran energy densty over the entre volume, we make use of the fact that the stored elastc (potental) energy n the sold must equal to the work done to t. By consderng a vrtual experment that transforms a stress-free sold nto a sold contanng an ncluson and accountng for the work done along the way, we can derve the total elastc energy (or the elastc energy stored wthn the ncluson or the matrx) usng consderably less math than before. To better llustrate ths method, let us consder a smple example. Consder a mass M attached to a lnear sprng wth stffness k. Let E 0 be the equlbrum state of the system under no appled force. Obvously E 0 = 0. Defne the orgn as the poston of the mass at ths state. Suppose we gradually apply a force to the mass untl the force reaches F 1. At ths pont the mass must have moved by a dstance x 1 = F 1 /k. Let the energy of ths state be E 1. The work done n movng the mass from 0 to x 1 equals to the average force F appled to the mass tmes the dstance travelled (x 1 ). Because the ntal force s 0 and the fnal force s F 1, the average force s F = F 1 /. Therefore, the work done n movng the mass from 0 to x 1 s, Hence W 01 = F x 1 = 1 F 1x 1 = 1 kx 1 () E 1 = E 0 + W 01 = 1 kx 1 (3) Suppose we further ncrease the force to F and the system reaches a new state at x = F /k wth energy E. Snce the ntal force durng ths transformaton s F 1 and the fnal force s F, the average force s F = (F 1 + F )/. The mass moves by a dstance of x x 1 under ths force. Therefore the work done s 5
6 Hence W 1 = 1 (F 1 + F )(x x 1 ) = 1 k(x x 1) (4) E = E 1 + W 1 = 1 kx (5) Now, lets apply ths method to Eshelby s ncluson problem. Let us consder the four steps n Eshelby s constructon of a sold contanng an ncluson. Recall that after step 1, the ncluson s outsde the matrx. The ncluson has undergone a deformaton due to ts egenstran. No forces are appled to ether the ncluson or the matrx. Obvously, the total elastc energy at ths state s E 1 = 0. In step, we apply a set of tracton forces on the ncluson surface. At the end of step, the tracton forces are T j = n j and the dsplacements on the surface are u j = e kj x k. Therefore, the work done n step s W 1 = 1 T j (x)u j (x) ds(x) = 1 n e kjx k ds(x) = 1 e kj x k n ds(x) (6) and usng Gauss s theorem W 1 = 1 e kj x k, d (x) = 1 e kj δ k d (x) = 1 e In step 3, the ncluson s put nsde the matrx wth the tracton force unchanged. No work s done n ths step,.e. W 3 = 0. In step 4, the tracton force T j s gradually reduced to zero. Both the ncluson and the matrx dsplace over a dstance of u c j. Snce the tracton force s T j at the begnnng of step 4 and 0 at the end of step 4, the average tracton force (7) 6
7 s, agan, T j /. The work done to the entre system (ncluson + matrx) s W 34 = 1 T j (x)u c j(x) ds(x) = 1 n u c j(x) ds(x) = 1 u c j,(x) ds(x) = 1 e c ds(x) = 1 e c (8) The total elastc energy at the end of step 4 s, E = E 1 + W 1 + W 3 + W 34 = e c e = 1 (c )e = 1 I e (9) whch s exactly the same as Eq. (19). The same approach can also be appled to obtan the elastc energy stored nsde the ncluson (E I ) or nsde the matrx (E M ). In step 4, the matrx also exert force on the ncluson, whch does work as the nterface moves. Ths leads to a transfer of elastc energy from the ncluson to the matrx. 3 Incluson energy n a fnte sold Let us now consder a problem wth an ncluson n a fnte sold. Agan, the stress-stran felds n ths case can be solved by superpostons. Suppose the fnte sold assumes the stress-stran felds of an nfnte sold contanng an ncluson, as solved prevously. Then we must apply a set of tracton forces T j to the outer surface S ext of the sold to mantan equlbrum. To obtan the soluton of a fnte sold wth zero tracton on ts outer surface, we need to remove T j on S ext. Ths s equvalent to apply a cancellng tracton force F j = T j on the outer surface S ext of the fnte sold. Let e m, m and u m be the stran, stress and dsplacement felds n response to the surface tracton F j. They are called mage felds. In ths case, the elastc felds nsde the matrx and the ncluson are, matrx e M M u M = e c + e m = c + m = u c + u m ncluson e I = e c e + e m I = c + m u I = u c e x j + u m 7
8 Obvously, the mage felds satsfy the condton, e m (x) = 1 (um,j(x) + u m j,(x)) (30) m (x) = C kl e m kl (x) (31) Notce that the mage felds are generally not unform wthn the ncluson. The free tracton boundary condton on the outer surface S ext can be expressed as, M n ext = 0 (on S ext ) (3) Smlar to Eq. (7), the total elastc energy n the sold can be expressed n terms of surface ntegrals, E = 1 ( u I I j M u M j )n out ds + M u M j n ext ds (33) S ext Because of Eq. (3), the second ntegral does not contrbute. Usng the tracton contnuty argument (n I out = M n out ) as before, we get E = 1 (u I I j u M j )n out ds (34) Agan, usng u I j u M j = e k x k, we get E = 1 e I jkx k n out ds (35) Ths s the same as Eq. (17) except that the stress feld nsde the ncluson now contans the mage component. Defne I, c (36) as the stress feld nsde the ncluson n an nfnte medum. Then I (x) = I, Smlarly, defne E 1 + m (x) (37) I, e jkx k n out ds = 1 I, e (38) as the ncluson energy n an nfnte sold. Then the ncluson energy n a fnte sold s, E = E 1 m e jkx k n out ds (39) 8
9 Convertng the second ntegral nto volume ntegral, we have E = E 1 ( m e jkx k ), d = E 1 m e d = E 1 e d where m m = E 1 m e = E + E m 1 m (x) d (x) (40) s the averaged mage stress nsde the ncluson. E m 1 m e s the mage contrbuton to the total ncluson energy. ncluson s, I I, + m (41) The average stress nsde the (4) Thus the total ncluson energy s stll related to the averaged stress nsde the ncluson as E = 1 I e (43) The results of the total ncluson energy for ellpsodal ncluson under varous boundary condtons are summarzed below. total elastc energy nfnte sold E = 1 I e fnte sold wth zero tracton E = 1 I e 4 Colonett s theorem I = I, + m We now wsh to study the energy of a sold contanng an ncluson and also subjectng to appled forces at ts outer surface. Before we do that, let us frst prove Colonett s theorem, whch s very helpful to study such problems. Colonett s theorem [1] states that There s no cross term n the total elastc energy of a sold, between the nternal stress feld and the appled stress feld. 9
10 However, there s an nteracton energy term between the nternal and appled felds when the energy of the appled loads s ncluded. Colonett s theorem can greatly smplfy the energy expressons when we apply stress to a fnte sold contanng an ncluson. To apprecate Colonett s theorem, we need to be specfc about the meanng of nternal and appled stress felds. Let us start wth a stress-free homogenous sold wth outer surface S ext. Defne nternal stress felds as the response to a heterogeneous feld of egenstran nsde the sold wth zero tracton on S ext. Defne appled stress felds as the response to a set of tractons on S ext when there s no egenstran nsde the sold. Let us consder two states of stress. State 1 s purely nternal, and state s appled. The total elastc energy nsde the sold for these two states are, E (1) = 1 (1) e(1) d E () = 1 () e() d Now consder a state 1+ whch s the superposton of state 1 and. Its total elastc energy should be, E (1+) = 1 ( (1) + () )(e(1) + e () ) d = 1 ( (1) e(1) + (1) e() + () e(1) + () e() ) d where E (1 ) 1 = E (1) + E () + E (1 ) ( (1) e() + () e(1) ) d (44) s the nteracton term between state 1 and state. Colonett s theorem states that E (1 ) must be zero, whch we wll prove below. Frst, we note that so that () e(1) = C kl e () kl e(1) (1) e() = C kl e (1) kl e() () e(1) = (1) e() E (1 ) = (1) e() d = (45) () e(1) d (46) Snce there s no body force, (1), = 0. Therefore, E (1 ) = ( (1) u() j ), d (47) 10
11 Now, we wsh to apply Gauss s theorem to convert the volume ntegral nto a surface ntegral. However, to use Gauss s theorem, the ntegrand must be contnuous nsde the entre volume. However, ths s not necessarly the case f the egenstran feld e (x) s not suffcently smooth. For example, n Eshelby s transformed ncluson problem, e (x) s not contnuous at the ncluson surface. As a result, the nternal stress feld (1) (x) s not contnuous at the ncluson surface ether. However, for clarty, let us assume for the moment that the egenstran feld e (x) and the nternal stress feld (1) (x) are suffcently smooth for the Gauss s theorem to apply. Ths corresponds to the case of thermal stran nduced by a smooth varaton of temperature nsde the sold. In ths case, E (1 ) = ( (1) ), d = u() j S ext n ext (1) u() j ds (48) Snce (1) s a purely nternal stress state, Hence (1) next = 0 (on S ext ) (49) E (1 ) = 0 (50) whch s Colonett s theorem. Let us now consder the case where the egenstran feld e (x) and the nternal stress feld (1) (x) are pecewse smooth nsde varous ncluson volumes K as well as n the matrx K K. Let n out,k be the outward normal vector of ncluson volume K. We can apply Gauss s theorem n each ncluson and the matrx separately, whch gves, E (1 ) = (1) u() j n ext ds + ( (1),K (1) )u() j n out,k ds (51) S ext K S K where (1),K s the stress nsde the volume k and (1) s the stress n matrx. The tracton force across the ncluson nterface must be contnuous,.e., ( K (1) )nout,k = 0 for any K, (5) where the summaton s not mpled over K. Therefore, agan we have, E (1 ) = n ext ds = 0 (53) (1) u() j S ext whch s Colonett s theorem. Colonett s theorem only deals wth the elastc stran energy that s stored nsde the sold (nternal energy). When the system s under appled load, ts evoluton proceeds 11
12 towards mnmzng ts enthalpy, whch s the nternal energy subtractng the work done by the external force. For example, the enthalpy of a system under external pressure p s H = E + p. The enthalpy for the sold under study s, H = E (1+) W W s the work done by the loadng mechansm, W = j + u () j S ext () (u(1) whch can be wrtten as where W (1 ) W = W (1 ) W (1 ) = W () = (54) )n ext ds (55) + W () (56) () u(1) j n ext S () u() j n ext S ds ds can be regarded as the cross term between the two stress states n the total enthalpy. 5 Fnte sold wth appled tractons We now apply Colonett s theorem to our problem of an ncluson n a fnte sold under a set of appled tractons. We wll use superscrpt A to denote the felds n response to the appled tractons when the egenstran vanshes (no ncluson). Let superscrpt F denote the felds of an ncluson n a fnte sold under zero external tractons (as n secton 3). From Colonett s theorem, where E = E A + E F E A = 1 E F = 1 (I, The enthalpy of the system s H = E W (57) A e A d (58) + m )e (59) where the A and F felds do have nteracton terms n the work term W,.e., W = W A + W A F (61) W A = u A A j n ext ds = e A A d = E A (6) S ext W A F = u A F j n ext ds (63) S ext 1 (60)
13 We would lke to express W A F result s, where n terms of an ntegral over the ncluson volume. The W A F = e A (64) A 1 A d (65) To show that ths s the case, frst note that n F ext S ext. Thus W A F = (u A F,M j S ext F,M u A j ) n ext ds = 0 on the surface where the superscrpt M denotes the felds n the matrx. Consder a volume ntegral of the same ntegrand over the matrx volume M =, 0 = (e A F,M F,M e A ) d M = (u A F,M j F,M u A j )n ext ds (u A F,M j F,M u A j )n out ds S ext whch means (u A F,M j S ext F,M u A j )n ext ds = (u A F,M j F,M u A j )n out ds (66) Hence W A F = (u A F,M j F,M u A j )n out ds (67) Notce that the ntegral s on the matrx sde of the ncluson nterface. We can smlarly wrte out the volume ntegral nsde the ncluson 0 = (e A F,I F,I e A ) d 0 = (u A F,I j F,I u A j )n out ds whch means that u A F,I j n out ds = F,I u A j n out ds (68) 13
14 Substtutng ths nto Eq. (67) and notng the tracton contnuty condton F,I n out = F,M n out, we have, W A F = (u A F,M j u F,I j )n out ds S 0 = e A jkx k n out ds (69) Usng Gauss s theorem, we get, W A F = e A d = e A d = e A (70) The major results of ths lecture note are summarzed below. total elastc energy nfnte sold E = 1 I e fnte sold wth zero tracton E = 1 I e E = E A fnte sold + E F E A = 1 wth tracton E F = 1 I e References A e A d I = I, I = I, + m + m total enthalpy H = E W = E A + E F W A A F W = E F E A Ae [1] G. Colonett, Per una teora generale delle coazon elastche, Att R. Acad. Sc. Torno: Cl. Sc. fs. mat. natur., 56, , (191). 14
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 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 informationPES 1120 Spring 2014, Spendier Lecture 6/Page 1
PES 110 Sprng 014, Spender Lecture 6/Page 1 Lecture today: Chapter 1) Electrc feld due to charge dstrbutons -> charged rod -> charged rng We ntroduced the electrc feld, E. I defned t as an nvsble aura
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 informationOne Dimensional Axial Deformations
One Dmensonal al Deformatons In ths secton, a specfc smple geometr s consdered, that of a long and thn straght component loaded n such a wa that t deforms n the aal drecton onl. The -as s taken as the
More information1 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 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 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 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 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 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 information1 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations
Physcs 171/271 -Davd Klenfeld - Fall 2005 (revsed Wnter 2011) 1 Dervaton of Rate Equatons from Sngle-Cell Conductance (Hodgkn-Huxley-lke) Equatons We consder a network of many neurons, each of whch obeys
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 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 informationCelestial Mechanics. Basic Orbits. Why circles? Tycho Brahe. PHY celestial-mechanics - J. Hedberg
PHY 454 - celestal-mechancs - J. Hedberg - 207 Celestal Mechancs. Basc Orbts. Why crcles? 2. Tycho Brahe 3. Kepler 4. 3 laws of orbtng bodes 2. Newtonan Mechancs 3. Newton's Laws. Law of Gravtaton 2. The
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 information1 (1 + ( )) = 1 8 ( ) = (c) Carrying out the Taylor expansion, in this case, the series truncates at second order:
68A Solutons to Exercses March 05 (a) Usng a Taylor expanson, and notng that n 0 for all n >, ( + ) ( + ( ) + ) We can t nvert / because there s no Taylor expanson around 0 Lets try to calculate the nverse
More informationEPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski
EPR Paradox and the Physcal Meanng of an Experment n Quantum Mechancs Vesseln C Nonnsk vesselnnonnsk@verzonnet Abstract It s shown that there s one purely determnstc outcome when measurement s made on
More informationPhysics 5153 Classical Mechanics. Principle of Virtual Work-1
P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal
More informationOpen Systems: Chemical Potential and Partial Molar Quantities Chemical Potential
Open Systems: Chemcal Potental and Partal Molar Quanttes Chemcal Potental For closed systems, we have derved the followng relatonshps: du = TdS pdv dh = TdS + Vdp da = SdT pdv dg = VdP SdT For open systems,
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 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 informationFrom 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 informationMath1110 (Spring 2009) Prelim 3 - Solutions
Math 1110 (Sprng 2009) Solutons to Prelm 3 (04/21/2009) 1 Queston 1. (16 ponts) Short answer. Math1110 (Sprng 2009) Prelm 3 - Solutons x a 1 (a) (4 ponts) Please evaluate lm, where a and b are postve numbers.
More information9 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations
Physcs 171/271 - Chapter 9R -Davd Klenfeld - Fall 2005 9 Dervaton of Rate Equatons from Sngle-Cell Conductance (Hodgkn-Huxley-lke) Equatons We consder a network of many neurons, each of whch obeys a set
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 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 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 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 information8.6 The Complex Number System
8.6 The Complex Number System Earler n the chapter, we mentoned that we cannot have a negatve under a square root, snce the square of any postve or negatve number s always postve. In ths secton we want
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 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 informationLectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix
Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could
More 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 informationMore 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 informationModule 1 : The equation of continuity. Lecture 1: Equation of Continuity
1 Module 1 : The equaton of contnuty Lecture 1: Equaton of Contnuty 2 Advanced Heat and Mass Transfer: Modules 1. THE EQUATION OF CONTINUITY : Lectures 1-6 () () () (v) (v) Overall Mass Balance Momentum
More informationRate of Absorption and Stimulated Emission
MIT Department of Chemstry 5.74, Sprng 005: Introductory Quantum Mechancs II Instructor: Professor Andre Tokmakoff p. 81 Rate of Absorpton and Stmulated Emsson The rate of absorpton nduced by the feld
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 information8.022 (E&M) Lecture 4
Topcs: 8.0 (E&M) Lecture 4 More applcatons of vector calculus to electrostatcs: Laplacan: Posson and Laplace equaton url: concept and applcatons to electrostatcs Introducton to conductors Last tme Electrc
More informationExercise Solutions to Real Analysis
xercse Solutons to Real Analyss Note: References refer to H. L. Royden, Real Analyss xersze 1. Gven any set A any ɛ > 0, there s an open set O such that A O m O m A + ɛ. Soluton 1. If m A =, then there
More informationAPPENDIX A Some Linear Algebra
APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,
More informationStanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011
Stanford Unversty CS359G: Graph Parttonng and Expanders Handout 4 Luca Trevsan January 3, 0 Lecture 4 In whch we prove the dffcult drecton of Cheeger s nequalty. As n the past lectures, consder an undrected
More informationIntroduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:
CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and
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 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 12. Ordinary Differential Equation Boundary Value (BV) Problems
Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(
More informationAERODYNAMICS I LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY
LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY The Bot-Savart Law The velocty nduced by the sngular vortex lne wth the crculaton can be determned by means of the Bot- Savart formula
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 informationON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION
Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION
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 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 informationCIVL 8/7117 Chapter 10 - Isoparametric Formulation 42/56
CIVL 8/77 Chapter 0 - Isoparametrc Formulaton 4/56 Newton-Cotes Example Usng the Newton-Cotes method wth = ntervals (n = 3 samplng ponts), evaluate the ntegrals: x x cos dx 3 x x dx 3 x x dx 4.3333333
More information1 NOTATIONS 2 NUMERICAL SIMULATIONS 3 TWO DIMENSIONAL CIRCULAR INCLUSION. Supplementary Notes
Supplementary Notes 1 NOTATIONS In the followng notes, we wll reserve the Greek letters α,β,γ,... for tensor subscrpts, whle the Latn letters,j,k,... wll be used for partcles and Eshelby nclusons. NUMERICAL
More information3.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 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 informationCalculus of Variations Basics
Chapter 1 Calculus of Varatons Bascs 1.1 Varaton of a General Functonal In ths chapter, we derve the general formula for the varaton of a functonal of the form J [y 1,y 2,,y n ] F x,y 1,y 2,,y n,y 1,y
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 informationΔ x. u(x,t) Fig. Schematic view of elastic bar undergoing axial motions
ME67 - Handout 4 Vbratons of Contnuous Systems Axal vbratons of elastc bars The fgure shows a unform elastc bar of length and cross secton A. The bar materal propertes are ts densty ρ and elastc modulus
More information8.323 Relativistic Quantum Field Theory I
MI OpenCourseWare http://ocw.mt.edu 8.323 Relatvstc Quantum Feld heory I Sprng 2008 For nformaton about ctng these materals or our erms of Use, vst: http://ocw.mt.edu/terms. MASSACHUSES INSIUE OF ECHNOLOGY
More informationAffine 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 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 informationMAE140 - Linear Circuits - Winter 16 Midterm, February 5
Instructons ME140 - Lnear Crcuts - Wnter 16 Mdterm, February 5 () Ths exam s open book. You may use whatever wrtten materals you choose, ncludng your class notes and textbook. You may use a hand calculator
More informationGeorgia Tech PHYS 6124 Mathematical Methods of Physics I
Georga Tech PHYS 624 Mathematcal Methods of Physcs I Instructor: Predrag Cvtanovć Fall semester 202 Homework Set #7 due October 30 202 == show all your work for maxmum credt == put labels ttle legends
More information763622S 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 information2.3 Nilpotent endomorphisms
s a block dagonal matrx, wth A Mat dm U (C) In fact, we can assume that B = B 1 B k, wth B an ordered bass of U, and that A = [f U ] B, where f U : U U s the restrcton of f to U 40 23 Nlpotent endomorphsms
More informationOne-sided finite-difference approximations suitable for use with Richardson extrapolation
Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,
More informationIntegrals and Invariants of Euler-Lagrange Equations
Lecture 16 Integrals and Invarants of Euler-Lagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,
More informationA boundary element method with analytical integration for deformation of inhomogeneous elastic materials
Journal of Physcs: Conference Seres PAPER OPEN ACCESS A boundary element method wth analytcal ntegraton for deformaton of nhomogeneous elastc materals To cte ths artcle: Moh. Ivan Azs et al 2018 J. Phys.:
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 informationSTATIC ANALYSIS OF TWO-LAYERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION
STATIC ANALYSIS OF TWO-LERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION Ákos József Lengyel István Ecsed Assstant Lecturer Emertus Professor Insttute of Appled Mechancs Unversty of Mskolc Mskolc-Egyetemváros
More informationG4023 Mid-Term Exam #1 Solutions
Exam1Solutons.nb 1 G03 Md-Term Exam #1 Solutons 1-Oct-0, 1:10 p.m to :5 p.m n 1 Pupn Ths exam s open-book, open-notes. You may also use prnt-outs of the homework solutons and a calculator. 1 (30 ponts,
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 informationn α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0
MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector
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 informationLecture 6/7 (February 10/12, 2014) DIRAC EQUATION. The non-relativistic Schrödinger equation was obtained by noting that the Hamiltonian 2
P470 Lecture 6/7 (February 10/1, 014) DIRAC EQUATION The non-relatvstc Schrödnger equaton was obtaned by notng that the Hamltonan H = P (1) m can be transformed nto an operator form wth the substtutons
More informationSolution Thermodynamics
Soluton hermodynamcs usng Wagner Notaton by Stanley. Howard Department of aterals and etallurgcal Engneerng South Dakota School of nes and echnology Rapd Cty, SD 57701 January 7, 001 Soluton hermodynamcs
More informationMATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018
MATH 5630: Dscrete Tme-Space Model Hung Phan, UMass Lowell March, 08 Newton s Law of Coolng Consder the coolng of a well strred coffee so that the temperature does not depend on space Newton s law of collng
More informationLinear Momentum. Center of Mass.
Lecture 6 Chapter 9 Physcs I 03.3.04 Lnear omentum. Center of ass. Course webste: http://faculty.uml.edu/ndry_danylov/teachng/physcsi Lecture Capture: http://echo360.uml.edu/danylov03/physcssprng.html
More informationImplicit 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 informationElectrical double layer: revisit based on boundary conditions
Electrcal double layer: revst based on boundary condtons Jong U. Km Department of Electrcal and Computer Engneerng, Texas A&M Unversty College Staton, TX 77843-318, USA Abstract The electrcal double layer
More informationA Quantum Gauss-Bonnet Theorem
A Quantum Gauss-Bonnet Theorem Tyler Fresen November 13, 2014 Curvature n the plane Let Γ be a smooth curve wth orentaton n R 2, parametrzed by arc length. The curvature k of Γ s ± Γ, where the sgn s postve
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 informationSpring Force and Power
Lecture 13 Chapter 9 Sprng Force and Power Yeah, energy s better than orces. What s net? Course webste: http://aculty.uml.edu/andry_danylov/teachng/physcsi IN THIS CHAPTER, you wll learn how to solve problems
More information4DVAR, according to the name, is a four-dimensional variational method.
4D-Varatonal Data Assmlaton (4D-Var) 4DVAR, accordng to the name, s a four-dmensonal varatonal method. 4D-Var s actually a drect generalzaton of 3D-Var to handle observatons that are dstrbuted n tme. The
More information8 Derivation of Network Rate Equations from Single- Cell Conductance Equations
Physcs 178/278 - Davd Klenfeld - Wnter 2015 8 Dervaton of Network Rate Equatons from Sngle- Cell Conductance Equatons We consder a network of many neurons, each of whch obeys a set of conductancebased,
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 information1. Review of Mechanics Newton s Laws
. Revew of Mechancs.. Newton s Laws Moton of partcles. Let the poston of the partcle be gven by r. We can always express ths n Cartesan coordnates: r = xˆx + yŷ + zẑ, () where we wll always use ˆ (crcumflex)
More informationPHYS 705: Classical Mechanics. Canonical Transformation II
1 PHYS 705: Classcal Mechancs Canoncal Transformaton II Example: Harmonc Oscllator f ( x) x m 0 x U( x) x mx x LT U m Defne or L p p mx x x m mx x H px L px p m p x m m H p 1 x m p m 1 m H x p m x m m
More informationFormal solvers of the RT equation
Formal solvers of the RT equaton Formal RT solvers Runge- Kutta (reference solver) Pskunov N.: 979, Master Thess Long characterstcs (Feautrer scheme) Cannon C.J.: 970, ApJ 6, 55 Short characterstcs (Hermtan
More informationVisco-Rubber Elastic Model for Pressure Sensitive Adhesive
Vsco-Rubber Elastc Model for Pressure Senstve Adhesve Kazuhsa Maeda, Shgenobu Okazawa, Koj Nshgch and Takash Iwamoto Abstract A materal model to descrbe large deformaton of pressure senstve adhesve (PSA
More informationForce = F Piston area = A
CHAPTER III Ths chapter s an mportant transton between the propertes o pure substances and the most mportant chapter whch s: the rst law o thermodynamcs In ths chapter, we wll ntroduce the notons o heat,
More informationprinceton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg
prnceton unv. F 17 cos 521: Advanced Algorthm Desgn Lecture 7: LP Dualty Lecturer: Matt Wenberg Scrbe: LP Dualty s an extremely useful tool for analyzng structural propertes of lnear programs. Whle there
More informationSolutions to Problem Set 6
Solutons to Problem Set 6 Problem 6. (Resdue theory) a) Problem 4.7.7 Boas. n ths problem we wll solve ths ntegral: x sn x x + 4x + 5 dx: To solve ths usng the resdue theorem, we study ths complex ntegral:
More informationFrequency dependence of the permittivity
Frequency dependence of the permttvty February 7, 016 In materals, the delectrc constant and permeablty are actually frequency dependent. Ths does not affect our results for sngle frequency modes, but
More information8 Derivation of Network Rate Equations from Single- Cell Conductance Equations
Physcs 178/278 - Davd Klenfeld - Wnter 2019 8 Dervaton of Network Rate Equatons from Sngle- Cell Conductance Equatons Our goal to derve the form of the abstract quanttes n rate equatons, such as synaptc
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 informationHow Differential Equations Arise. Newton s Second Law of Motion
page 1 CHAPTER 1 Frst-Order Dfferental Equatons Among all of the mathematcal dscplnes the theory of dfferental equatons s the most mportant. It furnshes the explanaton of all those elementary manfestatons
More informationChapter 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 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 informationPHYS 1101 Practice problem set 12, Chapter 32: 21, 22, 24, 57, 61, 83 Chapter 33: 7, 12, 32, 38, 44, 49, 76
PHYS 1101 Practce problem set 1, Chapter 3: 1,, 4, 57, 61, 83 Chapter 33: 7, 1, 3, 38, 44, 49, 76 3.1. Vsualze: Please reer to Fgure Ex3.1. Solve: Because B s n the same drecton as the ntegraton path s
More information