Inverse quasi static unsteady state thermal stresses in a thick circular plate

Transcription

1 APPLIED & INTERDISCIPLINARY MATHEMATICS RESEARCH ARTICLE Inverse quasi static unsteady state thermal stresses in a thick circular plate Ishaque Khan, Lalsingh Khalsa and Vinod Varghese Cogent Mathematics (217, 4: Page 1 of 1

2 Khan et al., Cogent Mathematics (217, 4: APPLIED & INTERDISCIPLINARY MATHEMATICS RESEARCH ARTICLE Inverse quasi static unsteady state thermal stresses in a thick circular plate Ishaque Khan 1, Lalsingh Khalsa 1 and Vinod Varghese 1 * Received: 21 March 216 Accepted: 16 January 217 First Published: 19 January 217 *Corresponding author: Vinod Varghese, Department of Mathematics, M. G. College, Armori, Gadchiroli, Maharashtra, India vino7997@gmail.com Reviewing editor: Yong Hong Wu, Rashtrasant Tukadoji Maharaj Nagpur University, India Additional information is available at the end of the article Abstract: In this article, an attempt has been made to discuss the inverse thermoelastic problem of a thick circular plate defined as r a, h z h subjected to the arbitrary heat supply at interior point while circular edge of the thick circular plate at the outer surface and at the lower surface is maintained at zero temperature. The conductivity equation and the corresponding initial and boundary conditions have been solved using finite Hankel and Laplace integral transform techniques. Goodier s and Michell s functions are used to obtain the displacement components and its associated stresses. The results are obtained in a form in terms of Bessel s function. The results for unknown temperature, displacement, and stresses have been computed numerically considering special functions and illustrated graphically. Subjects: Applied Mathematics; Applied Mechanics; Inverse Problems; Mathematical Modeling; Mathematics & Statistics; Science Keywords: thick circular plate; thermoelasticity; unsteady-state; integral transform AMS subject classifications: 35B7; 35G3; 35K5; 44A1 1. Introduction As a result of the increased usage of industrial and construction materials, the interest in the inverse thermal stress problems has grown considerably, typified by main shaft of lathe and the role of the rolling mill due to the elementary geometry involved. As a result of this, a number of theoretical studies concerning them have been reported so far. However, to simplify this, almost all the studies were conducted on the assumption that the upper and lower surfaces of the circular plate are insulated or that the heat is dissipated with uniform heat transfer coefficients throughout the surfaces as direct problems. For example, Sabherwal (1965 investigated the inverse problem of transient Lalsingh Khalsa ABOUT THE AUTHORS Ishaque Khan is an assistant professor in Mathematics at M. G. College, Armori, Gadchiroli (MS, India and pursuing research work in the field of Thermoelasticity BVP. Lalsingh Khalsa has published 21 research papers and 5 books at international level, besides being an excellent supervisor. Vinod Varghese has also extended his full support in assisting the work and has lent a helping hand in providing it with a complete touch. PUBLIC INTEREST STATEMENT The results are useful to scientists working on two-dimensional thick plate theories, especially when circumstances are of inverse quasi-static transient conditions. Inverse problems are very crucial, especially in view of its relevance to various industrial machines and it can be the theoretical base for several geophysical and industrial applications. The physical applications are encountered in the context of problems such as ground explosions and oil industries, where the interest is in various phenomena occurring during the measurement of displacements, stresses, and temperature field due to the presence of certain sources. 217 The Author(s. This open access article is distributed under a Creative Commons Attribution (CC-BY 4. license. Page 2 of 1

3 Khan et al., Cogent Mathematics (217, 4: heat conduction in circular plate. Grysa, Ciałkowski, and Kamiński (1981 discussed on an inverse temperature filled problem of the theory of thermal stresses. Noda (1989 discussed an analytical method for an inverse problem of three-dimensional transient thermoelasticity in a transversely isotropic solid by integral transform technique with newly designed potential function and illustrated practical application of the method in engineering problem. Ashida, Choi, and Noda (1996 investigated an inverse thermoelastic problem in an isotropic structural plate onto which a piezoelectric ceramic plate is perfectly bonded. When unknown heating temperature acts on the free surface of the isotropic structural plate, an electric potential is induced in the piezoelectric ceramic plate. Deshmukh and Wanhede (1997, 1998a, 1998b discussed the inverse transient problem of quasistatic thermal deflection in these clamped circular plates and axisymmetric inverse steady-state problem of thermoelastic deformation of finite length hollow cylinder and inverse quasi-static transient thermoelastic problem in a thin annular disk. Again Ashida et al. (22 emphasized on the inverse transient thermoelastic problem for a composite circular disk. Yang, Chen, Chang (22 studied inverse boundary value problem of coupled thermoelasticity in an infinitely long annular cylinder using simulated exact and inexact measurements. Patil and Prasad (213 studied inverse steady-state thermoelastic problem of a thin rectangular plate using operational methods. From the previous literatures regarding thick plate as considered, it was observed by the author that no analytical procedure has been established for thick circular plate, considering inverse quasi-static thermoelastic analysis. In this problem, we consider some new interesting results of the inverse heat conduction problem of thick circular plate occupying the space D {(x, y, z R 3 : (x 2 + y 2 1/2 a, h z h}, where r (x 2 + y In a condition wherein a thick circular plate is subjected to arbitrary heat supply at interior point while the circular edge of the thick circular plate at the outer surface and at the lower surface is maintained at zero temperature, the governing heat conduction equation has been solved using integral transform method. The results are obtained in series form in terms of Bessel s functions. The mathematical model of final thick circular plate has been constructed with the help of numerical illustrations. 2. Formulation of the problem Consider a thick circular plate of thickness 2 h occupying space D defined by r a, h z h, as shown in Figure 1. Let the plate be subjected to an arbitrary known interior temperature f(r, t within the region h z h. With lower surface and circular surface, r a at zero temperature. Under this more realistic prescribed condition, the unknown temperature g(r, t which is at the upper surface of the plate z h and quasi-static thermal stresses due to unknown temperature g(r, t are to be determined. Figure 1. Geometrical configuration of the problem. Page 3 of 1

4 Khan et al., Cogent Mathematics (217, 4: Temperature distribution The transient heat conduction equation of the plate is given as follows: 2 T r + 1 T 2 r r + 2 T z 1 T 2 κ t (1 where κ is the thermal diffusivity of the material of the disk (which is assumed to be constant, subjected to the initial and boundary conditions: T at t (2 T at r a (3 T at z h T f (r, t at z ξ T g(r, t at z h r a (Unknown (4 (5 ( Thermal displacements and thermal stresses Following Noda et al. (23, we assume that the Navier s equations in the absence of body forces for axisymmetric two-dimensional thermoelastic problem as: 2 U r U r r + 1 e 2(1+υ 1 2υ r 2 U z 1 e 2(1+υ 1 2υ z α T 1 2υ t r α T 1 2υ t z (7 where U r and U z are the displacement components in the radial and axial directions, respectively, and the dilatation e as: e U r r + U r r + U z z (8 The displacement function in the cylindrical coordinate system is represented by Goodier s thermoelastic displacement potential ϕ and Michel s function M: U r φ r 2 M r z, U z φ z + 2(1 υ 2 M 2 M 2 z (9 (1 in which Goodier s thermoelastic potential must satisfy: 2 φ Kτ with φ at t. (11 and the Michel s function M must satisfy: 2 2 M (12 in which K is the restraint coefficient and temperature change τ T T i, T i is the initial temperature, and: 2 2 r r r + 2 z. 2 The component of the stresses is represented as: Page 4 of 1

5 Khan et al., Cogent Mathematics (217, 4: [ 2 φ σ rr 2G r kτ + ( ] ν 2 M 2 M 2 z r 2 σ θθ 2G [ 1 φ r r kτ + ( ν 2 M 1 z r ] M r [ 2 ( ] φ σ zz 2G kτ + (2 v 2 M 2 M z 2 z z 2 [ 2 φ σ rz 2G r z + ( ] (1 v 2 M 2 M r z 2 (13 (14 (15 (16 where G and ν are shear modulus and Poisson s ratio, respectively; for the traction-free surface, the stress function is: σ rr σ rz at r a (17 Equations (1 to (17 constitute the mathematical formulation of the problem. 3. Solution of the problem 3.1. Solution for temperature distribution In order to solve Equation (1 under the boundary condition (3, we firstly introduce the finite Hankel transform of order m over the variable r; the integral transform and its inversion theorem (Sneddon, 1972 can be written as: a T, z, t rj, rt(r, z, tdr, ( 2J (α T(r, z, t m, r T(α m1 J 2 (α a m, z, t 1 m (18 in which α 1, α 2, α 3, are the roots of the transcendental equation J a. Applying the finite Hankel transform and Laplace integral transform, and its inversion theorems, results in the final temperature distribution as: 4( 1 n+1 nπ kj T(r, z, t r sin(β n (ξ + h 2 J 2 (α a 1 m t exp( k u f, t u du (19 The function given in Equation (23 represents the temperature at every instant and at all points of the circular thick plate of finite height. The unknown temperature g(r, t can be obtained by substituting z h in the equation. g(r, t T (r z, t at z h, one obtains: 2( 1 n+1 nπ kj g(r, t r sin(β n (ξ + h 2 J 2 (α a 1 m in which: t exp( k u f, t u du (2 Page 5 of 1

6 Khan et al., Cogent Mathematics (217, 4: β n nπ ( ( 2 z + h, ξ + h α 2 nπ + m ξ + h 3.2. Solution for thermal stresses Goodier thermoelastic displacement potential (φ Referring to the fundamental Equation (1 and its solution (19 for the heat conduction problem, the solution for the displacement function is represented by Goodier s thermoelastic displacement potential φ governed by Equation (11; this is represented by: 4kJ φ(r, z, t r ( 1 n sin(β n J 2 (α a 1 m (ξ + h 2 t [exp( k u f ( α m, t u ] du ( Michell s function M Similarly, the solution for Michell s function M is assumed so as to satisfy the governed condition of Equation (12 as: M 4k ( J αm r ( αm a { ( ( } H mn sinh[α m z + h ]+Rmn α m (z + h cosh[α m z + h ] J 2 in which H mn and R mn are arbitrary functions. ( Displacement and thermal stresses In this manner, two displacement functions in the cylindrical coordinate system ϕ and M are fully formulated. Now in order to obtain the displacement components, we substitute the values of thermoelastic displacement potential φ and Michell s function M in Equations (9 and (1, which results in: ( 4k ( α U r m ( 1 n { t nπ J 1 r sin(β n (ξ + h 2 J 2 (α a [exp( k u f, t u]du 1 m + H mn α 2 cosh [α (z + h] + R m m mn α2 m cosh[αm (z + h] + α m (z + h sinh[α m (z + h] } (23 ( 4k 2 ( 1 n n 2 π 2 J U z r (ξ + h 2 J 2 (α a 1 m { cos(β n H mn α 2 sinh[α m m (z + h] + R mn α 2 m 2 ( 1 2v sinh[α m (z + h] α m (z + h cosh h[α m + h] } ( { 4kG 1 σ rr J 2 (α a m 2 J r r 1 J 1 r ( 1n+1 t nπ sin(β n [exp( k (ξ + h 2 u f, t u]du 4k ( 1 n+1 t nπ J k r sin(β n J 2 (α [exp( k 1 ma(ξ + h2 u f, t u]du +H mn α 2 m J r r 1 J 1 r cosh[α m (z + h] +R mn α 2 m[ 2υαm, J ( αm r cosh [ α m (z + h] + J r r 1 J 1 r cosh[2α m (z + h] sinh[α m (z + h] } (24 (25 Page 6 of 1

7 Khan et al., Cogent Mathematics (217, 4: ( { 4k J σ θθ 1, r 2 ( 1 n+1 n πα m sin(β n F(t J 2 (α 1 m a (ξ + h 2, 4k( 1 n+1 t nπj k r sin(β n J 2 (α [exp( k 1 ma(ξ + h u f, t u]du +H mn r 1 α 2 J (α r cosh[α (z + h] + R m 1 m m mn m α2 2v αm J r ( cosh[α m (z + h] + r 1 J 1 r(cosh[α m (z + h] +α m z + h sinh[αm (z + h] } (26 ( 4k σ zz k 2( 1 n+1 n 2 π 3 J r sin(β n ( 4J ξ + h 2 (α a 1 m ( k ( 1 n+1 nπ J αm r t (ξ + h 2 J 2 1( αm a t [exp( k u f, t u]du [exp( k u f, t u]du H mn α 3 m cosh[α m (z + h] + R mn α3 m (1 2v cosh[αm (z + h] α m (z + h sinh[α m (z + h] } ( { 4k J σ rz 1 r ( 1 n+1 2n 2 π 2 t cos(β n J 2 (α a [exp( k m (ξ + h 3 u f, t u]du +H mn α 3 m sinh[α m (z + h] + R mn α3 m 2v sinh[αm (z + h] +α m (z + h cosh[α m [α m (z + h]] } (27 ( Determination of unknown arbitrary function H mn and R mn Applying boundary condition Equation (17 to Equations (25 and (28, one obtains: R mn ( 1n 2n 2 π 2 t cos(β n (ξ + h 3 α 3 [exp( k u f, t u]du m 1 sinh[α m (z + h] + α m (z + h cosh[α m (z + h] (29 H mn ( 1n+1 2n 2 π 2 t cos(β n (ξ + h 3 α 3 [exp( k u f, t u]du m 1 2υ sinh[α m (z + h] + α m (z + h cosh[α m (z + h] (3 4. Special case and numerical calculations Setting f (r, t (r 2 2 (1 e t (31 Applying finite Hankel transform to the Equation (31, one attains: f ( α m, t a r (r 2 2 J, r(1 e t dr (32 8a {(8 α 2 m J 1 a 4aα m J a} (1 e t α 5 m Page 7 of 1

8 Khan et al., Cogent Mathematics (217, 4: Figure 2. Temperature distribution. Figure 3. Radial displacement profile. Figure 4. Axial displacement profile. Figure 5. Radial stress distribution. 5. Numerical calculations The numerical calculation has been carried out for steel (SN 5C plate with parameters a 1 m, h.2 m, thermal diffusivity k (m 2 s 1, and Poisson s ration v.281, with α , Page 8 of 1

9 Khan et al., Cogent Mathematics (217, 4: Figure 6. Tangential stress distribution. α , α , α , α , α , α , α , α , and α being the roots of transcendental equation J a. For convenience, set A a, B 16K/1 2 a, and C 32GK/1 2 a in the expressions for obtaining the unknown temperature, displacement, and stress components. In order to examine the influence of unknown temperature on the upper surface of circular plate, the numerical calculation z h/2, r,.2, 4., 6., 8., & 1 and ξ.2,.1 and.1 was performed. Numerical variations in radial directions have been illustrated in the figure with the help of a computer program. 6. Concluding remarks In this problem, a thick circular plate is considered which is kept traction-free as well as subjected to arbitrary known interior temperature and determined for the expressions of unknown temperature, displacements, and stress functions due to the unknown temperature. As a special case, mathematical model is constructed for f (r (r 2 2 (1 e t and numerical calculations were performed. The thermoelastic behaviors such as temperature, displacements, and stresses are examined with the help of arbitrary known interior temperature along the radial direction as a.2, b.1, c, d.1. Figure 2 indicates that the unknown temperature decreases from r to r.3 and increases from.3 to 1 with the thickness of the circular plate. As the source of known temperature varies from a negative to positive value, the unknown temperature decreases its magnitude along the radial direction. As shown in Figure 3, the source of known temperature varies from bottom to top, the radial displacement decreases at r, and the radial displacement vanishes; otherwise, its existence would have been visible. As shown in Figure 4, the source of known temperature varies from bottom to top; the axial displacement increases along radial direction; and it shows its existence. Figure 5 shows that the radial stress decreasing from bottom to (lower surface to upper surface top. Stress at r and r a is zero; otherwise, it shows its existence. Figure 6 indicates that the stress function σ θθ decreases with the thickness of the circular plate. It shows the existence for small thickness. Also, it develops tensile stresses in the radial direction. In this article, we analyzed an inverse thermoelastic problem of a thick circular plate and determined the expressions of unknown temperature, displacement, and thermal stresses. The heat conduction differential equation is solved using finite Hankel and Laplace integral transform techniques, and their inversion theorems. Goodier s and Michell s functions are used to obtain the displacement components. As a special case, a mathematical model is constructed for steel (SN 5C thick plate, Page 9 of 1

10 Khan et al., Cogent Mathematics (217, 4: with the material properties specified as above, and examined for the thermoelastic behaviors in unsteady-state field for unknown temperature change, displacement, and thermal stresses. We conclude that the displacement and stress components occur near the heat source region. As the temperature increases, the circular plate will tend to expand in the radial direction as well as in the axial direction. Also any particular case of special interest can be derived by assigning values to the parameters and functions in the expressions (19 (28. Acknowledgments The author offers his grateful thanks to the Editor and the referee for their kind help and active guidance in the preparation of this revised paper. Funding The author is thankful to University Grants Commission, New Delhi, India, for providing partial financial assistance under Minor Research Project Scheme [File No /12(WRO]. Author details Ishaque Khan 1 iakhan_get@rediffmail.com Lalsingh Khalsa 1 lalsinghkhalsa@yahoo.com Vinod Varghese 1 vino7997@gmail.com ORCID ID: 1 Department of Mathematics, M. G. College, Armori, Gadchiroli, Maharashtra, India. Citation information Cite this article as: Inverse quasi static unsteady state thermal stresses in a thick circular plate, Ishaque Khan, Lalsingh Khalsa & Vinod Varghese, Cogent Mathematics (217, 4: Cover image Source: Authors. References Ashida, F., Choi, J. C., & Noda, N. (1996. An inverse thermoelastic problem in an isotropic plate associated with a piezoelectric ceramic plate. Journal of Thermal Stresses, 19, Ashida, F., Sakata, S., Tauchert, T. R., & Yamashita, Y. (22. Inverse transient thermoelastic problem for a composite circular disk. Journal of Thermal Stresses, 25, Deshmukh, K. C., & Wankhede, P. C. (1997. An axisymmetric inverse steady state problem of thermoelastic deformation of a finite length hollow cylinder. Far East Journal of Applied Mathematics, 1, Deshmukh, K. C., & Wankhede, P. C. (1998a. An inverse transient problem of quasi-static thermal deflection of a thin clamped circular plate. Bulletin of Pure and Applied Sciences, 17, Deshmukh, K. C., & Wankhede, P. C. (1998b. An Inverse quasistatic transient thermoelastic problem in a thin annular disc. Far East Journal of Applied Mathematics, 2, Grysa, K., Ciałkowski, M. J., & Kamiński, H. (1981. An inverse temperature field problem of the theory of thermal stresses. Nuclear Engineering and Design, 64, Noda, N. (1989. An inverse problem of coupled thermal stress filled in a circular cylinder. JSME, International Journal, 32, Noda, N., Hetnarski, R. B., & Tanigawa, Y. (23. Thermal stresses (2 ed., pp New York, NY: Taylor and Francis. Patil, S., & Prasad, J. S. V. R. K. (213. Some steady-state thermoelastic problems of rectangular plates. Applied Mathematical Sciences, 7, Sabherwal, K. C. (1965. An inverse problem of transient heat conduction. Indian Journal of Pure Applied Physics, 3E, Sneddon, I. N. (1972. The use of integral transform (pp New York, NY: McGraw Hill. Yang, Y. C., Chen, U. C., & Chang, W. J. (22. An inverse problem of coupled thermoelasticity in predicting heat flux and thermal stresses by strain measurement. Journal of Thermal Stresses, 25, The Author(s. This open access article is distributed under a Creative Commons Attribution (CC-BY 4. license. You are free to: Share copy and redistribute the material in any medium or format Adapt remix, transform, and build upon the material for any purpose, even commercially. The licensor cannot revoke these freedoms as long as you follow the license terms. Under the following terms: Attribution You must give appropriate credit, provide a link to the license, and indicate if changes were made. You may do so in any reasonable manner, but not in any way that suggests the licensor endorses you or your use. No additional restrictions You may not apply legal terms or technological measures that legally restrict others from doing anything the license permits. Page 1 of 1