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.     

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

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

Preliminary Study on the Translunar Halo Orbits of the Real Earth–Moon System

The computation of translunar Halo orbits of the real Earth–Moon system (REMS) has been an open problem for a long time, but now, it is possible to compute Halo orbits of the REMS in a systematic way. In this paper, we describe the method used for the numerical computation of Halo orbits for a time span longer than 41 years. Halo orbits of the REMS are computed from quasi-periodic Halo orbits of the quasi-bicircular problem (QBCP). The QBCP is a model for the dynamics of a spacecraft in the Earth–Moon–Sun system. It is a Hamiltonian system with three degrees of freedom and depending periodically on time. In this model, Earth, Moon and Sun are moving in a self-consistent motion close to bicircular. The computed Halo orbits of the REMS are compared with the family of Halo orbits of the QBCP. The results show that the QBCP is a good model to understand the main features of the Halo family of the REMS.     

长江三峡工程地壳形变监测网络

长江三峡工程地壳形变监测网络，采用了当今高精度GPS、INSAR空间大地测量技术，并与精密水准测量，精密重力测量，精密激光测距和峒体连续形变监测等技术相结合，构成一个空间上点、线、面结合，时间上长、中、短兼顾的高精度，高时空分辨率的地壳形变监测网络。该网络既可获取三峡库区特别是库首区区域形变场和区域应变场的动态变化，监测库区主要断层活动，为水库诱发地震预测及研究服务，又可用于气象、滑坡地质灾害监测等，该监测网络于1997年底开始建设，2001年6月建成，到目前为止，已获得大量宝贵的观测资料，必将产生显著的经济效益和社会效益。     

关于t-模的生成

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

银河系运动学参数的确定

利用 2 0 2个太阳附近疏散星团的视向速度和自行观测资料 ,对太阳的运动和银河系的运动学参数进行了研究。其中 ,距离在 0 .5kpc到 2kpc之间的 12 8个疏散星团对平均太阳运动分量的解算结果是 (u0 ,v0 ,w0 ) =(- 13.8± 1.4 ,- 5 .0± 1.6 ,- 11.6± 2 .9)km/s ;Oort常数和银河系径向运动参数的解算结果分别为 (A ,B) =(16 .9± 1.1,- 11.6± 2 .6 )km·s- 1·kpc- 1及 (C ,D) =(2 .5± 1.1,- 2 .1± 0 .9)km·s- 1·kpc- 1。     

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

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

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.     

Effects of Different Diets on the Dietary Attractability and Selectivity of Chinese Shrimp, Fenneropenaeus chinensis

Attractabilities of different diets and dietary selectivity of Chinese shrimp, Fenneropenaeus chinensis were studied through behavior observation and feeding experiment, respectively. The five diets used in the experiment are: Fish Flesh (FF), Shrimp Flesh (SF), Clam Foot (CF), Polychaete Worm (PW), and Formulated Diet (FD). No significant differences of attractability exist between any two different diets when every two natural diets or all five diets are provided simultaneously. On the other hand, significant differences of attractability exist between FD and every single natural diet when they are provided simultaneously. Results of behavioral observation indicate that natural diets are more attractive than FD. In feeding experiment, Chinese shrimp has distinct selectivity on different diets. It positively selects CF and PW, negatively selects FF and SF, and excludes FD absolutely. The results of the present studies indicate that the dietary selectivity of shrimp was based not only on the attractabilities of the diets, but also on the responses such as growth and food conversion.     

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.