Next Article in Journal
Performance Comparison of CdTe:Na, CdTe:As, and CdTe:P Single Crystals for Solar Cell Applications
Next Article in Special Issue
Nonlinear Inverse Analysis for Predicting the Tensile Properties of Strain-Softening and Strain-Hardening UHPFRC
Previous Article in Journal
Fly Ash-Added, Seawater-Mixed Pervious Concrete: Compressive Strength, Permeability, and Phosphorus Removal
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimental and Numerical Study on the Failure Characteristics of Brittle Solids with a Circular Hole and Internal Cracks

1
College of Water Conservancy and Hydropower Engineering, Hohai University, Nanjing 210098, China
2
Sichuan Woneng Investment Group Co., Ltd. Chengdu 610031, China
3
Nanjing Hydraulic Research Institute, Nanjing 210098, China
*
Author to whom correspondence should be addressed.
Materials 2022, 15(4), 1406; https://doi.org/10.3390/ma15041406
Submission received: 11 November 2021 / Revised: 9 February 2022 / Accepted: 10 February 2022 / Published: 14 February 2022

Abstract

:
A stress analysis of a circular hole is one of the classical problems in mechanics. Internal cracks are inherent properties of materials, and they are mostly three-dimensional in form. However, studies on hole problems with three-dimensional internal cracks are still lacking. In this paper, internal cracks were generated in brittle materials containing circular holes based on 3D internal laser-engraved crack technology. Then, uniaxial compression tests were performed. The experimental results were compared with the existing literature, and theoretical and numerical simulation studies were carried out. The results show that: (1) The main crack shapes are the primary cracks and remote cracks. (2) The dynamic fracture characteristics existed in the formation of primary cracks and the surface of remote cracks. The tips of primary cracks were arc-shaped, and the surfaces of the remote cracks were curved. Remote cracks were tangential to the orifice where type III spear-like characteristics appeared. (3) The stress birefringence technology can be combined with 3D internal laser-engraved crack technology for internal crack stress information monitoring, the moire around the orifice was “flamboyant”, and the moire at the tip of the prefabricated crack was “petallike”. (4) The existence of internal cracks reduced the cracking and breaking load of the specimen, and compared with the intact orifice specimen, the upper primary crack, the lower primary crack, the remote crack and the failure load were reduced by 41.2%, 31.7%, 15.9%, and 32.3%, respectively. (5) The results of qualitative stress analysis of the orifice specimen were consistent with the initiation law of primary cracks and remote cracks. The K distribution based on M integral and the numerical simulation of crack propagation process based on the maximum tensile stress criterion were consistent with the law of primary crack growth. Compared with the current mainstream method of transparent rock research, 3D internal laser-engraved crack technology has certain advantages in terms of brittleness, crack authenticity, stress field visualization, and fracture characteristics, and the result will provide experimental and theoretical references for research on three-dimensional problems and internal cracks in fracture mechanics.

1. Introduction

The stress and failure analysis of circular holes is one of the classical problems in the mechanics, for example, the orifice fatigue fracture problems in aviation field [1,2,3]; the stress concentration problems of orifice bolts in mechanical engineering [4,5,6]; the problems of perforating hydraulic fracturing in oil and shale gas extraction [7,8,9]; the stability of surrounding rock of tunnels in hydraulic and geotechnical fields [10,11,12]. Research on the stress and failure analysis of circular holes mainly concentrates on three aspects: (1) Theoretical research. For example, Lamé et al. [13] gave a theoretical solution of a ring subjected to uniform internal and external loads based on the linear elasticity theory, which later became the basis of inelastic orifice problems; Muskhelishvili et al. [14] used the plane elastic complex function method to give the expression for the stress field and displacement field of the orifice problem; Lv et al. [15] and Fan et al. [16] obtained the analytical expression of an orifice with arbitrary shape by deriving the orifice mapping functions; Mendelson et al. [17] extended the elasticity problem of planar orifices to elastoplasticity. (2) Experimental studies, for example, Zhang et al. [18] conducted experimental studies on mechanical properties and failure models of marble under different numbers and diameters of orifices; Gong et al. [19] analyzed the failure processes of granite specimens with rectangular holes under triaxial stress conditions; Zhu et al. [20] discussed the failure processes of sandstone with double round holes under uniaxial compression. (3) Numerical simulation, for example, Li et al. [21] used Flac3d software to simulate the stress, strain, and damage process of plate samples with bilateral orifice; Li et al. [22] used particle flow discrete element software PFC to conduct numerical simulation research on the uniaxial compression failure process of orifices containing circular, rectangular, and gate shapes; Xie et al. [23] used the damage and fracture software RFPA-3D to simulate the three-dimensional fracture of cuboid samples with orifices.
With the continuous advances in fracture mechanics research, the mechanisms of interactions of pre-existing cracks and orifice have also been further studied. For example, Wu et al. [24] conducted experimental and numerical simulation studies on cement mortar samples with circular holes and cracks; Zhu et al. [25] discussed the strength characteristics, deformation characteristics, and fracture evolution processes of tabular sandstone with prefabricated orifice and fissures; Zhang et al. [26] made an in-depth analysis of the mechanisms of rock mass with orifice and fracture after grouting. However, most of the previous studies on cracks are limited to two dimensions, and the research on the orifice with 3D crack is quite rare. In fact, simply transforming 3D problems into 2D cannot fully reflect the mechanical properties of defective materials [27]. The reasons contributing to the lack of research on 3D internal cracks are mainly reflected in the following aspects: (1) Internal crack visualization problems, for example, Li et al. [28] used the cement mortar pouring method to bury three-dimensional cracks, but they could not observe the propagation process of three-dimensional cracks during the test process. Although CT scanning [29] and acoustic emission technology [30] have been introduced, they still cannot solve the problems of scanning accuracy, accurate positioning, and real-time observation. (2) The generation of internal cracks, for example, with the continuous development of research, transparent brittle materials have been developed to facilitate the observation of internal cracks. There are three main methods to generate internal cracks in transparent brittle materials: (1) The cut and paste method [31]. This method was firstly proposed by Adams [31] in 1978. This is carried out by cutting the semi-circular surface crack into two pieces of plexiglass and forming an “inner crack” by almosaics, which obviously destroys the integrity of the original test when using this method. (2) Embedded casting method [32]. This method was firstly proposed by Dyskin [32] in 1994. It is formed by pouring epoxy resin into a metal or mica sheet that is prepositioned in a mold (regarded to be an internal crack). This method solves the problem of sample integrity of the cutting and pasting method and has been regarded as a mainstream research method for three-dimensional internal crack by scholars all over the world. However, this method requires casting molding samples at very low temperature (minus 20~30 degrees) to maintain brittleness, and the highest brittleness can only reach 1/3~1/7. At the same time, the test discreteness under low temperature is very large, the test success rate is low. In addition, regarding the heterogeneous sheet as an internal crack is not consistent with the actual situation, and it is difficult to imagine that the crack surface is a hard metal or mica sheet. (3) 3D printing method. This method is considered to be the most promising new technique for solving 3D internal crack generation problems. The advantage lies in that it can be modeled digitally and can produce arbitrary internal cracks. For example, Jv et al. [33] used 3D printing to prepare coal and rock models with complex fracture networks; Liu et al. [34] prepared transparent rocks using 3D printing. However, 3D printing technology still has many problems to be solved: firstly, the method still uses resin material for printing, which cannot guarantee enough brittleness; secondly, 3D printing internal cracks still needs to use heterogeneous support materials, essentially similar to the above embedded pouring method.
Therefore, this paper adopted a series of methods to solve the above problems. First, 3D-ILC technology proposed by the author [35] was used to generate the 3D internal cracks in a brittle solid without any heterogeneous support materials and, at the same time, did not have any impact on the surface of the brittle solid; The transparent brittle material glass was selected to solve the problems of observing the process of crack propagation. At the same time, the brittleness was higher, which provided a reference for similar research on natural brittle materials such as rock. The stress birefringence of glass lays a foundation for stress visualization. Based on this, uniaxial compression tests were carried out on the specimens with different embedded depths of internal cracks. The fracture mechanism and mechanical characteristics of the specimens with three-dimensional internal cracks were revealed through stress moire, failure pattern analysis, and numerical simulation. The results provide a physical test basis for the corresponding theoretical research and correction.

2. Research Scheme and Sample Preparation

2.1. 3D-ILC

3D-ILC was proposed by Wang et al. [35] in 2017, which realized the fabrication of arbitrary internal cracks in a material without any impact on the surface and be called the 3D internal laser-engraved crack technology.
It is difficult to make a purely closed internal crack. It can be imagined that it is impossible to make a “surgically” precise structural change in the interior of a complete object without affecting other parts. However, electromagnetic and wave physical fields can penetrate objects and be used in the interior. Wang et al. [35] finally proposed that 3D-ILC acts on the interior of materials through electromagnetic fields to form plasma (the fourth state of matter) blasting, thus producing macroscopic pure closed internal cracks. Due to the penetration of electromagnetic field, the internal crack can be made without any influence on the surface by controlling the parameters.

2.2. Sample Material

The glass was selected for testing due to the following advantages: (1) due to the homogeneity and isotropy, glass strictly follows Hooke’s law before crack propagation, and is a classic material in brittle solid fracture mechanics; (2) glass is a stress birefringent photosensitive material, which provides a basis for direct observation of stress distribution; (3) as an index to measure the brittleness of materials, the pull-to-pressure ratio of glass can reach 1/13~1/33, which is similar to that of rock materials (1/12~1/35), and can better simulate the brittleness of rock.
As is pointed out by Li [36] in his book, All Pioneering Tests of Solid Mechanics Use Glass as the Basic Sample Material, and Fracture Mechanics is No Exception, the Griffith fracture criterion, the basis of modern fracture mechanics, uses glass as a test material. Subsequently, many scholars have carried out classical tests based on glass. For example, Roesler et al. [37] conducted three-dimensional Hertz contact fracture tests on glass materials and observed conical crack morphology; Knauss et al. [38] studied the morphology of pure type Ⅲ crack with glass material and observed the “quasi-spiral” failure characteristics; Sommer et al. [39] used glass rods to conduct tensile and torsion tests and observed the classic fracture form of I–III “double spearhead”, which became the permanent cover of the international journal, Engineering Fracture Mechanics.

2.3. Experimental Conditions

For comparison, two types of brittle solid samples were used to carry out the tests: (1) intact samples (group A); (2) samples with internal crack (group B and group C), and the design, number, test purpose, and research content of each sample are shown in Table 1.
(1)
Group A (complete orifice sample without internal cracks). The control group with no internal cracks, the samples were loaded to crushing to provide comparison for the orifice samples (Group B and Group C) with internal cracks of different buried depths, which is marked as A0.
(2)
Group B (orifice samples with internal cracks of different buried depths). Samples buried at depths (d) 0.5, 1, 1.5, and 2 cm, which are marked B0.5–B2.5. The samples were loaded to crushing to observe the propagation of internal crack, the characteristic loads and other information during the loading process, and the fracture morphology was observed after the specimen was crushed.
(3)
Group C. The samples were the same as group A and B, but the loading stopped when the main crack and the remote crack appeared, which were recorded accordingly. The samples were marked Ci-j, where i ranges from 0 to 2, representing samples without internal cracks and samples with cracks in different buried depths. The values of j range from 1 to 2, representing the occurrence of main cracks and remote cracks.

2.4. Sample Preparations

(1)
Group A (samples without internal cracks): The sample was a cubic sample of 80 × 80 × 40 mm3 with a circular hole in the center, with a diameter that was 20 mm as shown in Figure 1a.
(2)
Group B and C (samples with internal cracks): The sample was a cubic sample of 80 ×80 × 40 mm3, with a circular hole in the center and a circular inner crack, with a diameter that was 2c = 13 mm. The d between the inner crack and the top surface of the sample was 0.5, 1, 1.5, and 2 cm, respectively, as shown in Figure 1b.
Figure 1. Sample design: (a) samples without internal cracks; (b) samples with internal cracks.
Figure 1. Sample design: (a) samples without internal cracks; (b) samples with internal cracks.
Materials 15 01406 g001

2.5. Test Instruments

The SUNS-650W loading equipment (Tinius kuli, Suzhou, China) was applied in this experiment, which had two loading modes. The maximum loading pressure was 600 kN, and the displacement range was 0–200 mm. In our experiment, the force control loading mode was utilized, and the loading rate was 0.5 kN/s. The equipment is shown in Figure 2, and the loading mode is shown in Figure 3. Stress concentration is more sensitive in brittleness test, and the unevenness or unsmoothness of the end will lead to cracks formed by stress concentration at the end of the sample, which will affect the test quality. To avoid this and reduce the end friction effect, smooth PMMA(plexiglass) is placed at the top and bottom.

2.6. Birefringence Effect Observation System

Glass and other amorphous media usually exhibit optical isotropy in a stress-free state. When subjected to stress, the refractive index characteristics change, showing optical anisotropy. When a beam of light passes through glass with internal stress, it will produce two beams of light with different propagation speed, ordinary light o following the law of refraction, and extraordinary light e not following the law of refraction, which is stress birefringence. According to the stress-optical law, when the incident light incident on the test object due to the birefringence effect, the relationship between the principal stress and the corresponding refractive index is as follows:
n1n2 = (C1C2)(σ1σ2)
where σ1 and σ2 are the principal stresses of the sample under uniaxial compression stress state, respectively; n1 and n2 are the refractive indices of the direction of σ1 and σ2, respectively; C1 and C2 are the material stress-optical coefficients. The optical path difference (Δ) generated by polarized light passing through the sample is:
Δ = (n1n2)h
where h is the thickness of the medium. According to Equations (1) and (2), the corresponding relationship between stress and optical quantity is established:
(σ1σ2) = Δ/hC
Figure 2. Loading system.
Figure 2. Loading system.
Materials 15 01406 g002
Figure 3. Uniaxial compression loading mode.
Figure 3. Uniaxial compression loading mode.
Materials 15 01406 g003
The arrangement of orthogonal plane-polarized light field is shown in Figure 4.

3. Experimental Results

3.1. Failure Process Analysis

3.1.1. Orifice Sample without Internal Cracks (Group A)

As can be seen from Figure 5a, the typical failure process of the sample under uniaxial compression can be divided into four stages: (1) stable stage; (2) main crack propagation stage; (3) remote crack propagation stage; (4) failure stage.
(1)
Stable stage: This stage persists until the loading force to 131.16 kN. No crack propagation was observed in this stage;
(2)
Main crack propagation stage: when loading to 131.16 kN, the upper and lower main cracks initiated almost simultaneously, and the main crack’s growth rate gradually decreased after the crack initiation.
Figure 5. Crack propagation processes of the orifice specimen. (a) A0; (b) B0.5; (c) B1; (d) B1.5; (e) B2.
Figure 5. Crack propagation processes of the orifice specimen. (a) A0; (b) B0.5; (c) B1; (d) B1.5; (e) B2.
Materials 15 01406 g005aMaterials 15 01406 g005b
The glass can be regarded as a completely elastic body before cracking. According to The Ramet solution (Equation (4)), the annular stress distribution at the orifice is similar to a “peanut shell”, which is shown in Figure 6; that is, the maximum annular tensile stress concentration occurs at the top of the orifice. Therefore, under the action of the maximum annular tensile stress, a mode I tensile crack first appears above and below the orifice, and the direction is perpendicular to the annular tensile stress.
{ σ r = 0 σ θ = p + 2 pcos 2 θ τ r θ = 0
It is noteworthy that the crack tip is circular in the main crack propagation process. Assume that the length and length radius of the crack tip are ellipses of a and c, respectively, the vertical stress is σA, and the following equations can be obtained:
K Ι = φ σ A c 0.5
φ   ( a c ,   β ) = π 0.5 [ cos 2 β + ( c a ) 2 sin 2 β ] 0.5 /   E   ( a c )
E   ( a c ) = 0 π 2 [ 1 ( 1 c 2 a 2 )   sin 2 ϕ ] 0.5 d ϕ
where KΙ is the mode I stress intensity factor; c is the crack length; φ is the geometry item; E(a/c) is the elliptic integral; ϕ is the dummy variable. When a/c ≠ 1, then:
K Ι ( β = π 2 ) / K Ι ( β = 0 ) = ( c a ) 0.5 1
Therefore, without external interference, crack propagation always tends to form a circular tip with a/c = 1.
(3)
Remote crack propagation stage: The orifice sample with no internal cracks did not continue to exhibit “bisection” fracture failure along the main crack, but it produced two pairs of remote cracks at the symmetrical positions above and below the orifice when the load was 438.25 kN. The surface was curved and tangent to the orifice. The line between the tangent point and the center of the circular hole was at an angle of approximately 30° with the horizontal axis. The occurrence time of the upper and lower remote cracks is similar;
(4)
Failure stage: When the failure load reached 587.13 kN, the sample showed violent “explosive” failure, and the sample was detrital with a long range of avalanche. This characteristic is typical of brittle material failure under ballast loads of a rigid testing machine.

3.1.2. Orifice Sample with Internal Cracks (Group B)

As can be seen from Figure 5b–d, the typical failure process of the sample with three-dimensional internal cracks under uniaxial compression was similar to that of the complete sample without cracks, which can be divided into four stages: (1) stable stage; (2) main crack propagation stage; (3) remote crack propagation stage; (4) failure stage.
(1)
Stable stage: This stage persisted until a loading force of 77.12 kN. No crack propagation was observed during this stage.
(2)
Main crack propagation stage: When loading to 77.12 kN, the upper main crack first initiated, and when the load reached 89.61 kN, the lower main crack initiated. The main crack initiation load of samples containing internal cracks was less than that of samples not containing internal cracks. Meanwhile, the internal cracks also influenced the stress field of the sample, which led to the upper main crack becoming easier to initiate than the lower main crack. What should be noticed is that for the samples with lower buried depths (d = 0.5, 1, and 1.5 cm), the main crack would not “penetrate” the prefabricated internal crack. However, when the buried depth is higher (d = 2 cm), the main crack would “penetrate” the prefabricated internal crack;
(3)
Remote crack propagation stage: When loaded to 368.72 kN, remote cracks appeared, and the crack form was similar to that of the sample without internal cracks, but the initiation load of remote cracks was smaller than that of the sample without internal cracks;
(4)
Failure stage: When the load reached 397.61 kN, there was a violent “explosive” failure, and the failure load was significantly lower than that of the complete orifice sample without crack.

3.2. Stress Moire Laws

The changes of stress moire in the test during the loading process is shown in Figure 7, which mainly shows the following rules:
(1)
When the pressure was relatively lower, the stress moire presented the shape of a “flame” in the top and bottom of the orifice for the condition of the sample without internal cracks. However, for the condition of the sample with internal cracks, the “petal shape” stress moire was observed in the tips of internal crack;
(2)
With the increase in load, the petal moire at the tip of the prefabricated crack gradually became lighter, while the flame moire at the upper and lower part of the orifice gradually highlighted;
(3)
The appearance of the main crack in the upper and lower part of the orifice changed the moire distribution of the “flame shape” and made the moire appear similar to an “inverted triangle” along the direction of the main crack propagation.
Due to the difficulty of observing internal cracks, visualization of the stress field inside rocks is a difficult problem. The traditional stress birefringence is basically used in plate tests without cracks or with penetrating cracks. In addition to the 3D-ILC technology, Jv et al. [33] carried out effective research on visualization of the internal stress field by combining 3D printing transparent rock materials with stress freezing and slicing technology.
This method using stress freezing technology to keep the strain unchanged, then one can observe the cut sliced section with stress moire, achieving three-dimensional stress field distributions inside the sample observations. This method cannot be used for real-time monitoring. However, it is still a creative and effective method for the viewing the “black box” of internal cracks in rocks.
In this paper, 3D-ILC combined with stress birefringence technology was used to try to visualize the stress field of the internal crack. It is a different idea. The stress mottling shows high quality and the “petal” distribution at the crack tip indicating that it can reflect the stress distribution characteristics of the internal crack tips. At present, the authors are also carrying out in-depth quantitative theoretical analysis and research and compiling corresponding procedures, and they plan to present thematic results in the future. Here, only qualitative analysis is carried out.

3.3. Crack Initiation and Failure Load

Figure 7 shows the characteristic loads of the orifice sample without internal cracks and the orifice sample with internal cracks. It should be noted that the characteristic load of the cracked orifice samples in the figure was the average load of samples with different buried depths.
Figure 7. Characteristic load distributions.
Figure 7. Characteristic load distributions.
Materials 15 01406 g007
As can be seen from the figure, the characteristic loads of each sample showed the following rules:
(1)
For samples with internal cracks, the proportion of crack initiation load to final failure load was 19.4% for the upper main crack, 22.5% for the lower main crack, and 92.7% for the remote crack. For the sample without internal cracks, the upper main crack was 22.3%, the lower main crack was 22.3%, and the remote crack was 74.6%;
(2)
The existence of prefabricated cracks greatly reduced the values of each characteristic loads. Compared with the orifice sample without internal cracks, the average load of the upper main crack, lower main crack, remote crack, and final failure of the specimen with internal cracks were reduced by 41.2%, 31.7%, 15.9%, and 32.3%, respectively.

4. Fractography Characteristics

As an independent discipline, “fractography” analyzes the process, type, nature, and mechanism behind material failure by studying the morphology and properties of macro and micro fractography. It can be said that “fractography is a crime scene of mechanics”. Fractography also has many applications in the field of rocks such as scanning electron microscopic (SEM) fracture analysis and fractal rock mechanics. However, in the current research field of transparent rock-like materials, most of the research is focused on the analysis of propagation path and failure modes, and few studies are combined with fractography. In this paper, significant fracture characteristics and mechanisms behind them are presented.

4.1. Main Crack Surface Characteristics

4.1.1. Circular Arc Characteristics of Crack Tip

No matter the sample containing the internal crack or not, the crack tip is circular in shape, and its mechanism can be explained in Section 3.1.1.
As shown in Figure 8, the blue line represents the crack tip when the specimen initiates, and the red line represents the crack tip when the specimen unloads. It can be seen that there were obvious arc-shaped irregular features on the fracture. The rectangular coordinate system, as shown in Figure 8, was defined to fit the arc of the crack surface at the fracture of sample C0-1. The arc equation of crack tip 1 (Equation (9)) and crack tip 2 (Equation (10)):
( x 20 ) 2 + ( y 1.81 ) 2 = 28.54 2
( x 20 ) 2 + ( y 5.99 ) 2 = 31.82 2
According to Equations (9) and (10), when the main crack expands, the circular radius of the crack front increases continuously, and the center of the circle is located on the axial center line of the orifice.
Figure 8. Arc propagation characteristics at crack tip.
Figure 8. Arc propagation characteristics at crack tip.
Materials 15 01406 g008

4.1.2. Dynamic Fracture Characteristics

Dynamic fracture occurs at the crack initiation moment of the main crack, that is, the “mirror zone”, “atomization zone”, and “feather zone” are presented spreading outward from the crack initiation point as shown in Figure 9a, and the results of uniaxial tensile test in previous experiments [40] are shown in Figure 9c.
Under the action of increasing load, these defects in the solid begin to connect and form cracks and expand outwards. When the crack growth reaches a certain speed, the equilibrium state of static crack growth is broken. At this point, the volume elements in the area near the crack tip inside the sample are affected by non-equilibrium forces, and these volume elements are accelerated to obtain kinetic energy accordingly. The original static solution cannot be established if kinetic energy is ignored, and the static fracture system becomes a dynamic fracture system.
In the dynamic system, the crack growth rate increased continuously until it reached the limit value, and this process showed a radial mirror region, atomization region, and feather region on the fracture surface. These regions also recorded different stages of kinetic energy dissipation: the mirror region was a semicircular plane, the fracture was bright and smooth, and the reflectivity was high, which was formed when the crack propagation rate was low. The atomization zone is the transition zone, which was the semi-ring region outside the mirror area and is the transition zone from slow crack growth to fast crack growth. It is basically in the same plane with the mirror area, and its appearance is characterized by “atomization”; that is, the reflectivity decreases. The reason is that the single smooth crack surface cannot release the kinetic energy of the crack completely under high-speed crack growth, and small-scale fracture occurs on the crack surface, resulting in the increase of surface roughness. The feather area radiated from the atomization area to the surrounding three-dimensional space. The surface of the fracture was very rough when viewed from the plane of the mirror and atomization area. There was still crack propagation outside the plane, and its mechanical mechanism was the dynamic crack bifurcation.
According to the experimental fracture characteristics, this paper analyzed the dynamic crack tip field distortion theory. The state of the near field at the tip of an expanding crack varies with the rate of crack growth. Definition of a “dynamic stress intensity factor”:
K = K ( v t / v ) 0.5
Then, the angle function, fθθ(θ), has the same form as the static crack tip stress solution:
fθθ(θ) = σθθ(2πr)0.5
σij = K(2πr)−0.5fij(θ)
Figure 9. Dynamic crack bifurcation and the characteristic of the Wallner line. (a) Main crack fractography; (b) Main crack characteristics; (c) Wallner line in previous studies; (d) Wallner line formation.
Figure 9. Dynamic crack bifurcation and the characteristic of the Wallner line. (a) Main crack fractography; (b) Main crack characteristics; (c) Wallner line in previous studies; (d) Wallner line formation.
Materials 15 01406 g009
Figure 10 shows the component, σθθ, calculated for some selected relative crack growth rates, vt. The instability propagation of the crack occurs when the propagation velocity approaches the theoretical limit rate (the rate of Rayleigh surface waves). As can be seem from Figure 11, when the velocity approached the theoretical limit rate 0.58vt (vt is the velocity of shear wave propagating in solid), the maximum local tensile stress shifts from 0° on the original crack surface to 25°~35° on the inclined surface. Similarly, the angular distribution of mechanical energy release rate also showed a similar trend. Therefore, in isotropic materials, the crack bifurcation was the inevitable result of the dynamic characteristics of the crack field.
In this paper, the dynamic fracture characteristic occurred near the orifice area. In the center of the axis of the orifice and the crack, the source can clearly be observed. The dynamic fracture was along the horizontal direction. The authors speculated that when the crack along the axis of the orifice rate of speed and along the main direction of crack propagation rate, it gradually became smaller, presenting the static fracture characteristics.

4.1.3. Secondary Wallner Lines

The Primary Wallner Lines were described in detail in the three-point bending test carried out by the author’s research group [41]. The primary Wallner lines are caused by the elastic stress waves generated by the crack tip during crack propagation, which are reflected at corner points and interact with the expanding crack tip. The propagation direction of the crack tip deviates from the normal path momentarily, and the crack oscillates, leading to the fluctuation of the propagation path of the crack. The reflection of light from these fluctuations makes the crack tip’s interaction with shear waves visible on the surface as continuous arcs. Secondary Wallner lines are caused by the dynamic fracture characteristics of cracks (usually the crack growth reaches its limit rate), and the elastic stress wave generated by the expanded crack front is reflected by the rough surface of the dynamic fracture characteristics of cracks, and finally interacts with the expanded crack front to form, as can be seen in Figure 9d.

4.2. Remote Crack Surface Characteristics

4.2.1. Dynamic Fracture Characteristics

The remote crack surface had the same dynamic fracture characteristics as the main crack surface, but the dynamic fracture characteristics of the main crack surface appeared close to the orifice, while the dynamic fracture characteristics of the remote crack surface appeared in the middle of the remote crack surface, and the crack propagation direction was in the mirror zone, atomization zone, and feather zone. The comparison with the literature [28] is shown in Figure 11b.
Figure 11. Surface features of remote crack. (a) Remote crack surface characteristics; (b) Dynamic fracture characteristics; (c) Spear torn Characteristics.
Figure 11. Surface features of remote crack. (a) Remote crack surface characteristics; (b) Dynamic fracture characteristics; (c) Spear torn Characteristics.
Materials 15 01406 g011

4.2.2. Characteristics of Tearing III—Shaped Crack

The characteristics of mode III crack in the tearing zone appear at the point where the remote crack is tangent to the orifice. In this paper, combined with the numerical simulation in Section 5, it is judged that this part is mode III crack initiation.
The spear-shaped cracks in the tearing region are closely distributed in line near the orifice. As the cracks continue to extend outwards, 2–3 short lines will converge and merge and continue to extend. The side is stepped, and the linear shape near the orifice is thin, while the outer line shape is thick.
Friedel [42] gave a theoretical explanation for the characteristics of “tearing” cracks caused by type III shear dislocation: There is linear tension (T) at the front end of the staggered crack. This tension is the unit energy increased by each unit length of the crack front, and it is assumed that T is independent of the direction of the crack front and the speed of the crack. The step also has a drag force of 2γh, γ is the energy required to form a unit surface at the step, and H is the height of the step, which is shown in Figure 12. For the step with equal spacing, the crack expands symmetrically with the step as the center, θ′ = θ″, so that the whole front end of the crack moves forward under the following conditions:
2 γ h = 2 T sin θ
If the distance between steps is not equal, as shown in Figure 12c, asymmetric forces will appear at the steps, i.e., θ′ ≠ θ″, thus Tsinθ′ ≠ Tsinθ″. The asymmetric stress will lead to the steps with close spacing being close to each other in the process of crack propagation, and these steps will eventually merge into a new step. As the steps become closer to each other, the θ′/θ″ ratio increases and the new and old steps merge further, resulting in a “binary tree” shape as shown in Figure 13.

4.3. Crack “Penetration”

For the specimen of B2, in the process of main crack propagation, the prefabricated crack would be penetrated. The fractography of the sample is shown in Figure 14.
When the primary crack penetrated the prefabricated crack, it cut off the prefabricated crack and, at the same time, produced mode III spear fracture features on the upper and lower sides of the prefabricated crack (see Section 4.2.1 for its cause). What should be noticed is that the spear crack at the upper part of the prefabricated crack was coarser, while the lower part is thinner, and there is a smaller spear crack extending at the distal end, showing a similar “tentacle shape” shape.
Figure 14. Main crack “penetrating” prefabricated crack. (a) Experimental results; (b) Schematic diagram.
Figure 14. Main crack “penetrating” prefabricated crack. (a) Experimental results; (b) Schematic diagram.
Materials 15 01406 g014

5. Numerical Simulations

Three-dimensional crack propagation in fracture mechanics is a difficult problem both theoretically and numerically. At present, there are almost no numerical simulation results involving three-dimensional crack and hole interactions as well as crack penetration under complex stress states. Limited by the existing theoretical and numerical simulation level, qualitative analysis will be carried out in this section, which mainly includes numerical simulation from the following three aspects: (1) the influence of different prefabricated crack burial depths on main crack initiation; (2) study on the propagation process of the main crack; (3) the influence of main crack generation and propagation on remote crack.

5.1. Stress Intensity Factor Calculation

The stress intensity factor is calculated by means of m-integral, assuming that the interaction J-integral of the superposition of the two equilibrium states is:
J ¯ = Γ ( σ i j ( 1 ) + σ i j ( 2 ) ) u i ( 1 ) + u i ( 2 ) x ( W ( 1 ) + W ( 2 ) ) δ i j q x j d s = J ¯ ( 1 ) + J ¯ ( 2 ) + M ¯ ( 1 , 2 )
For linear elastic materials, strain energy density can be expressed as W = σ ij ε ij / 2 . According to Betti reciprocity theorem:
W ( 1 , 2 ) = σ ij ( 1 ) ε ij ( 2 ) = σ ij ( 2 ) ε ij ( 1 )
where:
J ¯ ( 1 ) = Γ ( σ i j ( 1 ) u i ( 1 ) x 1 W ( 1 ) δ i j ) q x j d s
J ¯ ( 2 ) = Γ ( σ i j ( 2 ) u i ( 2 ) x 1 W ( 2 ) δ 1 j ) q x j d s
M ¯ ( 1 , 2 ) = [ σ i j ( 1 ) u i ( 2 ) x 1 + σ i j ( 2 ) u i ( 1 ) x 1 1 2 ( σ i j ( 1 ) ε i j ( 2 ) + σ i j ( 2 ) ε i j ( 1 ) ) δ 1 j ] q x j d s
Using the relationship between energy release rate, G, and stress intensity factor, K, under two equilibrium states:
G = J = 1 ν 2 E ( K I ( 1 ) + K I ( 2 ) ) 2 + 1 ν 2 E ( K II ( 1 ) + K II ( 2 ) ) 2 + 1 + ν E ( K III ( 1 ) + K III ( 2 ) ) 2
J ( 1 ) = 1 ν 2 E ( K I ( 1 ) ) 2 + 1 ν 2 E ( K II ( 1 ) ) 2 + 1 + ν E ( K III ( 1 ) ) 2
J ( 2 ) = 1 ν 2 E ( K I ( 2 ) ) 2 + 1 ν 2 E ( K II ( 2 ) ) 2 + 1 + ν E ( K III ( 2 ) ) 2
M ( 1 , 2 ) = 2 [ 1 ν 2 E K I ( 1 ) K I ( 2 ) + 1 ν 2 E K I I ( 1 ) K I I ( 2 ) + 1 + ν E K I I I ( 1 ) K I I I ( 2 ) ] d s
By selecting a specific auxiliary field, the solution of the stress intensity factors KI, KII, and KIII at the crack tips can be obtained by using the interaction M-integral of the far field through Equation (23).
State 1 is selected as the real state of the problem, state 2 is represented by the asymptotic solution of pure mode I crack:
K I ( 2 ) = 1 ,   K II ( 2 ) = 0 ,   K III ( 2 ) = 0
We can obtain the equations from Equations (23) and (24):
M ( 1 , 2 ) = 2 1 ν 2 E K I
Therefore:
K I = E 2 ( 1 ν 2 ) M ( 1 , 2 ) = E 2 ( 1 ν 2 ) M ¯ ( 1 , 2 )
Similarly, state 1 is selected as the real state of the problem, state 2 is represented by the asymptotic solution of pure type II crack, which can then be written as:
K II ( 1 ) = K II , K II ( 2 ) = 1 , K I ( 2 ) = 0 , K III ( 2 ) = 0
We can obtain from (16) and (20) that:
M ( 1 ,   2 ) = 2 1 ν 2 E K II
K II = E 2 ( 1 ν 2 ) M ( 1 , 2 ) = E 2 ( 1 ν 2 ) M ¯ ( 1 , 2 )
Similarly, state 1 is selected as the real state of the problem, state 2 is represented by the asymptotic solution of pure type II crack, and we can obtain:
K III ( 1 ) = K III , K III ( 2 ) = 1 , K I ( 2 ) = 0 , K II ( 2 ) = 0
We can then obtain the following equations from Equations (15) and (22):
M ( 1 , 2 ) = 2 1 + ν E K III
K III = E 2 ( 1 + ν ) M ( 1 , 2 ) = E 2 ( 1 + ν ) M ¯ ( 1 , 2 )
Then the stress intensity factor can be obtained from the above equations.

5.2. Crack Propagation Criterion

In this paper, the maximum tensile stress criterion (MTS) is adopted as the crack growth criterion. The crack will propagate along the direction of the maximum tensile stress σθmax, while the circumferential stress is related to mode I stress intensity factor, which can be expressed as:
K I r ( θ ) = σ θ θ 2 π r = cos θ 2 [ K I cos 2 θ 2 3 2 K II sin θ ]
The cracking angle θ0 can then be obtained:
θ 0 = arccos 3 K 2 II ± K 4 I + 8 K 2 I K 2 II K 2 I + 9 K 2 II

5.3. Numerical Models

Three numerical models are established as shown in Figure 15: Model 1: complete orifice model without internal crack; Model 2: the main crack initiation model of the orifice (the size of the main crack is defined according to the experiment results in this paper); Model 3: the model with internal cracks of different buried depths (taking 2 cm as an example). The boundary condition is that the stress boundary is imposed on the upper part of the model and the displacement constraint of fixed three directions is imposed on the lower part of the model.

5.4. Numerical Results

5.4.1. Effects of Different Buried Depths on the Initiation of Main Crack

Figure 16 shows the distributions of the maximum tensile stress in Model 1 and Model 3 and the variations of the maximum tensile stress with the prefabricated cracks at different buried depths. The existence of prefabricated cracks greatly changed the stress distributions around the orifice. Under the same load, the maximum tensile stress of the orifice in Model 1 increased by 63.74% on average compared with Model 3. Therefore, the main crack of the orifice sample with internal crack initiates first, which is consistent with the experimental laws. Secondly, by comparing the main cracks at different buried depths in Model 3, it can be found that the greater the buried depth is, the greater the maximum tensile stress around the orifice is. Therefore, in this paper, the main crack penetration phenomenon of group B2 samples is explained as follows: (1) the maximum tensile stress at the orifice of this group of specimens is the largest, and the initiation and propagation of the main crack are more sufficient, so the main crack length is longer; (2) the prefabricated crack is closest to the orifice and is most easily penetrated by the main crack.

5.4.2. Main Crack Propagation Process

The main crack propagation process of Model 1 is shown in Figure 17. In order to explain the general rule, the stress intensity factor was normalized, and the relative stress intensity factor was defined as Ki/|KImax|. The variation rule of the relative stress intensity factor with the extended time step is shown in Figure 18.
The propagation laws of the main crack can be seen from the figure and listed below:
(1)
After initiation, the main crack presents self-similar propagation along the crack surface, which is consistent with the experimental phenomenon. According to the stress intensity factor distributions, the Mode 2 and Mode 3 stress intensity factors of the main crack growth process are 0, which can determine that the main crack is pure mode I fracture;
(2)
In the process of main crack propagation, the crack tip presented a circular arc feature. According to the distribution of stress intensity factor, it can be seen that the part near the main crack to the sample surface had a larger KI, while the stress intensity factor in the middle of the main crack is smaller. Therefore, the part near the sample surface had a faster crack growth, while the middle part had a slower crack growth. Therefore, the radius of the arc gradually increased with the growth of the main crack, which is consistent with the fracture characteristics of the main crack in Section 4.1.1 of this paper;
(3)
It can be seen from the pattern of Mode I stress intensity factor over time that the mode I stress intensity factor decreases gradually during the growth of the main crack, which is consistent with the experimental phenomenon that the growth rate of the main crack decreases gradually after initiation in this paper.

5.4.3. Effects of Main Crack Generation and Propagation on Remote Crack

The distributions of torsional shear stress (S13) between Model 1 and Model 2 are shown in Figure 19. The value of shear stress at the shear concentration of the complete sample (Model 1) (170.1 kpa) was 40.97% of that of the sample with main crack (Model 2) (415.2 kpa). It can be seen that the initiation of the main crack leads to the concentration of the torsional shear stress at 45º, resulting in the initiation of the remote crack tangent to the orifice. Therefore, the remote crack initiation mode was the Mode 3 crack initiation mode, which was also verified by the characteristics of the Mode 3 spear fracture at the Section 4.2 of remote crack tangent to the orifice.

6. Conclusions

(1)
No matter whether the orifice sample contains internal crack or not, the crack morphology is main crack and remote crack. At the initiation of the main crack, there is a dynamic fracture feature near the orifice, and the tip is arc-shaped. The radius of the arc gradually increases in the process of propagation, while the growth rate becomes slower and slower. When the prefabricated crack buried depth is relatively larger, it “penetrates” the prefabricated crack, and presents a "tentacle-like" shape at the surface of the main crack. The surface of the remote crack is curved and tangential to the orifice. The tangential part presents mode III spear feature, and dynamic fracture occurs on the crack surface;
(2)
Stress birefringence can be used to qualitatively monitor the dynamic change of stress in the sample with 3D-ILC internal crack. The existence of internal crack changes the original “flame” moire distribution of the orifice sample, and the moire at the prefabricated crack tip presents "petal" characteristics;
(3)
The presence of prefabricated cracks reduced the characteristic load of the sample, the upper main crack decreased 41.2%, the lower main crack decreased 31.7%, the remote crack decreased 15.9%, and the failure load decreased 32.3%;
(4)
The results of qualitative stress analysis of the orifice sample were consistent with the crack initiation laws of the main crack and remote crack, and the numerical simulation of the crack propagation process based on the K distributions of M integral and MTS was consistent with the main crack propagation laws;
(5)
3D-ILC can be applied to the study of crack propagation in brittle solids. Compared with the current mainstream method of transparent rock research, the embedded casting method has certain progress in approaching rock brittleness, crack authenticity, stress field visualization, fracture characteristics and other aspects.

Author Contributions

Conceptualization, H.W. and X.R.; methodology, C.L.; software, S.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Nature Science Found (Aging characteristics and long-term safety assessment and control method of surrounding rock-support system in deep buried tunnel under internal water pressure), grant number U1765204.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the data, models, or codes that support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ghalib, M.A.; Makky, S.M. Fatigue fracture of rotating circular disks. Eng. Fract. Mech. 1977, 9, 397–401. [Google Scholar] [CrossRef]
  2. Arutyunyan, A.R. A Fatigue-Fracture Criterion for Composite Materials. Dokl. Phys. 2019, 64, 394–396. [Google Scholar] [CrossRef]
  3. Molina, C.; Araujo, A.; Bell, K.; Mendez, P.F.; Chapetti, M. Fatigue life of laser additive manufacturing repaired steel component. Eng. Fract. Mech. 2020, 241, 107417. [Google Scholar] [CrossRef]
  4. Bolt, R.H.; Labate, S. Some Observations on the Acoustic Impedance of Small Orifices. J. Acoust. Soc. Am. 1948, 20, 596. [Google Scholar] [CrossRef] [Green Version]
  5. Tei, C.; Chin, K.; Vijayaraghavan, G.; Boltwood, C.M.; Shah, P.M. Paradoxical deterioration of left ventricular asynergy after administration of nitroglycerin. Am. J. Cardiol. 1985, 55, 677–679. [Google Scholar] [CrossRef]
  6. Romaniv, O.N.; Simin’Kovich, V.N.; Stepanov, V.G. Effectiveness of using a surface plastic strain to prevent crack development during cyclic loading. Mater. Sci. 1979, 15, 104–109. [Google Scholar] [CrossRef]
  7. Kanaun, S. On the hydraulic fracture of poroelastic media. Int. J. Eng. Sci. 2020, 155, 103366. [Google Scholar] [CrossRef]
  8. Zhao, K.; Stead, D.; Kang, H.; Gao, F.; Donati, D. Three-dimensional simulation of hydraulic fracture propagation height in layered formations. Environ. Earth Sci. 2021, 80, 435. [Google Scholar] [CrossRef]
  9. Kaurin, D.; Arroyo, M. Surface Tension Controls the Hydraulic Fracture of Adhesive Interfaces Bridged by Molecular Bonds. Phys. Rev. Lett. 2019, 123, 228102. [Google Scholar] [CrossRef] [Green Version]
  10. Zeng, K.H.; Xu, J.X. Unified semi-analytical solution for elastic-plastic stress of deep circular hydraulic tunnel with support yielding. J. Cent. South Univ. 2013, 20, 1742–1749. [Google Scholar] [CrossRef]
  11. Zhou, L.; Kai, S.U.; Zhou, Y.; Wen, X. Hydro-mechanical coupling analysis of pervious lining in high pressure hydraulic tunnel. J. Hydraul. Eng. 2018, 49, 313–322. [Google Scholar]
  12. Yang, T.; Zhu, W.; Yu, Q.; Liu, H. The role of pore pressure during hydraulic fracturing and implications for groundwater outbursts in mining and tunnelling. Hydrogeol. J. 2011, 19, 995–1008. [Google Scholar] [CrossRef]
  13. Lamé, G.; Lequien, M.; Pionnier, P.A. Interpretation and limits of sustainability tests in public finance. Appl. Econ. 2014, 46, 616–628. [Google Scholar] [CrossRef]
  14. Muskhelishvili, L.; Latendresse, J.R.; Kodell, R.L.; Henderson, E.B. Evaluation of Cell Proliferation in Rat Tissues with BrdU, PCNA, Ki-67(MIB-5) Immunohistochemistry and In Situ Hybridization for Histone mRNA. J. Histochem. Cytochem. 2003, 51, 1681–1688. [Google Scholar] [CrossRef] [Green Version]
  15. Lv, A.; Wang, Q. A new method for solving mapping function of roadway with arbitrary section shape by using optimization technique. J. Rock Mech. Eng. 1995, 03, 269–274. [Google Scholar]
  16. Fan, G.; Tang, C. The outfield mapping function of non-circular cavity is solved by the multiplication of three absolutely convergent series. J. Rock Mech. Eng. 1993, 03, 255–264. [Google Scholar]
  17. Mendelson, S.D.; Gorzalka, B.B. A facilitatory role for serotonin in the sexual behavior of the female rat. Pharmacol. Biochem. Behav. 1985, 22, 1025–1033. [Google Scholar] [CrossRef]
  18. Zhang, C.; Tang, J.; Teng, J. Experimental study on the effect of number and aperture of holes on mechanical properties of marble. Rock Soil Mech. 2017, 38, 41–50. [Google Scholar]
  19. Gong, F.; Wu, W.; Li, T. Experimental simulation of slab crack failure in deep rectangular hard rock tunnel. Rock Soil Mech. 2019, 40, 2085–2098. [Google Scholar]
  20. Zhu, T.; Jing, H.; Su, H. Experimental study on mechanical properties of sandstone with double circular holes under uniaxial compression. J. Geotech. Eng. 2015, 37, 1047–1056. [Google Scholar]
  21. Li, D.; Li, X.; Li, C. Experimental study on uniaxial compression mechanical properties of sandstone with double circular orifices. J. Rock Mech. Eng. 2015, 37, 1047–1056. [Google Scholar]
  22. Li, C.; Li, X.; Li, D. Particle flow analysis of failure characteristics of marble with orifice. J. Eng. Sci. 2017, 39, 1791–1801. [Google Scholar]
  23. Xie, L.; Zhu, W.; Wang, S. Parallel calculation and analysis of 3d fracture process of rock sample with orifice. J. Geotech. Eng. 2011, 33, 1447–1455. [Google Scholar]
  24. Wu, T.; Zhou, Y.; Wang, L. Mesoscopic study of pore-gap interaction mechanism under uniaxial compression. Rock Soil Mech. 2018, 39, 463–472. [Google Scholar]
  25. Zhu, T.; Jing, H.; Su, H. Experimental study on mechanical properties of porous and fissure combined defective sandstone. J. Coal 2015, 40, 1518–1525. [Google Scholar]
  26. Zhang, K.; Liu, T.; Yang, H. Physical test and numerical simulation of mechanical properties of rock mass with fissure after grouting. Hydrogeol. Eng. Geol. 2019, 46, 79–85+110. [Google Scholar]
  27. Fu, J.; Zhu, W. Study on fracture process of new material fracture body with 3d built-in fracture surface. J. Cent. South Univ. (Nat. Sci.) 2014, 45, 3257–3263. [Google Scholar]
  28. Li, S.; Yang, L.; Li, M. Effect of inclination Angle of 3d built-in fracture on tensile mechanical properties and fracture characteristics of rock-like materials. J. Rock Mech. Eng. 2009, 28, 281–290. [Google Scholar]
  29. Liu, J.; Yang, C.; Mao, H. Study on crack propagation and evolution of shale based on CT image processing. J. Zhejiang Univ. Technol. 2015, 43, 66–72. [Google Scholar]
  30. Guo, Y.; Zhu, W. Experimental study on three-dimensional propagation of embedded elliptic crack under compression and shear conditions. J. Solid Mech. 2011, 32, 64–73. [Google Scholar]
  31. Adams, M.; Sines, G. Crack extension from flaws in a brittle material subjected to compression. Tectonophysics 1978, 49, 97–118. [Google Scholar] [CrossRef]
  32. Dyskin, A.V.; Jewell, R.J.; Joer, H. Experiments on 3-D crack growth in uniaxial compression. Int. J. Fract. 1994, 65, R77–R83. [Google Scholar] [CrossRef]
  33. Ju, Y.; Xie, H.; Zheng, Z.; Lu, J.; Mao, L.; Gao, F.; Peng, R. Visualization of the complex structure and stress field inside rock by means of 3D printing technology. Chin. Sci. Bull. 2014, 59, 5354–5365. [Google Scholar] [CrossRef]
  34. Liu, H. Crack Propagation Law and Application of 3D Printing Technology in Prefabricated Cracked Rock-like Materials; China University of Mining and Technology: Beijing, China, 2018. [Google Scholar]
  35. Wang, H.; Zhang, J.; Ren, R. Implementation of 3D deep buried crack 3D-ILC based on electromagnetic field-medium damage. J. Geotech. Eng. 2018, 15, 1–13. [Google Scholar]
  36. Li, S.; Zhong, F. A review of 3d fracture research. Int. Seism. Dyn. 1993, 34, 1–4. [Google Scholar]
  37. Roesler, W.J.; Vandenbark, G.R.; Hanson, R.W. Cyclic AMP and the induction of eukaryotic gene transcription. J. Biol. Chem. 1988, 263, 9063–9066. [Google Scholar] [CrossRef]
  38. Knauss, K.G.; Wolery, T.J. Muscovite dissolution kinetics as a function of pH and time at 70 °C. Geochim. Cosmochim. Acta 1989, 53, 1493–1501. [Google Scholar] [CrossRef]
  39. Sommer, R.; Tautz, D. Minimal homology requirements for PCR primers. Nucl. Acids Res. 1989, 17, 6749. [Google Scholar] [CrossRef] [Green Version]
  40. Zheng, J.Y.; Deng, G.D.; Chen, Y.J. Experimental investigation on dynamic response and fracture characteristics of discrete multilayered thick-walled explosion containment vessels. Explos. Shock Waves 2005, 25, 506–511. [Google Scholar]
  41. Wang, H.; Le, C.; Tang, L. Based on 3D-ILC three point bending fracture characteristics of brittle solid with horizontal internal crack. Rock Soil Mech. 2021, 42, 2773–2784. [Google Scholar]
  42. Friedel, R.; Norman, D.A. The Psychology of Everyday Things. Technol. Cult. 1990, 31, 554. [Google Scholar] [CrossRef]
Figure 4. Arrangement of stress birefringence effect test.
Figure 4. Arrangement of stress birefringence effect test.
Materials 15 01406 g004
Figure 6. Stress distribution around the orifice.
Figure 6. Stress distribution around the orifice.
Materials 15 01406 g006
Figure 10. Variation law of angular function.
Figure 10. Variation law of angular function.
Materials 15 01406 g010
Figure 12. The diagramming of step cracks and the causes of confluence. (a) Step crack morphology; (b) Force of uniform step; (c) Causes of confluence.
Figure 12. The diagramming of step cracks and the causes of confluence. (a) Step crack morphology; (b) Force of uniform step; (c) Causes of confluence.
Materials 15 01406 g012
Figure 13. Ideal reverse binary tree pattern.
Figure 13. Ideal reverse binary tree pattern.
Materials 15 01406 g013
Figure 15. Model mesh division. (a) Model 1: Complete orifice model; (b) Model 2: Main crack initiation model; (c) Model 3: Model with internal cracks. A and B are the start and end point of the crack front.
Figure 15. Model mesh division. (a) Model 1: Complete orifice model; (b) Model 2: Main crack initiation model; (c) Model 3: Model with internal cracks. A and B are the start and end point of the crack front.
Materials 15 01406 g015
Figure 16. Distribution of maximum tensile stress.
Figure 16. Distribution of maximum tensile stress.
Materials 15 01406 g016
Figure 17. Main crack propagation process.
Figure 17. Main crack propagation process.
Materials 15 01406 g017
Figure 18. Variation of stress intensity factor at the tip of the main crack. A and B are the start and end point of the crack front.
Figure 18. Variation of stress intensity factor at the tip of the main crack. A and B are the start and end point of the crack front.
Materials 15 01406 g018
Figure 19. Torsional shear distribution. (a) Model 1; (b) Model 2.
Figure 19. Torsional shear distribution. (a) Model 1; (b) Model 2.
Materials 15 01406 g019
Table 1. Research scheme.
Table 1. Research scheme.
SamplesSchematic DiagramNumberBuried Depths/cmWhether to CrushPurpose
Samples without internal cracks Materials 15 01406 i001A0/YesObserve crack propagation
C0-1NoObserve main cracks
C0-2NoObserve remote cracks
Samples with internal cracks Materials 15 01406 i002B0.50.5YesObserve crack propagation
C0.5-10.5NoObserve main cracks
C0.5-20.5NoObserve remote cracks
Materials 15 01406 i003B11YesObserve crack propagation
C1-11NoObserve main cracks
C1-21NoObserve remote cracks
Materials 15 01406 i004B1.51.5YesObserve crack propagation
C1.5-11.5NoObserve main cracks
C1.5-21.5NoObserve remote cracks
Materials 15 01406 i005B22YesObserve crack propagation
C2-12NoObserve main cracks
C2-22NoObserve remote cracks
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Le, C.; Ren, X.; Wang, H.; Yu, S. Experimental and Numerical Study on the Failure Characteristics of Brittle Solids with a Circular Hole and Internal Cracks. Materials 2022, 15, 1406. https://doi.org/10.3390/ma15041406

AMA Style

Le C, Ren X, Wang H, Yu S. Experimental and Numerical Study on the Failure Characteristics of Brittle Solids with a Circular Hole and Internal Cracks. Materials. 2022; 15(4):1406. https://doi.org/10.3390/ma15041406

Chicago/Turabian Style

Le, Chengjun, Xuhua Ren, Haijun Wang, and Shuyang Yu. 2022. "Experimental and Numerical Study on the Failure Characteristics of Brittle Solids with a Circular Hole and Internal Cracks" Materials 15, no. 4: 1406. https://doi.org/10.3390/ma15041406

APA Style

Le, C., Ren, X., Wang, H., & Yu, S. (2022). Experimental and Numerical Study on the Failure Characteristics of Brittle Solids with a Circular Hole and Internal Cracks. Materials, 15(4), 1406. https://doi.org/10.3390/ma15041406

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