Plates on elastic foundation

Size: px
Start display at page:

Download "Plates on elastic foundation"

Transcription

1 Pltes on elstic foundtion Circulr elstic plte, xil-symmetric lod, Winkler soil (fter Timoshenko & Woinowsky-Krieger (1959) - Chpter 8) Prepred by Enzo Mrtinelli Drft version ( April 016)

2 Introduction Winkler model: modulus of subgrde 0 σ s = k0w N/mm kg f /cm

3 Introduction Circulr elstic (thin) pltes The Kirchhoff Love theory of pltes is two-dimensionl mthemticl model tht is used to determine the stresses nd deformtions in thin pltes subjected to forces nd moments. The following kinemtic ssumptions tht re mde in this theory: stright lines norml to the mid-surfce remin stright fter deformtion stright lines norml to the mid-surfce remin norml to the mid-surfce fter deformtion the thickness of the plte does not chnge during deformtion. (fter Wikipedi) Germin-Lgrnge eqution (crtesin coordintes) w w w q x x y y D Prepred by Enzo Mrtinelli Drft version ( April 016)

4 Introduction Circulr elstic (thin) pltes Germin-Lgrnge eqution (crtesin coordintes) q ww D x y Conversion to polr (cylindricl) coordintes y r x r cos y r sin x w w 1 w cos sin x r r w w 1 w sin cos y r r 1 1 r r r r Prepred by Enzo Mrtinelli Drft version ( April 016)

5 Introduction Circulr elstic (thin) pltes under xil-symmetric lod Germin-Lgrnge eqution (polr coordintes) q w w D 1 1 r r r r 1 1 w 1 w 1 w q w r r r r r r r r D d 1 d d w 1 dw q dr r dr dr r dr D Prepred by Enzo Mrtinelli Drft version ( April 016)

6 Introduction Circulr elstic (thin) pltes on Winkler soil under xil-symmetric lod d 1 d d w 1 dw q k0w dr r dr dr r dr D q k w d w d w 1 d w 1 dw 0 dr r dr r dr r dr D Since n uniformly distributed lod q pplied on n elstic pltes (or elstic bems) on Winkler soils with uniform k 0 does not determine ny stress within the plte, the sole effect of point lod P pplied in the centre of the plte is considered herefter. dy q dx k 0 w q=0 P R Prepred by Enzo Mrtinelli Drft version ( April 016)

7 Conversion to dimensionless form - Chrcteristic length: k0 1 D l [F ][L] 1 [F ][L] [L] dr r dr dr r dr l d 1 d d w 1 dw w 0 - Dimensionless quntities: w Dimensionless displcement z l r Dimensionless rdius x 1 d 1 d d z 1 dz z l 0 l dx x dx dx x dx l Dimensionless derivtives: dw l dz dz d 1 d d z 1 dz z 0 dr l dx dx dx x dx dx x dx d w d dw d dz 1 d z dr dr dr l dx dx l dx i i d w d dw 1 d w dr dr dr l dx i i1 i d z d z 1 d z 1 dz dx x dx x dx x dx z 0 Prepred by Enzo Mrtinelli Drft version ( April 016)

8 Generl solution The proposed mthemticl trnsformtions led to the following homogeneous th -order liner differentil eqution with vrible coefficients: d z d z 1 d z 1 dz dx x dx x dx x dx z 0 Therefore, in principle the generl solution of this differentil eqution might be written s follows: z A X x A X x A X x A X x 1 1 where X i re four independent solutions of the differentil eqution under considertion nd A i re four integrtion constnts depending on the ctul boundry conditions. Note: d 1 d 1 dx x dx dx x dx d Xi dxi i X 0 i 1.. Prepred by Enzo Mrtinelli Drft version ( April 016)

9 Approximtion by power series Since X i re unknown (s the differentil eqution under considertion includes vrible coefficients), n pproximtion bsed on power series my be firstly ssumed. Hence, the generl expression of Xi cn be tken s follows: NT n 1 n NT Xi nx 0 1x x... nx... N x T n0 where N T is the highest monomil in the series. Since X i should be solution of the generl differentil eqution, the nd order lplcin of ech term n x n should find similr monomil such tht: X i X i 0 x d 1 d d x 1 d x dx x dx dx x dx n n n n n n n n x x 0 n n n n n n d nx 1 dnx n n 1 x n x n x dx x dx dx n n n n n n d 1 d n n n n n n x n n n n x n n n x n n n x x dx n n n n Recursive definition of n

10 Approximtion by power series: first independed solution X 1 n n n n z A X x A X x A X x A X x 1 1 Bsed on the recursive definition of the n coefficients, first independent solution X 1 cn be built up by ssuming 0 s the first nonzero coefficient: n 0 : n : n 8 : n 1 : X1 x 1 x x x Prepred by Enzo Mrtinelli Drft version ( April 016)

11 Approximtion by power series: second independed solution X n n n n Before proceeding with constructing the second independent solution X of z, it should be noted tht the ssumption of 1 =1 is not dmissible. In fct, if one ssumes 1 =1 (nd, t the sme time, 0 = = =0) the resulting solution X of z would hve first term x: X x... which implies nonzero first derivtive of X (nd, hence, of the generl solution z) which cnnot be ccepted becuse of the fct tht rottions (j=dx /dx=1+ ) should be zero t x=0. Therefore, the second solution X should be built by ssuming =1 nd, t the sme time, 0 = 1 = =0: X x... this leds to n expression of X whose first nonzero monomil is of nd order nd, hence, the second derivtive of X is non zero (c=d X /dx =1+ ) ): this is cceptble becuse the curvture is nonzero t the plte centre (x=0). Prepred by Enzo Mrtinelli Drft version ( April 016)

12 Approximtion by power series: second independed solution X n n n n Bsed on the recursive definition of the n coefficients, first independent solution X cn be built up by ssuming s the first nonzero coefficient: n : n 6 : n 10 : n 1 : X x x x x x Prepred by Enzo Mrtinelli Drft version ( April 016)

13 Approximtion by power series: third independed solution X n n n n The ssumption =1 (nd, t the sme time, 0 = 1 = =0) is not even dmissible, s it would led to the function hving the following expression: X x... This is not dmissible becuse it would led to finite vlue of sher for x=0: d d X 1 dx V... dx dx x dx wheres sher should be divergent (it should tend to infinity) for x0, s finite force P should be in equilibrium with sher stresses on n infinitesiml circumference of rdius e. P x=e V Prepred by Enzo Mrtinelli Drft version ( April 016)

14 Approximtion by power series: third independed solution X Since, on the one hnd, sher V for x0, nd on the other one hnd, both X 1 nd X would led to finite vlue of sher for x=0, the solutions X nd X should be cpble to reproduce the condition V for x0. Therefore, different solution, including the term log x (which diverges for x0 long with ll its derivtives) d n new unknown function F, is preliminrily defined: X X1 log x F The function X should be, in turn, solution of the generl differentil eqution: X X 0 i First of ll, the nd lplcin of X is determined i d 1 d 1 dx x dx dx x dx d Xi dxi i X 0 d X1 1 X log x X F x dx Prepred by Enzo Mrtinelli Drft version ( April 016)

15 Approximtion by power series: third independed solution X Then: The following condition hs to be met: X X 0 d X1 1 1 X X log x X F X log x F 0 x dx nd, fter collecting the term log x, d X1 x dx the following reltionship is derived for F : log x X X F F d X1 F F x dx F F x x... Prepred by Enzo Mrtinelli Drft version ( April 016)

16 Approximtion by power series: third independed solution X F F x x n (n 1) (n ) F F ( 1) x The function F my lso be pproximted by the following power series: F b x b x b x whose nd lplcin cn be expressed s follows: n T n 1 n n (n ).. n b x n n b x n n n n n (n 1) (n ) n n b x b x ( 1) x n (n ).. n 1 n n n n n n n b 1 1 n (n 1) (n ) n n b ( 1) n n n (n ).. Prepred by Enzo Mrtinelli Drft version ( April 016)

17 Approximtion by power series: third independed solution X n b 1 1 n (n 1) (n ) n n b ( 1) n n n (n ).. Prepred by Enzo Mrtinelli Drft version ( April 016)

18 Approximtion by power series: fourth independed solution X A similr procedure leds to defining the fourth independent solution: X X log x F nd, similr considertions, cn be done in defining F for X is solution of the generl differentil eqution. F c x c x c x Prepred by Enzo Mrtinelli Drft version ( April 016)

19 Approximtion by power series: fourth independed solution X Prepred by Enzo Mrtinelli Drft version ( April 016)

20 Boundry conditions The generl solution, whose expression ws defined by th pproximte solutions, is reported below: z A X x A X x A X x A X x 1 1 Four boundry conditions should be written for determining the four constnts A 1 A. The first two ones re imposed for r=r (nmely, for x=r/l): Dimensionl expression Dimensionless expression P R r R x R l dr d w dw r dr rr 0 dx d z dz x dx R x l 0 d d w 1 dw dr dr r dr rr 0 d d z 1 dz dx dx x dx R x l 0 Prepred by Enzo Mrtinelli Drft version ( April 016)

21 Boundry conditions The generl solution, whose expression ws defined by th pproximte solutions, is reported below: r 0 dw dr rr z A X x A X x A X x A X x Other two boundry conditions re imposed for r=0 (nmely, for x=0): Dimensionl expression Dimensionless expression x R l dz dx x0 0 P R P el V d d w 1 dw lim e l D P 0 dr dr r dr re l e 0 d d z 1 dz lim e k0l P 0 dx dx x dx xe e 0 Prepred by Enzo Mrtinelli Drft version ( April 016)

22 Boundry conditions: expressions of constnts First of ll, let s consider the third boundry condition: Since X 1 nd X re polynomils with even order power terms (i.e. x 0, x, x 8 nd x, x 6, x 10, respectively), their first derivtive is zero for x=0. Moreover, the first derivtive of X is lso zero for x=0 (or, better, in the limit of x 0): dx d x lim X log x F lim x log x dx dx x e0 e0 x0 xe xe Therefore, the first derivtive of z in x=0 only tkes nonzero contribution from X : dz dx d 1 A A lim X log x F A lim... 0 dx dx dx x 1 e0 e0 x0 x0 xe xe which is only stisfied for A =0. dz dx x0 0 Prepred by Enzo Mrtinelli Drft version ( April 016)

23 Boundry conditions: expressions of constnts Then, the fourth condition cn be considered: nd, once gin, no contribution comes from X 1 nd X, but the only nonzero terms is given by the third derivtives of X : d d z 1 dz d d X 1 dx d d 1 d A A X log x F dx dx x dx dx dx x dx dx dx x dx xe d d z 1 dz lim e k0l P 0 dx dx x dx xe e 0 d d 1 d d d 1 d A X log x F A x log x... dx dx x dx dx dx x dx Therefore: A xe xe A x x e e x xe xe lim e k0l A P 0 e e 0 P 8 k0l A P 0 A 8 k l 0 Prepred by Enzo Mrtinelli Drft version ( April 016)

24 Boundry conditions: expressions of constnts Finlly, the first two equtions cn be considered for determining the two constnts A 1 nd A : dx d z dz x dx R x l 0 d d z 1 dz dx dx x dx R x l 0 z x A X x A X x A X x 1 1 w r r l z l d w Mr r D dr r dr dw M r D r dr dr 1 dw d w d d w 1 dw Vr r D dr dr r dr Prepred by Enzo Mrtinelli Drft version ( April 016)

25 Comprisons: proposed procedure vs FEM solution N T = N T = Prepred by Enzo Mrtinelli Drft version ( April 016)

Kirchhoff and Mindlin Plates

Kirchhoff and Mindlin Plates Kirchhoff nd Mindlin Pltes A plte significntly longer in two directions compred with the third, nd it crries lod perpendiculr to tht plne. The theory for pltes cn be regrded s n extension of bem theory,

More information

1 Bending of a beam with a rectangular section

1 Bending of a beam with a rectangular section 1 Bending of bem with rectngulr section x3 Episseur b M x 2 x x 1 2h M Figure 1 : Geometry of the bem nd pplied lod The bem in figure 1 hs rectngur section (thickness 2h, width b. The pplied lod is pure

More information

Plate Theory. Section 13: PLATE BENDING ELEMENTS

Plate Theory. Section 13: PLATE BENDING ELEMENTS Section : PLATE BENDING ELEENTS Wshkeic College of Engineering Plte Theor A plte is structurl element hose mid surfce lies in flt plne. The dimension in the direction norml to the plne is referred to s

More information

Plate Theory. Section 11: PLATE BENDING ELEMENTS

Plate Theory. Section 11: PLATE BENDING ELEMENTS Section : PLATE BENDING ELEMENTS Plte Theor A plte is structurl element hose mid surfce lies in flt plne. The dimension in the direction norml to the plne is referred to s the thickness of the plte. A

More information

Partial Differential Equations

Partial Differential Equations Prtil Differentil Equtions Notes by Robert Piché, Tmpere University of Technology reen s Functions. reen s Function for One-Dimensionl Eqution The reen s function provides complete solution to boundry

More information

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1 3 9. SEQUENCES AND SERIES 63. Representtion of functions s power series Consider power series x 2 + x 4 x 6 + x 8 + = ( ) n x 2n It is geometric series with q = x 2 nd therefore it converges for ll q =

More information

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007 A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H Thoms Shores Deprtment of Mthemtics University of Nebrsk Spring 2007 Contents Rtes of Chnge nd Derivtives 1 Dierentils 4 Are nd Integrls 5 Multivrite Clculus

More information

Summary: Method of Separation of Variables

Summary: Method of Separation of Variables Physics 246 Electricity nd Mgnetism I, Fll 26, Lecture 22 1 Summry: Method of Seprtion of Vribles 1. Seprtion of Vribles in Crtesin Coordintes 2. Fourier Series Suggested Reding: Griffiths: Chpter 3, Section

More information

Math 124A October 04, 2011

Math 124A October 04, 2011 Mth 4A October 04, 0 Viktor Grigoryn 4 Vibrtions nd het flow In this lecture we will derive the wve nd het equtions from physicl principles. These re second order constnt coefficient liner PEs, which model

More information

1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.

1.2. Linear Variable Coefficient Equations. y + b ! = a y + b  Remark: The case b = 0 and a non-constant can be solved with the same idea as above. 1 12 Liner Vrible Coefficient Equtions Section Objective(s): Review: Constnt Coefficient Equtions Solving Vrible Coefficient Equtions The Integrting Fctor Method The Bernoulli Eqution 121 Review: Constnt

More information

Partial Derivatives. Limits. For a single variable function f (x), the limit lim

Partial Derivatives. Limits. For a single variable function f (x), the limit lim Limits Prtil Derivtives For single vrible function f (x), the limit lim x f (x) exists only if the right-hnd side limit equls to the left-hnd side limit, i.e., lim f (x) = lim f (x). x x + For two vribles

More information

Section 6.1 INTRO to LAPLACE TRANSFORMS

Section 6.1 INTRO to LAPLACE TRANSFORMS Section 6. INTRO to LAPLACE TRANSFORMS Key terms: Improper Integrl; diverge, converge A A f(t)dt lim f(t)dt Piecewise Continuous Function; jump discontinuity Function of Exponentil Order Lplce Trnsform

More information

Space Curves. Recall the parametric equations of a curve in xy-plane and compare them with parametric equations of a curve in space.

Space Curves. Recall the parametric equations of a curve in xy-plane and compare them with parametric equations of a curve in space. Clculus 3 Li Vs Spce Curves Recll the prmetric equtions of curve in xy-plne nd compre them with prmetric equtions of curve in spce. Prmetric curve in plne x = x(t) y = y(t) Prmetric curve in spce x = x(t)

More information

Jackson 2.26 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell

Jackson 2.26 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell Jckson 2.26 Homework Problem Solution Dr. Christopher S. Bird University of Msschusetts Lowell PROBLEM: The two-dimensionl region, ρ, φ β, is bounded by conducting surfces t φ =, ρ =, nd φ = β held t zero

More information

A ROTATING DISC IN CONSTANT PURE SHEAR BY S. KUMAR AND C. V. JOGA RAO

A ROTATING DISC IN CONSTANT PURE SHEAR BY S. KUMAR AND C. V. JOGA RAO A ROTATING DISC IN CONSTANT PURE SHEAR BY S. KUMAR AND C. V. JOGA RAO (Deprtment of Aeronuticl Engineering, Indin Institute of Science, Bnglore-3) Received April 25, 1954 SUMMARY The disc of constnt pure

More information

Math 426: Probability Final Exam Practice

Math 426: Probability Final Exam Practice Mth 46: Probbility Finl Exm Prctice. Computtionl problems 4. Let T k (n) denote the number of prtitions of the set {,..., n} into k nonempty subsets, where k n. Argue tht T k (n) kt k (n ) + T k (n ) by

More information

KINEMATICS OF RIGID BODIES

KINEMATICS OF RIGID BODIES KINEMTICS OF RIGID ODIES Introduction In rigid body kinemtics, e use the reltionships governing the displcement, velocity nd ccelertion, but must lso ccount for the rottionl motion of the body. Description

More information

Solution Manual. for. Fracture Mechanics. C.T. Sun and Z.-H. Jin

Solution Manual. for. Fracture Mechanics. C.T. Sun and Z.-H. Jin Solution Mnul for Frcture Mechnics by C.T. Sun nd Z.-H. Jin Chpter rob.: ) 4 No lod is crried by rt nd rt 4. There is no strin energy stored in them. Constnt Force Boundry Condition The totl strin energy

More information

Consequently, the temperature must be the same at each point in the cross section at x. Let:

Consequently, the temperature must be the same at each point in the cross section at x. Let: HW 2 Comments: L1-3. Derive the het eqution for n inhomogeneous rod where the therml coefficients used in the derivtion of the het eqution for homogeneous rod now become functions of position x in the

More information

10 Vector Integral Calculus

10 Vector Integral Calculus Vector Integrl lculus Vector integrl clculus extends integrls s known from clculus to integrls over curves ("line integrls"), surfces ("surfce integrls") nd solids ("volume integrls"). These integrls hve

More information

Pressure Wave Analysis of a Cylindrical Drum

Pressure Wave Analysis of a Cylindrical Drum Pressure Wve Anlysis of Cylindricl Drum Chris Clrk, Brin Anderson, Brin Thoms, nd Josh Symonds Deprtment of Mthemtics The University of Rochester, Rochester, NY 4627 (Dted: December, 24 In this pper, hypotheticl

More information

Note 16. Stokes theorem Differential Geometry, 2005

Note 16. Stokes theorem Differential Geometry, 2005 Note 16. Stokes theorem ifferentil Geometry, 2005 Stokes theorem is the centrl result in the theory of integrtion on mnifolds. It gives the reltion between exterior differentition (see Note 14) nd integrtion

More information

MAC-solutions of the nonexistent solutions of mathematical physics

MAC-solutions of the nonexistent solutions of mathematical physics Proceedings of the 4th WSEAS Interntionl Conference on Finite Differences - Finite Elements - Finite Volumes - Boundry Elements MAC-solutions of the nonexistent solutions of mthemticl physics IGO NEYGEBAUE

More information

MA Handout 2: Notation and Background Concepts from Analysis

MA Handout 2: Notation and Background Concepts from Analysis MA350059 Hndout 2: Nottion nd Bckground Concepts from Anlysis This hndout summrises some nottion we will use nd lso gives recp of some concepts from other units (MA20023: PDEs nd CM, MA20218: Anlysis 2A,

More information

New Expansion and Infinite Series

New Expansion and Infinite Series Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University

More information

Equations of Motion. Figure 1.1.1: a differential element under the action of surface and body forces

Equations of Motion. Figure 1.1.1: a differential element under the action of surface and body forces Equtions of Motion In Prt I, lnce of forces nd moments cting on n component ws enforced in order to ensure tht the component ws in equilirium. Here, llownce is mde for stresses which vr continuousl throughout

More information

THREE-DIMENSIONAL KINEMATICS OF RIGID BODIES

THREE-DIMENSIONAL KINEMATICS OF RIGID BODIES THREE-DIMENSIONAL KINEMATICS OF RIGID BODIES 1. TRANSLATION Figure shows rigid body trnslting in three-dimensionl spce. Any two points in the body, such s A nd B, will move long prllel stright lines if

More information

Explain shortly the meaning of the following eight words in relation to shells structures.

Explain shortly the meaning of the following eight words in relation to shells structures. Delft University of Technology Fculty of Civil Engineering nd Geosciences Structurl Mechnics Section Write your nme nd study number t the top right-hnd of your work. Exm CIE4143 Shell Anlysis Tuesdy 15

More information

Chapter 3. Vector Spaces

Chapter 3. Vector Spaces 3.4 Liner Trnsformtions 1 Chpter 3. Vector Spces 3.4 Liner Trnsformtions Note. We hve lredy studied liner trnsformtions from R n into R m. Now we look t liner trnsformtions from one generl vector spce

More information

Math 113 Exam 1-Review

Math 113 Exam 1-Review Mth 113 Exm 1-Review September 26, 2016 Exm 1 covers 6.1-7.3 in the textbook. It is dvisble to lso review the mteril from 5.3 nd 5.5 s this will be helpful in solving some of the problems. 6.1 Are Between

More information

Eigen Values and Eigen Vectors of a given matrix

Eigen Values and Eigen Vectors of a given matrix Engineering Mthemtics 0 SUBJECT NAME SUBJECT CODE MATERIAL NAME MATERIAL CODE : Engineering Mthemtics I : 80/MA : Prolem Mteril : JM08AM00 (Scn the ove QR code for the direct downlod of this mteril) Nme

More information

Conservation Law. Chapter Goal. 5.2 Theory

Conservation Law. Chapter Goal. 5.2 Theory Chpter 5 Conservtion Lw 5.1 Gol Our long term gol is to understnd how mny mthemticl models re derived. We study how certin quntity chnges with time in given region (sptil domin). We first derive the very

More information

1 1D heat and wave equations on a finite interval

1 1D heat and wave equations on a finite interval 1 1D het nd wve equtions on finite intervl In this section we consider generl method of seprtion of vribles nd its pplictions to solving het eqution nd wve eqution on finite intervl ( 1, 2. Since by trnsltion

More information

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies Stte spce systems nlysis (continued) Stbility A. Definitions A system is sid to be Asymptoticlly Stble (AS) when it stisfies ut () = 0, t > 0 lim xt () 0. t A system is AS if nd only if the impulse response

More information

Definite integral. Mathematics FRDIS MENDELU

Definite integral. Mathematics FRDIS MENDELU Definite integrl Mthemtics FRDIS MENDELU Simon Fišnrová Brno 1 Motivtion - re under curve Suppose, for simplicity, tht y = f(x) is nonnegtive nd continuous function defined on [, b]. Wht is the re of the

More information

Taylor Polynomial Inequalities

Taylor Polynomial Inequalities Tylor Polynomil Inequlities Ben Glin September 17, 24 Abstrct There re instnces where we my wish to pproximte the vlue of complicted function round given point by constructing simpler function such s polynomil

More information

APPLICATIONS OF THE DEFINITE INTEGRAL

APPLICATIONS OF THE DEFINITE INTEGRAL APPLICATIONS OF THE DEFINITE INTEGRAL. Volume: Slicing, disks nd wshers.. Volumes by Slicing. Suppose solid object hs boundries extending from x =, to x = b, nd tht its cross-section in plne pssing through

More information

20 MATHEMATICS POLYNOMIALS

20 MATHEMATICS POLYNOMIALS 0 MATHEMATICS POLYNOMIALS.1 Introduction In Clss IX, you hve studied polynomils in one vrible nd their degrees. Recll tht if p(x) is polynomil in x, the highest power of x in p(x) is clled the degree of

More information

Problem set 1: Solutions Math 207B, Winter 2016

Problem set 1: Solutions Math 207B, Winter 2016 Problem set 1: Solutions Mth 27B, Winter 216 1. Define f : R 2 R by f(,) = nd f(x,y) = xy3 x 2 +y 6 if (x,y) (,). ()Show tht thedirectionl derivtives of f t (,)exist inevery direction. Wht is its Gâteux

More information

Definite integral. Mathematics FRDIS MENDELU. Simona Fišnarová (Mendel University) Definite integral MENDELU 1 / 30

Definite integral. Mathematics FRDIS MENDELU. Simona Fišnarová (Mendel University) Definite integral MENDELU 1 / 30 Definite integrl Mthemtics FRDIS MENDELU Simon Fišnrová (Mendel University) Definite integrl MENDELU / Motivtion - re under curve Suppose, for simplicity, tht y = f(x) is nonnegtive nd continuous function

More information

Indefinite Integral. Chapter Integration - reverse of differentiation

Indefinite Integral. Chapter Integration - reverse of differentiation Chpter Indefinite Integrl Most of the mthemticl opertions hve inverse opertions. The inverse opertion of differentition is clled integrtion. For exmple, describing process t the given moment knowing the

More information

f(x) dx, If one of these two conditions is not met, we call the integral improper. Our usual definition for the value for the definite integral

f(x) dx, If one of these two conditions is not met, we call the integral improper. Our usual definition for the value for the definite integral Improper Integrls Every time tht we hve evluted definite integrl such s f(x) dx, we hve mde two implicit ssumptions bout the integrl:. The intervl [, b] is finite, nd. f(x) is continuous on [, b]. If one

More information

Chapter 5 Bending Moments and Shear Force Diagrams for Beams

Chapter 5 Bending Moments and Shear Force Diagrams for Beams Chpter 5 ending Moments nd Sher Force Digrms for ems n ddition to illy loded brs/rods (e.g. truss) nd torsionl shfts, the structurl members my eperience some lods perpendiculr to the is of the bem nd will

More information

Math 100 Review Sheet

Math 100 Review Sheet Mth 100 Review Sheet Joseph H. Silvermn December 2010 This outline of Mth 100 is summry of the mteril covered in the course. It is designed to be study id, but it is only n outline nd should be used s

More information

(uv) = u v + uv, (1) u vdx + b [uv] b a = u vdx + u v dx. (8) u vds =

(uv) = u v + uv, (1) u vdx + b [uv] b a = u vdx + u v dx. (8) u vds = Integrtion by prts Integrting the both sides of yields (uv) u v + uv, (1) or b (uv) dx b u vdx + b uv dx, (2) or b [uv] b u vdx + Eqution (4) is the 1-D formul for integrtion by prts. Eqution (4) cn be

More information

Numerical Integration

Numerical Integration Chpter 5 Numericl Integrtion Numericl integrtion is the study of how the numericl vlue of n integrl cn be found. Methods of function pproximtion discussed in Chpter??, i.e., function pproximtion vi the

More information

Math 113 Exam 2 Practice

Math 113 Exam 2 Practice Mth 3 Exm Prctice Februry 8, 03 Exm will cover 7.4, 7.5, 7.7, 7.8, 8.-3 nd 8.5. Plese note tht integrtion skills lerned in erlier sections will still be needed for the mteril in 7.5, 7.8 nd chpter 8. This

More information

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions Physics 6C Solution of inhomogeneous ordinry differentil equtions using Green s functions Peter Young November 5, 29 Homogeneous Equtions We hve studied, especilly in long HW problem, second order liner

More information

Best Approximation. Chapter The General Case

Best Approximation. Chapter The General Case Chpter 4 Best Approximtion 4.1 The Generl Cse In the previous chpter, we hve seen how n interpolting polynomil cn be used s n pproximtion to given function. We now wnt to find the best pproximtion to given

More information

Multiple Integrals. Review of Single Integrals. Planar Area. Volume of Solid of Revolution

Multiple Integrals. Review of Single Integrals. Planar Area. Volume of Solid of Revolution Multiple Integrls eview of Single Integrls eding Trim 7.1 eview Appliction of Integrls: Are 7. eview Appliction of Integrls: olumes 7.3 eview Appliction of Integrls: Lengths of Curves Assignment web pge

More information

Method of Localisation and Controlled Ejection of Swarms of Likely Charged Particles

Method of Localisation and Controlled Ejection of Swarms of Likely Charged Particles Method of Loclistion nd Controlled Ejection of Swrms of Likely Chrged Prticles I. N. Tukev July 3, 17 Astrct This work considers Coulom forces cting on chrged point prticle locted etween the two coxil,

More information

4. Calculus of Variations

4. Calculus of Variations 4. Clculus of Vritions Introduction - Typicl Problems The clculus of vritions generlises the theory of mxim nd minim. Exmple (): Shortest distnce between two points. On given surfce (e.g. plne), nd the

More information

Math 520 Final Exam Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008

Math 520 Final Exam Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008 Mth 520 Finl Exm Topic Outline Sections 1 3 (Xio/Dums/Liw) Spring 2008 The finl exm will be held on Tuesdy, My 13, 2-5pm in 117 McMilln Wht will be covered The finl exm will cover the mteril from ll of

More information

Chapter 4 Contravariance, Covariance, and Spacetime Diagrams

Chapter 4 Contravariance, Covariance, and Spacetime Diagrams Chpter 4 Contrvrince, Covrince, nd Spcetime Digrms 4. The Components of Vector in Skewed Coordintes We hve seen in Chpter 3; figure 3.9, tht in order to show inertil motion tht is consistent with the Lorentz

More information

3 Mathematics of the Poisson Equation

3 Mathematics of the Poisson Equation 3 Mthemtics of the Poisson Eqution 3. Green functions nd the Poisson eqution () The Dirichlet Green function stisfies the Poisson eqution with delt-function chrge 2 G D (r, r o ) = δ 3 (r r o ) (3.) nd

More information

Kinematic Waves. These are waves which result from the conservation equation. t + I = 0. (2)

Kinematic Waves. These are waves which result from the conservation equation. t + I = 0. (2) Introduction Kinemtic Wves These re wves which result from the conservtion eqution E t + I = 0 (1) where E represents sclr density field nd I, its outer flux. The one-dimensionl form of (1) is E t + I

More information

Chapter 3 Polynomials

Chapter 3 Polynomials Dr M DRAIEF As described in the introduction of Chpter 1, pplictions of solving liner equtions rise in number of different settings In prticulr, we will in this chpter focus on the problem of modelling

More information

Calculus of Variations

Calculus of Variations Clculus of Vritions Com S 477/577 Notes) Yn-Bin Ji Dec 4, 2017 1 Introduction A functionl ssigns rel number to ech function or curve) in some clss. One might sy tht functionl is function of nother function

More information

Introduction to Finite Element Method

Introduction to Finite Element Method Introduction to Finite Element Method Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pn.pl/ tzielins/ Tble of Contents 1 Introduction 1 1.1 Motivtion nd generl concepts.............

More information

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by. NUMERICAL INTEGRATION 1 Introduction The inverse process to differentition in clculus is integrtion. Mthemticlly, integrtion is represented by f(x) dx which stnds for the integrl of the function f(x) with

More information

CHM Physical Chemistry I Chapter 1 - Supplementary Material

CHM Physical Chemistry I Chapter 1 - Supplementary Material CHM 3410 - Physicl Chemistry I Chpter 1 - Supplementry Mteril For review of some bsic concepts in mth, see Atkins "Mthemticl Bckground 1 (pp 59-6), nd "Mthemticl Bckground " (pp 109-111). 1. Derivtion

More information

Plane curvilinear motion is the motion of a particle along a curved path which lies in a single plane.

Plane curvilinear motion is the motion of a particle along a curved path which lies in a single plane. Plne curiliner motion is the motion of prticle long cured pth which lies in single plne. Before the description of plne curiliner motion in n specific set of coordintes, we will use ector nlsis to describe

More information

Theoretical foundations of Gaussian quadrature

Theoretical foundations of Gaussian quadrature Theoreticl foundtions of Gussin qudrture 1 Inner product vector spce Definition 1. A vector spce (or liner spce) is set V = {u, v, w,...} in which the following two opertions re defined: (A) Addition of

More information

Fundamental Theorem of Calculus

Fundamental Theorem of Calculus Fundmentl Theorem of Clculus Recll tht if f is nonnegtive nd continuous on [, ], then the re under its grph etween nd is the definite integrl A= f() d Now, for in the intervl [, ], let A() e the re under

More information

Lecture 13 - Linking E, ϕ, and ρ

Lecture 13 - Linking E, ϕ, and ρ Lecture 13 - Linking E, ϕ, nd ρ A Puzzle... Inner-Surfce Chrge Density A positive point chrge q is locted off-center inside neutrl conducting sphericl shell. We know from Guss s lw tht the totl chrge on

More information

SOLUTION OF QUADRATIC NONLINEAR PROBLEMS WITH MULTIPLE SCALES LINDSTEDT-POINCARE METHOD. Mehmet Pakdemirli and Gözde Sarı

SOLUTION OF QUADRATIC NONLINEAR PROBLEMS WITH MULTIPLE SCALES LINDSTEDT-POINCARE METHOD. Mehmet Pakdemirli and Gözde Sarı Mthemticl nd Computtionl Applictions, Vol., No., pp. 37-5, 5 http://dx.doi.org/.99/mc-5- SOLUTION OF QUADRATIC NONLINEAR PROBLEMS WITH MULTIPLE SCALES LINDSTEDT-POINCARE METHOD Mehmet Pkdemirli nd Gözde

More information

Best Approximation in the 2-norm

Best Approximation in the 2-norm Jim Lmbers MAT 77 Fll Semester 1-11 Lecture 1 Notes These notes correspond to Sections 9. nd 9.3 in the text. Best Approximtion in the -norm Suppose tht we wish to obtin function f n (x) tht is liner combintion

More information

Chapter 4. Additional Variational Concepts

Chapter 4. Additional Variational Concepts Chpter 4 Additionl Vritionl Concepts 137 In the previous chpter we considered clculus o vrition problems which hd ixed boundry conditions. Tht is, in one dimension the end point conditions were speciied.

More information

Section 6.1 INTRO to LAPLACE TRANSFORMS

Section 6.1 INTRO to LAPLACE TRANSFORMS Section 6. INTRO to LAPLACE TRANSFORMS Key terms: Improper Integrl; diverge, converge A A f(t)dt lim f(t)dt Piecewise Continuous Function; jump discontinuity Function of Exponentil Order Lplce Trnsform

More information

Green function and Eigenfunctions

Green function and Eigenfunctions Green function nd Eigenfunctions Let L e regulr Sturm-Liouville opertor on n intervl (, ) together with regulr oundry conditions. We denote y, φ ( n, x ) the eigenvlues nd corresponding normlized eigenfunctions

More information

Integration Techniques

Integration Techniques Integrtion Techniques. Integrtion of Trigonometric Functions Exmple. Evlute cos x. Recll tht cos x = cos x. Hence, cos x Exmple. Evlute = ( + cos x) = (x + sin x) + C = x + 4 sin x + C. cos 3 x. Let u

More information

KINEMATICS OF RIGID BODIES

KINEMATICS OF RIGID BODIES KINEMTICS OF RIGI OIES Introduction In rigid body kinemtics, e use the reltionships governing the displcement, velocity nd ccelertion, but must lso ccount for the rottionl motion of the body. escription

More information

Chapter 28. Fourier Series An Eigenvalue Problem.

Chapter 28. Fourier Series An Eigenvalue Problem. Chpter 28 Fourier Series Every time I close my eyes The noise inside me mplifies I cn t escpe I relive every moment of the dy Every misstep I hve mde Finds wy it cn invde My every thought And this is why

More information

Families of Solutions to Bernoulli ODEs

Families of Solutions to Bernoulli ODEs In the fmily of solutions to the differentil eqution y ry dx + = it is shown tht vrition of the initil condition y( 0 = cuses horizontl shift in the solution curve y = f ( x, rther thn the verticl shift

More information

Generalizations of the Basic Functional

Generalizations of the Basic Functional 3 Generliztions of the Bsic Functionl 3 1 Chpter 3: GENERALIZATIONS OF THE BASIC FUNCTIONAL TABLE OF CONTENTS Pge 3.1 Functionls with Higher Order Derivtives.......... 3 3 3.2 Severl Dependent Vribles...............

More information

(6.5) Length and area in polar coordinates

(6.5) Length and area in polar coordinates 86 Chpter 6 SLICING TECHNIQUES FURTHER APPLICATIONS Totl mss 6 x ρ(x)dx + x 6 x dx + 9 kg dx + 6 x dx oment bout origin 6 xρ(x)dx x x dx + x + x + ln x ( ) + ln 6 kg m x dx + 6 6 x x dx Centre of mss +

More information

HOMEWORK SOLUTIONS MATH 1910 Sections 7.9, 8.1 Fall 2016

HOMEWORK SOLUTIONS MATH 1910 Sections 7.9, 8.1 Fall 2016 HOMEWORK SOLUTIONS MATH 9 Sections 7.9, 8. Fll 6 Problem 7.9.33 Show tht for ny constnts M,, nd, the function yt) = )) t ) M + tnh stisfies the logistic eqution: y SOLUTION. Let Then nd Finlly, y = y M

More information

The Wave Equation I. MA 436 Kurt Bryan

The Wave Equation I. MA 436 Kurt Bryan 1 Introduction The Wve Eqution I MA 436 Kurt Bryn Consider string stretching long the x xis, of indeterminte (or even infinite!) length. We wnt to derive n eqution which models the motion of the string

More information

E S dition event Vector Mechanics for Engineers: Dynamics h Due, next Wednesday, 07/19/2006! 1-30

E S dition event Vector Mechanics for Engineers: Dynamics h Due, next Wednesday, 07/19/2006! 1-30 Vector Mechnics for Engineers: Dynmics nnouncement Reminders Wednesdy s clss will strt t 1:00PM. Summry of the chpter 11 ws posted on website nd ws sent you by emil. For the students, who needs hrdcopy,

More information

Jackson 2.7 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell

Jackson 2.7 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell Jckson.7 Homework Problem Solution Dr. Christopher S. Bird University of Msschusetts Lowell PROBLEM: Consider potentil problem in the hlf-spce defined by, with Dirichlet boundry conditions on the plne

More information

Solutions to Homework Set 3

Solutions to Homework Set 3 Solutions to Homework Set 3 1 Bent brs: First recll some elementry clculus results: the curvture κ, nd the rdius R of the osculting circle, t point (x, y on curve y(x re given by κ = 1 R = dψ = 1 d tn

More information

Review of Calculus, cont d

Review of Calculus, cont d Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some

More information

u t = k 2 u x 2 (1) a n sin nπx sin 2 L e k(nπ/l) t f(x) = sin nπx f(x) sin nπx dx (6) 2 L f(x 0 ) sin nπx 0 2 L sin nπx 0 nπx

u t = k 2 u x 2 (1) a n sin nπx sin 2 L e k(nπ/l) t f(x) = sin nπx f(x) sin nπx dx (6) 2 L f(x 0 ) sin nπx 0 2 L sin nπx 0 nπx Chpter 9: Green s functions for time-independent problems Introductory emples One-dimensionl het eqution Consider the one-dimensionl het eqution with boundry conditions nd initil condition We lredy know

More information

1 The Lagrange interpolation formula

1 The Lagrange interpolation formula Notes on Qudrture 1 The Lgrnge interpoltion formul We briefly recll the Lgrnge interpoltion formul. The strting point is collection of N + 1 rel points (x 0, y 0 ), (x 1, y 1 ),..., (x N, y N ), with x

More information

Polynomials and Division Theory

Polynomials and Division Theory Higher Checklist (Unit ) Higher Checklist (Unit ) Polynomils nd Division Theory Skill Achieved? Know tht polynomil (expression) is of the form: n x + n x n + n x n + + n x + x + 0 where the i R re the

More information

Problems for HW X. C. Gwinn. November 30, 2009

Problems for HW X. C. Gwinn. November 30, 2009 Problems for HW X C. Gwinn November 30, 2009 These problems will not be grded. 1 HWX Problem 1 Suppose thn n object is composed of liner dielectric mteril, with constnt reltive permittivity ɛ r. The object

More information

6.5 Plate Problems in Rectangular Coordinates

6.5 Plate Problems in Rectangular Coordinates 6.5 lte rolems in Rectngulr Coordintes In this section numer of importnt plte prolems ill e emined ug Crte coordintes. 6.5. Uniform ressure producing Bending in One irection Consider first the cse of plte

More information

CMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature

CMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature CMDA 4604: Intermedite Topics in Mthemticl Modeling Lecture 19: Interpoltion nd Qudrture In this lecture we mke brief diversion into the res of interpoltion nd qudrture. Given function f C[, b], we sy

More information

Theoretische Physik 2: Elektrodynamik (Prof. A.-S. Smith) Home assignment 4

Theoretische Physik 2: Elektrodynamik (Prof. A.-S. Smith) Home assignment 4 WiSe 1 8.1.1 Prof. Dr. A.-S. Smith Dipl.-Phys. Ellen Fischermeier Dipl.-Phys. Mtthis Sb m Lehrstuhl für Theoretische Physik I Deprtment für Physik Friedrich-Alexnder-Universität Erlngen-Nürnberg Theoretische

More information

Jim Lambers MAT 280 Spring Semester Lecture 26 and 27 Notes

Jim Lambers MAT 280 Spring Semester Lecture 26 and 27 Notes Jim Lmbers MAT 280 pring emester 2009-10 Lecture 26 nd 27 Notes These notes correspond to ection 8.6 in Mrsden nd Tromb. ifferentil Forms To dte, we hve lerned the following theorems concerning the evlution

More information

12 TRANSFORMING BIVARIATE DENSITY FUNCTIONS

12 TRANSFORMING BIVARIATE DENSITY FUNCTIONS 1 TRANSFORMING BIVARIATE DENSITY FUNCTIONS Hving seen how to trnsform the probbility density functions ssocited with single rndom vrible, the next logicl step is to see how to trnsform bivrite probbility

More information

3.4 Numerical integration

3.4 Numerical integration 3.4. Numericl integrtion 63 3.4 Numericl integrtion In mny economic pplictions it is necessry to compute the definite integrl of relvlued function f with respect to "weight" function w over n intervl [,

More information

2. VECTORS AND MATRICES IN 3 DIMENSIONS

2. VECTORS AND MATRICES IN 3 DIMENSIONS 2 VECTORS AND MATRICES IN 3 DIMENSIONS 21 Extending the Theory of 2-dimensionl Vectors x A point in 3-dimensionl spce cn e represented y column vector of the form y z z-xis y-xis z x y x-xis Most of the

More information

PHYSICS 116C Homework 4 Solutions

PHYSICS 116C Homework 4 Solutions PHYSICS 116C Homework 4 Solutions 1. ( Simple hrmonic oscilltor. Clerly the eqution is of the Sturm-Liouville (SL form with λ = n 2, A(x = 1, B(x =, w(x = 1. Legendre s eqution. Clerly the eqution is of

More information

Geometric and Mechanical Applications of Integrals

Geometric and Mechanical Applications of Integrals 5 Geometric nd Mechnicl Applictions of Integrls 5.1 Computing Are 5.1.1 Using Crtesin Coordintes Suppose curve is given by n eqution y = f(x), x b, where f : [, b] R is continuous function such tht f(x)

More information

Hung problem # 3 April 10, 2011 () [4 pts.] The electric field points rdilly inwrd [1 pt.]. Since the chrge distribution is cylindriclly symmetric, we pick cylinder of rdius r for our Gussin surfce S.

More information

5.7 Improper Integrals

5.7 Improper Integrals 458 pplictions of definite integrls 5.7 Improper Integrls In Section 5.4, we computed the work required to lift pylod of mss m from the surfce of moon of mss nd rdius R to height H bove the surfce of the

More information

Density of Energy Stored in the Electric Field

Density of Energy Stored in the Electric Field Density of Energy Stored in the Electric Field Deprtment of Physics, Cornell University c Tomás A. Aris October 14, 01 Figure 1: Digrm of Crtesin vortices from René Descrtes Principi philosophie, published

More information

BME 207 Introduction to Biomechanics Spring 2018

BME 207 Introduction to Biomechanics Spring 2018 April 6, 28 UNIVERSITY O RHODE ISAND Deprtment of Electricl, Computer nd Biomedicl Engineering BME 27 Introduction to Biomechnics Spring 28 Homework 8 Prolem 14.6 in the textook. In ddition to prts -e,

More information

Conducting Ellipsoid and Circular Disk

Conducting Ellipsoid and Circular Disk 1 Problem Conducting Ellipsoid nd Circulr Disk Kirk T. McDonld Joseph Henry Lbortories, Princeton University, Princeton, NJ 08544 (September 1, 00) Show tht the surfce chrge density σ on conducting ellipsoid,

More information