Fast iterative solution of large undrained soil‐structure interaction problems

In view of rapid developments in iterative solvers, it is timely to re‐examine the merits of using mixed formulation for incompressible problems. This paper presents extensive numerical studies to compare the accuracy of undrained solutions resulting from the standard displacement formulation with a penalty term and the two‐field mixed formulation. The standard displacement and two‐field mixed formulations are solved using both direct and iterative approaches to assess if it is cost‐effective to achieve more accurate solutions. Numerical studies of a simple footing problem show that the mixed formulation is able to solve the incompressible problem ‘exactly’, does not create pressure and stress instabilities, and obviate the need for an ad hoc penalty number. In addition, for large‐scale problems where it is not possible to perform direct solutions entirely within available random access memory, it turns out that the larger system of equations from mixed formulation also can be solved much more efficiently than the smaller system of equations arising from standard formulation by using the symmetric quasi‐minimal residual (SQMR) method with the generalized Jacobi (GJ) preconditioner. Iterative solution by SQMR with GJ preconditioning also is more elegant, faster, and more accurate than the popular Uzawa method. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical study of enhanced energy dissipation near a seamount

Using a two-dimensional primitive equation model, we examine nonlinear responses of a semidiurnal tidal flow impinging on a seamount with a background Garrett-Munk-like (GM-like) internal wavefield. It is found that horizontally elongated pancake-like structures of high vertical wavenumber near-inertial current shear are created both in the near-field (the region over the slope of the seamount) and far-field (the region over the flat bottom of the ocean). An important distinction is that the high vertical wavenumber near-inertial current shear is amplified only at mid-latitudes in the far-field (owing to a parametric subharmonic instability (PSI)), whereas it is amplified both at mid-and high-latitudes (above the latitude where PSI can occur) in the near-field. In order to clarify the generating mechanism for the strong shear in the near-field, additional numerical experiments are carried out with the GM-like background internal waves removed. The experiments show that the strong shear is also created, indicating that it is not caused by the interaction between the background GM-like internal waves and the semidiurnal internal tides. One possible explanation is proposed for the amplification of high vertical wavenumber near-inertial current shear in the near-field where tide residual flow resulting from tide-topography interaction plays an important role in transferring energy from high-mode internal tides to near-inertial internal waves.     

Influence of Cyclic Wetting Drying on Collapse Behaviour of Compacted Residual Soil

The climatic zones where residual soils occur are often characterized by alternate wet and dry seasons. Laboratory studies of earlier workers have established that the alternate wetting and drying process affects the swell-shrink potentials, water content, void ratio and particle cementation of expansive soils. The influence of cyclic wetting and drying on the collapse behaviour of residual soils has not been examined. This paper examines the influence of alternate wetting and drying on the collapse behaviour of compacted residual soil specimens from Bangalore District. Results of such a study are useful in anticipating changes in collapse behaviour of compacted residual soil fills. Experimental results indicated that the cyclic wetting and drying process increased the degree of expansiveness of the residual soils and reduced their collapse tendency. Changes in the swell/collapse behaviour of compacted residual soil specimens from wetting drying effects are attributed to reduction in water content, void ratio and possible growth of cementation bonds.     

Compressibility and shear strength of a residual soil

A detailed study was performed to investigate the compressibility and the shear strength of a residual soil in Omdurman, Sudan. The soil resulted from weathering of sandstone and mudstone. The soil is classified as sandy clayey silt (MH-CH) and is generally partially saturated. Consolidation tests were performed to study the compressibility of the soil and UU, CIU and insitu tests (CPT and SPT) were performed to study the shear strength characteristics of the soil. The characteristics of this residual soil can generally be predicted using correlations proposed in the literature for transported soils. The soil is of low to medium compressibility and exhibits slight apparent overconsolidation behavior. The soil dilates during shear and shows no long term effective cohesion. Problems associated with open cuts in this soil are attributed to loss of strength upon saturation.     

冲绳海槽玄武岩中中酸性残余熔体研究及其岩石学意义

与洋壳有关的酸性岩由于对了解幔源岩浆的演化以及判别古老蛇绿岩套及其构造位置的重要意义而倍受岩石学家和构造地质学家的关注。本对发现于冲绳海槽玄武岩基质中的中酸性到酸性残余熔体进行了详细研究，它们提供了幔源玄武岩浆结晶分异形成酸性岩浆的直接证据。在细小的基质矿物间分布有一种玻璃质的残余熔体，其成分随距冷凝边距离(L)的增加而越来越酸性。在SiO2对Na2O K2O图解上，残余熔体的投影点从玄武岩到英安岩均有分布，反映了一个连续的演化系列。在AFM图解上，残余熔体表现出与Thingmuli火山岩系列类似的拉斑玄武岩系列的演化趋势。我们的研究表明：残余熔体的演化受结晶分异作用控制。在早期结晶阶段，辉石的结晶起主导作用，结果造成残余熔体中SiO2、Al2O3，Na2O的迅速增加，FeO、MgO、CaO迅速降低。在晚期结晶阶段，斜长石成为主导结晶相，导致残余熔体中Al2O3，Na2O的迅速消减。Al2O3、Na2O从增加到降低的转变出现在SiO2=62％左右。在L=27．5mm处，85～90％的基质岩浆已发生了结晶作用，导致残余熔体中SiO2含量达到69～70％，而且此处还新出现了一种富FeTi的氧化物。该玄武岩中残余熔体和基质矿物的成分及演化特征分别与Thingmuli火山岩系列中酸性端元的组成相似，在也佐证了Thingmuli火山岩系列是幔源岩浆结晶分异的产物.     

京珠高速公路东西湖段软基处理数值分析

采用等效地基固结系数法,对京珠高速公路东西湖段软土地基,进行了比奥固结平面有限元分析。计算值与实测值比较接近。计算结果为选择地基处理措施和确定道面施工时间提供了依据,保证了竣工后该路段的工后沉降满足设计要求。     

利用相关分析划分重力异常

本文从理论上对相关分析法区分重力异常的可行性进行了论述，并从一元线性回归方程着手，对大地水准面上的重力值与相应的地形高程作了相关分析，通过分析，认为在区域性小比例主均衡补偿基本完善的条件下，重力值与高程有极好的线性相关关系；当计算窗口较大时，地下的局部不均匀体，不足以破坏重力值与高程的相头发一，所以，在实际工作中，利用相关分析处理重力资料的叠加场，可以收到较好的效果，且方法简单易行。     

模拟退火剩余静校正的一步算法

当地震资料的剩余静校正量大、信噪比低时，利用传统的线性反演方法不能有效地拾取静校正量。本文介绍一种非线性反演方法（一步模拟退火法）来估算炮点和检波点的剩余静校正量。通过对合成地震记录的试算表明，该算法可获得近似于真实的剩余静校正量。