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.     

迁客离忧楚地颜--略说贬谪文学与荆湘地域之关系及其特点

古代荆湘之地域的特殊性使其成为朝廷贬官吏的首选之地和贬官的中转之地。在此 ,这些被贬官首创了大量关于贬谪的文学作品。其作品物色鲜明 ,楚地风情突出 ,具有近乎凝固化了的悲怨传统 ,早期逐臣、作品及谪地成为凝定的贬谪代码和符号 ,生成并积淀了深厚的人文精神。     

Analytical solutions of one‐dimensional infiltration before and after ponding

This paper is a continuation of previous research, which obtained a convenient solution for arbitrary surface fluxes before ponding. By means of Fourier Transformation this has been extended to derive analytical solutions of a linearized Richards' equation for arbitrary input fluxes after surface saturation. Copyright © 2003 John Wiley & Sons, Ltd.     

关于t-模的生成

讨论了概率度量空间中ｔ－模生成的概率意义。对构造ｔ－模提供了具体方法，并以实例介绍了这种构造方法。     

Sea‐surface temperature and salinity changes in the northwest Pacific since the Last Glacial Maximum

The oxygen isotope records of both benthic and planktonic Foraminifera in five piston cores, collected from the region between the Oyashio and Kuroshio Currents near Japan, clearly show the marked latitudinal shifts of these two currents during the past 25 kyr. Under the present hydrographic condition, a clear relationship between the sea‐surface temperature (SST) and oxygen isotope differences from benthic to planktonic Foraminifera is observed in this region. Using this relationship, we find decreased SSTs of 12–13°C (maximum 15°C) in the southernmost core site at the Last Glacial Maximum (LGM), indicating the Oyashio Current shifted southward. The SSTs at the southern two core sites abruptly increased more than 10°C at 10–11 ka, suggesting the Kuroshio Current shifted northward over these sites at 10–11 ka. In contrast, the northern two core sites have remained under the influence of the cold Oyashio Current for the past 25 kyr. With the reasonable estimate of bottom‐water temperature decrease of 2.5°C at the LGM, the SSTs estimated by this new method give exactly the same SST values calculated from Mg/Ca ratio of planktonic Foraminifera, allowing palaeosea‐surface salinities to be reconstructed. The result suggests that the ice volume effect was 1.0 ± 0.1‰ at the LGM. Copyright © 2004 John Wiley & Sons, Ltd.     

山东省产业重心转移与可持续发展研究

 本文利用经济空间结构的重心研究方法,分析了山东省近八年来产业重心的转移,提出了山东省产业可持续发展战略     

Motive force of Qinghai-Tibet Plateau moving to east

A calculation formula on spherical pattern of Qinghai-Tibet plateau moving model is established. Tibet massif moves norward by east in speed of 28 mm/a, Ganshu-Qinghai massif moves to northeast in speed of 15 mm/a, Qomolangma Feng moves northward by a few east in speed of 35 ～42 mm/a. The low latitude perimeter is longer than the nigh latitude perimeter. When the Tibet massif moves northward, its latitude perimeter must be contracted and the Tibet massif must move eastward by Coriolis. Coriolis force is inertial in earth rotation. It makes the fall body turning to east and the rising block turning westward. In the Northern Hemisphere, it makes the northward body turning to east and the southward block turning to west.This is the reason why the tectonic zones of western Pacific are different from those of eastern Pacific.     

Seismic rays in an ellipsoidal earth model: theoretical analysis and estimation of deviations due to the ellipticity

Summary. A first-order form of the Euler's equations for rays in an ellipsoidal model of the Earth is obtained. The conditions affecting the velocity law for a monotonic increase, with respect to the arc length, in the angular distance to the epicentre, and in the angle of incidence, are the same in the ellipsoidal and spherical models. It is therefore possible to trace rays and to compute travel times directly in an ellipsoidal earth as in the spherical model. Thus comparison with the rays of the same coordinates in a spherical earth provides an estimate of the various deviations of these rays due to the Earth's flattening, and the corresponding travel-time differences, for mantle P -waves and for shallow earthquakes. All these deviations are functions both of the latitude and of the epicentral distance. The difference in the distance to the Earth's centre at points with the same geocentric latitude on rays in the ellipsoidal and in the spherical model may reach several kilometres. Directly related to the deformation of the isovelocity surfaces, this difference is the only cause of significant perturbation in travel times. Other differences, such as that corresponding to the ray torsion, are of the first order in ellipticity, and may exceed 1 km. They induce only small differences in travel time, less than 0.01s. Thus, we show that the ellipticity correction obtained by Jeffreys (1935) and Bullen (1937) by a perturbational method can be recovered by a direct evaluation of the travel times in an ellipsoidal model of the Earth. Moreover, as stated by Dziewonski & Gilbert (1976), we verify the non-dependence of this correction on the choice of the velocity law.     

REVIEWS OF GEOGRAPHIC SOFTWARE

Reviews of geographic software in this article: DEMO-GRAPHICS: WORLD POPULATIONS AND PROJECTIONS. ESP GAUSS. CEMODEL S. Damus LIMDEP. William H. Greene MICROSTAT 4.1 OTIS PCIPS. (Personal Computer Image Processing System) . H.J. Meyers and R. Bernstein. REGRESSION ANALYSIS OF TIME SERIES (RATS) SPSS/PC+ URBAN DATA MANAGEMENT SOFTWARE (UDMS)     

The importance of oxygen isotope provenance in relation to solute content of bulk meltwaters at Imersuaq Glacier,West Greenland

Stable oxygen isotope analysis and measurement of several dissolved cations and anions of bulk meltwater samples have provided information about the hydrochemical environment of the glacial hydrological system at Imersuaq Glacier, an outlet tongue from the Greenland ice‐sheet, West Greenland. The samples were collected at frequent intervals during the period 20–28 July 2000 in a small (<20 L s^?1) englacial meltwater outlet at the glacier margin. The results document the following findings: (i) a marked diurnal variation of δ¹⁸O is related to the composition of oxygen isotope provenances, mainly near‐marginal local superimposed ice and basal up‐sheared ice further up‐glacier; (ii) a relationship is seen between all base cations (Na⁺, K⁺, Ca²⁺, Mg²⁺), SO₄^2? and δ¹⁸O, indicating that solute acquisition is provided by solid–solution contact with the up‐sheared ice—as the relationship with Cl^? is weak the influence of seasalt‐derived solutes is small in the area; (iii) when the melt rate is high, two diurnal maxima of δ¹⁸O values and solute concentrations are measured, and it is suggested that a snow meltwater component is responsible for the second maximum of δ¹⁸O—a short residence time leads to a delayed decrease in ion concentrations. Copyright © 2003 John Wiley & Sons, Ltd.