ANALYSIS OF THERMOELASTIC STRESSES IN LAYERED PLATES

Similar documents
Analysis of composites with multiple rigid-line reinforcements by the BEM

DYNAMIC ANALYSIS OF BEAM-LIKE STRUCTURES SUBJECT TO MOVING LOADS

Atomic Physics 4. Name: Date: 1. The de Broglie wavelength associated with a car moving with a speed of 20 m s 1 is of the order of. A m.

Principle Of Superposition

DETERMINATION OF MECHANICAL PROPERTIES OF A NON- UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS.

METHOD OF FUNDAMENTAL SOLUTIONS FOR HELMHOLTZ EIGENVALUE PROBLEMS IN ELLIPTICAL DOMAINS

A NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS

Stability analysis of numerical methods for stochastic systems with additive noise

UNIVERSITY OF CALIFORNIA - SANTA CRUZ DEPARTMENT OF PHYSICS PHYS 116C. Problem Set 4. Benjamin Stahl. November 6, 2014

The axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c.

H05 A 3-D INVERSE PROBLEM IN ESTIMATING THE TIME-DEPENDENT HEAT TRANSFER COEFFICIENTS FOR PLATE FINS

Finite Difference Method for the Estimation of a Heat Source Dependent on Time Variable ABSTRACT

The Advection-Diffusion equation!

Damped Vibration of a Non-prismatic Beam with a Rotational Spring

Numerical Simulation of Thermomechanical Problems in Applied Mechanics: Application to Solidification Problem

EVALUATION OF GLASS FIBER/EPOXY INTERFACIAL STRENGTH BY THE CRUCIFORM SPECIMEN METHOD

DEGENERACY AND ALL THAT

On the convergence, consistence and stability of a standard finite difference scheme

PRELIM PROBLEM SOLUTIONS

University of Washington Department of Chemistry Chemistry 453 Winter Quarter 2015

Discrete Orthogonal Moment Features Using Chebyshev Polynomials

Finally, we show how to determine the moments of an impulse response based on the example of the dispersion model.

Comparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series

The Method of Least Squares. To understand least squares fitting of data.

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,

THE ENERGY BALANCE ERROR FOR CIRCUIT TRANSIENT ANALYSIS

Numerical Conformal Mapping via a Fredholm Integral Equation using Fourier Method ABSTRACT INTRODUCTION

Estimating the Population Mean using Stratified Double Ranked Set Sample

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j

Streamfunction-Vorticity Formulation

Scientific Research of the Institute of Mathematics and Computer Science

ANALYSIS OF DAMPING EFFECT ON BEAM VIBRATION

Similarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall

FAILURE CRITERIA: MOHR S CIRCLE AND PRINCIPAL STRESSES

62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +

Analytical solutions for multi-wave transfer matrices in layered structures

ON LOCAL LINEAR ESTIMATION IN NONPARAMETRIC ERRORS-IN-VARIABLES MODELS 1

x x x 2x x N ( ) p NUMERICAL METHODS UNIT-I-SOLUTION OF EQUATIONS AND EIGENVALUE PROBLEMS By Newton-Raphson formula

THE NUMERICAL SOLUTION OF THE NEWTONIAN FLUIDS FLOW DUE TO A STRETCHING CYLINDER BY SOR ITERATIVE PROCEDURE ABSTRACT

MATHEMATICAL MODELLING OF ARCH FORMATION IN GRANULAR MATERIALS

ME 501A Seminar in Engineering Analysis Page 1

Inelastic spherical caps with defects

Application of a Two-Step Third-Derivative Block Method for Starting Numerov Method

Numerical Method for Blasius Equation on an infinite Interval

Numerical Study on MHD Flow And Heat Transfer With The Effect Of Microrotational Parameter In The Porous Medium

NUMERICAL SOLUTIONS OF THE FRACTIONAL KdV-BURGERS-KURAMOTO EQUATION

Regression with an Evaporating Logarithmic Trend

A CONFINEMENT MODEL OF HIGH STRENGTH CONCRETE

1 1 2 = show that: over variables x and y. [2 marks] Write down necessary conditions involving first and second-order partial derivatives for ( x0, y

2C09 Design for seismic and climate changes

Recurrence Relations

On Exact Finite-Difference Scheme for Numerical Solution of Initial Value Problems in Ordinary Differential Equations.

Celestin Chameni Nembua University of YaoundéII, Cameroon. Abstract

CO-LOCATED DIFFUSE APPROXIMATION METHOD FOR TWO DIMENSIONAL INCOMPRESSIBLE CHANNEL FLOWS

Homogenization Model for the Folded Core Sandwich Plates under the Transverse Loading

Lecture 4 Conformal Mapping and Green s Theorem. 1. Let s try to solve the following problem by separation of variables

Poincaré Problem for Nonlinear Elliptic Equations of Second Order in Unbounded Domains

PHYS-3301 Lecture 9. CHAPTER 5 Wave Properties of Matter and Quantum Mechanics I. 5.3: Electron Scattering. Bohr s Quantization Condition

MATH 6101 Fall 2008 Newton and Differential Equations

Chapter 9: Numerical Differentiation

Measurement uncertainty of the sound absorption

Nonparametric regression: minimax upper and lower bounds

ON POINTWISE BINOMIAL APPROXIMATION

Dynamical stability and parametrical vibrations of the laminated plates with complex shape

11 Correlation and Regression

APPENDIX A EARLY MODELS OF OXIDE CMP

Partial Differential Equations

A Pseudo Spline Methods for Solving an Initial Value Problem of Ordinary Differential Equation

THE SYSTEMATIC AND THE RANDOM. ERRORS - DUE TO ELEMENT TOLERANCES OF ELECTRICAL NETWORKS

A New Hybrid in the Nonlinear Part of Adomian Decomposition Method for Initial Value Problem of Ordinary Differential Equation

Some Basic Diophantine Equations

Chapter 1. Complex Numbers. Dr. Pulak Sahoo

Exact Solutions for a Class of Nonlinear Singular Two-Point Boundary Value Problems: The Decomposition Method

The standard deviation of the mean

ALLOCATING SAMPLE TO STRATA PROPORTIONAL TO AGGREGATE MEASURE OF SIZE WITH BOTH UPPER AND LOWER BOUNDS ON THE NUMBER OF UNITS IN EACH STRATUM

GROUND MOTION OF NON-CIRCULAR ALLUVIAL VALLEY FOR INCIDENT PLANE SH-WAVE. Hui QI, Yong SHI, Jingfu NAN

Generating Functions for Laguerre Type Polynomials. Group Theoretic method

On forward improvement iteration for stopping problems

A Block Cipher Using Linear Congruences

A PROOF OF THE TWIN PRIME CONJECTURE AND OTHER POSSIBLE APPLICATIONS

7.1 Convergence of sequences of random variables

ON THE EIGENFUNCTION EXPANSION METHOD FOR THE CALCULATION OF GREEN S FUNCTIONS

Computation of Hahn Moments for Large Size Images

ON THE EIGENFUNCTION EXPANSION METHOD FOR THE CALCULATION OF GREEN S FUNCTIONS

Similarity between quantum mechanics and thermodynamics: Entropy, temperature, and Carnot cycle

DECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS. Park Road, Islamabad, Pakistan

Mathematical Modeling of Optimum 3 Step Stress Accelerated Life Testing for Generalized Pareto Distribution

Lecture 9: Regression: Regressogram and Kernel Regression

Properties and Hypothesis Testing

Experiment (2) Upsetting Test of Forgeability of Metal 1. Objective: 2. Background: 3. Theory

1 6 = 1 6 = + Factorials and Euler s Gamma function

WHAT IS THE PROBABILITY FUNCTION FOR LARGE TSUNAMI WAVES? ABSTRACT

18.01 Calculus Jason Starr Fall 2005

Orthogonal Gaussian Filters for Signal Processing

Asymptotic distribution of products of sums of independent random variables

a b c d e f g h Supplementary Information

Session 5. (1) Principal component analysis and Karhunen-Loève transformation

[412] A TEST FOR HOMOGENEITY OF THE MARGINAL DISTRIBUTIONS IN A TWO-WAY CLASSIFICATION

Analysis of MOS Capacitor Loaded Annular Ring MICROSTRIP Antenna

Chapter 4. Fourier Series

Transcription:

6t Iteratioal DAAAM Baltic Coferece "INDUSTRIAL ENGINEERING" 4-6 April 8, Talli, Estoia ANALYSIS OF THERMOELASTIC STRESSES IN LAYERED PLATES Kõo, J. & Valgur, J. Abstract: Tis paper deals wit termoelastic effects i deformatio of plates wit arbitrarily cagig elastic parameters ad temperature troug tickess. Usig te semi-iverse metod, a simple aalytical solutio is obtaied for te termoelastic problem of a oomogeeous plate wit a arbitrary cotour. Plates wit free, slidig ad fixed edges are cosidered. Te results are applied to te layered plate. As a example of applicatio, termoelastic stresses ad deformatios are determied for a copper plate wit a steel coatig. Te edge effects for tis plate are examied by te fiite elemet metod. Key words: semi-iverse metod; termoelastic stresses; layered plates; edge effects; fiite elemet metod.. INTRODUCTION I moder egieerig desig, oomogeeous layered structures are widely used due to teir superiority i stiffess-toweigt ad stregt-to-weigt ratios, as well as logevity. Te study of mecaical, especially termoelastic, caracteristics of layered beams ad plates as bee of icreasig iterest for egieers. Termo-mecaical aalysis of layered beams (arrow strips) is preseted i papers [, ]. I te moograp [ 3 ] te termoelastic problem is solved for a omogeeous free plate wit temperature variatio troug tickess. I te paper [ 4 ] tis solutio is exteded to te plate wit variatio of termoelastic parameters troug tickess. I preset paper te solutio uder cosideratio will be applied to termoelastic aalysis of layered plates. I coectio wit tis, te plate wit free, slidig, ad fixed edges will be studied. As a example of suc applicatio, termoelastic stresses ad deformatios will be determied for a quadratic copper plate wit a steel coatig. Te edge effects for tis plate will be examied usig te fiite elemet program ANSYS.. PLATE WITH ARBITRARY VARIATION OF THERMOELASTIC PARAMETERS THROUGH THICKNESS.. Plate wit free edges Cosider a plate (Fig. ) of a arbitrary form ad of costat tickess. Te plate is completely free of te surface load. Te rectagular coordiates x, y ad are used, were te lower surface of te plate is take as te referece surface (x, y) ad te coordiate is directed upwards. Fig.. Free plate wit a arbitrary sape

Te plate is fixed at te origi of te coordiates, were deflectio ad its derivatios are ero, i.e. u u = = =. x y Te temperature T, te coefficiet of liear termal expasio α, te modulus of elasticity (Youg s modulus) E ad te Poisso s ratio μ vary troug tickess oly, = E = E ad i.e. T T( ), α=α( ), μ = μ ( ). By aalogy wit te omogeeous plate [ 3 ], it is reasoable to assume tat te stress compoets will be of te followig form: x = y = ( ). () = τxy = τ y = τx = Tese stresses satisfy idetically te followig equilibrium equatios of te teory of elasticity [ 3, 5 ]: τ x yx τx + + = (x, y, ). () x y Here ad below te symbol (x, y, ) deotes permutatio by meas of wic we ca write two more equatios correspodig to te oter two axes, cagig x for y, y for, ad for x. From te termoelastic stress-strai relatios [ 3, 5 ] ε = ( + ) μ + α T (3) x E x y (x, y, ) we ave εx = εy = ε( ) = + αt E, (4) μ ε = + αt E were E E = (5) μ is te biaxial elastic modulus. For strais (4) te compatibility equatios [ 3, 5 ] ε γ y x x y are reduced to ε x y xy + = (x, y, ) (6) + T =. E α (7) Terefore = E α T + a+ b, (8) were te costats a ad b are to be determied from te boudary coditios of ero resultat force ad momet o te edges of te plate. Tese coditios are reduced to: d = d =. (9) Substitutio of expressio (8) yields te two equatios for te ukow costats a ad b. From tese equatios we fid: CNT BMT a = C BD, () CMT DNT b = C BD were te quatities N T, M T, B, C ad D are defied as follows: NT = E αtd ; () M = T E αtd B= E d C = E d. () D= E d Te displacemets u x, u y ad u i te x, y ad directios, respectively, are obtaied by usig equatios [ 3, 5 ] ux ε x = (x, y, ) (3) x ad (4) as well as (8). Te solutio is te foud to be: ux = x( a+ b) ; (4) uy = y( a+ b)

a u = ( x + y ) + (5) μ + μ + ( a + b) + α T. d μ μ For completeess of te solutio, we determie te deformatio caracteristics of te referece surface ( = ). Usig equatios (4) ad (3) we ave for strai o te referece surface: ε ε =b. (6) I te case of small deflectios we ca use for curvatures te approximatios: u u æx = ; æ y =. x y Te, usig equatio (5) we fid te curvatures of te referece surface æx = æy æ = a. (7) Tus, te plate subjected to a uiform temperature cage is bet to a sperical surface... Plate wit slidig edges For a plate wit slidig edges (Fig. ), we ca suppose tat u =, = = x y for every poit of te referece surface. Fig.. Plate wit slidig edges Te a = ad from equatios () it follows tat N b = T (8) B ad expressio (8) for stresses is reduced to = E α T + b. (9).3. Plate wit fixed edges For a plate wit fixed edges (Fig. 3), we ca start from te assumptios tat ux = uy = u =, = = x y o te referece surface. Fig. 3. Plate wit fixed edges Te a = b = ad equatio (8) for stresses is reduced to = E α T. () 3. LAYERED PLATE SUBJECTED TO A UNIFORM TEMPERATURE CHANGE 3.. Plate wit free edges Cosider a layered plate (Fig. 4) cosistig of a arbitrary umber of layers wit differet costat tickesses i. Fig. 4. Layered plate Te locatio of a arbitrary layer i of te plate is specified by te coordiate i, wic is te distace from te bottom plae of te plate to te top plae of te it layer. Assume tat te coefficiet of termal expasio α i, te modulus of elasticity E i ad te Poisso s ratio μ i do ot cage troug te tickess of te it layer. Te, assumig tat te plate is subjected to a uiform temperature cage ΔT, ad usig piecewise itegratio, we ca write te parameters N T, M T, B, C ad D i equatios () ad () as NT =ΔT Eiαi( i i ) =ΔT Eiαii i= i= ;() ΔT M = ( ) T Eiα i i i i=

B E E = i ( i i ) = i i i= i= = i ( i i ) i= 3 3 = i ( i i ) 3 i= C E D E. () For te costats () we ave a = aδt ; b= bδt, were CNT BMT a = C BD. (3) CMT DNT b = C BD Here NT = Eiαi( i i ) = Eiαii i= i=. (4) M = ( ) T Eiαi i i i= I accordace wit equatio (8) stresses i te it layer are = E ΔT α + a+ b. (5) i i i Deformatio parameters are ε = bδt. (6) æ = aδt 3.. Plate wit slidig edges I tis case a =, ad N b = T B. (7) Stresses i te it layer are = E ΔT α + b. (8) i i i Strai o te referece surface ε = bδt (9). 3.3. Plate wit fixed edges For tis case a = b =, ad stresses i te it layer are = E α ΔT. (3) i i i 4. APPLICATIONS: THERMAL STRESSES DUE TO UNIFORM TEMPERATURE CHANGE IN A COPPER PLATE WITH A GALVANIC STEEL COATING Usually, te coefficiets of termal expasio of te substrate ad coatig are differet ad te temperature of te coatig process differs from room temperature (+ C). Terefore, termal stresses are geerated i coated parts. I experimetal aalysis of residual stresses it may tur out tat coatig temperature differs from te temperature at wic te deformatio parameters of te substrate are measured. I suc cases measuremet results ave to be corrected takig ito accout termal stresses. As a umerical example, termal stresses ad displacemets for a quadratic copper plate 3 3 mm wit galvaic steel coatig deposited at 95 C are calculated. Fig. 5. Copper plate wit a steel coatig I additio, te followig data is used: = =. mm; E = GPa; μ =.34; α = = 7.5-6 / C for te substrate ad = =.6 mm; E = GPa; μ =.8; α = = 4.3-6 / C for te coatig.

4.. Plate wit free edges For uiform temperature cage ΔT = = 95 = 75 C, te values defied by equatios (5), () ad (4) are: E = 66.7 GPa; E = 8.6 GPa; N = 396 N/m; M =.638 N; T B = 39.66 6 N/m; C = 65.79 3 N; D = 49. N m. By meas of equatios (3) we ave a = -.98-3 /m; b = 8.635-6. For stresses i layers from equatio (5) it follows tat: 3 6 i = 75Ei ( αi +.98 8.635 ) (i =, ). (3) At typical poits of tickess stresses i MPa are: = 4.9 ( 3.5) 3 ( ) = 3.7 (3.6). (3) 3 ( ) = 8.5 ( 7.56) 3 (.6 ) =. (.3) From (6) for te strai ad curvature of te referece surface we fid 6 6 ε = 398 ( 393.6 ) 3 3 -. (33) æ = 6 (.53 ) m 4.. Plate wit slidig edges Accordig to (7), we ave b =6.58-6 ad for stresses i layers (8) we obtai 6 = 75 ( 6.58 i Ei α i ) (34) (i =, ). From ere for stresses i MPa we ave: =.5 (.6). (35) = 46.89 ( 46.77) Strai o te referece surface (9) as te value 6 6 ε = 4 ( 39.6 ). (36) 4.3. Plate wit fixed edges By meas of equatio (3) we calculate stresses i te layers i MPa: = 8.79 (8.7). (37) = 3.94 (3.8) T 4.4. FEM aalysis A similar problem was aalysed by te fiite elemet program ANSYS. Te obtaied distributio of ormal stress x = = y = i te middle regio of te plate is preseted i Fig. 6. Te values of stresses at typical poits of tickess are give i te pareteses at te give data sets (3), (35) ad (37). Te values of te deformatio parameters ε ad æ are preseted i te same maer at te data sets (33) ad (36). Compariso of te results sows tat discrepacy betwee data is less ta 5%. Some results of stress aalysis i te edge regio (.5 mm x 5 mm, y =,.6 mm) are preseted i Fig. 7. As we ca see, te ormal stresses x = y = i te edge regio are reduced. Te ormal stress, kow as peelig stress, ad te sear stress τ x are igly localied ear te edge. Tese facts oce agai cofirm te validity of te well-kow Sait-Veat s priciple. Fig. 6. Normal stresses x = y = () i te middle regio of a free copper plate wit a steel coatig

(a) (c) (b) (d) Fig. 7. Stress distributios i te edge regio: (a) ormal stress x, (b) ormal stress y, (c) ormal stress (peelig stress), (d) sear stress τ x (MPa) 5. CONCLUSIONS. Usig te semi-iverse metod, a aalytical solutio is obtaied for te termoelastic problem of a plate wit a arbitrary cotour oomogeeous troug tickess. Te results are applied for aalysis of a layered plate subjected to a uiform temperature cage.. As a applicatio, stresses due to uiform temperature cage i a copper plate wit a galvaic steel coatig are determied. 3. Usig te fiite elemet metod, te termo-mecaical state of te same plate wit te coatig is studied. For te middle regio, good agreemet is acieved betwee te data of te aalytical ad umerical metods. Te data obtaied wit te use of te fiite elemet metod for stresses i te edge regio of te plate cofirm te validity of te well-kow Sait- Veat s priciple. 6. REFERENCES. Feg, Z. C. ad Liu, H. D. Geeralied formula for curvature radius ad layer stresses caused by termal strai i semicoductor multilayer structures. J. Appl. Pys., 983, 54 (), 83 85.. Hsue, C. H. Termal stresses i elastic multilayer systems. Ti Solid Films,, 48, 8 88. 3. Boley, B. A. ad Weier, J. H. Teory of Termal Stresses. Jo Wiley ad Sos, Ic., New York, Lodo, 96. 4. Kõo, J. P. Problem of termoelasticity for a free oomogeeous plate wit temperature variatio troug tickess. I Heat Stresses i Structural Members. Naukova Dumka, Kyiv, 979, 9, 89 9 (i Russia). 5. Timoseko, S. P., Goodier, J. N. Teory of Elasticity. Tird Editio. McGraw-Hill, Sigapore, 97. 7. CORRESPONDING AUTHOR Prof. emer. Jakub Kõo Estoia Uiversity of Life Scieces Istitute of Forestry ad Rural Egieerig Kreutwaldi 5, Tartu 54, Estoia, E-mail: jakub.koo@emu.ee