Computational homogenization of uncoupled consolidation in micro‐heterogeneous porous media     

Fredrik Larsson  Kenneth Runesson  Fang Su 《国际地质力学数值与分析法杂志》2010,34(14):1431-1458

Variationally consistent homogenization is exploited for the analysis of transient uncoupled consolidation in micro‐heterogeneous porous solids, whereby the classical approach of first‐order homogenization for stationary problems is extended to transient problems. Homogenization is then carried out in the spatial domain on representative volume elements (RVE), which are introduced in quadrature points in standard fashion. Along with the classical averages, a higher‐order conservation quantity is obtained. An iterative FE²‐algorithm is devised for the case of nonlinear permeability and storage coefficients, and it is applied to pore pressure changes in asphalt‐concrete (particle composite). Various parametric studies are carried out, in particular, with respect to the influence of the ‘substructure length scale’ that is represented by the size of the RVE's. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

Analytical layer‐element method for non‐axisymmetric consolidation of multilayered soils     

Zhi Yong Ai  Wen Ze Zeng 《国际地质力学数值与分析法杂志》2012,36(5):533-545

A numerically efficient and stable method is developed to analyze Biot's consolidation of multilayered soils subjected to non‐axisymmetric loading in arbitrary depth. By the application of a Laplace–Hankel transform and a Fourier expansion, the governing equations are solved analytically. Then, the analytical layer‐element (i.e. a symmetric stiffness matrix) describing the relationship between generalized displacements and stresses of a layer is exactly derived in the transformed domain. Considering the continuity conditions between adjacent layers, the global stiffness matrix of multilayered soils is obtained by assembling the inter‐related layer‐elements. Once the solution in the Laplace–Hankel transformed domain that satisfies the boundary conditions has been obtained, the actual solution can be derived by the inversion of the Laplace–Hankel transform. Finally, numerical examples are presented to verify the theory and to study the influence of the layered soil properties and time history on the consolidation behavior. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

Statistical estimation of mineral age by K-Ar method     

A. B. Vistelius  E. R. Drubetzkoy  A. V. Faas 《Mathematical Geology》1989,21(8):905-919

Statistical estimation of age of⁴⁰Ar/⁴⁰K ratios may be considered a result of convolution of uniform and normal distributions with different weights for different minerals. Data from Gul'shad Massif (Nearbalkhash, Kazakhstan, USSR) indicate that⁴⁰Ar/⁴⁰K ratios reflecting the intensity of geochemical processes can be resolved using convolutions. Loss of⁴⁰Ar in bioties is shown whereas hornblende retained the original content of⁴⁰Ar throughout the geological history of the massif. Results demonstrate that different estimation methods must be used for different minerals and different rocks when radiometric ages are employed for dating.  相似文献   

Mountain building and earth rotation     

L. L. A. Vermeersen  R. Sabadini  G. Spada  N. J. Vlaar 《Geophysical Journal International》1994,117(3):610-624

  相似文献   

盆地模拟技术中的化学动力学生烃模型与应用     

赵永军  李汉林  查明 《地质论评》1994,40(7):84-89

油气地表化探数据的失真(特别是数据中的高值失真点)直接影响到化探异常圈定。本文根据油气地表化探数据的分布特点，利用概率统计原理，提出一种适合油气化探资料前期处理的失真数据的概率判定方法和处理手段，在对实际资料的应用中，取得了良好的效果。  相似文献   

Some comments on the 3-D impulse response of a Maxwell viscoelastic earth     

J. X. Mitrovica  J. L. Davis 《Geophysical Journal International》1995,120(1):227-234

  相似文献   

IMPLICIT NONLINEAR NORMAL MODE INITIALIZATION:A MULTIGRID APPROACH          下载免费PDF全文

D.R.C.Nair  B.Chakravarty  P.Niyogi 《Acta Meteorologica Sinica》1993,(1)

The implicit nonlinear normal mode initialization (INMI) is applied to a tropical limited area shallow water modelin spherical coordinates.The boundary condition for the INMI scheme is based on the boundary formulation of themodel.The INMI scheme is found to be very efficient in suppressing spurious gravity wave oscillation and providing awell balanced initial data set for the model.The INMI scheme involves solving a number of elliptic type equations withvarying complexity.and hence an efficient numerical technique is required for solving such equations.In order to makeINMI computationally more attractive,we are employing the multigrid method for solving all the elliptic type equationsin the INMI scheme.The numerical procedures for the development of such multigrid solvers are briefly described.  相似文献   

孔隙水压力圆锥静力触探及其发展现状     

张喜发 《吉林大学学报(地球科学版)》1990,(4)

介绍孔隙水压力圆锥静力触探这门新技术,着重评述近年来在探头优化设计方面的进展,以及在划分地层和鉴别土类、估计土的超固结比和固结系数方面应用的最新研究成果。  相似文献   

An examination of the use of the logratio transformation for the testing of endmember hypotheses     

Ross M. Renner 《Mathematical Geology》1991,23(4):549-563

Given a compositional dataset in the absence of any prior information on any mixing process which may have formed it, a complete analysis of mixtures determines three distinct types of estimates in order. These are: (i) the estimate of the number of endmembers or fixed source compositions, of which all the sample compositions of the dataset must be approximate mixtures; (ii) the estimated compositions for each of these chosen number of endmembers; and (iii) the estimated contributions of each of these endmember estimates to each sample. Traditionally, the estimate for the number of endmembers has been assessed either by mapping or by inspection of the coefficients of determination between the observed and estimated variables. Mapping entails the plotting on a map of the region from which the samples were taken, either the contours of the contributions of each endmember to each sample, or some other portrayal of the distribution of endmember abundances. Because it requires the complete analysis, assessment by this method is too elaborate except for final confirmation and display. Alternatively, choosing a number of endmembers, which result in suitability high coefficients of determination for all or most variables, may account for elements which are not part of the conjectured mixing process or, worse, may result in the identification of endmembers which may never in fact have existed. Such an error is similar to overspecifying a multiple regression model. So, the obvious starting point from which to assess the validity, or otherwise choice of endmember numbers, is to examine the matrix of residuals. The differences between the logratio-transformed observed and estimated data form an array of residual logratios. A linear combination of these may be formed for each sample, which, under a random perturbation assumption, should follow a univariate normal distribution. Whether or not this scalar is normal can be readily tested. It can also be examined graphically for such desirable qualities as symmetry when the test for normality may be too severe. This procedure is employed to assess the decompositions of the U.S.G.S. Mid-Pacific data and the Nazca Plate Surface sediments.This paper was presented at the 18th Geochautauqua, Newark, Delaware, 13–14 October 1989.  相似文献   

水库区应力－孔压耦联分析     

梁青槐  高士钧  曾心传 《华南地震》1994,(3)

基于Ｂｉｏｔ耦联固结理论［１］，推导出二维饱水弹性岩石介质的应力－孔压耦联微分方程，并用有限元法和有限差分法分别对其进行了空间和时间的离散化。  相似文献