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

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

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

求解以水位为变量的连续方程,并根据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.     

Paléoenvironnements et caractérisation des roches mères pétrolières des séries pré-salifères du bassin intérieur du Gabon

The Interior Basin of Gabon, created during the break-up between South America and Africa, displays thick Neoproterozoic to Aptian p.p. fluvio-lacustrine deposits overlain by Aptian to Albian marine facies. Rock–Eval analyses from outcrop and drillhole samples show high content in organic matter (up to 25%) related to types I and II. These intervals are encountered within Permian, Neocomian–Barremian as well as Aptian siliciclastic succession. They constitute fairly good to excellent potential petroleum source rocks, which are most probably at the origin of oil indices recognized both in drillholes and in surface.     

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

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

饱水黄土隧道变形规律研究

新松树湾隧道位于饱和的砂粘土中,土体的稳定性极差,施工难度较大。为此,在隧道施工期间,通过土工试验、隧道收敛和围岩内部位移的现场监测、隧道变形的三维有限元仿真计算等手段,对该隧道的变形规律进行了系统研究,获得了大量试验数据。研究成果对隧道的施工组织和支护参数的确定起到了重要作用,为此类隧道的施工积累了经验。