森林沼泽区矿产资源地球化学勘查

通过景观调查和中大比例尺化探方法试验,对我国东北森林沼泽区景观地球化学特征、元素存在形式、迁移富集机制、影响因素、采样介质、采样粒度、采样方法等理论和方法技术问题进行研究,提出该类地区元素表生分散富集的某些地球化学规律和中大比例尺化探工作方法.     

成都市浮尘物理与化学特征的初步研究

通过对成都市一环、二环路交叉处和纵贯市区人民路浮尘进行粒度、X射线衍射物相、元素组成等分析,研究了成都市浮尘的成分、元素含量、组合特征及其在空间的分布特征,对降低浮尘污染提出了初步建议.     

考虑几何非线性的强夯流固动力耦合分析

针对强夯法加固地基机理的流固动力耦合分析问题,本文考虑土体的几何非线性,建立土体非线性动力平衡方程和整体流固动力耦合方程,给出土体非线性动力平衡方程的迭代计算格式和流固动力耦合方程的解耦计算方法,并对一具体算例进行了数值分析,得到了地基位移、应力等在强夯作用下的变化规律和在空间上的分布特征,为精确模拟分析强夯加固机理提供了有效途径。     

浅水方程数值计算方法的研究

求解以水位为变量的连续方程,并根据Navier-Stokes方程压力修正算法的基本思想,建立了浅水方程的水位修正算法,放宽了对离散时间步长的限制.通过对离散方程系数矩阵的重新构造,建立了高分辨率有限元格式,该格式既具有较高的离散精度又避免了数值解的伪振荡.对动量方程的阻力项做负坡线性化处理,提高了露滩计算的稳定性.数值模拟结果与解析解吻合良好,表明所建立的数值计算方法是正确的和可靠的.     

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.     

岩溶地区桩基的有限元分析

岩溶地区长桩存在着遇洞率的问题。为了研究岩溶地区桩的受力特点,采用有限元法模拟岩溶地区桩基的工作方式,分析了岩层厚度对桩的影响以及桩的受力状态,并研究了去桩的可能性。计算结果表明:岩层厚度对桩的影响较大,岩层的位置对桩的应力不产生影响,仅对桩的位移产生影响。