MPEMathematical Problems in Engineering1563-51471024-123XHindawi Publishing Corporation28438010.1155/2009/284380284380Research ArticleNew Poisson's Type Integral Formula for Thermoelastic Half-SpaceSeremetVictor1BonnetGuy2SpeianuTatiana3RagabSaad A.1Department of Mathematics and EngineeringAgrarian State University of MoldovaMircesti Street 442049 ChisinauMoldovauasm.md2Laboratoire de Modélisation et Simulation Multi-Echelle FRE 3160 CNRSUniversité Paris Est 5 boulevard Descartes77454 Marne la Vallée CedexFranceuniv-paris-est.fr3Department of MathematicsAgrarian State University of MoldovaMircesti Street 44 2049 ChisinauMoldovauasm.md20090707200920090501200923042009070520092009Copyright © 2009This 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.

A new Green's function and a new Poisson's type integral formula for a boundary value problem (BVP) in thermoelasticity for a half-space with mixed boundary conditions are derived. The thermoelastic displacements are generated by a heat source, applied in the inner points of the half-space and by temperature, and prescribed on its boundary. All results are obtained in closed forms that are formulated in a special theorem. A closed form solution for a particular BVP of thermoelasticity for a half-space also is included. The main difficulties to obtain these results are in deriving of functions of influence of a unit concentrated force onto elastic volume dilatation Θ(k) and, also, in calculating of a volume integral of the product of function Θ(k) and Green's function in heat conduction. Using the proposed approach, it is possible to extend the obtained results not only for any canonical Cartesian domain, but also for any orthogonal one.

1. Introduction

The main objective of this paper is to prove a theorem (Section 2.2) about deriving a Poisson’s type integral formula for a thermoelastic half-space with the homogeneous locally mixed mechanical boundary conditions and with the nonhomogeneous Dirichlet’s boundary condition for temperature. To prove this theorem we need some equations for Green’s integral formula in stationary thermoelasticity (Section 1.1) and thermoelastic influence functions (Sections 1.2 and 1.3) for a solid body, suggested and published by the first-named author earlier. Needed special cases of these equations for a half-space are given in Section 2.1. An example of application of the obtained new Poisson’s type integral formula is given in Section 3.

1.1. General Green’s Integral Formula in Stationary Thermoelasticity

The Green’s function plays the leading role in finding the solutions in integrals for boundary value problems (BVPs) in different fields of mathematical physics. The theory of thermoelasticity, which is a synthesis of the theory of heat conduction and elasticity theory, is one of such fields. By now, a number of theories of thermoelasticity have been developed and described in classical scientific literature . The most developed theory, which is widely used in practical calculations, is the theory of thermal stresses, that is, the theory of uncoupled thermoelasticity, where the temperature field does not depend on the field of elastic displacements. According to this theory the formulation of a BVP consists in nonhomogeneous Lame’s equations μ2uk(ξ)+(λ+μ)θ,k(ξ)-γT,k(ξ)=0 with respective mechanical boundary conditions, where λ, μ are Lame’s constants of elasticity; γ=αt(2μ+3λ) is the thermoelastic constant; αt is the coefficient of the linear thermal expansion. The temperature field T in (1.1) has to be determined from the BVP in heat conduction that consists in the Poisson’s equation 2T(ξ)=-a-1F(ξ)ξV, with the respective boundary conditions for heat actions: temperature T(y), heat flux a(T(y)/ny), or certain law of heat exchange between exterior medium and surface of the body αT(y)+aT(y)/ny, where F is the heat source; a is thermal conductivity; α is the coefficient of convective heat conductivity. To solve the BVP of stationary thermoelasticity in (1.1) and (1.2) using traditional methods, at the first stage we need to solve the BVP of heat conduction in (1.2) with the given boundary conditions and to find the temperature field. At the second stage we need to solve BVP of thermoelasticity in (1.1) with the already known temperature field and with the given mechanical boundary conditions. On the base of the influence functions introduced in the works  the following new general Green’s type integral formula for determining the fields of displacements, described by BVP in (1.1)-(1.2), was suggested: uk(ξ)=a-1VF(x)Uk(x,ξ)dV(x)-ΓDT(y)Uk(y,ξ)nydΓD(y)+ΓN[T(y)ny]Uk(y,ξ)dΓN(y)+ΓM[T(y)+(aα-1)T(y)ny]Uk(y,ξ)dΓM(y), where ΓD, ΓN, and ΓM denote the surfaces on which the boundary conditions of Dirichlet’s, Neumann’s, or mixed type are prescribed. One of the advantages of this formula is that the desired thermoelastic displacements uk are determined in the integral form directly via the prescribed inner heat source and other thermal data, given on the boundary.

1.2. Main Thermoelastic Influence Functions

Thus, the introduced functions of influence of a unit heat source on thermoelastic displacements Uk(x,ξ) in (1.3) are determined by the following general integral formula : Uk(x,ξ)=γVG(x,z)Θ(k)(z,ξ)dV(z),x,z,ξV. This formula represents the Mayzel’s formula, when the temperature field (the Green’s function G in heat conduction problem) is generated by an inner unit point heat source. The volume dilatation Θ(k) in (1.4) has to be determined from the following Lame set of equations: μ2Ui(k)(x,ξ)+(λ+μ)Θ,i(k)(x,ξ)=-δikδ(x-ξ),i,k=1,2,3, with the respective homogeneous mechanical boundary conditions. In (1.5) δ(x-ξ) is the Dirac’s function and δik is the Kronecker’s symbol. The Green’s function G in (1.4) has to be determined from the following equation in heat conduction: x2G(x,ξ)=-δ(x-ξ),x,ξV, with the respective homogeneous boundary conditions for heat actions. A generalization of influence function and Green’s integral formula on the BVP in the classical theory of stationary thermoelasticity  was given for the first time in the papers . The considered influence functions Uk(x,ξ) have the physical sense as displacements at an inner point of observation x(x1,x2,x3), generated by a unit heat source, applied at an inner point ξ(ξ1,ξ2,ξ3), and described by the Dirac’s function. According to (1.4) they are determined by a convolution over the body V of two influence functions. The first influence function is the Green’s function G for the BVP in heat conduction. The other functions Θ(k) are functions of influence of inner concentrated body forces on elastic volume dilatation.

Finally, the influence functions Uk(x,ξ) are functions of double influence , which take in consideration both physical phenomena (heat conduction and elasticity) in a solid body.

Over the coordinates of the point of observation x(x1,x2,x3) for thermoelastic displacements, they satisfy the equations of the BVP for determining Green's functions in the theory of heat conduction (1.6). The only difference is that the unit heat source is replaced by the function of influence of the unit concentrated forces on the volume dilatation x2Uk(x,ξ)=-γ  Θ(k)(x,ξ)with the respective homogeneous boundary conditions for heat actions.

Over the coordinates of the point of application ξ(ξ1,ξ2,ξ3) of the unit point heat source, they satisfy the BVP for determining components of the Green's matrix (1.5). The only difference is that the unit concentrated body forces are replaced with the derivatives of Green's functions of the heat conduction problem μξ2Uk(x,ξ)+(λ+μ)Θ,ξk(x,ξ)-γG,ξk(x,ξ)=0with the respective homogeneous mechanical boundary conditions.

1.3. Other Thermoelastic Influence Functions

Take note that all influence functions in (1.3) are determined on the boundary independently or using the respective limits from the main influence functions Uk(x,ξ):

the formula for influence functions of a unit point heat flux a(T(y)/ny)=δ(x-y) on the surface ΓN on thermoelastic displacements is Uk(y,ξ)=limxyγVG(x,z)Θ(q)(z,ξ)dV(z)=limxyUk(x,ξ),x,ξ,zV,yΓN,

the formula for influence functions of a unit point temperature T(y)=-δ(x-y) on the surface ΓD on thermoelastic displacements is, Uk(y,ξ)ny=limxyγVG(x,z)nxΘ(q)(z,ξ)dV(z),x,ξ,zV=limxyUk(x,ξ)nx,yΓD,

the formula for influence functions of unit point heat exchange of the body with exterior medium αT(y)+aT(y)/ny=δ(x-y) through the surface ΓM on thermoelastic displacements is, Uk(y,ξ)=limxyγVG(x,z)Θ(k)(z,ξ)dV(z),x,ξ,zV=limxyUk(x,ξ),yΓM.

The formula in (1.3) can be treated also as a generalization of Mayzel’s formula  for those cases when the temperature field satisfies the BVP of heat conduction. Temperature field in this case is caused by the inner heat source and by the prescribed on the boundary temperature, heat flux, or certain law of heat exchange between exterior medium and surface of the body.

The advantage of the proposed integral formula in (1.3) is that it allows us to unite the two-staged process of solving the BVP in the theory of thermoelasticity (the first stage comprises finding temperature fields, and the second stage comprises finding thermoelastic displacements) in one stage. Also, the advantage of the integral formula in (1.3) in comparison with the well-known Mayzel’s integral formula is that the thermoelastic displacements are determined directly via given heat actions. Besides, for any concrete type of BVP we can obtain all possible solutions for different laws describing the above-mentioned heat actions. The main difficulties for practical realization of the integral formula in (1.3) and (1.4) are to derive the functions of influence of a unit concentrated force on elastic volume dilatation Θ(k) and of the Green’s functions in heat conduction G. In addition we need to compute some volume integrals of the product of the above-mentioned functions. These difficulties, especially deriving functions Θ(k), were overcome successfully for Cartesian canonical domains . For cylindrical and spherical domains only general integral representations for Θ(k) and Green’s matrices were proposed [16, 17].

2. Elastic Response of a Half-Space to a Unit Point Heat Source

In this section we give a theorem for determining the thermoelastic displacements for a half-space in the form of volume and surface integrals, which is a particular case of the general integral formula in (1.3). To do this, first, on the basis of the theory described above we construct the functions of influence of the inner unit point heat source on the thermoelastic displacements Uk(x,ξ). At the second step we have to calculate (on the basis of the main influence functions Uk(x,ξ)) the other influence functions Uk(y,ξ)/ny and to write the Poisson’s type integral formulas for respective BVP of thermoelasticity for half-space. At the last step it is necessary to prove that the obtained influence functions and Poisson’s type integral formulas satisfy the respective BVP.

2.1. Main Influence Functions and Green’s Integral Formula for a Thermoelastic Half-space

To obtain the functions of influence of an inner unit point heat source F(ξ)=δ(x-ξ) on thermoelastic displacements for a half-space we will use the general formula in (1.4) rewritten as follows: Uk(x,ξ)=γ0-G(x,z)  Θ(k)(z,ξ)dz1dz2dz3,|Uk(x,ξ)|<. Before starting to formulate the theorem we have to give the following Poisson’s type integral formula for half-space V(0x1<;    -<x2,x3<): uk(ξ)=a-10-F(x)Uk(x,ξ)dx1dx2dx3--T(y)Uk(y,ξ)ny1dy2dy3,|uk(ξ)|<, where the thermoelastic displacements uk(ξ) are generated by the inner heat source F(x) and given temperature T(y) on the boundary plane Γ(y1=0;    -<y2,y3<). The formula in (2.2) is obtained from the suggested general integral formula from (1.3), wherein it should be taken into account that on the all boundary plane Γonly Dirichlet’s condition (temperature) is given. Note that when we check if the boundary conditions for thermoelastic stresses σij(x,ξ) with respect to coordinates of the point ξ(ξ1,ξ2,ξ3) are satisfied, we have to use the Duhamel-Neumann law σij=μ(Ui,j+Uj,i)+δij(λUj,j-γG),i,j=1,2,3.

2.2. A Theorem about Deriving a Poisson’s Integral Formula for Thermoelastic Half-SpaceTheorem 2.1.

Let the field of displacements uk(ξ) at inner points ξ(ξ1,ξ2,ξ3) of the elastic half-space V(0x1<;-<x2,x3<) be determined by nonhomogeneous Lame equations μ2uk(ξ)+(λ+μ)θ,k(ξ)-γT,k(ξ)=0, and at the points y(0,y2,y3) of its boundary plane Γ(y1=0;-<y2,y3<) the following homogeneous mixed mechanical conditions are given: σ11(ξ1=0,ξ2=y2,,ξ3=y3)=0,u2(ξ1=0,ξ2=y2,,ξ3=y3)=0,u3(ξ1=0,ξ2=y2,,ξ3=y3)=0, where σ11 are normal stresses which are determined by the Duhamel-Neumann law (2.3).

Let also the temperature field T(ξ) in (2.4), generated by the inner heat source F(ξ) and temperature T(y)Γ (Dirichlet’s boundary condition), satisfy the following BVP in heat conduction: 2T(ξ)=-a-1F(ξ),ξV,T(ξ1=0,ξ2=y2,,ξ3=y3)=T(y),y(0,y2,y3)Γ. If the inner heat source and boundary temperature satisfies the conditions 0-|F(x)|dx1dx2dx3<,-|T(y)|dy2dy3<, then the solution of this BVP in (2.4)–(2.7) of thermoelasticity for searched displacements uk(ξ) for half-space exists, and it can be presented by the following Poisson’s type integral formula, written in the matrix form, u(ξ)=1a0-F(x)U(x,ξ)dx1dx2dx3--T(y)Q(y,ξ)dy2dy3,y(0,y2,y3), where |u(ξ)|< everywhere and vanish at infinity limξ1uk(ξ)0, lim|ξ2|uk(ξ)0, lim|ξ3|uk(ξ)0.

The matrices U(x,ξ) and Q(y,ξ)=U(y,ξ)/ny1 of influence of an inner unit point heat source F(ξ)=δ(x-ξ) and of a unit point temperature T(y)=δ(y-ξ) on the Γ onto thermoelastic displacements, also the matrix u(ξ) of searched displacements in (2.8), are determined as follows: U(x,ξ)=(U1(x,ξ)U2(x,ξ)U3(x,ξ))=γ8π(λ+2μ)(-(x1+ξ1R1+x1-ξ1R)(x2-ξ2)(1R1-1R)(x3-ξ3)(1R1-1R)),Q(y,ξ)=(Q1(y,ξ)Q2(y,ξ)Q3(y,ξ))=γ(λ+2μ)R(1-ξ12R2(y2-ξ2)ξ1R2(y3-ξ3)ξ1R2),u(ξ)=(u1(ξ)u2(ξ)u3(ξ)), where    R=R(x,ξ)=|x-ξ|; R1(x,ξ)=|x-ξ1*|; ξ(ξ1,ξ2,ξ3); ξ1*(-ξ1,ξ2,ξ3); x(x1,x2,x3)   in the matrix U(x,ξ), and     R=R(y,ξ)=|y-ξ|; R1(y,ξ)=|y-ξ1*|; ξ(ξ1,ξ2,ξ3); ξ1*(-ξ1,ξ2,ξ3); y(0,y2,y3) in the matrix Q(y,ξ).

Proof.

First, well-known Green’s function G for Poisson’s equation for a half-space is rewritten, and then, in the Subsection 2.2.1, we derive the volume dilatation Θ(k)(x,ξ). In Subsection 2.2.2 it is shown how to derive thermoelastic influence functions Uk(x,ξ). Checking the correctitude of the derived already thermoelastic influence functions Uk(x,ξ) is given in Subsection 2.2.3. Finally, in Subsection 2.2.4, on the base of the functions Uk(x,ξ), the Poisson’s type integral formula for stated BVP of thermoelasticity is derived and checked.

To obtain the matrix U(x,ξ) in (2.9) for the BVP in (2.4)–(2.6) we will use the integral formula in (2.1). The functions G(x,ξ) and Θ(k)(x,ξ) in this equation are the Green’s functions of Dirichlet problem in heat conduction and, respectively, the influence functions of a unit concentrated body force δikδ(x-ξ) onto volume dilatation in theory of elasticity for the half-space V. So, to get the Green’s function G(x,ξ) we have to solve the BVP, which consists of the equation in the heat conduction with the homogeneous boundary conditions similar to those in (2.6): x2G(x,ξ)=-δ(x-ξ),x,ξV,G=0,y1=0,-<y2,y3<. In the presented paper we recall this Green’s function from the handbook  (see problem and answer 15.P.1): G(x,ξ)=14π(R-1-R1-1),R=(x1-ξ1)2+(x2-ξ2)2+(x3-ξ3)2,R1=(x1+ξ1)2+(x2-ξ2)2+(x3-ξ3)2. This Green’s function is well known and is presented also in the books [18, 19].

2.2.1. Deriving the Volume Dilatation <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M121"><mml:mrow><mml:msup><mml:mrow><mml:mi mathvariant="normal">Θ</mml:mi></mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>x</mml:mi><mml:mo>,</mml:mo><mml:mi>ξ</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:math></inline-formula>

To get the influence functions Θ(k)(x,ξ), usually, we have to solve the following BVP, which consists of Lame’s equations in theory of elasticity and homogeneous boundary conditions as in (2.5) μ2Ui(k)(x,ξ)+(λ+μ)Θ,i(k)(x,ξ)=-δikδ(x-ξ),x,ξV,σ11(k)(y,ξ)=U2(k)(y,ξ)=U3(k)(y,ξ)=0,y1=0,-<y2,y3< and then, on the base of displacements Ui(k)(x,ξ), to compute volume dilatation Θ(k)(x,ξ)=Uj,j(k)(x,ξ). But as it will be shown below, in the case of the boundary condition in (2.12), we can derive the volume dilatation Θ(k)(x,ξ) using the equation 2Θ(k)(x,ξ)=-1λ+2μxkδ(x-ξ) only, and its integral representation via respective Green’s function GΘ(x,ξ): Θ(k)(x,ξ)=-1λ+2μξkGΘ(x,ξ)+Γ[Θ(k)(y,ξ)nΓ-Θ(k)(y,ξ)nΓ]GΘ(y,x)dΓ(y). To show this, first prove that the boundary conditions in (2.12) lead to the zero volume dilatation on the boundary plane Γ(y1=0;-<y2,y3<): Θ(k)(y,ξ)=0. To do this we write the Hooke’s law for the normal stresses σ11(k) in the form σ11(k)=2μU1,1(k)+λΘ(k)=(λ+2μ)Θ(k)-2μ(U2,2(k)+U3,3(k)),i,j=1,2,3. Next we know that the boundary conditions U2(k)=U3(k)=0 lead to the following zero tangential derivatives on the boundary plane Γ(y1=0;-<y2,y3<): U2,2(k)=U3,3(k)=0. Finally, substituting equalities (2.19) and boundary conditions σ11(k)=0 in the Hooke’s law (2.17) we come to the zero volume dilatation on the boundary Γ in (2.16).

So, as for the BVP for Θ(k), we have (2.14) and (2.16), then the respective BVP for Green’s function GΘ for considered half-space can be written in the form x2GΘ(x,ξ)=-δ(x-ξ),x,ξV,GΘ=0,y1=0,-<y2,y3<, which coincides with the BVP in (2.10). This fact leads to the conclusion that GΘG, and, as a result, the Green’s function GΘ is determined by expressions in (2.11).

Finally, if we introduce the expressions GΘG in (2.10), (2.11) and the equality Θ(k)(y,ξ)=0 in (2.16) into the representation (2.15), taking into account the equality [Θ(k)(y,ξ)nΓ-Θ(k)(y,ξ)nΓ]GΘ(y,x)=0, we obtain the searched volume dilatation of the elastic BVP in (2.12) for half-space, written in the form Θ(k)(x,ξ)=-14π(λ+2μ)ξk(R-1-R1-1).

2.2.2. Deriving the Thermoelastic Influence Functions <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M148"><mml:msub><mml:mrow><mml:mi>U</mml:mi></mml:mrow><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>x</mml:mi><mml:mo>,</mml:mo><mml:mi>ξ</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:math></inline-formula>

Now we have both functions: G(x,ξ) and Θ(k)(x,ξ) needed for deriving the thermoelastic influence functions Uk(x,ξ), using (2.1). So, substituting functions G(x,ξ) and Θ(k)(x,ξ) from (2.11) and (2.22) in (2.1), then, calculating the volume integral in the special way (see appendix), we will get the following expression for functions Uk(x,ξ): Uk(x,ξ)=γ0-G(x,z)Θ(k)(z,ξ)dz1dz2dz3=-γ0-14π(R-1-R1-1)14π(λ+2μ)ξk(R-1-R1-1)dz1dz2dz3,Uk(x,ξ)=γ(λ+2μ)-18πξk[R(x,ξ)-R1(x,ξ)], that being presented in the matrix form coincide with the matrix U(x,ξ) in (2.9). From the expressions Uk(x,ξ) in (2.23) we can see that |Uk(x,ξ)|< and that these displacements vanish at infinity: limξ1Uk(x,ξ)0, lim|ξ2|Uk(x,ξ)0, and lim|ξ3|Uk(x,ξ)0.

At the next step of the proof of the theorem we have to check the correctness of the functions Uk(x,ξ). To do this we have to use the fact that these functions satisfy both physical phenomena: elasticity in (1.5) and heat conduction in (1.4). So, they must satisfy the following BVP of thermoelasticity in respect to the coordinates of the point of application ξ(ξ1,ξ2,ξ3) of the inner unit heat source F(ξ)=δ(x-ξ): μξ2Uk(x,ξ)+(λ+μ)Θ,ξk(x,ξ)-γG,ξk(x,ξ)=0,σ11=U2=U3=0,ξ1=0,-<ξ2,ξ3<, which follows from the integral formula in (2.31) and boundary conditions in (2.12) (see also (1.5). Also, the functions Uk(x,ξ) have to satisfy the following fictive BVP in heat conduction with respect to the coordinates of the point of observation x(x1,x2,x3): x2Uk(x,ξ)=-γ  Θ(k)(x,ξ),x,ξV,Uk(y,ξ)=0,y(0,y2,y3)Γ, which follows from the integral formula in (2.23) and boundary conditions in (2.10) (see also (1.4). In addition, the functions Uk(x,ξ) in (2.23) have to vanish at infinity. In (2.25) the functions   Θ(k)(x,ξ) are determined by expression in (2.22), but in (2.24) the functions   Θ(x,ξ) are determined on the base of displacements Uk(x,ξ) in (2.23), using the rule  Θ(x,ξ)=Uk,k(x,ξ), where derivatives are taken with respect to coordinates of the point ξ(ξ1,ξ2,ξ3).

2.2.3. Checking the Correctitude of the Thermoelastic Influence Functions <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M175"><mml:msub><mml:mrow><mml:mi>U</mml:mi></mml:mrow><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>x</mml:mi><mml:mo>,</mml:mo><mml:mi>ξ</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:math></inline-formula>

To check (2.24) and (2.25) we need to compute the following values: ξ2R(x,ξ)=x2R(x,ξ)=2R-1(x,ξ),ξ2R1(x,ξ)=x2R1(x,ξ)=2R1-1(x,ξ),ξ2R-1(x,ξ)=ξ2R1-1(x,ξ)=x2R-1(x,ξ)=x2R1-1(x,ξ)=0. Check first (2.24). So, from (2.23), (2.24), and (2.26) it follows that μξ2Uk(x,ξ)=γμ4π(λ+2μ)ξk[R-1(x,ξ)-R1-1(x,ξ)],(λ+μ)  Θ,ξk(x,ξ)=γ(λ+μ)4π(λ+2μ)ξk[R-1(x,ξ)-R1-1(x,ξ)]. From (2.27) and (2.11) it follows that μξ2Uk(x,ξ)+(λ+μ)  Θ,k(x,ξ)=γ4πξk[R-1(x,ξ)-R1-1(x,ξ)]=γG,k(x,ξ), where the expression for Poisson’s function G(x,ξ) in (2.11) was used. So, due to (2.28), we come to the conclusion that the functions in (2.23) verify (2.24) at inner points. To check boundary conditions in (2.24) we need to calculate the derivatives in the expressions for displacements Uk(x,ξ) in (2.23) and calculate the values of the normal stresses for ξ1=0. So, using the Duhamel-Neumann law (2.3) for normal stresses σ11σ11=2μU1,1+(λΘ-γG), where the displacements U1, the volume dilatation Θ are determined on the base of expressions in (2.23), and the Green’s function G is determined by (2.11), we will obtain 2μU1,1(x,ξ)|ξ1=0=γμ(λ+2μ)-14π[1R(x,ξ)-1R1(x,ξ)-(x1-ξ1)2R(x,ξ)+(x1+ξ1)2R1(x,ξ)]|ξ1=0=0,λΘ(k)(x,ξ)|ξ1=0=λγ4π(λ+2μ)(R-1-R1-1)|ξ1=0=0,G(x,ξ)|ξ1=0=14π(R-1-R1-1)|ξ1=0=0. Substituting expressions in (2.30) in (2.29) we obtain σ11|ξ1=0=0, so the first boundary condition in (2.24) is satisfied. Using the expressions in (2.23) for U2 and U3, and we obtain Uk(x,ξ)|ξ1=0=γ(λ+2μ)-18πξk[R(x,ξ)-R1(x,ξ)]|ξ1=0=0,k=2,3. So, as we can see from (2.31) and from σ11|ξ1=0=0, the BVP in (2.24) is satisfied. Now we can confirm that the thermoelastic influence functions Uk(x,ξ) satisfy BVP of thermoelasticity in (2.24). Next we have to prove that the influence functions in (2.23) satisfy with respect to the coordinates of the point of observation x(x1,x2,x3) the fictive BVP in heat conduction, described by (2.25). Indeed, using (2.22), (2.26) and (2.23), we will obtain x2Uk(x,ξ)=γ4π(λ+2μ)ξk(R-1(x,ξ)-R1-1(x,ξ))=-γΘ(k)(x,ξ), which means that at inner points (2.25) is satisfied. Also, taking into account expressions in (2.23) for displacements Uk(x,ξ) we come to the conclusion that the boundary conditions in (2.25) on the boundary plane are also satisfied, which means Uk(y,ξ)=0; y(0,y2,y3)Γ.

2.2.4. Deriving and Checking the Poisson’s Type Integral Formula

The next step is to calculate the other influence functions in (1.10) Uk(y,ξ)/ny1 and to get Poisson’s type integral formula. So, to get the Poisson’s type integral formula we have to use formula (2.2) for the half-space. The influence functions Uk(x,ξ) in this formula are determined by (2.23). The functions Uk(y,ξ)/ny1, as functions of influence of unit point temperature, given on the boundary plane, according to the formula in (1.10) have to be determined on the basis of the functions Uk(x,ξ) in (2.23) as following   Uk(y,ξ)ny1=-limxyUk(x,ξ)x1=-limxyγ(λ+2μ)-18π2x1ξk[R(x,ξ)-R1(x,ξ)]=limxyγ(λ+2μ)-14π2ξ1ξkR(x,ξ),Uk(y,ξ)ny1=γ(λ+2μ)-14π2ξ1ξkR(y,ξ). As we see, the expressions for Uk(y,ξ)/ny1 in (2.33) vanish at infinity.

Introducing the influence functions in (2.23) and (2.33) into the formula in (2.2) we obtain the following Poisson’s type integral formula: Uk(ξ)=a-10-F(x)Uk(x,ξ)dx1dx2dx3--T(y)Uk(y,ξ)ny1dy2dy3=γ8π(λ+2μ)[a-10-F(x)ξk[R(x,ξ)-R1(x,ξ)]dx1dx2dx3,-2-T(y)2ξ1ξkR(y,ξ)dy2dy3],|uk(ξ)|<. Finally, to notice that integrals with unbounded intervals will exist, which means that the displacements |uk(ξ)|<, when the following conditions are satisfied: 0-|F(x)|dx1dx2dx3<,-|T(y)|dy2dy3<, because, as was mentioned above, the kernels in (2.34) vanish at infinity. The conditions in (2.35) will be satisfied in case the functions F(x) and T(y) are given on the bounded domains.

If in (2.23) and (2.33) we calculate the derivatives, and, if we present the obtained results in matrix form, then we can be sure that the influence matrices U(x,ξ), U(y,ξ)/ny1=Q(y,ξ), and Poisson’s type integral formula in (2.34) coincide with the results in (2.8) and (2.9).

The obtained displacements, described by Poisson’s type integral formula (2.34), must satisfy the stated in the theorem BVP of thermoelasticity, described by (2.4) and (2.5). To check this we substitute the integral formula from (2.34) into (2.4): μ2uk(ξ)+(λ+μ)θ,k(ξ)=a-10-F(x)[μξ2Uk(x,ξ)+(λ+μ)Θ,k(x,ξ)]dx1dx2dx3--T(y)ny1[μξ2Uk(y,ξ)+(λ+μ)Θ,k(y,ξ)]dy2dy3. Taking into account the already proved equality in (2.28) and the Poisson’s type integral formula for BVP in heat conduction, described by (2.6), T(ξ)=0-F(x)G(x,ξ)dx1dx2dx3--T(y)ny1G(y,ξ)  dy2dy3, the following is obtained: μ2uk(ξ)+(λ+μ)θ,k(ξ)=a-10-F(x)G,k(x,ξ)dx1dx2dx3--T(y)ny1G,k(y,ξ)  dy2dy3=T,k(ξ). So, the displacements, described by Poisson’s type integral formula (2.34), satisfy (2.4). Finally, we have to show that the displacements in (2.34) satisfy the mechanical boundary conditions in (2.5). Of course, these conditions are satisfied, because, as we have shown above, the kernels of Poisson’s type integral formula (2.34) satisfied the same mechanical boundary conditions in (2.24) that are analogical to (2.5). So, now we are sure that all items of the theorem are proved. It means that the theorem is proved.

3. An Example of Application of the Obtained Poisson’s Type Integral Formula

Let us solve (2.4)–(2.6) and to determine thermoelastic displacements uk(ξ) in half-space if on the segment [y1=0, -by2b; y3=0] of its boundary plane Γ(y1=0;<y2,y3<), the constant temperature

T(y)={T0=const,y[y1=0,-by2b;y3=0];0,y[y1=0,-by2b;y3=0] is given. To solve this BVP we use the integral formula in (2.8) and (2.9) for F(x)=0 and T(y)=T0=const: u(ξ)=-T0-bbQ(y;ξ)dy2,y(0,y2,0),-by2b, where Q(y;ξ) is determined by the matrix (2.9). The formulas in (3.1) and (2.9) can be presented in terms of components of thermoelastic displacements uk; k=1,2,3 as follows (see (2.9) and (2.34): uk(ξ)=-γT04π(λ+2μ)-bb2ξ1ξkR(y2;ξ)dy2,R(y2;ξ)=ξ12+(y2-ξ2)2+ξ32,ξ(ξ1,ξ2,ξ3). So, taking the respective integrals over the boundary segment [y1=0,-by2b;y3=0], we obtain the following thermoelastic displacements at an arbitrary inner point ξ(ξ1,ξ2,ξ3): uk(ξ)=-γT04π(λ+2μ)-bb2ξkξ1R(y2;ξ)dy2=-γT04π(λ+2μ)ξk-bbξ1R(y2;ξ)dy2=-γT04π(λ+2μ)ξk[ξ1ln|y2-ξ2+R(y2;ξ)|-bb],uk(ξ)=-γT04π(λ+2μ)ξk[ξ1L(ξ)], where L(ξ)=ln|b2-ξ2+R(b;ξ)-(b+ξ2)+R(-b;ξ)|,R(b;ξ)=ξ12+(b-ξ2)2+ξ32,R(-b;ξ)=ξ12+(b+ξ2)2+ξ32. Check whether thermoelastic displacements uk(ξ) in (3.4) and (3.5) satisfy the stated in this section BVP, described by (2.4)–(2.6) and (3.1). For this, first calculate μ2uk(ξ)=-μγT04π(λ+2μ)2ξk[ξ1L(ξ)]=-γμT02π(λ+2μ)2ξkξ1L(ξ),(λ+μ)ξkΘ(k)(ξ)=(λ+μ)2ξkξjuj(ξ)=-(λ+μ)γT04π(λ+2μ)3ξkξj2[ξ1L(ξ)]=-(λ+μ)γT04π(λ+2μ)ξk2[ξ1L(ξ)]=-(λ+μ)γT02π(λ+2μ)2ξkξ1L(ξ). In (3.6) and (3.7) the Euler’s formula was used: 2(f·φ)=φ2+f2φ+2f,jφ,j. Also it was taken into account that the function L(ξ), described by (3.5), is a harmonic function, so that 2L(ξ)=0. Next calculate the temperature T(ξ) using the well-known Poisson’s integral formula in (2.37) rewritten in the case when on the segment [y1=0,-by2b; y3=0] of its boundary plane Γ(y1=0;<y2,y3<) the constant temperature T=T0 is given: T(ξ)=-T0-bbny1G(y2,ξ)dy2=-T0ξ1-bbdy22πR(y2;ξ),T(ξ)=-T02πξ1L(ξ). It is easy to see that temperature in (3.9) satisfies the BVP in heat conduction, described by Poisson equation (2.6) [for F(x)=0] and boundary condition in (3.1). So, substituting the expressions in (3.6), (3.7), and (3.9) into (2.4) we can see that this equation is satisfied. It remains to show that the thermoelastic displacements uk(ξ) in (3.4) and (3.5) satisfy the boundary conditions in (2.5), and also that they vanish at infinity. From (3.4) it is easy to see that for k=2,3 the thermoelastic displacements u2(ξ) and u3(ξ) for ξ1=0 satisfy the second and third boundary conditions in (2.5). To check the first boundary condition in (2.5) calculate the normal stresses σ11(ξ) using the Duhamel-Neumann law, the displacements in (3.4) and the temperature in (3.9): σ11=2μU1,1+λUj,j-γT,          i,j=1,2,3,        σ11=(λ+2μ)U1,1-λ(U2,2+U3,3)-γT,σ11=-(λ+2μ)γT04π(λ+2μ)2ξ12[ξ1L(ξ)]-λ(U2,2+U3,3)-γT,σ11=-γT04π(2+ξ1ξ1)ξ1L(ξ)-λ(U2,2+U3,3)+T02πξ1L(ξ). As on the boundary plane (for ξ1=0) u2=u3=0 (which means also u2,2=u3,3=0), then from the last equation it follows σ11=0 [the first boundary condition in (2.5) is satisfied].

Finally, calculating the derivatives of expressions in (3.4) and (3.5), we obtain the thermoelastic displacements of the stated in this section BVP in the explicit form u1(ξ)=γT04π(λ+2μ){ln|b2-ξ2+R(b;ξ)-(b+ξ2)+R(-b;ξ)|+ξ12[R-1(b;ξ)(b2-ξ2+R(b;ξ))-1-R-1(-b;ξ)(-(b2+ξ2)+R(-b;ξ))-1]};u2(ξ)=γT04π(λ+2μ)ξ1[R-1(b;ξ)-R-1(-b;ξ)];u3(ξ)=γT04π(λ+2μ)ξ1ξ3[R-1(b;ξ)(b2-ξ2+R(b;ξ))-1-R-1(-b;ξ)(-(b2+ξ2)+R(-b;ξ))-1], which vanish at infinity.

4. Conclusions

The Poisson’s type integral formula, obtained in this paper, is new, useful, and completely ready to be efficiently applied for computing of the thermoelastic displacements uk(ξ) in half-space (see example of its application given in Section 3). The main advantage of the obtained integral formula is that the searched half-space thermoelastic displacements are expressed directly via given inside heat source, boundary temperature, and known kernels. So it is not necessary to determine intermediately the inner temperature field or to solve, as in traditional methods, additional BVP.

The most difficult problems in the proposed here method are the problems of deriving the Green’s functions G(x,ξ) in heat conduction and the functions of influence for volume dilatation Θ(k)(x,ξ) in the elasticity theory. Note that for canonical Cartesian domains, these problems were solved successfully in the handbook , where about 190 Green’s functions for Laplace’s equation and 250 influence functions for volume dilatation Θ(k)(x,ξ) are presented. At last, computing the integral over volume in (1.4) of the product of functions G(x,ξ) and Θ(k)(x,ξ) is also solved successfully [see the recommendations of the appendix]. So, for canonical Cartesian domains, the proposed in thermoelasticity method will work successfully. This means that the presented paper opened great possibilities for researchers to derive many new Poisson’s type integral formulas, not only for half-space, but also for many other canonical Cartesian domains. These new Poisson’s type integral formulas are very useful to solve effectively not only deterministic BVPs of thermoelasticity, but also the stochastic ones .

The approach presented in this paper in thermoelasticity for Cartesian canonical domains can be extended onto spherical [21, 22], polar , cylindrical, and other canonical domains of any orthogonal systems of coordinates. This extension will be done when the lists of the respective functions G(x,ξ) and Θ(k)(x,ξ) are completed.

The approach presented in this paper is valid also for other physical phenomena as electroelasticity, magnetoelasticity, and poroelasticity, described by the same BVP as in thermoelasticity.

Appendix

The improper integral in (2.23) was taken using:

the following equalities on the boundary plane of the half-space: [R(x,z)-R1(x,z)]|z=y(0,y2,y3)=0,[R-1(z,ξ)-R1-1(z,ξ)]|z=y(0,y2,y3)=0,

the relations z2R(x,z)=2R-1(x,z),z2R1(x,z)=2R1-1(x,z),

the following property of Dirac’s function: Vf(x)δ(x,ξ)dV(x)=f(ξ),

the Green’s formula inside the half-space V(φ2f-f2φ)dV=Γ[φ(fn)-f(φn)]dΓ, where f=[R(x,z)-R1(x,z)],φ=(4π)-1[R-1(z,ξ)-R1-1(z,ξ)]. So, using the recommended equations (A.1)–(A.5), we obtain Uk(x,ξ)=γ0-G(x,z)Θ(k)(z,ξ)dz1dz2dz3=-γ0-14π(R-1-R1-1)14π(λ+2μ)ξk(R-1-R1-1)dz1dz2dz3=-γ8π(λ+2μ)ξk0-z2[R(x,z)-R1(x,z)]14π[R-1(z,ξ)-R1-1(z,ξ)]dz1dz2dz3=-γ8π(λ+2μ)ξk0-[R(x,z)-R1(x,z)]14πz2[R-1(z,ξ)-R1-1(z,ξ)]dz1dz2dz3=γ8π(λ+2μ)ξk0-[R(x,z)-R1(x,z)]δ(z-ξ)dz1dz2dz3=γ8π(λ+2μ)ξk[R(x,ξ)-R1(x,ξ)], which coincides with the final expression in (2.23).

Acknowledgments

One of the authors, Dr V. Seremet, is grateful to the University Paris-Est, Marne-la-Vallee, France, for the support of his two research visits to the university. Also he expresses many thanks to the reviewers of this paper, whose comments have contributed substantially to its improvements.

BoleyB. A.WeinerJ. H.Theory of Thermal Stresses1960New York, NY, USAJohn Wiley & Sonsxvi+586MR0112414ZBL0095.18407KovalenkoA. D.The Basis of Thermoelasticity1970Kiev, UkraineNaukova DumkaMelanE.ParkusH.Thermoelastic Stresses Caused by the Stationary Heat Fields1958Moscow, RussiaFizmatgizMayzelV. M.The Temperature Problem of the Theory of Elasticity1951Kiev, UkraineAN SSSRNowackiW.Thermoelasticity1962Warszawa, PolandPergamon Press, Oxford, UK; PWNNowackiW.TheTheory of Elasticity1975Moscow, RussiaMirSheremetV. D.The Integral Equations and Green's Matrices of the Method of Elements of Influence in the Mechanics of Solids, Thesis of the doctor of physical1995Chisinau, MoldovaMathematical Sciences, Technical University of MoldovaŞeremetV. D.The modification of Maizeli's formula in the stationary thermoelasticityBulletin of Academy of Science of Republic of Moldova, Mathematics199731922MR1602050ŞeremetV. D.Some new influence functions and integral solutions in theory of thermal stressesProceedings of the IVth International Congress on Thermal StressesJune 2001Osaka, Japan423427ŞeremetV. D.New results in 3-D thermoelasticityProceedings of the 14th U.S. National Congress of Theoretical and Applied Mechanics, Virginia TechJune 2002Blacksburg, Va, USAŞeremetV. D.Handbook of Green's Functions and Matrices2003Southampton, UKWIT Pressii+295MR1945701SheremetV.New formulae for dynamical thermal stressesJournal of Thermal Stresses200225212315310.1080/014957302753384388SheremetV.Generalization of Green's formulae in thermoelasticityCollection: Multiscale Green's Functions for Nanostructures, National Science Digital Library of USA, NSF, p. 4, 2003ŞeremetV.Influence Elements Method2003Chisinau, MoldovaEdit. Center, Agrar. Univ. of MoldovaŞeremetV. D.VladI.ŞeremetA.New integral formulae in thermoelasticityProceedings of the 16th ASCE Engineering Mechanics Conference (EM '03)July 2003Seattle, Wash, USAWashington University110SheremetV. D.PrecupanD.VladI.SheremetA. V.The constructing of Green's matrices in cylindrical coordinatesProceedings of the 17th Engineering Mechanics Conference of the American Society of Civil EngineersJune 2004Newark, Del, USAUniversity of Delaware Newark19ŞeremetV.BonnetG.SpeianuT.New results in construction of the Green's matrices in spherical coordinatesProceedings of the Inaugural International Conference of the Engineering Mechanics Institute (EM '08)May 2008University of Minnesota, USA17BudacV. M.SamarskiiA. A.TihonovA. N.Set of Problems on Mathematical Physics1980Moscow, RussiaGostehnizdatKartashovA. D.Analytic Methods in the Theory of Heat Conductibility of Rigid Bodies1980Moscow, RussiaHigh SchoolSabelfeldK.Expansion of random boundary excitations for elliptic PDEsMonte Carlo Methods and Applications2007135-640545310.1515/mcma.2007.022MR2379265ZBL1158.60027ŞeremetV. D.VladI.ŞeremetA.New influence functions for thermoelastic spherical shellsProceedings of the 5th International Congress on Thermal Stresses (ICTS '03), Virginia Tech.June 2003Blacksburg, Va, USA14ŞeremetV.BonnetG.SpeianuT.Influence functions and integral formulae for spherical thermoelastic bodiesProceedings of the 22nd International Congress of Theoretical and Applied Mechanics (ICTAM '08)August 2008Adelaide, AustraliaAdelaide University12ŞeremetV.BonnetG.Encyclopedia of Domain Green's Functions (Thermo-Magnetoelectrostatics of Solids in Rectangular and Polar Coordinates)2008Chisinau, MoldovaEdit. Center, Agrar. Univ. of Moldova