Comparison of plane-stress, generalized-plane-strain and. and 3D FEM elasticplastic analyses of thick-walled cylinders subjected.
|
|
- Juliana Flora Griffith
- 5 years ago
- Views:
Transcription
1 Loughborough University Institutionl Repository Comprison of plne-stress, generlized-plne-strin nd 3D FM elsticplstic nlyses of thick-wlled cylinders subjected to rdil therml grdient This item ws submitted to Loughborough University's Institutionl Repository by the/n uthor. Cittion: KAMAL, S.M....et l., 07. Comprison of plne-stress, generlizedplne-strin nd 3D FM elsticplstic nlyses of thick-wlled cylinders subjected to rdil therml grdient. Interntionl Journl of Mechnicl Sciences, 3-3, pp Additionl Informtion: This is n open ccess rticle under the CC BY license. Metdt Record: Version: Published Publisher: c 07 The Authors. Published by lsevier Ltd. Rights: This work is mde vilble ccording to the conditions of the Cretive Commons Attribution 4.0 Interntionl CC BY 4.0) licence. Full detils of this licence re vilble t: Plese cite the published version.
2 Interntionl Journl of Mechnicl Sciences ) Contents lists vilble t ScienceDirect Interntionl Journl of Mechnicl Sciences journl homepge: Comprison of plne-stress, generlized-plne-strin nd 3D FM elstic plstic nlyses of thick-wlled cylinders subjected to rdil therml grdient S.M. Kml, U.S. Dixit b,, A. Roy c, Q. Liu c, Vdim V. Silberschmidt c Deprtment of Mechnicl ngineering, Tezpur University, Npm, Tezpur , Indi b Deprtment of Mechnicl ngineering, Indin Institute of Technology Guwhti, Guwhti , Indi c Wolfson School of Mechnicl, lectricl nd Mnufcturing ngineering, Loughborough University, Loughborough L 3TU, UK r t i c l e i n f o b s t r c t Keywords: Thick-wlled cylinder lstic plstic Plne stress Generlized plne strin Three-dimensionl finite element method Therml stress In mny industril pplictions, thick-wlled cylindricl components re subjected to high pressure nd/or temperture. During the opertion the cylinder wll my undergo elstic plstic deformtion. This pper presents plne-stress nd plne-strin thermo-elstic plstic stress nlyses of thick-wlled cylinders subjected to rdil therml grdient. A three-dimensionl finite element method 3D FM) nlysis of the thermo-elstic plstic stresses in thick-wlled cylinder is lso crried out. The 3D FM results re compred with the nlyticl plne stress nd the generlized plne strin nlyses in order to study the vlidity of these models on the bsis of length to wll-thickness rtio of cylinders. The plne stress nd generlized plne strin nlyses re bsed on the Tresc yield criterion nd ssocited flow rule. The strin hrdening behvior of the mteril of the cylinder is tken into ccount. It is observed tht for the length to wll thickness rtio of more thn 6, the generlized plne strin nlysis cn provide sufficiently ccurte results. Similrly, for the length to wll thickness rtio of less thn 0.5, plne stress nlysis cn be used. When the length to wll thickness rtio is more thn 0.5 but less thn 6, three-dimensionl nlysis is needed. 07 The Authors. Published by lsevier Ltd. This is n open ccess rticle under the CC BY license. ). Introduction The thick-wlled cylinders subjected to pressures nd tempertures find severl pplictions, e.g., in chemicl industries nd nucler power plnts. In most of the cses, the design ttempts to keep the stresses in the cylinders within elstic limits. However, it is lwys better to crry out n elstic plstic nlysis in order to get n ide bout sfety in cse of untowrd situtions. Moreover, in n utofrettge process, the plstic deformtion is delibertely produced to induce compressive stresses in the inner side of the cylinder. There re lso exmples where plstic deformtion is permitted by design. Thus, the elstic plstic nlysis of the thick-wlled cylinders is n ttrctive reserch re. lstic plstic deformtion of thick-wlled cylinders due to internl pressure loding is well recognized nd hs been investigted by mny reserchers. Mny nlyticl solutions re vilble for stress, strins nd displcement fields during the elstic plstic deformtion of cylinders. An erly nlysis of elstic plstic deformtion in thick tubes ws crried out by Hill et l. [] ssuming non-hrdening mteril, Tresc yield criterion nd plne strin condition. Go [] developed closed form nlyticl solution for the elstic plstic stress, strin nd displcements of n internlly pressurized thick-wlled cylinders under plne stress condition using Hencky s deformtion theory nd von Mises yield criterion. Durbn [3], Go [4], nd Bonn nd Hupt [5] investigted lrge elstic plstic deformtions of thick-wlled cylindricl tube under internl pressure. An nlyticl solution for stress, strin nd displcements in thick cylinder subjected to internl pressure ws developed by Go [6] using strin grdient plsticity theory. An nlyticl solution for elstic plstic stresses considering the Buschinger effect nd Tresc yield criterion in thick-wlled cylindricl vessel mde of elstic liner-hrdening mteril ws given by Drijni et l. [7]. The clssicl solution for thermo-elstic stresses in thick-wlled cylindricl bodies under stedy stte temperture distribution due to temperture grdient is well known. However, only few ppers tret the thermo-elstic plstic nlysis of thick-wlled cylinders nlyticlly. Blnd [8] crried out n elstic plstic nlysis of thick-wlled tubes of work hrdening mteril subjected to internl nd externl pressures with outer nd inner surfces mintined t different tempertures. He considered tht yielding took plce only on the inner side of Corresponding uthor. -mil ddresses: udy@iitg.c.in, usd008@yhoo.com U.S. Dixit). Received 4 Septembe06; Received in revised form June 07; Accepted 0 July 07 Avilble online 0 August / 07 The Authors. Published by lsevier Ltd. This is n open ccess rticle under the CC BY license. )
3 Nomenclture inner rdius of cylinder b outer rdius of cylinder c, d, e, f interfce rdii Young s modulus of elsticity K hrdening coefficient n strin hrdening exponent r rdius of cylinder T temperture t the inner wll of cylinder T b temperture t the outer wll of cylinder u rdil displcement α coefficient of therml expnsion ε r totl rdil strin ε θ totl hoop strin ε 0 constnt xil strin ε e r, ε e, ε e θ z elstic rdil, hoop nd xil strin r, p, ε θ z plstic rdil, hoop nd xil strin uivlent plstic strin ν Poisson s rtio σ r rdil stress σ θ hoop stress σ z xil stress σ uivlent stress σ Y yield stress the cylinder leding to the formtion of n inner plstic zone. Wong nd Simionescu [9] developed n elstic plstic nlyticl model of thickwlled tube subjected to internl heting nd pressure ssuming smll displcements, plne strin condition nd the yield criterion of Tresc without strin hrdening. They did not consider the cse in which there re outer nd inner plstic zones with n intermedite elstic zone. The thermo-elstic plstic deformtion of tubes with inner plstic nd outer elstic zones due to internl het genertion ws investigted by Orçn nsln [0] considering the temperture-dependent mechnicl nd therml properties. Sdeghin nd Toussi [] crried out n xisymmetric thermo-elstic plstic stress nlysis in cylindricl vessels mde of functionlly grded mteril. The stress nlysis of functionlly grded thick-wlled hyperelstic sphericl vessel subjected to internl nd externl pressure ws crried out by Anni nd Rhimi [], lbeit without considering therml stresses. Some ppers in the literture del with the thermo-elstic plstic nlysis of solid cylinders nd sphericl vessels. Orçn [3] crried out thermo-elstic plstic nlysis of elstic-perfectly plstic cylindricl rod with uniform internl het genertion for generlized plne strin condition bsed on Tresc yield criterion. He considered the cse of two plstic zones seprted by n elstic zone. Cowper [4], Johnson nd Mellor [5], nd Drijni et l. [6] nlyzed the elstic plstic stresses in thick-wlled hollow sphere under stedy stte rdil temperture grdient. In recent pper, Zre nd Drijni [7] studied elstic plstic stresses in thick-wlled cylinder rotting t very lrge ngulr velocity. In generl, thick-wlled cylinder my undergo simultneous plstic deformtions emnting from both inner nd outer wlls under the condition of high therml grdient. This cse ws nlyzed by Kml nd Dixit [8,9] for hollow disk nd cylinder ssuming plne stress nd generlized plne strin conditions, respectively. The nlyses re bsed on Tresc yield criterion nd its ssocited flow rule. The effect of strin hrdening ws considered in the plne stress model [8]. However, the generlized plne strin model [9] did not include the strin hrdening, lthough the theoreticl results mtched well with the experiments [0]. In this work, the generlized plne strin model of Kml nd Dixit [9] is extended to incorporte the effect of strin-hrdening. The strinhrdening during plstic deformtion is ssumed to follow Ludwik s hrdening lw. The plne stress nd generlized plne strin models cnnot predict the dimension-rnge in which the respective ssumptions re justified. Hence, three-dimensionl finite element method 3D FM) nlysis is crried out in order to define the pplicbility of the nlyticl models. A criterion bsed on the length to wll thickness rtio of the cylinders is developed to ssess the vlidity of the plne stress nd generlized plne strin ssumptions. During the nlysis, it is ssumed tht the mximum temperture induced in the cylinder is well below the recrystlliztion temperture of the mteril. Hence, it is pproprite to tret the mechnicl nd therml properties of the mteril s temperture-independent.. Problem definition A thick-wlled cylinder with inner rdius nd outer rdius b is considered. The inner wll is subjected to temperture T nd the outer wll is subjected to temperture T b such tht T b > T. The thermoelstic plstic stress nlyses of short cylinder hollow circulr disk) using the ssumption of plne stress σ z = 0) incorporting strin hrdening ws crried out in [8]. Another nlyticl model bsed on generlized plne-strin condition z = constnt) pplicble for long cylinders is vilble in [9]. However, the generlized plne-strin model did not include the effect of strin hrdening. In the present work, the thermo-elstic plstic stress nlysis in the cylinder under the generlized plne-strin ssumption is presented incorporting the strin hrdening in Section 3. The mteril is ssumed to follow the Ludwik s hrdening lw given by []: σ = σ Y + K ) n, ) where σ Y is the yield stress in unixil tension or compression, σ >σ Y ) is the uivlent stress, is the uivlent plstic strin, K is the hrdening coefficient nd n is the strin hrdening exponent. A stedy-stte temperture distribution in the cylinder under rdil temperture difference T b T ) is given by []: T = T + ln r T b T. ) ln b To test the vlidity of plne stress nd generlized plne strin ssumptions, 3D FM nlysis is crried out for different length to wll thickness rtios of cylinder. 3. Generlized-plne-strin nlyticl model considering strin hrdening Under the condition of generlized plne strin condition z = ε 0 = constnt), the thermo-elstic stresses in the cylinder re given by the utions provided in [3], when the cylinder is subjected to sufficiently low temperture. At the inner rdius the yielding of the mteril of the cylinder begins ccording to the Tresc criterion given by σ θ σ r = k σ, σ z σ r = k σ, 3) where σ is given by the Ludwik s hrdening lw q. )). Beyond the inner rdil position, the cylinder yields s per the Tresc yield criterion: σ z σ r = k σ. 4) By substituting the thermo-elstic stresses in q. 3) t r =, nd tking σ = σ Y, the temperture difference ruired for the initition of yielding t the inner rdius cn be obtined. When the mteril of the cylinder yields s per the Tresc yield criterion, during first stge of elstic plstic deformtion the wll of the cylinder consists of three zones plstic zone I r c ), plstic zone II c r d ), nd elstic zone d r b ). During the second stge of elstoplstic deformtion, the cylinder wll gets divided into five zones two 745
4 [ r = ν r k σ Y ln r + C 3 + k K k K r r ε 3 k K r ν 4 k σ Y + ν) α T b T { r + ln ln b { k K ) n + r } ] r d αt + 3 } ε 0 C 4 r. 8) Fig.. lstic nd plstic zones in cylindricl segment during elstic plstic deformtion generlized-plne-strin). inner plstic zones: plstic zone I r c ) nd plstic zone II c r d ), n intermedite elstic zone: d r e nd two outer plstic zones: plstic zone III f r b ) nd plstic zone IV e r f ). The two consecutive plstic zones in the cylinder correspond to two different sides of Tresc yield locus. The generl geometry of the cylinder during elstic plstic deformtion is shown schemticlly in Fig.. In the following subsections the stresses nd plstic strin fields in the first nd second stge of elstic plstic deformtion re obtined by incorporting strin hrdening using Ludwik s hrdening lw. 3.. First stge of elstic plstic deformtion The stresses nd plstic strins in ll the three zones during the first stge of elstic plstic deformtion re obtined in similr mnner s described in Kml nd Dixit [9]. The stresses in the elstic zone, d r b re provided in [9]. However, the stresses nd strins in the plstic zone get modified due to the incorportion of strin-hrdening following the Ludwik s hrdening lw. The resulting fields of stresses nd strins in the plstic zones I, r c nd II, c r e obtined s follows: Plstic zone I, r c : Using q. 3) in the stress uilibrium ution [9], the rdil, hoop nd xil stresses in the plstic zone I re obtined s σ r = k σ Y ln r + k K r + C 3, 5) σ θ = σ z = k σ Y + ln r ) + k K ) r n ε + k K r + C 3. 6) Using the condition of plstic incompressibility, the strindisplcement reltion nd the generlized Hooke s lw, the expression for the rdil displcement, the totl hoop ndil strins cn be obtined [9]. The plstic prts of the hoop, rdil nd xil strins re obtined by subtrcting the elstic prts from the totl strin components. They re given by θ = ν [ k k σ Y ln r + C 3 + K r k K r + 3 k { K r 3 ν 4 k σ Y ν k K ) n + αt α T b T { r + ln ln b r ε ) n } ] r d p 3 } ε 0 + C 4 r. 7) { z = ε 0 ν } k σ Y ln r + C 3 + k K r ν k σ Y ν k K ) n αt α ln r T b T. 9) ln b The uivlent plstic strin is given by = p p ε 3 ij ε ij, 0) where ij denotes the component of the plstic prt of the strin tensor. Plstic zone II, c r d : Using the Tresc yield criterion q. 4)) nd its ssocited flow rule long with the strin comptibility [9], the elstic plstic rdil nd hoop stresses re obtined s σ r = C 5 r + ν) + C 6 r ν) + k σ Y ν ) k K r + { ν) ν + ν } ν) r + ν) c k K r { ν) ν ν } ν) r ν) c + αt ν ) + α T b T ln r ν ) ln b α T b T α T b T ν ) ln b ν ) ln ν) r ) n ν) ε ε 0 b ν ), ) σ θ = C 5 ν) ν) C 6 ν) r ν) + k σ Y ν ) { k K ν) ν) ν + ν } ν) + ν) r ν) r ) n c k K ν) r { ν) ν ν } ν) + ν) r ν) ε + νk K ) n c + αt ν ) + α T b T ln r α T b T ε 0 ν ) ln b ν ) ln b ν ). ) 746
5 The expression for σ z cn be obtined from the Tresc criterion given by q. 4). Using the Tresc ssocited flow rule, the plstic prt of the hoop strin in the plstic zone is obtined s zero implying r = z. In this zone, the plstic strins re given by z = C { 5 + ν) ν ν } ν) + C { 6 ν) ν + ν } ν) }{ ν + ν } ν) p r = ε k K r + { ν) ν ν ν) + ν) r ν) r ) n c k K r { ν) ν + ν }{ ν) ν ν } ν) ν) r ν) ε c + ν ) k K α T b T + ν) + ln b ν ), 3) The constnts C 5 nd C 6 re determined by using the continuity of rdil stresses nd plstic strins t the elstic plstic interfce. The continuity of rdil stresses t the interfce of the plstic zones I nd II provides the constnt C 3 nd the continuity of plstic hoop strin provides the constnt C 4. The constnt xil strin ε 0 is obtined by using the free-end condition, i.e., mking the resultnt of the xil stresses zero. To evlute the unknown boundry rdii, c nd d, numericl procedure needs to be used. The procedure involves n itertive pproch to estimte the vlues of c nd d. The initil estimtes for c nd d cn be obtined from non-hrdening cse K = 0) by using the boundry conditions of vnishing rdil stress t the inner rdius nd plstic zone II ) plstic zone II ) σ θ = σz t r = c nd solving them using FSOLV function in MATLAB. The initil guess vlue of is tken s zero everywhere in the plstic zones I nd II. With these vlues of c nd d, the vlues of p, ε θ r nd re updted in qs. 7) 9) for plstic zone I. 0 q. 0) provides the updted vlues of t ny rdil position in the plstic zone I. Similrly, using q. 3), the vlues of r t different rdil positions in the plstic zone II re updted. The updted vlues of t different rdil positions in the plstic zone II re obtined from q. 0). These vlues of uivlent plstic strin re used to obtin the updted components of plstic strin. From these components, the vlues of re updted further for fixed c nd d in both the plstic zones. This procedure is repeted till the convergence in is chieved. The integrl terms involved in the expressions cn be evluted numericlly by using two-guss-point formul. Now, using the converged vlues of, the boundry conditions of vnishing rdil stress t the inner rdius plstic zone II ) plstic zone II ) nd σ θ = σz t r = c re solved gin to get the new estimtes of c nd d. If these new estimtes of c nd e sme s the previously estimted vlues of c nd d the procedure is stopped, otherwise the whole procedure is repeted till the convergence for c nd d is chieved. 3.. Second stge of elstic plstic deformtion In the second stge of elstic plstic deformtion, the outer plstic zone, i.e., plstic zone III f r b ), develops in the wll of the cylinder ccording to the Tresc yield criterion given by σ θ σ r = k σ, σ z σ r = k σ 4) Another plstic zone, i.e., plstic zone IV e r f ) develops simultneously long with the plstic zone III s per the Tresc yield criterion: σ θ σ r = k σ. 5) The stress nd strin expressions for the plstic zones I nd II during the second stge of elstic plstic deformtion re sme s in the cse of first stge of elstic plstic deformtion. However, the constnts C 3, C 4, C 5 nd C 6 chnge due to chnge of continuity conditions. The stress solutions in the intermedite elstic zone re vilble in Ref. [9]. The utions for stresses nd plstic strins incorporting the Ludwik s hrdening lw in the plstic zones III nd IV re obtined s follows: Plstic zone III, f r b : Using the Tresc yield criterion q. 4)) nd the stress uilibrium ution, the rdil, hoop nd xil stresses re given by σ r = k σ Y ln r k K f r + C 7, 6) σ θ = σ z = k σ Y + ln r ) k K ) r n ε k K f r + C 7. 7) The plstic prts of the hoop, rdil nd xil strins re obtined in similr wy s described in Section 3. ne given by θ = ν [ C 7 k σ Y ln r ) k K r r ε f + k K f r 3 k { K r f f + 3 ν 4 k σ Y + ν k K ) n + αt α T b T { r + ln ln b [ r = ν C 7 k σ Y ln r ) r k K f + k K } ] r d 3 } ε 0 + C 8 r, 8) r r ε + 3 k K r f f 7 6 ν 4 k σ Y ν) α T b T { r + ln ln b { k K ) n + r f } ] r d αt + 3 } ε 0 C 8 r, 9) { z = ε 0 ν k σ Y ln r + C 7 k K f r ) + ν k K } + ν k σ Y ) n αt α ln r T b T. 0) ln b Plstic zone IV, e r f : Using q. 5) in stress uilibrium ution, the rdil stress is given by α σ r = ν ) T + α [ { T b T d ln + ν) ln b ν ) e } ln d { + + e ) + } r ν d ln k e σ Y + ε r 0 ν) k K e ] e p ) n r. ) Knowing the expression for σ r, the expression for σ θ cn be obtined from q. 5). The Tresc ssocited flow rule indictes tht the xil strin in the plstic zone IV is wholly elstic. Hence, the xil stress 747
6 Tble Mteril properties of luminum. Yield stress, σ Y MP) Modulus of elsticity, GP) Poisson s rtio, ν Coefficient of therml expnsion, α / C) Mss density, ρ kg/m 3 ) Therml conductivity, k W/m. K) Specific het, c p J/kg. K) component cn be obtined by using the generlized Hooke s lw [9]. The resulting expression for xil stress is given by [ { α σ z = ν ) T + α T b T ν d ln + ν) ln b ν ) e } d ] r r + ν ln ln e { ν + + e ) } r ν ln k ν d e σ Y + ε 0 ν) νk K e r νk K ) n. ) θ = The plstic prt of the hoop strin, p = ε θ α T b T ν) ln + 4 ν k K r ν ) + ν ) k σ Y b { e + ν r r k K e ν k K e r e } r d r ) is given by r ε + ν ) k K r + C 9 r. 3) The constnts C 5 nd C 6 re obtined by employing the continuity condition of rdil stresses ndil plstic strins t the elstic plstic interfce rdius d. The constnt C 3 is obtined by using continuity of rdil stresses t r = c. The continuity condition of the plstic hoop strins t r = c provides the constnt C 4. To obtin the constnt C 7, continuity of the rdil stresses t r = f cn be employed. The constnts C 9 nd C 8 cn be obtined by using the boundry conditions of continuity of the plstic prt of the hoop strins t the interfce rdii e nd f, respectively. The constnt xil strin ε 0 in the second stge of elstic plstic deformtion cn be obtined by using the free-end condition. The interfce rdii c, d, e nd f re obtined in mnner similr to tht described in Section 3.. Fig.. lstic plstic stresses in luminum cylinder: ) plne-stress, b) generlizedplne-strin. 4. Numericl simultions In this section, numericl simultions with the plne-stress nd generlized-plne-strin models re crried out for typicl cylinder. The objective is to compre the solutions for stresses nd plstic strins of these two models under rdil therml grdient. An luminum cylinder with = 0 mm nd b = 0 mm is considered. The mteril properties of luminum re provided in Tble. The hrdening coefficient K nd strin hrdening exponent n for the cylinder re tken s 58.8 MP nd 0.48 [8]. The temperture difference ruired to initite the yielding t the inner wll of the cylinder is obtined s C for plne stress model nd C for generlized plne strin model. The simultion is crried out for temperture difference, T b T ) = 75 C. With this temperture grdient mechnicl nd therml properties do not vry much [4]. The inner wll is ssumed to be t T = 5 C. According to plne stress model, for T b T ) = 75 C, the cylinder undergoes first stge of elstic plstic deformtion with n inner plstic zone propgting outwrds up to rdius c =.0 mm. However, s per generlized plne strin model, the cylinder undergoes second stge of elstic plstic deformtion dividing the cylinder wll into two inner-plstic, two outer-plstic nd n intermedite-elstic zone. At the inner wll, plstic zone I propgtes outwrds to rdius c =.4973 mm nd plstic zone II propgtes outwrds to rdius d =.05 mm. A plstic zone III propgtes inwrds from the outer rdius to rdius f = mm nd plstic zone IV propgtes inwrds to rdius e = mm. The elstic plstic stresses generted in different zones of the cylinder for plne stress nd generlized plne strin conditions re shown long rdil pth in Fig.. It is observed from Fig. ) tht for plne stress cse, the mximum vlue of hoop stress in the cylinder is generted t the inner rdius of the cylinder nd is tensile in nture. Fig. b) shows tht the mximum vlue of hoop stress exists t the interfce rdius c nd the mximum vlue of xil stress exists t the interfce rdius d s per generlized plne strin model. The mgnitudes of rdil stresses long the rdil pth re smller for plne stress s well s generlized plne strin model. 748
7 Fig. 4. Residul stress distribution in luminum cylinder: ) plne-stress, b) generlizedplne-strin. Fig. 3. Plstic strin distribution in luminum cylinder: ) plne-stress, b) generlizedplne-strin. The plstic prts of rdil, hoop nd xil strins produced in the cylinder re obtined for both plne stress nd generlized plne strin models Fig. 3 ). It is observed from Fig. 3 ) tht the plstic prts of strins re numericlly very smll for plne stress cse. However, Fig. 3 b) shows tht for generlized plne strin cse, lthough the plstic prts of strins re lrger in mgnitude s compred to plne stress cse, their mgnitudes re still smll of the order of 0 3 ). Thus, the mgnitude of uivlent plstic strin generted in the cylinder is lso not very substntil. As the vlue of uivlent stress σ in the plstic zone for yielding depends on the uivlent plstic strin, here the vlue of σ will not devite much from σ Y. The mximum devition of σ from σ Y is only % for the plne stress condition nd 4.87% for the generlized plne strin condition. The mximum uivlent plstic strin occurs t the inner rdius. When the temperture difference T b T ) is removed, i.e., when the cylinder is cooled to room temperture, the residul stresses re set up in the wll of the cylinder. Considering the unloding process to be completely elstic, the residul stresses cn be obtined by subtrcting the thermo-elstic stresses [] from the respective stress utions of elstic plstic zones. The resulting residul stress distribution in the cylinder for the plne stress nd the generlized plne strin condition is shown in Fig. 4. It is observed from Fig. 4 ) nd b) tht the compressive residul hoop stresses re generted t nound the inner rdius of the cylinder. For the sme temperture difference, the gen- erlized plne strin model predicts lrger mgnitude of compressive residul hoop stress t the inner rdius. Fig. 4 b) shows tht for generlized plne strin condition the xil residul stresses generted t nd round the inner rdius of the cylinder re lso compressive. The compressive residul stresses t nound the inner rdius of the cylinder helps in reducing the net mximum stress produced in the cylinder in the next loding stge. This mounts to increse the lod crrying cpcity of the cylinder. The mgnitudes of tensile hoop nd xil stresses re smll t nound the outer wll of the cylinder. 5. Three-dimensionl finite-element modeling of elstic plstic stresses in thick cylinders due to temperture grdient From the numericl simultions of both the plne stress nd generlized plne strin nlyticl models, it is not cler which model is vlid for prticulr length of the cylinder. To ssess the vlidity of the nlyticl models, 3D FM nlysis using ABAQUS 6.0 pckge is crried out. A homogeneous thick-wlled hollow cylinder with inner rdius nd outer rdius b is considered. The inner wll is subjected to temperture T nd the outer wll is subjected to temperture T b such tht T b > T. The problem is solved using ABAQUS stndrd code. Under the rdil temperture difference, T b T ), the stedy stte condition is ssumed in ABAQUS stndrd simultion. The strin hrdening of the mteril during plstic deformtion is considered. For comprison of the FM results with the plne-stress nd generlized-plne-strin models, the nlyticl solution developed in [8] nd the stress solution devel- 749
8 Tble Mesh sensitivity nlysis. Mesh Number of elements in Mximum vlue of Screentime s) Rdil direction Circumferentil direction Axil direction Rdil stress MP) Hoop stress MP) Axil stress MP) Mesh Mesh Mesh Mesh Mesh Mesh Mesh Mesh Mesh Mesh Mesh 9 indicted in the boldfced is the optimum mesh. Tble 3 Comprison of stresses between 3D FM nd nlyticl models. L / b ) L norm of error in stresses L norm of error in stresses between between FM nd plne stress MP) FM nd generlized plne strin MP) σ r σ θ σ z σ r σ θ σ z The boldfced vlues correspond to the cses in which ll stress components show less thn 0% error. oped in Section 3 re used. The detiled nlysis of 3D FM is presented in the following subsections D finite-element modeling A 3D FM model for nlysing thermo-elstic plstic stresses in the cylinder under the rdil temperture difference ws developed using commercil finite element pckge ABAQUS 6.0. A coupled thermomechnicl pproch is used to obtin the stedy-stte therml stress solutions in ABAQUS/Stndrd. In therml utofrettge process, the cylinder is subjected to non-homogeneous elstic plstic deformtion due to rdil therml grdient between the outer nd inner wll of the cylinder. Hence, for the stress solution in therml utofrettge process, coupled temperture-displcement elements of ABAQUS/Stndrd [5] re used. The prt considered for the 3D FM nlysis is thick-wlled solid extruded cylinder with open ends. 5.. Mteril properties An luminum cylinder is considered. The inner nd outer rdii of the cylinder re tken s 0 mm nd 0 mm, respectively. The temperture difference cross the wll thickness of the cylinder is tken s 75 C. The therml stress nlysis during the therml loding nd unloding of cylinder for different lengths is crried out. The mechnicl nd therml properties of luminum re presented in Tble. The model used von Mises criterion with isotropic hrdening Boundry conditions nd mesh genertion The therml stresses in the cylinder re induced due to the rdil therml grdient cross the wll thickness. For nlyzing the therml stress, the Dirichlet temperture boundry conditions re specified t the inner nd outer surfces of the cylinder. The temperture of the inner surfce is prescribed s T = 5 C nd tht of the outer surfce is Fig. 5. Schemtic of 3D prt in ABAQUS long with the boundry conditions. prescribed s T b = 00 C. At ech node, there re three trnsltionl degrees of freedom, viz., rdil displcement u r, circumferentil displcement u θ nd xil displcement u z. In the nlysis, the circumferentil displcement, u θ vnishes. This implies tht the rottion in the circumferentil direction is constrined, but the cylinder is free to expnd xilly ndilly. A schemtic of the 3D prt in ABAQUS long with the therml nd displcement boundry conditions is shown in Fig. 5. An eight-node continuum C3D8T thermlly coupled brick, triliner displcement nd temperture element is used to generte the mesh on the cylinder. The meshed cylinder is shown in Fig. 6. Although the cylinder is modeled s open-ended, the xil stresses re produced in it due to the therml grdient. When the inner surfce of the cylinder is t temperture lower thn tht of the outer surfce, the tendency to expnd in the xil direction is less in the vicinity of the inner surfce thn in the vicinity of the outer one. This cuses tensile xil stresses on the inner side nd compressive xil stresses t the outer side. However, the resultnt xil force due to these stresses is zero. 750
9 Fig. 6. Cylindricl geometry with typicl C3D8T element Mesh sensitivity nlysis A mesh sensitivity nlysis is crried out in order to select n pproprite mesh size. The mesh sensitivity nlysis is presented in Tble bsed on the mximum vlues of rdil, hoop nd xil stresses generted in the cylinder. The length of the cylinder ws tken s 00 mm. It is observed from Tble tht the Mesh 9 with 4 elements in the rdil direction, 3 elements in the circumferentil direction nd 50 elements in the xil direction leds to the convergent mesh. This mesh hs totl of 38,400 elements. Mesh 9 provides less thn % devition in the solution in comprison to Mesh 8 nd Mesh 0. Hence, it is considered s the optimum mesh. The thermo-elstic plstic stresses generted in the cylinder re nlyzed vrying the length of the cylinder using Mesh 9. A comprtive study with the plne stress nd generlized plin strin models is presented for different length to wll-thickness rtios of the cylinder in Section Comprison of 3D finite element elstic plstic therml stresses with plne-stress nd generlized-plne-strin models for different length to wll thickness rtios of cylinder The 3D FM simultions re crried out for different length to wll thickness rtios of the cylinder nd the results re compred with nlyticl plne stress nd generlized plne strin model. The present FM nlysis is crried out in.40-ghz processor nd.60-gb rndom-ccess memory RAM) ACR PC nd it tkes bout h of screen time. However, the nlyticl models provide the solution in less thn 0 min. The comprison of rdil, hoop nd xil stresses between 3D FM nd nlyticl models long the rdil pth for different length to wll thickness rtio is shown in Tble 3 using L norm of error in stresses. The stresses long rdil pth from the 3D ABAQUS finite element model re obtined considering rdil pth t the mid-length of the cylinder. The stresses long rdil pth t the edges become sufficiently low in FM model ne not tken into ccount for the nlysis. It is observed from Tble 3 tht L norm of error in stresses between 3D FM nd generlized plne strin re resonbly smll when the length to wll thickness rtio of the cylinder is greter thn or ul to 6. Thus, the generlized plne strin model is well estblished for L / b ) 6. It is lso observed from Tble 3 tht when L / b ), the L norm of error between 3D FM nd plne stress becomes smller. This shows tht the plne stress nlyticl model is vlid for length to wll thickness rtio, L / b ). The comprison of 3D finite element solutions with the generlized plne strin model for L / b ) = 0 is shown in Fig. 7 ). Fig. 7 b) shows the comprison of 3D finite element Fig. 7. Comprison of elstic plstic stresses between 3D FM nd ) generlized-plnestrin for L/ b ) = 0, b) plne-stress for L/ b ) = 0.5 in luminum cylinder. results with the plne stress model for L / b ) = 0.5 in luminum cylinder Residul stress solutions The whole cylinder is cooled to room temperture 5 C) by imposing the cooling boundry condition in ABAQUS/Stndrd. This genertes the residul stresses in the cylinder. The finite element residul stresses in the cylinder re obtined long rdil pth t the mid-length. The results for long nd short cylinders re depicted in Fig. 8 long with the predictions of generlized plne strin model for L/ b ) = 0 nd plne stress model for L/ b ) = 0.5. It is observed tht the finite element predictions of the residul stresses re in close greement with the nlyticl models. 6. Conclusions The min objective of this work is to compre thermo-elstic plstic plne-stress nd generlized-plne-strin nlyses with 3D FM simultions. For the plne-stress nlysis, the model of Kml nd Dixit [8] is used. For the generlized-plne-strin nlysis, the model of Kml nd Dixit [9] is extended to incorporte the effect of strin hrdening. The FM results re compred with the nlyticl solutions for different length to wll thickness rtio, L / b ), of the cylinder. The comprison suggests the pplicbility of the developed nlyticl models on the bsis of L / b ) rtio of the cylinder. The generlized plne strin nlyticl model provides relistic solution for L / b ) 6. Thus, this 75
10 Acknowledgments Funding from the ngineering nd Physicl Sciences Reserch Council UK) through grnt P/K0836/ nd Deprtment of Science nd Technology Indi) through grnt DST/RC-UK/4-AM/0, project Modeling of Advnced Mterils for Simultion of Trnsformtive Mnufcturing Processes MAST) is grtefully cknowledged. The uthors re lso thnkful to the orgnizers of WCCM XII & APCOM VI 06 Congress for giving the opportunity of orl presenttion of this pper. References Fig. 8. Comprison of residul stresses between 3D FM nd ) generlized-plne-strin for L/ b ) = 0, b) plne-stress for L/ b ) = 0.5 in luminum cylinder. model is suitble for long cylinders such s gun brrels, pressure vessels nd thick pipes subjected to therml grdient. The plne stress nlyticl model provides gooesults for very short cylinders thin disks) such s fstener holes with L / b ). It is to be noted tht the CPU time ruired is much more in cse of FM nlysis compred to the solution obtined by nlyticl models. Thus, the nlyticl models hve their own importnce minly due to computtionl efficiency. The vlidity of the plne-stress nd the generlized-plne-strin ssumptions is inferred bsed on the sme mteril model nd the fixed b/ rtio. This study highlights tht for L / b ) rtio between nd 6, 3D nlysis is ruired. For other types of cylinder, the rnge of L / b ) ruiring 3D nlysis my differ. It will be interesting to crry out similr study for number of cylinder geometries, loding conditions nd mteril model formultions, nd develop generlized criterion for the vlidity of ech model. [] Hill R, Lee H, Tupper SJ. The theory of combined plstic nd elstic deformtion with prticulr reference to thick tube under internl pressure. Proc R Soc A 947;9November): [] Go XL. An exct elsto-plstic solution for n open-ended thick-wlled cylinder of strin-hrdening mteril. Int J Press Vessel Pip 99;5:9 44. [3] Durbn D. Lrge strin solution for pressurized elsto/plstic tubes. J Appl Mech 979;46:8 30. [4] Go XL. An exct elsto-plstic solution for closed-end thick-wlled cylinder of elstic liner-hrdening mteril with lrge strins. Int J Press Vessel Pip 993;56: [5] Bonn R, Hupt P. xct solutions for lrge elstoplstic deformtions of thick- -wlled tube under internl pressure. Int J Plst 995;:99 8. [6] Go XL. lsto-plstic nlysis of n internlly pressurized thick-wlled cylinder using strin grdient plsticity theory. Int J Solids Struct 003;40: [7] Drijni H, Krgrnovin MH, Nghdbdi R. Design of thick-wlled cylinders under internl pressure bsed on elsto-plstic pproch. Mter Des 009;30: [8] Blnd DR. lstoplstic thick-wlled tubes of work-hrdening mteril subject to internl nd externl pressures nd to temperture grdients. J Mech Phys Solids 956;4:09 9. [9] Wong H, Simionescu O. An nlyticl solution of thermoplstic thick-wlled tube subject to internl heting nd vrible pressure, tking into ccount corner flow nd nonzero initil stress. Int J ng Sci 996;34: [0] Orçn Y, rsln AN. Therml stresses in elstic-plstic tubes with temperture-dependent mechnicl nd therml properties. J Therm Stresses 00;4: [] Sdeghin M, Toussi H. lsto-plstic xisymmetric therml stress nlysis of functionlly grded cylindricl vessel. J Bsic Appl Sci Res 0;: [] Anni Y, Rhimi GH. Stress nlysis of thick pressure vessel composed of functionlly grded incompressible hyperelstic mterils. Int J Mech Sci 05;04: 7. [3] Orçn Y. Therml stresses in het generting elstic-plstic cylinder with free ends. Int J ng Sci 994;3: [4] Cowper GR. The elstoplstic thick-wlled sphere subjected to rdil temperture grdient. J Appl Mech 960: [5] Johnson W, Mellor PB. lstic-plstic behviour of thick-wlled spheres of non work-hrdening mteril subject to stedy-stte rdil temperture grdient. Int J Mech Sci 96;4: [6] Drijni H, Krgrnovin MH, Nghdbdi R. Design of sphericl vessels under stedy-stte therml loding using thermo-elsto plstic concept. Int J Press Vessel Pip 009;86:43 5. [7] Zre HR, Drijni H. Strengthening nd design of the liner hrdening thick- -wlled cylinders using the new method of rottionl utofrettge. Int J Mech Sci 07:4 5: 8. [8] Kml SM, Dixit US. Fesibility study of therml utofrettge process. In: Nrynn RG, Dixit US, editors. Advnces in mteril forming nd joining. New Delhi: Springer; 05. p [9] Kml SM, Dixit US. Fesibility study of therml utofrettge of thick-wlled cylinders. ASM J Press Vessel Technol 05;376): [0] Kml SM, Borsiki AC, Dixit US. xperimentl ssessment of residul stresses induced by the therml utofrettge of thick-wlled cylinders. J Strin Anl ng Des 06;5): [] Dixit PM, Dixit US. Modeling of metl forming nd mchining processes: by finite element nd soft computing methods. London: Springer; 008. [] Chkrbrty J. Theory of plsticity. 3rd ed. Burlington: Butterworth-Heinemnn; 006. [3] Nod N, Hetnrski RB, Tnigw Y. Therml stresses. nd ed. New York: Tylor nd Frncis; 003. [4] Zhu XK, Cho YJ. ffects of temperture-dependent mteril properties on welding simultion. Comput Struct 00;80: [5] ABAQUS User s Mnul, Version 6.7, Dssult Systèmes Simuli Corp., Providence, RI, USA,
Studies on Nuclear Fuel Rod Thermal Performance
Avilble online t www.sciencedirect.com Energy Procedi 1 (1) 1 17 Studies on Nucler Fuel od herml Performnce Eskndri, M.1; Bvndi, A ; Mihndoost, A3* 1 Deprtment of Physics, Islmic Azd University, Shirz
More informationAn inverse steady state thermal stresses in a thin clamped circular plate with internal heat generation
Americn Journl of Engineering Reserch (AJER) e-issn : 2320-0847 p-issn : 2320-0936 Volume-02, Issue-10, pp-276-281 www.jer.org Reserch Pper Open Access An inverse stedy stte therml stresses in thin clmped
More information1 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 informationNew 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 informationENERGY-BASED METHOD FOR GAS TURBINE ENGINE DISK BURST SPEED CALCULATION
28 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES ENERGY-BASED METHOD FOR GAS TURBINE ENGINE DISK BURST SPEED CALCULATION Anton N. Servetnik Centrl Institute of Avition Motors, Moscow, Russi servetnik@cim.ru
More informationHPLP310 Biblio_35 Fissures radial intern in a thick cylinder under pressure and thermal loading
defult Titre : HPLP310 - Biblio_35 Fissure rdile interne dns u[...] Dte : 0/09/011 Pge : 1/15 Responsble : TRAN Vn Xun Clé : V7.0.310 Révision : HPLP310 Biblio_35 Fissures rdil intern in thick cylinder
More informationShear and torsion interaction of hollow core slabs
Competitive nd Sustinble Growth Contrct Nº G6RD-CT--6 Sher nd torsion interction of hollow core slbs HOLCOTORS Technicl Report, Rev. Anlyses of hollow core floors December The content of the present publiction
More informationEffects of peripheral drilling moment on delamination using special drill bits
journl of mterils processing technology 01 (008 471 476 journl homepge: www.elsevier.com/locte/jmtprotec Effects of peripherl illing moment on delmintion using specil ill bits C.C. Tso,, H. Hocheng b Deprtment
More informationA 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 informationAvailable online at ScienceDirect. Procedia Engineering 190 (2017 )
Aville online t www.sciencedirect.com ScienceDirect Procedi Engineering 19 (217 ) 237 242 Structurl nd Physicl Aspects of Construction Engineering Numericl Anlysis of Steel Portl Frme Exposed to Fire Lenk
More informationEntropy ISSN
Entropy 006, 8[], 50-6 50 Entropy ISSN 099-4300 www.mdpi.org/entropy/ ENTROPY GENERATION IN PRESSURE GRADIENT ASSISTED COUETTE FLOW WITH DIFFERENT THERMAL BOUNDARY CONDITIONS Abdul Aziz Deprtment of Mechnicl
More informationA Brief Note on Quasi Static Thermal Stresses In A Thin Rectangular Plate With Internal Heat Generation
Americn Journl of Engineering Reserch (AJER) 13 Americn Journl of Engineering Reserch (AJER) e-issn : 3-847 p-issn : 3-936 Volume-, Issue-1, pp-388-393 www.jer.org Reserch Pper Open Access A Brief Note
More informationDETERMINATION OF MECHANICAL PROPERTIES OF NANOSTRUCTURES WITH COMPLEX CRYSTAL LATTICE USING MOMENT INTERACTION AT MICROSCALE
Determintion RevAdvMterSci of mechnicl 0(009) -7 properties of nnostructures with complex crystl lttice using DETERMINATION OF MECHANICAL PROPERTIES OF NANOSTRUCTURES WITH COMPLEX CRYSTAL LATTICE USING
More informationMath 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 informationSUPPLEMENTARY INFORMATION
DOI:.38/NMAT343 Hybrid Elstic olids Yun Li, Ying Wu, Ping heng, Zho-Qing Zhng* Deprtment of Physics, Hong Kong University of cience nd Technology Cler Wter By, Kowloon, Hong Kong, Chin E-mil: phzzhng@ust.hk
More informationConsequently, 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 informationNumerical Study Of Coated Electrical Contacts
Excerpt from the Proceedings of the COMSOL Conference 21 Pris Numericl Study Of Coted Electricl Contcts Per Lindholm Mchine Design KTH Brinellvägen 83 SE-144 Stockholm per@md.kth.se Abstrct: Electricl
More information- 5 - TEST 2. This test is on the final sections of this session's syllabus and. should be attempted by all students.
- 5 - TEST 2 This test is on the finl sections of this session's syllbus nd should be ttempted by ll students. Anything written here will not be mrked. - 6 - QUESTION 1 [Mrks 22] A thin non-conducting
More informationAN EXACT SOLUTION OF MECHANICAL BUCKLING FOR FUNCTIONALLY GRADED MATERIAL BIMORPH CIRCULAR PLATES
Assocition of Metllurgicl Engineers of Serbi AMES Scientific pper UDC: 6.:66.7/.8 AN EXACT SOLUTION OF MECHANICAL BUCKLING FOR FUNCTIONALLY GRADED MATERIAL BIMORPH CIRCULAR PLATES Jfr Eskndri Jm, Mhmood
More informationHT Module 2 Paper solution. Module 2. Q6.Discuss Electrical analogy of combined heat conduction and convection in a composite wall.
HT Module 2 Pper solution Qulity Solutions wwwqulitytutorilin Module 2 Q6Discuss Electricl nlogy of combined het conduction nd convection in composite wll M-16-Q1(c)-5m Ans: It is frequently convient to
More informationMAC-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 informationCBE 291b - Computation And Optimization For Engineers
The University of Western Ontrio Fculty of Engineering Science Deprtment of Chemicl nd Biochemicl Engineering CBE 9b - Computtion And Optimiztion For Engineers Mtlb Project Introduction Prof. A. Jutn Jn
More informationOn the Linear Stability of Compound Capillary Jets
ILASS Americs, th Annul Conference on Liquid Atomiztion nd Spry Systems, Chicgo, IL, My 7 On the Liner Stbility of Compound Cpillry Jets Mksud (Mx) Ismilov, Stephen D Heister School of Aeronutics nd Astronutics,
More informationProblems 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 informationTheoretische 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 information5.2 Volumes: Disks and Washers
4 pplictions of definite integrls 5. Volumes: Disks nd Wshers In the previous section, we computed volumes of solids for which we could determine the re of cross-section or slice. In this section, we restrict
More informationCAPACITORS AND DIELECTRICS
Importnt Definitions nd Units Cpcitnce: CAPACITORS AND DIELECTRICS The property of system of electricl conductors nd insultors which enbles it to store electric chrge when potentil difference exists between
More informationKINEMATICS 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 informationIntro to Nuclear and Particle Physics (5110)
Intro to Nucler nd Prticle Physics (5110) Feb, 009 The Nucler Mss Spectrum The Liquid Drop Model //009 1 E(MeV) n n(n-1)/ E/[ n(n-1)/] (MeV/pir) 1 C 16 O 0 Ne 4 Mg 7.7 14.44 19.17 8.48 4 5 6 6 10 15.4.41
More informationINTRODUCTION. The three general approaches to the solution of kinetics problems are:
INTRODUCTION According to Newton s lw, prticle will ccelerte when it is subjected to unblnced forces. Kinetics is the study of the reltions between unblnced forces nd the resulting chnges in motion. The
More informationPhysics 3323, Fall 2016 Problem Set 7 due Oct 14, 2016
Physics 333, Fll 16 Problem Set 7 due Oct 14, 16 Reding: Griffiths 4.1 through 4.4.1 1. Electric dipole An electric dipole with p = p ẑ is locted t the origin nd is sitting in n otherwise uniform electric
More informationTHE INTERVAL LATTICE BOLTZMANN METHOD FOR TRANSIENT HEAT TRANSFER IN A SILICON THIN FILM
ROMAI J., v.9, no.2(2013), 173 179 THE INTERVAL LATTICE BOLTZMANN METHOD FOR TRANSIENT HEAT TRANSFER IN A SILICON THIN FILM Alicj Piseck-Belkhyt, Ann Korczk Institute of Computtionl Mechnics nd Engineering,
More informationPhysics 201 Lab 3: Measurement of Earth s local gravitational field I Data Acquisition and Preliminary Analysis Dr. Timothy C. Black Summer I, 2018
Physics 201 Lb 3: Mesurement of Erth s locl grvittionl field I Dt Acquisition nd Preliminry Anlysis Dr. Timothy C. Blck Summer I, 2018 Theoreticl Discussion Grvity is one of the four known fundmentl forces.
More informationTHERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION
XX IMEKO World Congress Metrology for Green Growth September 9,, Busn, Republic of Kore THERMAL EXPANSION COEFFICIENT OF WATER FOR OLUMETRIC CALIBRATION Nieves Medin Hed of Mss Division, CEM, Spin, mnmedin@mityc.es
More informationADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS
ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS F. Tkeo 1 nd M. Sk 1 Hchinohe Ntionl College of Technology, Hchinohe, Jpn; Tohoku University, Sendi, Jpn Abstrct:
More informationProgressive failure analysis of compression-loaded composite flat panel with cutout
Interntionl Journl on Theoreticl nd Applied Reserch in Mechnicl Engineering (IJTARME) Progressive filure nlysis of compression-loded composite flt pnel with cutout 1 Guspir S. Mkndr, 2 N.K. Chhpkhne, 3
More informationSimulated Performance of Packed Bed Solar Energy Storage System having Storage Material Elements of Large Size - Part I
The Open Fuels & Energy Science Journl, 2008, 1, 91-96 91 Open Access Simulted Performnce of Pcked Bed Solr Energy Storge System hving Storge Mteril Elements of Lrge Size - Prt I Rnjit Singh *,1, R.P.
More informationAPPLICATIONS 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 informationTHE VISCOELASTIC CONTACT BETWEEN HIGH-ORDER POLYNOMIAL SURFACES
THE VISCOELASTIC CONTACT BETWEEN HIGH-ORDER POLYNOMIAL SURFACES Sergiu Spinu 1,, Deli Cerlinc 1, 1 Deprtment of Mechnics nd Technologies, Stefn cel Mre University of Sucev, 13th University Street, 79,
More informationCHM 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 informationPart I: Basic Concepts of Thermodynamics
Prt I: Bsic Concepts o Thermodynmics Lecture 4: Kinetic Theory o Gses Kinetic Theory or rel gses 4-1 Kinetic Theory or rel gses Recll tht or rel gses: (i The volume occupied by the molecules under ordinry
More informationTheoretical 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 informationRel Gses 1. Gses (N, CO ) which don t obey gs lws or gs eqution P=RT t ll pressure nd tempertures re clled rel gses.. Rel gses obey gs lws t extremely low pressure nd high temperture. Rel gses devited
More informationCONTRIBUTION TO THE EXTENDED DYNAMIC PLANE SOURCE METHOD
CONTRIBUTION TO THE EXTENDED DYNAMIC PLANE SOURCE METHOD Svetozár Mlinrič Deprtment of Physics, Fculty of Nturl Sciences, Constntine the Philosopher University, Tr. A. Hlinku, SK-949 74 Nitr, Slovki Emil:
More informationTravelling Profile Solutions For Nonlinear Degenerate Parabolic Equation And Contour Enhancement In Image Processing
Applied Mthemtics E-Notes 8(8) - c IN 67-5 Avilble free t mirror sites of http://www.mth.nthu.edu.tw/ men/ Trvelling Profile olutions For Nonliner Degenerte Prbolic Eqution And Contour Enhncement In Imge
More informationStudy on the Calculation of Magnetic Force Based on the Equivalent Magnetic Charge Method
Avilble online t www.sciencedirect.com Physics Procedi 4 () 9 97 Interntionl Conference on Applied Physics nd Industril Engineering Study on the Clcultion of Mgnetic Force Bsed on the Equivlent Mgnetic
More informationFig. 1. Open-Loop and Closed-Loop Systems with Plant Variations
ME 3600 Control ystems Chrcteristics of Open-Loop nd Closed-Loop ystems Importnt Control ystem Chrcteristics o ensitivity of system response to prmetric vritions cn be reduced o rnsient nd stedy-stte responses
More informationThin and Thick Cylinders and Spheres
CHAPTR 8 Thin nd Thick Cylinders nd Spheres Prolem. A shell.5 m long nd m dimeter, is sujected to n internl pressure of. N/mm. If the thickness of the shell is 0 mm find the circumferentil nd longitudinl
More informationMultiple 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 information63. 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 informationJack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah
1. Born-Oppenheimer pprox.- energy surfces 2. Men-field (Hrtree-Fock) theory- orbitls 3. Pros nd cons of HF- RHF, UHF 4. Beyond HF- why? 5. First, one usully does HF-how? 6. Bsis sets nd nottions 7. MPn,
More informationMeasuring Electron Work Function in Metal
n experiment of the Electron topic Mesuring Electron Work Function in Metl Instructor: 梁生 Office: 7-318 Emil: shling@bjtu.edu.cn Purposes 1. To understnd the concept of electron work function in metl nd
More informationNUMERICAL 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 informationA - 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 informationThermal Performance of Electrocaloric Refrigeration using Thermal Switches of Fluid Motion and Changing Contact Conductance
Americn Journl of Physics nd Applictions 2016; 4(5): 134-139 http://www.sciencepublishinggroup.com/j/jp doi:.11648/j.jp.20160405.12 ISSN: 2330-4286 (Print); ISSN: 2330-4308 (Online) Therml Performnce of
More informationMultiple 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 informationProf. Anchordoqui. Problems set # 4 Physics 169 March 3, 2015
Prof. Anchordoui Problems set # 4 Physics 169 Mrch 3, 15 1. (i) Eight eul chrges re locted t corners of cube of side s, s shown in Fig. 1. Find electric potentil t one corner, tking zero potentil to be
More informationSolution 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 informationTerminal Velocity and Raindrop Growth
Terminl Velocity nd Rindrop Growth Terminl velocity for rindrop represents blnce in which weight mss times grvity is equl to drg force. F 3 π3 ρ L g in which is drop rdius, g is grvittionl ccelertion,
More informationNew data structures to reduce data size and search time
New dt structures to reduce dt size nd serch time Tsuneo Kuwbr Deprtment of Informtion Sciences, Fculty of Science, Kngw University, Hirtsuk-shi, Jpn FIT2018 1D-1, No2, pp1-4 Copyright (c)2018 by The Institute
More informationLecture 6: Singular Integrals, Open Quadrature rules, and Gauss Quadrature
Lecture notes on Vritionl nd Approximte Methods in Applied Mthemtics - A Peirce UBC Lecture 6: Singulr Integrls, Open Qudrture rules, nd Guss Qudrture (Compiled 6 August 7) In this lecture we discuss the
More informationKINEMATICS 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 informationragsdale (zdr82) HW2 ditmire (58335) 1
rgsdle (zdr82) HW2 ditmire (58335) This print-out should hve 22 questions. Multiple-choice questions my continue on the next column or pge find ll choices before nswering. 00 0.0 points A chrge of 8. µc
More informationLecture 14: Quadrature
Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl
More informationThe International Association for the Properties of Water and Steam. Release on the Ionization Constant of H 2 O
IAPWS R-7 The Interntionl Assocition for the Properties of Wter nd Stem Lucerne, Sitzerlnd August 7 Relese on the Ioniztion Constnt of H O 7 The Interntionl Assocition for the Properties of Wter nd Stem
More informationA027 Uncertainties in Local Anisotropy Estimation from Multi-offset VSP Data
A07 Uncertinties in Locl Anisotropy Estimtion from Multi-offset VSP Dt M. Asghrzdeh* (Curtin University), A. Bon (Curtin University), R. Pevzner (Curtin University), M. Urosevic (Curtin University) & B.
More informationPhysics 202, Lecture 14
Physics 202, Lecture 14 Tody s Topics Sources of the Mgnetic Field (Ch. 28) Biot-Svrt Lw Ampere s Lw Mgnetism in Mtter Mxwell s Equtions Homework #7: due Tues 3/11 t 11 PM (4th problem optionl) Mgnetic
More informationThe Moving Center of Mass of a Leaking Bob
The Moving Center of Mss of Leking Bob rxiv:1002.956v1 [physics.pop-ph] 21 Feb 2010 P. Arun Deprtment of Electronics, S.G.T.B. Khls College University of Delhi, Delhi 110 007, Indi. Februry 2, 2010 Abstrct
More informationFEM ANALYSIS OF ROGOWSKI COILS COUPLED WITH BAR CONDUCTORS
XIX IMEKO orld Congress Fundmentl nd Applied Metrology September 6 11, 2009, Lisbon, Portugl FEM ANALYSIS OF ROGOSKI COILS COUPLED ITH BAR CONDUCTORS Mirko Mrrcci, Bernrdo Tellini, Crmine Zppcost University
More informationLesson 8. Thermomechanical Measurements for Energy Systems (MENR) Measurements for Mechanical Systems and Production (MMER)
Lesson 8 Thermomechnicl Mesurements for Energy Systems (MEN) Mesurements for Mechnicl Systems nd Production (MME) A.Y. 205-6 Zccri (ino ) Del Prete Mesurement of Mechnicl STAIN Strin mesurements re perhps
More informationKirchhoff 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 informationLecture 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 informationContact Analysis on Large Negative Clearance Four-point Contact Ball Bearing
Avilble online t www.sciencedirect.co rocedi ngineering 7 0 74 78 The Second SR Conference on ngineering Modelling nd Siultion CMS 0 Contct Anlysis on Lrge Negtive Clernce Four-point Contct Bll Bering
More informationBasic model for traffic interweave
Journl of Physics: Conference Series PAPER OPEN ACCESS Bsic model for trffic interweve To cite this rticle: Ding-wei Hung 25 J. Phys.: Conf. Ser. 633 227 Relted content - Bsic sciences gonize in Turkey!
More informationMath 1B, lecture 4: Error bounds for numerical methods
Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the
More informationReading from Young & Freedman: For this topic, read the introduction to chapter 24 and sections 24.1 to 24.5.
PHY1 Electricity Topic 5 (Lectures 7 & 8) pcitors nd Dielectrics In this topic, we will cover: 1) pcitors nd pcitnce ) omintions of pcitors Series nd Prllel 3) The energy stored in cpcitor 4) Dielectrics
More informationSolutions to Problems in Merzbacher, Quantum Mechanics, Third Edition. Chapter 7
Solutions to Problems in Merzbcher, Quntum Mechnics, Third Edition Homer Reid April 5, 200 Chpter 7 Before strting on these problems I found it useful to review how the WKB pproimtion works in the first
More informationChapter 6 Electrostatic Boundary Value Problems. Dr. Talal Skaik
Chpter 6 Electrosttic Boundry lue Problems Dr. Tll Skik 1 1 Introduction In previous chpters, E ws determined by coulombs lw or Guss lw when chrge distribution is known, or potentil is known throughout
More informationMath 8 Winter 2015 Applications of Integration
Mth 8 Winter 205 Applictions of Integrtion Here re few importnt pplictions of integrtion. The pplictions you my see on n exm in this course include only the Net Chnge Theorem (which is relly just the Fundmentl
More information1 Online Learning and Regret Minimization
2.997 Decision-Mking in Lrge-Scle Systems My 10 MIT, Spring 2004 Hndout #29 Lecture Note 24 1 Online Lerning nd Regret Minimiztion In this lecture, we consider the problem of sequentil decision mking in
More informationEnergy creation in a moving solenoid? Abstract
Energy cretion in moving solenoid? Nelson R. F. Brg nd Rnieri V. Nery Instituto de Físic, Universidde Federl do Rio de Jneiro, Cix Postl 68528, RJ 21941-972 Brzil Abstrct The electromgnetic energy U em
More informationPhysics 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 informationThermal Diffusivity. Paul Hughes. Department of Physics and Astronomy The University of Manchester Manchester M13 9PL. Second Year Laboratory Report
Therml iffusivity Pul Hughes eprtment of Physics nd Astronomy The University of nchester nchester 3 9PL Second Yer Lbortory Report Nov 4 Abstrct We investigted the therml diffusivity of cylindricl block
More informationSystem Identification with Noisy Data
System Identifiction with Noisy Dt U K Dewngn & S V Bri 2 Dept of Civil Engineering, Ntionl Institute of Technology, Ripur, Indi 2 Dept of Civil Engineering, IIT Khrgpur, Khrgpur -72 302, Indi E-mil :
More informationPrediction of tool-chip contact length using a new slip-line solution for orthogonal cutting
Interntionl Journl of Mchine Tools & Mnufcture 43 (2003) 1209 1215 Prediction of tool-chip contct length using new slip-line solution for orthogonl cutting Andrey Toropov, Sung-im Ko School of Mechnicl
More information221B Lecture Notes WKB Method
Clssicl Limit B Lecture Notes WKB Method Hmilton Jcobi Eqution We strt from the Schrödinger eqution for single prticle in potentil i h t ψ x, t = [ ] h m + V x ψ x, t. We cn rewrite this eqution by using
More informationPlates 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 informationJackson 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 informationMath 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 informationMotion of Electrons in Electric and Magnetic Fields & Measurement of the Charge to Mass Ratio of Electrons
n eperiment of the Electron topic Motion of Electrons in Electric nd Mgnetic Fields & Mesurement of the Chrge to Mss Rtio of Electrons Instructor: 梁生 Office: 7-318 Emil: shling@bjtu.edu.cn Purposes 1.
More informationAcceptance Sampling by Attributes
Introduction Acceptnce Smpling by Attributes Acceptnce smpling is concerned with inspection nd decision mking regrding products. Three spects of smpling re importnt: o Involves rndom smpling of n entire
More informationElectromagnetism Answers to Problem Set 10 Spring 2006
Electromgnetism 76 Answers to Problem Set 1 Spring 6 1. Jckson Prob. 5.15: Shielded Bifilr Circuit: Two wires crrying oppositely directed currents re surrounded by cylindricl shell of inner rdius, outer
More informationVadose Zone Hydrology
Objectives Vdose Zone Hydrology 1. Review bsic concepts nd terminology of soil physics. 2. Understnd the role of wter-tble dynmics in GW-SW interction. Drcy s lw is useful in region A. Some knowledge of
More informationHomework Assignment 3 Solution Set
Homework Assignment 3 Solution Set PHYCS 44 6 Ferury, 4 Prolem 1 (Griffiths.5(c The potentil due to ny continuous chrge distriution is the sum of the contriutions from ech infinitesiml chrge in the distriution.
More informationTANDEM QUEUE WITH THREE MULTISERVER UNITS AND BULK SERVICE WITH ACCESSIBLE AND NON ACCESSBLE BATCH IN UNIT III WITH VACATION
Indin Journl of Mthemtics nd Mthemticl Sciences Vol. 7, No., (June ) : 9-38 TANDEM QUEUE WITH THREE MULTISERVER UNITS AND BULK SERVICE WITH ACCESSIBLE AND NON ACCESSBLE BATCH IN UNIT III WITH VACATION
More informationFreely propagating jet
Freely propgting jet Introduction Gseous rectnts re frequently introduced into combustion chmbers s jets. Chemicl, therml nd flow processes tht re tking plce in the jets re so complex tht nlyticl description
More informationExplain 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 informationThe 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 informationLinear and Non-linear Feedback Control Strategies for a 4D Hyperchaotic System
Pure nd Applied Mthemtics Journl 017; 6(1): 5-13 http://www.sciencepublishinggroup.com/j/pmj doi: 10.11648/j.pmj.0170601.1 ISSN: 36-9790 (Print); ISSN: 36-981 (Online) Liner nd Non-liner Feedbck Control
More informationLinear static analysis of perforated plates with round and staggered holes under their self-weights
Liner sttic nlysis of perforted pltes with round nd stggered holes under their self-weights Mustf Hlûk Srçoğlu*, Uğur Albyrk Online Publiction Dte: 8 Nov 2015 URL: http://www.jresm.org/rchive/resm2015.25me0910.html
More information