# 2D Quasi-Static Accurate Solutions for Isotropic Thermoelastic Materials with Applications

Authors: Wen-Ying Xiao; Jie Tong; Ying-Jie Liu; Jiang Su; Jian-Ping Li
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

Thermally induced stress is an important scientific problem in engineering applications. In this paper, an accurate and efficient method for the two-dimensional quasi-static thermal elastic problem is established to explore the thermal stress problem. First, the compact quasi-static two-dimensional general solution is derived in terms of simple potential functions. The general solution is simple in form and can be derived for arbitrary boundary problems subjected to a line heat load. This is completely new to the literature. Second, Green’s function solutions of an infinite plane under the line pulse heat source based on the general solutions are presented to analyze the thermal stress field. Lastly, numerical results are taken into account to study the temperature and stress field induced by the dynamic heat source load. The corresponding analysis can constitute to reveal the mechanism of thermal elastic problems and provide some guidance for experiments or engineering structural design in the future work.


---

## Body

## 1. Introduction

Thermal stress occurs when the temperature changes for structures [1]. The applications of thermoelastic mechanics have great significance in engineering. The stiffness and strength problems caused by thermal stress are indescribably important in some engineering structures which has been an important factor for the design of engineering devices in the field of aerospace, mechanical, civil, electronic engineering, and some microstructures [2–6]. Thermal stresses are generated during high-temperature curing of the material, which induces tensile and shear stress in the structure. These stresses may cause cracks during the heat conduction process, especially during the thermal shock process [5, 6]. The thermodynamic analysis is very important to reveal the interaction mechanism between the thermal source and structures [7–10].Many scholars have done excellent work in thermodynamics [11–14]. Goodier proposed the integration expression of thermoelastic equations in 1937, which laid the foundation for the study of thermodynamics [11]. Boley and Weiner explored the thermal stress problems in curve beams, rings, trusses, frames, and built-up structures [15]. In 1961, Mindlin derived the governing equations for thermal media from the fundamental laws of continuum mechanics and thermodynamics [16]. Nowinski proved the laws of uniqueness for the governing equations [17]. Hetnarski introduced the thermodynamic theory systematically [13]. Therefore, it is true to use these classical elasticity theories to solve many basic thermoelastic problems. In recent years, many research studies on the thermal elastic problem can be found in the literature. Numerical method and analytical method are used to study the thermoelastic problems. The finite element method (FEM) and boundary element method (BEM) are often preferred to analyze the thermoelastic problems. For example, a thermal finite element model for the multicomponent system of high-speed metal cutting was presented [18]. The thermal stress and local damage of the composite were simulated by finite element and discrete element coupling methods [19], but it is hard to obtain the accurate solutions of the thermal elastic structure, especially near the heat source. The reason is that extremely high-gradient coupling field near the heat source requires a high-density grid, which usually causes ill-conditioned equations in the calculation. Thus, the calculation cannot converge to obtain an exact solution of the coupling field naturally by the traditional software platforms.The main analytical methods are the integral transformation method, complex variable method, and classical general solution method. Integral transformation is a common method to solve differential equations. The Laplace and Hankel transform techniques have been used for the transversely isotropic thermoelastic thin plate [20]. The method of Fourier transformation is used to obtain an analytical expansion for the temperature and thermal stresses [21]. In the context of generalized theory of thermoelasticity proposed by Green and Naghdi, the thermal stress and temperature fields are obtained by using Fourier transforms [22, 23]. Lata and Parveen investigated the two-dimensional deformation of a homogeneous isotropic thermoelastic solid using the same method [24–26]. For quasi-static thermal conduction problems, three-dimensional transient heat conduction analysis was conducted by using the Laplace transformation method [27]. The thermoelasticity problem for the thick-walled cylinder subjected to transient thermal loading was analytically studied [28]. With the complexity of the structure, the integral matrix becomes larger and less efficient. And the integral transformation condition is very harsh in many cases which is difficult to get accurate results by inverse transformation. For using the complex variable method, Hasebe et al. presented thermal stresses and temperature distribution for an infinite plate with circular holes and cracks under uniform heat flux in the arbitrary direction [29]. Gao et al. presented an exact solution for the problem of an elliptical hole or a crack in a thermo-piezoelectric material [30]. Mohammad studied thermal stress analysis of the orthotropic plate containing a rectangular hole using the complex variable method [31], but complex function methods are usually only applicable to solve plane problems, which is very difficult to obtain a simple form.The general solution method is a basic method to solve mechanical problems. It directly starts from the governing equations of mechanical models to find the general solution of partial differential equations and then combines the boundary conditions of specific problems to find the specific solution. The general solution under point loading is called Green’s function solution. Green’s function is an efficient tool for refined analysis of thermal stress of solid structures which can provide necessary preconditions to many high-precision numerical methods, such as boundary element method [18, 19]. For isotropic materials, some scholars gave the basic solution of the famous closed form [32–34]. Hou’s team did some work on Green’s function of thermal elastic structure based on the general solutions [35–38], but the thermal elastic analysis of the isotropic structure under dynamic heat load has not been seen. In this literature, we would introduce how to get a compact and complete thermoelastic general solution for quasi-static heat conduction problems.In this background, a harmonic function and a function which satisfies the quasi-static heat conduction equation for Green’s functions under the quasi-static line heat source will be constructed in this study. Firstly, we use the differential operator theory to derive three general solutions. Secondly, by virtue of Almansi’s theorem, the three general solutions can be transformed to the general solution which is expressed in terms of two harmonic functions and a function which satisfies the quasi-static heat conduction equation, respectively. Thirdly, the integral general solution expressed in one harmonic function and one function which satisfies the quasi-static heat conduction equation is obtained. As an example, the suitable functions for a quasi-static line heat source in the interior of the infinite steel plane are constructed, and the corresponding Green’s solutions are presented by virtue of the obtained general solutions. In this literature, all these harmonic functions are in the form of elementary functions, and all the thermal stress components are presented in a closed form at an arbitrary point in a two-phase infinite body. The fundamental solutions for the arbitrary distributed heat source are obtained by numerical integration using Green’s function under the point heat source. The method proposed in this study presents an accurate and efficient tool to analyze the thermal elastic problems, which have great engineering application value.


## 2. Two-Dimensional Governing Equations for the Isotropic Thermoelastic Material

If all components are independent of coordinatey in Cartesian coordinates x,y,z, one will obtain the so-called two-dimensional problem or plane problem in two-dimensional Cartesian coordinates x,z. In this case, the corresponding constitutive equations of isotropic thermoelastic materials are in the form of(1)σx=λ+2G∂u∂x+λ∂ω∂z−αθ,σz=λ∂u∂x+λ+2G∂ω∂z−αθ,τzx=G∂u∂z+∂ω∂x,where u and ω are the components of the mechanical displacement in x- and z-directions, respectively, σx,σz, and τzx are the components of stress, θ is the temperature increment, and α is the thermal modulus. λ and G are the Lame constant and shear modulus, which are in the following forms for the plane stress problem and plane strain problem, respectively.For the plane stress problem,(2a)λ=νE1−ν2,G=E21+ν.For the plane strain problem,(2b)λ=νE1+ν1−2ν,G=E21+ν,where E and v are Young’s modulus and Poisson’s ratio of the material.In the absence of body forces, the mechanical equilibrium equations are(3)∂σx∂x+∂τzx∂z=0,∂τzx∂x+∂σz∂z=0.The two-dimensional quasi-static heat conduction equation is(4)∂2∂x2+∂2∂z2θ=1κ∂θ∂t,where κ is the thermal diffusivity of heat conduction:(5)κ=βρce,where β,ρ,andce are the coefficient of thermal conductivity, density, and specific heat of the material, respectively.


## 3. Two-Dimensional Quasi-Static General Solutions for the Isotropic Thermoelastic Material

In this section, the quasi-static general solutions are derived using the governing equations1–(5) for the isotropic thermal elastic materials by virtue of the differential operator theory. Firstly, the equilibrium equations by substituting equation (1) into (3) yield(6a)λ+2G∂2∂x2+G∂2∂z2u+λ+G∂2ω∂x∂z−α∂θ∂x=0,(6b)λ+G∂2u∂x∂z+G∂2∂x2+λ+2G∂2∂z2ω−α∂θ∂z=0.A combination of the results of equations (4), (6a), and (6b) gives(7)Duωθ=0,where D is the differential operator matrix defined by(8)D=λ+2G∂2∂x2+G∂2∂z2λ+G∂2∂x∂z−α∂∂xλ+G∂2∂x∂zG∂2∂x2+λ+2G∂2∂z2−α∂∂z00∂2∂x2+∂2∂z2−1κ∂∂t.Based on the differential operator theory, three general solutions of homogeneous differential equation (7) can be obtained as follows:(9)u=Ai1f,ω=Ai2f,θ=Ai3fi=1,2,3.Aiji,j=1,2,3 are the algebraic cofactors of operator matrix D. i,j=1,2,3 represent row i and column j of operator matrix D, and i=1,2,3 indicates three cases in which row i=1 or row i=2 or row i=3. The function f should satisfy the following homogeneous equation:(10)Df=0.The determinant of derivative operator matrixD is(11)D=Gλ+2G∂2∂x2+∂2∂z22∂2∂x2+∂2∂z2−1κ∂∂t,where Aiji,j=1,2,3 in equation (9) are the algebraic complements of matrix D which are as follows:(12a)A11=G∂2∂x2+λ+2G∂2∂z2∂2∂x2+∂2∂z2−1κ∂∂t,A12=−λ+G∂2∂x∂z∂2∂x2+∂2∂z2−1κ∂∂t,A13=0,(12b)A21=−λ+G∂2∂x∂z∂2∂x2+∂2∂z2−1κ∂∂t,A22=λ+2G∂2∂x2+G∂2∂z2∂2∂x2+∂2∂z2−1κ∂∂t,A23=0,(12c)A31=αG∂∂x∂2∂x2+∂2∂z2,A32=αG∂∂z∂2∂x2+∂2∂z2,A33=Gλ+2G∂2∂x2+∂2∂z22.Wheni=1,2 and Ai3=0, we can obtain the general solutions in the isothermal case. The solution is the case i=3 considering the variable temperature. In this paper, the case i=3 is analyzed.By virtue of Almansi’s theorem, the functionf in equation (10) can be expressed in the following form:(13)f=f1+zf2+f3,where fj satisfy the following harmonic equation:(14a)∂2∂x2+∂2∂z2fj=0,j=1,2,(14b)∂2∂x2+∂2∂z2−1κ∂∂tf3=0.For the convenience of operation, let(15)F1=2∂f2∂z,F2=∂2∂x2+∂2∂z2f3=1κ∂f3∂t,and then(16)∂2∂x2+∂2∂z2f=∂2∂x2+∂2∂z2f1+zf2+f3=F1+F2.Based on equations (14a)–(15), one can find that the functions Fjj=1,2 in equation (16) still satisfy the following harmonic equation:(17)∂2∂x2+∂2∂z2F1=0,∂2∂x2+∂2∂z2−1κ∂∂tF2=0.Two general solutions are obtained in whichθ=0 when i=1,2. These solutions are in fact identical to those without the thermal effect which are not discussed here. Therefore, the solution should be taken in equation (12c) in which i=3. The following general solution is obtained:(18)u=αG∂F1∂x+∂F2∂x,ω=αG∂F1∂z+∂F2∂z,θ=Gλ+2G∂2∂x2+∂2∂z2F2=Gλ+2Gκ∂F2∂t.Let(19)ψ1=2αG2F1,ψ2=2αG2F2.Based on equations (17) and (19), one can find that the functions ψjj=1,2 still satisfy the following harmonic equation:(20)∂2∂x2+∂2∂z2ψ1=0,∂2∂x2+∂2∂z2−1κ∂∂tψ2=0.In combination with equations (18) and (20), the solutions are obtained:(21)2Gu=∂ψ1∂x+∂ψ2∂x,2Gω=∂ψ1∂z+∂ψ2∂z,2GΓθ=∂ψ2∂t,where(22)Γ=λ+2Gακ=1−νEακ1+ν1−2ν.Substituting general solution equation (21) into constitutive equation (1), the stress general solution can be simplified as follows:(23)σx=∂2ψ1∂x2+∂2ψ2∂x2−1κ∂ψ2∂t,σz=∂2ψ1∂z2+∂2ψ2∂z2−1κ∂ψ2∂t,τzx=∂2ψ1∂x∂z+∂2ψ2∂x∂z.The two-dimensional quasi-static general solutions for isotropic thermoelastic materials are obtained in equations (21) and (23). These general solutions are expressed in terms of two functions ψjj=1,2 in equation (20). The suitable functions ψjj=1,2 with undetermined constants should be constructed according to the configuration of the structures and the form of loadings. The undetermined constants will be determined by the corresponding boundary conditions, compatible conditions, and other requests. In addition, it is important that the function ψ1 satisfies the harmonic equation and the function ψ2 satisfies the quasi-static heat conduction equation. As an example, the suitable functions ψjj=1,2 under a dynamic line heat source will be constructed based on the obtained general solutions. For the arbitrary heat sources in actual engineering problems, their stress fields can be obtained based on the fundamental solutions for the line heat source using the superposition principle.


## 4. Fundamental Solution for an Infinite Isotropic Thermoelastic Plane under a Pulse Line Heat Source

Consider an infinite isotropic thermoelastic thin plane in two-dimensional Cartesian coordinatesx,z with the principle material direction parallel to the x- and z-axis (Figure 1). A pulse line heat source Hx,z,t=Hδt is applied at the origin 0,0 on the xoz plane. Based on general solutions (21) and (23), the coupled thermoelastic field in the infinite thermoelastic plane is derived in this section.Figure 1
An infinite isotropic thermoelastic plane applied by a pulse line heat sourceH.The following functions are introduced by the trial and error method:(24)ψ1=−2A1lnr,ψ2=A2Ei−r2τ,where Ajj=1,2 are arbitrary constants to be determined and(25)Eix=∫−∞xeηηdη,τ=4κt,r=x2+z2.Substituting function (24) into general solutions (21) and (23) yields(26a)2Gux=−A12xr2+A22xr2e−r2/τ,2Guz=−A12zr2+A22zr2e−r2/τ,2GΓθ=−4κA21τe−r2/τ,(26b)σx=−2A11r2−2x2r4+2A21r2−2x2r4−2x2r2τ+2τe−r2/τ,σz=−2A11r2−2z2r4+2A21r2−2z2r4−2z2r2τ+2τe−r2/τ,τzx=4xzA11r4−4xzA21r4+1r2τe−r2/τ.When the mechanical equilibrium is considered (Figure1), the following equation can be obtained:(27)∫−∞+∞∫−∞+∞ρceθdxdz=Hδxδzδt.Substituting equation (26a) into (27) yields(28)−ρce2A2κΓG1τ∫−∞+∞∫−∞+∞exp−r2τdxdz=Hδxδzδt.One useful integral is listed as follows:(29)1τ∫−∞+∞∫−∞+∞exp−r2τdxdz=1τ∫−∞+∞∫−∞+∞exp−x2+z2τdxdz=∫−∞+∞∫−∞+∞exp−x2+z2τdxτdzτ=π222π2∫−∞+∞exp−x2τdxτ∫−∞+∞exp−z2τdzτ=πerf∞2=π,where erfx is the error function. It is defined as follows:(30)erfx=2π∫0xexp−η2dη.Substituting equation (29) into (28) yields(31)A2=−HG2βπΓ.If the time tends to infinity, all the components are given as follows:(32a)2Gux=−A12xr2+A22xr2,2Guz=−A12zr2+A22zr2,θ=0,(32b)σx=−2A11r2−2x2r4+2A21r2−2x2r4,σz=−2A11r2−2z2r4+2A21r2−2z2r4,τzx=4xzA11r4−4xzA21r4.The shear stress would be zero at infinity, and the following equation can be obtained:(33)A1−A2=0.Substituting equation (31) into (33) yields(34)A1=A2=−HG2βπΓ.Hence, Green’s solution for a pulse line heat source in an infinite isotropic thermoelastic plane is determined by equation (34).


## 5. Results and Discussion

Based on the solutions derived above, we now present some numerical results. The following nondimensional components are used in the figures:(35)λ1=GE=ν1+ν1−2ν,G1=λE=121+ν,α1=αEα0,E1=EE=1,Γ1=2λ1+2G1α1,β1=ββ=1,τ′=τr02,κ′=κ4κ=14,ξ=xr0,ζ=zr0,uk=uir0α0T0,ϑ=θT0,σk=σiEα0T0,τkl=τijEα0T0i,j=x,z;k,l=ξ,ζ,where r0 is a nonzero dimension and α0 and T0 are the thermal expansion coefficient and reference temperature, respectively. Equations (32a) and (32b) should be rewritten in the following nondimensional form:(36)A2=−2δGπΓ1,where δ is a nondimensional heat source given as follows:(37)δ=Hr02T0ρce.Here, letδ=1 act at the point 0,0, so the heat source H=r02T0ρce. The material properties for steel are listed in Table 1.Table 1
Material properties.
MaterialsYoung’s modulusE109Nm−2Poisson’s ratioνDensityρ103kg⋅m−3Specific heatce103J⋅kg−1°CThermal expansion coefficientα010−6K−1Steel210.00.297.80.4612.0


### 5.1. Thermoelastic Field Contours under the Pulse Line Heat Source

Temperature and thermal-stress distribution contours of steel are given in Figures2–5 when a pulse line heat source acts on the center of the plane. On the basis of the contours plotted, one can find all components are singular at the origin in which the line heat source is located. The components change gradually with time τ and then tend to be stable which do not change over time. It can be found that there is a great thermal impact near the heat source at the moment of the heat source action. So, the stresses are relatively large which change quite dramatically. With time going on, the stress values decrease gradually with the diffusion of the heat source. We can clearly see the trend of the temperature field and stress field, which helps engineers to predict the engineering structure stress state under the action of thermal load in advance.(1)
Figure2 shows that the contours of temperature increment ϑθ in the infinite plane are all normal circles, which is because the thermal diffusivity κ along the x- and z-axis is equal for steel, respectively. In the initial stage of the heat source acting on the plane, the heat source is confined in a small region. And the temperature increment is very small which nearly approaches to be zero, far away from the heat source region. This also indicates that the temperature gradient is very large under the pulse heat source. A large stress field is easily generated near the operating point of the heat source which could cause structural damage. With the passage of time, the gradients of temperature increment decrease gradually and tend to be in a stable state which clearly reveals the change process of the temperature field.(2)
It can be seen from Figures3 and 4 that the radial component σξσx and axial component σςσz of thermal stress have positive and negative values which vary with the heat source acting time. The gradients of negative contours for the axial stress σςσz are much larger than those of positive contours at the time τ′=1. As time goes on, the zero contour moves outward, and the whole stress fields would be compressive stress. In the initial stage of the heat acting, the stresses rise and drop sharply because the heat energy is mainly concentrated near the heat source region.(3)
Figure5 shows that there is a zero common tangent for shear stress τςξτzx in an infinite steel plane. And the contours of shear stress τςξτzx are symmetric with the zero contours. The shear stress τςξτzx changes significantly between the interfaces of the zero contour at the beginning which means the stress failure of the materials is possible to occur at the time.Figure 2
Temperature increment contourϑ×102 under a line heat source δ=1.
(a)(b)(c)(d)(e)(f)Figure 3
Stress contourσξ×102 under a line heat source δ=1.
(a)(b)(c)(d)(e)(f)Figure 4
Stress contourσς×102 under a line heat source δ=1.
(a)(b)(c)(d)(e)(f)Figure 5
Stress contourτςξ×102 under a line heat source δ=1.
(a)(b)(c)(d)(e)(f)


### 5.2. Thermoelastic Field Contours under Arbitrary Distributed Pulse Heat Source Load

In this section, we calculate the temperature and thermal stress distributions under different pulse-type heat source loads. By virtue of the superposition principle and numerical integration method, the coupled field under arbitrary distributed pulse heat source can be obtained easily. In this section, we consider one kind of heat source loaded on the plane.For example, the dimensionless distributed pulse heat source loaded on the plane is in the form of(38)fξ,ς=ξ2+ς2,−1≤ξ≤1,−1≤ς≤1.The region of the distributed pulse heat source is(39)Aξ,ς:ξ2+ς2=1,−1≤ξ≤1,−1≤ς≤1.The contours of the stress components and temperature component for the distributed heat source are shown in Figures6–9. Figures 6–9 show that the temperature and stress components tend to be stable gradually which is the same as the variation trend of that in Section 5.1.Figure 6
Temperature increment contourϑ×102 under the distributed pulse heat source.
(a)(b)(c)Figure 7
Stress contourσξ×102 under the distributed pulse heat source.
(a)(b)(c)Figure 8
Stress contourσζ×102 under the distributed pulse heat source.
(a)(b)(c)Figure 9
Stress contourτςξ×102 under the distributed pulse heat source.
(a)(b)(c)


### 5.3. Variation Characteristics of the Thermal Stress Field Components

The distributions of temperature and thermal stress at different positions of the plane with different times are given in Figures10–13. Figure 10 depicts the variation of temperature distribution with t. It indicates the values of temperature increment ϑ have a sharp rise and then have a rapid decline when 0<t<20. It is noticed that the temperature values slightly decrease when t≥50 and then tend to be stable. Figures 11 and 12 exhibit the variations of radial and axial stress components with time t. The thermal stresses show compressive stress when they are stable. Figure 13 shows the variation of the shear stress with time t. It is noticed that the values of the shear stress τξζ become stable when t>20. When 0<t<20, the shear stress values are negative and increase rapidly.Figure 10
Temperature distributions at different locations.Figure 11
Radial component of thermal stress distributions at different locations.Figure 12
Axial component of thermal stress distributions at different locations.Figure 13
Shear stress component of thermal stress distributions at different locations.


## 5.1. Thermoelastic Field Contours under the Pulse Line Heat Source

Temperature and thermal-stress distribution contours of steel are given in Figures2–5 when a pulse line heat source acts on the center of the plane. On the basis of the contours plotted, one can find all components are singular at the origin in which the line heat source is located. The components change gradually with time τ and then tend to be stable which do not change over time. It can be found that there is a great thermal impact near the heat source at the moment of the heat source action. So, the stresses are relatively large which change quite dramatically. With time going on, the stress values decrease gradually with the diffusion of the heat source. We can clearly see the trend of the temperature field and stress field, which helps engineers to predict the engineering structure stress state under the action of thermal load in advance.(1)
Figure2 shows that the contours of temperature increment ϑθ in the infinite plane are all normal circles, which is because the thermal diffusivity κ along the x- and z-axis is equal for steel, respectively. In the initial stage of the heat source acting on the plane, the heat source is confined in a small region. And the temperature increment is very small which nearly approaches to be zero, far away from the heat source region. This also indicates that the temperature gradient is very large under the pulse heat source. A large stress field is easily generated near the operating point of the heat source which could cause structural damage. With the passage of time, the gradients of temperature increment decrease gradually and tend to be in a stable state which clearly reveals the change process of the temperature field.(2)
It can be seen from Figures3 and 4 that the radial component σξσx and axial component σςσz of thermal stress have positive and negative values which vary with the heat source acting time. The gradients of negative contours for the axial stress σςσz are much larger than those of positive contours at the time τ′=1. As time goes on, the zero contour moves outward, and the whole stress fields would be compressive stress. In the initial stage of the heat acting, the stresses rise and drop sharply because the heat energy is mainly concentrated near the heat source region.(3)
Figure5 shows that there is a zero common tangent for shear stress τςξτzx in an infinite steel plane. And the contours of shear stress τςξτzx are symmetric with the zero contours. The shear stress τςξτzx changes significantly between the interfaces of the zero contour at the beginning which means the stress failure of the materials is possible to occur at the time.Figure 2
Temperature increment contourϑ×102 under a line heat source δ=1.
(a)(b)(c)(d)(e)(f)Figure 3
Stress contourσξ×102 under a line heat source δ=1.
(a)(b)(c)(d)(e)(f)Figure 4
Stress contourσς×102 under a line heat source δ=1.
(a)(b)(c)(d)(e)(f)Figure 5
Stress contourτςξ×102 under a line heat source δ=1.
(a)(b)(c)(d)(e)(f)


## 5.2. Thermoelastic Field Contours under Arbitrary Distributed Pulse Heat Source Load

In this section, we calculate the temperature and thermal stress distributions under different pulse-type heat source loads. By virtue of the superposition principle and numerical integration method, the coupled field under arbitrary distributed pulse heat source can be obtained easily. In this section, we consider one kind of heat source loaded on the plane.For example, the dimensionless distributed pulse heat source loaded on the plane is in the form of(38)fξ,ς=ξ2+ς2,−1≤ξ≤1,−1≤ς≤1.The region of the distributed pulse heat source is(39)Aξ,ς:ξ2+ς2=1,−1≤ξ≤1,−1≤ς≤1.The contours of the stress components and temperature component for the distributed heat source are shown in Figures6–9. Figures 6–9 show that the temperature and stress components tend to be stable gradually which is the same as the variation trend of that in Section 5.1.Figure 6
Temperature increment contourϑ×102 under the distributed pulse heat source.
(a)(b)(c)Figure 7
Stress contourσξ×102 under the distributed pulse heat source.
(a)(b)(c)Figure 8
Stress contourσζ×102 under the distributed pulse heat source.
(a)(b)(c)Figure 9
Stress contourτςξ×102 under the distributed pulse heat source.
(a)(b)(c)


## 5.3. Variation Characteristics of the Thermal Stress Field Components

The distributions of temperature and thermal stress at different positions of the plane with different times are given in Figures10–13. Figure 10 depicts the variation of temperature distribution with t. It indicates the values of temperature increment ϑ have a sharp rise and then have a rapid decline when 0<t<20. It is noticed that the temperature values slightly decrease when t≥50 and then tend to be stable. Figures 11 and 12 exhibit the variations of radial and axial stress components with time t. The thermal stresses show compressive stress when they are stable. Figure 13 shows the variation of the shear stress with time t. It is noticed that the values of the shear stress τξζ become stable when t>20. When 0<t<20, the shear stress values are negative and increase rapidly.Figure 10
Temperature distributions at different locations.Figure 11
Radial component of thermal stress distributions at different locations.Figure 12
Axial component of thermal stress distributions at different locations.Figure 13
Shear stress component of thermal stress distributions at different locations.


## 6. Conclusions

In this paper, the compact quasi-static two-dimensional general solution for isotropic thermoelastic materials is presented by virtue of differential operator theory and Almansi’s theorem. The general solutions are expressed using simple explicit functions. Based on the general solutions, Green’s functions are obtained for a line heat source. Also, all Green’s components are expressed in simple functions, which are very convenient to use. The solution in this paper is obtained completely from the governing equation. No simplifying assumptions have been made. As an example, Green’s solutions for a pulse line heat source in the interior of an infinite steel plane are presented by virtue of the obtained general solution. Some numerical analyses are taken into account to study the distribution of the thermoelastic field under the action of the dynamic heat source. The numerical results of this study can give some theoretical basis for revealing the mechanism of the thermal elastic problem and provide some guidance for experiments or engineering structural design in the future work.



---

*Source: 8825226-2021-04-28.xml*