Non-Fourier Heat Conduction of a Functionally Graded Cylinder Containing a Cylindrical Crack

Fu, Jiawei; Hu, Keqiang; Qian, Linfang; Chen, Zengtao

doi:https://doi.org/10.1155/2020/8121295

Advances in Mathematical Physics

On this page

Abstract Introduction Conclusion Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 8121295 | https://doi.org/10.1155/2020/8121295

Non-Fourier Heat Conduction of a Functionally Graded Cylinder Containing a Cylindrical Crack

Jiawei Fu,¹Keqiang Hu,²Linfang Qian,¹and Zengtao Chen³

Academic Editor: P. Areias

Received27 Sept 2019

Revised13 Nov 2019

Accepted19 Nov 2019

Published01 Feb 2020

Abstract

The present work investigates the problem of a cylindrical crack in a functionally graded cylinder under thermal impact by using the non-Fourier heat conduction model. The theoretical derivation is performed by methods of Fourier integral transform, Laplace transform, and Cauchy singular integral equation. The concept of heat flux intensity factor is introduced to investigate the heat concentration degree around the crack tip quantitatively. The temperature field and the heat flux intensity factor in the time domain are obtained by transforming the corresponding quantities from the Laplace domain numerically. The effects of heat conduction model, functionally graded parameter, and thermal resistance of crack on the temperature distribution and heat flux intensity factor are studied. This work is beneficial for the thermal design of functionally graded cylinder containing a cylindrical crack.

1. Introduction

Cylindrical composite structures are widely used in engineering, and the interfaces in the composites are key parts but at the same time weak regions, considering that these weak regions are generally subjected to various damages such as debonding or cracking. When a cylindrical interface debonds, the crack may occupy the whole angular range which is called a cylindrical crack [1]. Functionally graded (FG) materials have received considerable attention in engineering fields due to the fact that FG materials possess gradual change in composition and microstructure and have continuously spatial variations in thermal and mechanical properties. The application of FG materials can improve the bonding strength and toughness and reduce the probability of thermomechanical problems caused by the mismatch of material properties. Due to their important characteristics, FG materials found various applications, such as: inner walls of nuclear reactor, pressure vessels, thermal barrier coatings for combustion chambers, heat exchanger tubes, etc., [2]. The common fabrication technologies for producing FG bulk materials include powder metallurgy, centrifugal casting, tape casting, and additive manufacturing [3]. And, it has been found that embedded microcracks or cracks are prone to occur during the manufacture process of structures or in their service [4–6]. For instance, cylindrical cracks may initiate along the interface of fiber and matrix in the fiber reinforced composite cylinders [7]. Thus it is a significant design issue to study the thermal or mechanical fracture behaviors of FG structures [8, 9].

Great efforts have been made to investigate the fracture behavior of FG materials under mechanical and/or thermal loading [10–15]. Jin [16] studied the singular crack-tip fields in nonhomogeneous body under thermal stress and pointed out that the stress intensity factors (SIFs) were still applicable in fracture problems of FG materials. Li and Weng [17] investigated the dynamic problem of a cylindrical crack in a functionally graded interlayer bonded by two dissimilar, homogeneous cylinders under torsional impact loading. Jin and Paulino [18] studied a crack in a viscoelastic strip of a FG material under tensile loading and obtained both mode I and mode II SIFs and the crack opening/sliding displacements. Feng et al. [19] considered the problem of a cylindrical crack in the nonhomogeneous, interfacial zone in a composite loaded by torsional impact. The plane, thermo-mechanical behavior of a crack in a viscoelastic FG material coating with arbitrary material properties bonded to homogeneous substrate has been studied by Cheng et al. [20]. The problem of a cylindrical interface crack in a hollow, functionally graded cylinder under static torsion has been studied by Shi [21], and the coupled effects of geometrical, physical, and functionally graded parameters on the interfacial fracture behavior have been demonstrated by numerical results. Zhang et al. [22] analyzed a thermal shock problem of an elastic strip made of functionally graded materials containing a crack parallel to the free surface based on the fractional heat conduction theory.

Considering that cylindrical composites under thermal loading conditions have significant engineering applications, the thermal analysis of the composite structures with cracks is of great importance for the safety design. The classical Fourier heat conduction model implies that the thermal wave can propagate instantaneously in the media at an infinite speed, which contradicts physical facts. For crack problems, the application of non-Fourier heat conduction models on homogeneous materials has been reported in some recent literatures [23–27]. Hu and Chen [25] investigated a partially, thermally insulated crack in an elastic strip under a thermal impact loading and studied the possible crack kinking phenomena. Wang [26] considered the transient thermal cracking associated with non-classical heat conduction in a cylindrical coordinate system. Hu and Chen [28] analyzed the transient heat conduction of a cracked half-plane using dual-phase-lag (DPL) theory. The fracture behavior of a thermoelastic cylinder subjected to a sudden temperature change on its outer surface has been studied by Fu et al. [27] using the hyperbolic heat conduction model. The fractional calculus has been introduced into the hyperbolic heat conduction model to investigate the fracture problem of a circumferential crack in a hollow cylinder under thermal shock by Zhang and Li [29]. Yang and Chen [30] established a thermoelastic model for a functionally graded half-plane containing a crack under thermal shock loading in the framework of hyperbolic heat conduction theory.

To the authors’ best knowledge, the thermal problem of a cylindrical crack in a functionally graded cylinder has not been solved due to the mathematical complexity involved. This work is dedicated to obtain the dynamic thermal responses of an FG cylinder containing a cylindrical crack under thermal shock using non-Fourier heat conduction models. The Fourier transform and Laplace transform are applied, and the crack problem is reduced to solving a singular integral equation in the Laplace domain. Then the heat flux intensity factor (HFIF) and temperature field are obtained by applying the numerical Laplace inversion technique. Finally, the effects of heat conduction model, functionally graded parameter, and thermal resistance of crack on the temperature field and HFIF are analyzed numerically.

2. Problem Formulation

As illustrated in Figure 1, an FG hollow cylinder located in the cylindrical coordinate is subjected to transient thermal loadings on its inner and outer surfaces. The infinitely long, rigid cylinder has an inner radius , and outer radius , while an axisymmetric, cylindrical crack spreads along the plane with the size of . For simplicity, the thermal conductivity and mass density of the cylinder are assumed to be radial-coordinate-dependent, while other material properties are homogeneous. The non-Fourier heat conduction model is adopted to analyze the thermal behavior of the cracked FG cylinder with uniform initial temperature .

Considering the axisymmetric nature of the cracked cylinder, the heat flux only exists in the radial, and axial directions, and can be expressed within the structure of the DPL model [28] as,

in which, is the heat flux vector, is time, is temperature, and is the thermal conductivity. The phase lag of heat flux and phase lag of temperature gradient are material properties. By omitting the quadratic term on the left and setting in Eq. (1), another frequently-used, non-Fourier heat conduction model proposed by Cattaneo [31] and Vernotte [32] (C–V model), is obtained. Also, setting and leads to the classical Fourier heat conduction model. The advancement of the DPL heat conduction model has been clarified by predicting the temperature results tested for some materials, such as, processed meat [33]. Phase lags of heat flux and temperature gradient are material-dependent properties and vary from picoseconds for most of the metals [34] to 10² s for nonhomogeneous materials [33, 35].

Without heat source, the energy conservation equation reads

where and are mass density and specific heat capacity, respectively.

By eliminating the heat flux from Eqs. (1) and (2), the 2D axisymmetric governing equation of temperature is obtained

where the thermal conductivity and mass density vary with radial coordinate in a power-law form

in which, and are reference material properties, and the power exponent can also be called as FG parameter.

Substitution of Eq. (4) into (3) leads to

where is the reference thermal diffusivity.

The boundary conditions for the partially insulated crack can be written as

where and with the superscripts “(1)” and “(2)” represent temperature and heat flux in the area and , respectively. Equations (8)–(9b) describe the thermal condition along the crack plane, while Eqs. (9a) and (9b) is termed as the mixed boundary condition which enhances the complexity of the problem remarkably. When the thermal resistance of the crack equals 0, the crack will not disturb the temperature field, and a perfectly conductive crack condition arises. When approaches infinity, there is no heat flux across the insulating crack.

3. Solution Procedure

The following non-dimensional parameters are introduced to simplify the solution procedure:

It is noted that other parameters with the same units as those in Eq. (10) are omitted here for brevity.

Using the parameters in (10), the governing Equation (5) can be normalized as

Applying the Laplace transform to the above equation by assuming the zero initial conditions leads to

in which, and is the Laplace variable.

Similar normalization and Laplace transformation processes applied to the boundary conditions give,

Using superposition, the temperature field resulting from differential Equation (12) under the boundaries (13)~(16b) can be the addition of solutions of two subproblems: () inhomogeneous boundary condition applied to the cylinder surfaces with the absence of crack; () homogeneous thermal load applied to the cracked cylinder. For the first problem, the temperature distribution along the axial direction will be uniform. Thus one could have.

:

and

:

The total temperature and heat flux can be easily obtained from and after solving the above two problems.

can be directly solved as,

where, the imaginary number , , and . and represent the th-order modified Bessel functions of the first and second kinds, respectively. The unknowns and can be obtained from , in which, is defined as

and the non-zero elements of the 4 × 4 matrix and 4 × 1 vector are

Now, let’s turn to solve the mixed-boundary-value problem . Applying the Fourier transform to Equation (22) with respect to leads to,

The temperature and the heat flux are then obtained as,

where, .

Using the boundary Equations (23) and (24), the unknowns in Eq. (32) are linked as

in which,

The only unknown, independent coefficient can be solved from the mixed boundary conditions (25a) and (25b). A function, defined as the gradient of temperature jump, is introduced

Equation (36), expressing the axial gradient of temperature jump across the crack faces, serves as a basic function to solve the mixed-boundary-condition problem using the method of singular integral equation. Due to the existence of crack, the temperature along the plane jumps in the radial direction within the range of , while it is continuous within the range of . Thus, the axial gradient of the temperature jump can describe the nature of thermal disturbance caused by the crack in a unified form.

Thus from Eq. (25a), one could have

Eq. (38) is the single-valuedness condition. With the help of Eq. (34), substituting Eq. (32) into Eq. (36) leads to

in which, is the integral variable and,

Substituting Eqs. (27), (32)–(34), and (39) into Eq. (25b) results in a singular integral equation

in which,

The fundamental solution of the singular integral equation (41) under the single-valuedness condition (38) is [36],

Eqs. (41) and (38) can be solved numerically by adopting the Gauss–Jacobi integration formulas in [37] as,

in which,

After calculating the function numerically, the four unknown coefficients and can be obtained from Eqs. (34), (39), and (48). Finally, problem is solved.

4. Temperature Field and Heat Flux Intensity Factor

The total temperature in the cylinder can be obtained by adding Equations (26)–(32) as,

in which, , , and are expressed in terms of , while can be numerically obtained from Eq. (39) using the Gauss–Chebyshev integration equation

The concept of heat flux intensity factor (HFIF), similar to the stress intensity factor, is introduced to describe the heat gathering degree around the crack tip

Using the parameter , the non-dimensional HFIF in the Laplace domain is written as

Substituting Eqs. (27) and (33) into the above equation, one could have

in which, is calculated from using the interpolation method.

Finally, the transient temperature and HFIF in the time domain are obtained by applying the numerical Laplace inversion technique described in [38] to Eqs. (53) and (57), respectively.

5. Numerical Results

In this section, a fiber reinforced composite cylinder containing a cylindrical crack is taken as an example to demonstrate the applicability of the mathematical model proposed in this paper. A sudden temperature change is applied on the inner surface of cylinder, and transient temperature field and HFIF are obtained based on the DPL heat conduction model. The effects of heat conduction model and phase lags, FG parameter, and thermal resistance of crack on the results are analyzed in detail. In the meanwhile, the solution procedure is verified in Section 5.1 with the temperature result presented in [39] which described the non-Fourier heat conduction in uncracked FG cylinders.

5.1. Validation of the Solution Procedure

In order to justify the solution process in this paper, a cracked FG cylinder with the same configuration and temperature loading as the uncracked cylinder in [39] is taken for illustration. The mass density and thermal conductivity vary with the radial coordinate gradually with , while other parameters are homogeneous. The crack spreads along the plane and has the nondimensional length . It is known from the Saint-Venant principle that the temperature field near the crack will be disturbed and not be affected by the crack at a position far away from it. Thus, the temperature field far from the crack should be the same as that for uncracked cylinders. The temperature distributions along the radial direction at are drawn in Figure 2 and compared with the results of uncracked cylinder at two different instants and . It is seen that the thermal wave propagates from the outer surface of the cylinder to the inner surface under the DPL heat conduction model. Evidently, the difference between the transient temperature field away from the cylindrical crack obtained in this paper and that in the reference for uncracked cylinder is negligible, which proves the correctness of the solution procedure.

5.2. Effect of Heat Conduction Model and Phase Lag Parameters

In this subsection and thereafter, a cylinder with the size of and containing a cylindrical crack of length 2c at is chosen for the parameter analyses. The cylinder is made of carbon fiber reinforced resin matrix composite, and the reference properties are chosen with 20% volume fraction of fiber, shown in Table 1. It should be noted that the phase lag of heat flux is chosen referring to [40], and the phase lag of temperature gradient is assumed correspondingly due to the absence of this value in literature. The crack length is . The inner surface of the cylinder is heated to suddenly, while the outer surface is kept at the initial temperature . The effects of heat conduction model and corresponding phase lags on the transient heat flux intensity factor and temperature field are analyzed in this subsection.

The histories of HFIF calculated based on different heat conduction models for a homogeneous cylinder with an insulated cylindrical crack are depicted in Figure 3. As the HFIF represents the degree of heat concentration around the crack tip, a higher value of HFIF means a more dangerous thermal situation. It can be seen from Figure 3 that the value of at stable state is the same for all the heat conduction theories, however, the maximum value of strongly depends on the heat conduction theory adopted. Compared to the HFIF history obtained by the Fourier model, oscillates first and has a higher maximum value for non-Fourier models. Also, the DPL model induces a higher than the C–V model, which declares the importance of DPL heat conduction model to predict the thermal concentration degree around the crack.

The effect of phase lag of heat flux on the HFIF history and temperature histories of the crack face midpoints of a homogeneous cylinder are shown in Figures 4 and 5, respectively. Three different values of within the structure of C–V model are chosen for the analysis. It can be seen from Figure 4 that the increase of raises the oscillation range of HFIF, while its influence on the maximum value of HFIF is indiscernible. Figure 5 clearly shows the temperature jump across the crack faces, and the temperature at the inner crack face is much higher than that at the outer crack face. Due to the existence of the insulated crack, the thermal wave will be reflected back when it propagates from the inner surface to the outer surface. Then, more heat is gathered around the inner crack face, and the temperature is elevated substantially to a value higher than the applied loading . Similar to the conclusion given by Akbarzadeh and Chen [39], the increase of delays the start signal of the thermal wave and increases the maximum temperature of the two crack faces.

Figures 6 and 7 illustrate the effect of phase lag of temperature gradient on the transient heat flux intensity factor and temperature field, respectively, within the structure of DPL heat conduction model. The cylinder is taken to be homogeneous, while the phase lag of heat flux is kept constant at and is assumed to vary from to . Different from the effect of on the maximum in Figure 4, the increase of generally elevates the maximum value of HFIF, which is clearly shown in Figure 6. It reveals that a cracked structure made of a material with a higher is more dangerous under sudden temperature change. By comparing Figures 5 and 7, the temperature at crack face midpoints obtained using the DPL model needs more time to get stable than that using the C–V model. Moreover, the influence of on the maximum temperature is not obvious.

5.3. Effect of FG Parameter

The thermal conductivity and mass density of cylinder vary with the radial coordinate in a power-law form as shown in Eq. (4), while other parameters are homogeneous. Different values of FG parameter (power exponent) correspond to a cylinder made of materials with different properties, thus affects the thermal responses remarkably. In this subsection, the effects of FG parameter on the HFIF history, temperature history of crack face midpoints, and temperature distribution along the crack plane are shown in Figures 8~10, respectively. It is worth mentioning that the DPL heat conduction model is taken for the thermal analyses.

It can be seen from Figure 8 that the rise of inhomogeneous degree increases the heat flux intensity factor entirely, including the oscillation range, maximum value, and the stable value. This reveals that the employment of inhomogeneous cylinders with positive FG parameter is more thermally dangerous than the homogeneous ones. Meanwhile, an inhomogeneous cylinder with negative FG parameter can serve as a nice thermal barrier structure as it reduces the thermal concentration degree around the crack remarkably. Different from the situation of HFIF history, the effect of FG parameter on the temperature history of midpoints of crack faces is weak, which is depicted in Figure 9. Due to the use of DPL heat conduction model, the temperature jump across the crack faces shows clear oscillation behavior. The temperature jump is more obvious in Figure 10, which illustrates the temperature distribution along the crack plane at a specific instant for different FG parameters. Clearly, the temperature at one crack face differs from another within the range of crack , and this phenomenon becomes continuously beyond this range. Also, the temperature at the inner crack face is much higher than that at the outer crack face due to the outwards thermal flow direction. Figure 10 also indicates that has sensible effect on the temperature of inner crack face and negligible effect on the temperature of the outer crack face, as the thermal wave has not arrived at the outer crack face at the instant due to the thermal reflection by the insulated crack face.

5.4. Effect of Thermal Resistance of Crack

The effect of thermal resistance of crack on the transient HFIF and temperature field of crack face midpoints are shown in Figures 11 and 12, respectively. The cylinder is taken to be inhomogeneous with , and the thermal results are obtained based on the DPL model. As stated previously, an infinitely large value of corresponds to an insulated crack, while two other partially insulated cracks with and are chosen for analysis.

With the decrease of , the crack tends to be thermally conductive. As shown in Figure 11, a smaller value of thermal resistance of crack leads to a lower HFIF, which indicates that a crack with higher conductivity leads to a lower heat concentration degree at the crack tip. According to the temperature history shown in Figure 12, the temperature jump across crack face shrinks with the decrease of . Clearly, the temperature difference between the two crack faces of the insulating crack () is much higher than that of the nearly conducting crack () for the whole dynamic equilibrium process. Similar conclusions have been reported by Hu and Chen [28]. It is also obvious that the temperature oscillates first and tends to be steady because of the significant non-Fourier thermal effect.

6. Conclusion

The transient heat conduction of a functionally graded cylinder containing a cylindrical crack is investigated in this paper using the non-Fourier heat conduction theory. The mixed-boundary-value problem is solved with the employment of the singular integral equation and Laplace transform methods. The concept of heat flux intensity factor is introduced to evaluate the heat concentration degree at the crack tip. Finally, the effects of heat conduction model, FG parameter, and thermal resistance of the crack on the transient HFIF and temperature field are analyzed numerically.

It is found that the temperature field across the two crack faces is discontinuous due to the existence of insulating crack. The heat conduction models adopted and related phase lags have great effects on the thermal responses. Due to the use of non-Fourier models, the HFIF are elevated remarkably, and the temperature shows significant oscillation behavior. Within the structure of FG profile used in this paper, an inhomogeneous cylinder with negative FG parameter can serve as a better thermal barrier structure than the homogeneous one. Also, the thermal resistance of crack is an important parameter which reduces the HFIF and temperature jump across the crack faces when the crack becomes more conductive.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work is partially supported by the National Natural Science Foundation of China (grant number 11702137), Natural Science Foundation of Jiangsu Province of China (grant number BK20170816), and the Fundamental Research Funds for the Central Universities (grant number 309171B8802).

References

Y.-D. Li, H. Zhao, and T. Xiong, “The cylindrical interface crack in a layered tubular composite of finite thickness under torsion,” European Journal of Mechanics–A/Solids, vol. 39, pp. 113–119, 2013.
View at: Publisher Site | Google Scholar
Y. Miyamoto, W. A. Kaysser, B. H. Rabin, A. Kawasaki, and R. G. Ford, Functionally Graded Materials: Design, Processing and Applications, Kluwer Academic Publishers, Boston, 1999.
R. M. Mahamood and E. T. Akinlabi, Functionally Graded Materials, Springer, Switzerland, 2017.
H. G. Fu, Q. Xiao, and J. D. Xing, “A study on the crack control of a high-speed steel roll fabricated by a centrifugal casting technique,” Materials Science and Engineering: A, vol. 474, no. 1-2, pp. 82–87, 2008.
View at: Publisher Site | Google Scholar
Z. P. Zhou, L. Huang, Y. J. Shang, Y. P. Li, L. Jiang, and Q. Lei, “Causes analysis on cracks in nickel-based single crystal superalloy fabricated by laser powder deposition additive manufacturing,” Materials & Design, vol. 160, pp. 1238–1249, 2018.
View at: Publisher Site | Google Scholar
Y. D. Li and K. Y. Lee, “Fracture analysis on the arc-shaped interface in a layered cylindrical piezoelectric sensor polarized along its axis,” Engineering Fracture Mechanics, vol. 76, no. 13, pp. 2065–2073, 2009.
View at: Publisher Site | Google Scholar
H. M. Zbib, J. P. Hirth, and I. Demir, “The stress intensity factor of cylindrical cracks,” International Journal of Engineering Science, vol. 33, no. 2, pp. 247–253, 1995.
View at: Publisher Site | Google Scholar
L. C. Guo, L. Z. Wu, T. Zeng, and L. Ma, “Mode I crack problem for a functionally graded orthotropic strip,” European Journal of Mechanics–A/Solids, vol. 23, no. 2, pp. 219–234, 2004.
View at: Publisher Site | Google Scholar
W. Z. Yang and Z. T. Chen, “Thermo-viscoelastic response of a cracked, functionally graded half-plane under a thermal shock,” Engineering Fracture Mechanics, vol. 206, pp. 267–277, 2019.
View at: Publisher Site | Google Scholar
G. Bao and L. Wang, “Multiple cracking in functionally graded ceramic/metal coatings,” International Journal of Solids and Structures, vol. 32, no. 19, pp. 2853–2871, 1995.
View at: Publisher Site | Google Scholar
H. J. Choi, T. E. Jin, and K. Y. Lee, “Collinear cracks in a layered half-plane with a graded nonhomogeneous interfacial zone – part II: thermal response,” International Journal of Fracture, vol. 94, no. 2, pp. 123–135, 1998.
View at: Google Scholar
H. J. Choi, K. Y. Lee, and T. E. Jin, “Collinear cracks in a layered half-plane with a graded nonhomogeneous interfacial zone – part I: mechanical response,” International Journal of Fracture, vol. 94, no. 2, pp. 103–122, 1998.
View at: Publisher Site | Google Scholar
Z. H. Jin and N. Noda, “Transient thermal stress intensity factors for a crack in a semi-infinite plate of a functionally gradient material,” International Journal of Solids and Structures, vol. 31, no. 2, pp. 203–218, 1994.
View at: Publisher Site | Google Scholar
F. Erdogan and B. H. Wu, “The surface crack problem for a plate with functionally graded properties,” Journal of Applied Mechanics, vol. 64, no. 3, pp. 449–456, 1997.
View at: Publisher Site | Google Scholar
P. Gu and R. J. Asaro, “Cracks in functionally graded materials,” International Journal of Solids and Structures, vol. 34, no. 1, pp. 1–17, 1997.
View at: Publisher Site | Google Scholar
Z. H. Jin, “Crack-tip singular fields in functionally graded materials,” in Encyclopedia of Thermal Stresses, R. B. Hetnarski, Ed., pp. 798–805, Springer, Dordrecht, 2014.
View at: Google Scholar
C. Y. Li and G. J. Weng, “Dynamic stress intensity factor of a cylindrical interface crack with a functionally graded interlayer,” Mechanics of Materials, vol. 33, no. 6, pp. 325–333, 2001.
View at: Publisher Site | Google Scholar
Z. H. Jin and G. H. Paulino, “A viscoelastic functionally graded strip containing a crack subjected to in-plane loading,” Engineering Fracture Mechanics, vol. 69, no. 14–16, pp. 1769–1790, 2002.
View at: Publisher Site | Google Scholar
W. J. Feng, R. K. L. Su, and Z. Q. Jiang, “Torsional impact response of a cylindrical interface crack between a functionally graded interlayer and a homogeneous cylinder,” Composite Structures, vol. 68, no. 2, pp. 203–209, 2005.
View at: Publisher Site | Google Scholar
Z. Q. Cheng, S. A. Meguid, and Z. Zhong, “Thermo-mechanical behavior of a viscoelastic FGMs coating containing an interface crack,” International Journal of Fracture, vol. 164, no. 1, pp. 15–29, 2010.
View at: Publisher Site | Google Scholar
P. P. Shi, “Cylindrical interface crack in a hollow layered functionally graded cylinder under static torsion,” Mechanics Based Design of Structures and Machines, vol. 44, no. 3, pp. 250–269, 2016.
View at: Publisher Site | Google Scholar
X. Y. Zhang, Y. J. Xie, and X. F. Li, “Transient thermoelastic response in a cracked strip of functionally graded materials via generalized fractional heat conduction,” Applied Mathematical Modelling, vol. 70, pp. 328–349, 2019.
View at: Publisher Site | Google Scholar
Z. T. Chen and K. Q. Hu, “Thermo-elastic analysis of a cracked half-plane under a thermal shock impact using the hyperbolic heat conduction theory,” Journal of Thermal Stresses, vol. 35, no. 4, pp. 342–362, 2012.
View at: Publisher Site | Google Scholar
Z. T. Chen and K. Q. Hu, “Thermoelastic analysis of a cracked substrate bonded to a coating using the hyperbolic heat conduction theory,” Journal of Thermal Stresses, vol. 37, no. 3, pp. 270–291, 2014.
View at: Publisher Site | Google Scholar
K. Q. Hu and Z. T. Chen, “Thermoelastic analysis of a partially insulated crack in a strip under thermal impact loading using the hyperbolic heat conduction theory,” International Journal of Engineering Science, vol. 51, pp. 144–160, 2012.
View at: Publisher Site | Google Scholar
B. L. Wang, “Transient thermal cracking associated with non-classical heat conduction in cylindrical coordinate system,” Acta Mechanica Sinica, vol. 29, no. 2, pp. 211–218, 2013.
View at: Publisher Site | Google Scholar
J. W. Fu, Z. T. Chen, L. F. Qian, and K. Q. Hu, “Transient thermoelastic analysis of a solid cylinder containing a circumferential crack using the C–V heat conduction model,” Journal of Thermal Stresses, vol. 37, no. 11, pp. 1324–1345, 2014.
View at: Publisher Site | Google Scholar
K. Q. Hu and Z. T. Chen, “Transient heat conduction analysis of a cracked half-plane using dual-phase-lag theory,” International Journal of Heat and Mass Transfer, vol. 62, pp. 445–451, 2013.
View at: Publisher Site | Google Scholar
X. Y. Zhang and X. F. Li, “Transient thermal stress intensity factors for a circumferential crack in a hollow cylinder based on generalized fractional heat conduction,” International Journal of Thermal Sciences, vol. 121, pp. 336–347, 2017.
View at: Publisher Site | Google Scholar
W. Z. Yang and Z. T. Chen, “Investigation of the thermal-elastic problem in cracked semi-infinite FGM under thermal shock using hyperbolic heat conduction theory,” Journal of Thermal Stresses, vol. 42, no. 8, pp. 993–1010, 2019.
View at: Publisher Site | Google Scholar
C. Cattaneo, “A form of heat conduction equation which eliminates the paradox of instantaneous propagation,” Compte Rendus, vol. 247, no. 4, pp. 431–433, 1958.
View at: Google Scholar
P. Vernotte, “Some possible complications in the phenomena of thermal conduction,” Compte Rendus, vol. 252, pp. 2190-2191, 1961.
View at: Google Scholar
P. J. Antaki, “New interpretation of non-Fourier heat conduction in processed meat,” Journal of Heat Transfer-Transactions of the ASME, vol. 127, no. 2, pp. 189–193, 2005.
View at: Publisher Site | Google Scholar
D. Y. Tzou, Macro- to Microscale Heat Transfer: The Lagging Behavior, John Wiley & Sons, United Kingdom, 2nd edition, 2015.
W. Roetzel, N. Putra, and S. K. Das, “Experiment and analysis for non-Fourier conduction in materials with non-homogeneous inner structure,” International Journal of Thermal Sciences, vol. 42, no. 6, pp. 541–552, 2003.
View at: Publisher Site | Google Scholar
F. Erdogan, G. D. Gupta, and T. S. Cook, “Numerical solution of singular integral equations,” in Mechanics of Fracture, Volume 1: Methods of Analysis and Solutions of Crack Problems, G. C. Sih, Ed., pp. 368–425, Noordhoff, Leyden, Netherlands, 1973.
View at: Google Scholar
J. W. Fu, Z. T. Chen, L. F. Qian, and Y. D. Xu, “Non-Fourier thermoelastic behavior of a hollow cylinder with an embedded or edge circumferential crack,” Engineering Fracture Mechanics, vol. 128, pp. 103–120, 2014.
View at: Publisher Site | Google Scholar
F. Durbin, “Numerical inversion of Laplace transforms: an efficient improvement to Dubner and Abate’s method,” The Computer Journal, vol. 17, no. 4, pp. 371–376, 1974.
View at: Publisher Site | Google Scholar
A. Akbarzadeh and Z. Chen, “Heat conduction in one-dimensional functionally graded media based on the dual-phase-lag theory,” Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, vol. 227, no. 4, pp. 744–759, 2013.
View at: Publisher Site | Google Scholar
L. C. Zhou, J. S. Lin, H. F. Lin, and G. H. Chen, “Electrical-thermal switching effect in high-density polyethylene/graphite nanosheets conducting composites,” Journal of Materials Science, vol. 43, no. 14, pp. 4886–4891, 2008.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Jiawei Fu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

493

Downloads

890

Citations