Laplace-transform analytic element solution of transient flow in porous media

A Laplace-transform analytic element method (LT-AEM) is described for the solution of transient flow problems in porous media. Following Laplace transformation of the original flow problem, the analytic element method (AEM) is used to solve the resultant time-independent modified Helmholtz equation, and the solution is inverted numerically back into the time domain. The solution is entirely general, retaining the mathematical elegance and computational efficiency of the AEM while being amenable to parallel computation. It is especially well suited for problems in which a solution is required at a limited number of points in space–time, and for problems involving materials with sharply contrasting hydraulic properties. We illustrate the LT-AEM on transient flow through a uniform confined aquifer with a circular inclusion of contrasting hydraulic conductivity and specific storage. Our results compare well with published analytical solutions in the special case of radial flow.     

Estimation of fundamental periods of shear‐wall dominant building structures

Shear‐wall dominant multistorey reinforced concrete structures, constructed by using a special tunnel form technique are commonly built in countries facing a substantial seismic risk, such as Chile, Japan, Italy and Turkey. In spite of their high resistance to earthquake excitations, current seismic code provisions including the Uniform Building Code (International Conference of Building Officials, Whittier, CA, 1997) and the Turkish Seismic Code (Specification for Structures to be Built in Disaster Areas, Ankara, Turkey, 1998) present limited information for their design criteria. In this study, consistency of equations in those seismic codes related to their dynamic properties are investigated and it is observed that the given empirical equations for prediction of fundamental periods of this specific type of structures yield inaccurate results. For that reason, a total of 80 different building configurations were analysed by using three‐dimensional finite‐element modelling and a set of new empirical equations was proposed. The results of the analyses demonstrate that given formulas including new parameters provide accurate predictions for the broad range of different architectural configurations, roof heights and shear‐wall distributions, and may be used as an efficient tool for the implicit design of these structures. Copyright © 2003 John Wiley & Sons, Ltd.     

Application of extended distinct element method with lattice model to collapse analysis of RC bridges

This paper proposes a simple lattice model for collapse analysis of RC bridges subjected to earthquakes by using the extended distinct element method (EDEM). In the model, a concrete element consists of lumped masses connected to one another by springs, and a reinforcement bar is represented by a discrete model or an integrated model. The proposed lattice model is simple but its parameters are reasonably defined. It has fewer element nodes and connecting springs, which will be of benefit by shortening the CPU time. The processes to determine the initial stiffness of concrete and steel springs, the parameters of the constitutive model and the fracture criteria for springs are described. A re‐contact spring model is also proposed to simulate the re‐contact of the concrete after fracture of springs; and a general grid searching method is used to decrease the CPU time for judging re‐contact after fracture. The lattice model is assessed by numerical simulations and experiments. As an application, a damaged single‐column pier subjected to the Kobe Earthquake in 1995 is analysed by EDEM with the proposed model. The simulation results indicate that the proposed model predicts well qualitatively the collapse process of RC bridges. Copyright © 2003 John Wiley & Sons, Ltd.     

Evaluating foundation mass,damping and stiffness by the least‐squares method

This paper discusses how to use the three‐dimensional (3D) time‐domain finite‐element method incorporating the least‐squares method to calculate the equivalent foundation mass, damping and stiffness matrices. Numerical simulations indicate that the accuracy of these equivalent matrices is acceptable when the applied harmonic force of 1+sine is used. Moreover, the accuracy of the least‐squares method using the 1+sine force is not sensitive to the first time step for inclusion of data. Since the finite‐element method can model problems flexibly, the equivalent mass, damping and stiffness matrices of very complicated soil profiles and foundations can be established without difficulty using this least‐squares method. Copyright © 2003 John Wiley & Sons, Ltd.     

A hybrid indirect boundary element—discrete wave number method applied to simulate the seismic response of stratified alluvial valleys

A hybrid indirect boundary element – discrete wavenumber method is presented and applied to model the ground motion on stratified alluvial valleys under incident plane SH waves from an elastic half-space. The method is based on the single-layer integral representation for diffracted waves. Refracted waves in the horizontally stratified region can be expressed as a linear superposition of solutions for a set of discrete wavenumbers. These solutions are obtained in terms of the Thomson–Haskell propagators formalism. Boundary conditions of continuity of displacements and tractions along the common boundary between the half-space and the stratified region lead to a system of equations for the sources strengths and the coefficients of the plane wave expansion. Although the regions share the boundary, the discretization schemes are different for both sides: for the exterior region, it is based on the numerical and analytical integration of exact Green's functions for displacements and tractions whereas for the layered part, a collocation approach is used. In order to validate this approach results are compared for well-known cases studied in the literature. A homogeneous trapezoidal valley and a parabolic stratified valley were studied and excellent agreement with previous computations was found. An example is given for a stratified inclusion model of an alluvial deposit with an irregular interface with the half-space. Results are displayed in both frequency and time domains. These results show the significant influence of lateral heterogeneity and the emergence of locally generated surface waves in the seismic response of alluvial valleys.     

高填方路堤互锚式薄壁挡土墙土压力试验研究

在山区修建高等级公路工程中，挡土墙结构广泛应用于高填方路堤工程。在提出一种新型互锚式薄壁挡土墙结构的基础上，通过埋设土压力盒的方法对该种挡土墙的结构型式进行了力学模型试验，并进行了模型土压力强度特性分析。试验及分析表明，该挡土墙结构成本低廉，施工方便，结构合理，可以应用于高填方路堤实际工程。     

济阳坳陷稀土元素特征及其在物源对比中的应用

物源是控制沉积物中REE组成最主要的因素，济阳坳陷区岩石稀土元素特征表明：①济阳坳陷古生界与华北地台其他地区同时代地层岩石的REE分布特征具有极大的相似性，这体现了晚古生代济阳坳陷区所处的整个华北地台区为一稳定的克拉通沉积盆地，地层横向分布稳定，具有一致的物源和构造背景；②济阳坳陷古生界为中生界的物源，反映了济阳坳陷区由古生代稳定地台型沉积到中生代山间盆地沉积的转变，中生代洼陷区的沉积主要来自对附近凸起区古生代地层的剥蚀；③新生界样品与中生界样品的REE分布模式具有很大的相似性，从一个侧面反映了济阳坳陷中生代与新生代的构造格局的转变，中生代接受沉积的部分洼陷区至古近纪成为供给物源的凸起区。     

Anomalous Gadolinium,Cerium, and Yttrium Contents in the Adige and Isarco River Waters and in the Water of Their Tributaries (Provinces Trento and Bolzano/Bozen,NE Italy)

Rare earth patterns of surface and groundwaters near big cities often show anthropogenic Gd (Gd_ant) anomalies in addition to geogenic Ce and Y anomalies. The Gd_ant anomaly is caused by very stable organic complexes, one of which is gadopentetic acid, Gd‐DTPA. Derivatives of this and similar compounds are used as contrast agents in magnetic resonance imaging (MRI) of the human blood system. The organic Gd complexes are stable enough to pass nearly unaffected through sewage treatment plants and are, thereafter, discharged into surface water systems. Water of the rapidly flowing Isarco/Eisack and Adige/Etsch rivers (Provinces of Trento and Bolzano/Bozen, NE Italy) and their tributaries show remarkable variations in anthropogenic Gd contents (Gd_ant). Low Gd_ant values are found on Monday and Tuesday, whereas high values are observed during the remaining weekdays. Reliable Gd_ant balances are calculated for the river system at the confluence of the Adige and its tributaries. At two places local decrease of Gd_ant indicates exfiltration of groundwater. It is demonstrated that Gd_ant can be used as a reliably conservative tracer to study the water budget in rapidly flowing alpine river systems. The studied different river waters show considerable negative Ce and positive Y anomalies. Negative Ce anomalies are caused by scavenging of Ce(III) by FeO(OH) precipitates and subsequent oxidation to CeO₂. Y anomalies are attributed to less sorption of Y than REE onto particulate matter. Thus, Y moves faster than REE. Both, Ce and Y anomalies, are of geogenic origin.     

利用卫星测高空间技术建立数字化海洋

论述了卫星测高空间技术的主要特征、科学意义及其在军事和经济建设领域的重要价值，阐明了国内外测高卫星系统的发展现状及利用卫星测高空间技术建立数字化海洋所取得的成果，分析卫星测高目前的研究方向，提出了发展我国卫星测高空间技术以及利用该技术建立数字化海洋的一些建议。     

Reliability based risk index for the design of reinforced earth structures

The design methods currently used for earth reinforcement are mostly based on deterministic properties of both the soil and the construction materials used. Nowadays, however, the general trend is designing at a specific degree of reliability. This is even more true where the raw data such as soil properties exhibit significant variation. Deterministic solutions, in this case, may not suffice. Therefore, this paper will attempt to use probabilistic formulations thereby modifying the existing design procedure of reinforced earth retaining walls to account for uncertainties and variabilities. Through a first order Taylor's series expansion about the mean, the mean and variance of the strip reinforcing components, namely width and length, are derived in terms of the variations in the soil properties. Design charts that enable estimation of both mean and variance are developed to avoid extensive partial differentiation involved in the computations. Using appropriate probability distributions along with the mean and variance, the final design outputs are determined for a selected failure probability by introducing what is refered to as 'risk index'. The results indicate that the risk index increases with an increase in the coefficient of variations and a decrease in failure probability. Furthermore, it is shown that in some cases, depending on the variabilities of the soil properties, the classical design technique produced a relatively high failure probability. This revised version was published online in July 2006 with corrections to the Cover Date.