Three-dimensional finite element simulation of fracture propagation in rock specimens with pre-existing fissure(s) under compression and their strength analysis

A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials.     

Effects of the geometry of two‐dimensional fractures on their hydraulic aperture and on the validity of the local cubic law

Flow through rough fractures is investigated numerically in order to assess the validity of the local cubic law for different fracture geometries. Two‐dimensional channels with sinusoidal walls having different geometrical properties defined by the aperture, the amplitude, and the wavelength of the walls' corrugations, the corrugations asymmetry, and the phase shift between the two walls are considered to represent different fracture geometries. First, it is analytically shown that the hydraulic aperture clearly deviates from the mean aperture when the walls' roughness, the phase shift, and/or the asymmetry between the fracture walls are relatively high. The continuity and the Navier–Stokes equations are then solved by means of the finite element method and the numerical solutions compared to the theoretical predictions of the local cubic law. Reynolds numbers ranging from 0.066 to 66.66 are investigated so as to focus more particularly on the effect of flow inertial effects on the validity of the local cubic law. For low Reynolds number, typically less than 15, the local cubic law properly describes the fracture flow, especially when the fracture walls have small corrugation amplitudes. For Reynolds numbers higher than 15, the local cubic law is valid under the conditions that the fracture presents a low aspect ratio, small corrugation amplitudes, and a moderate phase lag between its walls.     

Solution for a plane strain rough‐walled hydraulic fracture driven by turbulent fluid through impermeable rock

The impact of turbulent flow on plane strain fluid‐driven crack propagation is an important but still poorly understood consideration in hydraulic fracture modeling. The changes that hydraulic fracturing has experienced over the past decade, especially in the area of fracturing fluids, have played a major role in the transition of the typical fluid regime from laminar to turbulent flow. Motivated by the increasing preponderance of high‐rate, water‐driven hydraulic fractures with high Reynolds number, we present a semianalytical solution for the propagation of a plane strain hydraulic fracture driven by a turbulent fluid in an impermeable formation. The formulation uses a power law relationship between the Darcy‐Weisbach friction factor and the scale of the fracture roughness, where one specific manifestation of this generalized friction factor is the classical Gauckler‐Manning‐Strickler approximation for turbulent flow in a rough‐walled channel. Conservation of mass, elasticity, and crack propagation are also solved simultaneously. We obtain a semianalytical solution using an orthogonal polynomial series. An approximate closed‐form solution is enabled by a choice of orthogonal polynomials embedding the near‐tip asymptotic behavior and thus giving very rapid convergence; a precise solution is obtained with 2 terms of the series. By comparison with numerical simulations, we show that the transition region between the laminar and turbulent regimes can be relatively small so that full solutions can often be well approximated by either a fully laminar or fully turbulent solution.     

Micromechanical approach to effective viscoelastic properties of micro‐fractured geomaterials

The aim of this paper is to formulate a micromechanics‐based approach to non‐aging viscoelastic behavior of materials with randomly distributed micro‐fractures. Unlike cracks, fractures are discontinuities that are able to transfer stresses and can therefore be regarded from a mechanical viewpoint as interfaces endowed with a specific behavior under normal and shear loading. Making use of the elastic‐viscoelastic correspondence principle together with a Mori‐Tanka homogenization scheme, the effective viscoelastic behavior is assessed from properties of the material constituents and damage parameters related to density and size of fractures. It is notably shown that the homogenized behavior thus formulated can be described in most cases by means of a generalized Maxwell rheological model. For practical implementation in structural analyses, an approximate model for the isotropic homogenized fractured medium is formulated within the class of Burger models. Although the approximation is basically developed for short‐term and long‐term behaviors, numerical applications indicate that the approximate Burger model accurately reproduce the homogenized viscoelastic behavior also in the transient conditions.     

农村信用社产权制度创新研究

现行农村信用社存在的出资人所有权虚置、收益权削弱、法人财产权丧失等产权制度缺陷 ,严重阻碍着它的进一步发展。形成这些缺陷的主要原因在于农村信用社产权改革的方向不明和对产权范畴的片面理解 ,因此 ,农村信用社产权制度创新应在清晰界定产权的基础上、向股份制、股份合作制、合作制等多种产权形式的方向发展。     

A STUDY ON LAND SYSTEM INNOVATIONS IN THE CONSTRUCTION OF SMALL TOWNS IN CHINA

The economic development in China and the rising of the living standard need to speed up the urbanisa-tion.The development of small towns is an important way to Chinese urbanization.Land use plays a very important role in the development of small towns.However there are many problems in the development of small towns,esp.in land use.The paper first discusses the land problems in the development of small towns.Such as much cultivated land lies idle,under-utilization and waste of land,increasing illegal use of land,unstable contractural relationship for land use.The relationship between the development of small towns and land use is also studied.Then the guidelines for the land system innovations of small towns are put forward.Namely the sustainability of social and economic development,the bal-ance between land reservation and land utilization,the provision of service to village ,agriculture and farmer,the manage-ment of land resource and land assets,the parsimonious andlegal use of land.The basic framework of land system innova-tions of small towns is put forward finally.It include the land replacement policy for small towns and the permanent tenan-cy of farmland,the overall plan for land uses and other plans in harmony,the establishment of a flexible system of land supplies,using land with payments ,the transfer of agricultural land and refining the land law related to the construction of small towns.     

FEASIBLE STUDY ON THE INTEGRATION SYSTEM FOR THE SPACE MONITORING OF MAJOR EARTHQUAKES AND VOLCANOES IN TERRESTRIAL LAND

With the rapid development of space technology, earth observation technology and sky observatory technolo-gy, they have played a more and more important part in monitoring and predicting of earthquakes and volcanoes in the terres-trial land. In recent years, the related agencies have done the experiments and researches on monitoring and predicting ofearthquakes and volcanoes in the forewarning period by means of many approaches, such as satellite thermal infrared re-mote sensing (TIRS), Global Positioning System (GPS), differential interferometric synthesis aperture radar (D-INSAR),astronomical time-latitude residual anomaly, and Geographic Information Systems (GIS), etc. A quite large number of re-search foundation has been built in the fundamental theories and application methods. The experiments and researcheshave shown that these technology is efficient methods for high frequency crust movement. If the existed separate scientificforces and results are possibly assembled together to form a more complete integration monitoring system with the combina-tion of space, sky observation, ground, deep geology and macro anomaly, it will come into a new stage of monitoring andpredicting of earthquakes and volcanic eruptions.     

GDJ人卫目视观测数据采集和处理系统的研制

本文介绍了ＧＤＪ跟踪打印经纬仪在人造卫星的观测中如何应用计算机对多圈卫星的目视观测资料进行实时采集、处理等工作过程的原理和方法。     

流星暴的航天定义

定义ZHR^*以每小时穿过1000平方公里天顶截面的，能在铝质表面打出直径不小于1厘米的弹坑的流星体数，作为流星群对航天安全威胁强度的定量标志，给出了从ZHR到ZHR^*的换算公式．用ZHR^*取代ZHR来评价20世纪90年代以来的一些强流星雨，发现1998年的贾可比尼(天龙座)流星雨对航天安全来说，是比公认的1999、2001、2002年3次狮子座流星暴强数倍的更强流星暴．     

中国1∶100万景观生态制图设计

中国1∶100万景观生态图系借鉴国内外传统景观制图方法和制图规范,在遥感、地理信息系统先进技术的支持下,设计研制其制图方法、制图内容、样图和技术流程,同时提出初步的景观分类系统,为今后编制中国1∶100万景观生态图奠定基础。