Numerical Analysis of Stress Fields and Crack Growths in the Floor Strata of Coal Seam for Longwall Mining

This paper predicts the possibility of water inrush from a confined aquifer under the action of mining activities and water pressure. The study uses numerical analyses to evaluate stress redistribution and crack growth which result from coal extraction operations. Two models are presented in this study. By simplifying the distribution of the disturbed vertical stress on the coal seam and floor around a working face, a model is established to analyze the additional stresses in the floor strata induced by mining activities. And some distribution features of all the additional stress components are described. By using the superposition principle in fracture mechanics, another model is developed to analyze the crack growth in the floor strata under the action of disturbed stresses and water pressure. And the stress intensity factors at the crack tip are presented and the process of crack growth is obtained in the advancement of a working face. Because of discretizing only loading areas and crack surfaces, the present methods can obtain the accurate numerical results. Finally, some suggestions are made for preventing the water inrush from a confined aquifer.


Introduction
In the coalfields of North China, coal seams are underlain by Ordovician limestone, which is strongly karstified and highly permeable.In most cases, the limestone aquifer contains an abundant supply of water and is confined.The activities of mining these coal seams are threatened with the confined water from the Ordovician limestone.Water inrush accidents result in heavy losses and casualties [1,2].With coal mining being executed in deeper sites, the pressure of confined water becomes larger and this type of water inrushes may occur more frequently.
In order to avoid water inrushes and keep the mines safe, large-scale dewatering or depressurizing of the Ordovician aquifer was chosen.In practice, dewatering causes serious environmental problems such as surface subsidence and underground water pollution and reduces the amount of water available for supply in the area.To solve this dilemma, some researches were devoted to make full use of the uncontaminated karst water obtained from dewatering for water supply [3].Another measure to prevent water inrushes is to use the floor strata between the coal seam and the aquifer as a protective rock pillar.However, once the strength of the floor strata is difficult to resist the invasion of high-pressure water, the confined water floods into working faces.In some cases, grouting into the rock layers is used to improve the strength of floor strata to resist the invasion of confined water.By taking these measures, the water in the Ordovician limestone can be preserved to some extent and the coal seams can be extracted in a relatively safe environment.
As early as in 1979, the field investigations of the floor strata were executed in Fengfeng coalfield, China, where mining activities were threatened by confined aquifers [4].After that, more field investigations of floor strata were performed in different types of mines.In Jingxing coalfield, China, it was found that the joints in floor strata were impermeable before mining, and confined water passed through the joints and flooded into the mines after mining [5].This phenomenon may be explained by the joint growth induced by mining activities and water pressure from a confined aquifer.Additionally, the test of water injection in boreholes was used for understanding the permeability variations in floor strata induced by mining activities [6,7].These measured data were further used to determine the failure zones directly below the coal seam and above the confined aquifer.
In the past several decades, other types of analytical methods have been employed to study mining with water pressure for preventing water inrushes.Wang et al. [8] established both physical and numerical models to gain a better understanding into fracture initiation and propagation in floor strata during a mining process.Yang et al. [9] and Li et al. [10] used finite element methods to understand the mechanisms of water inrushes.Because different coal seams have different vulnerabilities for water inrushes, Wu et al. [11] assessed the differences in vulnerability between multiple coal seams.To evaluate a major water inflow hazard in a mine, Rafiqul et al. [12] used boundary element methods to model the mininginduced stress field.Zuo et al. [13] explored the spatial morphology of karst collapse and employed the discontinuous deformation analysis method to investigate karst collapse.
Although a great deal of field investigations and numerical analyses have been performed and some research results of water inrushes have been obtained, all the mechanisms of water inrushes are not known and the water inrushes are often reported.It is estimated that more than 80 percent of water inrushes may occur in the case of the floor rocks having faults or joints.This means that faults and joints in floor strata may induce water inrushes under the action of water pressure from aquifers and mining activities.Therefore, the water inrushes induced by faults or joints in longwall mining need to be further studied.
The principal objectives of this paper are to understand the stress fields in floor strata induced by coal mining and to predict the joint growths in the floor strata of the coal seam threatened by confined water in the advancement of a working face.From now on, the term crack will be used to describe a joint.In this paper, firstly, two numerical methods to be used are briefly introduced.Secondly, a mechanical model for simulating coal mining is established and the additional stress fields induced by coal mining are presented.Thirdly, another mechanical model is established for analyzing the crack problems under the action of water pressure and mining activities, and the stress intensity factors (SIFs) are obtained for the cracks at different positions of floor strata.Finally, the crack growths are predicted using SIF values and fracture criteria, and some suggestions are made for preventing the water inrushes into working faces from confined aquifers in floor strata.

Analytical Methods for Elastic Fields and
Crack Problems in a Layered Medium where  *  (, ) and  *  (, ) are, respectively, stresses and displacements of Yue's solution for the field point  due to the unit force along the  direction at the source point ;   () is the traction at the source point ; the integral domain  is the loading area.It should be noted that the subscript  is a dummy index.
The corresponding computer program LayerSmart3D was written in FORTRAN.The techniques adopted in Lay-erSmart3D primarily involve the discretization of the loading area into a finite number of quadrilateral elements.Values of the loads are inputted at the node points of the discretized area.Numerical verification of LayerSmart3D indicates that numerical solutions of high accuracy can be efficiently calculated for elastic fields induced by the distributed loads in a layered medium.

Yue's Solution
Based DDM for Analyzing Crack Problems in a Layered Medium.Xiao and Yue [16,17] developed a new displacement discontinuity method (DDM) for the analysis of crack problems in a layered medium of infinite extent.This approach is also based on Yue's solution [15] and the corresponding computer program LayerDDM3D was written in FORTRAN.
Assume that the crack surfaces consist of Γ + and Γ − .In the global coordinates (, , ), the discontinuous displacements of the crack surfaces are defined as Δ  ( Γ + ) =   ( Γ + ) −   ( Γ − ).  is a displacement along the direction  = , , .The traction on the crack surface is defined as   ( Γ + ).The basic equations of Yue's solution based DDM are as follows: where   ( Γ + ) is the unit outward normal of the positive side of the crack surface at the source point  Γ + and  *  ( Γ + , ) is a new kernel function which can be calculated by using the traction  *  of Yue's solution.In developing the DDM for (2), the nine-node continuous and discontinuous elements are used to discretize a crack surface.Kutt's numerical method is used to calculate the hypersingular integral in the integral equation.By solving the discretized equations of (2), the discontinuous displacements on the crack surface can be obtained.The numerical verification shows that the present results are in very good agreement with the existing ones.However, the proposed DDM is only suitable for the analysis of crack problems in layered media of infinite extent under the action of loadings on crack surfaces.Other types of crack problems in infinite media can be analyzed by using the two numerical methods together.

The Analytical Methods of Crack Problems in Layered
Media.In the following, the two numerical methods mentioned above, i.e., LayerSmart3D and LayerDDM3D, and the superposition principle in fracture mechanics are utilized to analyze crack problems in a layered medium.LayerSmart3D is firstly employed to obtain the stress fields of a layered medium without a crack under the action of distributed loadings within layered media.Using the superposition principle [18], the tractions, whose values are equal to the mininginduced stress and which have opposite directions, are then applied on the crack surfaces in the layered medium without the above distributed loadings.LayerDDM3D is finally employed to obtain the discontinuous displacements of the crack surfaces under the action of the above tractions.The SIF values at the crack tip can be calculated by using the discontinuous displacements of the crack surfaces and the relationship between SIFs and discontinuous displacements.working face in this mine was chosen for field measurements.Table 1 describes the roof and floor strata of the extracted coal seam.The coal seam is about 48 m away from the Ordovician limestone in which the water pressure ranges from 1.5 to 2 MPa.The working face is located at a depth of about 145 m and the longwall panel is 90 m long and 350 m wide.The coal seam is 1.5 m thick and the dip angle is about 1-3 ∘ .For the sake of simplicity, the rock strata are treated as horizontal media in developing a mechanical model.The working face was advanced along the strike direction of the coal seam.From August 1979 to February 1980, comprehensive measurements of floor deformation and failure were performed in the advancing process of the working face.In the field measurements, large amounts of data were obtained and will be further employed to analyze the water inrush from confined aquifers here.

The Mining-Disturbed Area on the Interface of Coal
Seam and Floor.150 m wide.The vertical stress on the interface of the coal seam and the floor strata, which are not disturbed by mining, can be calculated using the following formula: where  is the average unit weight of the overlying strata and ℎ is the depth of the extracted coal seam below the ground surface.
Mining activities induce the stress concentration on the coal pillars and floor around the working face.In the field investigations at No. 2710 working face, the disturbed area and stress concentration on the coal seams were obtained.Based on these data, we develop a model for the analysis of the influence of coal mining on floor strata.
The vertical stresses along different cross sections on the disturbed area are assumed by using the measured data.The concentration factors of the vertical stress ahead of and behind the working face along the cross section A 1 -A 2 are taken to be  1 = 3 and  2 = 1, respectively.And the stress concentration factors along the cross sections B 1 -B 2 and C 1 -C 2 are taken to be  1 = 3.On the floor strata near the roadways and the working face, the stresses are highly concentrated and the stress concentration factor is taken to be  3 = 3.5.The stress concentration factors along other cross sections can be obtained by using the methods of linear or quadratic interpolation.
Figure 2 illustrates the process of obtaining the additional vertical stress on the interface of coal seam and floor along the cross section A1-A2 shown in Figure 1.It can be found

Element Mesh of the Disturbed Area and Comparison of Numerical and Field
Results.In order to analyze the elastic fields using the software LayerSmart3D, the disturbed area ABCD is discretized into the mesh with 2610 nodes and 8037 elements, shown in Figure 4.The additional vertical stresses on the disturbed area ABCD (see Figure 3) are applied on the nodes of the element mesh.Figure 5 illustrates the distribution of the additional vertical stress on the interface of coal seam and floor.It can be found that the distribution of the additional vertical stress shown in Figure 5 is in good agreement with the one shown in Figure 3.
After presenting the mechanical parameters of floor strata in Table 1 and the additional vertical stress in Figure 5, the software LayerSmart3D is used to obtain the stress and displacement fields of floor strata.Table 2 lists the vertical displacements at some points in the floor strata from numerical analyses and in situ measurements.It can be found that the displacements from the two methods are much close to each other if the differences of the geological and mining conditions for two methods are considered.

Variations of Additional Stresses in
Floor Strata Induced by Mining Figures 6-8 show the variations of three normal stresses in the floor strata.In Figure 6, the rectangular area can be divided into three areas by considering the distribution features of   .Below the goaf area,   is tensile and has small values.Below the coal seam in the front of the working face,     is compressive and its minimum value is about -1.6 MPa.And below the coal pillar outside the goaf area,   is tensile and its maximum value is about 1.6 MPa.In Figure 7, the rectangular area can also be divided into three areas.Below the goaf area,   approaches zero.Below the coal seam in the front of the working face,   is tensile and its minimum value is about 1.6 MPa.And below the coal pillar outside the   goaf area,   is compressive and its minimum value is about 1.6 MPa.In Figure 8, the rectangular area can be divided into two areas.  is compressive below the goaf area and is tensile away from the goaf area.
Figures 9-11 show the variations of three shear stresses in the floor strata.In Figure 9,   has a maximum value of about -1.0 MPa at the point (, , ) = (94m, 46m, 15m) and   becomes small away from this point.In Figure 10,   is concentrated in the front of the working face and becomes small away from this zone and its maximum value is about 1.0 MPa.In Figure 11,   is concentrated below the coal pillar outside the goaf area and becomes small away from this zone.4, this plane is along the advancing direction of the working face and has a horizontal distance of 25 m from the coal pillar outside the goaf.For the sake of simplification, the stress variations are referred to as the ones with the depth increasing in the following analyses.From these figures, we can have the following observations:

Distribution of the Additional Stresses in the Plane
(i) At  = 0m, all the stresses have very small values because the position is away from the maximum loadings.In the first and second layers, the stresses vary obviously with depth.In the other layers, all the stresses are very small and vary weakly with depth.
(ii)   has small values in the first layer and the absolute values of   increase in the second layer.In the first layer,   increases at  = 100m whilst   decreases at the other  values.In the second layer,   is negative and decreases.At  = 40, 81m,   is positive and decreases.At  = 90m,   is positive in the first layer and negative in the other remaining layers and decreases.  is negative and decreases at  = 95m, and   is negative and increases at  = 100m.
(iii) At  = 0, 40m,   decreases in the first, second, and third layers and increases in the remaining two layers.The variation of   at  = 81, 90, 95, 100m is different from the one at  = 0, 40m.  is negative at  = 0, 40m whilst   is positive at  = 81, 90, 95, 100m.In the different layers,   first increases and then decreases for a given  value from  = 81 to 100 m.The absolute values of   first increase and then decrease.  is always positive for all the  values.For  = 0m,   increases.Except for  = 0m, in the first and second layers,   first increases and then decreases.And from the third layer to the fifth layer,   first increases and then decreases.

A Model of Analyzing Crack Problems in Floor Strata
As shown in Figure 18, a rectangular crack is located in floor strata and is subject to mining-disturbed stresses and water pressure.A global coordinate system Oxyz, which is completely the same as the one shown in Figure 1, is attached to the model to be analyzed and a local coordinate system         is attached on the crack surface.The rectangular crack is 2 m Mathematical Problems in Engineering wide (AB and CD sides) and 4 m long (AD and BC sides).Assume that the crack surfaces are parallel to the advancing direction of the working face and the crack center is located at the point (, , ) = ( 1 ,  2 , 46m).
Using the superposition principle in fracture mechanics [18], the disturbed stresses induced by mining activities can be transformed to the tractions on the crack surfaces.The disturbed stresses should be the summation of the horizontal normal stress  ℎ induced by rock weights and the additional stresses (  ,   ,   ) induced by mining activities.Additionally, the crack surfaces are subject to the water pressure   from confined aquifer.Thus, the crack surfaces are subjected to the following tractions:   = −  (4a) where   is taken to be 2 MPa.In (4b),  ℎ is calculated by the following formula: where ] and  are, respectively, the average values of Poisson's ratio and unit weight of rock strata and ℎ is the depth of the crack center below the ground surface.Take  = 25kN/m 2 , ] = 0.28, and ℎ = 192m.Thus,  ℎ ≈ 1.8667MPa.The crack surface is discretized into 435 nodes and 98 nine-node elements with 60 nine-node isoparametric elements, 34 discontinuous elements of type I, and four discontinuous elements of type II.The discontinuous displacements of the crack surfaces can be obtained using the element mesh, the tractions (  ,   ,   ) on the crack surface, and the computer code LayerDDM3D.

SIFs of Cracks under the Action of
Water Pressure and Disturbed Stresses 5.1.General.The stress intensity factors (SIF) values at the crack tip can be calculated by using the discontinuous displacements of crack surfaces.When the crack is not located at the interfaces of multilayered solids, the crack tip field singularities and angular distributions are the same as those in a homogeneous medium and the material parameters exert an influence on the magnitudes of the SIFs.Based on the relationship of displacements and the SIFs, the formulae of the SIFs along the crack side BC are as follows: where (, , ) is the spherical coordinates located at the crack front and  is an angle in the coordinate plane normal to the crack line of crack front.
The following two cases are analyzed: Case 1: a rectangular crack is located in the coordinate plane Oxz, i.e.,  2 = 0m; Case 2: a rectangular crack is located in the plane  = 40m, i.e.,  2 = 40m.

Case 1: The Crack Is Located in the Coordinate Plane Oxz.
In this case, the shear stresses   and   in the coordinate plane Oxz are equal to zero because of the symmetry of the model.As a result,  II and  III are equal to zero. Figure 19 illustrates the variations of  I for the crack at different positions along the advancing direction of the working face.
For the crack only subject to water pressure, the maximum value of  I is 0.1134 MPa⋅m 1/2 .However, for the crack subject to water pressure and disturbed stresses, the maximum value of  I is 0.1351 MPa ⋅ m 1/2 at  1 = 0m.With  1 increasing,  I increases, arrives at maximum values at  1 = 40, 50m, and then decreases.For  1 = 80m and   ≥ 0.5m,  I is equal to zero and  I along the whole BC side is equal to zero at  1 ≥ 90m.These findings mean that the crack below the coal seam is closed and the crack below the goaf is opened under the action of disturbed stresses and water pressure.When the crack is located at some horizontal distance of 40-50m after the coal wall of the working face, there are the maximum opening displacements.In this case, the crack becomes easy to open and induces the uplifting of the confined water.

Case 2:
The Crack Is Located in the Plane  = 40m, i.e.,  2 = 40m.In this case, the crack is close to the coal pillars outside the working face and has a horizontal distance of 5 m from the coal pillar.As discussed in Section 3, the shear stresses induced by mining exist at this position and as a result the crack has three types of fracture modes.show the variations of  I ,  II , and  III along the crack front BC, respectively.In Figure 20, the  I values for  1 = 0, 40m are very close to the ones of the crack only subject to water pressure whilst x  (m)  arrive at the minimum values at  1 = 50m, and then increase.For  1 = 110 − 150m, as the  1 values increase, the  II values first increase, arrive at the maximum values at  1 = 130m, and then decrease.It is believed that the  II values tend to zero when the absolute values of  1 tend to an infinite value.
In Figure 22, the  III values for  1 = 0m are positive.As the  1 values increase, the  III values decrease along the most part of the crack side BC, become negative values, and arrive at the minimum values at  1 = 90m.Then, the  III values increase along the most part of the crack side BC.It is believed that the  III values tend to zero when the absolute values of  1 tend to an infinite value.

Analysis of Crack Growths
Induced by Mining Sih [20].The S-criterion was developed on the basis of the strain energy density concept.It states that local instability is assumed to occur when the local minimum energy factor  min reaches a critical value  cr .Using the fracture criterion and the SIF values in Section 5.3, the crack growth for Case 2 will be discussed in the following.The strain energy density factor  of a three-dimensional crack can be defined as [21]  () =  11 () where where  is the shear modulus of elasticity, ] is Poisson's ratio, and  is an angle in the coordinate plane normal to the crack line of crack front.
Using the SIF values presented in Section 5.3 and (7), the local minimum energy factor  min can be calculated.Figure 23 illustrates the variations of  min values with the distance  1 .At  1 = 0, 90, 110, 130, 150m, the  min values are smaller than the ones at  1 = 30, 50, 60, 70m.From  1 = 0 to 90 m,  min first increases, arrives at the maximum value at  1 = 50m, and then decreases.From  1 = 90 to 150 m,  min first decreases, arrives at the minimum value at  1 = 110m, and then increases.It can be estimated that  min tends to zero as the floor strata are far from the coal wall of the working face.
The following conclusions can be drawn: (i) For  1 ≤ 90m,  min is the largest at  1 = 50m whilst for  1 > 90m,  min is the largest at  1 = 110m.This means that  min has two largest values before and after the coal pillar of the working face.Whether the crack grows or not depends on the  cr magnitude of the shale.
(ii) The maximum values of  min at  1 = 50m are larger than the ones at  1 = 110m.The crack growth becomes easy when the crack is located below the goaf.It has been found that the  I values for Case 1 are the largest at  1 = 40, 50m.This means that the growth of any cracks parallel to the coordinate plane Oxz is easy at some horizontal distance of 40-50m behind the coal wall of the working face.

Conclusions
Numerical methods and mechanical models have been proposed for analyzing the stress fields and crack growths in floor strata for longwall mining.The additional stresses in the floor strata induced by mining are calculated and the crack problems subject to mining-induced stresses and water pressure are investigated.The research results presented in this paper can be summarized as follows: (1) The additional stresses of the floor strata for longwall mining are described in detail and the distribution features of all the stress components with the depth below the floor and the distance from the working face are obtained.
(2) By analyzing the crack problems under the action of water pressure and disturbed stresses, the SIF values of the crack at different positions of floor strata are obtained.When the crack is located in the coordinate plane Oxz, only the mode I deformation of the crack exists.However, when the crack is not located in the coordinate plane Oxz, the mode I, II, and III deformations of the crack are coupled together.
(3) At some horizontal distance of 40-50m behind the working face, the vertical crack in the coordinate plane Oxz has the largest values of  I and the vertical crack, which is not located in the coordinate plane Oxz, has the largest values of  min .At these positions, the crack growths become easier.Therefore, the field measures should be taken in the advancement of the working face to prevent the water inrushes into the goaf behind working faces.
It is the authors' belief that the analytical methods presented in this paper can be a powerful numerical tool, which can apply to other types of cracks in the floor strata of the coal seam threatened by confined aquifers.Some of the related problems are currently under investigation by the authors for further understanding the complex mechanisms of water inrushes.

Data Availability
The [DATA TYPE] data used to support the findings of this study are included within the article.

3. 1 .
Geological and Mining Conditions of the Working Face.Jing et al. [4] investigated the deformation and failure of floor strata in mining the coal seam threatened by confined water in No. 2 Coal Mine of Fengfeng coalfield.No. 2710

Figure 1
illustrates the disturbed area ABCD on the interface between the mined coal seam and the floor.It is assumed that the domain out of the area ABCD is not disturbed by coal mining.According to the in situ investigations, the disturbed area ABCD is 150 m long and nterface Pos t on of the work ng face Vertical stress on the interface of coal seam and floor for the disturbed case nterface p 0 (b) Vertical stress on the interface of coal seam and floor for the undisturbed case Advanc ng d rect on nterface Pos t on of the work ng face (k 1 -1)p 0 -k 2 p 0 (c) Additional vertical stress on the interface of coal seam and floor for the disturbed case

Figure 2 :
Figure 2: Analytical method of additional vertical stress on the interface of coal seam and floor along the cross section A 1 -A 2 .
Additional vertical stress on the interface of coal seam and floor along the cross section A 1 -A 2 1 -1)p 0 (b) Additional vertical stress on the interface of coal seam and floor along the cross section B 1 -B 2 y L (k 1 -1)p 0 (k 1 -1)p 0 c 3 c 2 c 1 -k 2 p 0 c 3 c 2 c 1 (c) Additional vertical stress on the interface of coal seam and floor along the cross section C 1 -C 2

Figure 3 :
Figure 3: Additional vertical stress on the interface of coal seam and floor along different cross sections.

Figure 4 :Figure 5 :
Figure 4: Element mesh of the disturbed area ABCD, the projection abcd of the calculating area at  = 15m, and the plane at  = 20m for calculating the stresses.

Figure 7 :
Figure 7: Contours of the additional stress   at the depth  = 15m.

Figure 8 :
Figure 8: Contours of the additional stress   at the depth  = 15m.

Figure 9 :
Figure 9: Contours of the additional stress   at the depth  = 15m.

Figure 10 :
Figure 10: Contours of the additional stress   at the depth  = 15m.

Figure 11 :
Figure 11: Contours of the additional stress   at the depth  = 15m.
= 20m.Figures 12-17 illustrate the distribution of all the stress components in the plane  = 20m for different  and  values.As shown in Figure

Figure 12 :Figure 13 :
Figure 12: Variations of the additional stress   with depth at  = 20m.

Figure 14 :Figure 15 :
Figure 14: Variations of the additional stress   with depth at  = 20m.

Figure 16 :Figure 17 :
Figure 16: Variations of the additional stress   with depth at  = 20m.

1 Figure 18 :
Figure 18: Model of analyzing crack problems in floor strata under the action of distributed stresses and water pressure.

Figure 19 :
Figure 19: Variations of  I values at the crack side (BC) for Case 1.

Figure 20 :
Figure 20: Variations of  I values at the crack side (BC) for Case 2.

Figure 21 :Figure 22 :
Figure 21: Variations of  II values at the crack side (BC) for Case 2.

Figure 23 :
Figure 23: Variations of  min values at the crack side (BC) for Case 2.

Table 1 :
The roof and floor strata of the extracted coal seam in No. 2710 working face.

Table 2 :
Comparison of maximum vertical displacements from in situ measurements and numerical analyses along the advancing direction (unit: mm).
[19]Case 1: The Crack Is Located in the Coordinate PlaneOxz.As discussed in Section 5.2, the crack has only mode I deformation type in this case.The crack growth can be analyzed by using the fracture criterion  I =  IC .Chong et al.[19]investigated the fracture toughness  IC of layered rocks by using three-point-bend specimens and presented the  IC values of oil shale.Using an analogical method, the  IC value of the shale in Table1is estimated to be 0.74 MPa ⋅ m 1/2 .As shown in Figure19, the maximum  Imax values are about 0.2335 MPa ⋅ m 1/2 and  Imax <  IC .This means that this type of crack does not grow and the confined water does not uplift under the action of water pressure and disturbed stresses.
6.2.Case 2:The Crack Is Located in the Plane  = 40m, i.e.,  2 = 40m.As discussed in Section 5.3, the mode I, II, and III deformations of the crack are coupled together.Growth of a crack in elastic solids subject to complex stress states can be assessed using the so-called S-criterion proposed by 2 I + 2 12  I  II +  22 ()  2