Open access peer-reviewed chapter

# Probe on Rupture Theory of Soft-Matter Quasicrystals

Written By

Hui Cheng and Tian-You Fan

Submitted: 25 April 2019 Reviewed: 19 September 2019 Published: 23 September 2020

DOI: 10.5772/intechopen.89813

From the Edited Volume

## Fracture Mechanics Applications

Edited by Hayri Baytan Ozmen and H. Ersen Balcioglu

Chapter metrics overview

View Full Metrics

## Abstract

In this chapter, a probe on an important aspect, rupture theory of soft matter, is discussed, in which the soft matter and soft-matter quasicrystals are introduced at first. Then, we discuss the behaviour of the matter. For the soft-matter quasicrystals, there are basic equations describing their dynamics; due to the work of the second author of this chapter, this provides a fundamental for studying the rupture feature quantitatively. For general soft matter, there are no such equations so far, whose rupture theory quantitatively is not available at moment. The discussion on the soft-matter quasicrystals may provide a reference for other branches of soft matter.

### Keywords

• soft matter
• quasicrystals
• generalized dynamics
• equation of state
• crack
• stress intensity factor
• generalized Griffith criterion of rupture

## 1. Introduction

Soft-matter quasicrystals belong to a category of soft matter. It is well known that the soft matter is an intermediate phase between solid and liquid, which presents the behaviour of both solid and liquid, the first nature of which is fluidity, and behaves as a complex structure. Hence, the soft matter is named as a complex liquid or a structured liquid. Soft-matter quasicrystals are soft-matter with quasi-periodic symmetry. In this sense, they are a category of soft matter with highest symmetry so far. The high symmetry of the matter presents symmetry breaking and leads to importance of elementary excitations. This helps us to set up their dynamic equations and provide the fundamental for the rupture study.

## 2. Soft matter and the rupture problem

The liquid crystals, polymers, colloids, surfactants, and so on are in common named as soft matter. The 12-fold symmetry quasicrystals were observed most frequently in the soft matter [1, 2, 3, 4, 5, 6, 7]. At the same time, cracks in soft matter [8, 9] cannot be ignored and should be prevented [10] from the material safety. Other crack and rupture problems in soft matter can be referred from [11, 12, 13, 14, 15, 16, 17, 18]. This shows that the study on crack and rupture for soft matter including soft-matter quasicrystals is very important.

It is well known that the failure of brittle structural materials has been well studied. According to the work of Griffith [19], the existence and propagation of crack is the reason of the failure of these materials. Griffith used the exact solution of a crack in an infinite plate and calculated the crack strain energy. The energy is the function of crack size. He further calculated crack energy release rate and suggested that when the release rate equals to the surface tensile of the material, the solid will be a failure. His hypothesis was proved by experiments. This is the famous Griffith criterion. Afterward, the classical Griffith theory was developed by Irwin et al. for studying quasi-brittle failure, where the strain energy release rate was replaced by so-called stress intensity factor and the surface tensile was replaced by the fracture toughness of the material. This is the engineering approach of the Griffith theory, and played an important role in engineering application. The failure of ductile materials is also related with the existence and propagation of crack, but the problem has not been well studied due to the plastic deformation around the crack tip. The plastic deformation is a nonlinear irreducible process. The problem is extremely complex physically and mathematically. The failure of soft matter will be more complex than the ductile structural materials because of the existence of the fluidity. Especially, the experimental results are few of reported. As a most preliminary probe for studying the problem of soft matter, we try to draw from the idea of the classical Griffith work, i.e., to study the crack stability/instability, we can use the so-called following Griffith-Irwin criterion

K 1 a σ = K IC T f s E1

in which K 1 a σ represents elastic stress intensity factor, which a function of crack size a and applied external stress σ , can be determined by stress analysis of cracked materials, and the K IC T f s fracture toughness, a material constant but influenced by temperature and the structure factor f s , and the suffice I expresses mode I, i.e., the opening mode fracture (and the mode II is shearing mode, or slip mode fracture, and mode III is longitudinal shear mode, or tearing mode fracture, we here discuss only the opening mode). If the value of K 1 a σ is greater than that of K IC T f s , then the crack will propagate and the material will fracture.

Of course, the criterion (1) is only a reference for the soft matter, and a further analysis will be given in the subsequent sections.

## 3. Soft-matter quasicrystals and their generalized dynamics

There are a quite lot of references concerning the crack and rupture problems in soft matter; however, the quantitative analyses are not so much, because most branches of soft matter science are in qualitative stage so far. Either theoretical research or engineering application, the rupture problem of the soft matter needs a quantitative analysis.

Recently, generalized dynamics of soft-matter quasicrystals has been developed [20, 21, 22, 23, 24], which may become another quantitative branch in soft matter apart from the liquid crystals science. The generalized dynamics of soft-matter quasicrystals provides a tool for analyzing quantitatively rupture problem of the matter, whose result may be references of other categories of soft matter.

Soft-matter quasicrystals look like other categories of soft matter which belong to complex fluid; at mean time they present highly symmetry. For the high ordered phase, the symmetry breaking and elementary excitation principle are important. By using the Landau-Anderson [25, 26] symmetry breaking and elementary excitation principle, there are phonon and phason elementary excitations. As a class of soft matter, the fluidity is the substantive nature of the soft-matter quasicrystals; so Fan [20, 21, 22, 23, 24] introduced another elementary excitation—fluid phonon apart from phonons and phasons; of course, the concept of the fluid phonon is originated from the Landau School [27]. The introducing of fluid phonon requires a supplemented equation and the equation of state as well, which is also completed by Fan and co-worker [28]. With these bases, the generalized dynamics of soft-matter quasicrystals is set up.

For the need of the present chapter, we list the two-dimensional equations of the dynamics of soft-matter quasicrystals of 5- and 10-fold symmetry as follows:

ρ t + · ρ V = 0 ρ V x t + V x ρ V x x + V y ρ V x y = p x + η 2 V x + 1 3 η x · V + M 2 u x + L + M B x · u + R 1 2 w x x 2 + 2 2 w y x y 2 w x y 2 R 2 2 w y x 2 2 2 w x x y 2 w y y 2 A B 1 ρ 0 δρ x ρ V y t + V x ρ V y x + V y ρ V y y = p y + η 2 V y + 1 3 η y · V + M 2 u y + L + M B y · u + R 1 2 w y x 2 2 2 w x x y 2 w y y 2 + R 2 2 w x x 2 + 2 2 w y x y 2 w x y 2 A B 1 ρ 0 δρ y u x t + V x u x x + V y u x y = V x + Γ u [ M 2 u x + L + M x · u + R 1 2 w x x 2 + 2 2 w y x y w x y 2 R 2 2 w y x 2 2 2 w x x y 2 w y y 2 ] u y t + V x u y x + V y u y y = V y + Γ u [ M 2 u y + L + M y · u + R 1 2 w y x 2 2 2 w x x y 2 w y y 2 + R 2 2 w x x 2 + 2 2 w y x y 2 w x y 2 ] w x t + V x w x x + V y w x y = Γ w [ K 1 2 w x + R 1 2 u x x 2 2 2 u y x y 2 u x y 2 + R 2 2 u y x 2 + 2 2 u x x y 2 u y y 2 ] w y t + V x w y x + V y w y y = Γ w [ K 1 2 w y + R 1 2 u y x 2 + 2 2 u x x y 2 u y y 2 R 2 2 u x x 2 2 2 u y x y 2 u x y 2 ] p = 3 k B T l 3 ρ ρ 0 + ρ 2 ρ 0 2 + ρ 3 ρ 0 3 E2

in which u i denotes the phonon field, w i the phason field, V i the fluid velocity field, C ijkl is the phonon elastic constant tensor, K ijkl phason elastic constant tensor, and R ijkl and R klij are the phonon-phason coupling elastic constant tensors.

A and B are the constants describing density variation, g = ρ V , and η is the fluid viscosity, k B is the Boltzmann constant, T is the absolute temperature, and l 10 100 nm is the characteristic size of the soft matter, respectively.

## 4. An example of stress analysis of soft-matter quasicrystals

With Eq. (2), we can carry out a stress analysis of some fundamental specimens with crack of soft-matter quasicrystals; we here give only a computational example as shown in Figure 1.

If we want to obtain further information on deformation and motion of the material, we must solve the equations under appropriate initial and boundary conditions. To solve the problem, a specimen made by the matter should be optioned, which is subjected to some initial and boundary conditions. Here, the corresponding conditions of the specimen shown in Figure 1 are as follows:

t = 0 : V x = V y = 0 , u x = u y = 0 , w x = w y = 0 , p = p 0 ;
y = ± H , x < W : V x = V y = 0 , σ yy = σ 0 f t , σ yx = 0 , H yy = H yx = 0 , p = p 0 ; x = ± W , y < H : V x = V y = 0 , σ xx = σ xy = 0 , H xx = H xy = 0 , p = p 0 ; y = 0 , x < a : V x = V y = 0 , σ yy = σ yx = 0 , H yy = H yx = 0 , p = p 0 E3

In the present computation we take, 2 H = 0.01 m , 2 W = 0.01 m , 2 a = 0.0024 m , σ 0 = 0.01 MPa , ρ 0 = 1.5 × 10 3 kg / m 3 , η = 0.1 Pa s or 0.2 Pa s , ζ = 0 , L = 10 MPa , M = 4 MPa , K 1 = 0.5 L , R = R 1 = 0.04 M or 0.06 M , R 2 = 0 , Γ u = 4.8 × 10 17 m 3 s / kg , Γ w = 4.8 × 10 19 m 3 s / kg , A 0.2 MPa , B 0.2 MPa , and p 0 denotes 1 atm.

The initial and boundary value problem (3) of Eq. (2) can be solved by the finite difference method to solve the boundary value problem (Figure 2), e.g.,

2 u x x 2 = u x i + 1 j 2 u x i j + u x i 1 j h 2 , 2 u x y 2 = u x i j + 1 2 u x i j + u x i j 1 h 2 , u x t = u x k + 1 u x k τ .

so determine the phonon and phason displacement fields and fluid phonon velocity filed, then the phonon and phason strain tensors.

ε ij = 1 2 u i x j + u j x i , w ij = w i x j , E4

and the fluid phonon deformation rate tensor

ξ ̇ ij = 1 2 V i x j + V j x i E5

and according to constitutive laws

σ ij = C ijkl ε ik + R ijkl w kl , H ij = K ijkl w ij + R klij ε kl , p ij = p δ ij + σ ij ' = p δ ij + 2 η ξ ̇ kl , E6

we obtain the phonon stresses σ ij and fluid phonon stresses p ij , respectively, in which recall C ijkl the phonon elastic constant tensor, K ijkl phason elastic constant tensor, R ijkl and R klij the phonon-phason coupling elastic constant tensors, etc., refer to [24].

Due to the complexity of the equations, the computation is complex too. For the dynamic problems, i.e., there are manmade damping terms θρ , θρ V x , and so on, the iterative computation is easily stable, and for the static problems, we still take the iterative computation; however, the choosing of the manmade damping coefficient θ depends upon experience.

The computational results are mostly dependent on the ratio value of σ 0 / p 0 , i.e., the ratio of amplitudes of the phonon stress and fluid phonon stress apart from other factors.

## 5. Significance of fluid stress to crack initiation of growth and crack propagation

In the crack and fracture of brittle or quasibrittle failure of structural materials, the stress analysis is a basic task, from which one can determine the stress intensity factor, and then use the Griffith-Irwin criterion to analyze the crack stability/instability. For the Mode I crack, the tensile stress σ yy x y t is the most important, which leads to the crack surface opening and so the crack propagation.

After our computation, in soft matter including soft-matter quasicrystals, apart from the phonon stress σ yy x y t , there is the fluid phonon stress p yy x y t which is pressure and leads to crack closing (Figures 3 and 4).

In principle, the Griffith theory holds for describing brittle and quasi-brittle rupture of structural materials (or engineering materials). However, for soft matter including soft-matter quasicrystals, there is fundamental difference with the structural materials. The key reason about this is the existence of fluid effect; in terminology of soft-matter quasicrystal study, it is also called the existence of fluid phonon. According to our analysis, the effect of fluid stress intensity factor is opposite to the elastic stress intensity factor (or by using the terminology of soft-matter quasicrystals study, the elastic stress intensity factor is also called the phonon stress intensity factor). K 1 phonon denotes the phonon stress intensity factor and K 1 fluid denotes the fluid phonon stress intensity factor which can be defined as follows:

K 1 phonon = lim x a 0 + π x a 0 σ yy x 0 t E7
K 1 fluid = lim x a 0 + π x a 0 p yy x 0 t E8

The computational results are dependent upon the ratio σ 0 / p 0 , while the fluid stress intensity factor is expressed by K 1 fluid , which is also a function of time in general, because K 1 phonon > 0 leads to crack surface opening, and according to our computation, K 1 fluid < 0 leads to crack surface closing. The total stress intensity factor is denoted by

K 1 total = K 1 phonon + K 1 fluid . E9

Therefore, the value of K 1 total is smaller than K 1 phonon .

K 1 is the most often used in engineering design parameter in fracture mechanics and, hence, must be understood if we are to design fracture tolerant materials used in bridges, buildings, aircraft, or even bells. Polishing cannot detect a crack. Typically, if a crack can be seen, it is very close to the critical stress state predicted by the stress intensity factor. The magnitude of K ˜ 1 depends on sample geometry, the size and location of the crack, and the magnitude and the modal distribution of loads on the material.

The normalized elastic dynamic stress intensity factor is defined as

K ˜ 1 t = K 1 t / π a 0 σ 0 E10

## 6. Possible criterion of rupture

For soft matter, there is fracture criterion:

K 1 total = K IC T . E11

The Griffith theory can be extended to soft matter. If K 1 total t > K IC T then crack surface to be opening, and crack propagate leads to fracture; otherwise K 1 total t < K IC T then crack surface to be closing, and crack cannot propagate, the material is in safe. The conclusion is that the physical state around crack tip in soft matter quasicrystals is dependent on the competition of phonon stress and fluid stress.

## 7. The references for other soft matter

The example given in previous section is only the most preliminary discussion on the probe of rupture of soft-matter quasicrystals. It does not mean any in-depth and exact theory of the problem. The significance of the discussion lies in providing a possible drawing for other categories of soft matter in the quantitative study of rupture problems. The formulation of the present study may be referenced by other branches of soft matter, in which the constitutive laws for themselves are needed, and phason field can be excluded for non-quasicrystalline materials, and so on. With these considerations, the rupture theory for the other categories of soft matter out of the soft-matter quasicrystals may be set up.

## 8. Conclusion and discussion

From Figure 5, we also can see that the real value of fluid DSIF is negative; thus, we can say it causes the crack closed. When comparing (a) with (b), the absolute value of fluid DSIF is larger than the elastic ones; thus, the crack cannot propagate. Then, we can conclude that in the quasicrystal of soft matter, the fluid effect is important, and we should do the study under the hydrodynamics theory.

The equation set of generalized dynamics of soft-matter quasicrystals is complete; it is possible to provide the hydrodynamics study on the mass distribution, deformation, and motion, including rupture analysis of the new phase. In particular, the dynamics stress intensity factor (DSIF) was evaluated, and the computation presents highly stability. In the results, there are two regimes: one is elastic DSIF and the other is fluid DSIF. The former leads to crack surface opening and the latter leads to crack surface closing. We can conclude that whether or not the crack open or close depends on the competition between the elastic DSIF and fluid DSIF. Of course, the theoretical analysis is anxious for the experimental verification currently, which is our attempt work in future.

The above introducing is rough and preliminary, in fact, around the crack tip, due to the stress concentration and high stress grad, nonelastic stress appears, the Griffith theory will not be effective, i.e., the theory based on the stress intensity factor will not be valid, one must carry out nonlinear analysis, e.g. the theory based on the crack tip opening displacement, in which the size of the plastic zone near the crack tip is also considered, refer to Fan [24], which may help us to overcome the difficulty mentioned above. However this is a very difficult question which has not been solved so far.

## Acknowledgments

This work is supported by the National Natural Science Foundation of China through grant 11272053.

## Conflict of interest

The authors declare no conflict of interest.

## A.1 The final governing equations of generalized dynamics of soft-matter quasicrystals with 10-fold symmetry in three dimension

In the text, Eq. (2) is the planar field form only. In book [24], it provides an explicit form of the equations and lists below soft-matter quasicrystals with 10-fold symmetry:

ρ t + · ρ V = 0 ρ V x t + V x ρ V x x + V y ρ V x y + V z ρ V x z = p x + η 2 V x + 1 3 η x · V + C 66 2 y 2 + C 44 2 z 2 u x + C 12 + C 66 2 u y x y + C 13 + C 44 C 11 2 u z x z + C 11 B x · u + R x 1 · w R y w x y w y x A B 1 ρ 0 δρ x ρ V y t + V x ρ V y x + V y ρ V y y + V z ρ V y z = p y + η 2 V y + 1 3 η y · V + C 12 + C 66 2 u x x y + C 66 2 x 2 + C 11 2 y 2 + C 44 2 z 2 u y + C 13 + C 44 2 u z y z + C 11 B y · u R x w x y w y x R y 1 · w A B 1 ρ 0 δρ y ρ V z t + V x ρ V z x + V y ρ V z y + V z ρ V z z = p z + η 2 V z + 1 3 η z · V + C 44 2 x 2 + C 44 2 y 2 + C 33 C 13 C 44 2 z 2 u z + C 13 + C 44 B z · u A B 1 ρ 0 δρ z u x t + V x u x x + V y u x y + V z u x z = V x + Γ u [ C 11 2 x 2 + C 66 2 y 2 + C 44 2 z 2 u x + C 12 + C 66 2 u y x y + C 13 + C 44 2 u z x z + R x 1 · w R y w x y w y x ] u y t + V x u y x + V y u y y + V z u y z = V y + Γ u [ C 12 + C 66 2 u x x y + C 66 2 x 2 + C 11 2 y 2 + C 44 2 z 2 u y + C 13 + C 44 2 u z y z R x w x y w y x R y 1 · w ] u z t + V x u z x + V y u z y + V z u z z = V z + Γ u [ C 13 + C 44 2 u x x z + 2 u y y z + C 44 2 x 2 + C 44 2 y 2 + C 33 2 z 2 u z ] w x t + V x w x x + V y w x y + V z w x z = Γ w K 1 1 2 w x + K 4 2 w x z 2 + K 2 y x y w y + R x u x x u y y R y u x y + u y x w y t + V x w y x + V y w y y V z w y z = Γ w K 1 1 2 w y + K 4 2 w y z 2 + R x u x y + u y x + R y u x x u y y p = f ρ = 3 k B T l 3 ρ 0 3 ρ 0 2 ρ + ρ 0 ρ 2 + ρ 3 EA1

in which 2 = 2 x 2 + 2 y 2 + 2 z 2 , 1 2 = 2 x 2 + 2 y 2 , = i x + j y + k z , V = i V x + j V y + k V z , u = i u x + j u y + k u z , and C 11 , C 12 , C 13 , C 33 , C 44 , C 66 = C 11 C 12 / 2 are the phonon elastic constants, K 1 , K 2 , K 3 , K 4 are the phason elastic constants, R is the phonon-phason coupling constant, η is the fluid dynamic viscosity, Γ u and Γ w the phonon and phason dissipation coefficients, and A and B the material constants due to variation of mass density, respectively.

Equation (A1) is the final governing equations of dynamics of soft-matter quasicrystals of 10-fold symmetry in three-dimensional case with field variables u x , u y , u z , w x , w y , V x , V y , V z , ρ and p ; the amount of the field variables is 10, and the amount of field equations is 10 too; among them are: (A1a) (the first of (A1)) is the mass conservation equation, (A1b)–(A1d) (the second to fourth of (A1)) are the momentum conservation equations or generalized Navier-Stokes equations, (A1e)–(A1g) (the fifth to seventh of (A1)) are the equations of motion of phonons due to the symmetry breaking, (A1h) and (A1i) (the eighth to ninth of (A1)) are the phason dissipation equations, and (A1j) (the tenth of (A1)) is the equation of state, respectively. The equations are consistent to be mathematical solvability, if there is lack of the equation of state, the equation system is not closed, and has no meaning mathematically and physically. This shows that the equation of state is necessary.

The equation set (A1) is the three-dimensional form of the equations, and in our solution, we computed only the two-dimensional form, i.e., the plane field form; in this special case, the 5- and 10-fold symmetry quasicrystals have the same governing equations, but in the three-dimensional case, the equation set (A1) is valid only for the 10-fold symmetry quasicrystals.

## References

1. 1. Zeng X, Ungar G, Liu Y, Percec V, Dulcey AE, Hobbs JK. Supramolecular dendritic liquid quasicrystals. Nature. 2004;428(6979):157-160. DOI: 10.1038/nature02368
2. 2. Takano A, Kawashima W, Noro A, Isono Y, Tanaka N, Dotera T, et al. A mesoscopic Archimedean tiling having a new complexity in an ABC star polymer. Journal of Polymer Science Part B: Polymer Physics. 2005;43(18):2427-2432. DOI: 10.1002/polb.20537
3. 3. Hayashida K, Dotera T, Takano A, Matsushita Y. Polymeric quasicrystal: Mesoscopic quasicrystalline tiling in ABC star polymers. Physical Review Letters. 2007;98(19):195502. DOI: 10.1103/PhysRevLett.98.195502
4. 4. Talapin DV, Shevchenko EV, Bodnarchuk MI, Ye X, Chen J, Murray CB. Quasicrystalline order in self-assembled binary nanoparticle superlattices. Nature. 2009;461(7266):964-967. DOI: 10.1038/nature08439
5. 5. Fischer S, Exner A, Zielske K, Perlich J, Deloudi S, Steurer W, et al. Colloidal quasicrystals with 12-fold and 18-fold diffraction symmetry. Proceedings of the National Academy of Sciences. 2011;108(5):1810-1814. DOI: 10.1073/pnas.1008695108
6. 6. Cheng SZD. Giant surfactants based on precisely functionalized POSS nano-atoms: Tuning from crystals to Frank-Kasper phases and quasicrystals. In: First Annual Symposium on Frontiers of Soft Matter Science and Engineering, Beijing; 2015. Bulletin of the American Physical Society, 2016
7. 7. Yue K, Huang MJ, Marson R, He JL, Huang JH, Zhou Z, et al. Geometry induced sequence of nanoscale Frank-Kasper and quasicrystal mesophases in giant surfactants. Proceedings of the National Academy of Sciences. 2016;113:1392-1400
8. 8. Goehring L, Conroy R, Akhter A, et al. Evolution of mud-crack patterns during repeated drying cycles. Soft Matter. 2010;6(15):3562-3567. DOI: 10.1039/b922206e
9. 9. Carbone G, Pierro E, Gorb SN. Origin of the superior adhesive performance of mushroom-shaped microstructured surfaces. Soft Matter. 2011;7(12):5545-5552. DOI: 10.1039/c0sm01482f
10. 10. Xuan Y, Guo X, Cui Y, et al. Crack-free controlled wrinkling of a bilayer film with a gradient interface. Soft Matter. 2012;8(37):9603-9609. DOI: 10.1039/c2sm25487e
11. 11. Haque MA, Kurokawa T, Kamita G, et al. Lamellar bilayers as reversible sacrificial bonds to toughen hydrogel: Hysteresis, self-recovery, fatigue resistance, and crack blunting. Macromolecules. 2011;44(22):8916-8924. DOI: 10.1021/ma201653t
12. 12. Hatton B, Mishchenko L, Davis S, et al. Assembly of large-area, highly ordered, crack-free inverse opal films. Proceedings of the National Academy of Sciences. 2010;107(23):10354-10359. DOI: 10.1073/pnas.1000954107
13. 13. Lazarus V, Pauchard L. From craquelures to spiral crack patterns: Influence of layer thickness on the crack patterns induced by desiccation. Soft Matter. 2011;7(6):2552-2559. DOI: 10.1039/c0sm00900h
14. 14. Tirumkudulu MS, Russel WB. Cracking in drying latex films. Langmuir. 2005;21(11):4938-4948. DOI: 10.1021/la048298k
15. 15. van der Kooij HM, Sprakel J. Watching paint dry; more exciting than it seems. Soft Matter. 2015;11(32):6353-6359. DOI: 10.1039/C5SM01505G
16. 16. Martinez CJ, Lewis JA. Shape evolution and stress development during latex-silica film formation. Langmuir. 2002;18(12):4689-4698. DOI: 10.1021/la0114833
17. 17. Tirumkudulu MS, Russel WB. Role of capillary stresses in film formation. Langmuir. 2004;20(7):2947-2961. DOI: 10.1021/la0356250
18. 18. Bohn S, Pauchard L, Couder Y. Hierarchical crack pattern as formed by successive domain divisions. Physical Review E. 2005;71(4):046214. DOI: 10.1103/PhysRevE.71.046215
19. 19. Griffith AA. The phenomena of rupture and flow in solids. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences. 1921;221:582-593. DOI: 10.2307/91192
20. 20. Fan TY. Mathematical Theory of Elasticity of Quasicrystals and Its Applications. 1st ed./2nd ed. Beijing/Heidelberg: Science Press/Springer-Verlag; 2010/2016. DOI: 10.1007/978-3-642-14643-5
21. 21. Fan TY. Poisson bracket method and its applications to quasicrystals, liquid crystals and a class of soft matter. Acta Mechanica Sinica. 2013;45(4):548-559. DOI: 10.6052/0459-1879-12-346 (in Chinese)
22. 22. Fan TY. Equation system of generalized hydrodynamics of soft-matter quasicrystals. Applied Mathematics and Mechanics. 2016;37(4):331-344 (in Chinese)
23. 23. Fan TY. Generalized hydrodynamics for second 2D soft-matter quasicrystals. Applied Mathematics and Mechanics. 2017;38(2):189-199 (in Chinese)
24. 24. Fan TY. Generalized Dynamics of Soft-Matter Quasicrystals—Mathematical Models and Solutions. Beijing/Heidelberg: Beijing Institute of Technology Press/Springer-Verlag; 2017
25. 25. Landau LD, Lifshitz EM. Theoretical Physics V, Statistical Physics, Part I. 3rd ed. New York: Pergamon Press; 1980
26. 26. Anderson PW. Basic Notations of Condensed Matter Physics. Menlo Park: Benjamin Cummings; 1984
27. 27. Lifshitz EM, Pitaevskii LP. Theoretical Physics V, Statistical Physics, Part II. New York: Pergamon Press; 1980
28. 28. Fan L, Fan TY. Equation of state of structured liquid. In: First Annual Symposium on Frontiers of Soft Matter Science and Engineering; December 12, 2015; Beijing, China; 2015

Written By

Hui Cheng and Tian-You Fan

Submitted: 25 April 2019 Reviewed: 19 September 2019 Published: 23 September 2020