A Brief Note on Quasi Static Thermal Stresses In A Thin Rectangular Plate With Internal Heat Generation

Similar documents
An inverse steady state thermal stresses in a thin clamped circular plate with internal heat generation

Solutions of Klein - Gordan equations, using Finite Fourier Sine Transform

Summary: Method of Separation of Variables

MAC-solutions of the nonexistent solutions of mathematical physics

THE INTERVAL LATTICE BOLTZMANN METHOD FOR TRANSIENT HEAT TRANSFER IN A SILICON THIN FILM

Module 9: The Method of Green s Functions

arxiv: v1 [math.gm] 30 Dec 2015

M344 - ADVANCED ENGINEERING MATHEMATICS

The Shortest Confidence Interval for the Mean of a Normal Distribution

Lecture 24: Laplace s Equation

Hamiltonian Cycle in Complete Multipartite Graphs

The Dirichlet Problem in a Two Dimensional Rectangle. Section 13.5

Studies on Nuclear Fuel Rod Thermal Performance

Differential Equations 2 Homework 5 Solutions to the Assigned Exercises

Green function and Eigenfunctions

1.1. Linear Constant Coefficient Equations. Remark: A differential equation is an equation

Thin and Thick Cylinders and Spheres

New Expansion and Infinite Series

1 Bending of a beam with a rectangular section

Creating A New Planck s Formula of Spectral Density of Black-body Radiation by Means of AF(A) Diagram

Adomian Decomposition Method with Green s. Function for Solving Twelfth-Order Boundary. Value Problems

Consequently, the temperature must be the same at each point in the cross section at x. Let:

Entropy ISSN

Math 124A October 04, 2011

INEQUALITIES FOR TWO SPECIFIC CLASSES OF FUNCTIONS USING CHEBYSHEV FUNCTIONAL. Mohammad Masjed-Jamei

Edexcel GCE Core Mathematics (C2) Required Knowledge Information Sheet. Daniel Hammocks

(PDE) u t k(u xx + u yy ) = 0 (x, y) in Ω, t > 0, (BC) u(x, y, t) = 0 (x, y) on Γ, t > 0, (IC) u(x, y, 0) = f(x, y) (x, y) in Ω.

Modification Adomian Decomposition Method for solving Seventh OrderIntegro-Differential Equations

ANALYSIS OF MECHANICAL PROPERTIES OF COMPOSITE SANDWICH PANELS WITH FILLERS

Review of Gaussian Quadrature method

MA 201: Partial Differential Equations Lecture - 12

1 E3102: a study guide and review, Version 1.0

Calculus of variations with fractional derivatives and fractional integrals

Deteriorating Inventory Model for Waiting. Time Partial Backlogging

Plate Bending Analysis by using a Modified Plate Theory

Method of Localisation and Controlled Ejection of Swarms of Likely Charged Particles

AMPERE CONGRESS AMPERE on Magnetic Resonance and Related Phenomena. Under the auspices of The GROUPEMENT AMPERE

Derivations for maximum likelihood estimation of particle size distribution using in situ video imaging

Generalizations of the Basic Functional

Harmonic Mean Derivative - Based Closed Newton Cotes Quadrature

PLEASE SCROLL DOWN FOR ARTICLE

An iterative method for solving nonlinear functional equations

Fredholm Integral Equations of the First Kind Solved by Using the Homotopy Perturbation Method

3 Mathematics of the Poisson Equation

7.1 Integral as Net Change and 7.2 Areas in the Plane Calculus

OBSERVATIONS ON TERNARY QUADRATIC EQUATION z 82 x

The Existence of the Moments of the Cauchy Distribution under a Simple Transformation of Dividing with a Constant

QUADRATURE is an old-fashioned word that refers to

Research Article Analytical Solution of the Fractional Fredholm Integrodifferential Equation Using the Fractional Residual Power Series Method

ORDER REDUCTION USING POLE CLUSTERING AND FACTOR DIVISION METHOD

HT Module 2 Paper solution. Module 2. Q6.Discuss Electrical analogy of combined heat conduction and convection in a composite wall.

Research Article On Existence and Uniqueness of Solutions of a Nonlinear Integral Equation

Lecture 4 Notes, Electromagnetic Theory I Dr. Christopher S. Baird University of Massachusetts Lowell

13.4 Work done by Constant Forces

Section 4: Integration ECO4112F 2011

Section 6.1 INTRO to LAPLACE TRANSFORMS

10 Vector Integral Calculus

AN INEQUALITY OF OSTROWSKI TYPE AND ITS APPLICATIONS FOR SIMPSON S RULE AND SPECIAL MEANS. I. Fedotov and S. S. Dragomir

A Modified ADM for Solving Systems of Linear Fredholm Integral Equations of the Second Kind

EFFECT OF RADIATION ON NATURAL CONVECTION FLOW FROM A POROUS VERTICAL PLATE IN PRESENCE OF HEAT GENERATION

Evaluation of Reliability and Availability Characteristics of 2-out of -3 Standby System under a Perfect Repair Condition

Eigen Values and Eigen Vectors of a given matrix

Torsion in Groups of Integral Triangles

Calculus of Variations

2. THE HEAT EQUATION (Joseph FOURIER ( ) in 1807; Théorie analytique de la chaleur, 1822).

International ejournals

CHM Physical Chemistry I Chapter 1 - Supplementary Material

Research Article Fejér and Hermite-Hadamard Type Inequalities for Harmonically Convex Functions

1 The fundamental theorems of calculus.

Numerical simulation of ohmic heating in idealized thin-layer electrodeposition cells

Final Exam - Review MATH Spring 2017

Partial Derivatives. Limits. For a single variable function f (x), the limit lim

u t = k 2 u x 2 (1) a n sin nπx sin 2 L e k(nπ/l) t f(x) = sin nπx f(x) sin nπx dx (6) 2 L f(x 0 ) sin nπx 0 2 L sin nπx 0 nπx

APM346H1 Differential Equations. = u x, u = u. y, and u x, y =?. = 2 u t and u xx= 2 u. x,t, where u t. x, y, z,t u zz. x, y, z,t u yy.

Thomas Whitham Sixth Form

In the next lecture... Tutorial on ideal cycles and component performance.

Homogeneous Bi-Quadratic Equation With Five Unknowns

PDE Notes. Paul Carnig. January ODE s vs PDE s 1

(uv) = u v + uv, (1) u vdx + b [uv] b a = u vdx + u v dx. (8) u vds =

Research Article Harmonic Deformation of Planar Curves

POSITIVE SOLUTIONS FOR SINGULAR THREE-POINT BOUNDARY-VALUE PROBLEMS

Time Dependent Slip MHD Flow Due to Stretching or Shrinking Vertical Sheet with Thermal Radiation

We partition C into n small arcs by forming a partition of [a, b] by picking s i as follows: a = s 0 < s 1 < < s n = b.

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1

When a force f(t) is applied to a mass in a system, we recall that Newton s law says that. f(t) = ma = m d dt v,

DISTRIBUTION OF SUB AND SUPER HARMONIC SOLUTION OF MATHIEU EQUATION WITHIN STABLE ZONES

Section 3.2 Maximum Principle and Uniqueness

Some Improvements of Hölder s Inequality on Time Scales

Wave Equation on a Two Dimensional Rectangle

Chapter 6 Techniques of Integration

THE VISCOELASTIC CONTACT BETWEEN HIGH-ORDER POLYNOMIAL SURFACES

u(x, y, t) = T(t)Φ(x, y) 0. (THE EQUATIONS FOR PRODUCT SOLUTIONS) Plugging u = T(t)Φ(x, y) in (PDE)-(BC) we see: There is a constant λ such that

ISSN: International Journal of Science, Engineering and Technology Research (IJSETR) Volume 2, Issue 4, April 2013

( ) + ( ) = ( ), [, ]

Solving the (3+1)-dimensional potential YTSF equation with Exp-function method

Thin plates with large deflection - Buckling instability and response

Problem 1. Solution: a) The coordinate of a point on the disc is given by r r cos,sin,0. The potential at P is then given by. r z 2 rcos 2 rsin 2

7.1 Integral as Net Change Calculus. What is the total distance traveled? What is the total displacement?

Variational Techniques for Sturm-Liouville Eigenvalue Problems

Theoretische Physik 2: Elektrodynamik (Prof. A.-S. Smith) Home assignment 4

Transcription:

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 on Qusi Sttic Therml Stresses In A Thin Rectngulr Plte With Internl Het Genertion C. M. Bhongde, nd M. H. Durge 1 Deprtment of Mthemtics, Shri. Shivji College, Rjur, Mhrshtr, Indi Deprtment of Mthemtics, Annd Niketn College, Wror, Mhrshtr, Indi Astrct: - The present pper dels with the determintion of qusi sttic therml stresses in thin rectngulr plte with internl het genertion. A thin rectngulr plte is considered hving zero initil temperture nd sujected to ritrry het supply t x = nd x = where s thin rectngulr plte is insulted t y = nd y =. Here we modify Kulkrni (7). The governing het conduction eqution hs een solved y the method of integrl trnsform technique. The results re otined in series form in terms of Bessel s functions. The results for temperture, displcement nd stresses hve een computed numericlly nd illustrted grphiclly. Keywords: - Qusi sttic therml stresses, thermoelstic prolem, internl het genertion, thin rectngulr plte. I. INTRODUCTION Gogulwr nd Deshmukh (4) determined the therml stresses in rectngulr plte due to prtilly distriuted het supply. Nsser M. et l. (4) solved two dimensionl prolem of thick plte with het sources in generlized thermoelsticity. Khndit nd Deshmukh (1) studied thermoelstic prolem in rectngulr plte with het genertion. Recently Ptil et l. (13) determined the therml stresses in rectngulr sl with internl het source, now here thin rectngulr plte is considered hving zero initil temperture nd nd sujected to ritrry het supply t x = nd x = where s the plte is insulted t y = nd y =. Here we modify Kulkrni (7). To otin the temperture distriution, cosine integrl trnsform nd Lplce trnsform re pplied. The results re otined in series form in terms of Bessel s functions nd the temperture chnge, displcement nd stresses hve een computed numericlly nd illustrted grphiclly. A mthemticl model hs een constructed of thin rectngulr plte with the help of numericl illustrtion y considering steel (.5% cron) rectngulr plte. No one previously studied such type of prolem. This is new contriution to the field. The direct prolem is very importnt in view of its relevnce to vrious industril mechnics sujected to heting such s the min shft of lthe, turines nd the role of rolling mill, se of furnce of oiler of therml power plnt nd gs power plnt. II. FORMULATION OF THE PROBLEM A thin rectngulr plte occupying the spce D: x, y, is considered. A thin rectngulr plte is considered hving zero initil temperture nd sujected to ritrry het supply t x = nd x = where s the plte is insulted t y = nd y =. Here the plte is ssumed sufficiently thin nd considered free from trction. Since the plte is in plne stress stte without ending. Airy stress function method is pplicle to the nlyticl development of the thermoelstic field. The eqution is given y the reltion + U = x y t E + x y T (1) where t, E nd U re liner coefficient of the therml expnsion, Young s modulus elsticity of the mteril of the plte nd Airys stress functions respectively. w w w. j e r. o r g Pge 388

Americn Journl of Engineering Reserch (AJER) 13 The displcement components u x nd u y in the X nd Y direction re represented in the integrl form nd the stress components in terms of U re given y u x = 1 U U E y x tt dx () u y = 1 U U E x y tt dy (3) ς xx = U y (4) ς yy = U x nd ς xy = U x y (6) ς xy = t y =. (7) where v is the Poisson s rtio of the mteril of the rectngulr plte. The temperture of the thin rectngulr plte t time t stisfying het conduction eqution s follows, T T T x y k α t (8) T x, y, t = t t = x, y (9) T(x, y, t) = f 1 y, t t x =, y (1) T x, y, t = f y, t t x =, y (11) T y T y = t y =, x (1) = t y =, x (13) q x, y, t = δ y y sin β m x + 1 e t, < y < (14) where α is the therml diffusivity of the mteril of the plte, k is the therml conductivity of the mteril of the plte, q is the internl het genertion nd δ r is well known dirc delt function of rgument r. Eq. (1) to Eq. (14) constitute mthemticl formultion of the prolem. III. SOLUTION To otin the expression for temperture T(x, y, t), we introduce the cosine integrl trnsform nd its inverse trnsform re T x, β m, t = K β m, y T x, y, t dy (15) T(x, y, t) = m =1 K β m, y T x, β m, t (16) where the kernel K β m, y = cos β m y (17) where β m is the m th root of the trnscendentl eqution sin β m =, β m = mπ, m = 1,,. On pplying the cosine integrl trnsform defined in the Eq. (15), its inverse trnsform defined in Eq. (16), pplying Lplce trnsform nd its inverse y residue method successively to the Eq. (1), one otins the expression for temperture s T x, y, t = where n=1 m=1 cosβ m y + α 1 n sin x Q t + α k Q 1 t = t e α β m + π n (t u) 1 + e u + e αβ m u αβ m 1 αβ m αβ m (αβm 1) Q t = t e α β m + π n nd (t u) α sin ( 1) n Q 1 t cosβ m y sin β m x + Q 3 t (18) cosβ m y α sin β m k F β m, u 1 + e u + e αβ m u αβ m 1 αβ m αβ m (αβm 1) Q 3 t = 1 + e t e + αβ m t αβ m 1 αβ m αβ m (αβm 1) F 1 β m, u du cosβ m y α sin β m k du w w w. j e r. o r g Pge 389

Americn Journl of Engineering Reserch (AJER) 13 Airy stress function U Using Eq.(18) in Eq.(1), one otins the expression for Airy s stress function U s U = t E m =1 n=1 cos β m y β m + π n α sin ( 1) n Q 1 t + α 1 n sin x Q t + α β m k cos(β m y ) sin β m x + Q 3 t β m + π n (19) Displcement nd Stresses Now using Eqs. (18) nd (19) in Eqs. () to (6) one otins the expressions for displcement nd stresses s u x = t m =1 n=1 cosβ m y cos α k 1+v β m u y = sin β m y t m=1 n=1 β m sin απ n 3 ( 1) n 1+v β m + π n Q 1 t + cos x Q t cos(β m y ) cos β m x + Q 3 t () α ( 1) n β m 1+v π β n Q 1 t + sin x Q t + α k 1 + v cos(β m y ) sin β m x + Q 3 t (1) ς xx = t E n=1 sin ( β m )cos β m y m=1 π β n α ( 1) n Q 1 t + sin x Q t + α kβ m β m + π n cos(β m y ) sin β m x + Q 3 t () ( π n )cos β m y m =1 π β n ς yy = t E n=1 sin + ς xy = t E cos + α β kπ n m + π n β m sin β m y m =1 n=1 π β n α ( 1) n Q 1 t + sin x Q t cos(β m y ) sin β m x + Q 3 t (3) n π α 3 ( 1) n Q 1 t + cos x Q t α β kβ π n cos(β m m y ) cos β m x + Q 3 t (4) IV. SPECIAL CASE AND NUMERICAL CALCULATIONS Setting f 1 y, t = f y, t = δ y y 1 δ t t, y 1, < t < F 1 β m, t = F β m, t = cos(β m y 1 ) δ t t = 1m, = m, t =,, 4, 6, 8 sec nd y = y 1 = 1m. 4.1 Mteril Properties The numericl clcultion hs een crried out for steel (.5% cron) rectngulr plte with the mteril properties defined s Therml diffusivity α = 14.74 1 6 m s 1, Specific het c ρ = 465 J/kg, Therml conductivity k = 53.6 W/m K, Poisson rtio θ =.35, w w w. j e r. o r g Pge 39

Americn Journl of Engineering Reserch (AJER) 13 Young s modulus E = 13 G p, Lme constnt μ = 6.67, Coefficient of liner therml expnsion t = 13 1 6 1 K 4. Roots of Trnscendentl Eqution The β 1 = 3.1414, β = 6.88, β 3 = 9.44, β 4 = 1.5656, β 5 = 15.77, β 6 = 18.8484 re the roots of trnscendentl eqution sin β m =. The numericl clcultion nd the grph hs een crried out with the help of mthemticl softwre Mt l. V. DISCUSSION In this pper thin rectngulr plte is considered which is free from frction nd determined the expressions for temperture, displcement nd stresses due to ritrry het supply on the edges x = nd x = of plte wheres the plte is insulted t y = nd y =. A mthemticl model is constructed y considering steel (.5% cron) rectngulr plte with the mteril properties specified ove. Fig. 1 Temperture T in Y- direction. Fig. The displcement u x in Y- direction. w w w. j e r. o r g Pge 391

Americn Journl of Engineering Reserch (AJER) 13 Fig. 3 The displcement u y in Y- direction. Fig. 4 Therml stresses ς xx in Y-direction Fig. 5 Therml stresses ς yy in Y-direction. w w w. j e r. o r g Pge 39

Americn Journl of Engineering Reserch (AJER) 13 Fig. 6 Therml stresses ς xy in Y-direction. From figure 1, it is oserved tht temperture T decreses s the time increses. The overll ehvior of temperture is decresing nd it is symmetric out y= 1 in thin rectngulr plte with internl het genertion long Y-direction. From figure, it is oserved tht displcement u x decreses s the time increses. The overll ehvior of temperture is decresing nd it is symmetric out y= 1 in thin rectngulr plte with internl het genertion long Y-direction. From figure 3, it is oserved tht the displcement u y is incresing for y.5, 1.5 y nd decresing for.5 y 1.5. The overll ehvior of displcement u y is decresing long Y-direction nd it is ntisymmetric out y= 1 in thin rectngulr plte with internl het genertion long Y-direction. From figure 4 nd 5, it is oserved tht therml stresses ς xx, ς yy increses s the time increses. Mximum vlue of stresses ς xx, ς yy occure ner het source nd it is tensile in nture in thin rectngulr plte with internl het genertion long Y-direction. From figure 6, it is oserved tht the therml stresses ς xy is decresing for x.5, 1.5 x nd incresing for.5 x 1.5. The overll ehvior of stresses ς xy is incresing nd it is ntisymmetric out y= 1 in rectngulr plte with internl het genertion long Y-direction. VI. CONCLUSION We cn conclude tht temperture T, displcement u x nd u y re decresing with time in thin rectngulr plte with internl het genertion long Y-direction. The therml stresses ς xx, ς yy re tensile in nture in thin rectngulr plte with internl het genertion long Y-direction. The therml stresses ς xy is incresing nd it is ntisymmetric out y = 1 in thin rectngulr plte with internl het genertion long Y- direction. The results otined here re useful in engineering prolems prticulrly in the determintion of stte of stress in thin rectngulr plte, se of furnce of oiler of therml power plnt nd gs power plnt. REFERENCES [1] V. S. Kulkrni nd K. C. Deshmukh, A rief note on qusi- sttic therml stresses in rectngulr plte, Fr Est J. Appl. Mth.6(3),7, 349-36. [] V. S. Gogulwr nd K. C. Deshmukh, Therml stresses in rectngulr plte due to prtilly distriuted het supply, Fr Est J. Appl. Mth.16(),4, 197-1. [3] M. Nsser, EI-Mghry, Two dimensionl prolem in generlized thermoelsticity with het sources, Journl of Therml Stresses, 7, 4, 7-39. [4] M. V. Khndit nd K. C. Deshmukh, Thermoelstic prolem in rectngulr plte with het genertion, Presented in Indin Mthemticl Society conference, Pune, Indi, 7. [5] V. B. Ptil, B. R. Ahirro nd N. W. Khorgde, Therml stresses of semi infinite rectngulr sl with internl het source, IOSR Journl of Mthemtics, 8(6), 13, 57-61. [6] M. N. Ozisik, Boundry Vlue Prolems of Het Conduction, Interntionl Text Book Compny, Scrnton, Pennsylvni, 1968. w w w. j e r. o r g Pge 393