Abstract
This paper studies the problem of a cracked orthotropic solid subject to linear thermal flux and linear mechanical load. The proposed extended partially insulated crack model is employed to simulate two collinear cracks. Taking advantage of Fourier transform technique and superposition theory, the closed form of some physical quantities and fracture parameters is obtained. Some simple examples are employed to demonstrate dimensionless thermal conductivity () between the upper and below crack regions, and the proposed coefficient () has great effects on some physical quantities and fracture parameters.
1. Introduction
Multicomponent composite materials are widely used in the material industry. However, considering the complex factors involving working environment, internal and external loads, and production process, it is inevitable to contain a series of various cracks in these solids. The appearances of different kinds of cracks will reduce the capacity of cracked structures and even bring about severe accidents. Therefore, it is vital to do some research on fracture analysis of a cracked solid by utilizing the theory of thermal elasticity for the purpose of safety [1–3]. With the rapid growth of thermoelasticity theory, a great deal of treatises and papers was published to investigate fracture characters of cracked solids [4–6]. The fracture parameters of an orthotropic material containing a central crack under heat flow were obtained by Tsai [7]. The closed form of fracture parameters of cracked orthotropic solids was calculated by Ju and Rowlands [8]. The closed form of some physical quantities of two collinear cracks was studied by Chen and Zhang [9]. The transient thermal problem of a cracked orthotropic plate was taken into account by Noda [10]. Some physical quantities of a cracked orthotropic semi-infinite medium were given by Rizk [11]. On the other hand, the thermoelastic problems of orthotropic functionally graded solids brought about the widespread attention. For example, the fracture parameters of orthotropic functionally graded solids under mechanical load were given explicitly by Kim and Paulino [12]. The problem of a cracked solid subject to plane temperature-step waves was investigated by Brock [13]. The equivalent domain integral was formulated to study the fracture problems subject to thermal stresses by Dag [14]. The problems of cracked orthotropic solids subject to symmetrical thermomechanical loads with application of Fourier transform technique (FTT) and superposition principle were studied by Wu et al. [15].
Subject to thermal load, the analysis of fracture behavior for cracked solids which were often regarded as orthotropic or isotropic has generated enormous publicity [16–20]. To simulate two collinear cracks, a partially insulated crack model prevailed [21–23]. where the definitions of , , and have been given in detail [3]. The case of or denotes a fully thermally impermeable or permeable state.
The following extended partially insulated crack model is also put forward by virtue of mathematical intuition. where presents initial heat flux. The coefficient is considered a constant. Whether it is negative or positive is mainly relies on the portion of thermal flux and mechanical load. Clearly, the crack model proposed in (2) returns to (1) when .
The reasons of introducing constant in (2) are as follows. First, the value of does not precisely address the cracks with thermal resistance. Second, the constant , which is introduced as an adjustment factor, conforms to the complex situation and meets the abnormal state of crack surface.
This paper employs an extended partially insulated crack model to discuss two collinear cracks under linear thermal flux and linear mechanical load. The thermoelastic field is given in explicit form based on the proposed extended partially insulated crack model, Fourier transform, and superposition theory. The results show the effects of dimensionless thermal conductivity () between the upper and below crack regions and the proposed coefficient () on and and . It is revealed the boundary conditions of crack surface, thermal properties of crack, and the raised coefficient should be paid attention to the analysis of crack growth under thermal load in numerical results.
2. Problem Statement
Two collinear cracks in an orthotropic solid are taken into account as shown in Figure 1. They are located at .
Making use of the state of plane stress [3], we obtain where where the definitions of , , , , , , , , , , and have been given in [3].
One obtains
Making use of the Fourier heat conduction leads to where the definitions of , , , and have been given in [3]. Furthermore, based on the equilibrium equation, one has
Taking advantage of the thermal equilibrium equations brings out where
The crack face boundary conditions are depicted as
Hereafter, the subscript ‘’ or ‘’ denotes the physical quantity of the upper () or below () part. , , , and stand for the prescribed constants. Based on (14) and (15), thermal flux is composed of only one part () and mechanical loading is divided into two parts ( and ). As linear thermal flux and linear mechanical load are antisymmetrical and symmetrical, respectively, the thermoelastic field of the region () is only dealt with. The crack-surface boundary conditions are expressed with the application of the improved partially insulated crack model.
According to Equations (17) and (18), the solutions under thermal flux () and mechanical loading () have been given explicitly in [24, 25]. Next, we depict the boundary conditions of crack-surface subject to linear mechanical load (). where
Besides, some physical quantities conform to the following conditions:
3. Solution Procedure
3.1. Temperature Field
According to [25], one obtains the explicit form of temperature difference on crack faces as
3.2. Elastic Field
To achieve the goal of explicit form in Equations (8) and (9), and are expressed according to [26]. where the definitions of and have been given in [26].
Hereafter, or denotes or . need to solve. The definitions of have been given in [26]. where
Furthermore, and are chosen as
Taking advantage of Equations (29), (30), (8), and (9), we have where
By the aid of Equations (3)–(5), (25), (26), (29), and (30), the components of stress are in the form of the following expressions:
In order to get the explicit solution of this considered problem, we depict the dual integral equations as
Using Equations (36) and (37), one gets
Applying Equations (20) and (37), one obtains where
In order to solve Equations (39) and (40), the auxiliary function is introduced as
Applying inverse Fourier transform leads to
Inserting Equation (43) into (40), one has
Recalling the known result [27],
Based on Equation (45), Equation (44) can be expressed as
It is convenient to introduce , and . Equation (46) is rewritten as
According to the singular integral containing the Cauchy kernel4 [28], the solution of Equation (47) is obtained
In the application of Equation (42), the closed form of elastic displacement is obtained
Inserting Equation (48) into (34), the stress field is obtained as
and denote the first and second kinds of complete elliptical integrals, respectively, where
For simplicity, the detailed procedure of reduction under thermal flux is omitted. The shear stresses are obtained according to [25]. where
By superposition theory, the exact solutions of the physical quantities are obtained subject to linear thermal flux () and linear mechanical load ().
3.3. Crack-Tip Field
Using Equations (2) and (23), one obtains the closed form of heat flux to the crack surface.
We define the value of to stand for the dimensionless thermal resistance between crack faces. It is easily found Equation (54) is different from that in [25]. When and , one obtains or , meaning partially thermally insulated or fully conductive cracks. When and , one has , meaning fully thermally insulated cracks.
4. Fracture Parameters
It is important that the stress intensity factors including the mode-I and mode-II should be defined as the analysis of cracked growth.
Based on Equation (55), one can obtain
When , it means the mode-I stress intensity factor of a single crack with . It is easily found that
According to Equation (56) and Reference [25], one has
The importance of strain energy in a unit volume of the solid is illustrated for nonisothermal [29, 30]. where the definitions of and have been given in [24]. For the orthotropic solid, Equation (60) can also be given based on the above concepts of energy density function
The following strain energy density factor on the crack line is defined to study crack growth in fracture mechanics [24]. where
5. Numerical Results
For the sake of simplicity, some numerical examples are employed to demonstrate and have great effects on subject to linear thermal flux() and linear mechanical load (). The orthotropic material like Tyrannohex is selected as in [31] (Table 1).
Figure 2 shows versus with for and . increases with an increase of for a fixed as shown in Figure 2. The case of implies the heat conduction between the upper and blow crack faces is the same as that of external material of crack. The cases of imply to the heat conduction between the upper and blow crack faces are twice, triple, and quadruple as much as thermal conductivities of the external material of crack. Figure 3 displays versus with for and . As the dimensionless thermal resistance increases, increases. The constant is considered to be an adjustment quantity. The bigger the value of constant , the greater the heat flux per thickness through crack. It means making use of the extended partially insulated crack model involving the greater or will overestimate the heat flux per thickness to the crack surface. Furthermore, it is suitable to decease or increase stress field near outer and inn cracks tip by the way of filling certain materials into the region between the upper and below crack faces according to Figures 2 and 3.
Figure 4 displays or versus with where and denote and for , respectively. or decreases when increases for a fixed . Figure 5 displays or versus with for and . As the dimensionless thermal resistance increases, the mode-II stress intensity factors decrease. The bigger the value of constant , the smaller the mode-II stress intensity factors. It means making use of the extended partially insulated crack model involving the greater and will underestimate the mode-II stress intensity factors. The obtained results reveal that the crack face boundary conditions, the thermal properties of crack, and the raised coefficients have great influences on the heat flux per thickness to the crack surface and the mode-II stress intensity factors.
In order to present the influence of the thermal properties of crack on , the value of is easily defined. which denotes the strain energy density factor of a crack with under mechanical load . Figures 6 and 7 show and versus with for , MPa, and . Figures 6 and 7 respond to the two cases of strain energy density factor near outer and inn cracks for partially thermally insulated cracks. It is easily seen that and are made up of the mode-II stress intensity factor under thermal flux () and the mode-I stress intensity factor induced by mechanical load (). The corresponding and increase with an increase of thermal flux and mechanical load. The strain energy density factor on the crack line is greatly influenced by the adjustment quantity . The bigger the value of constant , the smaller and . So, applying the bigger value of constant will underestimate and .
From the above figures, it is revealed and have significant impacts on the analysis of a cracked solid. In other words, some physical quantities (i.e., and ) should be given enough attention to the analysis of the thermoelastic field.
6. Conclusions
This paper addresses two collinear cracks in an orthotropic solid under linear thermal flux and linear mechanical load in this paper. Some physical quantities and fracture parameters are obtained in explicit forms with application of the proposed extended partially insulated crack model, Fourier transform, and superposition theory. The results show that and have vital effects on Q1c and some fracture parameters. The obtained results reveal the boundary conditions of crack face, thermal properties of crack, and the raised coefficients should be concerned about the analysis of a cracked solid under the thermal load.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The work was supported by the Hebei University Scientific Research Foundation for higher-level talent (No. 521100221019), the Project Funded by China Postdoctoral Science Foundation (No. 2019M662946), and the Ningxia Key Research and Development Program (Introduction of Talents Project) (No. 2020BEB04039).