Introduction of Finite Elemény neh Ana, ie a6) for the simplest ¢ stiffness matrix axisym, opment of the section isa plane triangle. We' begin with the developmen a element, the triangular torus, tion 2, Axisymmetric Formulati neato Praae a typical axisymmetric triangular element wil 3a iting i Fig.4.5. Fig. 4.5. Typical axisymmetric element In two dimensional problems, the displacements and distributed body force ‘lie’: we indicated by x-y plane. But in case of axisymmetric problems, these values are a r-z plane as shown in Fig.4.5: For two dimensional problem, the displacement vector w is given by, Gey) = {"} i where, wand v are the x and y components of u respectively. In case of axisymmetric problems, the displacement vector uw is given by, > u(r,2) = {“} where, u and w are r and z components of u respectively, The stresses and strains for two dimensional element are given by, Stress, {0} = Strain, {e} = & scanned with OKEN Scannerproblems — Vector Variable Problems call PU exisymmettic element, stresses and strains are given by, of i og? Stress, {0} = Tre where, 6, — Radial stress 6, —> Longitudinal stress 69 —> Circumferential stress Tz — Shear stress Strain, {e} = Yes where, e, -> Radial strain e, —» Longitudinal strain ¢ > Circumferential strain Yq > Shear strain For two dimensional problem, body force is given by, St (pt F =-18,f AG Incase of axisymmetric problem, Fe {E ic Element .3. Derivation of Shape Function for Axisymmetri : 4 i shown in Fig.4.3., ce Cana a veri ‘element with nodes 1, 2 and 3 as in Fig. jider an axis ieee the nodal displacements be m, Wi ta "2 and ts Ws; Wi “1 wy Displacement, (#} = Vw, us : Wy ‘ - ait: & scanned with OKEN Scanner‘AN Introduction of Fini te Ele ment ince the triangular element has two degrees of freedom at each nod a co-ordinates. arte te, Displacement functions, u = aytagr+ayz w= aytasrt+agz ie (yy where, a}, ay a, ay as and a ae global of generalized co-ordinates, “ig Let my = ay tar +a32 th = ay tayh +4427 ty = ay tay ty +323 Write the above equations in matrix form, 1 ln zy a yal A el ie U3 1 ry z3 % 4 a ly zy uy Qe =| 1 ry Zz, % i a 1 73 23 Uy " ae $27, Sa se sgt hs ene Let, D lin x t+ o-t living. 2, ct pi= & 7 een (Al Find the co-factors of matrix D, mt a ny 7 23] — 223-7322) Va, 1 2,| 7 ~~) = @-2) ly, 1 | = +03-7) i m2 Cy == 5°23, 1731321) = yz) —71 25 & scanned with OKEN Scannerproblems Vector Variable Problems ond! 2 1 “| = Goa) Caz ln Ca =~ \nergll Ter = triad mi 2) : Cy = * me = 2712 Lz Cy = |y y = ~(@)-2)) = @-z,)- ln Cao * | yopyaGaren. (1.231322) @2-23) (35.9) = c= | G321-14)) @-2) God) (22-22) 1-2) Ma) My Zy—Pg Zo MeFi 23MIN) che 29-23 aieet fos gh ret AB mon We know that, {08.0} 10 ; ln 4 p-|in pel Lory 23 , ry 4 2 \in® pic i 1 13 %3 ) vee (419) py-n@sAd tO \DI= 10223773 jn equation 41D» 1 x Tag Fai 3-72) t= Try" 3 =ryzy 17277271 41) > D (2737-73 23-972 371 ne wie ; Aaa a a non e non 2 ants Lett) . & scanned with OKEN ScannerSubstitute D-! value in equation (4.16), a ln 7 (416) 3) 2p =] 1 2 a 1 ry 23 1 x - — SS (23-73 22) 1) 3-22) $2) (3-19) 1923-1322 1321-123 My 2.-1221 uy ra a ae 1 a day 13-1 "N-"s 277 us The area of the triangle can be expressed as a function of the r, z co-ordinates of the Nodes 1,2 and3. ; ory 2 1 == | bay Z A 2 2 22 Lory 23 1 : A = 511 0423-7329) =) (3-2) #2; (5-7) => [2A = 225575 22)~"1 G3 722) $21.5 —7) Substitute equation (4.21) in equation (4.20), tion (4.14), we know that, a 1223—%322 13 21—M 23-1 29572] [My 1 G20) >) 2p = Fr) 2725 a3 eit 74g |x % TT ae) RTn a Oy G2 ay uy sly 1 { 27 = 35K] Br Bo By |x 4 wm 42) % Ub % 4 where, o., = T223—T3 298 Oy = ry zporpzgpi os ajv= B) = z)-z5; By = 24-255 B = Y= 3-195 h = 1-713 = a t+a,r+a3z & scanned with OKEN Scannerpimensional Problems ~ Vector Variable Problems i two. an write this equation in matrix form, weet My we tra ye ay i Gy Og) Oy uy = [lr z]xzq| Br Bo Bs |x) & fab saclay YN Ye 13 uy My Oy Oy uy 1 = 3,4 Ur zl] Bi Bo Bs |x 4 nN bh % us ~ Ys / ih = Felon BF erie gt Byrty)z) Gy +Byreyzbe ym PO? [Note: (1x3) xB x3) = (1x3)] uy ay t+Byrtyiz Oot Bort 22 wate tne , 4 s|—~. TSA tel nar’ -[ 2A 2A 2K : ‘The above equation is in the form of : 1 w= NNN) % ca G29 ; ! z : wr)e i ao of wee, INNa Nal) foros te ven (4.25) Similarly, os ; : i ft ‘ay t+Birty2 function, Ny = 2a here, Shape funct - ii gt Bor +2 " No DA 7 yh HE 03+ By r+¥32~ “aN 2A & scanned with OKEN Scannera Introduction of Finite Element Analysis We can write equations (4.24) and (4.25) as follows: u = Ny uy +N my +N343 (426) w = N,w, +N), +N; w3 ws (421) Assembling the equations (4.26) and (4.27) in matrix form, Displacement function, uy Wy | _ ap Ni (0 Ny 0 N;. 0 uw, (428) we.z)} [| 0 N, 0 N, 0 Nj] )™{ ~ ty % u(r,z) = { 4.3.4. Strain-Displacement Matrix [B] for Axisvmmetric Trianaular Element & scanned with OKEN Scannerai Introduction of Finite Elejn, idl Tere Ana [4:16. if Matrix [K]_. im 43.6. Assemplgn fine men SHITE ‘ = “We know that, i > t took : i av = [By i srittness matrix, [KJ = i [D](B] av = By oe fav 7 tel" [D] [B] V iene (X= 27 ABD ay [eV =2nray where, iq aN " pb e nytt Co-ordinate, 7 = 3 3 |A = Areaofthe triangularelement = 4 (Bx) ? [B] = Strain-Displacement matrix, 7 a 4a ; oe By 0 Ss Tihs =hl tht 9 u n By where, oy = 723 — 73235 By Rony s2y N= yn; Co-ordinate z= = & scanned with OKEN ScannerIntroduction of Finite El ey, ie 4.18} ie. =m fot eco)” da { fw )" ey" yt * e(p) xr dA Internal virtual work = 2n A ($0810) = (4) pp ae f (TBI) xe) * dA Ee DPI A an f (e)TBI"D] x [B] x tv } xr da A [~ e(6) = [B] [y] where, y = Arbitrary nodal displacemeny an x (uJ? (BIT DIB) (y}r faa = 2nx {uw}? [BJ [D] [B] {y}r x A i Internal virtual work, W, = {}" x 2nrA [BJ" [D] [B] x fy} vs (4.46) ‘The above equation is in the form of, W. = (ux [K]x{y} (447 ; : aan, ela BL) Comparing equations (4.46) and (4.47), . oe Stiffness matrix, [K] = 2nrA [B]" [D] [B] [Zaumpte 41) The nodal co-ordinates for an axisymmetric triangular ele given below, 1, =10mm; 2,=10mm 1 =30 mm; 2)= 10mm 13=30mm; 23=40mm Evaluate [B] matrix for that element, Given: Co-ordinates: "=10mm; 2,=10mm )=30mm; z= 10mm 753=30mm; 2)=40mm To find: Strain-Displacement mat y @ scanned with OKEN Scanner1 problems ~ Vector Variable Problems [pe p ‘We know that, oss splacement matrix sui By 0 By 0 8, 0 a m2 9 M0 he ce Sept 0 +B o epee 9 i [B] = 2A 0 1 0 h 6 7 ee Nn B, Ya By % By [From equation no. (4.36)] where, A = Area of the triangular element ln 3 1 ry zp [From equation no. (4.21)] 1 ryezg Nie (ra 23-15 22)" @3-22) 421 -7D] [0 x 40) — (30 x 10)=10 (40 ~ 10) + 10 GO-30)] nie ie = ; x (600) nytt ts 10+30+30 Co-ordinates: Tinerras: 3 «aot ene. : P= 23.334mm zy tz,+z3°, 10+10+40 tnt nets 3 z= 20mm Gay =P ry 2y'— 75227 G0 x 40) —G0% 19) = 900 mm? a4 ty Pirsz 723 = BOX 10)+ (10 x 40) ——— i [oo ome oscil gy =r F2—7221 F (10x 10) — GO x 10) : Oo Gd fag = 5500 men] & scanned with OKEN ScannerIntroduction of Finite Elemen, ‘Ar 4.20) “ui B= 2-437 10-40 0 mm by = ay-81 = 40-10 8 By = 30mm {a by = 2)-# = 10-10 By = 0 ; ye yin = 30236 y= o i ay = ory = 10-30 = 20mm] pie % = rar S30-10 : PE ei] PO- Ors ODIs a uz _ 900 : = eet = hig + C30 4 eS = 8571 mm S44, e _ -100 20% 20) _ 5 +t = Bayt on +! 3394 3571 mm 3 2 _ 200 20x20 2. 7 Bat "2 = 33334 + 0 + 23.334 7°8:571.mm 5 oy PON ZO iS oz as Substitute, AB Boy By» G+ BATS Pape BE = and 7 values in equation (1), 3 30 0 30 __1_ | 8571 0 asm (= Bl= 75300! 9g 00s 9 0.0142 0 B=] 0 & scanned with OKEN Scanneranal Problems ~ Vector: Variable Problems ai 421 pid all strain-Displacement matrix, esi -0.05 0 0.05. 0 0 0 ee 0.0142 0 0.0142" 0 0.0142 0 t 0 oo 0 -0.0333 0 0.0333, 0 -0.05 -0.0333 0.05 0.0333 0 The nodal co-ordinates for an axisymmetric triangular element are givert below: 7, =20 mm; 2,=40mm 1) =40mm; 72 = 40mm 1;=30 mm; 73= 60mm Evaluate [B] matrix for that element. » Given: Co-ordinates: r= 20mm; 2, =40.mm , r= 40mm; 0 7 =30 mm; To find: Strain-Displacement matrix (By © Solution: We know that, Strain-Displacement matrix B, 0 B> s Bai | 0 @. nz op 8594 BF 1 | Sag et 0 Se ome a -il7 r B] = 24 ° ‘i oe fap 0S 2 Of] ont 8 Yagn « 623 2% Bs _ “ [rom equation no. (4:36)] where, A= ‘Area of the triangular element 1 Peapl eter eo “at =} Lory #2 | cui) > {From equation no. (4.21)] 1 et q 3-22) +41 0s- 1] 4 1 . 3 nanan L 2 [(40 x 60)- (60 = 40) +40 (20-409) & scanned with OKEN ScannerIntroduction of Finite Elemeny 4 = £ (1200-400 - 400] [A= 200mm?) rt thy — 20440430 Co-ordinates: rer 3 [£= 30mm, zytzp +z) _ 40+40+60 3 3 2 = = 46.667 mm Oy = 723-1322 = (40 x 60)—G0.x 40) \ @ = 1200 mm? Gy = 752,723 = (30% 40)-(20x 60)” a = 0 = 729-721 = (20% 40) 40 x 40)... 03 = ~ 800mm? are B, = 22-75 = 40-60 inemaosiquiCenig By = —20mm apis By = 23-z, = 60-40 La ee By = 20mm 4 | By = 2)-z) = 40-40 By = 0 YN = Y3-7 = 30-40 ? yy = -10mm Y2 = "—r3 = 20-30 5 t2-7 = 40-20 a & scanned with OKEN Scanner. sional ‘problems ~ Vector Variable Problems iy 7 = 2 ne e hs Bt p= Gq FE 20) +Cx 6.67 a ° 30 Sap eRe = 4444 mm 7 122 0 2 +Byt 7 = 39t -10 x 46.667, }? Br = 39720 + s Se | 8 Ys? _ -800 an | 2 +Byt7 = 39 FOF es sn | a, te A, By, Ba» Bs. G+ z 2 Yaz Oy substitute A, By Ba & 7 Bt GH tRt Go > tBt 1 snaandys values in equation (1), 20.0 2 0 0 0 1 | 4444 0 4.444 0 4.444 0 | as (Bl = 3x200| ‘o'-10 0 -10 0 20 10-2010 °26"""20"" 0 | 0.05. 0 | 8} 0:05: za Ole rmelT “HOr Ooi 0 ool 00.0111 0 | BB fg 9.038270 a.025 “0” 00 9.095 “Uo.os “6.025 © 0.05 005 0 Result: Strain-Displacement matrix MomnA = pos 0-5 gyons, eh. re oul o. , ool o 0.011 0 Bi- |} 9 | Loos, jos 0025 9 0.05 19.025 0.05 (0.05 0 | 9.025 | “0.05 the stiffness matrix. Take For the element shown in Fig.» E=200GPaand vane determine | h b b E p & scanned with OKEN ScannerIntroduction of Finite Elemeny ~ [424 The co-ordinates shown in Fig.() are in millimeters. Given: Co-ordinates: r=O0mm; 7, =0mm 7 =50mm; 7) = 0mm ry =50mm; 242 50mm Young’s modulus, E = 200 GPa ; = 200 x 109 Pa = 200 x 109 Nim? (0,0) (60, 0 ' (4.24) (t,22) = 200 x 103 N/mm? = |——__. E =.2x 108 N/mm? Poisson’s ratio, v = 0.25 To find: Element stiffness matrix [K]. © © Solution: For axisymmetric triangular element, stiffness matrix [K ] is given by, (K] = 2r AB] [D) Bl (1) age [From equ: where, A = Area of the triangular element =~ Try zy Saat tt 2 | : 1 rs 23 or We can calculate as, ; 1 ae Area ofthe triangular element, A= + x Breadth x Height = $ x 50x50 A = 1250 mm2 = Q Co-ordinate, r = “tain O+50+50 ee ? = 33.333 mm Co-ordinate, 7 = 21422423 040450 ! 3 ee 2 = 16.666 mm & scanned with OKEN Scannerote we 0 a tcebattig( (1S) slau : pip matrix, [D] = petationship ™ Dl=Gayd-2y| YoY inv 0 0 1-2 _ y [From equation no. (4.42)] 1-025 025 025 0 Sic! 025 1-025 025 0 = D> (4095) 1-@x 025) 025 © 025.1-025 0 12025 0 2-02 38 » ° 2 075 025 025 0 J, . 3110 : , | 025.0275 025 0 ee oats 310 = 320x10| 95 9.25 0750 | peclDh xO 13 or 007 035 Een 0 0001 3110 ix | 1 Oo | pps sox] 13.0 @) 00 01- Weknow that, aol Strain-Displacemtent matrix: geet . é 8, Bp 0 By ; : ay te ; oa Seg oe 0 Sth HE oP tty se usifeaps ir ® a 0 " e y ‘ ; { Ong a i fh Dae | +: [From equation no. (436)] = (30% 50)— (60 « 0) = "3500 mm? oHnen se (50x 0)- (0 50) | = 00)— (0x 0) Ye Wea & scanned with OKEN ScannerIntroduction of Finite Remon, 426 ea [e= 9) By = F-#3 = 0-50 By = 50mm Pa @ ag-ay = 50-0 By = 50 mm B= -H = 0-0 Bs = 0 th = 1ynr = 50-50 f 1 =. 0 et = 0-50 yi=ih-n = 50-0... p | Ol x O88 % = 50mm - = Sep4t- Faayyt C50) +0 = 25 mm pa + 22. o4s0 + £502 16.666) 25 mm : 33.333 itn s S +p+ 2 = o+04 50 « 16.666) = 25mm» 33.333 Yaz Oy a. ZO. Ys? Substitute B,, By, B3, > + B+ ne t Bet 3p, t Bat tp prerms Substitute Area A, value Values in equation no. (6) r 50 0 50 000 25 0 25 025 0 0 0 0-50 0 50 wont 0 -30 -50 50 50 0 500 50 0.0 0 1 25 0 25 025 0 2*1250]} 0 09 | 9 =50° 0 50 0 -50 -50 50 50 0 (B] = on ae & scanned with OKEN Scannery onal Problems = Vector Variable Problems im We 4.27) ae es, 6 450 0, 0 05] 0 1 2 0 25 = 103 pata 025 0 oe 80x10} 1 13.0 ]%2x1250) 0 0 0-50 0 50 0001 0-50 -50 50 50 04 a4 0 40 2000 Tks 1 0 25 101 010 a 3 80x10?] 1 1 3.9,,% 2x1250| 0 0 0-202 0001 0-2-2 220 3110 99 250100 13 toi} ¢1.20 81h olr 0 = 800], 130 || 0.0 0-202 0001 022 22.0 _gei40r0 O+04040 G+1404000040-240 0414040 0+0+240 Tpese0e0 404040 2434040 0+02240 (0734010 0+0+2+0 = 800 | 5414040 OFOHOHO 2414040 4026400414040 0+0+6+0 qeroe0 OF0H0-2 0+040-2-OFOHOH? ovoroi2 o}0H0#0 005-232 ; @ pie) = 800] 0 We know that, ° & 3 ° g 3 1 BI = 21250 29; 0. 25 10 Bie es) g 0 0-24) + Q) & scanned with OKEN ScannerS| Introduction of Finite Element sn sig 4.28 21 0 0 00 0-2 5 07-212 “i 10 5-3 » 21.02 1. an 32 => BI DIB) = O01}. 9 9 -2 2 Mylo 3 6-16 o1 0 2 0.22 229 00 2.0 21 0 0 f 00 02/75 0 7212 2160 2ffol 05232 = 001x800] 5.9 22 |) at) 0 3-616 0.19052 [L)02 22 220 060 2 90. 1OHI+0+0 0+0+0+0 \~14+5+0+0'4=2+040 -243+0+0 4424040" O+04040" 001014 “OF0F0HE 0040-4 0400-4. O+0+0+0 VOHI+O+O 0F0+0H4 | 14454044) -4-240-4-124340-4 4424040 = O+04240 OF0H0-4 0+0-6-40 04041244 100-244 040-1240 OFHO+0 OF0H0-4 —OFS+H0-4 O-2H0H4,. 0434044 0424040 040-240 0404040 — OF046HO 040-1240, 0404240 O+0H12H H 0-9 [28 1 aye oe hae 04 4 44 04 & J 9 4 2B-10 1 6 BFPIB)= 8) 3/4 to “16 {3 a |? li al a2207 92 | 2.0 6 41292 4 Substitute [B]" [D] [B] value in equation (1), 1 O62 a 0 44-4 06 9 23 <1 1 *é ad [K] = 2nrAx8 & scanned with OKEN Scanner5 pimensional Problems — Vector Variable Problems 5 4.29 Mi 9s sgPitgleser?, 94 44-4 0 = 2x x33.333x1250x8| 2 4 3-0 1 6 24-10 16 2 -12 14 10°97 2 2 0 6-12 2 12 Hoo 9 21 04 4 4-4 0 a ie 9 4 23-10 1 -eeness matrix [K] = 2.094 x 106 ig ai A sine 2-4-1016 2-12 | Nant peng 1 (AplCa7 Bigh | 7S Loh ym PMID? 212 tmay be noted that stiffness matrix [K.] is symmetric. graisd a0 om Ho 2 1 - 4 44450 9423.19 1 6 | am Result: stiffness matrix, [K] = 2.094% 106] 9 4 19 16 2-2 —Sssr[e abe V2 ah 2 0 6-12 2°12 in Fig.(), determine the ve For the axisymmetric elements shown i i = nani aid v = 0.25. The co-ordinates shown in Fig.) stiffness matrix. Let E 2.1 x 105 Nim anes¥ en ae are in millimeters. & scanned with OKEN Scanner(430 Introduction of Finite Element Angy abysig Given: Co-ordinates: ry = 0 mms Omm; rs Young’s modulus, E Poisson's ratio, v = 0.25 To find: Stiffness matrix, [K]}- © Solution: For axisymmetric triangular ele- ment, stiffness matrix [K ] is given by, [K] = 2xr AB D]B) [From equation no. (4.43)] Area of the triangular element, A 0mm; z= 0mm 2.1 x 108 N/mm? y= 50mm (0,0) (60, 0) (2) (2.2): (I) Fig. (ii) 1 } x Breadth x Height =" >-* 50 %50 Co-ordinate, r = We know that, Stress-Strain relationship matrix, [D] 2.1 x 105 = (D] ~ (F025) (1-@x 023) ‘A =. 1250 mm? : _intntn _o+s0+0 3 = 3 oe oer 0) cw ct tda* es _ OF 0450 adits Senna P= 16667 mm | t-yew vee ea Ewe Tovey 2 (l#v)(- vod sy, , 1=2¥ / oo ee [From equation no. (442) 1-0.25...0.25' 0.25 0 0.25 1-025 0.25 025 0.25 1-0.25 0 0 0 & scanned with OKEN Scannervines Problems = Vector Variable Problems fie 075 025 o25 = 336x103 | 975 975 025 9 025 028 0.75 [D] = 84x 103 We know that, Strain-Displacement matrix a. ae: hd B= alr Pr ak fare aE 2) orpgtfoaps #1 [From equation no. (4.36)) where, ay = mz30 87: ay = 2500? : oy MAGNET a = 0 j é e 8 = (0 x 0) +50 * 0) 0, = 22 & scanned with OKEN Scanner= Inredtion of Fite lene as By = 0 My = yom = 0-50 ‘ 1 = 50mm th = ners = 0-04 =o = on = 90-0 a oy nz | -50 x 16.667) _ ‘ & Se p+tt + (50) +e = 50mm Oy nz _ _0 0 x 16.667 _ > SPB ETT FO 16.66T Oa ty wt 0 5, 50x 16.667 _ > 7p tBt > = 76.667 + 9+” 16.667 = Om Siibsitute By By Py T+ B+ =, +p, BF, Fp + BE roandn values in equation 'no. (6), 7 =} 50 050 0 0 0 yf 50" 0 5005 0.50 0 fe BI= 34) 9 -50 (0:70; y0°50" |; -50 50 0 50°50 0° Substitute Area, A value, . = 50 050 0 0 0 1 50. 0°50. 0°50° 0 £3) = [Bl = 751250] 0 -so 0 0 0 50] ry 50 -50 0 50 50 0 -1 1 = 50. 1 1 2x 1250 0-10 -1 -10 0 0 0 0 0 1 Hoo o! 0 co 1 0. & scanned with OKEN Scannerlems = Vector Variable Prob sional PP OO erro -1 01000 1 o101d pc 002} 9 10001 wD) a -1 0110 : 3110 +1 01000 ; 1310 lo1010 2 4x 10 x 0.02 | ppm : Imre 7 oeoo0r E ooo 1, ee a) 31107F-1 01000 1310]| 1'01010 = 1680) 7130}! 010001 ooo wba Horio 3414040 OHO-140 3414040 O+F04040 O+1+0+0 0+0+140 1434040 OH0-140 143+0+0 OF0HOHO 0434040 +0440 = 1680 | 1414040 0+0-3+0 1414040 OFOHO+0 OF1+0+0 0404340 O+0+0-1 \0+0+0-1 OFOHOHO O+0+0H 0+0+0+1 0+04040 214011 2-14031]" 8) _ piel = 1680) 9 -3° 2:04 3 -1-10110 Weknow that, £ 8 i 1 01000 ‘ 101010 B= 00) 9 ;ooor -1-10110 -11 0-1 *) 90-1 -1 0 11 0. Og] cast wo [BI = 002]: 90 0 1 o1 01 00 1 0 | <1 0-1 2-14011 a9 ott 2-14031 ee 80 [BY () (B] = 0.02] 9 00 1 xd@ wee Qadye OF er or oo 10 & scanned with OKEN Scannera Introduction of Finite Element dng sis 434 -11 0-1 = DO rh 4! 2-14014 11.0 0 214034 = 002x160] 99 0 1/*%| 032013 o1 01 -1-101 00 10 a DIOHL LeLHOHL AH44040 OFOHO-1 -14340-1 ~1414049 OL0H0HT OFOFSH 040-240 O+0+0-1 O+0-I-1 040-345 “3424040 -I-14040 4444040 OFOHOHO 1434040 1414049 5336 0+0+0-1 0+0+0-1 0+0+0+0 O+0+0+1 OF0+O+1 0404049 OF240-1 O-140-1 0+4+0+0, 0+0+0+1 0+3+0+1 OF1+049 0+0+040 040-340 O+0+2+0 O+0+040 O+F0FI+0 0404340, 5.1 0-1/1 0 142-123 0-2 8 0 4 2, = BPDIBI= 36), 9 1 1 o w+ (10) 12.4 1.4 1 03 20 1 3 E Substitute [B]" [D] [B] value in equation (1), 5 1 0 2 1 4 2 2nr Ax 33.6 1 8 = [K] 0 12 4 0-3 2 0 2x m x 16.667 x 1250 x 33.6 & scanned with OKEN Scannerjonal Problems = Vector Variable Prog pines lems 4.35 $1 ofp alo tg 04 [K] = 439x196] 92 8 1 0 Nimm 12 4 032 woeon LA 14 O14 may be noted that stiffness matrix [ K ]is symmetric 5 10410 E 42A2a24 -gesuit: Stiffness matrix, [K] = 439x196] © 2 8 0 4 2 : Va wae 124144 03201 3. For the axisymmetric elements shown in Fig.(i), determine the element “stresses. Let E = 210 GPa and v= 0.25, The co-ordinates (in millimeters) are shown in Nimm ‘The nodal displacements are: 4, = 0.05 mm; w, = 0.03 mm 4, = 0.02 mm; w, = 0.02 mm a a =Omm; w;=Omm [Anna University, Dec. 2006, M.E. CAD/CAM] & scanned with OKEN ScannerIntroduction of Finite Elem, ent Analy, i 4.36, Young's modulus, E = (60, 0) (ra 22) Fig. (i) 210 GPa = 210 x 10° Pa = 210x 10° Nim? = 210 x 103 Nimm? E = 21 x 105N/mm? Poisson’s ratio, v = 0.25 Nodal displacements are* “1 =0.05 mm; w, = 0.03 mm t= 0.02 mm; Ww = 0.02 mm i= 0mm; w;=0mm To find: Element stresses (i) Radial stress, 6, (ii) Circumferential stress, 9 (iti) Longitudinal stress, o. (iv) Shear stress, 1,. © Solution: [Note: In this problem, the given axisymmetric triangular element, co ordinates, Young’s modulus and Poisson’s ratio are same as Example 4.4, So, we can take [D] [B] value from Example 4.4.] We know that, 2-1 " [DI[B] = 1680 i Stress (6 } = (D] [B] {wv} | Uy » ay | = [D) [B] & scanned with OKEN Scanner= 1680 = 1680 = 1680 => Radial stress, 9, Circumferential stress, 9 Longitudinal stress, . Shear stress, Tz oO, = <, [Example 46 ] For the axisymmetric Take E = 2.1 x 105 Nima and ¥= 038 4.37 2 = 14044 0.05 22-1403 4] | 2% A0110 0 0 (2x 0.05) +e 1 x 0.03) +4 : 0.02) +0+0+ (2 0.05) + (1 <00)+4x0M +0060 0+ (3 x 0.03) +(2x 0.02) +0+0+0 = 0.05 + 0.03) +0+0.02+0+0 (Note: (4 x 6) x (6x 1) = (4 -0.05 a 0.15 -0.05 -0.06 -84 1 1 -100.8 = -84N/mm? 252 N/mm? -84 Nimm? = -100,8 N/mm? -84 N/mm? 252 Nimm? 34 Nimmo sp et (gdetermint the element 100.8 N/mm? clement shown i FIS- The co-ordinates shown in Fig. (0.076 © a 4, = 0.05 mm; 7 00 | y= 0.02 mm; mcnn, = O mm; wy ne & scanned with OKEN Scanner