Thermoelastic Problem of a Long Annular Multilayered Cylinder

Similar documents
A Study about One-Dimensional Steady State. Heat Transfer in Cylindrical and. Spherical Coordinates

COMPLEMENTARY ENERGY METHOD FOR CURVED COMPOSITE BEAMS

5-99C The Taylor series expansion of the temperature at a specified nodal point m about time t i is

Rotating Variable-Thickness Inhomogeneous Cylinders: Part II Viscoelastic Solutions and Applications

Chapter 8. Linear Momentum, Impulse, and Collisions

Analytical and Numerical Solutions for a Rotating Annular Disk of Variable Thickness

GENERALIZATION OF AN IDENTITY INVOLVING THE GENERALIZED FIBONACCI NUMBERS AND ITS APPLICATIONS

gravity r2,1 r2 r1 by m 2,1

Optimal Design of Step Stress Partially Accelerated Life Test under Progressive Type-II Censored Data with Random Removal for Gompertz Distribution

ANALYSIS OF AXIAL LOADED PILE IN MULTILAYERED SOIL USING NODAL EXACT FINITE ELEMENT MODEL

THE EQUIVALENCE OF GRAM-SCHMIDT AND QR FACTORIZATION (page 227) Gram-Schmidt provides another way to compute a QR decomposition: n

Multistage Median Ranked Set Sampling for Estimating the Population Median

PHYS 705: Classical Mechanics. Derivation of Lagrange Equations from D Alembert s Principle

On Accurate Stress Determination in Laminated Finite Length Cylinders Subjected to Thermo Elastic Load

PHYS 1443 Section 003 Lecture #21

Design and Flow Parameters Calculation of the Turbomachine Channels

Thermal behavior of friction clutch disc based on uniform pressure and uniform wear assumptions

BOUNDARY-ONLY INTEGRAL EQUATION APPROACH BASED ON POLYNOMIAL EXPANSION OF PLASMA CURRENT PROFILE TO SOLVE THE GRAD-SHAFRANOV EQUATION

24-2: Electric Potential Energy. 24-1: What is physics

The Finite Strip Method (FSM) 1. Introduction

V. Principles of Irreversible Thermodynamics. s = S - S 0 (7.3) s = = - g i, k. "Flux": = da i. "Force": = -Â g a ik k = X i. Â J i X i (7.

Fatigue equivalent loads for visualization of multimodal dynamic simulations

Set of square-integrable function 2 L : function space F

UNIT10 PLANE OF REGRESSION

Some Approximate Analytical Steady-State Solutions for Cylindrical Fin

2/24/2014. The point mass. Impulse for a single collision The impulse of a force is a vector. The Center of Mass. System of particles

Chapter 23: Electric Potential

Molecular Dynamic Simulations of Nickel Nanowires at Various Temperatures

CSJM University Class: B.Sc.-II Sub:Physics Paper-II Title: Electromagnetics Unit-1: Electrostatics Lecture: 1 to 4

Thermodynamics of solids 4. Statistical thermodynamics and the 3 rd law. Kwangheon Park Kyung Hee University Department of Nuclear Engineering

Energy in Closed Systems

Optimal System for Warm Standby Components in the Presence of Standby Switching Failures, Two Types of Failures and General Repair Time

A New Approach for Deriving the Instability Potential for Plates Based on Rigid Body and Force Equilibrium Considerations

Asymptotic Solutions of the Kinetic Boltzmann Equation and Multicomponent Non-Equilibrium Gas Dynamics

Potential Theory. Copyright 2004

Rigid Bodies: Equivalent Systems of Forces

Fundamental principles

Please initial the statement below to show that you have read it

A. Proofs for learning guarantees

A Study of C-Reducible Finsler Space. With Special Cases

MHD Oscillatory Flow in a Porous Plate

7/1/2008. Adhi Harmoko S. a c = v 2 /r. F c = m x a c = m x v 2 /r. Ontang Anting Moment of Inertia. Energy

User-friendly model of heat transfer in. preheating, cool down and casting

Generating Functions, Weighted and Non-Weighted Sums for Powers of Second-Order Recurrence Sequences

Transport Coefficients For A GaAs Hydro dynamic Model Extracted From Inhomogeneous Monte Carlo Calculations

Properties of the magnetotelluric frequency-normalised impedance over a layered medium

Integral Vector Operations and Related Theorems Applications in Mechanics and E&M

Chapter 3 Waves in an Elastic Whole Space. Equation of Motion of a Solid

Test 1 phy What mass of a material with density ρ is required to make a hollow spherical shell having inner radius r i and outer radius r o?

A. Thicknesses and Densities

If there are k binding constraints at x then re-label these constraints so that they are the first k constraints.

Physics 207 Lecture 16

Stellar Astrophysics. dt dr. GM r. The current model for treating convection in stellar interiors is called mixing length theory:

Physics 1501 Lecture 19

Exact Three-Dimensional Elasticity Solution for Buckling of Beam of a Circular Cross-Section *

Slobodan Lakić. Communicated by R. Van Keer

Mechanics Physics 151

COMPUTATIONAL METHODS AND ALGORITHMS Vol. I - Methods of Potential Theory - V.I. Agoshkov, P.B. Dubovski

4.4 Continuum Thermomechanics

An Approach to Inverse Fuzzy Arithmetic

Asymptotic Waves for a Non Linear System

Integral Transforms and Dual Integral Equations to Solve Heat Equation with Mixed Conditions

Journal of Naval Science and Engineering 2015, Vol.11, No.1, pp FINITE DIFFERENCE MODEL OF A CIRCULAR FIN WITH RECTANGULAR PROFILE.

Chapter 12 Equilibrium and Elasticity

ON THE FRESNEL SINE INTEGRAL AND THE CONVOLUTION

Distinct 8-QAM+ Perfect Arrays Fanxin Zeng 1, a, Zhenyu Zhang 2,1, b, Linjie Qian 1, c

3.1 Electrostatic Potential Energy and Potential Difference

Engineering Mechanics. Force resultants, Torques, Scalar Products, Equivalent Force systems

Chapter Fifiteen. Surfaces Revisited

Modeling and Analysis of Laminated Composite Plates Embedded with Shape Memory Alloy Wires

Budding yeast colony growth study based on circular granular cell

ALL QUESTIONS ARE WORTH 20 POINTS. WORK OUT FIVE PROBLEMS.

10/15/2013. PHY 113 C General Physics I 11 AM-12:15 PM MWF Olin 101

Chapter I Matrices, Vectors, & Vector Calculus 1-1, 1-9, 1-10, 1-11, 1-17, 1-18, 1-25, 1-27, 1-36, 1-37, 1-41.

Dynamic Performance, System Identification and Sensitivity Analysis of the Ladder Tracks. Ontario, Canada

DRBEM Applied to the 3D Helmholtz Equation and Its Particular Solutions with Various Radial Basis Functions

TEST-03 TOPIC: MAGNETISM AND MAGNETIC EFFECT OF CURRENT Q.1 Find the magnetic field intensity due to a thin wire carrying current I in the Fig.

BALANCING OF ROTATING MASSES

General Variance Covariance Structures in Two-Way Random Effects Models

Analysis of the magnetic field, force, and torque for two-dimensional Halbach cylinders

ISSN: [Ramalakshmi* et al., 6(1): January, 2017] Impact Factor: 4.116

Quantum Particle Motion in Physical Space

The Forming Theory and the NC Machining for The Rotary Burs with the Spectral Edge Distribution

Groupoid and Topological Quotient Group

XII.3 The EM (Expectation-Maximization) Algorithm

Optimization Methods: Linear Programming- Revised Simplex Method. Module 3 Lecture Notes 5. Revised Simplex Method, Duality and Sensitivity analysis

APPLICATIONS OF SEMIGENERALIZED -CLOSED SETS

1. Starting with the local version of the first law of thermodynamics q. derive the statement of the first law of thermodynamics for a control volume

1. A body will remain in a state of rest, or of uniform motion in a straight line unless it

MODFLOW and More 2003: Understanding through Modeling Conference Proceedings, Poeter, Zheng, Hill & Doherty

International Journal of Solids and Structures

Physics 2A Chapter 11 - Universal Gravitation Fall 2017

BALANCING OF ROTATING MASSES

N = N t ; t 0. N is the number of claims paid by the

Efficiency of the principal component Liu-type estimator in logistic

Large scale magnetic field generation by accelerated particles in galactic medium

Implementation in the ANSYS Finite Element Code of the Electric Vector Potential T-Ω,Ω Formulation

Scalars and Vectors Scalar

Unconventional double-current circuit accuracy measures and application in twoparameter

Design and Simulation of a Three-Phase Electrostatic Cylindrical Rotary Micromotor

Transcription:

Wold Jounal of Mechancs, 3, 3, 6- http://dx.do.og/.436/w.3.35a Publshed Onlne August 3 (http://www.scp.og/ounal/w) Theoelastc Poble of a Long Annula Multlayeed Cylnde Y Hsen Wu *, Kuo-Chang Jane Depatent of Infoaton Manageent, Oental Insttute of Technology, Tape, Chnese Tape Depatent of Appled Matheatcs, atonal Chung Hsng Unvesty, Tachung, Chnese Tape Eal: * yhwu@al.ot.edu.tw Receved May, 3; evsed June, 3; accepted June 9, 3 Copyght 3 Y Hsen Wu, Kuo-Chang Jane. Ths s an open access atcle dstbuted unde the Ceatve Coons Attbuton Lcense, whch pets unestcted use, dstbuton, and epoducton n any edu, povded the ognal wok s popely cted. ABSTRACT Theoelastc tansent esponse of ultlayeed annula cylndes of nfnte lengths subected to known nne pessue and oute sufaces coolng ae consdeed. A ethod based on the Laplace tansfoaton and fnte dffeence ethod has been developed to analyze the theoelastcty poble. Usng the Laplace tansfo wth espect to te, the geneal solutons of the govenng equatons ae obtaned n tansfo doan. The soluton s obtaned by usng the atx slaty tansfoaton and nvese Laplace tansfo. Solutons fo the tepeatue and theal stess dstbutons n a tansent state wee obtaned. It was found that the tepeatue dstbuton, the dsplaceent and the theal stesses change slghtly as te nceases. Keywods: Theoelastc; Multlayeed Annula Cylndes; Laplace Tansfoaton; Fnte Dffeence Method. Intoducton * Coespondng autho. A theal poble ases when the coposed ateals ae geneated by a sudden change n tepeatue. Shell stuctues ae wdely used n contepoay ndustes, so we ust take cae of the theal poble. The shell stuctues ay be affected due to the pessue change o the vaous tepeatue dstbutons. It s necessay to solve fo tepeatue o pessue at fst. The dynac theoelastc esponse of ccula shell apdly change of theal envonents s potant fo the desgn of any engneeng stuctues. Due to the coplexty of the govenng equatons and the atheatcal dffcultes assocated wth the soluton, seveal splfcatons have been used. Fo exaple, Sheef and Anwa [] dscussed the poble of an annula nfntely long elastc ccula. They have neglected both the neta tes and the elaxaton effects of the poble. Sheef and Anwa [] consdeed the theoelastcty poble of an nfntely long annula cylnde coposed of two dffeent ateals wth axal syety. The soluton was obtaned n the Laplace tansfo doan by usng the potental functon appoach. The pesent wok deals wth the one-densonal quasstatc coupled theoelastc pobles of an nfntely long annula ultlayeed cylnde coposed of ultlayeed dffeent ateals. The edu has a pessue at the nne laye, the tepeatue to be heated at the oute laye, wthout body foces and ntenal heat geneaton. Devatves ae appoxated by cental dffeences esultng n an algebac epesentaton of the patal dffeental equaton. By takng the Laplace tansfo wth espect to te, the geneal solutons n the tansfo doan ae fst obtaned. The fnal solutons n the eal doan can be obtaned by nvetng the Laplace tansfo.. Foulaton Ths wok deals wth the one-densonal, quas-statc coupled, theoelastc pobles of an nfntely long annula cylnde coposed of ultlayeed lanated ateals wth axal syety unde the followng assuptons: ) Mateals of each laye ae assued to be non-hoogeneous; ) Defoaton and stan satsfy the Hooke s law and sall stan theoy; 3) The coposte cylnde s constucted of ultlayeed lanates bonded togethe pefectly; 4) The edu s ntally undstubed, and wthout body foces and ntenal heat souces; 5) The edu s appled by a foce, whch s the functon of te; 6) The tepeatue at nne laye and oute laye ae the functons of te. We now consde an nfntely long annula cylnde Copyght 3 ScRes.

Y. H. WU, K.-C. JAE 7 ade of ultple layes of dffeent ateals. The nne and oute ad of the cylnde ae denoted by and, o espectvely. The ultlayeed coposte s assued to be heated suddenly at the nne and oute suface unde tepeatues f and f espectvely. The tansent heat conducton equaton fo the th laye n densonal fo can be wtten as (see Equaton () below) whee E and E n whch U s the adal coponent of dsplaceent, s adus, Cv and ae specfc heat and densty of ateal, and ae the Posson s ato k, k ae adal, ccufeental theal conductvty,, ae adal and ccufeental theal expanson coeffcent, E, E ae adal and ccufeental Young s odulus,, ae the tepeatue, efeence tepeatue, and s te, espectvely. If the body foces ae absent, the equaton of equlbu fo a cylnde along the adal decton can be wtten as U E U E U E E E E () The stess-dsplaceent elatons ae E U E U (3) E U E U, (4) whee, ae Lae s constant, adal and ccufeental stesses espectvely. Let the bounday condtons of ultlayeed cylnde be at at R t U, e ct t P f at R o, e ct out f t whee f, f, P, ae nne and oute suoundng tepeatues, ntal nne pessue, the ntal tepeatue at the oute laye espectvely. The non-densonal vaables ae defned as follows: T k cos k sn k cos k sn a Cv Cv ksn k cos k cos k sn b Cv Cv sn cos cos sn w Cv Cv k cos k sn t Cv R cos sn u U R Cv R f E E e E E cos sn g Cv cos sn E Cv E h E cos sn Q Cv E cos sn Q Cv Q 3 E cos sn R Cv E cos sn R Cv U k cos k sn k sn k cos Cv cos sn () U sn cos Copyght 3 ScRes.

8 Y. H. WU, K.-C. JAE 3R whee T, t,, u,, ae non-densonal tepeatue, te, adus, dsplaceent, adal stess and ccufeental stess fo the th laye espectvely. Substtutng the nondensonal quanttes nto the govenng Equatons ()-(4), the tansent heat conducton equaton and stess-dsplaceent elatons have the followng nondensonal fo: b u w u a T (5) t t t u e u u T T u u Q Q 3QT u u R R 3RT f g h 3. Coputatonal Pocedues (6) (7) (8) Applyng cental dffeence n Equatons (5)-(8), we ave at the followng dscetzed equatons: u u T T T T T T w u t t a b t t (9) u u u u u T T e f u g T () h whee u u u Q Q 3 u u u R R 3 QT () RT () ae de- and,,,. The Laplace tansfo of a functon fned by st s L t t e d t t Take the Laplace tansfo fo Equatons (9)-(), we obtan the followng equatons: T T T T T w a b T st ( u su ) n, n, u, n su u, n su (3) u u u u u T T T (4) e f u g h Let the suface of the cylndcal nne suface be stess fee and subect to a te-dependent tepeatue. Afte takng Laplace tansfoaton, the bounday condtons n tansfoed doan becoe P s, T fs s c at ; at out. s, T s s c and the nteface condtons ae as follows: u, s u, s whee s, s,,, q s q, s T s T, s,,,. Substtutng the bounday condtons and the nteface condtons nto Equatons (3), (4), we obtan the fol- Copyght 3 ScRes.

Y. H. WU, K.-C. JAE 9 lowng equaton n atx fo (see Equatons (5) below) whee 3Q a B Q a B 3Q Q w E Q Q C a F D X Y Z G E A b D 3Q a b Q k B w E 3Q a Q F 3Q w Q Q k Q k X Y Z X Q Q P Q 3 Y 3Q a b Q k 3Q a b Q k G Z G 3Q b a P T Q Q whee denotes the last laye, k the last pont, and denotes t h laye fo,3,,, (see Equaton (6) below) whee: g h H I K L M J e f e,3,, g B C T X Y Z A B C T X Y Z si s c s A B C T X Y Z A B T X Y Z E F u G D E F u G s D E F u G D E u G I J T L M u H I J K L M u T H I J T K L M u H I T K L u (5) (6) Copyght 3 ScRes.

Y. H. WU, K.-C. JAE Equatons (5) and (6) can be ewtten n the followng atx fos M sit su X s c (7) Y Z G s RT Qu (8) whee the atx M,, R and Q ae the coespondng atx n Equatons (5) and (6). Substtutng Equaton (7) nto (8), we have AsIT B s c (9) C DF s whee A Q R M B Q R X C Q R Y D Q R Z F Q R G Snce the atx A s a nonsngula eal atx, the atx A possesses a set of lnealy ndependent egenvectos, hence the atx A s dagonalzable. Thee exst a nonsngula tanston atx P such that P A P daga, that s, the atces A and dag A ae sla, whee the atx dag A s a dagonal atx wth eleents,,,, whee s the egenvalue of atx A. The equaton can be obtaned as P APsIP T P B s c () P C P DP F s Equaton () can be ewtten as s c s whee dag A si T B C D F () and T P T, B P B C P C, D P D F P F Fo Equaton (), the followng solutons can be obtaned edately. F B C D T s s sc s s s () By applyng the nvese Laplace tansfos to Equaton (), we get the soluton T. The egenvalue, egenvecto and nvese Laplace tansfo of atx [A] can be solved by applyng the IMSL MATH/LIBRARY suboutnes. Afte we have obtaned T, then we can use Equatons (3) and (4) to obtan the solutons T and u T PT u Q RT (3) (4) Substtutng T and u nto Equatons () and (), we obtan the adal and ccufeental stesses. 4. uecal Results and Dscussons In ths secton, we pesent soe nuecal esults of the tepeatue dstbuton n a long ultlayeed coposte hollow cylnde, and dsplaceent and theal stesses unde tepeatue changes. The nne and oute ad of the cylnde ae assued to be. and 4.5 espectvely. Fo an nfntely long annula ultlayeed cylnde, the geoety and ateal quanttes of the cylnde (n the case of thee layes, laye :E = 58E6, k =, =., =.8E 6, =.95, C v =.3 and laye :E = 3E6, k =, =.35, =.3E 6, =.53, C v =.5 and laye 3 : E = E6, k = 7, =., =.8E 6, =.9, C v =.7 ; n the case of fve layes, laye : E = 58E6, k =, =., =.8E 6, =.95, C v =.3 and laye : E = 3E6, k =, =.35, =.3E 6, =.53, C v =.5 and laye 3:E = E6, k = 7, =., =.8E 6, =.9, C v =.7 and laye 4:E = 3E6, k =, =.35, =.3E 6, =.53, C v =.5 and laye 5 : E = E6, k = 7, =., =.8E 6, =.9, C v =.7). Each laye s assued to have a dffeent thckness (n the case of thee layes, =.5, =.5 and 3 =.5; n the case of fve layes, =., =.5, 3 =., 4 =.5 and 5 =.5). The pessue of the nne suface s assued to be P =.5E6. The constant coeffcent c = c =.. The tepeatue at nne suface s assued to be 3, at oute Copyght 3 ScRes.

Y. H. WU, K.-C. JAE suface whch s a functon of te s assued to be to. Fgues -4 show soe nuecal esults of thee and fve layeed cylndes at te step t =.5,,, 5 and. Fgues and show the tepeatue dstbutons along adal decton fo 3 and 5 layes case. Because of the dffeence n theal conductvty and the effect of the oute laye s to be heated. As te s sall, say t =.5, the oute laye tepeatue whch s to be heated s not so oe, so the dstbuton deceasng at fst and then nceasng. Fgues 3 and 4 show the dsplaceent along the adal decton. The axu dsplaceent occued at the nteface of fst and second layes. Fgues 5 and 6 show the adal stess dstbuton along the adal decton. Fgues 7 and 8 show the ccufeental stess along the ccufeental decton. Fgue 3. Radal dsplaceent dstbuton along adal decton fo 3 layes case. Fgue. Tepeatue dstbuton along adal decton fo 3 layes case. Fgue 4. Radal dsplaceent dstbuton along adal decton fo 5 layes case. Fgue. Tepeatue dstbuton along adal decton fo 5 layes case. Fgue 5. Radal stess dstbuton along adal decton fo 3 layes case. Copyght 3 ScRes.

Y. H. WU, K.-C. JAE Fgue 6. Radal stess dstbuton along adal decton fo 5 layes case. Fgue 8. Ccufeental stess dstbuton along adal decton fo 5 layes case. dstbutons have been obtaned, all of whch can be used to desgn useful stuctues o achnes fo engneeng applcatons. Thee s no lt to the nube of annula layes n a cylnde. Exeplfyng nuecal esults fo thee- and fve- layeed cylndes at dffeent te steps have been pesented. The dscontnuty n ccufeental stess at each nteface was found. It was found that the tepeatue dstbuton, the dsplaceent and the theal stesses vay slghtly as the te nceases. Fgue 7. Ccufeental stess dstbuton along adal decton fo 3 layes case. A ethod based on the fnte dffeence and Laplace tansfoaton has been developed to obtan nuecal esults. The tepeatue, dsplaceent and theal stess REFERECES [] H. H. Sheef and M.. Anwa, Poble n Genealzed Theoelastcty, Jounal of Theal Stesses, Vol. 9, o., 986, pp. 65-8. do:.8/4957386896895 [] H. H. Sheef and M.. Anwa, A Poble n Genealzed Theoelastcty fo an Infntely Long Annula Cylnde Coposed of Two Dffeent Mateals, Acta Mechanca, Vol. 8, o. -, 989, pp. 37-49. do:.7/bf7885 Copyght 3 ScRes.