Recent advances of digitization in rock mechanics and rock engineering

In recent years, there are growing demands of representing rock mechanics and rock engineering in a digital format that can be easily managed, manipulated, analyzed and shared. The objective of this paper is to give a comprehensive review of the status quo and future trends of digitization in rock mechanics and rock engineering. Research topics essential to the process of digitization are firstly discussed, including data acquisition, data standardization, geological modeling, visualization and digital-numerical integration. New techniques that will play an important role in digitization process but require further improvement are then briefly proposed. Finally, achievements of present methods and techniques for digitization in substantial rock mechanics and rock engineering are presented.

评价饱和砂土液化过程中小应变到大应变的本构模型(英文)

在RambergOsgood模型的基础上,建立了一个实用的能描述饱和砂土从液化前小应变到初始液化后大应变范围的非线性本构模型。该模型主要有三个特点:①采用作者曾提出的移动相态转换线(MPTL)和移动临界状态线(MCSL)概念来较为合理地描述变形过程中的有效应力路径和临界状态;②采用一个应力软化模型来描述剪胀和定义超静孔压变化引起的骨架曲线的软化与硬化;③基于笔者等提出的液化后大变形理论来描述初始液化后产生的大剪切变形。文中通过对室内循环剪切试验结果的模拟,验证了模型的预测能力。该模型适用于水平饱和砂土地基的有效应力地震反应分析。

Seismic responses of high concrete face rockfill dams:A case study

Seismic responses of the Zipingpu concrete face rockfill dam were analyzed using the finite element method. The dynamic behavior of rockfill materials was modeled with a viscoelastic model and an empirical permanent strain model. The relevant parameters were obtained either by back analysis using the field observations or by reference to parameters of similar rockfill materials. The acceleration responses of the dam,the distribution of earthquake-induced settlement, and the gap propagation under the concrete slabs caused by the settlement of the dam were analyzed and compared with site investigations or relevant studies. The mechanism of failure of horizontal construction joints was also analyzed based on numerical results and site observations. Numerical results show that the input accelerations were considerably amplified near the top of the dam, and the strong shaking resulted in considerable settlement of the rockfill materials, with a maximum value exceeding 90 cm at the crest.As a result of the settlement of rockfill materials, the third-stage concrete slabs were separated from the cushion layer. The rotation of the cantilever slabs about the contacting regions, under the combined action of gravity and seismic inertial forces, led to the failure of the construction joints and tensile cracks appeared above the construction joints. The effectiveness and limitations of the so-called equivalent linear method are also discussed.

Centrifugal and field studies on water infiltration characteristics below canals under wetting-drying-freezing-thawing cycles

Seepage is one of the main causes for the deformation and instability of canal slopes in Xinjiang,China.In this study,centrifugal model tests under wetting-drying(WD)and wetting-drying-freezing-thawing(WDFT)cycles were performed to investigate the water infiltration characteristics below a canal.The results show that the shallow soil of the canal models was fully saturated in the wetting process.Compared with the canal model under the WD cycles,the canal model under the WDFT cycles had larger saturated areas and a higher degree of saturation below the canal top after each cycle,indicating that the freezing-thawing(FT)process in the WDFT cycles promoted the water infiltration behavior below the canal slope.The cracks on the surface of the canal model under the cyclic action of WDFT developed further and had a higher connectivity,which provided the conditions for slope instability from a transverse tensile crack running through the canal top.On this basis,a field test was conducted to understand the water infiltration distribution below a typical canal in Xinjiang,China,which also verified the accuracy of the centrifugal results.This study provides a preliminary basis for the maintenance and seepage treatment of canals in Xinjiang,China.

Influences affecting the soil-water characteristic curve

The soil-water characteristic curve (SWCC) is the primary partially saturated soil information as its behavior and properties can be derived from it. Although there have been many studies of unsaturated soils and the SWCC, there is still no combined constitutive model that can simulate soil characteristics accurately. In cases when hydraulic hysteresis is dominant (e.g. under cyclic loading) it is particularly important to use the SWCC. In the past decades, several mathematical expressions have been proposed to model the curve. There are various influences on the SWCC as a source of information, so the curves obtained from conventional tests often cannot be directly applied; and the mathematical expressions from one scenario cannot be used to simulate another situation. The effects of void ratio, initial water content, stress state and high suction were studied in this work revealing that water content and stress state are more important than the other effects; but that the influences tend to decrease when suction increases. The van Genuchten model was modified to simulate better the changes in the degree of saturation at low values of suction. Predictions were compared with experimental results to determine the simulation capability of the model.

Numerical simulation of dynamic response of operating metro tunnel induced by ground explosion

To evaluate the effects of possible ground explosion on a shallow-buried metro tunnel, this paper attempts to analyze the dynamic responses of the operating metro tunnel in soft soil, using a widely applied explicit dynamic nonlinear finite element software ANSYS/LS-DYNA. The blast induced wave propagation in the soil and the tunnel, and the von Mises effective stress and acceleration of the tunnel lining were presented, and the safety of the tunnel lining was evaluated based on the failure criterion. Besides, the parametric study of the soil was also carried out. The numerical results indicate that the upper part of the tunnel lining cross-section with directions ranging from 0° to 22.5° and horizontal distances 0 to 7 m away from the explosive center are the vulnerable areas, and the metro tunnel might be safe when tunnel depth is more than 7 m and TNT charge on the ground is no more than 500 kg, and the selection of soil parameters should be paid more attentions to conduct a more precise analysis.

Numerical Analysis of Interaction Between Pile-Supported Pier and Bank Slope

Two and three-dimensional finite element analysis programs for pile-soil interaction are compiled. Duncan-Chang's Model is used. The construction sequence of the pier is modeled. The pile-soil interface element is used. The influence of the combination type of piles on the deformation of bank slope and pile behaviour is analyzed. Different designs of a pile-supported pier are compared thoroughly. Calculation results show that the stresses and displacements of the pile are directly related to the distance from the bank slope and the direction of inclination. An inclined prop pile set in the rear platform would remarkably reduce the stresses of piles and the displacement of the pier.

Research on shape optimization of CSG dams

The multi-objective optimization method was used for shape optimization of cement sand and gravel （CSG） dams in this study. The economic efficiency, the sensitivities of maximum horizontal displacement and maximum settlement of the dam to water level changes, the overall stability, and the overall strength security were taken into account during the optimization process. Three weight coefficient selection schemes were adopted to conduct shape optimization of a dam, and the case studies lead to the conclusion that both the upstream-and downstream dam slope ratios for the optimal cross-section equal 1：0.7, which is consistent with the empirically observed range of 1：0.6 to 1;0.8 for the upstream and downstream dam slope ratios of CSG dams. Therefore, the present study is of certain reference value for designing CSG dams.

Fly Ash/Paraffin Composite Phase Change Material Used to Treat Thermal and Mechanical Properties of Expansive Soil in Cold Regions

Phase change materials(PCMs)can store large amounts of energy in latent heat and release it during phase changes,which could be used to improve the freeze-thaw performance of soil.The composite phase change material was prepared with paraffin as the PCM and 8%Class C fly ash(CFA)as the supporting material.Laboratory tests were conducted to reveal the influence of phase change paraffin composite Class C fly ash(CFA-PCM)on the thermal properties,volume changes and mechanical properties of expansive soil.The results show that PCM failed to establish a good improvement effect due to leakage.CFA can effectively adsorb phase change materials,and the two have good compatibility.CFA-PCM reduces the volume change and strength attenuation of the soil,and 8 wt.%PCM is the optimal content.CFA-PCM turns the phase change latent heat down of the soil and improves its thermal stability.CFA-PCM makes the impact small of freeze-thaw on soil pore structure damage and improves soil volume change and mechanical properties on a macroscopic scale.In addition,CFA-8 wt.%PCM treated expansive soil has apparent advantages in resisting repeated freeze-thaw cycles,providing a reference for actual engineering design.

Blast-induced ground vibration prediction in granite quarries:An application of gene expression programming,ANFIS,and sine cosine algorithm optimized ANN

Blasting of rocks has intrinsic environmental impacts such as ground vibration,which can interfere with the safety of lives and property.Hence,accurate prediction of the environmental impacts of blasting is imperative as the empirical models are not accurate as evident in the literature.Therefore,there is need to consider some robust predictive models for accurate prediction results.Gene expression programming(GEP),adaptive neuro-fuzzy inference system(ANFIS),and sine cosine algorithm optimized artificial neural network(SCA-ANN)models are proposed for predicting the blast-initiated ground vibration in five granite quarries.The input parameters into the models are the distance from the point of blasting to the point of measurement(D),the weight of charge per delay(W),rock density(q),and the Schmidt rebound hardness(SRH)value while peak particle velocity(PPV)is the targeted output.100 datasets were used in developing the proposed models.The performance of the proposed models was examined using the coefficient of determination(R2)and error analysis.The R2 values obtained for the GEP,ANFIS,and SCA-ANN models are 0.989,0.997,and 0.999,respectively,while their errors are close to zero.The proposed models are compared with an empirical model and are found to outperform the empirical model.

Gradation equation of coarse-grained soil and its applicability

Gradation equation is one way to describe the gradation of coarse-grained soil conveniently,exactly and quantitatively.With the gradation equation,the influence of gradation on the mechanical behaviors of coarse-grained soil can be expressed quantitatively.A new gradation equation with a parameter is proposed.The basic properties and applicability of the new equation are studied.The results show that the proposed equation has the applicability to express coarse-grained soil gradation(CSG),and the range of the parameter β is found to be 0<β<1.The value ofbdetermines the gradation curve shape.If β>0.5,the gradation curve is sigmoidal,otherwise the gradation curve is hyperbolic.For well graded gradations,the parameter has the value of 0.13<β<1.Several CSGs used in domestic and foreign earth-rockfill dams are probed,and the value of the parameterbfalls in the range of 0.18 to 0.97.The investigation of the range of β is of value to guide the design for CSG of earth-rockfill dam.

Undrained monotonic and cyclic behavior of sand-ground rubber mixtures

In this study the stress–strain characteristics of sand-ground rubber mixtures are investigated in the sandlike zone,at different confining pressures,using hollow cylinder specimens subjected to torsional monotonic and cyclic loading.Under monotonic loading a mixture of sand-ground rubber with 10% and 25% rubber content show more contraction behaviour than that observed in a pure sand specimen.Phase transformation point in these mixtures are located on a larger shear strain.As expected,the shear strength of specimens decreases with increase of ground rubber content.However,with increasing of effective confining pressure,the loss in shear strength of the mixture is decreased.In addition,a mixture with 25% ground rubber shows a smaller loss in shear strength compared to a mixture with 10% ground rubber mixture.Under cyclic loading mixtures with 10% and 25% ground rubber have similar liquefaction resistance,especially at confining pressures of 110 k Pa and 260 k Pa.Therefore,by using of the mixture with 25% ground rubber,a larger volume of scrap tires could be recycled.The addition of ground rubber to sand would affect the shear strain variation and excess pore water pressure trends,and this effect was further intensified with increasing ground rubber percentage.

Evaluation of the impact of commodity price change on mine plan of underground mining

Fluctuations in commodity prices should influence mining operations to continually update and adjust their mine plans in order to capture additional value under new market conditions. One of the adjustments is the change in production sequencing. This paper seeks to present a method for quantifying the net present value(NPV) that may be directly attributed to the change in commodity prices. The evaluation is conducted across ten copper price scenarios. Discrete event simulation combined with mixed integer programming was used to attain a viable production strategy and to generate optimal mine plans. The analysis indicates that an increase in prices results in an increased in the NPV from$96.57M to $755.65M. In an environment where mining operations must be striving to gain as much value as possible from the rights to exploit a finite resource, it is not appropriate to keep operating under the same mine plan if commodity prices alter during the course of operations.

Rock borehole shear tests in dam foundation of Xiangjiaba hydropower station

Xiangjiaba hydropowcr station is one of the complicated geological conditions of its dam foundation, parameters of rock masses are very important issues. To cascade power stations on the Jinsha River, China. Due to the evaluating the rock mass quality and determining the mechanical address these issues, several groups of rock borehole shear tests （RBSTs） were conducted on the black mudstone in the dam foundation of Xiangjiaba hydropower station in the second construction phase. Forty three groups of shear strengths of black mudstone samples were obtained from RBSTs, and the shear strength parameters （c and f） were calculated using the least squares method. In addition, the limitations and merits of RBST employed in the Xiangjiaba hydropower station were discussed. Test results indicate that the shear strength parameters obtained from RBST have a good correlation with the results from sotmd wave test in borehole. It is believed that RBST has a good adaptability and applicability in geotechnical engineering.

Accurate detection method of high-density polyethylene fi lm leakage location based on secondary electric potential difference

High-density polyethylene(HDPE)film leakage location detection is frequently accomplished using the double-electrode technique.The electric potential and potential difference are the main physical parameters in the double-electrode approach.Due to the impact of the complex geoelectric environment,the electric potential and the electric potential difference are not sensitive enough to respond to minimal leakage.The tiny leaking area cannot be precisely located using the electric potential and electric potential difference.Using the COMSOL Multiphysics software,this study created a standard geoelectric model of the double-electrode method.We calculated a new parameter—the G parameter through secondary electric potential difference—based on the response characteristics of the electric potential and the electric potential difference while the HDPEfilm is leaking.The experiment demonstrates that the G parameter is more sensitive than the electric potential and electric potential difference for detecting the leaking area of HDPE film.The G parameter is more effective at detecting leakage than the electric potential and electric potential difference.The results of this study can be used to locate HDPEfilm leakage areas in a landfill.

Performance analysis of empirical models for predicting rock mass deformation modulus using regression and Bayesian methods

Deformation modulus of rock mass is one of the input parameters to most rock engineering designs and constructions.The field tests for determination of deformation modulus are cumbersome,expensive and time-consuming.This has prompted the development of various regression equations to estimate deformation modulus from results of rock mass classifications,with rock mass rating(RMR)being one of the frequently used classifications.The regression equations are of different types ranging from linear to nonlinear functions like power and exponential.Bayesian method has recently been developed to incorporate regression equations into a Bayesian framework to provide better estimates of geotechnical properties.The question of whether Bayesian method improves the estimation of geotechnical properties in all circumstances remains open.Therefore,a comparative study was conducted to assess the performances of regression and Bayesian methods when they are used to characterize deformation modulus from the same set of RMR data obtained from two project sites.The study also investigated the performance of different types of regression equations in estimation of the deformation modulus.Statistics,probability distributions and prediction indicators were used to assess the performances of regression and Bayesian methods and different types of regression equations.It was found that power and exponential types of regression equations provide a better estimate than linear regression equations.In addition,it was discovered that the ability of the Bayesian method to provide better estimates of deformation modulus than regression method depends on the quality and quantity of input data as well as the type of the regression equation.

Wave-current-induced soil response in marine sediments

This letter presents a new analytical approximation for wave （current）-induced dynamic soil response in marine sediments. In the model, the third-order approximation for wave-current interactions is employed for the flow model, while Biot＇s dynamic poro-elastic model is used to simulate the porous flow in a seabed. The newly analytical solution is validated with the field observations. Based on the solution, effects of currents and wave-nonlinearity on soil response are examined and a parametric study will be carried out to examine the influence of currents on the liquefaction potential. C 2013 The Chinese Society of Theoretical and Applied Mechanics. [doi： 10.1063/2.1301202]

Caving performance through the integration of microseismic activity and numerical modeling at DOZ-PT Freeport,Indonesia

This article describes an undergoing research at PT Freeport, Indonesia, in which the main goal is to use the microseismic information recorded as a result of mining to analyze cave propagation and stress performance on the actual production and fixed infrastructure. At the moment, several numerical experiments have been conducted to correlate the mining activity with the microseismic events using the data collected during year 2005 and 2006. As a result of the preliminary analysis a micro- and a macrocracking envelop were proposed on the basis of computation of stress behavior at the location of the events. Stresses have been computed using standard elastic continuous boundary element models. The correlation between the average source radius and the stress performance has provided a method to propose a macrocracking criterion. Several techniques have been tested to nucleate the microseismic activity around different geological features. This last attempt was aimed to look at potential overstresses induced over the undercut and extraction level drifts. A method was devised to integrate the microseismicity into a 3-dimensional ride distribution model. This model has shown to be very effective to quantify the overstress induced as a result of computing volumetric microseismicity density. The volumetric microseismic model showed to induce overstress up to 10 MPa over a period of two months. The future work will concentrate on the calibration of the integrated model with actual damage observations made at the current mining infrastructure.

Laboratory-scale investigation of the material distribution characteristics of landslide dams in U-shaped valleys

Material distribution characteristics during sliding and depositing is particularly significative to investigate the internal structure and spatial variation of landslide dams,which are fundamentally determining the mechanical and hydraulic behavior and the susceptibility to cause dam failure.However,limited by longevity shortages and special geographic environments,the material distribution characteristics and their formation mechanisms are difficult to observe in the field.Therefore,an experimental apparatus modeling a landslide dam was developed in this paper,designing three sampling methods with two valley states.The internal deposit characteristics,void ratio variation and relative content of the particle size range(PSR)were analyzed,and the mechanics of deposit structure were also delicately ascertained.The results indicate that granular material deposited in valley shows a structure of inverse grain size accumulation in both vertical and horizontal directions,exhibiting spatial variability of particle gradation and void ratio.The characteristic PSR decreases from 22-30 mm in the two-dimensional state to 10-14 mm in the threedimensional state.Vibration excitation and vibration sieve are the intrinsic mechanisms of granular flow segregation,intrinsically inducing the formation of inverse grading deposit structures.Consequently,spatial variability in size is mainly trig gered by segregation,whereas coarse particle content and deposition boundaries merely exacerbate the difference degree.

Numerical Simulation of Vacuum Preloading for Chemically Conditioned Municipal Sludge

Municipal sludge is a sedimentation waste produced during the wastewater process in sewage treatment plants.Among recent studies,pilot and field tests showed that chemical conditioning combined with vacuum preloading can effectively treat municipal sludge.To further understand the drainage and consolidation characteristics of the conditioning sludge during vacuum preloading,a large deformation nonlinear numerical simulation model based on the equal strain condition was developed to simulate and analyze the pilot and field tests,whereas the simulation results were not satisfactory.The results of the numerical analysis of the pilot test showed that the predicted consolidation degree was greater than that measured by the field tests,which is attributed to the relatively low permeability layer formed during the preloading process of the prefabricated vertical drain.To better reflect the consolidation process of the conditioned sludge,a simplified analysis method considering the low permeability layer around the prefabricated vertical drain was proposed.The initial permeability coefficient of the low permeability layer is determined via numerical simulations using finite difference method.The predicted settlement curve was in good agreement with the measured results,which indicated that the numerical simulation based on the equal strain condition considering the relatively low permeability layer can better analyze the consolidation process of ferric chloride-conditioning sludge with vacuum preloading.