Fiber-reinforced materials have widespread applications, which prompt the study of the effect of fiber reinforcement. Research studies have indicated that thermal conductivity cannot be considered as a constant, which is closely related to temperature change. Based on those studies, we investigate the fiber-reinforced generalized thermoelasticity problem under thermal stress, with the consideration of the effect of temperature-dependent variable thermal conductivity. The problem is assessed according to the L-S theory. A fiber-reinforced anisotropic half-space is selected as the research model, and a region of its surface is subjected to a transient thermal shock. The time-domain finite element method is applied to analyze the nonlinear problem and derives the governing equations. The nondimensional displacement, stress, and temperature of the material are obtained and illustrated graphically. The numerical results reveal that the variable conductivity significantly influences the distribution of the field quantities under the fiber-reinforced effect. And also, the boundary point of thermal shock is the most affected. The obtained results in this paper can be applied to design the fiber-reinforced anisotropic composites under thermal load to satisfy some particular engineering requirements.

1. Introduction

Fiber reinforcement is an inherent property of materials that is considered an effect rather than a form of inclusion in it [1]. In an elastic state, the components of fiber-reinforced composites act as a single anisotropic unit without relative displacement [26]. Fiber-reinforced material produces higher specific strength and a larger specific modulus in the direction of fiber reinforcement. The material performance of fiber-reinforced composites is designable, and its corrosion resistance and durability are good. Those outstanding features inherent in fiber-reinforced composites have led to their widespread applications in aerospace, building engineering, automotive industries, and so on [7, 8]. Therefore, the fiber-reinforced effect of materials should be considered when studying mechanical behavior.

The assumption of infinite propagation speed of the thermal signal in classical thermoelasticity theory is inconsistent with the real phenomenon. Several generalized thermoelasticity theories have been developed to eliminate this paradox [912], such as L-S theory. It firstly used the Maxwell–Cattaneo law of heat conduction instead of the conventional Fourier’s law and presented the generalized thermoelastic theory with one relaxation time, which has been proved to be well investigated and well established. For fiber-reinforced generalized thermoelasticity problems, Othman and Said [4] investigated the thermal shock of 2D fiber-reinforced materials and found that the temperature, displacement, and stress components change drastically at the front of the heat wave. Othman and Lotfy [13] compared coupling theory, G-L theory, and L-S theory to prove that the fiber reinforcement and the magnetic field significantly influenced the physical quantities, such as stress and strain. And also, the results from the three theories are in accordance with each other. Abouelregal and Zenkour [14] analyzed the effects of the fractional parameter, reinforcement, and rotation on the variations of different field quantities inside the elastic medium and found that fiber reinforcement plays an important role in the distributions of the field quantities. Abbas [15] investigated the generalized thermoelastic interaction of an infinite fiber-reinforced anisotropic plate containing a circular hole and found that field quantities are significantly varied in the presence and absence of reinforcement.

Thermal conductivity is an important parameter of a material which is typically considered constant. However, several experimental and theoretical studies have indicated that thermal conductivity is closely related to temperature change [1622]. Xiong and Guo [23] validated the effects of variable temperature-dependent properties on field quantities based on a one-dimensional generalized magneto-thermoelastic problem. Wang et al. [24] studied generalized thermoelasticity with variable thermal material properties and found that variable thermal material properties significantly affect the thermoelastic behaviors, particularly the magnitude of thermoelastic response. Ezzat and El-Bary [25] examined the effects of variable thermal conductivity in a problem of a thermo-viscoelastic infinitely long hollow cylinder and discovered that all functions for the generalized theory with a variable thermal conductivity distinctly differ from those obtained for the generalized theory with a constant thermal conductivity. Abo-Dahab and Abbas [26] evaluated the thermal shock problem of generalized magneto-thermoelasticity and concluded that as the variable thermal conductivity increases, the temperature increases, whereas the radial and hoop stresses decrease. These studies demonstrated that variable thermal conductivity significantly influences the material properties and the distribution of field quantities [27]. Instantaneous changes of temperature can markedly change the thermal conductivity of a material. Therefore, the influence of temperature-dependent variable thermal conductivity must be considered in solving fiber-reinforced generalized thermoelasticity problems suffered from thermal stress.

Normal mode analysis is applicable only for solving the steady state problem, whereas integral transformation and time-domain finite element method are suitable for solving dynamic problems. Integral transforms, such as Fourier and Laplace transforms, have been widely used for processing the related generalized thermoelastic problems [2830]. In considering the effect of temperature-dependent variable conductivity in the fiber-reinforced generalized thermoelasticity problem, the governing equation has a nonlinear form. Given that the governing equations contain higher-order terms and nonlinear coupling terms, integral transformation is difficult to perform, and inverse transformation is needed. This process inevitably produces truncation and discrete errors [3133]. The governing equations can be directly solved in the time domain by using the finite element method, thereby avoiding the tedious processes in integral transformation. This method may be more efficient and may have higher precision. Furthermore, the time history of variables in constitutive relations can be reflected. Tian et al. [34] solved 2D generalized thermoelasticity problems by using a direct finite element method, which decreases the solving difficulty in the 2D model due to integral transformation. Li et al. [35] analyzed the nonlinear transient response under the generalized thermal diffusion theory based on the time-domain finite element method and obtained a good effect.

This paper investigates the transient thermal shock problem for fiber-reinforced materials with a time-dependent variable thermal conductivity according to the L-S theory. Time-domain finite element method is applied to derive the nonlinear governing equations. Numerical examples are presented to clarify the transient thermal shock response on a half-space. Field quantities are obtained for different thermal conductivities and illustrated graphically.

2. Governing Equations

Belfield et al. [5] proposed a constitutive equation for a fiber-reinforced linearly thermoelastic anisotropic medium in studying the deformation of fiber-reinforced composites. Given the reinforced direction , , and with consideration of variable thermal conductivity, the constitutive equations can be expressed aswhere is the stress tensor; is the Kronecker delta; is the strain tensor; , , , and are reinforcement parameters, with , in which are Lame constants; is the coefficient of linear thermal expansion; is the reference temperature; , in which is the temperature difference; is the entropy density; is the mass density; is the specific heat at constant strain; and .

The equation of motion (in the context of L-S theory) iswhere is the displacement vector.

The equation of energy conservation iswhere is the heat flux vector.

The geometrical equation is

The fiber-reinforced direction is defined as , and . Thus, equation (1) can be written aswhere

From equations (5) and (7), equation (3) yieldswhere

The equation of heat conduction iswhere is the relaxation time.

The thermal conductivity is temperature dependent and assumed to have the following linear form:where is the initial thermal conductivity and is the small quantity for measuring the influence of temperature on thermal conductivity.

From equations (2) and (4), equation (11) then yields

3. Finite Element Formulations

The finite element method is an approximate method for solving differential equations. The first step in solving the problem is to establish the governing equations, followed by defining the boundary conditions based on the specific problems and then performing the the structural discrete, unit analysis and overall analysis to obtain the numerical solutions. For nonlinear problems, the finite element expression is obtained using the finite element method, which can eliminate the influence of truncation errors and avoid the tediousness of integral transformation. In addition, the time history of the variables in the constitutive relation can be determined to better reflect the wavefront characteristics. FlexPDE is a useful tool for solving partial differential equations, which can form Galerkin finite element integrals, derivatives, and dependencies aiming at the problem description and then build a coupling matrix and solve it. Therefore, FlexPDE is employed to deal with the related partial differential equations generated by the finite element method. For convenience, the constitutive equations of equations (1) and (2) can be expressed in the matrix form as follows:

The heat conduction equation of equation (11) can be written as

The basic variables in this study include displacement and temperature. After the elements are divided, the variables are represented by shape functions in each element as follows:where and are the nodal displacement and temperature, respectively. and are shape functions:where denotes the number of nodes in the grid.

According to equation (5) and given , equation (15) yieldswhere and are the first-order derivative of the components in and with respect to the material coordinates, respectively.

Then, the variational forms of equation (5) are

According to virtual displacement principles, the fiber-reinforced generalized thermoelasticity problem with variable thermal conductivity can be formulated aswhere is the traction vector.

According to equations (13)–(15), and (19), equation (20) yields

These expressions can be summed as the following matrix form:where

4. Numerical Results and Discussion

4.1. Verification

To check the validity of the proposed method, reference [35] is chosen for comparison. Li et al. [35] had investigated the generalized diffusion-thermoelasticity problems with variable thermal conductivity by using the finite element method and had verified its effectiveness. This comparison research is conducted without the consideration of diffusion. In addition to this, the numerical model, initial conditions, and boundary conditions are same with reference [35]. The distribution of the temperature profile at the dimensionless time t = 0.06 is shown graphically in Figure 1, from which a trend consistency can be observed. This guarantees the validity and accuracy of the present method.

4.2. Results and Discussion

Consider the problem of a fiber-reinforced anisotropic elastic half-space () with variable thermal conductivity. As shown in Figure 2(a), the boundary surface is assumed to be without traction, and the banded area on is subjected to a time-dependent transient thermal shock.Initial conditions:Boundary conditions:where is the Heaviside unit step function.

Under the given conditions, the half-space model can be simplified as a xy plane model, as shown in Figure 2(b). The components of displacement and temperature can be simplified as follows:

Copper material is selected for numerical evaluation, and the parameters are presented in Table 1 [36].

For convenience, the following nondimensional quantities are introduced:

According to equation (28), equations (7), (9), and (12) can be written as follows (the asterisk is removed for brevity):where

Under the given conditions of equations (24)–(27), the nonlinear governing equations (29)–(31) can be solved directly in the time domain. Under the given conditions, the research model is symmetrical about the x-axis. Thus, the rectangular area is used for the subsequent analysis for simplification.

The dimensionless distributions of temperature, displacement, and stress are illustrated graphically in Figures 38. The K1 value is expressed in four cases to discuss the effects of variable thermal conductivity as follows:(1)Case 1: K1 = −0.6(2)Case 2: K1 = −0.3(3)Case 3: K1 = 0.

Figures 3 and 4 show the distributions of nondimensional temperature along and . The dimensionless temperature is equal to 1 at and region, which agrees with the boundary conditions that were previously assigned. As the thermal conductivity increases, the temperatures along and increase. This result indicates that the coefficient of thermal conductivity is positively correlated with the temperature change, and that variable thermal conductivity significantly affects the temperature distribution. Moreover, the wavefront effect is more pronounced as the thermal conductivity increases. There is a diverse trend at in Figure 3, which is mainly caused by the effect of wavefront and the calculation errors.

Given that the research model is symmetrical about the x-axis, the vertical displacement along is zero and need not be considered. Figure 5 depicts the distribution of the nondimensional horizontal displacement , which shows that the closest section to the origin undergoes expansion, the next undergoes compression, and the rest away from the origin is undisturbed. Figure 5 demonstrates that variable thermal conductivity is positively correlated with the distribution of the horizontal displacement along .

Figures 6 and 7 show the distribution of the nondimensional horizontal displacement and the vertical displacement , respectively. As shown, the variable thermal conductivity obviously affected the distribution of displacement. In addition, is the most affected thermal shock boundary point. The negative value in Figure 6 indicates that the particles tend to move toward the unconstrained direction.

Figure 8 shows the distribution of the nondimensional stress . The stress values show a violent oscillation in the zone, and variable thermal conductivity obviously affects the distribution of stress along , particularly in the zone.

5. Concluding Remarks

It is well known that instantaneous changes of temperature can markedly change the thermal conductivity of a material. Therefore, this article investigated the effect of temperature-dependent variable thermal conductivity on a fiber-reinforced generalized thermoelastic half-space. Given the reinforced direction , a region of its surface is subjected to a transient thermal shock. The problem is studied in the context of L-S theory. The time-domain finite element method is proposed to analyze the nonlinear response.

Based on the simulation, we can draw that the time-finite element method is very effective for analyzing nonlinear problems with given initial and boundary conditions, and by which we can capture a pronounced wavefront effect. In consideration of the fiber-reinforced effect, variable thermal conductivity positively affects the distributions of temperature, displacement, and stress. In addition, the boundary point of thermal shock is affected the most.


:Components of stress tensor
:Kronecker delta
:Components of strain tensor
, , and :Reinforcement parameters
, :Lame constants
:Coefficient of linear thermal expansion
:Reference temperature
:Temperature difference
:Entropy density
:Mass density
:Specific heat at constant strain
:Displacement vector
:Heat flux vector
:Relaxation time
:Initial thermal conductivity
:Small quantity for measuring the influence of temperature on thermal conductivity
:Number of nodes in the grid
:Heaviside unit step function
and :Shape functions
and :First-order derivative of and with respect to the material coordinates
:Traction vector.

Data Availability

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

Conflicts of Interest

There are no conflicts of interest regarding the publication of this paper.