## 2 Methods  
### Method  
Fracture mechanics were widely used in crack analysis for their high accuracy, and the development of computer technology offers the convenience of calculating. Thus fracture mechanics are commonly used in pavement engineering. Using the finite element method for digital modeling can calculate the behavior of cracks. Responses of cracks could be summarized by modifying mechanical parameters of finite element model such as working conditions, structure combination, material composition. Fracture mechanics divide cracks into three species (open crack--type I, sliding crack--type II, tear-open crack--type III) and stress intensity factors (\(K_{I}\), \(K_{II}\), \(K_{III}\)) were used to describe the tendency of cracking [12]. Linear superposition of three kinds of cracks can be used to describe the stress field at the tip of cracks. Normal stress and shear stress of the tip of cracks can be calculated using the following equations. \[\sigma_{x}=\frac{K_{I}}{\sqrt{2\pi r}}cos\frac{\theta}{2}\bigg{(}1-sin\frac{ \theta}{2}sin\frac{3\theta}{2}\bigg{)}+\frac{K_{II}}{\sqrt{2\pi r}}sin\frac{ \theta}{2}\bigg{(}-2-cos\frac{\theta}{2}cos\frac{3\theta}{2}\bigg{)} \tag{1}\]  
\[\sigma_{y}=\frac{K_{I}}{\sqrt{2\pi r}}cos\frac{\theta}{2}\bigg{(}1+sin\frac{ \theta}{2}sin\frac{3\theta}{2}\bigg{)}+\frac{K_{II}}{\sqrt{2\pi r}}sin\frac{ \theta}{2}cos\frac{3\theta}{2} \tag{2}\]  
\[\tau_{xy}=\frac{K_{I}}{\sqrt{2\pi r}}sin\frac{\theta}{2}cos\frac{3\theta}{2} +\frac{K_{II}}{\sqrt{2\pi r}}cos\frac{\theta}{2}\bigg{(}1-sin\frac{\theta}{2 }sin\frac{3\theta}{2}\bigg{)} \tag{3}\]  
\[\tau_{zx}=-\frac{K_{III}}{\sqrt{2\pi r}}\sin\frac{\theta}{2} \tag{4}\]  
\[\tau_{yz}=\frac{K_{III}}{\sqrt{2\pi r}}cos\frac{\theta}{2} \tag{5}\]  
where \(r\) represents the distance from the calculation point to the crack tip. \(\theta\) represents the angle between the calculation point and the crack tip in the polar coordinate system and the polar axis. Analysis of the above formula shows that the key to determining the stress field at the crack tip is the stress intensity factor. These factors are mainly determined by the properties of external load, the shape of cracks, and the geometry of elastomer, which can be expressed by the equation below. \[K_{m}=\lim_{|r|\to 0}\sqrt{2\pi r}Z_{m}(r) \tag{6}\]  
where \(Z_{m}\) (\(r\)) is an analytic function related to the boundary conditions, which can be calculated using finite element method software ABAQUS. The stress intensity factor at the ultimate fracture of a material is called fracture toughness, usually expressed by \(K_{IC}\), which can be determined by experiments. Compared with the stress intensity factor, J-integral is more suitable for analyzing fracture problems under elastoplastic conditions [13]. The J-integral is used to solve the situation that a certain range of plastic zone appears at the tip of the crack under the elastoplastic condition, which makes the problem very complicated. The J-integral can describe the stress and strain field strength in the crack tip area, and it can be easily determined by experiment. Use stress intensity factor and J-integral can reveal the crack behavior in pavement engineering. A plane crack problem was assumed when calculating J-integral. The integral is performed around the crack tip, starting from the lower crack surface, and stopping counterclockwise to the crack's upper surface, thus forming an integral loop. This integral value obviously has nothing to do with the integral path, so J-integral can reflect the intensity of the stress field near the crack tip. Figure 1 is a schematic diagram of the calculation of J-integral. J-integral is defined as:  
\[J=\int_{\Gamma}\biggl{(}\omega dy-\overset{\rightarrow}{T}\cdot\frac{ \partial\overset{\rightarrow}{u}}{\partial x}ds\biggr{)} \tag{7}\]  
where \(\omega\) represents the strain energy density on the loop. \(\overset{\rightarrow}{T}\) is the stress vector at any point on the loop. \(\overset{\rightarrow}{u}\) is the displacement vector at any point on the loop.