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.     

试论新世纪地理信息科学的发展

本文重点分析了当前国内外地理信息科学的形成与发展趋势和地理信息科学不断发展的历程 ,提出近期发展的前沿 ,以及当前地理信息产业化面临的一些重要问题。最后提出对地理信息科学发展的意见。     

Canonical theory of perturbations using eccentric anomaly as independent variable

It is shown in this paper how to build a canonical transformation of variables, so that the eccentric anomaly becomes the new independent variable. In the case of eccentric elliptical orbits it changes the equations of motion so, that they can be integrated analytically to any order of approximation comparatively easy.     

On the expansions of intermediate orbit for general planetary theory

Intermediate orbit for general planetary theory is constructed in the form of multivariate Fourier series with numerical coefficients. The structure and efficiency of the derived series are illustrated by giving various statistical properties of the coefficients.The ability of the recently proposed elliptic function approach to compress the Fourier series representing the intermediate orbit is investigated. Our results confirm that when mutual perturbations of a pair of planets are considered the elliptic function approach is quite efficient and allows one to compress the series substantially. However, when perturbations of three or more planets are under study the elliptic function approach does not give any advantages.