Next Article in Journal
Compensation Method for Correcting the Topography Convolution of the 3D AFM Profile Image of a Diffraction Grating
Previous Article in Journal
Electromyography-Based Biomechanical Cybernetic Control of a Robotic Fish Avatar
Previous Article in Special Issue
Unveiling Inertia Constants by Exploring Mass Distribution in Wind Turbine Blades and Review of the Drive Train Parameters
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Digital Twins to Predict Crack Propagation of Sustainable Engineering Materials under Different Loads

1
School of Intelligent Manufacturing, Chengdu Technological University, Chengdu 610031, China
2
Department of Civil and Mechanical Engineering, Purdue University Fort Wayne, Fort Wayne, IN 46805, USA
*
Author to whom correspondence should be addressed.
Machines 2024, 12(2), 125; https://doi.org/10.3390/machines12020125
Submission received: 9 January 2024 / Revised: 31 January 2024 / Accepted: 8 February 2024 / Published: 10 February 2024
(This article belongs to the Special Issue Estimation and Mitigation of Fatigue Damage for Wind Turbines)

Abstract

:
Computer-aided engineering (CAE) is an essential tool in a digital twin not only to verify and validate a virtual twin before it is transformed into a physical twin, but also to monitor the use of the physical twin for enhanced sustainability. This paper aims to develop a CAE model for a digital twin to predict the fatigue life of materials. Fatigue damage is represented by the size of a macro-crack that grows with a cluster of micro-cracks subjected to three different loads. The growth angle is related to the maximum circumferential tensile stress, and the growth rate is determined by the stress intensity factor (SIF) at the crack tip. The prediction model takes into consideration the main factors, including micro-cracks, crack closures, and initial configurations. Simulations are developed for the growth of macro-cracks with radially distributed micro-cracks and randomly distributed micro-cracks, and we find that (1) the macro-crack in the second case grows faster than that in the first case; (2) a pure shear load affects the macro-crack propagation more than a combined shear and tensile load or a tensional load; (3) the external stresses required to propagate are reduced when the inclination angle of the micro-crack is small and within (−25° < β < 25°); (4) micro-cracks affect the propagating path of the macro-crack and generally guide the direction of propagation. The developed model has been verified and validated experimentally for its effectiveness in predicting the fracture or fatigue damage of a structure.

1. Introduction

Pursuing sustainability in design and manufacturing implies that the capability of an engineering material in a product, structure, or system must be fully utilized to reduce or avoid any waste. This cannot be achieved without developing a digital twin for a physical twin so that the physical system (1) is built when its virtual model has been fully verified and validated and (2) is constantly monitored to respond to any conditional changes for maintenance or remedies as soon as possible [1,2,3]. Since physical products or systems are mainly made from engineering materials, we are especially interested in developing a generic CAE tool as a part of a digital twin to predict the fatigue damage of engineering materials. This can be an integral model of a digital twin where fatigue failure can be a major source of product failure, especially for aerospace products [4].
This paper focuses on fatigue prediction, which is associated with the accumulated damage of crack, and can be measure based on crack size. It is rare that fatigue damage consists of a single crack; in fact, numerous micro-cracks have been found near the tip and the path of a macro-crack using scanning electron microscopy in 7075-T6 aluminum alloy, X65 steel, and U71Mn steel [5,6,7]. Depending on their relative sizes, cracks are referred to as macro-cracks and micro-cracks. We are especially interested in the effects of micro-cracks on macro-cracks that are subject to various loading patterns.
Micro-cracks can enhance or reduce the fatigue strength of materials depending on a number of factors. To investigate the effects of fatigue strength on micro-cracks, the complex potential method and pseudo-traction method can be used to model the interactions between macro- and micro- cracks [8,9,10]. Kachanov [11,12] developed an interaction matrix and represented the effects of micro-cracks on a macro-crack. Tamuzs and Petrova [13,14] used the small parameter method to study the effect of a cluster of cracks on a macro-crack. Li et al. [15,16,17] assumed that micro-cracks were arbitrarily distributed and adopted the distributed dislocation technique to estimate their impacts on a macro-crack. The extended finite element method (XFEM) was used by Loehnert and Belytschko [18] and Sutula et al. [19] to analyze the effects of micro-cracks.
While the afore-mentioned works discussed the interactions between a macro-crack and micro-cracks, none of them considered the relationship and propagation of micro-cracks. We argue linear elastic fracture mechanics (LEFM) should be adopted to investigate how micro-cracks are propagated.
The direction and growth rate are two critical factors in propagating cracks [20,21,22,23]. Kahn [24] analyzed the impact of the initial crack angle in the mixed crack mode subject to different loads and concluded that the extension of cracks depended on the direction of stress with respect to the crack face. Wang [25] developed analytical and numerical solutions to the stress distribution in the surrounding area of a crack and found that the stress distribution was affected by the loadings and geometrical variations. Finite element modelling (FEM) was used by Ghaffari [26] to study the initiation and propagation of cracks in gears. These works ignored the conditions in which micro-cracks pre-existed. In general, a fracture in the real world is mostly caused by the propagation of a single macro-crack and is affected by pre-existing micro-cracks [7,27,28,29,30].
Many researchers have contributed to the study of initiation, propagation, and coalescence of cracks. It has been generally agreed upon that crack propagation is closely associated with fatigue fractures. For example, Li et al. [31] applied the distributed dislocation technique to variations in plastic zones (PZ) and the locations of micro-cracks. They found that the micro-cracks in front of a propagated macro-crack did not affect the growth of the macro-crack. Basoglu [32], Yates [33], and Liu [34] investigated the changes in a micro-crack when multiple micro-cracks were presented. The experimental methods were also used to investigate the initiation of and changes in macro-cracks [35,36,37,38]. There was a need to model the effects of micro-cracks on the fatigue growth of macro-cracks.
The rest of the paper is organized as follows. Section 2 presents a model of an infinite plane that consists of a macro-crack and multiple micro-cracks. The Muskhelishvili’s complex potential method is used to find the solution, subject to various loading conditions. Section 3 describes the propagations of cracks, including a macro-crack and multiple micro-cracks. Section 4 verifies the results of the proposed solutions. Section 5 summarizes the effects of micro-cracks on a macro-crack and provides a summary with conclusions.

2. Crack Propagation Modelling and Solution

In this section, a generalized representation of an engineering material with accumulated fatigue damage is presented, and the Muskhelishvili’s complex potential method is used to find the solution of crack propagation subject to various loading conditions.

2.1. Problem Description

To maintain generality, this paper models an engineering material as an elastic plane with one macro-crack and a set of arbitrarily oriented and distributed micro-cracks subjected to a uniaxial tensional stress of σ and a shear stress of τ, as shown in Figure 1.
To align the propagation of the macro-crack to the set of micro-cracks shown in Figure 1, we decompose the original problem in Figure 1 into three sub-problems, as shown in Figure 2, and look into the corresponding solutions.
Figure 2 shows that the original problem is decomposed into three sub-problems. The macro-crack is sized based on its length (2a). The micro-crack k is sized based on its length (2ak), where k is the index of the micro-crack. The global coordinate system (o-x-y) is established with its origin (o) at the center of the macro-crack, and x and y are aligned with the length and perpendicular direction of the macro crack, respectively. A local coordinate system (o-xk-yk) is established with its origin (ok) at the center of micro-crack k, and xk and yk are aligned with the length and perpendicular direction of micro-crack k, respectively. Moreover, βk and θk denote the inclined orientational angle of micro-crack k, respectively, and dk is the distance from o to ok.
The Muskhelishvili’s complex potential method is adopted to solve the formulated problems, and the solutions must satisfy all the boundary conditions applied to the cracks. For sub-problem (i), the stress distribution surrounding micro-crack j is determined based on the stress function available in ref. [30]. For sub-problem (ii), the resulting stresses apply to micro-crack j in the opposite direction. Due to the coupling of multiple micro-cracks, the stress on micro-crack j exists, but with a small magnitude. For sub-problem (iii), the solution includes both the redundant stress of micro-crack j from sub-problem (ii) and that at the position of non-existent micro-crack j. For sub-problem (iv), the stress of non-existent micro-crack j obtained from the solution to sub-problem (iii) is applied to real micro-crack j.

2.2. Solving Procedure

Assume the material is isotropic and elastic and its deformation on the plane of interest is within the elastic range of the material, the stress distribution around a crack is expressed by two complex potentials, ϕ(z) and ψ(z) [39], as follows:
σ x x + σ y y = 2 ϕ ( z ) + ϕ ¯ ( z ) σ x x σ y y + 2 i τ x y = 2 z ¯ ϕ ( z ) + ψ ( z )
where z = x + i y ; σ x x and σ y y are normal components along the horizontal and vertical directions, respectively; τ x y is the shear stress; and (′) and (″) are the operations for the first and second derivatives with respective to z, respectively.

2.2.1. Solution to Sub-Problem (i) for the Macro-Crack

The stress is distributed as follows [39]:
ϕ z = i τ 2 z z 2 a 2 ,   ψ z = i τ a 2 2 z a 2 z 2 a 2
substituting Equation (2) into Equation (1) to find the stress around the tip of a macro-crack. For micro-crack k, the stress components σ x x ( k ) ( i ) , σ y y ( k ) ( i ) , and τ x y ( k ) ( i ) can be obtained as follows:
σ x x ( k ) ( i ) = σ x x σ + σ x x τ ,   σ y y ( k ) ( i ) = σ y y σ + σ y y τ ,   τ x y ( k ) ( i ) = τ x y σ + τ x y τ
Further conducing coordinate transformation, the normal and tangential components of the stress around micro-crack k in the local coordinate system are as follows:
σ N k M I i = σ x x k i + σ y y k i 2 + σ x x k i σ y y k i 2 c o s 2 β k τ x y k i s i n 2 β k τ T k M I i = σ x x k i σ y y k i 2 s i n 2 β k + τ x y k i c o s 2 β k   k = 1 , 2 , n
where subscripts N and T denote the normal and tangential direction of the crack, respectively.

2.2.2. Solution to Sub-Problem (ii) for Micro-Cracks

We find the sub-solution one micro-crack at a time; if there are k micro-cracks, the solving procedure will be repeated k times. In each step, it is assumed that an infinite plane contains a single micro-crack j, one virtual macro-crack, and k − 1 virtual micro-cracks. It is further assumed that micro-crack j is subjected to a non-uniform distributed load of σ N ( j ) M A ( i ) and τ T ( j ) M A ( i ) from the solution in Section 2.2.1, but with the opposite direction.
Firstly, the load is treated as an assembly of numerous small, concentrated forces F, and the stress due to F applied to the upper and lower boundary of the crack becomes [39]:
ϕ ζ = F 2 π i ζ 0 ζ ( ζ 0 2 1 ) κ κ + 1 l n ζ ln ζ 0 ζ + ζ 0 2 ( a ζ + 1 ζ 2 z 1 ) a ζ 0 ζ ( ζ 0 2 1 ) i ( ζ 0 2 1 ) ζ 0 ζ 0 ζ ¯ 0 ψ ζ = F 2 π i 2 ζ ζ 0 ζ 0 2 ( ζ 2 + 1 ) ( ζ 2 1 ) ( ζ 0 2 1 ) + l n ζ κ + 1 + 2 κ + 1 ( ζ 2 1 ) ln ζ 0 ζ 2 z ¯ 1 ζ 0 2 a ζ 0 ζ ( ζ 0 2 1 ) i ( ζ 0 2 1 ) ζ 0 ζ 0 ζ ¯ 0  
The concentrated forces F in the local coordinate system O j x j y j correspond to the distributed load [39]:
F = P + i Q = σ N ( j ) M I ( i ) d x F + i   τ T ( j ) M I ( i ) d x F  
where P and Q are the components in the normal and tangential directions, respectively. Equation (6) can be rewritten using the conformal mapping method to transform ϕ(z) and ψ(z).
The stress components σ x x F M A ( i i ) , σ y y F M A ( i i ) , τ x y F M A ( i i ) , σ x x F M I ( i i ) , σ y y F M I ( i i ) , and τ x y F M I ( i i ) acting on the virtual macro-crack and micro-cracks due to F can be obtained using Equations (1), (4), and (5). The coordination transformation can be applied to find these stress components in the local coordinate system O j x j y j as follows:
σ x x j M A i i , σ y y j M A i i , τ x y j M A i i = a k a k σ x x F M A i i , σ y y F M A i i , τ x y F M A i i d x F σ x x k M I i i , σ y y k M I i i , τ x y k M I i i = a k a k σ x x F M I i i , σ y y F M I i i , τ x y F M I i i d x F σ x x k M I i i , σ y y k M I i i , τ x y k M I i i = 0 k = j   k = 1 , 2 , n , k j
For the macro-crack, subscript (j) shows the contribution of micro-crack j to the macro-crack. For a micro-crack, subscript (k) shows the contribution of micro-crack j to the k − 1 virtual micro-crack.
Next, transforming the results from O j x j y j to the global coordinate system o-x-y and local coordinate systems for a micro-crack as O k x k y k ( k , 1 , 2 , , n , ) , σ N ( j ) M I ( i i ) , τ T ( j ) M I ( i i ) , σ N ( k ) M I ( i i ) , and σ T ( k ) M I ( i i ) can be obtained when the stress is applied to micro-crack j. The same procedure applies to any of other k − 1 micro-cracks as follows:
σ N M A ( i i ) = j = 1 n σ N ( j ) M A ( i i ) τ T M A ( i i ) = j = 1 n σ T ( j ) M A ( i i )  
σ N M I ( i i ) = k = 1 n σ N ( k ) M I ( i i ) τ T M I ( i i ) = k = 1 n σ T ( k ) M I ( i i )   ( k = 1 , 2 , n )  
Taking into consideration of effects of other micro-cracks, the stress in Equation (8) is not zero. However, the stress magnitude in each micro-crack is insignificant in comparison with σ N ( j ) M I ( i ) and τ T ( j ) M I ( i ) , which were solution to sub-problem (i).

2.2.3. Solution to Sub-Problem (iii) for the Macro-Crack

Here, the effects of one macro-crack in the plane are investigated. The macro-crack is subjected to a non-uniformly distributed stress whose magnitude and direction are determined in Equation (7). The sample procedure in Section 2.2.2 is applied. When a concentrated force F is applied to the macro-crack, σ x x F M I ( i i i ) , σ y y F M I ( i i i ) , and τ x y F M I ( i i i ) on a virtual micro-crack are obtained as follows [39]:
F = P + i Q = σ N M A ( i i ) d x F + i   τ T M A ( i i ) d x F  
The stress components in micro-crack k are as follows:
σ x x ( k ) M I ( i i i ) , σ y y ( k ) M I ( i i i ) , τ x y ( k ) M I ( i i i ) = a k a k σ x x F M I ( i i i ) , σ y y F M I ( i i i ) , τ x y F M I ( i i i ) d x F ( k = 1 , 2 , n )
The normal and shear stresses σ N ( k ) M I ( i i i ) and τ T ( k ) M I ( i i i ) in O k x k y k   ( k , 1 , 2 , , n ,   ) can be obtained through coordinate transformation, and the stresses of the virtual micro-cracks are obtained in this step.

2.2.4. Solution to Sub-Problem (iv) from Micro-Cracks

The stresses of micro-cracks have been found based on the solutions to sub-problems (ii) and (iii); therefore, the loads by the micro-cracks become σ N ( k ) M I ( i i ) + σ N ( k ) M I ( i i i ) and τ T ( k ) M I ( i i ) + τ T ( k ) M I ( i i i ) , and, similar to Equations (7) and (8), the stress components of σ N M A ( i v ) ,   τ T M A ( i v ) ,   σ N ( k ) M I ( i v ) and τ T ( k ) M I ( i v )   ( k = 1 , 2 , , n ) can be obtained.
The aforementioned procedure shows that solving sub-problem (i) leads to the stress components of σ N ( k ) M I ( i ) and τ T ( k ) M I ( i ) . Solving the sub-problem leads to the stress components of σ N ( k ) M I ( i ) and τ T ( k ) M I ( i ) on micro-cracks; moreover, σ N ( k ) M I ( i i ) and τ T ( k ) M I ( i i ) in Equation (8) are negligeable in comparison with σ N ( j ) M I ( i ) and τ T ( j ) M I ( i ) from the solution to sub-problem (i). Solving sub-problem (iii) leads to the stress components of σ N ( k ) M I ( i i i ) and τ T ( k ) M I ( i i i ) .
To define the constraints in sub-problem (iv), σ N ( k ) M I ( i i ) ,   τ T ( k ) M I ( i i ) , σ N ( k ) M I ( i i i ) , and τ T ( k ) M I ( i i i ) are summed as the stresses due to micro-cracks. Note that the stress components σ N ( k ) M I ( i v ) and τ T ( k ) M I ( i v ) are negligible in comparison with σ N ( k ) M I ( i i i ) and τ T ( k ) M I ( i i i ) . The sub-problems can be further decomposed to improve the accuracy of the developed models.
In the proposed model, the coupling among the cracks is considered for all sub-problems. Sub-problem (i) has modelled the effects of the macro-crack on the micro-cracks; sub-problem (ii) has modelled the effects of the micro-cracks on the macro-crack; sub-problem (iii) has re-refined the effects of the macro-crack on the micro-cracks, and this procedure can be iteratively performed until the accuracy becomes satisfactory.

2.3. Solution to KI and KII

Traction-free boundary conditions on the crack surface can be satisfied using a step-by-step procedure of problem decomposition, as discussed in Section 2.2, since the solution to each sub-problem can be obtained using LEFM independently.

2.3.1. Solution of LEFM

Let the loads on the macro-crack surface be σ N M A ( i ) = σ and τ N M A ( i ) = τ in sub-problem (i), and the SIFs at the tip of crack are as follows [39]:
K I M A ( i ) = a a σ π a a + x a x 1 / 2 d x = σ π a K I I M A ( i ) = a a τ π a a + x a x 1 / 2 d x = τ π a  
Using the solutions to sub-problems (i), (iii), and (v) in the previous section, the SIFs at the tip can be further determined as follows:
K I , e l a s t i c M A = σ π a + a a σ N M A ( i i ) + σ N M A ( i v ) + σ N M A ( v i ) . π a a + x a x 1 / 2 d x K I I , e l a s t i c M A = τ π a + a a τ T M A ( i i ) + τ T M A ( i v ) + τ T M A ( v i ) . π a a + x a x 1 / 2 d x  
The stress components on micro-crack k include (1) σ N ( k ) M I ( i ) and τ T ( k ) M I ( i ) from the solution to sub-problem (ii), (2) σ N ( k ) M I ( i i ) + σ N ( k ) M I ( i i i ) and τ T ( k ) M I ( i i ) + τ T ( k ) M I ( i i i ) from the solution to sub-problem (iv), and (3) σ N ( k ) M I ( i v ) + σ N ( k ) M I ( v ) and τ T ( k ) M I ( i v ) + τ T ( k ) M I ( v ) from the solution to sub-problem (vi). Therefore, the SIFs at the tip of micro-crack k become:
K I ( k ) , e l a s t i c M I = a k a k σ N ( k ) M I ( i ) + σ N ( k ) M I ( i i ) + σ N ( k ) M I ( i i i ) + σ N ( k ) M I ( i v ) . π a k a k + x a k x 1 / 2 d x K I I ( k ) , e l a s t i c M I = a k a k τ T ( k ) M I ( i ) + τ T ( k ) M I ( i i ) + τ T ( k ) M I ( i i i ) + τ T ( k ) M I ( i v ) . π a k a k + x a k x 1 / 2 d x ( k = 1 , 2 , n )  

2.3.2. Correction of SIF for Ductile Materials

Crack propagation will be affected when a yield occurs to a ductile material. Here, the effect of plastic strain is modelled as a correction to SIF by a plastic zone (PZ). It is assumed that the plastic strain is relatively small (a >> r) in comparison with the size of the visible crack [4].
Firstly, for a crack with a mixed mode, the stress components in three directions are determined for failure mode I and mode II as follows [39]:
σ x = K I 2 π r 1 , σ y = K I 2 π r 1 , τ x y = K I I 2 π r 1 , σ z = τ x z = τ z y = 0
where KI and KII are evaluated in Equations (12) and (13) for macro- and micro cracks, respectively. Accordingly, the von Mises stress can be found as follows [39]:
σ e q = σ x σ y 2 + σ x σ z 2 + σ z σ y 2 + 6 ( τ x y 2 + τ y z 2 + τ x z 2 ) 2 1 / 2
where σeq is the equivalent von Mises stress. The characteristic size r1 of the plastic zone meets the condition of ref. [39]:
σ e q = K I 2 + 3 K I I 2 π r 1 = σ s
When the material is ideally elastoplastic, r1 = ry is small and is added as the extended length of a crack (see Figure 3) as follows [39]:
σ e f f = a + r y
Figure 3 shows that ry is doubled due to the effects of stress relaxation. The distance from the crack tip to the tip of an imaginary crack is as follows [39]:
r y = K I 2 + 3 K I I 2 π σ s 2
The resultant stress at the position of the extended length ry is equal to the yielding strength and is distributed uniformly in Figure 3. Therefore, the normal and tangential stresses of the crack become:
σ n , p l a s = K I 2 K I 2 + 3 K I I 2 1 / 2 σ s τ n , p l a s = 3 K I I 2 K I 2 + 3 K I I 2 1 / 2 σ s
With the consideration of PZ, the SIFs are corrected as follows:
K I , p l a s t i c = ( a + r 1 ) a σ n , p l a s π a c a c + x a c x 1 2 d x + a a + r 1 + σ n , p l a s π a c a c + x a c x 1 2 d x K I I , p l a s t i c = ( a + r 1 ) a τ n , p l a s π a c a c + x a c x 1 2 d x + a a + r 1 + τ n , p l a s π a c a c + x a c x 1 2 d x  
where a c = a + r 1 + r 1 + is the modified length of the macro-crack.

3. Discussion on Crack Propagation

Fatigue damage is processed at three phases, i.e., initiation, steady propagation (SP), and rapid propagation (RP) [4]. When the size of a crack is less than a threshold, it experiences SP; otherwise, it experiences RP, which rapidly leads to a fracture. It is worth looking into how a crack is propagated in SP and RP.

3.1. Steady Propagation (SP)

Based on their size, cracks can generally can be divided into three categories: microstructural short cracks, physical short cracks, and macroscopic cracks. A short crack becomes a macroscopic crack and SP is transformed into RP when the size of a crack reaches a critical length. Crack propagation is affected by many factors such as mean and fluctuated stresses, temperature, and microstructures. Pook [40] and Milella [41] commented that a critical length for metals is usually in a range of (0.25, 0.3) mm. Miller [42] and Stephens et al. [43] indicated that a critical size related to grain size is typically a depth of three to five times the grain size.
Here, we assume that the length of a macro-crack is fixed and above the critical length of a crack while the sizes of micro-cracks are randomly selected and below the critical length of a crack. The threshold SIF of the micro-crack is found using the El Haddad model [44] as follows:
Δ K t h = a k a k + a 0
where ak is the micro-crack length, Δ K t h is the threshold SIF of the micro-crack, and a0 is the critical length of the material, determined by [44]:
a 0 = 1 π Δ K t h Δ σ o 2  
where Δσ0 is the fully reversed fatigue strength. Taking an example, Table 1 shows the material properties of U71Mn steel.

3.2. Rapid Propagation (RP)

The level of a crack closure (U) is characterized by:
U = 1 K o p / K I , m a x 1 R  
where R is the stress ratio, Kop is the SIF corresponding to an opened crack, and KI,max is the intensity factor of the mode-I crack.
It is known that the residual PZ at the tip of a crack in the last loading cycle causes a crack closure. Wolf [46] provided an empirical equation to calculate the U of an aluminum alloy as U = 0.55 + 0.33R + 0.12R2. Note that Equation (24) takes into consideration the plastic-hardening effect and Kop = 0.45KI,max when U = 0.55 and R = 0. Therefore, the crack was closure free, and the closure level was U = 1.
For RP, the closure level was affected significantly by the stress ratio and propagation modes. Moreover, we took into consideration the impact of mode I and mode II on the SIF in Equation (25) [47] as follows:
Δ K = Δ K I 2 + 2 Δ K I I 2

3.3. Propagation Direction

Here, the fracture was justified based on the maximum circumferential tensile stress, and it was used to predict the direction of propagation [48] as follows:
γ = a r c c o s 3 K I I 2 + K I K I 2 + 8 K I I 2 K I 2 + 9 K I I 2  
where γ is the direction of crack propagation.
For radially distributed micro-cracks, the propagation angle of macro-crack γ was −1.89°, −70.26°, and −38.86° under pure tension, pure shear, and tension-shear loads, respectively. For randomly distributed micro-cracks, γ was 5.65°, −70.53°, and −45.68° under pure tension, pure shear, and tension-shear loads, respectively. When the micro-cracks were virtual around the tip of the macro-crack, γ was 0°, (−70°, −100°) and −40° [24], respectively. Thus, the propagation angle was unique to a specified loading condition and damage configuration, and the micro-cracks directly affected the propagation path of the macro-crack. This phenomenon was attributed to the attraction effects of the micro-cracks, inclusion, and dislocations on crack propagation [49,50]. In other words, the micro-cracks guided the propagation of the macro-crack, ultimately leading to a fracture. The macro-crack would coalesce with a favorable micro-crack near the tip of macro-crack to reform a new macro-crack.

4. Experiment

The experiment was conducted in this section to verify the results of the developed models based on observation, shown in Figure 4 and Figure 5. The specimen and loading conditions are shown in Figure 4, and U71Mn steel was selected as the material. The specimen was rectangular and was 120 mm in length, 15 mm in width, and 1 mm in height. The macro crack was designed on the edge. It was located in the middle of an edge of the specimen, and its initial length was 2.846 mm. A fluctuating load with a constant amplitude was applied through three-point bending at a loading frequency of 1 Hz. The maximum applied load was 1000 N and the stress ratio was 0.1, recorded using a digital force gauge.
There were numerous micro-cracks near the propagation path of the macro-crack, as shown in Figure 5a. The micro-cracks on the plane are highlighted with red arrows, and the coalescence of the macro-crack and micro-cracks are highlighted with a red circle. It was observed that the macro-crack was deflected and propagated in a zigzag manner; this was validated by the observations in ref. [27]. The zigzag path of propagation was greatly affected by the presence and propagation of micro-cracks; the micro-cracks guided the propagation of the macro-crack, since it would coalesce with a favorable micro-crack to form an extended macro-crack. According to Section 3.3, we found that (1) for radially distributed micro-cracks, the influence of the micro-cracks on macro-crack propagation was relatively small, the changes in the propagation angle of the macro-crack were small; (2) for randomly distributed micro-cracks, the influence of the micro-cracks on macro-crack propagation was relatively large. The macro-crack path would deflect and coalescence with the micro-cracks that were in the propagation direction.

5. Summary and Conclusions

This paper aimed to develop a CAE model for a digital twin to predict the fatigue life of selected engineering materials. Fatigue failure is a major concern in applications involving metal products [4]. In this paper, the solutions of crack propagation for a plane that contained a macro-crack and a number of micro-cracks under various loading conditions were obtained based on highly efficient theoretical models. The effects of micro-cracks on macro-crack growth were analyzed, and the following findings were obtained:
(1)
The propagation of a macro-crack and micro-cracks depended mainly on the failure configurations of the micro-cracks in front of the macro-crack as well as the loading conditions.
(2)
Micro-cracks with a small inclination angle (−25° < β < 25°) required less external stress to activate propagation.
(3)
Randomly distributed micro-cracks affected the propagation of the micro-cracks and the macro-crack more significantly than radially distributed micro-cracks.
(4)
Micro-cracks affected the propagation of the macro-crack primarily due to the pure shear load, followed by the tension-shear load and the pure tension load.
(5)
The presence of micro-cracks guided the propagation path of the macro-crack. The macro-crack would coalesce with a favorable micro-crack to form an extended macro-crack.
The influence of micro-cracks on macro-crack propagation is an important subject in fracture mechanics; however, due to the complex calculations of the mutual superposition of micro-cracks, researchers may choose to ignore or approximate the relationship among micro-cracks to predigest the calculation. This approach may create a bottleneck when the number of micro-cracks increases significantly, which can cause an inaccurate result. In the future, we will continue studying the prediction of micro-cracks on macro-crack paths, considering a large number of micro-cracks near the macro-crack tip based on the developed modal, neural network models, or a machine learning method. We also plan to investigate the impacts of micro-crack parameters (size, initial angle, distance between cracks, etc.) on propagation behavior. We expect that the results will be utilized to optimize anti-fracture designs by predicting the fracture behaviors of engineering materials.

Author Contributions

Conceptualization, X.L. and G.L.; methodology, X.L. and G.L.; software, X.L.; validation, X.L., G.L. and Z.B.; formal analysis, X.L.; resources, G.L.; writing—original draft preparation, X.L. and Z.B.; writing—review and editing, X.L. and Z.B.; supervision, G.L.; project administration, G.L.; funding acquisition, X.L. and G.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the Talent Introduction Project of Chengdu Technological University (2022RC005).

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Bi, Z.M.; Jin, Y.; Markopoulos, P.; Zhang, W.J.; Wang, L. Internet of Things (IoT) and Big Data Analytics (BDA) for Digital Manufacturing. Int. J. Prod. Res. 2023, 62, 4004–4021. [Google Scholar] [CrossRef]
  2. Bi, Z.M.; Zhang, W.J.; Wu, C.; Luo, C.; Xu, L. Generic Design Methodology for Smart Manufacturing Systems for a Practical Perspective, Part II-Systematic Designs of Smart Manufacturing Systems. Machines 2021, 9, 208. [Google Scholar] [CrossRef]
  3. Bi, Z.M.; Zhang, W.J.; Wu, C.; Luo, C.; Xu, L. Generic Design Methodology for Smart Manufacturing Systems for a Practical Perspective, Part I–Digital Triad Concept and Its Application as a System Reference Model. Machines 2021, 9, 207. [Google Scholar] [CrossRef]
  4. Bi, Z.M. Petersons Stress Concentration Factors, 4th ed.; John Wiley & Sons: New York, NY, USA, 2021. [Google Scholar]
  5. Zhang, W.; Liu, Y. Investigation of incremental fatigue crack growth mechanisms using in situ SEM testing. Int. J. Fatigue 2012, 42, 14–23. [Google Scholar] [CrossRef]
  6. Ronevich, J.A.; Somerday, B.P.; San Marchi, C.W. Effects of microstructure banding on hydrogen assisted fatigue crack growth in X65 pipeline steels. Int. J. Fatigue 2016, 82, 497–504. [Google Scholar] [CrossRef]
  7. Li, X.T.; Jiang, X.Y. Effect of a Micro-crack on the Edge Macro-crack Propagation Rate and Path Under Mixed Loads. Acta Mech. Solida Sin. 2019, 32, 517–532. [Google Scholar] [CrossRef]
  8. Petrova, V.; Tamuzs, V.; Romalis, N. A Survey of Macro-Microcrack Interaction Problems. Appl. Mech. Rev. 2000, 53, 117–146. [Google Scholar] [CrossRef]
  9. Petrova, V.; Schmauder, S. Mathematical modelling and thermal stress intensity factors evaluation for an interface crack in the presence of a system of cracks in functionally graded/homogeneous biomaterials. Comput. Mater. Sci. 2012, 52, 171–177. [Google Scholar] [CrossRef]
  10. Petrova, V.; Schmauder, S. Thermal fracture of a functionally graded/homogeneous bimaterial with system of cracks. Theor. Appl. Fract. Mech. 2011, 55, 148–157. [Google Scholar] [CrossRef]
  11. Kachanov, M. A simple technique of stress analysis in elastic solids with many cracks. Int. J. Fract. 1985, 28, 11–19. [Google Scholar] [CrossRef]
  12. Kachanov, M. Elastic Solids with Many Cracks: A Simple Method of Analysis. Int. J. Solids Struct. 1987, 23, 23–43. [Google Scholar] [CrossRef]
  13. Tamuzs, V.; Romalis, N.; Petrova, V. Fracture of Solids with Microdefects; Nijhoff: Leiden, The Netherlands, 2000. [Google Scholar]
  14. Tamuzs, V.P.; Petrova, V.E. On macrocrack–microdefect interaction. Int. Appl. Mech. 2002, 38, 1157–1177. [Google Scholar] [CrossRef]
  15. Li, X.; Jiang, X.; Xu, L.; Yang, H.; Zhang, Y. Solution of an inclined crack in a finite plane and a new criterion to predict fatigue crack propagation. Int. J. Mech. Sci. 2016, 119, 29–36. [Google Scholar] [CrossRef]
  16. Li, X.; Xu, L.; Jiang, X. Influence of a micro-crack on the finite macro-crack. Eng. Fract. Mech. 2017, 177, 95–103. [Google Scholar] [CrossRef]
  17. Li, X.; Xu, L.; Yang, H.; Jiang, X. Effect of micro-cracks on plastic zone ahead of the macro-crack tip. J. Mater. Sci. 2017, 52, 13490–13503. [Google Scholar]
  18. Loehnert, S.; Belytschko, T. Crack shielding and amplification due to multiple microcracks interacting with a macrocrack. Int. J. Fract. 2017, 145, 1–8. [Google Scholar] [CrossRef]
  19. Sutula, D.; Kerfriden, P.; Dam, T.V.; Bordas, S.P.A. Minimum energy multiple crack propagation, Part III: XFEM computer implementation and applications. Eng. Fract. Mech. 2018, 191, 257–276. [Google Scholar] [CrossRef]
  20. Gong, S.X.; Horii, H. General solution to the problem of microcracks near the tip of a main crack. J. Mech. Phys. Solids 1989, 37, 27–46. [Google Scholar] [CrossRef]
  21. Roberts, R.; Erdogan, F. Closure to “Discussion of ‘The Effect of Mean Stress on Fatigue Crack Propagation in Plates Under Extension and Bending’”. J. Basic Eng. 1968, 90, 414–417. [Google Scholar] [CrossRef]
  22. El Haddad, M.H.; Topper, T.H.; Smith, K.N. Prediction of non propagating cracks. Eng. Fract. Mech. 1979, 11, 573–584. [Google Scholar] [CrossRef]
  23. Liu, J.Y.; Bao, W.J.; Zhao, J.Y.; Zhou, C.Y. Fatigue Crack Growth Behavior of CP-Ti Cruciform Specimens with Mixed Mode I-II Crack under Biaxial Loading. Materials 2022, 15, 1926. [Google Scholar] [CrossRef] [PubMed]
  24. Khan, S.M.A.; Khraisheh, M.K. Analysis of mixed mode crack initiation angles under various loading conditions. Eng. Fract. Mech. 2000, 67, 397–419. [Google Scholar] [CrossRef]
  25. Wang, L.Y.; Chen, W.Z.; Tan, X.Y.; Yang, J.P. The impact of various crack geometrical parameters on stress field over tip under different mixed loading conditions and inclination angle. Theor. Appl. Fract. Mech. 2019, 102, 239–254. [Google Scholar] [CrossRef]
  26. Ghaffari, M.A.; Pahl, E.; Xiao, E. Three-dimensional fatigue crack initiation and propagation analysis of a gear tooth under various load conditions and fatigue life extension with boron/epoxy patches. Eng. Fract. Mech. 2015, 135, 126–146. [Google Scholar] [CrossRef]
  27. Petel, A.; Jager, A.; Babai, D.; Jopp, J. Fatigue Crack Growth in a Monocrystal and Its Similarity to Short-Crack Propagation in a Polycrystal of Nickel. Metals 2023, 13, 790. [Google Scholar] [CrossRef]
  28. Liu, Y.; Jiang, P.; Duan, G.; Wang, J. Investigation on the Fatigue Crack Propagation of Medium-Entropy Alloys with Heterogeneous Microstructures. Materials 2022, 15, 6081. [Google Scholar] [CrossRef] [PubMed]
  29. Reymer, P.; Leski, A.; Dziendzikowski, M. Fatigue Crack Propagation Estimation Based on Direct Strain Measurement during a Full-Scale Fatigue Test. Sensors 2022, 22, 2019. [Google Scholar] [CrossRef]
  30. Dong, Q.; Xu, G.; Hu, Y.; Peng, Z. Research on the Residual Strength of Cracked Plate Considering Fatigue Crack Propagation under Cyclic Load. J. Mar. Sci. Eng. 2023, 11, 706. [Google Scholar] [CrossRef]
  31. Li, X.; Yang, H.; Zan, X.; Li, X.; Jiang, X. Effect of a micro-crack on the kinked macro-crack. Theor. Appl. Fract. Mech. 2018, 96, 102–110. [Google Scholar] [CrossRef]
  32. Basoglu, M.; Zerin, Z.; Kefal, A.; Oterkus, E. A computational model of peridynamic theory for deflecting behavior of crack propagation with micro-cracks. Comput. Mater. Sci. 2019, 162, 33–46. [Google Scholar] [CrossRef]
  33. Yates, J.; Zanganeh, M.; Tomlinson, R.; Brown, M.; Garrido, F.D. Crack paths under mixed mode loading. Eng. Fract. Mech. 2018, 75, 319–330. [Google Scholar] [CrossRef]
  34. Liu, G.; Zhou, D.; Guo, J.; Bao, Y.; Han, Z.; Lu, J. Numerical simulation of fatigue crack propagation interacting with micro-defects using multiscale XFEM. Int. J. Fatigue 2018, 109, 70–82. [Google Scholar] [CrossRef]
  35. Zhu, Y.; Yang, J.; Pan, H. Three-Dimension Crack Propagation Behavior of Conical-Cylindrical Shell. Metals 2023, 13, 698. [Google Scholar] [CrossRef]
  36. Ohr, S. An electron microscope study of crack tip deformation and its impact on the dislocation theory of fracture. Mater. Sci. Eng. 1985, 72, 1–35. [Google Scholar] [CrossRef]
  37. Gao, Y.; Zhang, C.; Xiong, X.; Zheng, Z.; Zhu, M. Intergranular corrosion susceptibility of a novel Super304H stainless steel. Eng. Fail. Anal. 2012, 24, 26–32. [Google Scholar] [CrossRef]
  38. Bai, Y.; Wang, H.; Xia, M.; Ke, F. Statistical meso-mechanics of solid, linking coupled multiple space and time scales. Appl. Mech. Rev. 2015, 58, 372–388. [Google Scholar] [CrossRef]
  39. Muskhelishvili, N.I. Some Basic Problems of the Mathematical Theory of Elasticity; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  40. Pook, L. Why Metal Fatigue Matters; Springer: Berlin/Heidelberg, Germany; Dordrecht, The Netherlands, 2003. [Google Scholar]
  41. Milella, P. Fatigue and Corrosion in Metals; Springer: Berlin/Heidelberg, Germany; Dordrecht, The Netherlands, 2013. [Google Scholar]
  42. Miller, K. Short Fatigue Cracks; Springer: Berlin/Heidelberg, Germany; Dordrecht, The Netherlands, 1989. [Google Scholar]
  43. Stephens, I. Metal Fatigue in Engineering; Wiley: Hoboken, NJ, USA, 1980. [Google Scholar]
  44. Haddad, M.H.; Smith, K.N.; Topper, T.H. Fatigue crack propagation of short cracks. J. Eng. Mater. Technol. Trans. Asme 1979, 101, 42–46. [Google Scholar] [CrossRef]
  45. Wang, W.; Liu, Q. Study on fatigue crack growth rate of PD3 and U71Mn rail steel. J. Mech. Strength 2007, 29, 1026–1029. [Google Scholar]
  46. Wolf, E. Fatigue crack closure under cyclic tension. Eng. Fract. Mech. 1970, 2, 37–45. [Google Scholar] [CrossRef]
  47. Tanaka, K. Fatigue crack propagation from a crack inclined to the cyclic tensile axis. Eng. Fract. Mech. 1974, 6, 493–507. [Google Scholar] [CrossRef]
  48. Erdogan, F.; Sih, G. On the crack extension in plates under plane loading and transverse shear. J. Basic Eng. 1963, 85, 519–525. [Google Scholar] [CrossRef]
  49. Bouchard, P.O.; Bay, F.; Chastel, Y. Numerical modelling of crack propagation: Automatic remeshing and comparison of different criteria. Comput. Methods Appl. Mech. Eng. 2003, 192, 3887–3908. [Google Scholar] [CrossRef]
  50. Miranda, A.; Meggiolaro, M.; Castro, J.; Martha, L.; Bittencourt, T. Fatigue life and crack path predictions in generic 2D structural components. Eng. Fract. Mech. 2003, 70, 1259–1279. [Google Scholar] [CrossRef]
Figure 1. A plane with one macro-crack and a set of micro-cracks subject to mixed loads.
Figure 1. A plane with one macro-crack and a set of micro-cracks subject to mixed loads.
Machines 12 00125 g001
Figure 2. (ac) Decomposition of the original problem into a set of sub-problems, where the dashed line represents the absence of cracks and the solid line represents the real crack.
Figure 2. (ac) Decomposition of the original problem into a set of sub-problems, where the dashed line represents the absence of cracks and the solid line represents the real crack.
Machines 12 00125 g002
Figure 3. The PZ of the Irwin model for mode I.
Figure 3. The PZ of the Irwin model for mode I.
Machines 12 00125 g003
Figure 4. U71Mn steel under a three-point bending load.
Figure 4. U71Mn steel under a three-point bending load.
Machines 12 00125 g004
Figure 5. Scanning electron microscope photograph of crack propagation. (a) Numerous micro-cracks near the macro-crack path; (b) coalescence of the macro-crack and micro-cracks.
Figure 5. Scanning electron microscope photograph of crack propagation. (a) Numerous micro-cracks near the macro-crack path; (b) coalescence of the macro-crack and micro-cracks.
Machines 12 00125 g005
Table 1. Material properties of U71Mn steel [45].
Table 1. Material properties of U71Mn steel [45].
Fatigue Limits Δσ0CmΔKthDYield Limits σsa0
96 MPa4.597 × 10−132.882.2 MPa·m1/2100 μm550 MPa167 μm
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, X.; Li, G.; Bi, Z. Digital Twins to Predict Crack Propagation of Sustainable Engineering Materials under Different Loads. Machines 2024, 12, 125. https://doi.org/10.3390/machines12020125

AMA Style

Li X, Li G, Bi Z. Digital Twins to Predict Crack Propagation of Sustainable Engineering Materials under Different Loads. Machines. 2024; 12(2):125. https://doi.org/10.3390/machines12020125

Chicago/Turabian Style

Li, Xu, Gangjun Li, and Zhuming Bi. 2024. "Digital Twins to Predict Crack Propagation of Sustainable Engineering Materials under Different Loads" Machines 12, no. 2: 125. https://doi.org/10.3390/machines12020125

APA Style

Li, X., Li, G., & Bi, Z. (2024). Digital Twins to Predict Crack Propagation of Sustainable Engineering Materials under Different Loads. Machines, 12(2), 125. https://doi.org/10.3390/machines12020125

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