Introduction to Finite Elements in Engineering. Tirupathi R. Chandrupatla Rowan University Glassboro, New Jersey

Size: px
Start display at page:

Download "Introduction to Finite Elements in Engineering. Tirupathi R. Chandrupatla Rowan University Glassboro, New Jersey"

Transcription

1 Introduction to Finite Elements in Engineering Tirupthi R. Chndruptl Rown Universit Glssboro, New Jerse Ashok D. Belegundu The Pennslvni Stte Universit Universit Prk, Pennslvni Solutions Mnul Prentice Hll, Upper Sddle River, New Jerse 758 Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

2 CONTENTS Prefce Chpter Fundmentl Concepts Chpter Mtri Algebr nd Gussin Elimintion 8 Chpter One-Dimensionl Problems 6 Chpter Trusses 6 Chpter 5 Bems nd Frmes 86 Chpter 6 Two-Dimensionl Problems Using Constnt Strin Tringles Chpter 7 Aismmetric Solids Subjected to 5 Aismmetric Loding Chpter 8 Two-Dimensionl Isoprmetric Elements 8 nd Numericl Integrtion Chpter 9 Three-Dimensionl Problems in Stress 7 Anlsis Chpter Sclr Field Problems 8 Chpter Dnmic Considertions 6 Chpter Preprocessing nd Postprocessing 8 Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

3 PREFACE This solutions mnul serves s n id to professors in teching from the book Introduction to Finite Elements in Engineering, th Edition. The problems in the book fll into the following ctegories:. Simple problems to understnd the concepts. Derivtions nd direct solutions. Solutions requiring computer runs. Solutions requiring progrm modifictions Our bsic philosoph in the development of this mnul is to provide complete guidnce to the techer in formulting, modeling, nd solving the problems. Complete solutions re given for problems in ll ctegories stted. For some lrger problems such s those in three dimensionl stress nlsis, complete formultion nd modeling spects re discussed. The students should be ble to proceed from the guidelines provided. For problems involving distributed nd other tpes of loding, the nodl lods re to be clculted for the input dt. The progrms do not generte the lods. This clcultion nd the boundr condition decisions enble the student to develop phsicl sense for the problems. The students m be encourged to modif the progrms to clculte the lods utomticll. The students should be introduced to the progrms in Chpter right from the point of solving problems in Chpter 6. This will enble the students to solve lrger problems with ese. The input dt file for ech progrm hs been provided. Dt for problem should follow this formt. The best strteg is to cop the emple file nd edit it for the problem under considertion. The dt from progrm MESHGEN will need some editing to complete the informtion on boundr conditions, lods, nd mteril properties. We thnk ou for our enthusistic response to our first three editions of the book. We look forwrd to receive our feedbck of our eperiences, comments, nd suggestions for mking improvements to the book nd this mnul. Tirupthi R. Chndruptl P.E., CMfgE Deprtment of Mechnicl Engineering Rown Universit, Glssboro, NJ 88 e-mil: chndruptl@rown.edu Ashok D. Belegundu Deprtment of Mechnicl nd Nucler Engineering The Pennslvni Stte Universit Universit Prk, PA 68 e-mil: db@psu.edu Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

4 CHAPTER FUNDAMENTAL CONCEPTS. We use the first three steps of Eq.. Adding the bove, we get ε ε ε σ σ σ ν ν E E E σ σ σ ν ν E E E σ σ σ ν ν E E E ε ε Adding nd subtrcting ( σ σ σ ) ν ε E σ ν from the first eqution, E ( σ σ σ ) ν ν ε σ E E Similr epressions cn be obtined for ε, nd ε. From the reltionship for γ nd Eq.., E τ γ etc. ( ν) Above reltions cn be written in the form σ Dε where D is the mteril propert mtri defined in Eq..5.. Note tht u () stisfies the ero slope boundr condition t the support. Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

5 . Plne strin condition implies tht σ σ ε ν ν E E which gives σ ν σ σ ( ) σ E 6 We hve, σ psi σ psi E psi ν. On substituting the vlues,. σ psi. Displcement field u v u v ε u t, ε ( 6) ( 6 ) ( 6) ( 6) u v v 6 9 u v ( 6 ).5 On inspection, we note tht the displcements u nd v re given b It is then es to see tht u. v Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

6 ε ε γ u v.6 The displcement field is given s u v. u 6 v 7 () The strins re then given b u ε 6 v ε u v γ (b) In order to drw the contours of the strin field using MATLAB, we need to crete script file, which m be edited s tet file nd sve with.m etension. The file for plotting ε is given below file probp5b.m [X,Y] meshgrid(-:.:,-:.:); Z..*X.^6.*Y.^; [C,h] contour(x,y,z); clbel(c,h); On running the progrm, the contour mp is shown s follows: Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

7 Contours of ε Contours of ε nd γ re obtined b chnging Z in the script file. The numbers on the contours show the function vlues. (c) The mimum vlue of ε is t n of the corners of the squre region. The mimum vlue is..7 (, ) (u, v) ). u u. v u v u v b) ε ε γ. Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

8 .8. σ τ σ MP MP n From Eq..8 we get T T T n σ 5.67 MP τ 6. MP τ.7 MP T n n n n τ σ τ σ τ.9 MP T n n n n MP 5 MP T τ τ σ T n.9 From the derivtion mde in P., we hve n n n σ τ MP MP σ σ E ( )( ) ( ν ) ν ν E [ ε νε νε ] which cn be written in the form ( )( ) [( ν ) ε ] νε v ν ν nd E τ γ ( ν) Lme s constnts λ nd µ re defined in the epressions σ τ λ λε µγ ( ν)( ν) E µ Eν ( ν) µε On inspection, v Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

9 µ is sme s the sher modulus G... ε ε. T E GP α ε α T σ E δ C du d L L 5-6 /.6 ( ε ε ) 69.6 MP du d C d L L. Following the steps of Emple., we hve ( 8 5) 8 8 q 8 q 6 5 Above mtri form is sme s the set of equtions: 7 q 8 q 6 8 q 8 q 5 Solving for q nd q, we get q. mm q.87 mm Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

10 . When the wll is smooth, σ. T is the temperture rise. ) When the block is thin in the direction, it corresponds to plne stress condition. The rigid wlls in the direction require ε. The generlied Hooke s lw ields the equtions σ ε ν α T E σ ε α T E From the second eqution, setting ε, we get σ Eα T. ε is then clculted using the first eqution s ( ν) α T. b) When the block is ver thick in the direction, plin strin condition previls. Now we hve ε, in ddition to ε. σ is not ero. σ σ ε ν ν α T E E σ σ ε ν α T E E σ σ ε ν α T E E From the lst two equtions, we get Eα T ν σ σ Eα T ν ν ε is now obtined from the first eqution. Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

11 . For thin block, it is plne stress condition. Treting the nominl sie s, we m set the. initil strin ε α T in prt () of problem.. Thus σ.e..5 The potentil energ Π is given b Π du d EA d ugad Consider the polnomil from Emple., u du d ( ) ( ) ( ) g E A On substituting the bove epressions nd integrting, the first term of becomes nd the second term Thus ugad ud ( ) Π Π u this gives ( ). 5.6 E A f Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

12 We use the displcement field defined b u. u t u t We then hve u ( ), nd du/d ( ). The potentil energ is now written s Π 6 du d ( ) d ( ) 5 ( ) d ( ) 5 6 d fud d d Π This ields,. Displcemen u.( ) Stress σ E du/d.( ).7 Let u be the displcement t mm. Piecewise liner displcement tht is continuous in the intervl 5 is represented s shown in the figure. u u u 5 Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

13 u t u u t u / u (u /) du/d u / 5 u t 5 5 u u t u u / (5/)u u (5/)u (u /) du/d u / Π E l du A d d u u Π El A Est A El A Est A u u 5 E st du A d d u u Π u E A E A u l st Note tht using the units MP (N/mm ) for modulus of elsticit nd mm for re nd mm for length will result in displcement in mm, nd stress in MP. Thus, E l 7 MP, E st, nd A 9 mm, A mm. On substituting these vlues into the bove eqution, we get u.9 mm This is precisel the solution obtined from strength of mterils pproch.8 In the Glerkin method, we strt from the equilibrium eqution d du EA g d d Following the steps of Emple., we get Introducing du EA d dφ d d gφd Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

14 ( ) ( ) nd, φ φ u u where u nd φ re the vlues of u nd φ t respectivel, ( ) ( ) φ d d u On integrting, we get 8 φ u This is to be stisfied for ever φ, which gives the solution u.5.9 We use t t u u u This implies tht 8 nd ( ) ( ) ( ) ( ) d du u nd re considered s independent vribles in ( ) ( ) [ ] ( ) Π d on epnding nd integrting the terms, we get Π We differentite with respect to the vribles nd equte to ero. Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

15 Π Π On solving, we get.7856 nd.5. On substituting in the epression for u, t, u.79 This pproimtion is close to the vlue obtined in the emple problem.. () Π L L T σ εad T ( )ud σ E ε nd ε On substitution, Π Π 6 du EA d du d d du ( 6 ) d ud ud d Tud 6 Tud (b) Since u t nd 6, nd u, we hve u du d ( 6) ( 6) On substituting nd integrting, Π Setting dπ/d gives Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

16 .5 du σ E 6.95( 6) d Plots of displcement nd stress re given below: Displcement u Stress.. t 6 implies tht 6 ( 8 ) 6, which ields Substituting for k, h, L, nd in I, we get I I I ( ) d ( 5) [ ( 8 ) 8] ( ) d 5( 8 9) Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

17 di d di d On solving, Substituting into the epression for, we get Since u t, the displcement stisfing the boundr condition is u. Also the coordintes re, nd. The potentil energ for the problem is du π EA d P u Pu d du We hve u, u, E, A, nd d. Thus π ( ) d. dπ For sttionr vlue, setting, we get d, which gives.75. The pproimte solution is u.75. Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

18 . Use Glerkin pproch with pproimtion u b c to solve du u d u ( ) The week form is obtined b multipling b φ stisfing φ ( ). du φ u d d We now set u b c into the bove integrl, stisfing u ( ) nd φ. On introducing these ( )( b c b c ) d ( ) ( ) b b c c d b b c c d On integrting, we get b c c b b c c b b c b c 6 This must be stisfied for ever nd. Thus the equtions to be solved re 7 7 b c 6 b c The solution is b.957, c.8. Thus u Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

19 P The deflection nd slope t due to P re nd EI nd slope t L due to lod P re ( ) P L P v EI EI P v EI The deflection nd slope due to lod P re PL v EI PL v EI We then get v v v v v v P EI. Using this the deflection.5 (,) (,) (,) (,) () The displcement of B is given b (.,.) nd A, C, nd D remin in their originl position. Consider displcement field of the tpe u v b b b b The four constnts cn be evluted using the known displcements Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

20 At A (, ) At B (, ) At C (, ) At D (, ) b. b b. b b b b b b The solution is,.,,. b, b., b, b u.. This gives v.. (b) The sher strin t B is u v γ... B ( ) ( ) γ.... Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

21 CHAPTER MATRIX ALGEBRA AND GAUSSIAN ELIMINATION., 8 d A () [ ] T d d I (b) det A 8 [ ()() ( ) ( ) ] ( ) [( )() ()( )] (c) The chrcteristic eqution is det (A - λ I ), or 8 det λ λ λ which ields 55 5 λ λ λ Hndbooks (e.g., CRC Mthemticl Hndbook) give eplicit solutions to cubic equtions. Here equtions given in Chpter 9 re used, which give formuls for finding the eigenvlues of the () smmetric stress tensor. Referring to Section 9. in the tet, we hve I A A A 5, I 55, I Thus,, b, c 5.6, θ 7. whence λ.5, λ 5.665, λ 9. Note: Since ll λ i >, A is positive definite. Now, eigenvector i corresponding to eigenvlue λ i is obtined from ( - λ i I ) i, i,, Thus, is obtined s Thus, Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

22 Onl two of the bove three equtions re independent. We hve Letting, we get [.7,.9, ] T The length of the vector is T. 8. Normliing to be unit vector ields [.7,.668,.7] T. Similrl, [.95,.577, -.65] T, [.85, -.7,.] T. (d) Solution to A b using Algorithm for generl mtri: n First Step (k ) i ( nd row) 8 A c / -/8 -/ () (-/)(-) 7/, () -, d () --(-/)() -/ i ( rd row) c () -, (), d () 8 () Thus () A 7 /, d / Second Step (k ) i ( rd row) c -6/7, () (-6/7)(-) /7, d () (-6/7)(-/)8/7 Thus A () 8 7 /, d / 7 () / 8/ 7 Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

23 Bck-Substitution 8 7 / / / 7 8/ 7 Third row gives: /7 8/7 whence 6 nd row then gives 5, st row gives.5 Thus, solution is [.5, 5, 6 ] T. Solution to A b using Algorithm for smmetric, bnded mtri 8 A is stored s n, nbw First Step (k ) nbk min (,) nd row (i ): i c / -/ j j (-/)(-) 7/ 8 Thus A () 7 / Second Step (k ) nbk rd row (i ): c -6/7 j j /7 8 Thus A () 7 / / 7 reduction of right-hnd-side vector d nd bck-substitution is sme s in Algorithm bove, resulting in the the sme solution [.5, 5, 6 ] T. Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

24 . N [ ξ ξ ] () (b) N dξ ξ ξ ( ξ ) T / N N dξ ( ) ( ) 6/ 5 ξ ξ ξ. q ( ) ( ) 6.5 T Q c T where.5 Q nd c.5 6. The detiled lgorithm is given in the tet. This is n ecercise in computer progrmming. The solutions re () (.5,.5,.5) (b) (.55, 5., 6.).5 A M M M det det, M i j C ij ( ) M ij The minors re, M X X X X.., M X X.... X X X. X., M, M, M Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

25 Thus, the co fctor mtri is C A nd A det A.6 Are det.7 A det A det A det A A A A.65 A / A.8, A / A.5, A / A.7 C T ields.8 A i,j B i, j-i for j i Thus, A, corresponds to B, nd B 6, A 6,6.9 Full () mtri BANDED n, nbw No. of storge loctions (n) (nbw) SKYLINE Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

26 No. of storge loctions no. of column entries... () () / 55. A 6 Cholesk decomposition follows the steps s per Eq. (.): k ( st row) k ( nd row), 6 (/ ). 969 k ( rd row) (/ )(/ ), Thus, L We m see tht A L L T. Also, progrm implementtion is given in Subroutine CHOLESKI within Progrm GENEIGEN.. B epnding the mtri equtions, we obtin set of simultneous equtions for the nd u coefficients, resulting in the solution 5 L., U Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

27 . D (, 7) C (6, 6) B (7, ) A (, ) We split the qudrilterl into two tringles s shown. Using the epression from Problem.6, the re is then given b Are det 7 det det 6 det ( 5) ( ).5 At the second step, we subtrct the first row from second row nd third row, which does not chnge the determinnt vlue.. Setting cosθ sinθ T sinθ cosθ, it is es to show tht T TT.. A 5 Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

28 () Minors M 7 M 9 M 5 5 M M 7 M M M 5 M (b) Cofctors C M 7 C M 9 C M C M C M 7 C M 7 C M C M 5 C M (c) Adjoint 7 7 T AdjA C (d) Determinnt deta (7)-(9)() (e) Inverse A 7 AdjA det A 7 Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

29 CHAPTER ONE-DIMENSIONAL PROBLEMS. q.in P q.5in 5in in in A. in E 6 psi () At point P, (b) N N u Nq N q ε in du d ε Bq q q B N q [ ] B [ ] 8. ε [ ] σ Eε 875 psi 8.75 ksi (c) k EA L 6.5 (d) U e ½ q T kq. U [..5].5 6 e 56.5 lb-in.5. NBW m [(-), (-), (5-), (5-)] Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

30 . (i) Plotting the displcement is stright forwrd.6 u mm -. 5 mm (ii) Strins re constnt in ech of the three elements Since the displcements re liner within the elements, the strins re constnt over ech element. Strin in Element -: L - 5 mm ε 5 (.). Strin in Element -: L - 8 mm.6 ε.75 8 Strin in Element -: L - mm (.6). ε.7 The plotted curve represents the slope of the curve shown in section (i). (iii) The B mtri for element - is given b (iv) B The stiffness mtri of element - is given b [ ] [.5.5] k EA L ( )( ) 5 N/mm Introduction to Finite Elements in Engineering, Fourth Edition, b T. R. Chndruptl nd A. D. Belegundu. ISBN Person Eduction, Inc., Upper Sddle River, NJ. All rights reserved. This publiction is protected b Copright write to: Rights nd Permissions Deprtment, Person Eduction, Inc., Upper Sddle River, NJ 758.

Introduction to Finite Elements in Engineering. Tirupathi R. Chandrupatla Rowan University Glassboro, New Jersey

Introduction to Finite Elements in Engineering. Tirupathi R. Chandrupatla Rowan University Glassboro, New Jersey Introduction to Finite lements in ngineering Tirupthi R. Chndruptl Rown Universit Glssboro, New Jerse Ashok D. Belegundu The Pennslvni Stte Universit Universit Prk, Pennslvni Solutions Mnul Prentice Hll,

More information

A - INTRODUCTION AND OVERVIEW

A - INTRODUCTION AND OVERVIEW MMJ5 COMPUTATIONAL METHOD IN SOLID MECHANICS A - INTRODUCTION AND OVERVIEW INTRODUCTION AND OVERVIEW M.N. Tmin, CSMLb, UTM MMJ5 COMPUTATIONAL METHOD IN SOLID MECHANICS Course Content: A INTRODUCTION AND

More information

The Algebra (al-jabr) of Matrices

The Algebra (al-jabr) of Matrices Section : Mtri lgebr nd Clculus Wshkewicz College of Engineering he lgebr (l-jbr) of Mtrices lgebr s brnch of mthemtics is much broder thn elementry lgebr ll of us studied in our high school dys. In sense

More information

The Islamic University of Gaza Faculty of Engineering Civil Engineering Department. Numerical Analysis ECIV Chapter 11

The Islamic University of Gaza Faculty of Engineering Civil Engineering Department. Numerical Analysis ECIV Chapter 11 The Islmic University of Gz Fculty of Engineering Civil Engineering Deprtment Numericl Anlysis ECIV 6 Chpter Specil Mtrices nd Guss-Siedel Associte Prof Mzen Abultyef Civil Engineering Deprtment, The Islmic

More information

Matrices and Determinants

Matrices and Determinants Nme Chpter 8 Mtrices nd Determinnts Section 8.1 Mtrices nd Systems of Equtions Objective: In this lesson you lerned how to use mtrices, Gussin elimintion, nd Guss-Jordn elimintion to solve systems of liner

More information

Matrix Eigenvalues and Eigenvectors September 13, 2017

Matrix Eigenvalues and Eigenvectors September 13, 2017 Mtri Eigenvlues nd Eigenvectors September, 7 Mtri Eigenvlues nd Eigenvectors Lrry Cretto Mechnicl Engineering 5A Seminr in Engineering Anlysis September, 7 Outline Review lst lecture Definition of eigenvlues

More information

INTRODUCTION TO LINEAR ALGEBRA

INTRODUCTION TO LINEAR ALGEBRA ME Applied Mthemtics for Mechnicl Engineers INTRODUCTION TO INEAR AGEBRA Mtrices nd Vectors Prof. Dr. Bülent E. Pltin Spring Sections & / ME Applied Mthemtics for Mechnicl Engineers INTRODUCTION TO INEAR

More information

Calculus - Activity 1 Rate of change of a function at a point.

Calculus - Activity 1 Rate of change of a function at a point. Nme: Clss: p 77 Mths Helper Plus Resource Set. Copright 00 Bruce A. Vughn, Techers Choice Softwre Clculus - Activit Rte of chnge of function t point. ) Strt Mths Helper Plus, then lod the file: Clculus

More information

LINEAR ALGEBRA APPLIED

LINEAR ALGEBRA APPLIED 5.5 Applictions of Inner Product Spces 5.5 Applictions of Inner Product Spces 7 Find the cross product of two vectors in R. Find the liner or qudrtic lest squres pproimtion of function. Find the nth-order

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

M344 - ADVANCED ENGINEERING MATHEMATICS

M344 - ADVANCED ENGINEERING MATHEMATICS M3 - ADVANCED ENGINEERING MATHEMATICS Lecture 18: Lplce s Eqution, Anltic nd Numericl Solution Our emple of n elliptic prtil differentil eqution is Lplce s eqution, lso clled the Diffusion Eqution. If

More information

Stress distribution in elastic isotropic semi-space with concentrated vertical force

Stress distribution in elastic isotropic semi-space with concentrated vertical force Bulgrin Chemicl Communictions Volume Specil Issue pp. 4 9 Stress distribution in elstic isotropic semi-spce with concentrted verticl force L. B. Petrov Deprtment of Mechnics Todor Kbleshkov Universit of

More information

1 Part II: Numerical Integration

1 Part II: Numerical Integration Mth 4 Lb 1 Prt II: Numericl Integrtion This section includes severl techniques for getting pproimte numericl vlues for definite integrls without using ntiderivtives. Mthemticll, ect nswers re preferble

More information

Some Methods in the Calculus of Variations

Some Methods in the Calculus of Variations CHAPTER 6 Some Methods in the Clculus of Vritions 6-. If we use the vried function ( α, ) α sin( ) + () Then d α cos ( ) () d Thus, the totl length of the pth is d S + d d α cos ( ) + α cos ( ) d Setting

More information

Chapter 6 Notes, Larson/Hostetler 3e

Chapter 6 Notes, Larson/Hostetler 3e Contents 6. Antiderivtives nd the Rules of Integrtion.......................... 6. Are nd the Definite Integrl.................................. 6.. Are............................................ 6. Reimnn

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

Section - 2 MORE PROPERTIES

Section - 2 MORE PROPERTIES LOCUS Section - MORE PROPERTES n section -, we delt with some sic properties tht definite integrls stisf. This section continues with the development of some more properties tht re not so trivil, nd, when

More information

Engineering Analysis ENG 3420 Fall Dan C. Marinescu Office: HEC 439 B Office hours: Tu-Th 11:00-12:00

Engineering Analysis ENG 3420 Fall Dan C. Marinescu Office: HEC 439 B Office hours: Tu-Th 11:00-12:00 Engineering Anlysis ENG 3420 Fll 2009 Dn C. Mrinescu Office: HEC 439 B Office hours: Tu-Th 11:00-12:00 Lecture 13 Lst time: Problem solving in preprtion for the quiz Liner Algebr Concepts Vector Spces,

More information

September 13 Homework Solutions

September 13 Homework Solutions College of Engineering nd Computer Science Mechnicl Engineering Deprtment Mechnicl Engineering 5A Seminr in Engineering Anlysis Fll Ticket: 5966 Instructor: Lrry Cretto Septemer Homework Solutions. Are

More information

MATH , Calculus 2, Fall 2018

MATH , Calculus 2, Fall 2018 MATH 36-2, 36-3 Clculus 2, Fll 28 The FUNdmentl Theorem of Clculus Sections 5.4 nd 5.5 This worksheet focuses on the most importnt theorem in clculus. In fct, the Fundmentl Theorem of Clculus (FTC is rgubly

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: Volumes 7.3 eview Appliction of Integrls: Lengths of Curves Assignment web pge

More information

ES.182A Topic 32 Notes Jeremy Orloff

ES.182A Topic 32 Notes Jeremy Orloff ES.8A Topic 3 Notes Jerem Orloff 3 Polr coordintes nd double integrls 3. Polr Coordintes (, ) = (r cos(θ), r sin(θ)) r θ Stndrd,, r, θ tringle Polr coordintes re just stndrd trigonometric reltions. In

More information

Plates on elastic foundation

Plates on elastic foundation 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) Introduction Winkler

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

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

AQA Further Pure 2. Hyperbolic Functions. Section 2: The inverse hyperbolic functions

AQA Further Pure 2. Hyperbolic Functions. Section 2: The inverse hyperbolic functions Hperbolic Functions Section : The inverse hperbolic functions Notes nd Emples These notes contin subsections on The inverse hperbolic functions Integrtion using the inverse hperbolic functions Logrithmic

More information

Module 1. Energy Methods in Structural Analysis

Module 1. Energy Methods in Structural Analysis Module 1 Energy Methods in Structurl Anlysis Lesson 4 Theorem of Lest Work Instructionl Objectives After reding this lesson, the reder will be ble to: 1. Stte nd prove theorem of Lest Work.. Anlyse stticlly

More information

DETERMINANTS. All Mathematical truths are relative and conditional. C.P. STEINMETZ

DETERMINANTS. All Mathematical truths are relative and conditional. C.P. STEINMETZ All Mthemticl truths re reltive nd conditionl. C.P. STEINMETZ 4. Introduction DETERMINANTS In the previous chpter, we hve studied bout mtrices nd lgebr of mtrices. We hve lso lernt tht system of lgebric

More information

Chapter 5 Determinants

Chapter 5 Determinants hpter 5 Determinnts 5. Introduction Every squre mtri hs ssocited with it sclr clled its determinnt. Given mtri, we use det() or to designte its determinnt. We cn lso designte the determinnt of mtri by

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

13.3. The Area Bounded by a Curve. Introduction. Prerequisites. Learning Outcomes

13.3. The Area Bounded by a Curve. Introduction. Prerequisites. Learning Outcomes The Are Bounded b Curve 3.3 Introduction One of the importnt pplictions of integrtion is to find the re bounded b curve. Often such n re cn hve phsicl significnce like the work done b motor, or the distnce

More information

Chapter 3 MATRIX. In this chapter: 3.1 MATRIX NOTATION AND TERMINOLOGY

Chapter 3 MATRIX. In this chapter: 3.1 MATRIX NOTATION AND TERMINOLOGY Chpter 3 MTRIX In this chpter: Definition nd terms Specil Mtrices Mtrix Opertion: Trnspose, Equlity, Sum, Difference, Sclr Multipliction, Mtrix Multipliction, Determinnt, Inverse ppliction of Mtrix in

More information

13.3. The Area Bounded by a Curve. Introduction. Prerequisites. Learning Outcomes

13.3. The Area Bounded by a Curve. Introduction. Prerequisites. Learning Outcomes The Are Bounded b Curve 3.3 Introduction One of the importnt pplictions of integrtion is to find the re bounded b curve. Often such n re cn hve phsicl significnce like the work done b motor, or the distnce

More information

Chapter Direct Method of Interpolation More Examples Electrical Engineering

Chapter Direct Method of Interpolation More Examples Electrical Engineering Chpter. Direct Method of Interpoltion More Emples Electricl Engineering Emple hermistors re used to mesure the temperture of bodies. hermistors re bsed on mterils chnge in resistnce with temperture. o

More information

Mathematics. Area under Curve.

Mathematics. Area under Curve. Mthemtics Are under Curve www.testprepkrt.com Tle of Content 1. Introduction.. Procedure of Curve Sketching. 3. Sketching of Some common Curves. 4. Are of Bounded Regions. 5. Sign convention for finding

More information

approaches as n becomes larger and larger. Since e > 1, the graph of the natural exponential function is as below

approaches as n becomes larger and larger. Since e > 1, the graph of the natural exponential function is as below . Eponentil nd rithmic functions.1 Eponentil Functions A function of the form f() =, > 0, 1 is clled n eponentil function. Its domin is the set of ll rel f ( 1) numbers. For n eponentil function f we hve.

More information

Operations with Polynomials

Operations with Polynomials 38 Chpter P Prerequisites P.4 Opertions with Polynomils Wht you should lern: How to identify the leding coefficients nd degrees of polynomils How to dd nd subtrct polynomils How to multiply polynomils

More information

ES.182A Topic 30 Notes Jeremy Orloff

ES.182A Topic 30 Notes Jeremy Orloff ES82A opic 3 Notes Jerem Orloff 3 Non-independent vribles: chin rule Recll the chin rule: If w = f, ; nd = r, t, = r, t then = + r t r t r t = + t t t r nfortuntel, sometimes there re more complicted reltions

More information

Quantum Physics II (8.05) Fall 2013 Assignment 2

Quantum Physics II (8.05) Fall 2013 Assignment 2 Quntum Physics II (8.05) Fll 2013 Assignment 2 Msschusetts Institute of Technology Physics Deprtment Due Fridy September 20, 2013 September 13, 2013 3:00 pm Suggested Reding Continued from lst week: 1.

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

We divide the interval [a, b] into subintervals of equal length x = b a n

We divide the interval [a, b] into subintervals of equal length x = b a n Arc Length Given curve C defined by function f(x), we wnt to find the length of this curve between nd b. We do this by using process similr to wht we did in defining the Riemnn Sum of definite integrl:

More information

Reference : Croft & Davison, Chapter 12, Blocks 1,2. A matrix ti is a rectangular array or block of numbers usually enclosed in brackets.

Reference : Croft & Davison, Chapter 12, Blocks 1,2. A matrix ti is a rectangular array or block of numbers usually enclosed in brackets. I MATRIX ALGEBRA INTRODUCTION TO MATRICES Referene : Croft & Dvison, Chpter, Blos, A mtri ti is retngulr rr or lo of numers usull enlosed in rets. A m n mtri hs m rows nd n olumns. Mtri Alger Pge If the

More information

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 2

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 2 CS434/54: Pttern Recognition Prof. Olg Veksler Lecture Outline Review of Liner Algebr vectors nd mtrices products nd norms vector spces nd liner trnsformtions eigenvlues nd eigenvectors Introduction to

More information

DETERMINANTS. All Mathematical truths are relative and conditional. C.P. STEINMETZ

DETERMINANTS. All Mathematical truths are relative and conditional. C.P. STEINMETZ DETERMINANTS Chpter 4 All Mthemticl truths re reltive nd conditionl. C.P. STEINMETZ 4. Introduction In the previous chpter, we hve studied bout mtrices nd lgebr of mtrices. We hve lso lernt tht sstem of

More information

3.1 Exponential Functions and Their Graphs

3.1 Exponential Functions and Their Graphs . Eponentil Functions nd Their Grphs Sllbus Objective: 9. The student will sketch the grph of eponentil, logistic, or logrithmic function. 9. The student will evlute eponentil or logrithmic epressions.

More information

Math 113 Exam 2 Practice

Math 113 Exam 2 Practice Mth Em Prctice Februry, 8 Em will cover sections 6.5, 7.-7.5 nd 7.8. This sheet hs three sections. The first section will remind you bout techniques nd formuls tht you should know. The second gives number

More information

The First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a).

The First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a). The Fundmentl Theorems of Clculus Mth 4, Section 0, Spring 009 We now know enough bout definite integrls to give precise formultions of the Fundmentl Theorems of Clculus. We will lso look t some bsic emples

More information

Now, given the derivative, can we find the function back? Can we antidifferenitate it?

Now, given the derivative, can we find the function back? Can we antidifferenitate it? Fundmentl Theorem of Clculus. Prt I Connection between integrtion nd differentition. Tody we will discuss reltionship between two mjor concepts of Clculus: integrtion nd differentition. We will show tht

More information

Lesson 1: Quadratic Equations

Lesson 1: Quadratic Equations Lesson 1: Qudrtic Equtions Qudrtic Eqution: The qudrtic eqution in form is. In this section, we will review 4 methods of qudrtic equtions, nd when it is most to use ech method. 1. 3.. 4. Method 1: Fctoring

More information

MATRICES AND VECTORS SPACE

MATRICES AND VECTORS SPACE MATRICES AND VECTORS SPACE MATRICES AND MATRIX OPERATIONS SYSTEM OF LINEAR EQUATIONS DETERMINANTS VECTORS IN -SPACE AND -SPACE GENERAL VECTOR SPACES INNER PRODUCT SPACES EIGENVALUES, EIGENVECTORS LINEAR

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

6.2 CONCEPTS FOR ADVANCED MATHEMATICS, C2 (4752) AS

6.2 CONCEPTS FOR ADVANCED MATHEMATICS, C2 (4752) AS 6. CONCEPTS FOR ADVANCED MATHEMATICS, C (475) AS Objectives To introduce students to number of topics which re fundmentl to the dvnced study of mthemtics. Assessment Emintion (7 mrks) 1 hour 30 minutes.

More information

Student Handbook for MATH 3300

Student Handbook for MATH 3300 Student Hndbook for MATH 3300 0.8 0.6 0.4 0.2 0 0.2 0.4 0.6 0.8 0.5 0 0.5 0.5 0 0.5 If people do not believe tht mthemtics is simple, it is only becuse they do not relize how complicted life is. John Louis

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

x = b a n x 2 e x dx. cdx = c(b a), where c is any constant. a b

x = b a n x 2 e x dx. cdx = c(b a), where c is any constant. a b CHAPTER 5. INTEGRALS 61 where nd x = b n x i = 1 (x i 1 + x i ) = midpoint of [x i 1, x i ]. Problem 168 (Exercise 1, pge 377). Use the Midpoint Rule with the n = 4 to pproximte 5 1 x e x dx. Some quick

More information

5.5 The Substitution Rule

5.5 The Substitution Rule 5.5 The Substitution Rule Given the usefulness of the Fundmentl Theorem, we wnt some helpful methods for finding ntiderivtives. At the moment, if n nti-derivtive is not esily recognizble, then we re in

More information

SCHOOL OF ENGINEERING & BUILT ENVIRONMENT. Mathematics

SCHOOL OF ENGINEERING & BUILT ENVIRONMENT. Mathematics SCHOOL OF ENGINEERING & BUIL ENVIRONMEN Mthemtics An Introduction to Mtrices Definition of Mtri Size of Mtri Rows nd Columns of Mtri Mtri Addition Sclr Multipliction of Mtri Mtri Multipliction 7 rnspose

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

Algebra Of Matrices & Determinants

Algebra Of Matrices & Determinants lgebr Of Mtrices & Determinnts Importnt erms Definitions & Formule 0 Mtrix - bsic introduction: mtrix hving m rows nd n columns is clled mtrix of order m n (red s m b n mtrix) nd mtrix of order lso in

More information

The Trapezoidal Rule

The Trapezoidal Rule _.qd // : PM Pge 9 SECTION. Numericl Integrtion 9 f Section. The re of the region cn e pproimted using four trpezoids. Figure. = f( ) f( ) n The re of the first trpezoid is f f n. Figure. = Numericl Integrtion

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

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

Topic 1 Notes Jeremy Orloff

Topic 1 Notes Jeremy Orloff Topic 1 Notes Jerem Orloff 1 Introduction to differentil equtions 1.1 Gols 1. Know the definition of differentil eqution. 2. Know our first nd second most importnt equtions nd their solutions. 3. Be ble

More information

1.1. Linear Constant Coefficient Equations. Remark: A differential equation is an equation

1.1. Linear Constant Coefficient Equations. Remark: A differential equation is an equation 1 1.1. Liner Constnt Coefficient Equtions Section Objective(s): Overview of Differentil Equtions. Liner Differentil Equtions. Solving Liner Differentil Equtions. The Initil Vlue Problem. 1.1.1. Overview

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

4.6 Numerical Integration

4.6 Numerical Integration .6 Numericl Integrtion 5.6 Numericl Integrtion Approimte definite integrl using the Trpezoidl Rule. Approimte definite integrl using Simpson s Rule. Anlze the pproimte errors in the Trpezoidl Rule nd Simpson

More information

PROPERTIES OF AREAS In general, and for an irregular shape, the definition of the centroid at position ( x, y) is given by

PROPERTIES OF AREAS In general, and for an irregular shape, the definition of the centroid at position ( x, y) is given by PROPERTES OF RES Centroid The concept of the centroid is prol lred fmilir to ou For plne shpe with n ovious geometric centre, (rectngle, circle) the centroid is t the centre f n re hs n is of smmetr, the

More information

ME 501A Seminar in Engineering Analysis Page 1

ME 501A Seminar in Engineering Analysis Page 1 Phse-plne Anlsis of Ordinr November, 7 Phse-plne Anlsis of Ordinr Lrr Cretto Mechnicl Engineering 5A Seminr in Engineering Anlsis November, 7 Outline Mierm exm two weeks from tonight covering ODEs nd Lplce

More information

STRAND J: TRANSFORMATIONS, VECTORS and MATRICES

STRAND J: TRANSFORMATIONS, VECTORS and MATRICES Mthemtics SKE: STRN J STRN J: TRNSFORMTIONS, VETORS nd MTRIES J3 Vectors Text ontents Section J3.1 Vectors nd Sclrs * J3. Vectors nd Geometry Mthemtics SKE: STRN J J3 Vectors J3.1 Vectors nd Sclrs Vectors

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

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

Duality # Second iteration for HW problem. Recall our LP example problem we have been working on, in equality form, is given below.

Duality # Second iteration for HW problem. Recall our LP example problem we have been working on, in equality form, is given below. Dulity #. Second itertion for HW problem Recll our LP emple problem we hve been working on, in equlity form, is given below.,,,, 8 m F which, when written in slightly different form, is 8 F Recll tht we

More information

Session Trimester 2. Module Code: MATH08001 MATHEMATICS FOR DESIGN

Session Trimester 2. Module Code: MATH08001 MATHEMATICS FOR DESIGN School of Science & Sport Pisley Cmpus Session 05-6 Trimester Module Code: MATH0800 MATHEMATICS FOR DESIGN Dte: 0 th My 06 Time: 0.00.00 Instructions to Cndidtes:. Answer ALL questions in Section A. Section

More information

Calculus 2: Integration. Differentiation. Integration

Calculus 2: Integration. Differentiation. Integration Clculus 2: Integrtion The reverse process to differentition is known s integrtion. Differentition f() f () Integrtion As it is the opposite of finding the derivtive, the function obtined b integrtion is

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

10.2 The Ellipse and the Hyperbola

10.2 The Ellipse and the Hyperbola CHAPTER 0 Conic Sections Solve. 97. Two surveors need to find the distnce cross lke. The plce reference pole t point A in the digrm. Point B is meters est nd meter north of the reference point A. Point

More information

Higher Checklist (Unit 3) Higher Checklist (Unit 3) Vectors

Higher Checklist (Unit 3) Higher Checklist (Unit 3) Vectors Vectors Skill Achieved? Know tht sclr is quntity tht hs only size (no direction) Identify rel-life exmples of sclrs such s, temperture, mss, distnce, time, speed, energy nd electric chrge Know tht vector

More information

Review Exercises for Chapter 4

Review Exercises for Chapter 4 _R.qd // : PM Pge CHAPTER Integrtion Review Eercises for Chpter In Eercises nd, use the grph of to sketch grph of f. To print n enlrged cop of the grph, go to the wesite www.mthgrphs.com... In Eercises

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

SCHOOL OF ENGINEERING & BUILT ENVIRONMENT

SCHOOL OF ENGINEERING & BUILT ENVIRONMENT SCHOOL OF ENGINEERING & BUIL ENVIRONMEN MARICES FOR ENGINEERING Dr Clum Mcdonld Contents Introduction Definitions Wht is mtri? Rows nd columns of mtri Order of mtri Element of mtri Equlity of mtrices Opertions

More information

SCHOOL OF ENGINEERING & BUILT ENVIRONMENT. Mathematics. Basic Algebra

SCHOOL OF ENGINEERING & BUILT ENVIRONMENT. Mathematics. Basic Algebra SCHOOL OF ENGINEERING & BUILT ENVIRONMENT Mthemtics Bsic Algebr Opertions nd Epressions Common Mistkes Division of Algebric Epressions Eponentil Functions nd Logrithms Opertions nd their Inverses Mnipulting

More information

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

(4.1) D r v(t) ω(t, v(t))

(4.1) D r v(t) ω(t, v(t)) 1.4. Differentil inequlities. Let D r denote the right hnd derivtive of function. If ω(t, u) is sclr function of the sclrs t, u in some open connected set Ω, we sy tht function v(t), t < b, is solution

More information

A-Level Mathematics Transition Task (compulsory for all maths students and all further maths student)

A-Level Mathematics Transition Task (compulsory for all maths students and all further maths student) A-Level Mthemtics Trnsition Tsk (compulsory for ll mths students nd ll further mths student) Due: st Lesson of the yer. Length: - hours work (depending on prior knowledge) This trnsition tsk provides revision

More information

Maths in Motion. Theo de Haan. Order now: 29,95 euro

Maths in Motion. Theo de Haan. Order now:   29,95 euro Mths in Motion Theo de Hn Order now: www.mthsinmotion.org 9,95 euro Cover Design: Drwings: Photogrph: Printing: Niko Spelbrink Lr Wgterveld Mrijke Spelbrink Rddrier, Amsterdm Preview: Prts of Chpter 6,

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

Quotient Rule: am a n = am n (a 0) Negative Exponents: a n = 1 (a 0) an Power Rules: (a m ) n = a m n (ab) m = a m b m

Quotient Rule: am a n = am n (a 0) Negative Exponents: a n = 1 (a 0) an Power Rules: (a m ) n = a m n (ab) m = a m b m Formuls nd Concepts MAT 099: Intermedite Algebr repring for Tests: The formuls nd concepts here m not be inclusive. You should first tke our prctice test with no notes or help to see wht mteril ou re comfortble

More information

Section 4: Integration ECO4112F 2011

Section 4: Integration ECO4112F 2011 Reding: Ching Chpter Section : Integrtion ECOF Note: These notes do not fully cover the mteril in Ching, ut re ment to supplement your reding in Ching. Thus fr the optimistion you hve covered hs een sttic

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

1. Extend QR downwards to meet the x-axis at U(6, 0). y

1. Extend QR downwards to meet the x-axis at U(6, 0). y In the digrm, two stright lines re to be drwn through so tht the lines divide the figure OPQRST into pieces of equl re Find the sum of the slopes of the lines R(6, ) S(, ) T(, 0) Determine ll liner functions

More information

50. Use symmetry to evaluate xx D is the region bounded by the square with vertices 5, Prove Property 11. y y CAS

50. Use symmetry to evaluate xx D is the region bounded by the square with vertices 5, Prove Property 11. y y CAS 68 CHAPTE MULTIPLE INTEGALS 46. e da, 49. Evlute tn 3 4 da, where,. [Hint: Eploit the fct tht is the disk with center the origin nd rdius is smmetric with respect to both es.] 5. Use smmetr to evlute 3

More information

Counting intersections of spirals on a torus

Counting intersections of spirals on a torus Counting intersections of spirls on torus 1 The problem Consider unit squre with opposite sides identified. For emple, if we leve the centre of the squre trveling long line of slope 2 (s shown in the first

More information

Chapter 5. , r = r 1 r 2 (1) µ = m 1 m 2. r, r 2 = R µ m 2. R(m 1 + m 2 ) + m 2 r = r 1. m 2. r = r 1. R + µ m 1

Chapter 5. , r = r 1 r 2 (1) µ = m 1 m 2. r, r 2 = R µ m 2. R(m 1 + m 2 ) + m 2 r = r 1. m 2. r = r 1. R + µ m 1 Tor Kjellsson Stockholm University Chpter 5 5. Strting with the following informtion: R = m r + m r m + m, r = r r we wnt to derive: µ = m m m + m r = R + µ m r, r = R µ m r 3 = µ m R + r, = µ m R r. 4

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

In-Class Problems 2 and 3: Projectile Motion Solutions. In-Class Problem 2: Throwing a Stone Down a Hill

In-Class Problems 2 and 3: Projectile Motion Solutions. In-Class Problem 2: Throwing a Stone Down a Hill MASSACHUSETTS INSTITUTE OF TECHNOLOGY Deprtment of Physics Physics 8T Fll Term 4 In-Clss Problems nd 3: Projectile Motion Solutions We would like ech group to pply the problem solving strtegy with the

More information

ARITHMETIC OPERATIONS. The real numbers have the following properties: a b c ab ac

ARITHMETIC OPERATIONS. The real numbers have the following properties: a b c ab ac REVIEW OF ALGEBRA Here we review the bsic rules nd procedures of lgebr tht you need to know in order to be successful in clculus. ARITHMETIC OPERATIONS The rel numbers hve the following properties: b b

More information

CET MATHEMATICS 2013

CET MATHEMATICS 2013 CET MATHEMATICS VERSION CODE: C. If sin is the cute ngle between the curves + nd + 8 t (, ), then () () () Ans: () Slope of first curve m ; slope of second curve m - therefore ngle is o A sin o (). The

More information

Chapter 3 Exponential and Logarithmic Functions Section 3.1

Chapter 3 Exponential and Logarithmic Functions Section 3.1 Chpter 3 Eponentil nd Logrithmic Functions Section 3. EXPONENTIAL FUNCTIONS AND THEIR GRAPHS Eponentil Functions Eponentil functions re non-lgebric functions. The re clled trnscendentl functions. The eponentil

More information

Calculus Module C21. Areas by Integration. Copyright This publication The Northern Alberta Institute of Technology All Rights Reserved.

Calculus Module C21. Areas by Integration. Copyright This publication The Northern Alberta Institute of Technology All Rights Reserved. Clculus Module C Ares Integrtion Copright This puliction The Northern Alert Institute of Technolog 7. All Rights Reserved. LAST REVISED Mrch, 9 Introduction to Ares Integrtion Sttement of Prerequisite

More information

SECTION A STUDENT MATERIAL. Part 1. What and Why.?

SECTION A STUDENT MATERIAL. Part 1. What and Why.? SECTION A STUDENT MATERIAL Prt Wht nd Wh.? Student Mteril Prt Prolem n > 0 n > 0 Is the onverse true? Prolem If n is even then n is even. If n is even then n is even. Wht nd Wh? Eploring Pure Mths Are

More information