Research Article  Open Access
Fracture of Parallel Strand Bamboo Composite under Mode I Loading: DCB Test Investigation
Abstract
This paper describes the experimental studies on Mode I fracture of parallel strand bamboo (PSB) by the double cantilever beam (DCB) test. Rcurves based on the elementary beam theory and specimen compliance are proposed in order to overcome the difficulties to monitor the crack propagation during experiments. The results demonstrate that the energy release rate (ERR) is influenced by the specimen geometry, i.e., the specimen width and initial crack length. The ERR at the plateau level is similar for the range of the analyzed widths (B = 20, 40, and 60 mm), while it decreases with width increasing up to 80 mm and 100 mm. The energy release rate for PSB specimens would verge to a stable value with the width increasing up to a specific value, while the value of the energy release rate will be influenced by the initial crack length. Consequently, the DCB tests also show that the obtained Rcurve in this study is not a material property.
1. Introduction
Parallel strand bamboo (PSB) is manufactured by parallelly gluing bamboo strands together under controlled temperature and pressure. Because bamboo strips are parallelly glued along the longitudinal direction and uniformly distributed over the transverse direction, PSB can be considered as a unidirectional and orthotropic fibrous composite, as shown in Figure 1 [1]. More recently, PSB wins growing interests as an alternative of wood composites for construction engineering in China due to its fastgrowing feature and excellent structural performance. Aimed at the structural use, the failure modes of combined compression and bending PSB members were studied by Huang et al. [2–4]. It was found that the fracture along the fiber interfaces was one of the major failure modes of PSB composites. Once a PSB member subjected to an increasing external load, the initial deflects, such as microvoids and microcracks, could be advancing or growing into macrocracks and consequently result in catastrophic failure of structures. Therefore, the fracture of PSB composite is one of the major concerns in PSB structural design.
(a)
(b)
The objective of the present study is to investigate the Mode I fracture properties and to measure the fracture toughness of the PSB composite through DCB experiments, which is one of the necessities for establishing design allowable values used in damage tolerance analysis of PSB structures.
The double cantilever beam (DCB) test [5–7], a standard test procedure prescribed in ASTM D552813 [8], was employed as the test procedure in this study. The test uses a rectangular specimen with different widths, different initial crack lengths, and constant thickness. It contains a preimplanted, nonadhesive insert or crack as an initial delamination. The opening load is applied perpendicular to the initial crack surface to induce the Mode I crack. Such a test analysis is based on the beam theory, and the fracture toughness can be measured through the energy release rate (ERR) for orthotropic fibrous composites. The ERR can be obtained by various methods, including direct area integration of loadingunloading curves or by means of compliance measurements [9]. Since the data reduction of DCB test is easier than the other test approaches [10, 11], it is now the most popular method used to determine Mode I fracture toughness.
Recommended by ASTM and ISO standards [9, 12, 13], the energy release rate can be estimated through the change of compliance. Theoretically, it can be obtained by differentiating the specimen compliance, C, with respect to the crack length a, i.e., [14–18], where P is the applied load and B the width of the specimen. In the frame of linear fracture mechanics, the root condition is supposed to be fully builtin; hence, the compliance equation can be given by simple beam theory, i.e., [9]. Nevertheless, because of the fuzzy boundary conditions at the end of opening arms, the compliance equation C(a) cannot be determined without controversy. On the other hand, the crack length cannot be measured either with desired precision. Therefore, various data reduction methods were developed to overcome this disadvantage, such as the area method, in which the compliance equation was built through data fitting for tests [19, 20]. Hashemi et al. [21] compared the different data reduction methods for obtaining the energy release rate of fibrous composites. They proposed that the correction of crack length was necessary in the beam theorybased method because the end of opening arms was not perfectly clamped. They developed the corrected beam theory (CBT), providing a correction in the crack length based on a compliance calibration. However, the fracture mechanism of fibrous composite is much complicated than that of brittle materials. The major aspect is that the microvoids coalescing, fine cracks extending, and fiber bridging lead to a large fracture process zone (FPZ) in front of the crack tip, which makes the clear location of the crack tip impracticable [22–25]. Furthermore, the development of the large FPZ in the cracktip front delays the fracture of specimen, which consequently makes the builtin assumption no longer valid. It has been well recognized that the use of elementary beam theory with the builtin assumption at the crack tip has a significant error for calculating the energy release rate of fibrous composites [26]. To date, there is no method available to exactly identify the crack tip for the DCB test of fibrous composites. For this reason, extensive researches have been carried out to avoid locating the crack tip [27–32]. Among these researches, the concept of equivalent linear elastic fracture mechanics (LEFM) was widely accepted to deal with the fracture with a large FPZ [31]. According to the concept of equivalent LEFM, the increasing of compliance, owing to the development of FPZ or main crack propagation with FPZ, is attributed to the elastic crack which equals to the actual crack with its FPZ [33]. Therefore, the crack length can be estimated according to the associated compliance obtained by experiments. Furthermore, the Rcurve can also be measured using only a monotonic loaddisplacement record.
In the present study, the data reduction method was based on the theory of equivalent LEFM. Compliance calibration was introduced to determine the crack length. The direct integral approach was adopted to calculate the energy release rate.
2. Data Reduction Scheme
Figure 2 illustrates the principle of evaluating energy release rate based on the theory of equivalent LEFM. In the schematic loaddisplacement curve , two consecutive points and corresponding to two different crack lengths and , respectively, are selected. Suppose that the crack extends from length to with a small amount of propagation increment and that the loaddisplacement trajectory goes from M_{1} to M_{2}; thus, the shadow area circulated by the dashed lines OM_{1}, and OM_{2} and the segment of loaddisplacement M_{1}M_{2} must be equal to the energy released as the crack growth of the extension [19]. Therefore, the energy release rate can be calculated bywhere and stand for the applied load and associated displacement at loading position, respectively. The crack length can be determined from the test compliance . According to the beam theory, the compliance can be calculated by [34]where E_{L} is the longitudinal elastic modulus, G_{LT} is the shear modulus in LT plane, i.e., the 12 or 13 plane as shown in Figure 1(b), I is the inertia moment of the cross section, and A is the area of cross section. Once the compliance is obtained by experiment and the crack length can be determined by resolving equation (2) or by numerical iteration in terms of the following equation:where is the calibration parameter to eliminate the system error of the experiment, in which is the test compliance corresponding to the initial crack length a_{0} and stands for the theoretical compliance corresponding to the initial crack length calculated by equation (2).
3. Experimental Investigation
3.1. Materials
The test PSB was provided by Feiyu Bamboo Products, Jianghxi, China. The material density was 1.26 g/cm^{3}, and the moisture content was 11%. Mechanical properties in paralleltograin direction, which are involved in test analysis, were pretested following the method recommended in ASTM D14314 [35], and the results are collected in Table 1, where the subscripts L and T are the paralleltograin direction and perpendiculartograin direction, respectively, and represents Poisson’s ratio of PSB.

3.2. Specimen Preparation
The geometry of DCB specimen is shown in Figure 3. The dimensions of DCB specimen were determined referring to the standard test procedure addressed in ASTM D552813 [8]. In this method, the DCB dimension was designed to ensure the damage zone or nonlinear deformation developed along the delamination front and the stable crack growth can be achieved.
Totally four groups of specimens with different initial crack lengths and different widths were prepared, as illustrated in Table 2, in which each group consists of the same initial crack length and 5 types of specimens with different widths. A label system was designed to identify the specimen. In this system, the latter A, B, C, and D correspond to the initial crack length of 34 mm, 67 mm, 100 mm, and 134 mm, respectively. The number prior to the letter of the group name stands for the width of the specimen, as shown in Table 2. The initial crack was firstly introduced by a 1.5 mm thick saw kerf; afterwards, precrack with length about 1 mm was extended by using a cutting blade [10]. At the end of the DCB test, two bolt holes of 8 mm in diameter were drilled for the sake of joining DCB specimen to the actuator of the test machine, as illustrated in Figure 3.

3.3. Test Procedure
The fracture tests were performed on a servohydraulic universal test machine of 20 kN capacity in room ambient circumstance. The test setup is illustrated in Figure 4(a). The specimen was joined to the load actuator through two steel rods of 7.5 mm diameter, as shown in Figure 4(b). Loading was controlled by the displacement of a moveable actuator at the speed of 1.0 mm/min. A microscope digital camera was mounted in front of the specimen to monitor the crack propagation and take the images of the crack tip. The applied load and the displacement at loading position were simultaneously recorded at a frequency of 10 Hz. The opening displacement at the initial crack tip was measured by using a clipon gauge (COD gauge) symmetrically fixed at the two sides of the crack tip through embedded aluminum flakes (Figure 4(b)). The crack length during propagation can be observed by using the microscopic camera, as shown in Figure 4(b).
(a)
(b)
4. Test Results
4.1. LoadDisplacement Curve
Figure 5 illustrates a typical loaddisplacement curve. Roughly three stages can be observed. The first one (stage I) is nodamage stage from origin to the proportional limit point where the curve deviates from its original direction. The loaddisplacement curve exhibits perfect linearity in this stage. Once the load exceeds the point of proportional limit (stage II), damage onset extends and coalesces to form the fracture process zone (FPZ) and consequently leads to the compliance augment and the curve continuously deviates from its original direction as the load increases. Hence, this stage can be understood as the stage of FPZ development. When the FPZ fully developed, the loaddisplacement curve would reach its critical point, where the load reached its maximum value and the crack propagation began. Hence, the third stage (stage III) is the crack propagation stage.
Figure 6 presents the loaddisplacement profiles of the specimens. For the specimens of Group A, the loaddisplacement curve sharply declines once the load reaches the maximum value, which indicates that the crack propagation is unstable. For the specimens in Groups B, C, and D with longer initial crack length, the loaddisplacement profile slowly declines in postsummit segments. Furthermore, the loaddisplacement profile declined more slowly with longer initial crack length, which indicated that the stability of crack propagation is sensitive to the initial crack length. Therefore, enough initial crack length is needed for obtaining stable crack propagation. In the present study, the ratios of initial crack length to ligament length of Groups B, C, and D are 0.191, 0.286, and 0.383, respectively, and their crack propagates in a stable and ductile manner, which means specimens with the ratio of Groups B, C, and D can be used to determine the fracture properties.
(a)
(b)
(c)
(d)
The second factor is the width scaling which has a significant influence on fiber bridging in the wake of the crack, consequently affecting the delamination behavior. As shown in Figure 7, the loaddisplacement curves for partial specimens are presented with the same initial crack length but different widths. It can be concluded that the crack grows in a “thumbnail” shape [36] under constant or decreasing load, while the overall displacement is increased. The “thumbnail” theory means that the crack extends mainly from the center of the thickness of specimens, while the edge regions are plastically deformed. Therefore, with the width increasing, the energy release rate of DCB specimens for PSB composite would converge to a stable value.
(a)
(b)
It is worth mentioning that only partial specimens exhibit their loaddisplacement curves or Rcurves in later analysis due to the similar tendency of these curves. Moreover, because of the manufacturing defects for some PSB specimens, the data cannot be used to analyze the energy release rate in the DCB test. Therefore, some specimens present two lines in partial loaddisplacement curves or Rcurves.
4.2. Fracture Process Zone
Damage begins when the external load exceeds its proportional limit, which is characterized by local microcracks at the crack tip front zone. This damaged zone is comprised of microcracks between the grains or through the grains. The microcracks are consequently extended and coalesced to form macrocracks and local microcracks at the advancing crack tip front as the external load increases. This fracture process and damaged zone can be called fracture process zone (FPZ). Before the load reaches its maximum value, the crack interfaces are bonded by fiber bridges which either rupture or peel off from the crack surfaces. This fracture mechanism leads to a large FPZ which can be observed at the crack tip front before crack propagation, as shown in Figure 8. Due to the restriction of fiber bridges, the strain energy is smoothly consumed through the fibrous rupture or pulloff, which leads to stable crack propagation and a rising Rcurve.
Several micrographs were taken from the tested specimens at the fracture surfaces using the scanning electron microscope (SEM), as shown in Figure 9. Figure 9(a) exhibits the typical fracture surface with a limited amount of fiber pullout in the debonding plies, which is presented with higher magnitude in Figure 9(b). The pullout process could have created a fibrebridged zone in the wake of the advancing crack tip. During the whole experiment, the growth and eventual stabilization of this fibrebridged zone could account for the tendency observed in later Rcurves.
4.3. RCurve Measurement
Energy release rate for each test specimen was measured by the method addressed in Section 2. The Rcurves of partial DCB specimens with different widths and different initial crack lengths corresponding to specific thickness are presented in Figures 10 and 11. It was observed that the energy release rate increases up to a steadystate toughness after the initiation of delamination. As illustrated in Section 4.2, the fiber bridging effect during the whole crack propagation led to eventual stabilization of crack growths for all tested specimens. Meanwhile, the results also implied that the Rcurve obtained in this study was not a material property, i.e., it depended on specimen geometry. Consequently, when the crack starts to propagate, the bridging zone is created. By developing the bridging zone length, the strain energy release rate increases up to the steadystate fracture toughness. The length between the initial crack length and crack length corresponding to the steadystate fracture toughness is called as the steadystate bridging zone length. From Figure 10, the value of the energy release rate was always affected by the initial crack length. The major reason may be that the steadystate bridging zone of PSB specimens was always influenced by the preimbedded crack.
(a)
(b)
(a)
(b)
From Figure 11, it can also be concluded that the energy release rate of DCB specimens for the PSB composite would verge to a stable value with the width increasing up to a specific value, which will be identified in the future DCB tests for PSB specimens. As illustrated in the abovementioned “thumbnail” theory, the width scaling has a significant effect on fiber bridging in the wake of the crack, which make the crack grow in a “thumbnail” shape, i.e., the crack extends mainly from the center of specimens, while the edge regions are plastically deformed. Therefore, the energy release rate of DCB specimens for the PSB composite would converge to a stable value with the increasing width. Table 3 gives the average fracture toughness for partial test samples.

5. Conclusion
This paper has investigated an experimental procedure to obtain the energy release rate (ERR) using DCB tests for parallel strand bamboo (PSB) composites. The effects of specimen width and initial crack length on largescale fiber bridging in Mode I fracture of unidirectional PSB composites were investigated. The Rcurve measurements, for the range of the analyzed widths (B = 20, 40, 60, 80, and 100 mm) and initial crack lengths ( = 34, 67, 100, and 134 mm), show that the energy release rate is influenced by both the width and initial crack length of specimens, in that the ERR at the plateau level decreased with an increasing width and initial crack length. Consequently, the DCB tests imply that the Rcurve obtained in the present study is not a material property because it depends on the specimen geometry. However, the ERR at the plateau level is decreased to similar values as the width increases up to 80 mm and 100 mm. It may be concluded that the energy release rate would verge to a stable value with the width increasing up to more than a specific value, which can be clarified in the future DCB tests for PSB composites.
Nomenclature
:  Specimen width 
:  Specimen compliance 
:  Crack length 
:  Applied load 
:  Energy release rate 
:  Young’s modulus 
:  Height of the cantilever portion, H/2 
:  Specimen thickness 
:  Loadingline displacement 
:  Crack extension 
:  Elastic modulus in longitudinal direction 
:  Shear modulus in LT plane 
:  Inertia moment of cross section 
:  Area of cross section 
:  Calibration parameter to eliminate errors from experiments 
:  Theoretical compliance corresponding to initial crack length 
:  Test compliance corresponding to the initial crack length 
:  Initial crack length 
:  Poisson’s ratio 
:  Specimen length. 
Data Availability
The data of the DCB test for PSB composites used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The research was supported by the National Natural Science Foundation of China (no. 51778299) and the Priority Academic Program Development of Jiangsu Higher Education Institutions.
References
 D. Huang, Y. Bian, A. Zhou, and B. Sheng, “Experimental study on stressstrain relationships and failure mechanisms of parallel strand bamboo made from phyllostachys,” Construction and Building Materials, vol. 77, pp. 130–138, 2015. View at: Publisher Site  Google Scholar
 D. S. Huang, A. P. Zhou, H. T. Li, Y. Su, and G. Chen, “Experimental study on the tensile properties of bamboo related to its distribution of vascular bundles,” Key Engineering Materials, vol. 517, pp. 112–117, 2012. View at: Publisher Site  Google Scholar
 D. Huang, A. Zhou, and Y. Bian, “Experimental and analytical study on the nonlinear bending of parallel strand bamboo beams,” Construction and Building Materials, vol. 44, pp. 585–592, 2013. View at: Publisher Site  Google Scholar
 D. Huang, Y. Bian, D. Huang, A. Zhou, and B. Sheng, “An ultimatestatebasedmodel for inelastic analysis of intermediate slenderness PSB columns under eccentrically compressive load,” Construction and Building Materials, vol. 94, pp. 306–314, 2015. View at: Publisher Site  Google Scholar
 M. F. Kanninen, “An augmented double cantilever beam model for studying crack propagation and arrest,” International Journal of Fracture, vol. 9, pp. 83–92, 1973. View at: Google Scholar
 R. Olsson, “A simplified improved beam analysis of the DCB specimen,” Composites Science and Technology, vol. 43, no. 4, pp. 329–338, 1992. View at: Publisher Site  Google Scholar
 F. Ozdil and L. A. Carlsson, “Beam analysis of angleply laminate DCB specimens,” Composites Science & Technology, vol. 59, no. 2, pp. 305–315, 1999. View at: Publisher Site  Google Scholar
 American Society for Test Materials (ASTM), ASTM D552813, Standard Test Method for Mode I Interlaminar Fracture Toughness of Unidirectional FiberReinforced Polymer Matrix Composite, ASTM, West Conshohocken, PA, USA, 2013.
 H. Yoshihara and T. Kawamura, “Mode I fracture toughness estimation of wood by DCB test,” Composites Part A: Applied Science and Manufacturing, vol. 37, no. 11, pp. 2105–2113, 2006. View at: Publisher Site  Google Scholar
 M. F. S. F. de Moura, J. J. L. Morais, and N. Dourado, “A new data reduction scheme for mode I wood fracture characterization using the double cantilever beam test,” Engineering Fracture Mechanics, vol. 75, no. 13, pp. 3852–3865, 2008. View at: Publisher Site  Google Scholar
 J. De Gracia, A. Boyano, A. Arrese, and F. Mujika, “A new approach for determining the Rcurve in DCB tests without optical measurements,” Engineering Fracture Mechanics, vol. 135, pp. 274–285, 2015. View at: Publisher Site  Google Scholar
 American Society for Test Materials (ASTM), ASTM D6671M19, Standard Test Method for Mixed Mode IMode II Interlaminar Fracture Toughness of Unidirectional Fiber Reinforced Polymer Matrix Composites, ASTM, West Conshohocken, PA, USA, 2019.
 International Organization for Standardization (ISO), ISO 150242001, Fiberreinforced Composite  Determination of Mode I Interlaminar Fracture Toughness, GIc, for Unidirectional Reinforced Materials, ISO, Geneva, Switzerland, 2001.
 W. Xu and Z. Z. Guo, “A simple method for determining the mode I interlaminar fracture toughness of composite without measuring the growing crack length,” Engineering Fracture Mechanics, vol. 191, pp. 476–485, 2018. View at: Publisher Site  Google Scholar
 D. Huang, B. Sheng, Y. Shen, and Y.H. Chui, “An analytical solution for double cantilever beam based on elasticplastic bilinear cohesive law: analysis for mode I fracture of fibrous composites,” Engineering Fracture Mechanics, vol. 193, pp. 66–76, 2018. View at: Publisher Site  Google Scholar
 S. Hofmann, “Mode I delamination onset in carbon fibre reinforced SiC: double cantilever beam testing and cohesive zone modelling,” Engineering Fracture Mechanics, vol. 182, pp. 506–520, 2017. View at: Publisher Site  Google Scholar
 E. FarmandAshtiani, J. Cugnoni, and J. Botsis, “Specimen thickness dependence of large scale fiber bridging in mode I interlaminar fracture of carbon epoxy composite,” International Journal of Solids and Structures, vol. 55, pp. 58–65, 2015. View at: Publisher Site  Google Scholar
 A. Jyoti, R. F. Gibson, and G. M. Newaz, “Experimental studies of Mode I energy release rate in adhesively bonded width tapered composite DCB specimens,” Composites Science and Technology, vol. 65, no. 1, pp. 9–18, 2005. View at: Publisher Site  Google Scholar
 N. A. Plan, S. Morel, and M. Chaplain, “Mixedmode fracture in a quasibrittle material: Rcurve and fracture criterion—application to wood,” Engineering Fracture Mechanics, vol. 156, pp. 96–113, 2016. View at: Publisher Site  Google Scholar
 S. Bennati and P. S. Valvo, “An experimental compliance calibration strategy for mixedmode bending tests,” Procedia Materials Science, vol. 3, pp. 1988–1993, 2014. View at: Publisher Site  Google Scholar
 S. Hashemi, A. J. Kinloch, and J. G. Williams, “Corrections needed in doublecantilever beam tests for assessing the interlaminar failure of fibrouscomposites,” Journal of Material Science Letters, vol. 8, no. 2, pp. 125–129, 1989. View at: Publisher Site  Google Scholar
 B. F. Sørensen and T. K. Jacobsen, “Large scale bridging in composites: Rcurve and bridging laws,” Composites Part A, vol. 29, no. 11, pp. 1443–1451, 1998. View at: Publisher Site  Google Scholar
 B. F. Sørensen, E. K. Gamstedt, R. C. Østergaard, and S. Goutianous, “Micromechanical model of crossover fibre bridgingprediction of mixed mode bridging laws,” Mechanics of Materials, vol. 40, no. 45, pp. 220–234, 2008. View at: Publisher Site  Google Scholar
 K. R. Pradeep, B. N. Rao, S. M. Srinivasan, and K. Balasubramaniam, “Interface fracture assessment on honeycomb sandwich composite DCB specimens,” Engineering Fracture Mechanics, vol. 93, pp. 108–118, 2012. View at: Publisher Site  Google Scholar
 U. Stigh, “Damage and crack growth analysis of the double cantilever beam specimen,” International Journal of Fracture, vol. 37, no. 1, pp. R13–R18, 1988. View at: Publisher Site  Google Scholar
 J. G. Williams, “End corrects for orthotropic DCB specimens,” Composite Science and Technology, vol. 35, no. 4, pp. 367–376, 1989. View at: Publisher Site  Google Scholar
 L. P. Canal, M. Alfano, and J. Botsis, “A multiscale based cohesive zone model for the analysis of thickness scaling effect in fiber bridging,” Composites Science and Technology, vol. 139, pp. 90–98, 2017. View at: Publisher Site  Google Scholar
 B. D. Manshadi, E. FarmandAshtiani, J. Botsis, and A. P. Vassilopoulos, “An iterative analytical/experimental study of bridging in delamination of the double cantilever beam specimen,” Composites Part A: Applied Science and Manufacturing, vol. 61, pp. 43–50, 2014. View at: Publisher Site  Google Scholar
 H. Danielsson and P. J. Gustafsson, “A three dimensional plasticity model for perpendicular to grain cohesive fracture in wood,” Engineering Fracture Mechanics, vol. 98, pp. 137–152, 2013. View at: Publisher Site  Google Scholar
 A. B. de Morais, “A new fibre bridging based analysis of the double cantilever beam (DCB) test,” Composites: Part A, vol. 42, no. 10, pp. 1361–1368, 2011. View at: Publisher Site  Google Scholar
 S. Morel, C. Lespine, J. L. Coureau, J. Planas, and N. Dourado, “Bilinear softening parameters and equivalent LEFM Rcurve in quasibrittle failure,” International Journal of Solids and Structures, vol. 47, no. 6, pp. 837–850, 2010. View at: Publisher Site  Google Scholar
 J.L. Coureau, S. Morel, and N. Dourado, “Cohesive zone model and quasibrittle failure of wood: a new light on the adapted specimen geometries for fracture tests,” Engineering Fracture Mechanics, vol. 109, pp. 328–340, 2013. View at: Publisher Site  Google Scholar
 Z. P. Bazant and M. T. Kazemi, “Size effect in fracture of ceramics and its use to determine fracture energy and effective process zone length,” Journal of the American Ceramic Society, vol. 73, no. 7, pp. 1841–1851, 1990. View at: Publisher Site  Google Scholar
 F. A. M. Pereira, M. F. S. F. de Moura, D. Dourado, J. J. L. Morais, J. Xavier, and M. I. R. Dias, “Direct and inverse method to the determination of cohesive law of bovine cortical bone using the DCB test,” International Journal of Solids and Structures, vol. 128, pp. 210–220, 2017. View at: Publisher Site  Google Scholar
 American Society for Test Materials (ASTM), ASTM D14314, Standard Test Methods for Small Clear Specimens of Timber, ASTM, West Conshohocken, PA, USA, 2014.
 H. D. Bui, Fracture Mechanics, Springer, Berlin, Germany, 2006.
Copyright
Copyright © 2019 Yurong Shen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.