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

Similar documents
COMPLEMENTARY ENERGY METHOD FOR CURVED COMPOSITE BEAMS

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

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

Multistage Median Ranked Set Sampling for Estimating the Population Median

Molecular Dynamic Simulations of Nickel Nanowires at Various Temperatures

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

Thermoelastic Problem of a Long Annular Multilayered Cylinder

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

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

Scalars and Vectors Scalar

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

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

A. Thicknesses and Densities

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

Energy in Closed Systems

P 365. r r r )...(1 365

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

APPLICATIONS OF SEMIGENERALIZED -CLOSED SETS

UNIT10 PLANE OF REGRESSION

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

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?

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

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.

8 Baire Category Theorem and Uniform Boundedness

Some Approximate Analytical Steady-State Solutions for Cylindrical Fin

19 The Born-Oppenheimer Approximation

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

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

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

MHD Oscillatory Flow in a Porous Plate

Rigid Bodies: Equivalent Systems of Forces

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

Correspondence Analysis & Related Methods

An Approach to Inverse Fuzzy Arithmetic

The Greatest Deviation Correlation Coefficient and its Geometrical Interpretation

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

3. A Review of Some Existing AW (BT, CT) Algorithms

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

Mechanics Physics 151

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

Khintchine-Type Inequalities and Their Applications in Optimization

On the Distribution of the Product and Ratio of Independent Central and Doubly Non-central Generalized Gamma Ratio random variables

Part V: Velocity and Acceleration Analysis of Mechanisms

EVALUATION OF SPRINGBACK PREDICTION CAPABILITY USING UNIFORM PURE BENDING. A Thesis by. Kunal Indravadan Patel

Groupoid and Topological Quotient Group

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

Chapter Fifiteen. Surfaces Revisited

Asymptotic Waves for a Non Linear System

Exact Simplification of Support Vector Solutions

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

A Brief Guide to Recognizing and Coping With Failures of the Classical Regression Assumptions

9/12/2013. Microelectronics Circuit Analysis and Design. Modes of Operation. Cross Section of Integrated Circuit npn Transistor

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

Experimental study on parameter choices in norm-r support vector regression machines with noisy input

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

PARAMETER ESTIMATION FOR TWO WEIBULL POPULATIONS UNDER JOINT TYPE II CENSORED SCHEME

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

Chapter 8. Linear Momentum, Impulse, and Collisions

On Maneuvering Target Tracking with Online Observed Colored Glint Noise Parameter Estimation

Observer Design for Takagi-Sugeno Descriptor System with Lipschitz Constraints

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

Physics 1501 Lecture 19

A NOTE ON ELASTICITY ESTIMATION OF CENSORED DEMAND

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

Rotational Kinematics. Rigid Object about a Fixed Axis Western HS AP Physics 1

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

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

A Hybrid Variational Iteration Method for Blasius Equation

A Micro-Doppler Modulation of Spin Projectile on CW Radar

A Queuing Model for an Automated Workstation Receiving Jobs from an Automated Workstation

New Condition of Stabilization of Uncertain Continuous Takagi-Sugeno Fuzzy System based on Fuzzy Lyapunov Function

Efficiency of the principal component Liu-type estimator in logistic

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

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.

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

Bayesian Assessment of Availabilities and Unavailabilities of Multistate Monotone Systems

Review of Vector Algebra and Vector Calculus Operations

CSU ATS601 Fall Other reading: Vallis 2.1, 2.2; Marshall and Plumb Ch. 6; Holton Ch. 2; Schubert Ch r or v i = v r + r (3.

COST EVALUATION OF A TWO-ECHELON INVENTORY SYSTEM WITH LOST SALES AND NON-IDENTICAL RETAILERS

Density Functional Theory I

Lie Subalgebras and Invariant Solutions to the Equation of Fluid Flows in Toroidal Field. Lang Xia

q-bernstein polynomials and Bézier curves

PHY126 Summer Session I, 2008

Journal of Physics & Astronomy

Physics 11b Lecture #2. Electric Field Electric Flux Gauss s Law

4 Recursive Linear Predictor

A Characterization Theorem in Magneto-Thermal Convection in Rivlin-Ericksen Viscoelastic Fluid

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

Event Shape Update. T. Doyle S. Hanlon I. Skillicorn. A. Everett A. Savin. Event Shapes, A. Everett, U. Wisconsin ZEUS Meeting, October 15,

Chapter 12 Equilibrium and Elasticity

Rotating Disk Electrode -a hydrodynamic method

A Method of Reliability Target Setting for Electric Power Distribution Systems Using Data Envelopment Analysis

Vibration Input Identification using Dynamic Strain Measurement

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

Hybrid lattice Boltzmann finite-difference simulation of axisymmetric swirling and rotating flows

Dynamics of Rigid Bodies

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

EFFICIENT COMPUTATION OF THE GENERALIZED INERTIAL TENSOR OF ROBOTS BY USING THE GIBBS- APPELL EQUATIONS

Pattern Analyses (EOF Analysis) Introduction Definition of EOFs Estimation of EOFs Inference Rotated EOFs

Transcription:

Appled Mathematcs 00 43-438 do:0.436/am.00.5057 Publshed Onlne Novembe 00 (http://www.scrp.og/jounal/am) Analytcal and Numecal Solutons fo a Rotatng Annula Ds of Vaable Thcness Abstact Ashaf M. Zenou Daoud S. Mashat Depatment of Mathematcs Faculty of Scence Kng AbdulAzz Unvesty Jeddah Saud Aaba Depatment of Mathematcs Faculty of Scence Kafelsheh Unvesty Kaf l-sheh gypt -mal: zenou@sc.fs.edu.eg Receved July 9 00; evsed Octobe 9 00; accepted Octobe 00 In ths pape the analytcal and numecal solutons fo otatng vaable-thcness sold ds and numecal soluton fo otatng vaable-thcness annula ds ae pesented. The oute edge of the sold ds and the nne and oute edges of the annula ds ae consdeed to have clamped bounday condtons. Two dffeent cases fo the adally vayng thcness of the sold and annula dss ae gven. The numecal soluton as well as the analytcal soluton s avalable fo the fst case of the sold ds whle the analytcal soluton s not avalable fo the second case of the annula ds. Both analytcal and numecal esults fo dsplacement and stesses wll be nvestgated fo the fst case of adally vayng thcness. The accuacy of the pesent numecal soluton s dscussed and ts ablty of use fo the second case of adally vayng thcness s nvestgated. Fnally the dstbutons of dsplacement and stesses wll be pesented and the appopate compasons and dscussons ae made at the same angula velocty. Keywods: Rotatng Annula Ds Sold Ds Fnte Dffeence Method. Intoducton The poblems of otatng sold and annula dss have been pefomed unde vaous nteestng assumptons and the topc can be easly found n most of the standad elastcty boos [-5]. Most of the eseach wos ae concentated on the analytcal solutons of otatng dss wth smple coss-secton geometes of unfom thcness and especally vaable thcness [6-4]. The analytcal elastcty solutons of such otatng dss ae avalable n many boos of elastcty. As many otatng components n use have complex coss-sectonal geometes they cannot be dealt wth usng the exstng analytcal methods. Numecal methods such as the fnte element method [5] the bounday element method [6] and Runge-Kutta s algothm [7] can be appled to cope wth these otatng components. In ths pape we wll pesent the analytcal soluton fo the otatng sold ds wth abtay coss-secton of contnuously vaable thcness. In the followng a unfed govenng equaton wll be fst deved fom the basc equatons of the otatng dss and the poposed stess-stan elatonshp. Next fnte dffeence method (FDM) s ntoduced to solve the govenng equaton. A compason between both analytcal and numecal solutons s made. The accuacy of the numecal soluton s used to fnd the dsplacement and stesses of otatng vaable-thcness annula ds whose analytcal soluton s not avalable. Fnally a numbe of numecal examples ae gven to demonstate the valdty of the poposed method.. Basc quatons As the effect of thcness vaaton of otatng dss can be taen nto account n the equaton of moton the theoy of the vaable-thcness dss can gve good esults as that of the unfom-thcness dss as long as they meet the assumpton of plane stess. Afte consdeng ths effect the equaton of moton of otatng dss wth vaable thcness can be wtten as whee d d and h h h 0 () ae the adal and ccumfeental Copyght 00 ScRes.

43 A. M. ZNKOUR T AL. stesses s the adal coodnate s the densty of the otatng ds s the constant angula velocty and h s the thcness whch s functon of the adal coodnate. The elatons between the adal dsplacement u and the stans ae espectve of the thcness of the otatng ds. They can be wtten as du u () d whee and ae the adal and ccumfeental stans espectvely. Fo the elastc defomaton the consttutve equatons fo the otatng ds can be descbed wth Hooe s law. (3) Usng () nto quaton (3) one can obtan the consttutve equatons fo and as: du u d u du. d Let us consde a symmetc thn ds wth espect to the md-plane ts pofle vayng n the adal decton accodng to the fomula see Fgue : h h 0 n (5) b whee h 0 s the thcness at the axs of the ds n and ae geometc paametes h () (0n 0) and b s the oute adus of the sold ds. A unfom-thcness ds s obtaned by settng n = 0 and a lnealy deceasng thcness s obtaned by settng =. Futhemoe f < the pofle s concave and f > t s convex. The thcness of the ds s assumed to be suffcently small compaed to ts damete so that 3. Fomulaton and lastc Soluton fo Sold Ds The substtuton of quatons (4) and (5) nto quaton () poduces the followng confluent hypegeometc dffeental equaton fo the adal dsplacement u (): d u d ( ) u n d b d n b (6) n 3 b u 0. n b (4) (a) (c) Fgue. Vaable-thcness sold ds pofles fo (a) =.5 and n = 0.8; = 0.5 and n = 0.8 and (c) =.5 and n = 0.5. R b b UR u b. (7) Then quaton (6) may be wtten n the followng Copyght 00 ScRes.

A. M. ZNKOUR T AL. 433 smple fom d U R dr nr n R 3 U R nr n R R du d R 0. (8) The geneal soluton of the above equaton can be wtten as U R C F R C F R P R (9) whee C and C ae abtay constants and F and F ae gven by: F R RH nr (0) F R H nr () R n whch whee 4 () 4. The functons H ([ ] [ ] z) ae the genealzed hype-geometc functons q q q H z z q0 q! (3) q 0 z whee ( ) q s the Pochhamme symbol gven by... q q (4) q n whch epesents Gamma functon. The tem PR ( ) n quaton (9) s the patcula soluton of quaton (8) whch can be wtten as RFR RFR (5) PR4 FR drf R d R n R FF 4 F Fn R F 3 (6) n whch F3 R RH nr (7) F4 R H 3 nr (8) R The substtuton of quaton (9) nto quaton (4) wth the ad of the dmensonless foms gven n quaton (7) gves the adal and ccumfeental stesses n the foms of (9) and (0): nr nr dp P RC F F3C F F4 R R dr R nr nr d P P R C F F3C F F4. R R dr R (9) (0) 4. Analytcal Soluton fo the Rotatng Sold Ds The analytcal elastc soluton fo the sold ds wth vaable-thcness s completed by the applcaton of the bounday condtons. Snce the adal dsplacement should be vanshed and the stesses should be fnte at the cente of the ds then the constant C vanshes. The adal dsplacement s vanshed at the oute edge of the ds = b o (R = ) hence P C. () F So one can easly obtan the soluton fo the pesent otatng vaable-thcness sold ds by the substtuton of quaton () nto quatons (9) (9) and (0). 5. Fnte Dffeence Algothm fo Sold Ds The esoluton of the elastc poblem of otatng ds wth vaable thcness s to solve a second-ode dffeental equaton quaton (8) unde the gven bounday condtons. Ths equaton can be wtten n the followng geneal fom: U p R U q R U s R 0 R U 0 U 0 () whee the pme (') denotes dffeentaton wth espect to R and Copyght 00 ScRes.

434 A. M. ZNKOUR T AL. n R pr nr R n R qr nr R s R R.. (3) It s clea that the above poblem has a unque soluton because pr ( ) qr ( ) and s( R ) ae contnuous on the gven nteval ]0] and qr ( ) 0 on ]0]. The lnea second-ode bounday value poblem gven n quaton () eques that dffeence-quotent appoxmatons be used fo appoxmatng U and U. Fst we select an ntege N 0 and dvded the nteval ]0] nto ( N ) equal subntevals whose end ponts ae the mesh ponts R R fo 0... N whee R /( N ). At the nteo mesh ponts R... N the dffeental equaton to the appoxmated s U R p R U R q R U R s R (4) If we apply the centeed dffeence appoxmatons of U ( R ) and U ( R ) to quaton (4) we ave at the system: R p R ( ) U R q R U (5) R p R U R s R fo each... N. The N equatons togethe wth the bounday condtons U0 0 (6) U 0 N Ae suffcent to detemne the unnowns U 0... N. The esultng system of quatons (5) s expesses n the t-dagonal N N -matx fom: whee AU B (7) A R q R... N R A pr... N R A pr 3... N A j Aj 0... N j 3 4... N j B R s R... N. (8) The soluton of the fnte dffeence dscetzaton of the two-pont lnea bounday value poblem can theefoe be found easly even fo vey small mesh szes. 6. Numecal xamples and Dscusson fo Sold Ds Some numecal examples fo the otatng vaablethcness sold dss wll be gven accodng the analytcal and numecal solutons ( 0.3). Accodng to quaton (7) the followng dmensonless esponse chaactestcs u U R u b (9). R R detemned as pe the analytcal soluton ae compaed wth those obtaned by the numecal FDM soluton. The esults of the pesent nvestgatons fo the adal dsplacement u ae epoted n Table. Fo ths example N = 9 9 39 and 79 so R has the coespondng values 0. 0.05 0.05 and 0.05 espectvely. If we use the Rchadson extapolaton method wth R 0. 0.05 0.05 and 0.05 we obtan esults lsted n Table. The fst extapolaton s 4U R 0.05 U R 0. xt ; (30) 3 the second extapolaton s 4U R 0.05U R 0.05 xt ; (3) 3 the thd extapolaton s 4U R 0.05U R 0.05 xt 3 ; (3) 3 the foth extapolaton s 6xt xt xt 4 ; (33) 5 and the fnal extapolaton s 6 xt3 xt xt 5. (34) 5 All of the esults of xt4 and xt5 ae coect to the decmal places lsted. In fact f suffcent dgts ae mantaned these appoxmatons gve esults that agee wth the exact soluton wth maxmum eo of.9 0 8 and 4.8 0 8 espectvely at the mesh ponts. Copyght 00 ScRes.

A. M. ZNKOUR T AL. 435 The dstbuton of the adal dsplacement adal and ccumfeental stesses ae pesented n Fgue. The numecal FDM soluton s compaed wth the analytcal soluton fo the otatng vaable-thcness sold ds wth = 3 and vaous values of n. It s clea that the FDM gves dsplacement and stesses wth excellent accuacy wth the exact analytcal soluton. It can be seen fom Fgue that the FDM can descbe the dsplacement and stesses though-the-thcness vey well enough. 7. Fomulaton and Numecal Soluton fo Annula Ds Hee we consde a thn annula ds vaes contnuously n the fom of a fom of a geneal paabolc functon (see Fgue 3): a h h 0 n (35) b a (a) (a) (c) Fgue. Dmensonless adal dsplacement u adal stess σ and ccumfeental stess σ fo the vaable-thcness sold ds ( = 3): (a) n = 0. n = 0.5 and (c) n = 0.8. (c) Fgue 3. Vaable-thcness annula ds pofles fo (a) =.5 and n = 0.8 = 0.5 and n = 0.8 and (c) =.5 and n = 0.5. Copyght 00 ScRes.

436 A. M. ZNKOUR T AL. whee a s the nne adus of the annula ds. The equlbum equaton coespondng to quaton (35) may be easly gven but ts analytcal soluton s not. The substtuton of quatons (4) and (35) nto quaton () wth the help of the dmensonless foms gven n quaton (7) poduces the followng confluent hypegeometc dffeental equaton fo the adal dsplacement U( R ) of the annula ds: Fgue 5. Dmensonless adal dsplacement u of the vaable-thcness annula ds fo d eent values of and n. (a) Fgue 4. Dmensonless adal dsplacement u adal stess σ and ccumfeental stess σ fo the vaable-thcness annula ds ( = 3): (a) n = 0. n = 0.5 and (c) n = 0.8. whee d U nrr du R R d R R A nr d R 4 n RR R U 0 RAnR R A (36) a R A R A. (37) b A Mang analogous steps as gven fo the sold ds quaton (36) can be wtten n the followng geneal fom: U pruqru sr (38) A R U A U 0 whee the pme (') denotes dffeentaton wth espect to R and nrr pr R R A nr n RR qr R R AnR R sr. R A (39) So FDM gves easly the adal dsplacement of the otatng vaable-thcness annula ds. Usng the cuve fttng and least squae method one can obtan the adal Copyght 00 ScRes.

A. M. ZNKOUR T AL. 437 and ccumfeental stesses. Tang n = 0. 0.5 and 0.8 and = 0.5.5.5 and 3 n the vaable thcness functon gven n quaton (35) a otatng annula ds wth such vaable thcness s studed. The nne and oute ad of the ds ae taen to be a = 0. b (R = A = 0.) and b (R = ) and the esults ae gven n tems of the otatng angula velocty. The esults calculated wth the FDM fo dsplacement and stesses of a otatng vaable-thcness annula ds ae gven n Fgues 3-7. 8. Conclusons Ths pape pesents a unfed numecal method fo the elastc calculaton of otatng dss wth a geneal abtay confguaton. The govenng equaton was deved fom the equlbum equaton and the stess- stan elatonshp. The analytcal soluton was gven fo the otatng vaable-thcness sold ds. The calculaton of the otatng sold and annula dss was tuned nto fndng the soluton of a second-ode dffeental equaton unde the gven condtons at two bounday fxed ponts. Fnte dffeence method algothm was ntoduced to solve the govenng equaton fo both sold and annula dss and a numbe of numecal examples wee studed. The esults fom the analytcal and FDM solutons wee compaed. The poposed FDM appoach gves vey ageeable esults to the analytcal soluton. 9. Acnowledgements The nvestgatos would le to expess the appecaton to the Deanshp of Scentfc Reseach at Kng AbdulAzz Unvesty fo ts fnancal suppot of ths. 68/48 /ع No. study Gant 0. Refeences Fgue 6. Dmensonless adal stess σ n the vaablethcness annula ds fo dffeent values of and n. Fgue 7. Dmensonless ccumfeental stess σ n the vaable-thcness annula ds fo d eent values of and n. [] S. P. Tmosheno and J. N. Goode Theoy of lastcty McGaw-Hll New Yo 970. [] S. C. Ugal and S. K. Fenste Advanced Stength and Appled lastcty lseve New Yo 987. [3] U. Game Tesca s Yeld Condton and the Rotatng Ds ASM Jounal of Appled Mechancs Vol. 50 No. 983 pp. 676-678. [4] U. Game lastc-plastc Defomaton of the Rotatng Sold Ds Ingeneu-Achv Vol. 45 No. 4 984 pp. 345-354. [5] A. M. Zenou Analytcal Solutons fo Rotatng xponentally-gaded Annula Dss wth Vaous Bounday Condtons Intenatonal Jounal of Stuctual Stablty and Dynamcs Vol. 5 No. 4 005 pp. 557-577. [6] U. Güven lastc-plastc Stesses n a Rotatng Annula Ds of Vaable Thcness and Vaable Densty Intenatonal Jounal of Mechancal Scences Vol. 34 No. 99 pp. 33-38. [7] U. Güven On the Stess n the lastc-plastc Annula Ds of Vaable Thcness unde xtenal Pessue Intenatonal Jounal of Solds and Stuctues Vol. 30 No. 5 993 pp. 65-658. [8] U. Güven Stess Dstbuton n a Lnea Hadenng Annula Ds of Vaable Thcness Subjected to xtenal Pessue Intenatonal Jounal of Mechancal Scences Vol. 40 No. 6 998 pp. 589-60. [9] U. Güven lastc-plastc Stesses Dstbuton n a Rotatng Hypebolc Ds wth Rgd Incluson Intenatonal Jounal of Mechancal Scences Vol. 40 No. 998 pp. 97-09. Copyght 00 ScRes.

438 A. M. ZNKOUR T AL. [0] A. N. aslan Inelastc Defomaton of Rotatng Vaable Thcness Sold Dss by Tesca and Von Mses Ctea Intenatonal Jounal fo Computatonal Methods n ngneeng Scence Vol. 3 000 pp. 89-0. [] A. N. aslan and Y. Ocan On the Rotatng lastc-plastc Sold Dss of Vaable Thcness Havng Concave Pofles Intenatonal Jounal of Mechancal Scences Vol. 44 No. 7 00 pp. 445-466. [] A. N. aslan Stess Dstbutons n lastc-plastc Rotatng Dss wth llptcal Thcness Pofles Usng Tesca and von Mses Ctea Zuch Altenatve Asset Management Vol. 85 005 pp. 5-66. [3] A. M. Zenou and M. N. M. Allam On the Rotatng Fbe-Renfoced Vscoelastc Composte Sold and Annula Dss of Vaable Thcness Intenatonal Jounal fo Computatonal Methods n ngneeng Scence and Mechancs Vol. 7 No. 006 pp. -3. [4] A. M. Zenou Themoelastc Solutons fo Annula Dss wth Abtay Vaable Thcness Stuctual ngneeng and Mechancs Vol. 4 No. 5 006 pp. 55-58. [5] O. C. Zenewcz The Fnte lement Method n ngneeng Scence McGaw-Hll London 97. [6] P. K. Banejee and R. Buttefeld Bounday lement Methods n ngneeng Scence McGaw-Hll New Yo 98. [7] L. H. You Y. Y. Tang J. J. Zhang and C. Y. Zheng Numecal Analyss of lastc-plastc Rotatng Dss wth Abtay Vaable Thcness and Densty Intenatonal Jounal of Solds and Stuctues Vol. 37 No. 5 000 pp. 7809-780. Copyght 00 ScRes.