Fully coupled elastoplastic hydro-mechanical analysis of unsaturated porous media using a meshfree method

This paper presents the first application of an advanced meshfree method, ie, the edge-based smoothed point interpolation method (ESPIM), in simulation of the coupled hydro-mechanical behaviour of unsaturated porous media. In the proposed technique, the problem domain is spatially discretised using a triangular background mesh, and the polynomial point interpolation method combined with a simple node selection scheme is adopted for creating nodal shape functions. Smoothing domains are formed on top of the background mesh, and a constant smoothed strain, created by applying the smoothing operation over the smoothing domains, is assigned to each smoothing domain. The deformation and flow models are developed based on the equilibrium equation of the mixture, and linear momentum and mass balance equations of the fluid phases, respectively. The effective stress approach is followed to account for the coupling between the flow and deformation models. Further coupling among the phases is captured through a hysteretic soil water retention model that evolves with changes in void ratio. An advanced elastoplastic constitutive model within the context of the bounding surface plasticity theory is employed for predicting the nonlinear behaviour of soil skeleton. Time discretisation is performed by adopting a three-point discretisation method with growing time steps to avoid temporal instabilities. A modified Newton-Raphson framework is designed for dealing with nonlinearities of the discretised system of equations. The performance of the numerical model is examined through a number of numerical examples. The state-of-the-art computational scheme developed is useful for simulation of geotechnical engineering problems involving unsaturated soils.     

用双三次样条函数和GPS资料反演现今中国大陆构造形变场

将中国大陆现今构造变动视为一种连续的地壳变形,利用双三次样条函数模拟了近期GPS测定的大陆内部及周边地区412个测站速率,反演大陆地区自洽的构造变动速度场和应变率场.模拟结果显示:印度板块与欧亚板块的碰撞、挤压是构成中国大陆内部岩石层水平形变的主要驱动力.印度板块在东喜马拉雅构造结深深插入青藏高原,造成地壳大规模的缩短和抬升.青藏高原东南部的喜马拉雅带、拉萨和羌塘地块以及青藏高原东南边的川滇地区,内部构造活动强烈,其内部的构造变形包含地壳碎片的冲断、褶皱和侧向逃逸.大陆地壳(或岩石圈)的增厚,尤其是喜马拉雅山脉南北向的快速缩短和青藏高原东西向的缓慢拉张,大约吸收了印欧板块会聚量的85%,西藏中东地区东西向的拉张速率达到了(16±2.0)mm/a,且顺时针方向扭转明显.印度板块相对欧亚板块运动的欧拉极为(29.7°N, 19.3°E, 0.392°/Ma);华南地块相对于欧亚大陆向东(102°±7.4°)南的运动速率是(11±1.54)mm/a,华南块体相对欧亚板块运动的欧拉极为(62.25°N, 126.56°E, 0.141°/Ma);塔里木地块相对较稳定,其西部运动速度高于东部运动速度,作顺时针方向旋转.总体上讲,中国大陆运动方向为北偏东呈辐射状,从西部近南北方向的运动转向东部地区东南方向的运动,绕东喜马拉雅构造结有一顺时针方向的旋转.横穿喜马拉雅构造带及青藏内部的南北向压缩速率为(19±2.0)mm/a,横穿西天山构造带的南北向压缩平均速率为(13±1.5)mm/a,横穿东天山构造带的南北向压缩平均速率为(6.0±1.4)mm/a.阿尔金断裂带的左旋走滑速率为(6±1.2)mm/a.     

Evolutionary Processes in Low-Mass Binary Systems

We consider the evolution of certain low-mass binaries, incorporating models of (a) internal evolution, (b) tidal friction, (c) dynamo activity driven by an elementary α,Ω dynamo, (d) stellar wind driven by the activity, and (e) magnetic braking as a consequence of wind and poloidal dynamo-generated magnetic field. In some circumstances the stellar wind is found to remove mass on a nuclear timescale, as is necessary to explain some observed systems. We can hope that various uncertainties in the model may be clarified by a careful comparison of the models with such observed quantities as rotation periods. These are modified by processes (a), (b) and (e). Assuming that stellar evolution is slow, rotation rate should in some circumstances represent a balance between magnetic braking trying to slow the star down and tidal friction trying to spin it up. Preliminary attempts are promising, but indicate that some fine tuning is necessary. When there is a third body present, in an orbit which is inclined but not necessarily of short period, the eccentricity of a close binary can be strongly modified by ‘Kozai cycles’. We show that this may complicate attempts to account for spin rates of stars in close binaries.     

利用Lageos 1 SLR资料解算1990-2001年的地球定向参数

利用 12年的Lageos 1卫星激光测距资料 (1990 - 2 0 0 1)解算得到了地球定向参数 (EOP) ,将该序列的结果与同期的EOP(IERS)C0 4进行比较 ,其外符精度为 :极移XP— 0 .4 0mas,YP— 0 .4 2mas ,日长变化Dr— 0 .0 35ms。     

SAMPLING REMOTELY SENSED IMAGERY FOR STORAGE,RETRIEVAL, AND RECONSTRUCTION*

A major difficulty in remote sensing is handling the many data from sensors aboard aircraft and satellites. In this paper we identify an optimal procedure for sampling remotely sensed data before their storage or on their retrieval. The procedure depends on spatial correlation in the scene and uses kriging to estimate values that have been lost. An example in which data from an airborne multispectral scanner could be diminished to only about one tenth without serious loss of precision illustrates the method.     

International program of observations of selected minor planets: first results

International programs of observations of selected minor planets have lasted about 50 years and the last one comes to the end in 2000. The main aim of these observational programs consists in obtaining the orientation of the stellar reference frame with respect to the dynamical one using observations of the bright minor planets. The observations are also useful for the orbital improvements of the asteroids themselves. They are available from the author via e-mail at the address .During the above mentioned period more than 23 000 observations of minor planets, referred to different reference star catalogues, have been obtained. The reduction procedure of observations to the PPM star catalogue is described. The orientation parameters are given and discussed.     

重力插值方法研究

为了根据离散观测数据构制连续空间重力变化图像，分析和讨论了３种数值插值方法，计算结果表明多面函数方法插值精度最高。由于逐步回归分析筛选核函数中心点的计算繁琐，文中提出根据分形理论和Ｓｈａｎｎｏｎ取样定理来确定核函数中心点。对滇西试验场进行模拟试算，插值精度可达到４～５（１０－８ｍｓ－２）。     

Simulation of double cold cores of the 35°N section in the Yellow Sea with a wave-tide-circulation coupled model

Based on the MASNUM wave-tide-circulation coupled numerical model, the temperature structure along 35°N in the Yellow Sea was simulated and compared with the observations. One of the notable features of the temperature structure along 35°N section is the double cold cores phenomena during spring and summer. The double cold cores refer to the two cold water centers located near 122°E and 125°E from the depth of 30m to bottom. The formation, maintenance and disappearance of the double cold cores are discussed. At least two reasons make the temperature in the center (near 123°E) of the section higher than that near the west and east shores in winter. One reason is that the water there is deeper than the west and east sides so its heat content is higher. The other is invasion of the warm water brought by the Yellow Sea Warm Current (YSWC) during winter. This temperature pattern of the lower layer (from 30m to bottom) is maintained through spring and summer when the upper layer (0 to 30m) is heated and strong thermocline is formed. Large zonal span of the 35°N section (about 600 km) makes the cold cores have more opportunity to survive. The double cold cores phenomena disappears in early autumn when the west cold core vanishes first with the dropping of the thermocline position. Supported by the National Basic Research Program of China (No. G1999043809) and the National Science Foundation of China (No. 49736190).     

存储式岩心定向取心技术应用研究

采用岩心定向取心技术不仅可以迅速而方便地确定地下矿层和结构面的产状,而且还可以判断矿层的走向,了解地层倾向,预测钻孔弯曲的趋势。介绍了存储式岩心定向取心原理、关键技术问题、器具构成、技术参数以及现场应用情况。     

Possible Density Change of the Origin of North Pacific Intermediate Water Due to Double Diffusion in the Mixed Water Region

In this study we test Talley's hypothesis that Oyashio winter mixed-layer water (26.5–26.6σ _θ) increases its density to produce the North Pacific Intermediate Water (NPIW) salinity minimum (26.7– 26.8σ_θ) in the Mixed Water Region, assuming a combination of cabbeling and double diffusion. The possible density change of Oyashio winter mixed-layer water is discussed using an instantaneous ratio of the change of temperature and salinity along any particular intrusion (R _l ). We estimate the range of R _l ^DD required to convert Oyashio winter mixed-layer water to the NPIW salinity minimum due to double diffusion, and then assume double-diffusive intrusions as this conversion mechanism. A double-diffusive intrusion model is used to estimate R _l ^DD in a situation where salt fingering dominates vertical mixing, as well as to determine whether Oyashio winter mixed-layer water can become the NPIW salinity minimum. Possible density changes are estimated from the model R _l ^DD by assuming the amount of density change due to cabbeling. From these results, we conclude that Oyashio winter mixed-layer water contributes to a freshening of the lighter layer of the NPIW salinity minimum (around 26.70σ_θ) in the MWR.