Study on Mechanics Parameters of Unloaded Rock Mass about the Underground Powerhouse of Jin-ping Hydroelectric Station

After the excavation of Jinping underground cavern, a strong expansion appears along the unloading direction of the rock mass, mainly in the type of tensile rupture, accompanied by shear destruction, unloading resulted in significant deterioration of mechanical properties of rock. Based on the in-site investigation of rock mass structure, via analyzing the acoustic testing data, we identify the unloading range of the side walls and the division of rock types, and carry out with the solution of rock mechanical parameters about different unloading zone, providing geological foundation for the supporting design of the following design of the side walls, at the same time, providing reference for the selection of mechanical parameters of other underground excavation engineering with similar geological conditions.

Effect of Rock Mass Structure and Block Size on the Slope Stability-Physical Modeling and Discrete Element Simulation

This paper studies the stability of jointed rock slopes by using our improved three-dimensional discrete element methods (DEM) and physical modeling. Results show that the DEM can simulate all failure modes of rock slopes with different joint configurations. The stress in each rock block is not homogeneous and blocks rotate in failure development. Failure modes depend on the configuration of joints. Toppling failure is observed for the slope with straight joints and sliding failure is observed for the slope with staged joints. The DEM results are also compared with those of limit equilibrium method (LEM). Without considering the joints in rock masses, the LEM predicts much higher factor of safety than physical modeling and DEM. The failure mode and factor of safety predicted by the DEM are in good agreement with laboratory tests for any jointed rock slope.

Mechanical behavior of rock joints: influence of joint roughness on its closure and shear behavior

En esta tesis se estudia la fricción en juntas rocosas utilizando el Método de Elementos Discretos (DEM). En particular, se estudia la influencia de la rugosidad de las superficies de la junta, la elasticidad, la fractura, y el nivel de carga de compresión sobre el comportamiento de cierre y de cizalla de las juntas rocosas. Por primera vez la rugosidad de las juntas considerada como auto-afín es utilizada para estudiar la fricción de juntas rocosas, la rugosidad se describe mediante tres parámetros: el exponente de rugosidad, la longitud de correlación auto-afín y la varianza de alturas. Mediante un algoritmo de computadora basado en métodos espectrales, ocho superficies autoafines isotrópicas con diferente rugosidad fueron creadas. Posteriormente, las ocho superficies fueron utilizadas como moldes para generar las juntas utilizando elementos discretos. Antes de realizar las simulaciones de compresión y cizallaura, se calibraron las propiedades elásticas y de fractura (criterio de fractura elíptico basado en esfuerzos) de las juntas numéricas a los datos experimentales (obtenidos previamente) de unas muestras de mortero mediante la utilización de un volumen elemental representativo (REV). Una vez que las propiedades mecánicas de las juntas se obtuvieron mediante la calibración del REV, se realizaron las pruebas de cierre (prueba de compresión) de las ocho juntas DEM. Se utilizaron dos niveles de esfuerzo de compresión para las pruebas de cierre: 14 MPa y 21 MPa. Después, las ocho juntas DEM fueron cizalladas en dos direcciones mutuamente perpendiculares. Para cada dirección de cizalla y cada nivel de esfuerzo de compresión (14 y 21 MPa), las juntas fueron cizalladas usando uno de los tres modelos mecánicos siguientes: 1) un modelo rígido, en el que las juntas no se pueden deformar, excepto en su superficie, 2) un modelo puramente elástico, en el que las juntas se pueden deformar en todo su volumen y 3) un modelo elástico con fractura en el que las juntas se pueden deformar en su volumen y, si el esfuerzo sobre las uniones entre partículas excede cierto nivel de esfuerzo máximo, las uniones se rompen de una manera irreversible. El uso de estos tres modelos mecánicos nos permitirá estudiar de manera sistemática: la influencia de la rugosidad (modelo rígido), la influencia de la elasticidad y rugosidad (modelo puramente elástico) y, finalmente, el efecto combinado de la rugosidad de las juntas, la elasticidad y la fractura (modelo elástico con fractura). El estudio de los resultados obtenidos de las simulaciones DEM es seguido por una análisis energético el cual permite estudiar la evolución de los diferentes tipos de energía en función del desplazamiento de cizalla: energía elástica almacenada en el sistema, energía de fricción entre elementos discretos, el trabajo relacionado con la dilatación de la junta y la energía disipada por el amortiguamiento interno del DEM.

锦屏电站地下厂房卸荷岩体参数研究

It is a basic work to ascertain the parameters of rock mass for evaluation about stability of the engineering. Anisotropism、inhomogeneity and discontinuity characters of the rock mass arise from the existing of the structural plane. Subjected to water、weathering effect、off-loading, mechanical characters of the rock mass are greatly different from rock itself, Determining mechanical parameters of the rock mass becomes so difficult because of structure effect、dimension effect、rheological character, ‘Can’t give a proper parameter’ becomes one of big problems for theoretic analysis and numerical simulation. With the increment of project scale, appraising the project rock mass and ascertaining the parameters of rock mass becomes more and more important and strict. Consequently, researching the parameters of rock mass has important theoretical significance and actual meaning. The Jin-ping hydroelectric station is the first highest hyperbolic arch dam in the world under construction, the height of the dam is about 305m, it is the biggest hydroelectric station at lower reaches of Yalong river. The length of underground factory building is 204.52m, the total height of it is 68.83m, the maximum of span clearance is 28.90m. Large-scale excavation in the underground factory of Jin-ping hydroelectric station has brought many kinds of destructive phenomenon, such as relaxation、spilling, providing a precious chance for study of unloading parameter about rock mass. As we all know, Southwest is the most important hydroelectric power base in China, the construction of the hydroelectric station mostly concentrate at high mountain and gorge area, basically and importantly, we must be familiar with the physical and mechanical character of the rock mass to guarantee to exploit safely、efficiently、quickly, in other words, we must understand the strength and deformation character of the rock mass. Based on enough fieldwork of geological investigation, we study the parameter of unloading rock mass on condition that we obtain abundant information, which is not only important for the construction of Jin-ping hydroelectric station, but also for the construction of other big hydroelectric station similar with Jin-ping. This paper adopt geological analysis、test data analysis、experience analysis、theory research and Artificial Neural Networks （ANN） brainpower analysis to evaluate the mechanical parameter, the major production is as follows: (1)Through the excavation of upper 5-layer of the underground powerhouse and the statistical classification of the main joints fractures exposed, We believe that there are three sets of joints, the first group is lay fracture, the second group and the fourth group are steep fracture. These provide a strong foundation for the following calculation of and analysis; (2)According to the in-situ measurement about sound wave velocity、displacement and anchor stress, we analyses the effects of rock unloading effect，the results show a obvious time-related character and localization features of rock deformation. We determine the depth of excavation unloading of underground factory wall based on this. Determining the rock mass parameters according to the measurement about sound wave velocity with characters of low- disturbing、dynamic on the spot, the result can really reflect the original state, this chapter approximately the mechanical parameters about rock mass at each unloading area; (3)Based on Hoek-Brown experienced formula with geological strength index GSI and RMR method to evaluate the mechanical parameters of different degree weathering and unloading rock mass about underground factory, Both of evaluation result are more satisfied; (4)From the perspective of far-field stress, based on the stress field distribution ideas of two-crack at any load conditions proposed by Fazil Erdogan (1962)，using the strain energy density factor criterion (S criterion) proposed by Xue changming(1972)，we establish the corresponding relationship between far-field stress and crack tip stress field, derive the integrated intensity criterion formula under the conditions of pure tensile stress among two line coplanar intermittent jointed rock，and establish the corresponding intensity criterion for the exploratory attempt; (5)With artificial neural network, the paper focuses on the mechanical parameters of rock mass that we concerned about and the whole process of prediction of deformation parameters, discusses the prospect of applying in assessment about the parameters of rock mass，and rely on the catalog information of underground powerhouse of Jinping I Hydropower Station, identifying the rock mechanics parameters intellectually，discusses the sample selection, network design, values of basic parameters and error analysis comprehensively. There is a certain significance for us to set up a set of parameters evaluation system，which is in construction of large-scale hydropower among a group of marble mass.

A new approach for semi-automatic rock mass joints recognition from 3D point clouds

Rock mass characterization requires a deep geometric understanding of the discontinuity sets affecting rock exposures. Recent advances in Light Detection and Ranging (LiDAR) instrumentation currently allow quick and accurate 3D data acquisition, yielding on the development of new methodologies for the automatic characterization of rock mass discontinuities. This paper presents a methodology for the identification and analysis of flat surfaces outcropping in a rocky slope using the 3D data obtained with LiDAR. This method identifies and defines the algebraic equations of the different planes of the rock slope surface by applying an analysis based on a neighbouring points coplanarity test, finding principal orientations by Kernel Density Estimation and identifying clusters by the Density-Based Scan Algorithm with Noise. Different sources of information —synthetic and 3D scanned data— were employed, performing a complete sensitivity analysis of the parameters in order to identify the optimal value of the variables of the proposed method. In addition, raw source files and obtained results are freely provided in order to allow to a more straightforward method comparison aiming to a more reproducible research.

Stochastic Structural Model of Rock and Soil Aggregates by Continuum-Based Discrete Element Method

This paper first presents a stochastic structural model to describe the random geometrical features of rock and soil aggregates. The stochastic structural model uses mixture ratio, rock size and rock shape to construct the microstructures of aggregates,and introduces two types of structural elements (block element and jointed element) and three types of material elements (rock element, soil element, and weaker jointed element)for this microstructure. Then, continuum-based discrete element method is used to study the deformation and failure mechanism of rock and soil aggregate through a series of loading tests. It is found that the stress-strain curve of rock and soil aggregates is nonlinear, and the failure is usually initialized from weaker jointed elements. Finally, some factors such as mixture ratio, rock size and rock shape are studied in detail. The numerical results are in good agreement with in situ test. Therefore, current model is effective for simulating the mechanical behaviors of rock and soil aggregates.

Influence of stochastic mesoscopic structure on macroscopic mechanical behavior of brittle material

A newly developed numerical code, MFPA(2D) (Material Failure Process Analysis), is applied to study the influence of stochastic mesoscopic structure on macroscopic mechanical behavior of rock-like materials. A set of uniaxial compression tests has been numerically studied with numerical specimens containing pre-existing crack-like flaw. The numerical results reveal the influence of random mesoscopic structure on failure process of brittle material, which indicates that the variation of failure mode is strongly sensitive to the local disorder feature of the specimen. And the patterns of the crack evolution in the specimens are very different from each other due to the random mesoscopic structure in material. The results give a good explanation for various kinds of fracture modes and peak strength variation observed in laboratory studies with specimens made from the same rock block being statistically homogenous in macro scale. In addition, the evolution of crack is more complicated in heterogeneous cases than in homogeneous cases.

小湾坝基卸荷岩体参数研究

The Xiao-wan Power Station is the second highest arch dam in the world under construction. The height of the dam is about 292m. Large-scale excavation in the dam foundation of Xiao-wan Power Station has brought intensive unloading phenomenon. We collected a large number of firsthand data on unloaded rock mass in dam foundation, which supplies a natural testing ground for researching unloaded rock mass after excavation. Detailed study was carried out on the parameters of unloaded rock mass in the dam foundation of Xiao-wan Power Station. The study is not only importance for the Xiao-wan Power Station, but also has important instruction significance to similar projects in the Southwest of China. In order to study the mechanical parameters of unloaded rock mass, large field and laboratory tests were carried out. The test results showed the size effect of the sample is obvious. The change of deformation modulus of rock is not obvious. However, the Poisson's ratio of rock is increased under unloaded condition, its value is even more than 0.5. The theoretical forecasted results is accordance to the field tests including sound wave data and deformation monitoring data, which shows the forecasted results were reasonable. The soften yield criterion was adopted in the thesis to study the characteristics of the brittle rock mass in order to simulate their brittle failures. Based on the study results above, the transform of the structural plane network model to the numerical one was carried out, which made it feasible to consider the influence of large amount joints on the mechanical characteristic of rock mass in the numerical analyses. Using a factor, the degree of the damage or strengthen of rock mass can be determined rapidly, which proposed a rapid and feasible method for the determination of the parameters of rock mass.

基于窗口精细结构描述的岩体质量评价及力学参数综合优选

This thesis bases on horizontal research project “The research about the fine structure and mechanical parameters of abutment jointed rock mass of high arch dam on Jinping Ⅰ Hydropower Station, Yalong River” and “The research about the fine structure and mechanical parameters of the columnar basalt rock mass on Baihetan Hydropower Station, Jinsha River”. A rounded system about the fine structure description and rock mass classification is established. This research mainly contains six aspects as follow: (1) Methods about fine structure description of the window rock mass; (2) The window rock mass classification about the fine structure; (3) Model test study of intermittent joints; (4) Window rock mass strength theory; (5) Numerical experimentations about window rock mass; (6) The multi-source fusion of mechanical parameters based on Bayes principle. Variation of intact rock strength and joint conditions with the weathering and relaxation degree is studied through the description of window rock mass. And four principal parameters: intact rock point load strength, integration degree of window rock mass, joint conditions, and groundwater condition is selected to assess the window rock mass. Window rock mass is classified into three types using the results of window rock mass fine structure description combined with joints develop model. Scores about intact rock strength, integrality condition, divisional plane condition and groundwater conditions are given based on window rock mass fine structure description. Then quality evaluation about two different types of rock mass: general joint structure and columnar jointing structure are carried out to use this window rock mass classification system. Application results show that the window rock mass classification system is effective and applicable. Aimed at structural features of window structure of “the rock mass damaged by recessive fracture”, model tests and numerical models are designed about intermittent joints. By conducting model tests we get shear strength under different normal stress in integrated samples, through samples and intermittent joints samples. Also, the changing trends of shear strength in various connectivity rates are analyzed. We numerically simulate the entire process of direct shear tests by using PFC2D. In order to tally the stress-strain curve of numerical simulation with experimental tests about both integrated samples and through samples, we adjust mechanical factors between particles. Through adopting the same particle geometric parameter, the numerical sample of intermittent joints in different connective condition is re-built. At the same time, we endow the rock bridges and joints in testing samples with the fixed particle contacting parameters, and conduct a series of direct shear tests. Then the destructive process and mechanical parameters in both micro-prospective and macro-prospective are obtained. By synthesizing the results of numerical and sample tests and analyzing the evolutionary changes of stress and strain on intermittent joints plane, we conclude that the centralization of compressive stress on rock bridges increase the shear strength of it. We discuss the destructive mechanics of intermittent joints rock under direct shear condition, meanwhile, divide the whole shear process into five phases, which are elasticity phase, fracture initiation phase, peak value phase, after-peak phase and residual phase. In development of strength theory, the shear strength mechanisms of joint and rock bridge are analyzed respectively. In order to apply the deducted formulation conveniently in the real projects, a relationship between these formulations and Mohr-Coulomb hypothesis is built up. Some sets of numerical simulation methods, i.e. the distinct element method (UDEC) based on in-situ geology mapping are developed and introduced. The working methods about determining mechanical parameters of intact rock and joints in numerical model are studied. The operation process and analysis results are demonstrated detailed from the research on parameters of rock mass based on numerical test in the Jinping Ⅰ Hydropower Station and Baihetan Hydropower Station. By comparison，the advantages and disadvantages are discussed. Results about numerical simulation study show that we can get the shear strength mechanical parameters by changing the load conditions. The multi-source rock mass mechanical parameters can be fused by the Bayes theory, which are test value, empirical value and theoretical value. Then the value range and its confidence probability of different rock mass grade are induced and these data supports the reliability design.

复杂应力环境下断续节理岩体真三轴实验研究及其工程应用

The great deal of joints and faults , existing in the rock mass , are the leading cause of discontinuous rock mass. Structural planes not only destroy the integrality of rock mass, but also lead nonlinearity、heterogeneity、anisotropy and failure mode on mechanical properties of rock mass. Therefore the selection of strength and deformation parameters was very difficult. In practical rock mass engineering, equivalent parameters of rock mass were selected by the method of expert experience and engineering analogy. Based on the fine description of discontinuous joints in the type Ⅳ and Ⅴ rock mass and geological survey datum in situ, models was obtained by generalizing the structure of rock mass by the method of statistical analysis. Model intensity and deformation test were carried out on the true triaxial apparatus. Intermediate principle stress effect, anisotropy and dimension effect of discontinuous rock mass were considered in the model test. 3-D correction to Hoek-Brown empirical criterion was done by analysed the test datum. Detailed works were listed as follows: (1) The factors influenced intensity and deformation of discontinuous joints rock mass were the value of 、continuity, density and included angle of joints and anisotropy of joint plane. True triaxial intensity and deformation tests were carried out by considering above factors. The influence rule was obtained and corresponding relation formulary was established； (2) Based on the true triaxial tests under different stress path and load modes, we obtain intensity and deformation rule of rock mass； (3) Based on a great deal of true triaxial tests and other test datum, correction to the Hoek-Brown empirical criterion was done in the chapter 4. The intermediate principle stress was considered in the corrected formulary. It indicated that the formulary was applicable under a certain condition. In addition, the yield plane form of corrected Hoek-Brown empirical criterion under principle stress space was described in the paper. And the question of corner of yield plane was discussed； (4) Based on the single discontinuity theory, the three-dimensional intensity formulary of discontinuous joint rock mass was established. Correction to the intensity formulary was done considering intermediate principle stress effect. We may obtain the conclusion that the intensity of the discontinuous joint rock mass was influenced on compositive factors. They were 、 、continuity、internal frictional angle and cohesiveness of joint plane and rock； (5) The results of the true triaxial model test was applied into parameters evaluation of dam foundation rock mass of JinPing hydropower station. For there were abundant ophicalcite in the dam foundation, the interval of intensity and formation parameters influenced on continuity were determined based on test datum. (6) Especial mould for prismatic jointing model was designed. True triaxial intensity and deformation tests by Basalt with prismatic jointing were carried out. The influence of intermediate principle stress, stress path, anisotropy effect and dimensional effect to intensity and deformation was discussed in the chapter 6. The work of (3)、(4)、(6) was significative supplement and innovation to current test and theory.

岩体力学参数及稳定性能量判据研究

There are two major problems that have been concerned all the times, which are the mechanics characters of joint rock mass and the criterion for stability of engineering rock. Aim at the two problems, several works were conducted as follow: (1) Firstly, the mechanics characters of rock mass was studied by means of the Distinct Element Code. Subsequently, it was studied that the sensibility of joint surface roughness, strength of joint wall, joint stiffness ( i.e. tangential and normal stiffness) on the rock mass strength. (2) Based on the experimental rock mass classification methods of RMR and GSI, the program of “Parameters Calculation of the Rock Mass ” was developed. It has realized the rapid choice of rock mass parameters. (3) The concept of Representive Element Volume was induced based on the study of dimensional effect of rock mass. The Representive Element Volume of the horizontal and vertical pillar (ab. Two Pillars ) in the 2nd zone of Jinchuan mine were gained by the Geology Statistic Method and the Distinct Element Code. And then, the strength and deformatiom parameters of rock mass of the Two Pillars were obtained through numerical experiment. (4) From the confining depressure after thriaxial compression test of rock sample, it was concluded that the failure of rock is caused mainly by the lateral deformation and energy release happened during the confining depressure processure. The criterion of plastic energy catastrophe of rock engineering failure was proposed and validated. Subsquently, the stability of the horizontal pillar and Qianjiangping landslide in Three Gorges was judged by means of above-mentioned method. (5) Based on the fact there is a phenomenon of increasing energy concentration while the rock mass was compressed, rock information entropy (i.e. energy distribution entropy) was proposed. And it was revealed that there was change of energy distribution entropy while the rock mass was compressed to failure.

基于精细结构描述及数值试验的断续节理岩体力学参数智能选取

Evaluating the mechanical properties of rock masses is the base of rock engineering design and construction. It has great influence on the safety and cost of rock project. The recognition is inevitable consequence of new engineering activities in rock, including high-rise building, super bridge, complex underground installations, hydraulic project and etc. During the constructions, lots of engineering accidents happened, which bring great damage to people. According to the investigation, many failures are due to choosing improper mechanical properties. ‘Can’t give the proper properties’ becomes one of big problems for theoretic analysis and numerical simulation. Selecting the properties reasonably and effectively is very significant for the planning, design and construction of rock engineering works. A multiple method based on site investigation, theoretic analysis, model test, numerical test and back analysis by artificial neural network is conducted to determine and optimize the mechanical properties for engineering design. The following outcomes are obtained: (1) Mapping of the rock mass structure Detailed geological investigation is the soul of the fine structure description. Based on statistical window，geological sketch and digital photography，a new method for rock mass fine structure in-situ mapping is developed. It has already been taken into practice and received good comments in Baihetan Hydropower Station. (2) Theoretic analysis of rock mass containing intermittent joints The shear strength mechanisms of joint and rock bridge are analyzed respectively. And the multiple modes of failure on different stress condition are summarized and supplied. Then, through introducing deformation compatibility equation in normal direction, the direct shear strength formulation and compression shear strength formulation for coplanar intermittent joints, as well as compression shear strength formulation for ladderlike intermittent joints are deducted respectively. In order to apply the deducted formulation conveniently in the real projects, a relationship between these formulations and Mohr-Coulomb hypothesis is built up. (3) Model test of rock mass containing intermittent joints Model tests are adopted to study the mechanical mechanism of joints to rock masses. The failure modes of rock mass containing intermittent joints are summarized from the model test. Six typical failure modes are found in the test, and brittle failures are the main failure mode. The evolvement processes of shear stress, shear displacement, normal stress and normal displacement are monitored by using rigid servo test machine. And the deformation and failure character during the loading process is analyzed. According to the model test, the failure modes quite depend on the joint distribution, connectivity and stress states. According to the contrastive analysis of complete stress strain curve, different failure developing stages are found in the intact rock, across jointed rock mass and intermittent jointed rock mass. There are four typical stages in the stress strain curve of intact rock, namely shear contraction stage, linear elastic stage, failure stage and residual strength stage. There are three typical stages in the across jointed rock mass, namely linear elastic stage, transition zone and sliding failure stage. Correspondingly, five typical stages are found in the intermittent jointed rock mass, namely linear elastic stage, sliding of joint, steady growth of post-crack, joint coalescence failure, and residual strength. According to strength analysis, the failure envelopes of intact rock and across jointed rock mass are the upper bound and lower bound separately. The strength of intermittent jointed rock mass can be evaluated by reducing the bandwidth of the failure envelope with geo-mechanics analysis. (4) Numerical test of rock mass Two sets of methods, i.e. the distinct element method (DEC) based on in-situ geology mapping and the realistic failure process analysis (RFPA) based on high-definition digital imaging, are developed and introduced. The operation process and analysis results are demonstrated detailedly from the research on parameters of rock mass based on numerical test in the Jinping First Stage Hydropower Station and Baihetan Hydropower Station. By comparison，the advantages and disadvantages are discussed. Then the applicable fields are figured out respectively. (5) Intelligent evaluation based on artificial neural network (ANN) The characters of both ANN and parameter evaluation of rock mass are discussed and summarized. According to the investigations, ANN has a bright application future in the field of parameter evaluation of rock mass. Intelligent evaluation of mechanical parameters in the Jinping First Stage Hydropower Station is taken as an example to demonstrate the analysis process. The problems in five aspects, i. e. sample selection, network design, initial value selection, learning rate and expected error, are discussed detailedly.

复杂裂隙岩体统计断裂本构关系与变形参数研究

Rock mass is widely recognized as a kind of geologic body which consists of rock blocks and discontinuities. The deformation and failure of rock mass is not only determined by rock block，but also by discontinuity which is virtually more important. Mutual cutting and combination of discontinuities controlled mechanical property of rock mass. The complex cutting of discontinuities determine the intense anisotropy on mechanical property of rock mass，especially under the effect of ground stress. Engineering practice has show that the brittle failure of hard rock always occurs when its working stress is far lower than the yield strength and compressive strength，the failure always directly related to the fracture propagation of discontinuities. Fracture propagation of discontinuities is the virtue of hard rock’s failure. We can research the rock mass discontinuous mechanical properties precisely by the methods of statistical analysis of discontinuities and Fracture Mechanics. According to Superposition Principle in Fracture Mechanics，A Problem or C Problem could be chosen to research. Problem A mainly calculates the crack-tip stress field and displacement field on internal discontinuities by numerical method. Problem C calculate the crack-tip stress field and displacement field under the assumption of that the mainly rock mass stress field has been known. So the Problem C avoid the complex mutual interference of stress fields of discontinuities，which is called crack system problem in Fracture Mechanics. To solve Problem C, field test on stress field in the rock mass is needed. The linear Superposition of discontinuities strain energies are Scientific and Rational. The difference of Fracture Mechanics between rock mass and other materials can mostly expression as：other materials Fracture Mechanics mostly face the problem A，and can’t avoid multi-crack puzzle, while the Rock mass Fracture Mechanics answer to the Problem C. Problem C can avoid multi-discontinuities mutual interference puzzle via the ground stress test. On the basis of Problem C, Fracture Mechanics could be used conveniently in rock mass. The rock mass statistics fracture constitutive relations, which introduced in this article, are based on the Problem C and the Discontinuity Strain Energy linear superposition. This constitutive relation has several merits: first, it is physical constitutive relation rather than empirical; second, it is very fit to describe the rock mass anisotropy properties; third, it elaborates the exogenous factors such as ground stress. The rock mass statistics fracture constitutive relation is the available approach to answer to the physical, anisotropic and ground stress impacted rock mass problems. This article stand on the foundation of predecessor’s statistics fractures constitutive relation, and improved the discontinuity distributive function. This article had derived the limitation of negative exponential distribution in the course of regression analysis, and advocated to using the two parameter negative exponential distribution for instead. In order to solve the problems of two-dimension stability on engineering key cross-sectional view in rock mass, this article derived the rock mass planar flexibility tensor, and established rock mass two-dimension penetrate statistics fracture constitutive relation on the basis of penetrate fracture mechanics. Based on the crack tip plasticity research production of penetrate fracture, for example the Irwin plasticity equifinality crack, this article established the way to deal with the discontinuity stress singularity and plastic yielding problem at discontinuity tip. The research on deformation parameters is always the high light region of rock mass mechanics field. After the dam foundation excavation of XiaoWan hydroelectric power station, dam foundation rock mass upgrowthed a great deal of unload cracks, rock mass mechanical property gotten intricacy and strong anisotropy. The dam foundation rock mass mostly upgrowthed three group discontinuities： the decantation discontinuity, the steep pitch discontinuity, and the schistosity plane. Most of the discontinuities have got partial unload looseness. In accordance with ground stress field data, the dam foundation stress field greatly non-uniform, which felled under the great impaction of tectonic stress field, self-weight stress field, excavation geometric boundary condition, and excavation, unload. The discontinuity complexity and stress field heterogeneity, created the rock mass mechanical property of dam foundation intricacy and levity. The research on the rock mass mechanics, if not take every respected influencing factor into consideration as best as we can, major errors likely to be created. This article calculated the rock mass elastic modulus that after Xiao Wan hydroelectric power station dam foundation gutter excavation finished. The calculation region covered possession monolith of Xiao Wan concrete double-curvature arch dam. Different monolith were adopted the penetrate fracture statistics constitutive relation or bury fracture statistics constitutive relation selectively. Statistics fracture constitutive relation is fit for the intensity anisotropy and heterogeneity rock mass of Xiao Wan hydroelectric power station dam foundation. This article had contrastive analysis the statistics fracture constitutive relation result with the inclined plane load test actual measurement elastic modulus and RMR method estimated elastic modulus, and find that the three methods elastic modulus have got greatly comparability. So, the statistics fracture constitutive relations are qualified for trust. Generally speaking，this article had finished following works based on predecessors job: “Argumentation the C Problems of superposition principle in Fracture Mechanics, establish two-dimension penetrate statistics fracture constitutive relation of rock mass, argue the negative exponential distribution limitation and improve it, improve of the three-dimension berry statistics fracture constitutive relation of rock mass, discontinuity-tip plastic zone isoeffect calculation, calculate the rock mass elastic modulus on two-dimension cross-sectional view”. The whole research clue of this article inherited from the “statistics rock mass mechanics” of Wu Faquan（1992）.