Three‐dimensional nonlinear finite element analysis of pile groups in saturated porous media using a new transmitting boundary

The dynamic behaviour of pile groups subjected to an earthquake base shaking is analysed. An analysis is formulated in the time domain and the effects of material nonlinearity of soil, pile–soil–pile kinematic interaction and the superstructure–foundation inertial interaction on seismic response are investigated. Prediction of response of pile group–soil system during a large earthquake requires consideration of various aspects such as the nonlinear and elasto‐plastic behaviour of soil, pore water pressure generation in soil, radiation of energy away from the pile, etc. A fully explicit dynamic finite element scheme is developed for saturated porous media, based on the extension of the original formulation by Biot having solid displacement (u) and relative fluid displacement (w) as primary variables (u–w formulation). All linear relative fluid acceleration terms are included in this formulation. A new three‐dimensional transmitting boundary that was developed in cartesian co‐ordinate system for dynamic response analysis of fluid‐saturated porous media is implemented to avoid wave reflections towards the structure. In contrast to traditional methods, this boundary is able to absorb surface waves as well as body waves. The pile–soil interaction problem is analysed and it is shown that the results from the fully coupled procedure, using the advanced transmitting boundary, compare reasonably well with centrifuge data. Copyright © 2007 John Wiley & Sons, Ltd.     

A New Method of Significance Testing for Correlation-Coefficient Fields and Its Application

Correlation-coefficient fields are widely used in short-term climate prediction research. The most frequently used significance test method for the correlation-coefficient field was proposed by Livezey, in which the number of significantcorrelation lattice(station) points on the correlation coherence map is used as the statistic. However, the method is based on two assumptions:(1) the spatial distribution of the lattice(station) points is uniform;and(2) there is no correlation between the physical quantities in the correlation-coefficient field. However, in reality, the above two assumptions are not valid.Therefore, we designed a more reasonable method for significance testing of the correlation-coefficient field. Specifically, a new statistic, the significant-correlation area, is introduced to eliminate the inhomogeneity of the grid(station)-point distribution, and an empirical Monte Carlo method is employed to eliminate the spatial correlation of the matrix.Subsequently, the new significance test was used for simultaneous correlation-coefficient fields between intensities of the atmospheric activity center in the Northern Hemisphere and temperature/precipitation in China. The results show that the new method is more reasonable than the Livezey method.     

Estimation of Nash's IUH parameters using stochastic search algorithms

In rainfall–runoff studies, it is often necessary to change the duration of a given unit hydrograph. Nash's Instantaneous Unit Hydrograph (IUH) is an ideal method that eliminates the hydrograph duration. This paper presents the results of the application of search algorithms, namely a genetic algorithm and hill climbing, to develop the IUH that minimizes the error between the observed and generated hydrographs. Also the performance of these methods has been compared with that of the classical method used for estimation of IUH, namely the method of moments. The genetic algorithm is a popular search procedure for function optimization that applies the mechanics of natural genetics and natural selection to explore a given search space. Hill climbing is an optimization technique that belongs to the family of local search and algorithms can be used to solve problems that have many solutions, with some solutions better than others. The results obtained from both the genetic algorithm and hill climbing algorithm for estimation of Nash's IUH parameters were compared with the results obtained by the method of moments for storms from two river basins that are located in different climatic regions. It was found that both the genetic algorithm and hill climbing provided improved and consistent results. Copyright © 2008 John Wiley & Sons, Ltd.     

Earthquake magnitude estimation using the τc and Pd method for earthquake early warning systems

Earthquake early warning （EEW） systems are one of the most effective ways to reduce earthquake disaster. Earthquake magnitude estimation is one of the most important and also the most difficult parts of the entire EEW system. In this paper, based on 142 earthquake events and 253 seismic records that were recorded by the KiK-net in Japan, and aftershocks of the large Wenchuan earthquake in Sichuan, we obtained earthquake magnitude estimation relationships using the τe and Pa methods. The standard variances of magnitude calculation of these two formulas are ±0.65 and ±0.56, respectively. The Pd value can also be used to estimate the peak ground motion of velocity, then warning information can be released to the public rapidly, according to the estimation results. In order to insure the stability and reliability of magnitude estimation results, we propose a compatibility test according to the natures of these two parameters. The reliability of the early warning information is significantly improved though this test.     

Numerical simulation of migration–accumulation of oil resources

Numerical simulation of oil migration and accumulation is to describe the history of oil migration and accumulation in basin evolution. It is of great value to the evaluation of oil resources and to the determination of the location and amount of oil deposits. This thesis discusses the characteristics of petroleum geology and permeation fluid mechanics. For the three-dimensional problems of Dongying hollow of Shengli Petroleum Oil Field, it puts forward a new model and a kind of modified method of upwind finite difference fractional steps implicit interactive scheme. For the famous hydraulic experiment of secondary migration–accumulation, the numerical simulation test has been done, and both the computational and experimental results are basically identical. For the actual problem of Dongying hollow, the numerical simulation test and the actual conditions are basically coincident. Thus, the well-known problem has been solved.     

TBM在不良地质洞条件下的施工技术

TBM掘进机以其快速、高效、安全、优质等优点越来越广泛地被应用于隧洞开挖施工中,尤其更适用于深埋超长隧洞,然而在不良地质洞段中TBM掘进缓慢,甚至有卡刀可能,反而不如钻爆法灵活,这就需要根据围岩性状采取特殊技术处理措施,辅以监控量测手段对支护方案进行验证、调整支护措施、修正设计参数等。结合辽宁大伙房输水隧洞工程,总结了在不良地质洞条件下的超前地质预报方法、不良地质段处理措施以及围岩变形监测方法。     

自适应自然单元法研究——自适应细化

对h型自适应自然单元法的自适应细化方案进行了初步研究。在ZZ误差分析的基础上实现了节点的自动加密,使得随着自适应细化的进行,求解误差减小,而且误差分布趋于均匀。在数值分析中,主要有两种误差来源--插值误差和积分误差。随着节点的加密,Delaunay三角形的尺寸随之减小,三角形内的应力场趋于线性分布,那么插值误差和积分误差也都会随之减小。因而,h型自适应分析可以同时减小上述两种误差而达到不断提高求解精度的目的。由于自然单元法求解依赖于求解域内离散节点的Voronoi结构,建议的细化方案中新节点的引入只需局部调整Delaunay结构,算法的实现极为容易,程序实现简单、高效。研究表明,建议的自适应方案是可行的,自然单元法特别适合进行h型自适应分析。     

凸包插值算法的改进与实现

研究了凸包插值算法在n维欧拉空间构造Delaunay三角网的应用,提出了其三角剖分不惟一的重大缺陷,采用了限制区域生长法解决了这一缺陷。     

2008年初地面观测数据资料审核的异常分析及处理

结合2008年1~2月地面气象观测A文件数据的审核情况,针对台站气压突跳的反常记录,通过理论分析,判断是因天气现象变化所致,同时对特殊天气现象如固态降水、积雪的观测记录,提出处理方法。     

Nuclear Waste Disposal Simulations: Couplex Test Cases

We give some results obtained for the Couplex test cases proposed by the ANDRA. In this paper our aim is twofold. Firstly, to compute the release of nuclides out of the repository by concentrating on the 3D near field (Couplex 2). The simulation of the transport phenomena takes into account the dissolution of the glass containers and congruent emissions of the radio-nuclides including filiation chains and some simplified chemistry. Secondly, it is to use the near field computations in order to simulate the nuclide migrations in a 2D far field (Couplex 3). Coupling in between the two simulations takes into consideration the periodicity of the disposal modules and the geometry of the repository described in Couplex 1. The mixed finite element and discontinuous Galerkin methods are used to solve the convection–diffusion equations. In order to handle the nonlinear precipitation/dissolution term, we developed a new iterative technique that combines Picard and Newton–Raphson methods.