The failure process of notched specimens with pre-existed cracks or hard granules under external tensile loading are simulated with RFPA2D(Rock Failure Process Analysis)code. The influence of different heterogeneity on the coalescence of cracks and the load-displacement curve is analyzed. The simulating results indicate that it is the local stress concentration and deformation localization due to the heterogeneity that causes the complex interaction of cracks and different patterns of propagation and coalescence.

Based on the concept of equivalent stiffness of rock masses，the estimation of joint stiffness in distinct element model of stiff blocks is discussed. The propagation of P-wave in one-dimensional rock masses is analyzed according to plane to plane contact model. The P-wave shape is presented. The effects of damping，soft interlayer and tension property on the propagation characteristics are studied.

The middle-isolated pier is one of the key structures of permanent shiplock of the Three Gorges Project on Yangtze River，China. The rockmass deformation of middle-isolated pier of the second shiplock under various operation conditions is studied. The numerical simulations using 2-D distinct element method and 3-D finite element method are applied to analyze rockmass deformation law. Simultaneously，the time series analysis and grey system theory are used to predict the tendency of rock displacements based on in-situ measured data. The comprehensive results can be summarized as follows. (1) The horizontal displacement perpendicular to axial line of shiplock varies with process of excavation，and is in steady state after excavation. During operation of shiplock，the displacement is slightly decreased. (2) The horizontal displacement parallel to axial line of shiplock is of tendency of developing to downstream. (3) The vertical displacement occurred as resilient deformation during excavation，and then subsidence with slightly various values takes place in different part of the middle-isolated pier under different operation cases.

The relation between the falling ability and damage of top coal is studied by damage mechanics. The model of damage parameter of top coal is found using meso-statistical damage mechanics. The damage parameter of top coal is regarded as the index of falling ability. Some factors，such as the mining depth，the strength of coal seam, initial crack and horizontal stress，which affect the falling ability of top coal，are discussed. The theoretical results are well agreeable with the experimental results.

Two-dimensional FEM，which can consider rheology of soil，is used to simulate the construction of excavation. Making use of this method，the rules of ground deformation in the bottom of foundation pit and round it are studied in details for different steps and different supporting conditions during excavation. The correlation between displacements and settlements behind supporting wall is also gained. The results can be used to consult in practice.

Bolt-grouting support is active method to combine bolt-spouting and grouting. In light of its support mechanism，a calculating model about its deformation is derived by means of action of bolt represented by volume-force and action of grouting regarded as improvement of properties of rock masses. The action of the support in preserving stability of surrounding rock masses is analyzed and its effect is proved being remarkable with using Kelvin-Voigt model and visco-elastic theory.

Tieshan tunnel of the Bapeng road is located in Daxiang of Shichuan province with the total length of 2 099 m. The tunnel passes through coal mined-out area K17-1，K21 and K24. The deformation，cracking，alteration and falling occur at the secondary lining at the station of K140+527～608 due to the influence of mined-out area，underground water and other factors，which do harm to the traffic safety，and must be cured. The stability of roof rocks of mined-out area is analysed according to engineering geological conditions and the reinforcement technology of mined-out area is presented.

By field observation in three longwall faces with different shallow overburden，the ground pressure behavior and the basic law of roof weighting are revealed. It is found that the strata breaking runs though whole roof，and there are big roof steps pushing in face during roof weighting. Shallow seam roof can not form stable roof block structure，and it will cause obvious roof pressure. It is also found that there are only two zones forming after strata collapse，the caving zone and the cracked zone. Speeding up face advancing is propitious to roof stability. Finally，the shallow seam definition for strata control is put forward by use of three key index，the key stratum，the thickness ratio of rock to soft overburden and the depth of coal seams.

With the help of theoretical analysis，numerical calculation and in-situ observations，the structure around the tailentry in soft rocks and its deformation and failure regularity in pillarless fully-mechanized sublevel caving are mainly introduced. The in-seam distressing entry layout is combined with concussion blasting technique，so that the tailentry in soft rocks can be effectively maintained. Meanwhile，the general distressing principle and practical distressing program are discussed. The theory and practice prove that this is a scientific strategy and an effective deci-sion-making to bring the tailentry in soft rocks under control.

A case-based reasoning approach to slope stability evaluation based on fuzzy analogy preferred ratio is presented，and a relative slope base case indexing model is set up. In this model，the fuzzy analogy preferred relationship between slope target case and slope base case is determined，and the similar series between the two cases with different affecting factor on slope stability is given. Determining the factor weights，the comprehensive similar series between the two cases is determined，and the most similar base case to the slope target case in the base cases of slope is found out. Finally，the slope stability of target case is evaluated. It is shown from examples that the result of stability evaluation of the slope is the same as the practical state，and the approach is simple，visual，practical，and convenient to use.

A new laboratory apparatus for coupled shear-flow of jointed rock was developed and used to investigate the coupled characteristic of a single rock fracture under shear and normal effective stresses. Based on shearing test of jointed marble sample under constant normal loads，the relationships among normal stress，shear stress，water pressure and flow quantity are presented, and the seepage charactristic of jointed rock is analysed. The relation equation between flow quantity and shear stress is deduced to fit the curves of testing data. By using an approximate regression equation，the relationship between the hydraulic conductivity and shear stress is estabilished eventually. The testing results reveal that the hydaulic conductivity of a sand-filled joint is reduced slightly with increasing shear stress，and the relationship between hydraulic conductivity and shear stress is linear approximately.

Following the development of resources exploitation and engineering construction，the excavation depth is increased. The special mechanic character of surrounding rockmass in deep engineering appears and it is the challenge to the traditional theory. The current research state of surrounding rockmass character in deep engineering and design method is summarized and its developing tendency is put forward.

The probability distribution model of discontinuities is established through the three dimensional network modeling. A computational program applying the model is developed. With the program，all the finite rock blocks can be found out with free face exposed on excavation face. The geometrical parameters of the rock block are presented and the areas of the block faces together with the volume and center of block are calculated with this program.

The expressing method of synthesis knowledge by frame plus rule and neural network plus mathematics model is presented. Through corresponding transform，the expressing method will become reasoning rule. According to the possible modes of stabilization of gangue dam slope，the hardness degree of rock，joined state and development degree of structure surface，the character of flabbiness surface of dam bed，the fragmentation state of dam shoulder and seepage are considered. The illation and evaluation on the stability of the 2nd gangue dam in Dexing copper mine are made.

On the basis of Kelvin-solution，the analytical solution of influence coefficients of elastic analysis is achieved in indirect three-dimensional boundary element method based on FSM. The numerical initial system of the method is set up. For elastic analysis of large region，the calculating speed is faster and the precision is higher in comparision with tranditional BEM.

The search angle is an imporment factor to be considered to improve the prodiction distance and accuracy in the gedogical prodiction with TSP-202 seismic system. With the theory and method of probability，tectonic geology and material mechanics，the search angle obtained by TSP-202 seismic system is analysed. An adjustment method is proposed for the search angle to improve the prediction accuraey，and applied in Xinluona tunnel with good result.

The fault effect of strata displacement in underground engineering is studied. Through some preferred index，the preferred fault can be determined and the intensity of special strata displacement is forecasted by preferred index method. Finally，the cause of preferred fault effect of strata displacement is analyzed with preferred plane theory .

Under the condition of pumping water from the deep well point，the viscoelastic flow model in deep excavation is proposed to consider the three-dimensional deformation of the viscoelastic compaction of the soil skeleton of aquifer due to pumping and coupling influence of soil-water. Furthermore，on the basis of the model，numerical simulation with FEM is used. In order to improve computing precision，unsaturated zone and unconfined flow are taken into consideration. In addition，an object-oriented FEM software based on Microsoft Windows 95/98/NT is developed. The basic program frame and main codes are given.

According to the distribution regularity of axial force and frictional resisting force of cement treated composite foundation，the force acting on soft soil-cement soil pile is considered simply as two kinds of friction forces uniformly and triangularly acting along the pile. This result is well fitted to Mindlin solution for the given cases. As a result，a simple calculation method is put forward for engineering design estimation beforehand.