Next Article in Journal
A Compact Octa-Band Frequency Reconfigurable Antenna for Wireless Applications
Next Article in Special Issue
A New Class of Plane Curves with Arc Length Parametrization and Its Application to Linear Analysis of Curved Beams
Previous Article in Journal
Selecting Bloggers for Hotels via an Innovative Mixed MCDM Model
Previous Article in Special Issue
A Study of Yielding and Plasticity of Rapid Prototyped ABS
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Two-Dimensional Thermoelasticity Solution for Bimodular Material Beams under the Combination Action of Thermal and Mechanical Loads

1
School of Civil Engineering, Chongqing University, Chongqing 400045, China
2
Key Laboratory of New Technology for Construction of Cities in Mountain Area, Chongqing University, Ministry of Education, Chongqing 400045, China
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(13), 1556; https://doi.org/10.3390/math9131556
Submission received: 25 May 2021 / Revised: 25 June 2021 / Accepted: 29 June 2021 / Published: 2 July 2021
(This article belongs to the Special Issue Applied Mathematics and Solid Mechanics)

Abstract

:
A typical characteristic of bimodular material beams is that when bending, the neutral layer of the beam does not coincide with its geometric middle surface since the mechanical properties of materials in tension and compression are different. In the classical theory of elasticity, however, this characteristic has not been considered. In this study, a bimodular simply-supported beam under the combination action of thermal and mechanical loads is theoretically analyzed. First, a simplified mechanical model concerning the neutral layer is established. Based on this mechanical model, Duhamel’s theorem is used to transform the thermoelastical problem into a pure elasticity problem with imaginary body force and surface force. In solving the governing equation expressed in terms of displacement, a special solution of the displacement equation is found first, and then by utilizing the stress function method based on subarea in tension and compression, a supplement solution for the displacement governing equation without the thermal effect is derived. Lastly, the special solution and supplement solution are superimposed to satisfy boundary conditions, thus obtaining a two-dimensional thermoelasticity solution. In addition, the bimodular effect and temperature effect on the thermoelasticity solution are illustrated by computational examples.

1. Introduction

In the classical theory of elasticity [1], it is generally assumed that materials exhibit the same elastic properties in tension and compression, but this is only a simplified result and does not account for the nonlinear characteristics of materials. Many studies have indicated that most materials [2,3], including concrete, ceramics, graphite, rubber and some biomedical materials, exhibit different tensile and compressive strains when they are subjected to tensile and compressive stresses with the same magnitude. Thus, these materials exhibit different elastic moduli in tension and compression, and are called bimodular materials [4]. Generally, there are two basic material models widely used in theoretical analysis within the engineering profession. One is the criterion of positive–negative signs in the longitudinal strain of fibers proposed by Bert [5]. This model is mainly applicable to orthotropic materials and is therefore widely used in the analysis of laminated composites [6,7,8,9]. Another model is the criterion of positive–negative signs of principal stress proposed by Ambartsumyan [10], which is mainly applicable to isotropic materials. In structural analysis, this model is of particular significance, since it is this factor that determines whether the point is in tension or in compression. This paper will focus on the study of the latter model based on principal direction.
Ambartsumyan [10] linearized the bimodular materials model by two straight lines whose tangents at the origin are discontinuous, as shown in Figure 1, in which σ is the principal stress and ε is the principal strain. The basic assumptions of this model are as follows: (1) The studied body is continuous, elastic, homogeneous, and isotropic; (2) Small deformation is satisfied; (3)Young’s modulus of elasticity and Poisson’s ratio of materials are E + and μ + when the materials are in tension along a certain principal direction; and they are E and μ when the materials are in compression; (4) When the three principal stresses are uniformly positive or uniformly negative, the three basic equations are essentially the same as those of the classical theory of elasticity; when the signs of the three principal stresses are different, the differential equations of equilibrium and the geometrical equations are the same as those of the classical theory of elasticity, with the exception of the physical equations; (5) μ + / E + = μ / E is introduced and this will ensure symmetry of the flexibility matrix.
Due to the fact that this bimodular theory defines its constitutive model based on principal directions, the principal stress is generally obtained as a final result but not as a known condition before solving; this inevitably incurs difficulties of description of the stress state of a point. In addition, this model lacks the ability to describe experimental results of elastic coefficients in the complex states of stress. Analytical solutions are available in a few simple cases, although they only concern beams and plates [11,12,13,14]. In some complex problems, it is necessary to resort to the finite element method (FEM) based on an iterative technique [15,16,17,18,19]. According to our literature collection, there are only a few works on the application of a bimodular materials model in the analysis of thermal stress, and even if there are several studies, most of them are based on the Bert model, not on the Ambartsumyan model.
In the theory of thermoelasticity [20], the influence of temperature field in the governing equations is presented through the constitutive law, in which the modulus of elasticity of materials is an important constant. The theory of linear thermoelasticity is established on linear supplement of thermal strains to mechanical strains. Generally, problems of thermoelasticity have been solved by finding solutions of the Lamé displacement equations when a body is acted upon by arbitrary mass forces. Thus, many basic thermoelasticity problems have been considered within the classical theory of elasticity. This is the classical body force analogy, which may date back to Duhamel, who made great contribution to this field in history.
With the development of the theory of thermoelasticity, some generalized thermoelastic models have been proposed for transient responses in many applications (for example, ultra-fast lasers heating and low temperatures) where the classical theory of thermoelasticity fails. Some representative theories in this regard can be found in [21,22,23,24]. It should be emphasized that Green and Lindsay’s theory has been addressed for many types of media, in which Marin et al. [25] used this theory in the context of the dipolar thermoelastic bodies. On the other hand, besides the development of the theory itself, it is also very important to apply the theory to the analysis of the engineering components, or more specially, it is equally important to analyze the thermoelastic behavior of engineering structures, for example, the thermoelasticity behaviors of nanobeams [26], microbeams [27], composite beams [28], and laminated beams [29]. Of course, the bimodular material beams should also be included in the thermoelastical analysis of structures.
At present, the bimodular problem is mostly not considered in the theory of thermoelasticity. Once the bimodular problem is introduced into the theory of thermoelasticity, it is predictable that there will be changes in the governing equation used for the solution of the bimodular thermal problem and in the solving method, and what are the corresponding changes? This study will focus on the establishment of the displacement governing equation and the corresponding solving process, with the emphasis on the differences introduced by bimodular characteristics of materials.
This study is devoted to obtaining a two-dimensional thermoelasticity solution for the bimodular beam under the combination action of thermal and mechanical loads. For this purpose, this paper is organized as follows. First, the mechanical model based on the neutral layer is established in Section 2, thus realizing the regional segmentation of tension and compression. The governing equation expressed in terms of displacement components is given in Section 3, and the composition of the solution is presented in Section 4. Aiming at the governing equation, the special solution and corresponding supplement solution are derived in Section 5, in which the supplement solution is obtained by applying the stress function method as well as the de Saint-Venant’s Principle. In Section 6, some typical computational examples are given, in which the bimodular effect and temperature effect on the thermoelasticity solution are discussed. Section 7 is the conclusions.

2. Mechanical Model Based on Neutral Layer

A rectangular section beam with the height h and thickness t is subjected to the bending moment M at its two ends, and we isolated a part from the whole beam as our studied object, as shown in Figure 2. Obviously, the beam will deflect downward to resist the external bending moment, thus resulting in tension in the lower part of the beam and compression in the upper part. For convenience, we established the x-axis of the coordinate system (xoy) on the unknown neutral axis, which may be determined thereafter. The height of the tensile part of the beam is h 1 and the corresponding modulus in this part is E + ; while at the same time, the compressive height and the modulus in the compressive part is h 2 and E , as shown in Figure 2, in which y 0 is the offset distance between the neutral axis and the geometrical middle layer.
According to the plane section assumption, the curvature ρ of the neutral axis may be expressed as
1 ρ = d θ d x ,
in which, d x and d θ denote the length and rotation angle of a segment on the neutral axis, respectively, as shown in Figure 2. The relative elongation of the segment AB with distance y from the neutral axis, ε x , is equal to
ε x = ( y + ρ ) d θ ρ d θ d x = y d θ d x = y ρ .
Thus the tensile normal stress and compressive normal one may be expressed as follows
σ x + = E + y ρ ,   0 y h 1 ;   σ x = E y ρ ,   h 2 y 0 .
The condition of equilibrium gives
0 h 1 σ x + t d y + h 2 0 σ x t d y = 0 ,
and
0 h 1 σ x + y t d y + h 2 0 σ x y t d y = M .
Combining Equations (3) and (4) yields
E + h 1 2 = E h 2 2 .
After considering h 1 + h 2 = h , we have
h 1 = E E + + E h , h 2 = E + E + + E h
which determines the location of the unknown neutral axis. Similarly, combining Equations (3) and (5) yields
1 ρ ( E + t h 1 3 3 + E t h 2 3 3 ) = M .
If D is defined as the bending stiffness of the bimodular problem, that is,
D = t 3 ( E + h 1 3 + E h 2 3 ) .
Equation (8) will return to a familiar form
1 ρ = M D .
When E + = E = E , we have h 1 = h 2 = h / 2 , the above expressions may be reduced to the classical forms, in which D = E t h 3 / 12 . Thus, we establish the so-called mechanical model on subarea in tension and compression under mechanical loads. Next, based on this model, we will analyze the thermal stress of a bimodular beam under the combination action of mechanical and thermal loads.
It should be noted here that if the bimodular materials beam is subjected to uniformly distributed loads, not only the bending moment but also the shearing forces will develop in the beam, that is to say, the beam is under lateral bending but under pure bending. In this case, the mechanical model based on the neutral layer is still applicable since the shearing stress has no influence on the determination of the neutral axis, which depends only on the bending stress. This conclusion was demonstrated by Yao and Ye [11].

3. Displacement Governing Equation

The problem we consider here is a typical plane stress problem concerning thermal effect. If we let ε x , ε y and γ x y be the strain components of a plane stress problem and σ x , σ y and τ x y be the corresponding stress components, α represents the linear thermal expansion coefficient and T represents the temperature variation, the physical equation of the two-dimensional theory of thermoelasticity may be given as follows
{ ε x = 1 E + / ( σ x μ σ y ) + α T ε y = 1 E + / ( σ y μ σ x ) + α T γ x y = 2 ( 1 + μ ) E + / τ x y ,
in which the tensile-compressive Young’s modulus is denoted by E + / . In our study, however, there is no way to specify that, from the very beginning, the stress state of any point is tensile or compressive, thus we have to express the modulus in the form of E + / . Once the final result is obtained, with the stress state becoming clear, we can differentiate the tension from the compression. This fact is quite different from the case that is only under mechanical load shown in Section 2, in which the so-called subarea in tension and compression, that is, the neutral axis, is obvious from the beginning and may also be determined later. However, under the combination action of mechanical and thermal loads, not only the bending from the mechanical load, but also the axial tension or compression from the temperature rise, may have influences on the final stress results, thus the so-called neutral layer becomes obscure at present. Therefore, the form E + / has to be adopted to express the physical equation for the time being.
For the physical equation, we may have another form as follows,
{ σ x = E + / 1 μ 2 ( ε x + μ ε y ) α T E + / 1 μ σ y = E + / 1 μ 2 ( ε y + μ ε x ) α T E + / 1 μ τ x y = E + / 2 ( 1 + μ ) γ x y .
If we let u and v be the displacement along x and y direction, respectively, the geometrical equation of the plane problem gives [1]
ε x = u x , ε y = v y , γ x y = v x + u y .
Substituting Equation (13) into Equation (12), we may obtain the stress components expressed in terms of the displacement u and v and the temperature variation T as
{ σ x = E + / 1 μ 2 ( u x + μ v y ) α T E + / 1 μ σ y = E + / 1 μ 2 ( v y + μ u x ) α T E + / 1 μ τ x y = E + / 2 ( 1 + μ ) ( v x + u y ) .
Substituting Equation (14) into the equation of equilibrium of the plane problem
{ σ x x + τ x y y + X = 0 σ y y + τ x y x + Y = 0 ,
in which X and Y are the body forces along x and y directions, respectively, we will have the Lamé equation, which may be used in the displacement method (noting X = Y = 0 ).
{ 2 u x 2 + 1 μ 2 2 u y 2 + 1 + μ 2 2 v x y ( 1 + μ ) α T x = 0 2 v y 2 + 1 μ 2 2 v x 2 + 1 + μ 2 2 u x y ( 1 + μ ) α T y = 0 .
In addition, substituting Equation (14) into the following stress boundary condition of the plane problem
{ l ( σ x ) s + m ( τ x y ) s = X ¯ m ( σ y ) s + l ( τ x y ) s = Y ¯ ,
in which, l and m are direction cosines, and X ¯ and Y ¯ are the surfaces forces, we will have the stress boundary condition used in the displacement method (noting X ¯ = Y ¯ = 0 ).
{ l ( u x + μ v y ) s + m 1 μ 2 ( u y + v x ) s = l ( 1 + μ ) α T m ( v y + μ u x ) s + l 1 μ 2 ( v x + u y ) s = m ( 1 + μ ) α T .
Therefore, the Lamé equation and corresponding boundary conditions of a two-dimensional thermoelasticity problem are established, as shown in Equations (16) and (18), also including displacement boundary condition not listed here.
Comparing Equations (16) and (18) with the counterparts in the plane stress problem without thermal effect, it is found that the original X and Y are now replaced by α E + / ( T / x ) / ( 1 μ ) and α E + / ( T / y ) / ( 1 μ ) , respectively, and the original X ¯ and Y ¯ are now replaced by l α T E + / / ( 1 μ ) and m α T E + / / ( 1 μ ) , respectively. This fact suggests to us that under certain displacement boundary conditions, the displacement due to temperature variation T is equal to the displacement of the elastic body without temperature variation, which is subjected to the imaginary body forces, α E + / ( T / x ) / ( 1 μ ) and α E + / ( T / y ) / ( 1 μ ) , and the imaginary surface forces, l α T E + / / ( 1 μ ) and m α T E + / / ( 1 μ ) . Therefore, the thermoelasticity plane stress problem is transformed into a pure elasticity problem under the action of the known body forces and known surface forces; this is Duhamel’s theorem that we are familiar with.
It should be noted here that there exists no term of modulus of elasticity in Equations (16) and (18) since this term is eliminated during the derivation of the two equations. For the classical thermoelasticity problem with singular modulus, this conclusion is true, which may be easily found in many textbooks. However, once the bimodular effect is introduced, will the form of the equations change? From Equations (14) and (15), we may see that whether the tensile modulus or the compressive modulus may be eliminated will depend on the fact that σ x and τ x y have the same positive or negative sign, or σ y and τ x y have the same sign. First, σ x and σ y are two normal stresses along the directions of x and y, respectively, and they may be positive or negative; in this case, the positive or negative sign of τ x y is a relatively key factor, that is to say, whether the sign of τ x y may change as the sign of the corresponding σ x and σ y during the derivation. In the bimodular theory proposed by Ambartsumyan, there is a basic assumption that μ + / E + = μ / E , which ensures the symmetry of the flexibility matrix. In this study, we ignore the difference of μ + and μ , regarding them as a constant μ , thus when only τ x y is considered, according to μ + / E + = μ / E , we may have E + = E only in this case. That is to say, the positive of negative sign of τ x y may change as the positive of negative sign of σ x and σ y . Therefore, even if the bimodular effect is introduced, the basic forms of Equations (16) and (18) remain unchanged.

4. Composition of Solution

Generally speaking, under the stress boundary condition (18) and displacement boundary condition, it is hard to obtain the analytical solution of the governing Equation (16). In an actual solving process, it is a common practice to follow the next two steps. First, any special solution of Equation (16) is found that does not necessarily satisfy the boundary conditions. Second, another supplement solution is found, not considering temperature changes, and then by superimposing the supplement solution and the previous special solution, the total stress is obtained to satisfy the boundary conditions.
To find the special solution of Equation (16), a potential function of displacement, ψ ( x , y ) , is introduced and the special solution of displacement may be taken as
u = ψ x , v = ψ y .
Regarding the u and v as u and v, respectively, and substituting them into Equation (16), we have
{ x 2 ψ = ( 1 + μ ) α T x y 2 ψ = ( 1 + μ ) α T y .
Note that μ and α are both constants in this problem, this shows if the function ψ ( x , y ) can satisfy the following differential equation
2 ψ = ( 1 + μ ) α T ,
then the ψ ( x , y ) can satisfy Equation (20), thus also satisfying Equation (16). Finally, the u and v in Equation (19) may be selected as a special solution of displacement. Substituting Equations (19) and (21) into Equation (14), the stress component corresponding to the special solution of displacement may be obtained as
{ σ x = E + / 1 + μ 2 ψ y 2 σ y = E + / 1 + μ 2 ψ x 2 τ x y = E + / 1 + μ 2 ψ x y .
On the other hand, the supplement solution of displacement, u and v , need to satisfy the homogeneous form of Equation (16), this gives
{ 2 u x 2 + 1 μ 2 2 u y 2 + 1 + μ 2 2 v x y = 0 2 v y 2 + 1 μ 2 2 v x 2 + 1 + μ 2 2 u x y = 0 .
By using Equation (14) and also letting T = 0 , we have the stress components corresponding to the supplement solution of displacement,
{ σ x = E + / 1 μ 2 ( u x + μ v y ) σ y = E + / 1 μ 2 ( v y + μ u x ) τ x y = E + / 2 ( 1 + μ ) ( v x + u y ) .
Lastly, the total displacement is
u = u + u , v = v + v ,
which needs to satisfy displacement boundary conditions. In addition, the total stress is
σ x = σ x + σ x , σ y = σ y + σ y , τ x y = τ x y + τ x y   ,
which needs to satisfy stress boundary conditions.

5. Obtainment of Solution

5.1. Special Solution

If the temperature variation is always a function of y, that is, T = T ( y ) , it is easy to obtain the special solution according to the description above. Substituting Equation (21) into Equation (22), the following stress components may be obtained
{ σ x = E + / α T ( y ) = E α T ( y ) σ y = 0 τ x y = 0 ,
which corresponds to the special solution of displacement. Note that the σ x in the special solution always corresponds to the compressive stress since the linear thermal expansion coefficient α is positive and T ( y ) also represents the temperature rise. Therefore, the E + / in the term σ x should be modified as E , according to the bimodular characteristics of tension and compression described above.

5.2. Supplement Solution and Superposition with Special Solution

In an elasticity problem only with stress boundary conditions, it is hard to find the supplement solution of displacement. To overcome the difficulties, it is a common practice to adopt the stress function method, that is, introduce a stress function, φ ( x , y ) , and then express the stress components corresponding to the supplement solution of displacement as follows
σ x = 2 φ x 2 , σ y = 2 φ y 2 , τ x y = 2 φ x y   ,
in which the selection for φ ( x , y ) may be carried out according to the specific stress boundary conditions, which has been demonstrated in many real examples.
As shown in Figure 3, a simply-supported beam in a thermal environment is subjected to the uniformly distributed loads, in which the beam consists of certain bimodular materials, and it will deflect downward to resist the external loads, thus resulting in tension in the lower part of the beam and compression in the upper part. In Figure 3, the span length of the beam is 2l, the height is h and the thickness is t and is taken as unit 1; q is the magnitude of the uniformly distributed loads and ql denotes the reaction of two supports of the beam. Note that in such a bimodular beam, the tensile height of the beam is h 1 , the compressive height h 2 ; accordingly, the tensile modulus is E + in the tensile part and the compressive modulus is E in the compressive part. For convenience, we establish the x-axis of the coordinate system (xyz) on the known neutral axis, which was determined in Section 2.
If the two plane dimensions of the beam, 2l and h, are in the same order of magnitude, it is hard to obtain the analytical solution of this problem. However, if 2 l > > h , the left and right sides of the beam become second boundary, thus making the application of the de Saint-Venant’s Principle rational. For this purpose, the following stress function is adopted
φ ( x , y ) = x 2 2 f ( y ) + x f 1 ( y ) + f 2 ( y ) ,
in which, f ( y ) , f 1 ( y ) and f 2 ( y ) are three unknown function of y. The obtainment of this stress function is based on the so-called semi-inverse method, which may be simply described as follows. First, according to the real stress condition, the beam is subjected to the uniformly distributed loads along the direction of span length, and it is assumed that the stress component along y direction does not change with x, thus beginning with σ y = f ( y ) . Using σ y = 2 φ / 2 x and integrating twice with respect to x, Equation (29) may be derived. By satisfying the consistency equation 4 φ ( x , y ) = 0 , in which 4 denotes the dual Laplace operator, as well as after considering the symmetry of this problem, that is, f 1 ( y ) = 0 , lastly the stress function expressed in terms of unknown constants may be obtained,
φ ( x , y ) = x 2 2 ( A y 3 + B y 2 + C y + D ) A 10 y 5 B 6 y 4 + H y 3 + K y 2 ,
in which, A, B, C, D, H and K are six unknown constants that may be determind by boundary conditions. Note that the mechanical model on subarea in tension and compression is adopted here, the stress components are given in the tensile part and the compressive one, respectively. Substituting Equation (30) into Equation (28), we have the stress components expressed in terms of unknown constants, for the tensile part, 0 y h 1
{ σ x + = x 2 2 ( 6 A + y + 2 B + ) 2 A + y 3 2 B + y 2 + 6 H + y + 2 K + σ y + = A + y 3 + B + y 2 + C + y + D + τ x y + = x ( 3 A + y 2 + 2 B + y + C + ) ,
and for the compressive part, h 2 y 0 ,
{ σ x = x 2 2 ( 6 A y + 2 B ) 2 A y 3 2 B y 2 + 6 H y + 2 K σ y = A y 3 + B y 2 + C y + D τ x y = x ( 3 A y 2 + 2 B y + C ) .  
Now, superposing with the special solution (27), we have, for 0 y h 1 ,
{ σ x + = α E T ( y ) + x 2 2 ( 6 A + y + 2 B + ) 2 A + y 3 2 B + y 2 + 6 H + y + 2 K + σ y + = A + y 3 + B + y 2 + C + y + D + τ x y + = x ( 3 A + y 2 + 2 B + y + C + ) ,
and for h 2 y 0 ,
{ σ x = α E T ( y ) + x 2 2 ( 6 A y + 2 B ) 2 A y 3 2 B y 2 + 6 H y + 2 K σ y = A y 3 + B y 2 + C y + D τ x y = x ( 3 A y 2 + 2 B y + C ) ,
in which the quantities with superscript ‘+’ denote the quantities relating to tension and the quantities with superscript ‘−’ denote the quantities relating to compression. Note that after the subarea in tension and compression, the number of undetermined constants is double, from the original six to twelve. Obviously, not only stress boundary conditions but also certain continuity conditions need to be combined to determine these unknown constants.
These total stresses should satisfy the following boundary conditions: first, at the bottom of the beam,
σ y + = 0   τ x y + = 0   at   y = h 1 ,
and at the top of the beam,
σ y = q   τ x y = 0   at   y = h 2 .
Due to 2 l > > h indicated above, the left and right sides become the secondary conditions, thus the de Saint-Venant’s Principle may be used to formulate the boundary conditions at two ends of the beam. At x = ± l , three boundary conditions with integrals are as follows
h 2 0 σ x t d y + 0 h 1 σ x + t d y = 0 ,
h 2 0 σ x t y d y + 0 h 1 σ x + t y d y = 0 ,
h 2 0 τ x y t d y + 0 h 1 τ x y + t d y = q l .
Obviously, the above boundary conditions are insufficient for solving the twelve unknown constants. Therefore, the continuity conditions at the neutral axis must be used, which gives
σ y + = σ y   τ x y + = τ x y   at   y = 0 ,
and the normal stress along x direction should be zero, that is,
σ x + = σ x = 0   at   y = 0 .
Eventually, the twelve unknown constants may be determined by the boundary conditions (35)–(39), in combination with the continuity conditions (40) and (41).
Now, let us begin with the application of continuity conditions. First, substituting the last two expressions of Equations (33) and (34) into Equation (40), we have
C + = C ,   D + = D .
Substituting the first expression of Equations (33) and (34) into Equation (41), we have
σ x + = B + x 2 + 2 K + = 0 ,   σ x = B x 2 + 2 K = 0 .
It should be noted here that the neutral axis of the beam is determined only under the action of mechanical load, which has nothing to do with the thermal load, thus the term α E T ( y ) in σ x + and σ x should not be included here. Thus, to satisfy (43) for all x, we obtain
B + = B = K + = K = 0 .
It is found that after the application of all continuity conditions, we obtain six relations concerning the twelve unknown constants, as shown in Equations (42) and (44).
Substituting the last two expressions of Equations (33) and (34) into Equations (35) and (36), we have
A + h 1 3 + B + h 1 2 + C + h 1 + D + = 0 ,
3 A + h 1 2 + 2 B + h 1 + C + = 0 ,
A h 2 3 + B h 2 2 C h 2 + D = q ,
and
3 A h 2 2 2 B h 2 + C = 0 .
Substituting the last two expressions of Equations (33) and (34) into Equations (37) and (38), we have
h 2 0 [ α E T ( y ) + l 2 2 ( 6 A y + 2 B ) 2 A y 3 2 B y 2 + 6 H y + 2 K ] d y + 0 h 1 [ α E T ( y ) + l 2 2 ( 6 A + y + 2 B + ) 2 A + y 3 2 B + y 2 + 6 H + y + 2 K + ] d y = 0 ,
and
h 2 0 [ α E T ( y ) + l 2 2 ( 6 A y + 2 B ) 2 A y 3 2 B y 2 + 6 H y + 2 K ] y d y + 0 h 1 [ α E T ( y ) + l 2 2 ( 6 A + y + 2 B + ) 2 A + y 3 2 B + y 2 + 6 H + y + 2 K + ] y d y = 0 .
Using Equations (45)–(50) in combination of Equations (42) and (44), we obtain
A + = q 2 h h 1 2 , A = q 2 h h 2 2 , B + = B = 0 , C + = C = 3 q 2 h , D + = D = q h 1 h , H + = J 1 h 2 3 h h 1 2 + J 2 2 h h 1 2 + q l 2 4 h h 1 2 q N 1 10 h , H = J 1 h 1 3 h h 2 2 + J 1 2 h h 2 2 + q l 2 4 h h 2 2 q N 2 10 h , K + = K = 0 ,
in which, J 1 and J 2 are two integral terms about T ( y ) , N 1 and N 2 are two known constants, they are, respectively,
J 1 = h 2 0 [ α E T ( y ) ] d y + 0 h 1 [ α E T ( y ) ] d y J 2 = h 2 0 [ α E T ( y ) ] y d y + 0 h 1 [ α E T ( y ) ] y d y ,
and
N 1 = 1 1 6 h 2 h 1 + 1 6 ( h 2 h 1 ) 2 , N 2 = 1 1 6 h 1 h 2 + 1 6 ( h 1 h 2 ) 2 .
Up to now, Equation (39) has not been used. Substituting the known solution into Equation (39), it is found that it is naturally satisfied. Lastly, substituting the known constants into Equations (33) and (34), we have the final stress for the tensile part 0 y h 1 ,
{ σ x + = α E T ( y ) + 3 q 2 h h 1 2 ( l 2 x 2 ) y + q h y ( y 2 h 1 2 3 5 N 1 ) + ( 2 J 1 h 2 h h 1 2 + 3 J 2 h h 1 2 ) y σ y + = q 2 h h 1 2 y 3 + 3 q 2 h y q h 1 h τ x y + = 3 q 2 h h 1 2 x y 2 3 q 2 h x
and for the compressive part h 2 y 0 ,
{ σ x = α E T ( y ) + 3 q 2 h h 2 2 ( l 2 x 2 ) y + q h y ( y 2 h 2 2 3 5 N 2 ) + ( 2 J 1 h 1 h h 2 2 + 3 J 2 h h 2 2 ) y σ y = q 2 h h 2 2 y 3 + 3 q 2 h y q h 1 h τ x y = 3 q 2 h h 2 2 x y 2 3 q 2 h x

6. Results and Discussions

6.1. A Computational Example

If the function of the temperature rise is known, for example, we assume the highest temperature rise takes place at the middle of the beam and the lowest temperature rise at the up and down edges, this temperature rise may be simulated with a function as
T ( y ) = T 0 cos y y 0 h π ,
in which T 0 is the temperature at the center and y 0 is the distance of the geometric center of the beam from the neutral axis, see Figure 2. Via Equation (7) and h / 2 + y 0 = h 1 (or h / 2 y 0 = h 2 ), the magnitude of y 0 may be determined as
y 0 = E E + E + E + h 2 .
Substituting the known T ( y ) into Equation (52), the two integrals may be computed as
J 1 = α E T 0 h π ( sin h 2 + y 0 h π + sin h 1 y 0 h π ) J 2 = α E T 0 h π 2 ( π h 2 sin h 2 + y 0 h π π h 1 sin h 1 y 0 h π + h cos h 2 + y 0 h π h cos h 1 y 0 h π ) .
Thus, the normal stress along x direction is, for 0 y h 1 ,
σ x + = α E T 0 cos y y 0 h π + 3 q 2 h h 1 2 ( l 2 x 2 ) y + q h y ( y 2 h 1 2 3 5 N 1 ) + α E T 0 [ 2 h 2 π h 1 2 ( sin h 2 + y 0 h π + sin h 1 y 0 h π ) + 3 π 2 h 1 2 ( π h 2 sin h 2 + y 0 h π π h 1 sin h 1 y 0 h π + h cos h 2 + y 0 h π h cos h 1 y 0 h π ) ] y .
and for h 2 y 0 ,
σ x = α E T 0 cos y y 0 h π + 3 q 2 h h 2 2 ( l 2 x 2 ) y + q h y ( y 2 h 2 2 3 5 N 2 ) + α E T 0 [ 2 h 1 π h 2 2 ( sin h 2 + y 0 h π + sin h 1 y 0 h π ) + 3 π 2 h 2 2 ( π h 2 sin h 2 + y 0 h π π h 1 sin h 1 y 0 h π + h cos h 2 + y 0 h π h cos h 1 y 0 h π ) ] y .
Note that there are four terms in the above stress expression, in which the first and the fourth are from the temperature stress, and the second and the third are from mechanical stress. In the mechanical stress, the second stands for the stress from the bending moment in a one-dimensional solution, while the third stands for the rectified term from a two-dimensional elasticity solution.
When E + = E = E , we have h 1 = h 2 = h / 2 , y 0 = 0 , J 1 = 2 E α T 0 h / π , J 2 = 0 and N 1 = N 2 = 1 . Substituting them into Equations (59) and (60), the stress expression without bimodular effect is obtained as follows
σ x = α E T 0 ( 8 π y h + cos y h π ) + 6 q h 3 ( l 2 x 2 ) y + q y h ( 4 y 2 h 2 3 5 ) .
Thus, the regression of the solution is satisfied, which also verifies the correctness of the derivation from the side.

6.2. Bimodular Effect on Thermoelasticity Solution

For analyzing the bimodular effect on the final stress σ x , the following three typical cases relating to modulus are adopted: (1) E + = 2 E , (2) E + = E and (3) E + = 0.5 E , which correspond to, respectively, the tensile modulus is greater than, equal to and less than the compressive modulus. Thus, the tensile modulus and the compressive one, as well as the tensile section height, the compressive one and the distance between the geometrical middle surface and the neutral layer, are given in Table 1 in dimensionless form.
For a better understanding of the effects of the parameters on the stress, all figures in this section are plotted in dimensionless form, in which the x-axis stands for the dimensionless stress σ x / α E T 0 and the y-axis for the dimensionless section height y / h . On the cross section of the beam midspan, the bending moment achieves the most, thus we take this section x = 0 as our analysis object. In addition, it is also necessary to give the relative magnitudes between α E T 0 and the intensity q of the uniformly distributed loads, thus after considering the real case, we take the following three ratios: (1) α E T 0 / q = 20 , (2) α E T 0 / q = 30 and (3) α E T 0 / q = 40 . At the same time, the ratio of beam span to beam height is assumed to be 6, that is, 2 l / h = 6 , to agree with the geometrical characteristic of a two-dimensional deep beam.
Figure 4, Figure 5 and Figure 6 show the stress composition from the temperature stress, the mechanical stress, and total stress under three bimodular cases, in which only the curves when α E T 0 / q = 20 are given here, and the other two cases, that is, α E T 0 / q = 30 and α E T 0 / q = 40 , are not presented here since they are basically the same as Figure 4, Figure 5 and Figure 6.
From Figure 4, Figure 5 and Figure 6, it is easy to see that
(i) Under three cases of different moduli, the mechanical stresses all present as obviously linear, although the stress σ x is the cubic relation with respect to y; this indicates that the second term in the stress expressions is dominant, compared to the third, whose influence on the mechanical stress is relatively small.
(ii) When E + = E , the slopes of the straight line are consistent in the tensile and compressive parts since the stress growth ratio is constant with the strain (which is proportional to the distance along the y-axis); when E + = 2 E , the slope of the straight line in the compressive part is greater than that in the tensile part, showing that when E + > E , the stress growth has slowed down; when E + = 0.5 E , the slope of the straight line in the compressive part is smaller than that in the tensile part, showing that when E + < E , the stress growth has accelerated.
(iii) The temperature stress curves basically present the shape of a cosine function, having no up and down symmetry. After superposition with the mechanical stress curves, which are linear, the total stress curves still present the shape of a cosine function, only having their vertices somewhat offset.
Figure 7, Figure 8 and Figure 9 show that the total stress when α E T 0 / q = 20 , α E T 0 / q = 30 and α E T 0 / q = 40 varies with the section height direction, in which three bimodular cases are given.
From Figure 7, Figure 8 and Figure 9, it is easy to see that
(i) Under three cases of different moduli, and also for α E T 0 / q = 20 , α E T 0 / q = 30 and α E T 0 / q = 40 , the total stress curves all present the shape of a cosine function, having no up and down symmetry, and their vertices are at different places.
(ii) With the decrease of q, or alternatively, with the increase of α E T 0 , that is, from Figure 7, Figure 8 and Figure 9, we may find that when α E T 0 / q = 20 , there are some tensile stresses in the lower part of the beam (see Figure 7), when α E T 0 / q = 30 , there are only small quantities of tensile stresses remaining (see Figure 8), and when α E T 0 / q = 40 , there are no any tensile stresses in the lower part (see Figure 9). The tensile stress situation in the upper part of the beam is just the opposite.
(iii) Under three cases of different moduli, and also for α E T 0 / q = 20 , α E T 0 / q = 30 and α E T 0 / q = 40 , most of the total stresses are in compression, this reveals that after the introduction of the temperature stress, the original stress state with tension in the lower part and compression in the upper part, now gradually transits to the status almost in compression. This phenomenon is especially obvious for the case E + = 2 E .

6.3. Temperature Effect on Thermoelasticity Solution

For analyzing the temperature effect on the thermoelasticity solution, the three modes of temperature rise, including linear, square, and cubic relations, are taken to discuss the different modes of temperature rise on the thermal stress, in which the maximum temperature rise keeps T 0 in the three modes. First, for the linear temperature rise, we have
T ( y ) = T 0 2 [ 1 + 2 ( y y 0 ) h ] .
Substituting the known T ( y ) into Equation (52), the two integrals give
J 1 = α E T 0 2 ( h 1 h 2 2 + h 2 y 0 ) J 2 = α E T 0 2 [ 2 ( h 1 3 + h 2 3 ) 3 h + h 1 2 h 2 2 2 ( 1 2 y 0 h ) ] .
Thus, the normal stress along x direction is, for 0 y h 1
σ x + = α E T 0 2 [ 1 + 2 ( y y 0 ) h ] + 3 q 2 h h 1 2 ( l 2 x 2 ) y + q h y ( y 2 h 1 2 3 5 N 1 ) α E T 0 2 { 2 h 2 h h 1 2 ( h 1 h 2 2 + h 2 y 0 ) + 3 h h 1 2 [ 2 ( h 1 3 + h 2 3 ) 3 h + h 1 2 h 2 2 2 ( 1 2 y 0 h ) ] } y ,
and for h 2 y 0
σ x = α E T 0 2 [ 1 + 2 ( y y 0 ) h ] + 3 q 2 h h 2 2 ( l 2 x 2 ) y + q h y ( y 2 h 2 2 3 5 N 2 ) α E T 0 2 { 2 h 1 h h 2 2 ( h 1 h 2 2 + h 2 y 0 ) + 3 h h 2 2 [ 2 ( h 1 3 + h 2 3 ) 3 h + h 1 2 h 2 2 2 ( 1 2 y 0 h ) ] } y .
Similarly, for the square temperature rise, we have
T ( y ) = T 0 [ 1 - 4 ( y y 0 ) 2 h 2 ] ,
J 1 = α E T 0 [ 4 ( h 1 2 + h 2 3 ) 3 h 2 + y 0 ( h 1 h 2 ) h + h 4 y 0 2 h ] J 2 = α E T 0 [ h 1 4 h 2 4 h 2 + 8 y 0 ( h 1 3 + h 2 3 ) 3 h 2 + h 1 2 h 2 2 2 ( 1 4 y 0 2 h 2 ) ] .
Thus, the normal stress along x direction is, for 0 y h 1
σ x + = α E T 0 [ 1 4 ( y y 0 ) 2 h 2 ] + 3 q 2 h h 1 2 ( l 2 x 2 ) y + q h y ( y 2 h 1 2 3 5 N 1 ) α E T 0 { 2 h 2 h h 1 2 [ 4 ( h 1 3 + h 2 3 ) 3 h 2 + y 0 ( h 1 h 2 ) h + h 4 y 0 2 h ] + 3 h h 1 2 [ h 1 4 h 2 4 h 2 + 8 y 0 ( h 1 3 + h 2 3 ) 3 h 2 + h 1 2 h 2 2 2 ( 1 4 y 0 2 h 2 ) ] } y ,
and for h 2 y 0 ,
σ x = α E T 0 [ 1 4 ( y y 0 ) 2 h 2 ] + 3 q 2 h h 2 2 ( l 2 x 2 ) y + q h y ( y 2 h 2 2 3 5 N 2 ) α E T 0 { 2 h 1 h h 2 2 [ 4 ( h 1 3 + h 2 3 ) 3 h 2 + y 0 ( h 1 h 2 ) h + h 4 y 0 2 h ] + 3 h h 2 2 [ h 1 4 h 2 4 h 2 + 8 y 0 ( h 1 3 + h 2 3 ) 3 h 2 + h 1 2 h 2 2 2 ( 1 4 y 0 2 h 2 ) ] } y .
For the cubic temperature rise, we have
T ( y ) = T 0 2 [ 1 + 8 ( y y 0 ) 3 h 3 ] ,
J 1 = α E T 0 2 [ 2 ( h 1 4 h 2 4 ) h 3 8 y 0 ( h 1 3 + h 2 3 ) h 3 + 12 y 0 2 ( h 1 h 2 ) h 2 + h 8 y 0 3 h 2 ] , J 2 = α E T 0 2 [ 8 ( h 1 5 + h 2 5 ) 5 h 3 6 y 0 ( h 1 4 h 2 4 ) h 3 + 8 y 0 2 ( h 1 3 + h 2 3 ) h 3 + h 1 2 h 2 2 2 ( 1 8 y 0 3 h 3 ) ] ,
Thus, the normal stress along x direction is, for 0 y h 1 ,
σ x + = α E T 0 2 [ 1 + 8 ( y y 0 ) 3 h 3 ] + 3 q 2 h h 1 2 ( l 2 x 2 ) y + q h y ( y 2 h 1 2 3 5 N 1 ) α E T 0 2 { 2 h 2 h h 1 2 [ 2 ( h 1 4 h 2 4 ) h 3 8 y 0 ( h 1 3 + h 2 3 ) h 3 + 12 y 0 2 ( h 1 h 2 ) h 2 + h 8 y 0 3 h 2 ] + 3 h h 1 2 [ 8 ( h 1 5 + h 2 5 ) 5 h 3 6 y 0 ( h 1 4 h 2 4 ) h 3 + 8 y 0 2 ( h 1 3 + h 2 3 ) h 3 + h 1 2 h 2 2 2 ( 1 8 y 0 3 h 3 ) ] } y
and for h 2 y 0
σ x = α E T 0 2 [ 1 + 8 ( y y 0 ) 3 h 3 ] + 3 q 2 h h 2 2 ( l 2 x 2 ) y + q h y ( y 2 h 2 2 3 5 N 2 ) α E T 0 2 { 2 h 1 h h 2 2 [ 2 ( h 1 4 h 2 4 ) h 3 8 y 0 ( h 1 3 + h 2 3 ) h 3 + 12 y 0 2 ( h 1 h 2 ) h 2 + h 8 y 0 3 h 2 ] + 3 h h 2 2 [ 8 ( h 1 5 + h 2 5 ) 5 h 3 6 y 0 ( h 1 4 h 2 4 ) h 3 + 8 y 0 2 ( h 1 3 + h 2 3 ) h 3 + h 1 2 h 2 2 2 ( 1 8 y 0 3 h 3 ) ] } y .
Taking the case α E T 0 / q = 20 as an example, the total thermal stress curves are plotted under the three modes of temperature rise, as shown in Figure 10, Figure 11 and Figure 12, which correspond to the three modulus cases, E + = 2 E , E + = E and E + = 0.5 E , respectively. Figure 10, Figure 11 and Figure 12 show that the variation trends of the three thermal stress curves are basically the same as those from the three modes of temperature rise, in which the curve from the square relation is basically consistent with the curve from the cosine relation, with a vertex but no up and down symmetry.
Table 2 lists dimensionless thermal stresses for the linear, square, and cubic rise modes as well as the variations of amplitude value throughout the beam. The computation of amplitude value is, for the linear and cubic relation, from the top to bottom of the beam; and for the square relation, is from the maximum tensile stress to maximum compressive one. It is easy to see from Table 2 that, under the temperature rise mode of square relation, the amplitude is 2.21 for E + = 2 E , versus 1.80 for E + = E , and decreases to 1.58 for E + = 0.5 E ; under the temperature rise mode of linear and cubic relations, the case is just the opposite. Under the linear rise, the amplitude is 0.28 for E + = 2 E , versus 1.28 for E + = E , and increases to 2.43 for E + = 0.5 E ; and under the cubic rise, the amplitude is 0.41 for E + = 2 E , versus 0.88 for E + = E , and increases to 1.47 for E + = 0.5 E . This indicates that if the highest temperature rise takes place in the middle of the beam, like the square relation, more attention should be given when E + > E , and if the highest temperature rise takes place at the bottom of the beam, like the linear and cubic relations, more attention should be given when E + < E .

7. Concluding Remarks

In this study, based on the simplified mechanical model on regional segmentation of tension and compression, we obtained a two-dimensional thermoelasticity solution of a bimodular rectangular section beam under the combination action of thermal and mechanical loads. The following conclusions can be drawn.
(i) After the introduction of the bimodular effect, there is no essential change in the form of displacement governing the equation and the composition of the solution since the modulus of elasticity is not included in Lamé’s displacement equation.
(ii) The stress superposition of special solution and supplement solution shows clearly the stress composition contributed by the thermal load and mechanical load, respectively. Which one is dominant will depend on the relative magnitude between the external load q and the thermal stress term α E T 0 .
(iii) Different temperature rise modes will influence the distribution of thermal stresses. The analysis of the amplitude values of thermal stresses indicates that if the highest temperature rise takes place in the middle of the beam, the amplitude varies from 2.21 ( E + = 2 E ) to 1.58 ( E + = 0.5 E ); while if the highest temperature rise takes place at the bottom of the beam, the case is just the opposite. Specifically, the amplitude varies from 0.28 to 2.43 for the linear rise and from 0.41 to 1.47 for the cubic rise.
(iv) The introduction of bimodular characteristics will bring great change in the total stress. More specifically, when the tensile modulus is greater than the compressive one, most of the total stresses are in compression, which should be give more consideration in the analysis and design of bimodular material beams.
In addition, this study may be further extended into the field of functionally graded materials. In this case, however, the obtained solution of the Lamé displacement equation will be somewhat challenging, due to the fact that not only the modulus of elasticity but also some redundant strain terms will appear in this equation. The changes of the governing equation and the investigation into the corresponding solving methods (for example, the application of the multi-parameter perturbation method [30]) will present us with huge challenges. At the same time, in a thermal environment, the investigation into bimodular functionally graded piezoelectric beams [31] as well as into carbon nanotube-reinforced composite curved sandwich nanobeams [32] points out our future research direction. The relevant work is in progress.

Author Contributions

Conceptualization, X.-T.H. and J.-Y.S.; methodology, X.-T.H. and S.-R.W.; formal analysis, X.-T.H. and S.-R.W.; writing—original draft preparation, X.-T.H. and S.-R.W.; writing—review and editing, H.C. and J.-Y.S.; visualization, H.C. and S.-R.W.; funding acquisition, X.-T.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 11572061).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Timoshenko, S.P.; Goodier, J.N. Theory of Elasticity, 3rd ed.; McGraw Hill: New York, NY, USA, 1970. [Google Scholar]
  2. Barak, M.M.; Currey, J.D.; Weiner, S.; Shahar, R. Are tensile and compressive Young’s moduli of compact bone different? J. Mech. Behav. Biomed. Mater. 2009, 2, 51–60. [Google Scholar] [CrossRef] [PubMed]
  3. Destrade, M.; Gilchrist, M.D.; Motherway, J.A.; Murphy, J.G. Bimodular rubber buckles early in bending. Mech. Mater. 2010, 42, 469–476. [Google Scholar] [CrossRef] [Green Version]
  4. Jones, R.M. Stress-strain relations for materials with different moduli in tension and compression. AIAA J. 1977, 15, 16–23. [Google Scholar] [CrossRef]
  5. Bert, C.W. Models for Fibrous Composites with Different Properties in Tension and Compression. J. Eng. Mater. Technol. 1977, 99, 344–349. [Google Scholar] [CrossRef]
  6. Bruno, D.; Lato, S.; Sacco, E. Nonlinear analysis of bimodular composite plates under compression. Comput. Mech. 1994, 14, 28–37. [Google Scholar] [CrossRef]
  7. Tseng, Y.P.; Lee, C.T. Bending analysis of bimodular laminates using a higher-order finite strip method. Compos. Struct. 1995, 30, 341–350. [Google Scholar] [CrossRef]
  8. Zinno, R.; Greco, F. Damage evolution in bimodular laminated composites under cyclic loading. Compos. Struct. 2001, 53, 381–402. [Google Scholar] [CrossRef]
  9. Khan, A.H.; Patel, B.P. Nonlinear periodic response of bimodular laminated composite annular sector plates. Compos. Part B Eng. 2019, 169, 96–108. [Google Scholar] [CrossRef]
  10. Ambartsumyan, S.A. Elasticity Theory of Different Moduli (Wu RF and Zhang YZ Trans.); China Railway Publishing House: Beijing, China, 1986. [Google Scholar]
  11. Wen-Juan, Y.; Zhi-ming, Y. Analytical solution for bending beam subject to lateral force with different modulus. Appl. Math. Mech. 2004, 25, 1107–1117. [Google Scholar] [CrossRef]
  12. He, X.T.; Chen, S.L.; Sun, J.Y. Applying the equivalent section method to solve beam subjected lateral force and bending-compression column with different moduli. Int. J. Mech. Sci. 2007, 49, 919–924. [Google Scholar] [CrossRef]
  13. He, X.T.; Sun, J.Y.; Wang, Z.X.; Chen, Q.; Zheng, Z.L. General perturbation solution of large-deflection circular plate with different moduli in tension and compression under various edge conditions. Int. J. Non-Linear Mech. 2013, 55, 110–119. [Google Scholar] [CrossRef]
  14. He, X.T.; Cao, L.; Wang, Y.Z.; Sun, J.Y.; Zheng, Z.L. A biparametric perturbation method for the Föppl-von Kármán equations of bimodular thin plates. J. Math. Anal. Appl. 2017, 455, 1688–1705. [Google Scholar] [CrossRef]
  15. Zhiming, Y.E.; Wenjuan, C.T.Y. Progresses in elasticity theory with different moduli in tension and compression and related FEM. Mech. Eng. 2004, 26, 9–14. [Google Scholar]
  16. Sun, J.Y.; Zhu, H.Q.; Qin, S.H.; Yang, D.L.; He, X.T. A review on the research of mechanical problems with different moduli in tension and compression. J. Mech. Sci. Technol. 2010, 24, 1845–1854. [Google Scholar] [CrossRef]
  17. Du, Z.; Zhang, Y.; Zhang, W.; Guo, X. A new computational framework for materials with different mechanical responses in tension and compression and its applications. Int. J. Solids Struct. 2016, 100, 54–73. [Google Scholar] [CrossRef]
  18. Gao, J.L.; Yao, W.J.; Liu, J.K. Temperature stress analysis for bi-modulus beam placed on Winkler foundation. Appl. Math. Mech. 2017, 38, 921–934. [Google Scholar] [CrossRef]
  19. Ma, J.W.; Fang, T.C.; Yao, W.J. Nonlinear large deflection buckling analysis of compression rod with different moduli. Mech. Adv. Mater. Struct. 2019, 26, 539–551. [Google Scholar] [CrossRef]
  20. Hetnarski, R.B.; Eslami, M.R. Thermal Stresses—Advanced Theory and Applications; Springer: Berlin, Germany, 2019. [Google Scholar]
  21. Green, A.E.; Lindsay, K.A. Thermoelasticity. J. Elast. 1972, 2, 1–7. [Google Scholar] [CrossRef]
  22. Green, A.E.; Naghdi, P.M. On Undamped Heat Waves in an Elastic Solid. J. Therm. Stress. 1992, 15, 253–264. [Google Scholar] [CrossRef]
  23. Choudhuri, S.K.R. On A Thermoelastic Three-Phase-Lag Model. J. Therm. Stress. 2007, 30, 231–238. [Google Scholar] [CrossRef]
  24. Svanadze, M.; Scalia, A. Mathematical problems in the coupled linear theory of bone poroelasticity. Comput. Math. Appl. 2013, 66, 1554–1566. [Google Scholar] [CrossRef]
  25. Marin, M.; Craciun, E.M.; Pop, N. Some Results in Green–Lindsay Thermoelasticity of Bodies with Dipolar Structure. Mathematics 2020, 8, 497. [Google Scholar] [CrossRef] [Green Version]
  26. Abouelregal, A.E.; Marin, M. The Size-Dependent Thermoelastic Vibrations of Nanobeams Subjected to Harmonic Excitation and Rectified Sine Wave Heating. Mathematics 2020, 8, 1128. [Google Scholar] [CrossRef]
  27. Abouelregal, A.E.; Zenkour, A.M. Thermoelastic problem of an axially moving microbeam subjected to an external transverse excitation. J. Theor. Appl. Mech. 2015, 53, 167. [Google Scholar] [CrossRef] [Green Version]
  28. Warminska, A.; Manoach, E.; Warminski, J. Vibrations of a Composite Beam under Thermal and Mechanical Loadings. Procedia Eng. 2016, 144, 959–966. [Google Scholar] [CrossRef] [Green Version]
  29. Tao, C.; Fu, Y.M.; Dai, H.L. Nonlinear dynamic analysis of fiber metal laminated beams subjected to moving loads in thermal environment. Compos. Struct. 2016, 140, 410–416. [Google Scholar] [CrossRef] [Green Version]
  30. He, X.T.; Yang, Z.X.; Li, Y.H.; Li, X.; Sun, J.Y. Application of Multi-Parameter Perturbation Method to Functionally-Graded, Thin, Circular Piezoelectric Plates. Mathematics 2020, 8, 342. [Google Scholar] [CrossRef] [Green Version]
  31. Jing, H.X.; He, X.T.; Du, D.W.; Peng, D.D.; Sun, J.Y. Vibration Analysis of Piezoelectric Cantilever Beams with Bimodular Functionally-Graded Properties. Appl. Sci. 2020, 10, 5557. [Google Scholar] [CrossRef]
  32. Daikh, A.A.; Houari, M.S.A.; Karami, B.; Eltaher, M.A.; Dimitri, R.; Tornabene, F. Buckling Analysis of CNTRC Curved Sandwich Nanobeams in Thermal Environment. Appl. Sci. 2021, 11, 3250. [Google Scholar] [CrossRef]
Figure 1. Bimodulus materials model proposed by Ambartsumyan.
Figure 1. Bimodulus materials model proposed by Ambartsumyan.
Mathematics 09 01556 g001
Figure 2. The mechanical model on subarea in tension and compression under mechanical load.
Figure 2. The mechanical model on subarea in tension and compression under mechanical load.
Mathematics 09 01556 g002
Figure 3. A bimodular beam under uniformly distributed loads in a thermal environment.
Figure 3. A bimodular beam under uniformly distributed loads in a thermal environment.
Mathematics 09 01556 g003
Figure 4. Stress composition when E + = 2 E and α E T 0 / q = 20 .
Figure 4. Stress composition when E + = 2 E and α E T 0 / q = 20 .
Mathematics 09 01556 g004
Figure 5. Stress composition when E + = E and α E T 0 / q = 20 .
Figure 5. Stress composition when E + = E and α E T 0 / q = 20 .
Mathematics 09 01556 g005
Figure 6. Stress composition when E + = 0.5 E and α E T 0 / q = 20 .
Figure 6. Stress composition when E + = 0.5 E and α E T 0 / q = 20 .
Mathematics 09 01556 g006
Figure 7. The total stress when α E T 0 / q = 20 .
Figure 7. The total stress when α E T 0 / q = 20 .
Mathematics 09 01556 g007
Figure 8. The total stress when α E T 0 / q = 30 .
Figure 8. The total stress when α E T 0 / q = 30 .
Mathematics 09 01556 g008
Figure 9. The total stress when α E T 0 / q = 40 .
Figure 9. The total stress when α E T 0 / q = 40 .
Mathematics 09 01556 g009
Figure 10. The total stress when E + = 2 E and α E T 0 / q = 20 .
Figure 10. The total stress when E + = 2 E and α E T 0 / q = 20 .
Mathematics 09 01556 g010
Figure 11. The total stress when E + = E and α E T 0 / q = 20 .
Figure 11. The total stress when E + = E and α E T 0 / q = 20 .
Mathematics 09 01556 g011
Figure 12. The total stress when E + = 0.5 E and α E T 0 / q = 20 .
Figure 12. The total stress when E + = 0.5 E and α E T 0 / q = 20 .
Mathematics 09 01556 g012
Table 1. Values for three bimodular cases.
Table 1. Values for three bimodular cases.
Bimodular CasesTensile ModulusCompressive Modulus h 1 / h h 2 / h y 0 / h
E + > E 2 E E 0.41420.5858−0.0858
E + = E E E 0.50.50
E + < E 0.5 E E 0.58580.41420.0858
Table 2. Dimensionless thermal stresses and amplitude values.
Table 2. Dimensionless thermal stresses and amplitude values.
Bimodular CasesLinearSquareCubic
Top
Stress
Bottom
Stress
Amplit.
Value 1
Max Tens. StressMax Compr. StressAmplit.
Value
Top
Stress
Bottom
Stress
Amplit.
Value
E + = 2 E −0.43−0.710.280.96−1.252.21−0.42−0.830.41
E + = E 0.14−1.141.280.69−1.111.80−0.06−0.940.88
E + = 0.5 E 1.06−1.372.430.58−1.001.580.45−1.021.47
1 If two stresses are the same positive–negative sign, take the absolute value after their subtraction; if they are the contrary sign, take the sum of the two abstract values of two stresses.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wen, S.-R.; He, X.-T.; Chang, H.; Sun, J.-Y. A Two-Dimensional Thermoelasticity Solution for Bimodular Material Beams under the Combination Action of Thermal and Mechanical Loads. Mathematics 2021, 9, 1556. https://doi.org/10.3390/math9131556

AMA Style

Wen S-R, He X-T, Chang H, Sun J-Y. A Two-Dimensional Thermoelasticity Solution for Bimodular Material Beams under the Combination Action of Thermal and Mechanical Loads. Mathematics. 2021; 9(13):1556. https://doi.org/10.3390/math9131556

Chicago/Turabian Style

Wen, Si-Rui, Xiao-Ting He, Hao Chang, and Jun-Yi Sun. 2021. "A Two-Dimensional Thermoelasticity Solution for Bimodular Material Beams under the Combination Action of Thermal and Mechanical Loads" Mathematics 9, no. 13: 1556. https://doi.org/10.3390/math9131556

APA Style

Wen, S. -R., He, X. -T., Chang, H., & Sun, J. -Y. (2021). A Two-Dimensional Thermoelasticity Solution for Bimodular Material Beams under the Combination Action of Thermal and Mechanical Loads. Mathematics, 9(13), 1556. https://doi.org/10.3390/math9131556

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop