Advances in Materials Science and Engineering

Volume 2019, Article ID 8625371, 9 pages

https://doi.org/10.1155/2019/8625371

## Effect of Variable Thermal Conductivity on the Generalized Thermoelasticity Problems in a Fiber-Reinforced Anisotropic Half-Space

School of Civil Engineering, Tianjin University, Tianjin 300072, China

Correspondence should be addressed to Li-Na Yu; nc.ude.ujt@aniluy

Received 23 April 2019; Accepted 28 July 2019; Published 3 September 2019

Academic Editor: Paweł Kłosowski

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

#### Abstract

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 [2–6]. 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 [9–12], 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 [16–22]. 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 [28–30]. 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 [31–33]. 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.