Next Article in Journal
Correlation between Magnetocaloric Properties and Magnetic Exchange Interaction in Gd54Fe36B10−xSix Amorphous Alloys
Next Article in Special Issue
Special Issue “Fracture Mechanics and Fatigue Damage of Materials and Structures”
Previous Article in Journal
Effect of Shear Stress during Processing on Structure, Morphology, and Properties of Isotactic Polypropylene Nucleated with Silsesquioxane-Based β-Nucleating Agent
Previous Article in Special Issue
Material’s Strength Analysis of the Coupling Node of Axle of the Truck Trailer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Phase Field Approach to Two-Dimensional Quasicrystals with Mixed Mode Cracks

1
State Key Laboratory of Structural Analysis, Optimization and CAE Software for Industrial Equipment, Department of Engineering Mechanics, International Center for Computational Mechanics, Dalian University of Technology, Dalian 116024, China
2
School of Mechanics, Civil Engineering and Architecture, Northwestern Polytechnical University, Xi’an 710072, China
*
Author to whom correspondence should be addressed.
Materials 2023, 16(10), 3628; https://doi.org/10.3390/ma16103628
Submission received: 31 March 2023 / Revised: 5 May 2023 / Accepted: 8 May 2023 / Published: 9 May 2023
(This article belongs to the Special Issue Fracture Mechanics and Fatigue Damage of Materials and Structures)

Abstract

:
Quasicrystals (QCs) are representatives of a novel kind of material exhibiting a large number of remarkable specific properties. However, QCs are usually brittle, and crack propagation inevitably occurs in such materials. Therefore, it is of great significance to study the crack growth behaviors in QCs. In this work, the crack propagation of two-dimensional (2D) decagonal QCs is investigated by a fracture phase field method. In this method, a phase field variable is introduced to evaluate the damage of QCs near the crack. Thus, the crack topology is described by the phase field variable and its gradient. In this manner, it is unnecessary to track the crack tip, and therefore remeshing is avoided during the crack propagation. In the numerical examples, the crack propagation paths of 2D QCs are simulated by the proposed method, and the effects of the phason field on the crack growth behaviors of QCs are studied in detail. Furthermore, the interaction of the double cracks in QCs is also discussed.

1. Introduction

Quasicrystals (QCs) are a new kind of material with perfect long-range order but no periodicity. Due to their unique atomic structures, QCs exhibit many excellent physical properties, such as high hardness, a low friction coefficient, and high resistance. Therefore, QCs are very promising materials for potential applications in corrosion-resistant coatings, hydrogen storage, photovoltaic solar cells, etc. [1]. However, QCs are usually brittle, and consequently, cracks, holes, and other defects will inevitably occur during daily use. If the crack propagates, the QC will fail, which may lead to a catastrophic accident. Therefore, the fracture analysis of QCs is of great importance, and it is very necessary to investigate crack propagation behaviors in cracked QCs.
Plenty of research work has been carried out on the fracture problems of QCs. Li et al. [2] extended the classical linear elasticity fracture mechanics to investigate a decagonal QC with a Griffith crack. The results indicated the phonon and phason stresses at the crack tip exhibit the well-known square root singularity, and the strain energy release rate was obtained. Zhou and Fan [3] developed the plane elasticity theory of two-dimensional (2D) octagonal QCs and obtained the exact analytic solution of a Mode I Griffith crack. Guo and Fan [4] studied the fracture problem of a Mode II Griffith crack in decagonal quasicrystals and obtained the corresponding stress intensity factors and strain energy release rate. Shen and Fan [5] calculated the stress intensity factors for an infinitely long strip of finite height containing two straight, semi-infinite collinear cracks. Li and Fan [6] obtained exact solutions for two semi-infinite collinear cracks in a strip of 1D hexagonal QCs. After that, they [7] further obtained an analytic solution for the elliptic notch problem of the material in icosahedral QCs by using the complex variable method. The solution can be reduced to that of a Griffith crack problem. Fan et al. [8] presented linear, nonlinear, and dynamic fracture problems for different QCs. Li and Liu [9] obtained closed-form expressions for the elastic displacement and stress fields induced by a dislocation in icosahedral QCs. Sladek et al. [10] developed a meshless method based on the local Petrov-Galerkin approach for fracture analysis of decagonal QCs; both static and transient dynamic boundary value problems were considered. After that, they [11] present path-independent integrals for accurate evaluation of energy release and stress intensity factors in decagonal QCs. Wang et al. [12] obtained the phonon-phason coupling field in the half-space, which can be expressed in terms of elementary functions. These solutions could have applications in 3D contact and crack problems in QCs. Li et al. [13,14] derived solutions for elliptical crack and planar crack problems in 2D hexagonal QCs. Li and Shi [15] employed the method of potential function theory to solve plane defect problems originating from two-dimensional decagonal QCs. Zhao et al. [16] derived the fundamental solutions for interface cracks in 1D hexagonal QC coatings under in-plane loads.
For crack propagation in QCs, some investigations have been conducted based on the atomistic model [17,18,19,20,21,22,23,24,25,26,27,28,29]. However, it is inconvenient to apply atomistic simulation to engineering. Therefore, Wang and Ricoeur [30] adopted the finite element method (FEM) to simulate the crack growth in 1D QCs and predicted the crack pattern for different boundary conditions. Fan et al. [31,32,33,34,35,36,37,38,39] investigated crack propagation based on the elastodynamic/hydrodynamic model.
From the above literature review, it is found that increasing attention has been paid to the fracture of QCs. However, most of them were concentrated on the derivation of exact solutions and determination of fracture parameters of cracked QCs, and the investigations on crack propagation were mainly based on the atomistic model [17,18,19,20,21,22,23,24,25,26,27,28,29]. The work on crack growth in QCs based on continuum mechanics is still insufficient. Additionally, in the traditional FEM, the crack topology is modeled by geometry. The remesh process is necessary when simulating crack propagation, which has a huge computational cost. In this paper, a phase field method is introduced to predict the crack propagation path in QCs. Unlike conventional discrete crack models (such as the FEM), the fracture phased field method employs diffusive cracks to avoid an explicit representation of kinematic discontinuities, and therefore the propagating cracks are tracked automatically without additional ad-hoc criteria in the classical Griffith fracture theory [40,41]. In this method, the fracture energy and degraded strain energy of QCs are formulated using the phase field variable. Subsequently, the total potential energy of QCs under the phase field framework is obtained, and the governing equations for the phase field model are derived by means of the Francfort-Marigo variational principle. The phase field evolution equation for QCs is constructed. Finally, the FEM is adopted to solve the phase field governing equations. The phase field variable and phonon/phason displacement of the entire model can be obtained, as well as the reaction force and the crack pattern.
This paper is organized as follows: The basic equations of 2D decagonal QCs are presented in Section 2.1. The phase field model for 2D QCs is formulated in Section 2.2. The finite element implementation for the phase field model is presented in Section 2.3. Several numerical examples are presented in Section 3. Conclusions are summarized in Section 4.

2. Phase Field Method for 2D Decagonal QCs

2.1. The Basic Equations

According to the elasticity of 2D decagonal QCs, the basic equations for the plane problem of QCs in the absence of body force and phason self-action are [42]:
σ x x + τ x y y = 0 τ y x x + σ y y = 0 ,   H x x + H x y y = 0 H y x x + H y y = 0
ε x ε y γ x y = x y y x u x u y ,   ω x ω y ω x y ω y x = x y y x w x w y
σ x σ y τ x y = C 11 C 12 C 12 C 22 C 66 ε x ε y γ x y + R 1 R 1 R 2 R 2 R 1 R 1 R 2 R 2 R 2 R 2 R 1 R 1 ω x ω y ω x y ω y x
H x H y H x y H y x = R 1 R 1 R 2 R 1 R 1 R 2 R 2 R 2 R 1 R 2 R 2 R 1 ε x ε y γ x y + K 1 K 2 K 2 K 1 K 1 K 2 K 2 K 1 ω x ω y ω x y ω y x
where u x and u y are the phonon displacements; w x and w y are the phason displacements; ε x , ε y , and γ x y are the phonon strains; ω x , ω y , ω x y , and ω y x are the phason strains; σ x , σ y , and σ x y are the phonon stresses; H x , H y , H x y , and H y x are the phason stresses; C 11 , C 12 , C 22 , and C 66 are the phonon moduli; K 1 and K 2 are the phason moduli; and R 1 and R 2 are the phonon-phason coupling coefficients.
The strain energy of QCs is:
II s = Ω ψ ε , ω d V = Ω 1 2 ε T σ + 1 2 ω T H d V
where σ = σ x , σ y , σ x y T and H = H x , H y , H x y , H y x T are the stress vectors; ε = ε x , ε y , γ x y T and ε = ω x , ω y , ω x y , ω y x T are the strain vectors; and ψ ε , ω is the strain energy density.

2.2. Phase Field Method

Consider a 1D QC strip with a centered crack, as shown in Figure 1. The fracture energy is:
II c = Ω c G c d S = G c A c
where G c is the critical energy release rate (CERR); Ω c is the crack surface; and A c is the cross-sectional area of the strip.
In a phase field model, the crack is supposed to be surrounded by a diffusive degraded zone, and a phase field variable d is introduced to describe the damage in the diffusive degraded zone:
d = e x l c
where lc is an internal length scale that controls the width of the diffusive zone. Here, it is noted that, when x approaches infinity (±∞), d converges to zero, which indicates that the material is intact; when x is zero, d equals one, which indicates that the material is totally broken and the crack surface is produced. Define a crack’s surface density as functional [43]:
γ d , d = 1 2 l c d 2 + l c 2 d 2
The crack topology can then be described by the phase field variable.
Substituting Equation (8) into Equation (6) yields the fracture energy in the phase field model:
II c = Ω G c γ d , d d V
where Ω is the domain of the overall model. As observed in Equation (9), the integral over the crack path is transformed into the integral over the model. The crack topology is implicitly expressed in the framework of the phase field method. Therefore, it is unnecessary to track the crack’s path while it is propagating.
Similarly, the 2D crack surface density functional can be defined as:
γ d , d = 1 2 l c d 2 + l c 2 d 2
Therefore, the fracture energy of QCs in a plane problem is:
II c = Ω G c γ d , d d V
In the phase field model, the material at the crack is assumed to be softened. Therefore, a degradation function is introduced to evaluate the damage to the material in the diffusive degraded zone [43]:
ω d = 1 d 2 + k
where k is a small positive number to ensure the non-singularity of the matrix. It can be observed that ω d satisfies ω 0 = 1 and ω 1 = 0 .
The strain energy for QCs is modified as follows:
II s = Ω ω d ψ ε , ω d V
It should be pointed out that the crack will not propagate if the crack face is under compression. Therefore, the strain energy density ψ should be decomposed into a tensile part and a compressed part, and the energy degradation only occurs on the tensile part [43]:
II s = Ω ω d ψ + + ψ d V
where ψ + = λ ε 1 + ε 2 + ε 3 + 2 / 2 + μ ε 1 + 2 + ε 2 + 2 + ε 3 + 2 and ψ = ψ ψ + are the tensile and compressed strain energy densities, respectively; + is defined as x + = x + x / 2 ; and λ and μ are the Lamé constants.
The potential energy of the overall model contains the strain energy, fracture energy, and potential energy of the external force:
II p = II s + II c II e   = Ω ω d ψ + + ψ d V + Ω G c γ d , d d V Ω u T t u d S Ω w T t w d S
where u = u x , u y T ; w = w x , w y T ; and t u and t w are the external forces in the phonon and phason fields, respectively.
The Francfort-Marigo variational principle states that the real displacement u and the phase field variable d will minimize the potential energy.
δ II p = Ω u ω d σ i j + n i δ u j + σ i j n i δ u j d S + Ω w ω d H i j + n i δ w j + H i j n i δ w j d S Ω ω d σ i j + x i δ u j + σ i j x i δ u j d V Ω ω d H i j + x i δ w j + H i j x i δ w j d V Ω ω d ψ + δ d d V + Ω G c l c d δ d d V + Ω G c l c d x i n i δ d d S Ω G c l c 2 d x i 2 δ d d V Ω δ u j t j u d S Ω δ w j t j w d S = 0
where σ i j + and H i j + are the stresses induced by stretch, while σ i j and H i j are induced by compression. Equation (16) is valid for arbitrary δ u i , δ w i , and δ d . Hence, the governing equations for the phase field model of QCs are:
ω d σ i j + + σ i j = 0 ,   Ω
ω d H i j + + H i j = 0 ,   Ω
2 d 2 ψ + + G c l c d G c l c Δ d = 0 ,   Ω
ω d σ i j + + σ i j n i = t j u ,   Ω
ω d H i j + + H i j n i = t j w ,   Ω
d , i n i = 0 ,   Ω
It should be mentioned that crack growth is an irreversible process. Therefore, Equation (19) should be modified by considering the history of the load [44]:
2 d 2 H + G c l c d G c l c Δ d = 0
where H = max 0 , t ψ + is a history variable that is the maximum strain energy during the crack propagation. This history variable ensures that the crack face does not close under compression. Equation (23) is the evolution law of the crack phase field for QCs. Cracks grow only if this equation is valid.

2.3. Finite Element Implementation

Due to the strong nonlinearity of the governing equations in the phase field model, the FEM is often adopted to solve the problem. In the FEM, the phonon/phason displacements and the phase field variable are approximated by the shape functions:
u x , u y T = i = 1 4 N i u u i
w x , w y T = i = 1 4 N i w w i
d = i = 1 4 N i d i
where Ni is the shape function; N i u = N i w = diag N i , N i ; and u i , w i , and d are the nodal phonon and phason displacements and phase field variables, respectively.
Therefore, the phonon and phason strains and the gradient of the phase field variable are, respectively:
ε x , ε y , γ x y T = i = 1 4 B i u u i
ω x , ω y , ω x y , ω y x T = i = 1 4 B i w w i
d = i = 1 4 B i d d i
where B i u = N i / x N i / y N i / y N i / x T , B i w = N i / x N i / y N i / y N i / x T , and B i d = N i / x N i / y T . Substituting Equations (24)–(29) into Equation (15) yields the residuals:
r a = Ω N T t d S Ω 1 d 2 B T D Q C B a d V = 0
r d = Ω 2 d 1 H u N d T d V Ω G c c 0 2 d l c N d T + 2 l c B d T d d V = 0
where a = u x 1 , u y 1 , , u x 4 , u y 4 , w x 1 , w y 1 , , w x 4 , w y 4 T ; N = N 0 0 2 × 8 0 2 × 8 N 0 is the shape function matrix where N 0 = N 1 0 N 2 0 N 3 0 N 4 0 0 N 1 0 N 2 0 N 3 0 N 4 ; B = L N is the strain matrix where L = x 0 y 0 0 0 0 0 y x 0 0 0 0 0 0 0 x 0 y 0 0 0 0 0 y 0 x T ; and t = t x u , t y u , t x w , t y w T is the load vector.
Equations (30) and (31) can be solved by this iteration method:
a n + 1 d n + 1 = a n d n + K n a a K n a d K n d a K n d d 1 r n a r n d
where K a a = Ω ω d B T D B d V , K a d = K d a T = Ω 2 d 1 B T σ N d V , and K d d = Ω 2 H + G c l N d T N d d V + Ω l G c B d T B d d V .
Due to the high nonlinearity of Equation (32), there is a convergence problem using the classic Newton iteration. A highly robust staggered algorithm is usually adopted to solve Equation (32) [41,44]. In this algorithm, one of the two unknowns (displacement and phase field variable) is assumed to be unchanged while solving the other unknown during one iteration, which yields:
a n + 1 = a n + K n a a 1 r n a
d n + 1 = d n + K n d d 1 r n d
Finally, the nodal phonon/phason displacements and the phase field variables can be calculated by Equations (33) and (34), respectively.

3. Numerical Results

In this section, a few numerical examples are presented to illustrate the application of the phase field method to the fracture of 2D decagonal QCs. The material parameters are tabulated in Table 1 [45,46,47].
According to Fan’s work [42], the expression of CERR is:
G c = λ K 1 + K 2 + 2 R 1 2 + R 2 2 8 λ + M c K I C 2
where M = C 11 C 12 / 2 and c = M K 1 + K 2 2 R 1 2 + R 2 2 . The fracture toughness of Al-Ni-Co QCs is K I C = 1 MPa m . Therefore, the CERR in this paper is selected as G c = 5.56 × 10 4 N / mm , according to the material constants in Table 1. Although no experimental method has been reported to apply a constant phason displacement on the surface of QCs, some investigations [48,49] reveal that some ways can cause the disorder of the phason field. Therefore, different phason displacement loads are considered in the following calculation to investigate the effect of the phason field on the fracture behavior of QCs.

3.1. The Rectangular QCs with Edge Crack

A rectangular QC model is considered, as shown in Figure 2. The geometry of the model is width W, height L, and crack length a. The upper edge is constrained to have the same displacement in phonon and phason fields, and the lower edge is only constrained in the vertical direction. A concentrated phonon force P σ is applied to the upper edge. The angle between P σ and the horizon is φ .
At the current stage, the crack growth of 2D QCs has not been reported in the open literature. To demonstrate the accuracy of the present method, an elastic material is selected by degenerating all the material constants in phason field, i.e., K 1 = K 2 = R 1 = R 2 = 0 . The geometrical parameters are taken as a = 0.5 mm and W = L = 1 mm . The elastic material constants come from Ref. [44]. The variation of the vertical concentrated force P σ ( φ = π / 2 ) versus the displacement u is plotted in Figure 3. As observed, the present results are in excellent agreement with the reference data from Ref. [44]. Furthermore, it is noted that the peak force P σ shows an increasing trend as the length scale lc decreases, which indicates the crack grows slowly as the length scale declines. In the previous study, to ensure the high resolution of the crack topology, the minimum size of the element was required to meet the condition of h < l c / 2 [43]. Therefore, in the following calculations, the length scale is selected as l c = 1 mm , and the element size is selected as h < l c / 5 .
Subsequently, the effect of the phason field on the crack growth is studied in Figure 2 by considering an initial phason displacement w 0 on the upper edge of the model. The parameters are selected as a = 50 mm and W = L = 100 mm . The variation of the applied force versus the displacement for different w 0 values is illustrated in Figure 4. The results show that the peak force decreases as the initial phason displacement increases, which indicates that the initial stretch phason load will undermine the strength of the model. In addition, as the initial phason load gets larger, the crack grows slower. The crack pattern with w 0 = 6 × 10 4 mm is plotted in Figure 5.
Finally, the QC model subjected to a concentrated shear load at its upper edge is considered. The variation of the shear phonon force P σ ( φ = π ) versus the horizontal displacement at the upper edge is illustrated in Figure 6. Two peak forces were observed, and they decrease as the initial vertical phason displacement w 0 increases. The crack patterns with different initial phason displacements are shown in Figure 7. Clearly, the increasing initial phason displacement will lead to a significant crack deflection. Therefore, it is concluded that the phason load is a key influencing factor for the peak force and crack patterns under the applied shear loads.

3.2. The Rectangular QCs with an Internal Crack

In this example, a rectangular QC with an internal crack is considered in Figure 8. The length of the crack is a, and the angle between the crack and the horizon is θ . The middle of the crack is centered in the model. The boundary conditions of the QC are the same as those in Section 3.1. The parameters are selected as a = 25 mm and W = L = 100 mm . The QC model is subjected to a vertical initial phason displacement w 0 at the upper edge.
Firstly, the QC model is subjected to a tensile load in the phonon field ( φ = π / 2 ). The peaks of the phonon forces for different angles θ are shown in Figure 9. It can be observed that the peak force monotonously increases with the increase in θ , while it shows an opposite trend as the initial phason displacement increases. The force-displacement relation for different angles θ is shown in Figure 10. Clearly, the slope of the force-displacement curve increases as the angle θ increases. The crack patterns with different angles are plotted in Figure 11. As observed, the crack grows along a straight line to the edge of the model when subjected to a tensile load.
Subsequently, the model is subjected to shear loads in the phonon field at the upper edge. The peak phonon forces for different angles θ are illustrated in Figure 12. As observed, the peak force first increases and then decreases as the angle θ increases. The crack patterns are plotted in Figure 13. Similar to the observation in Figure 7, the increasing initial phason displacement has a significant influence on the crack propagation path.

3.3. The Rectangular QCs with Double Cracks

As the last example, the QC model with double cracks is investigated in Figure 14. The parameters are selected as a 1 = a 2 = 25 mm and W = L = 100 mm . The left crack is horizontal. The angle between the right crack and the horizon is θ . The lengths of the left and right cracks are a1 and a2, respectively. Point B is located at the center of the right crack. The distance between the right edge of the model and point B is Δ 1 = 27.5 mm . The distance between the left crack and point B is Δ 2 . The model is subjected to a tensile phonon/phason displacements. The peak forces for different angles θ with Δ 2 = 32.5 mm are plotted in Figure 15. It can be found that, as θ increases, the peak force monotonously increases for w 0 = 0 mm and w 0 = 0.0004 mm , while it increases first and then decreases for w 0 = 0.0008 mm . The crack patterns for different angles θ and distances Δ 2 are illustrated in Figure 16 and Figure 17. As observed, the interaction of two cracks has a big contribution to their crack propagation path. The two cracks eventually merge together as the applied load increases.

4. Conclusions

In this paper, a phase field model is developed to predict crack propagation in 2D decagonal QCs. The contribution of the phason field to the potential energy is considered, and the evolution law of the crack phase field for QCs is established. Therefore, the crack topology in QCs is described by the phase field variable, which can be solved by the FEM. In this manner, the crack propagation in QCs can be accurately simulated without remeshing, and the evolution of the crack can be vividly illustrated. Numerical examples illustrate that the proposed model can predict accurate results for the crack propagation of QCs, and the phason field has a big contribution to both the force-displacement relation and the crack pattern.

Author Contributions

Conceptualization, Z.Z. and Z.Y.; methodology, Z.Y. and T.L.; data curation, T.L.; writing—preparation of original draft, T.L., Z.Y. and C.X.; writing—proofreading and editing, Z.Z. and X.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Fundamental Research Funds for the Central Universities, [grant number DUT21LK35, DUT22LK16]; State Key Laboratory of Structural Analysis, Optimization and CAE Software for Industrial Equipment, [grant number GZ22109].

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Maciá Barber, E. Quasicrystals: Fundamentals and Applications, 1st ed.; CRC Press: Boca Raton, MA, USA, 2020. [Google Scholar]
  2. Li, X.F.; Fan, T.Y.; Sun, Y.F. A decagonal quasicrystal with a Griffith crack. Philos. Mag. 1999, 79, 1943–1952. [Google Scholar]
  3. Zhou, W.M.; Fan, T.Y. Plane elasticity problem of two-dimensional octagonal quasicrystals and crack problem. Chin. Phys. 2001, 10, 743–747. [Google Scholar]
  4. Guo, Y.C.; Fan, T.Y. A mode-II Griffith crack in decagonal quasicrystals. Appl. Math. Mech. 2001, 22, 1311–1317. [Google Scholar] [CrossRef]
  5. Shen, D.W.; Fan, T.Y. Exact solutions of two semi-infinite collinear cracks in a strip. Eng. Fract. Mech. 2003, 70, 813–822. [Google Scholar] [CrossRef]
  6. Li, L.H.; Fan, T.Y. Exact solutions of two semi-infinite collinear cracks in a strip of one dimensional hexagonal quasicrystal. Appl. Math. Comput. 2008, 196, 1–5. [Google Scholar] [CrossRef]
  7. Li, L.H.; Fan, T.Y. Complex variable method for plane elasticity of icosahedral quasicrystals and elliptic notch problem. Sci. China Phys. Mech. 2008, 51, 773–780. [Google Scholar] [CrossRef]
  8. Fan, T.Y.; Tang, Z.Y.; Chen, W.Q. Theory of linear, nonlinear and dynamic fracture for quasicrystals. Eng. Fract. Mech. 2012, 82, 185–194. [Google Scholar] [CrossRef]
  9. Li, L.H.; Liu, G.T. Stroh formalism for icosahedral quasicrystal and its application. Phys. Lett. A 2012, 376, 987–990. [Google Scholar] [CrossRef]
  10. Sladek, J.; Sladek, V.; Krahulec, S.; Zhang, C.; Wünsche, M. Crack analysis in decagonal quasicrystals by the MLPG. Int. J. Fract. 2013, 181, 115–126. [Google Scholar] [CrossRef]
  11. Sladek, J.; Sladek, V.; Atluri, S.N. Path-independent integral in fracture mechanics of quasicrystals. Eng. Fract. Mech. 2015, 140, 61–71. [Google Scholar] [CrossRef]
  12. Wang, T.; Li, X.Y.; Zhang, X.; Müller, R. Fundamental solutions in a half space of two-dimensional hexagonal quasicrystal and their applications. J. Appl. Phys. 2015, 117, 154904. [Google Scholar] [CrossRef]
  13. Li, Y.; Fan, C.; Qin, Q.H.; Zhao, M. Closed-form solutions of an elliptical crack subjected to coupled phonon–phason loadings in two-dimensional hexagonal quasicrystal media. Math. Mech. Solids 2019, 24, 1821–1848. [Google Scholar] [CrossRef]
  14. Li, Y.; Zhao, M.; Qin, Q.H.; Fan, C. Analysis solution method for 3D planar crack problems of two-dimensional hexagonal quasicrystals with thermal effects. Appl. Math. Model. 2019, 69, 648–664. [Google Scholar] [CrossRef]
  15. Li, W.; Shi, Y.Q. Extension of elastic models to decagonal quasicrystals. Crystals 2020, 10, 469. [Google Scholar] [CrossRef]
  16. Zhao, M.H.; Fan, C.Y.; Lu, C.S.; Dang, H.Y. Analysis of interface cracks in one-dimensional hexagonal quasi-crystal coating under in-plane loads. Eng. Fract. Mech. 2021, 243, 107534. [Google Scholar] [CrossRef]
  17. Stadler, J.; Mikulla, R.; Trebin, H.R. IMD: A software package for molecular dynamics studies on parallel computers. Int. J. Mod. Phys. C 1997, 8, 1131–1140. [Google Scholar] [CrossRef]
  18. Trebin, H.R.; Mikulla, R.; Stadler, J.; Schaaf, G.; Gumbsch, P. Molecular dynamics simulations of crack propagation in quasicrystals. Comput. Phys. Commun. 1999, 121–122, 536–539. [Google Scholar] [CrossRef]
  19. Schaaf, G.D.; Roth, J.; Trebin, H.R.; Mikulla, R. Numerical simulation of dislocation motion in three-dimensional icosahedral quasicrystals. Philos. Mag. A 2000, 80, 1657–1668. [Google Scholar] [CrossRef]
  20. Krdzalic, G.; Brunelli, M.; Trebin, H.R. Temperature dependence of dislocation motion and crack propagation in a two-dimensional binary model quasicrystal. MRS Online Proc. Libr. 2001, 643, K7.1.1-6. [Google Scholar] [CrossRef]
  21. Rudhart, C.; Gumbsch, P.; Trebin, H.R. From crystalline to glassy: Crack propagation modes in decagonal quasicrystals. MRS Online Proc. Libr. 2003, 805, 312–317. [Google Scholar] [CrossRef]
  22. Rudhart, C.; Rösch, F.; Gähler, F.; Roth, J.; Trebin, H.R. Crack Propagation in Icosahedral Model Quasicrystals; Springer: Berlin/Heidelberg, Germany, 2003; pp. 107–116. [Google Scholar]
  23. Schaaf, G.D.; Roth, J.; Trebin, H.R. Dislocation motion in icosahedral quasicrystals at elevated temperatures: Numerical simulation. Philos. Mag. 2003, 83, 2449–2465. [Google Scholar] [CrossRef]
  24. Rosch, F.; Rudhart, C.; Gumbsch, P.; Trebin, H.R. Cleavage planes of icosahedral quasicrystals: A molecular dynamics study. Mat. Res. Soc. Symp. Proc. 2004, 805, 329–334. [Google Scholar] [CrossRef]
  25. Rudhart, C.; Trebin, H.R.; Gumbsch, P. Crack propagation in perfectly ordered and random tiling quasicrystals. J. Non-Cryst. Solids 2004, 334–335, 453–456. [Google Scholar] [CrossRef]
  26. Rösch, F.; Rudhart, C.; Roth, J.; Trebin, H.R.; Gumbsch, P. Dynamic fracture of icosahedral model quasicrystals: A molecular dynamics study. Phys. Rev. B 2005, 72, 014128. [Google Scholar] [CrossRef]
  27. Rudhart, C.; Gumbsch, P.; Trebin, H.R. Temperature dependence of crack propagation in a two-dimensional model quasicrystal. Philos. Mag. 2005, 85, 3259–3272. [Google Scholar] [CrossRef]
  28. Engel, M.; Umezaki, M.; Trebin, H.R.; Odagaki, T. Dynamics of particle flips in two-dimensional quasicrystals. Phys. Rev. B 2010, 82, 134206. [Google Scholar] [CrossRef]
  29. Lipp, H.; Engel, M.; Sonntag, S.; Trebin, H.R. Phason dynamics in one-dimensional lattices. Phys. Rev. B 2010, 81, 064302. [Google Scholar] [CrossRef]
  30. Wang, Z.; Ricoeur, A. Numerical crack path prediction under mixed-mode loading in 1D quasicrystals. Theor. Appl. Fract. Mech. 2017, 90, 122–132. [Google Scholar] [CrossRef]
  31. Sun, Z.F.; Fan, T.Y.; Wu, X.F. Convolution of the impact three-dimensional elasto-dynamics and dynamic stress intensity factor for an elliptic crack. Acta. Mech. Sin. 2002, 18, 302–308. [Google Scholar]
  32. Wang, X.F.; Fan, T.Y.; Zhu, A.Y. Dynamic behaviour of the icosahedral Al-Pd-Mn quasicrystal with a Griffith crack. Chin. Phys. B 2009, 18, 709–714. [Google Scholar]
  33. Zhu, A.Y.; Fan, T.Y. Dynamic crack propagation in decagonal Al-Ni-Co quasicrystal. J. Phys.-Condes. Matter 2008, 20, 295217. [Google Scholar] [CrossRef]
  34. Yin, Z.H.; Fan, T.Y.; Zhu, A.Y. Dynamic crack propagation in five-fold symmetry quasicrystals. Mod. Phys. Lett. B 2009, 23, 1509–1518. [Google Scholar] [CrossRef]
  35. Qiao, L.P.; Wu, L.; Fan, T.Y. Dynamic response of an icosahedral quasi-crystalline medium with a Griffith crack under mechanical loadings. Adv. Mech. Eng. 2017, 9, 1–12. [Google Scholar] [CrossRef]
  36. Li, W.; Fan, T.Y. Elasto-dynamics of quasicrystals. Crystals 2016, 6, 152. [Google Scholar] [CrossRef]
  37. Cheng, H.; Fan, T.Y.; Wei, H. Phonon–phason dynamics and hydrodynamics of fivefold and tenfold symmetry quasicrystals. Acta Mech. 2017, 228, 136–1372. [Google Scholar] [CrossRef]
  38. Tang, Z.Y.; Fan, T.Y. Three-dimensional equations of generalized dynamics of 18-fold symmetry soft-matter quasicrystals. Mod. Phys. Lett. B 2020, 34, 2050109. [Google Scholar] [CrossRef]
  39. Fan, T.Y.; Tang, Z.Y. Three-dimensional generalized dynamics of soft-matter quasicrystals. Adv. Mater. Sci. Eng. 2020, 2020, 4875854. [Google Scholar] [CrossRef]
  40. Ambati, M.; Gerasimov, T.; De Lorenzis, L. A review on phase-field models of brittle fracture and a new fast hybrid formulation. Comput. Mech. 2015, 55, 383–405. [Google Scholar] [CrossRef]
  41. Bui, T.Q.; Hu, X.F. A review of phase-field models, fundamentals and their applications to composite laminates. Eng. Fract. Mech. 2021, 248, 107705. [Google Scholar] [CrossRef]
  42. Fan, T.Y. Mathematical Theory of Elasticity of Quasicrystals and Its Applications; Springer: Singapore, 2016. [Google Scholar]
  43. Miehe, C.; Welschinger, F.; Hofacker, M. Thermodynamically consistent phase-field models of fracture: Variational principles and multi-field FE implementations. Int. J. Numer. Methods Eng. 2010, 83, 1273–1311. [Google Scholar] [CrossRef]
  44. Miehe, C.; Hofacker, M.; Welschinger, F. A phase field model for rate-independent crack propagation: Robust algorithmic implementation based on operator splits. Comput. Meth. Appl. Mech. Eng. 2010, 199, 2765–2778. [Google Scholar] [CrossRef]
  45. Chernikov, M.A.; Ott, H.R.; Bianchi, A.; Migliori, A.; Darling, T.W. Elastic moduli of a single quasicrystal of decagonal Al-Ni-Co: Evidence for transverse elastic isotropy. Phys. Rev. Lett. 1998, 80, 321–324. [Google Scholar] [CrossRef]
  46. Jeong, H.C.; Steinhardt, P.J. Finite-temperature elasticity phase transition in decagonal quasicrystals. Phys. Rev. B 1993, 48, 9394–9403. [Google Scholar] [CrossRef]
  47. Edagawa, K. Phonon–phason coupling in decagonal quasicrystals. Philos. Mag. 2007, 87, 2789–2798. [Google Scholar] [CrossRef]
  48. Jaric, M.V.; Nelson, D.R. Diffuse scattering from quasicrystals. Phys. Rev. B 1988, 37, 4458–4472. [Google Scholar] [CrossRef] [PubMed]
  49. Coddens, G.; Bellissent, R.; Calvayrac, Y.; Ambroise, J.P. Evidence for phason hopping in icosahedral AlFeCu quasi-crystals. Europhys. Lett. 1991, 16, 271–276. [Google Scholar] [CrossRef]
Figure 1. A QC strip with a centered crack.
Figure 1. A QC strip with a centered crack.
Materials 16 03628 g001
Figure 2. A rectangular QC with an edge crack.
Figure 2. A rectangular QC with an edge crack.
Materials 16 03628 g002
Figure 3. Variation of the force P σ versus the displacement u [44].
Figure 3. Variation of the force P σ versus the displacement u [44].
Materials 16 03628 g003
Figure 4. Force-displacement relations for different initial phason loads w 0 .
Figure 4. Force-displacement relations for different initial phason loads w 0 .
Materials 16 03628 g004
Figure 5. The crack pattern with w 0 = 6 × 10 4 mm .
Figure 5. The crack pattern with w 0 = 6 × 10 4 mm .
Materials 16 03628 g005
Figure 6. The shear force-displacement relations for different initial tensile phason loads.
Figure 6. The shear force-displacement relations for different initial tensile phason loads.
Materials 16 03628 g006
Figure 7. The crack pattern with shear load for different initial phason loads (edge crack): (a) w 0 = 0 mm ; (b) w 0 = 4 × 10 4 mm ; and (c) w 0 = 8 × 10 4 mm .
Figure 7. The crack pattern with shear load for different initial phason loads (edge crack): (a) w 0 = 0 mm ; (b) w 0 = 4 × 10 4 mm ; and (c) w 0 = 8 × 10 4 mm .
Materials 16 03628 g007
Figure 8. A rectangular QC with an internal crack.
Figure 8. A rectangular QC with an internal crack.
Materials 16 03628 g008
Figure 9. The peak phonon forces for different angles θ (edge crack).
Figure 9. The peak phonon forces for different angles θ (edge crack).
Materials 16 03628 g009
Figure 10. The force-displacement relation for different angles θ .
Figure 10. The force-displacement relation for different angles θ .
Materials 16 03628 g010
Figure 11. The crack patterns with w 0 = 6 × 10 4 mm for different angles θ : (a) θ = 0 ; (b) θ = π / 6 ; and (c) θ = π / 3 .
Figure 11. The crack patterns with w 0 = 6 × 10 4 mm for different angles θ : (a) θ = 0 ; (b) θ = π / 6 ; and (c) θ = π / 3 .
Materials 16 03628 g011
Figure 12. The peak forces for different angles θ (internal crack).
Figure 12. The peak forces for different angles θ (internal crack).
Materials 16 03628 g012
Figure 13. The crack pattern with shear load for different initial phason loads (internal crack): (a) w 0 = 8 × 10 4 mm ; (b) w 0 = 1.2 × 10 3 mm ; and (c) w 0 = 1.6 × 10 3 mm .
Figure 13. The crack pattern with shear load for different initial phason loads (internal crack): (a) w 0 = 8 × 10 4 mm ; (b) w 0 = 1.2 × 10 3 mm ; and (c) w 0 = 1.6 × 10 3 mm .
Materials 16 03628 g013
Figure 14. A rectangular QC with double internal cracks.
Figure 14. A rectangular QC with double internal cracks.
Materials 16 03628 g014
Figure 15. The peak forces for different angles θ (double cracks).
Figure 15. The peak forces for different angles θ (double cracks).
Materials 16 03628 g015
Figure 16. The crack pattern with Δ 2 = 32.5 mm for different angles (double cracks): (a) θ = π / 12 ; (b) θ = π / 4 ; (c) θ = 5 π / 12 .
Figure 16. The crack pattern with Δ 2 = 32.5 mm for different angles (double cracks): (a) θ = π / 12 ; (b) θ = π / 4 ; (c) θ = 5 π / 12 .
Materials 16 03628 g016
Figure 17. The crack pattern with θ = π / 4 for different distances Δ 2 : (a) Δ 2 = 32.5 mm ; (b) Δ 2 = 22.5 mm ; and (c) Δ 2 = 12.5 mm .
Figure 17. The crack pattern with θ = π / 4 for different distances Δ 2 : (a) Δ 2 = 32.5 mm ; (b) Δ 2 = 22.5 mm ; and (c) Δ 2 = 12.5 mm .
Materials 16 03628 g017
Table 1. Material properties of 2D decagonal QCs.
Table 1. Material properties of 2D decagonal QCs.
C 11 C 12 C 66 K 1 K 2 R 1 R 2
Al-Ni-Co (GPa)234.3057.3488.4512224−1.10.1
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, T.; Yang, Z.; Xu, C.; Xu, X.; Zhou, Z. A Phase Field Approach to Two-Dimensional Quasicrystals with Mixed Mode Cracks. Materials 2023, 16, 3628. https://doi.org/10.3390/ma16103628

AMA Style

Li T, Yang Z, Xu C, Xu X, Zhou Z. A Phase Field Approach to Two-Dimensional Quasicrystals with Mixed Mode Cracks. Materials. 2023; 16(10):3628. https://doi.org/10.3390/ma16103628

Chicago/Turabian Style

Li, Tong, Zhenting Yang, Chenghui Xu, Xinsheng Xu, and Zhenhuan Zhou. 2023. "A Phase Field Approach to Two-Dimensional Quasicrystals with Mixed Mode Cracks" Materials 16, no. 10: 3628. https://doi.org/10.3390/ma16103628

APA Style

Li, T., Yang, Z., Xu, C., Xu, X., & Zhou, Z. (2023). A Phase Field Approach to Two-Dimensional Quasicrystals with Mixed Mode Cracks. Materials, 16(10), 3628. https://doi.org/10.3390/ma16103628

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