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.     

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.     

重力插值方法研究

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

海上不同高度风速换算关系的研究

本文用海上石油平台上观测的风梯度资料,计算了在不同风况下海面摩擦速度和粗糙高度。结果表明,海面摩擦速度u_*与海上10m高度处风速u_(10)的关系为u_*=0.055u_(10)-0.058。在4—5级风时,海面平均摩擦速度为50cm/s,粗糙度为0.022m,它相当于陆上的平均粗糙高度。在风速为2—3级时,二者与海上通常采用的数据接近。文中最后给出了上述两种风况下海上不同高度的换算系数。     

JONSWAP谱和N谱的几个对应特征量间的关系

研究了风浪频谱ＪＯＮＳＷＡＰ谱的峰度因子γ、尖度因子Ｐｊ和Ｎ谱尖度因子ＰＮ，参变量ｐ，ｑ间的关系；并在此关系上给出了由这两种谱计算的几个对应参量间的关系，在使用不同的谱形进行特征量的计算时，应注意不同频谱间的成长因子致性，以免引起计算结果的错误。     

南海波高熵和风速熵

根据风速的统计分布，给出了有因次风速熵和无因次风速熵的定义及其计算方法，使用ＧＥＯＳＡＴ高度计１９８６年１１月－１９８９年２月的有效波高和风速的资料，计算，分析了南海海域上的波高熵，风速熵，给出它们的时间变化特征和空间变化特征，并对不同随机量的无因次熵，即随机度进行了比较。     

Triple diagram method for the prediction of wave height and period

Many formulations have been developed so far to predict the wave height and period from fetch length and wind blowing duration for a constant wind speed. This study aimed to predict wave parameters from fetch length and meteorological factors by using triple diagram methodology based on Kriging principles. Proposed model results were compared with Joint North Sea Wave Project (JONSWAP) model which is used so commonly in the ocean and coastal engineering studies. For the implementation of the methodology hourly wave and wind data were obtained from a buoy located in Lake Ontario. Numerical and graphical comparisons demonstrated that the proposed method outperforms the classical formulation.     

Atmospheric forcing of the eastern tropical Pacific: A review

The increase in marine, land surface, atmospheric and satellite data during recent decades has led to an improved understanding of the air–sea interaction processes in the eastern tropical Pacific. This is also thanks to extensive diagnoses from conceptual and coupled ocean–atmosphere numerical models. In this paper, mean fields of atmospheric variables, such as incoming solar radiation, sea level pressure, winds, wind stress curl, precipitation, evaporation, and surface energy fluxes, are derived from global atmospheric data sets in order to examine the dominant features of the low level atmospheric circulations of the region. The seasonal march of the atmospheric circulations is presented to depict the role of radiative forcing on atmospheric perturbations, especially those dominating the atmosphere at low levels.In the tropics, the trade winds constitute an important north–south energy and moisture exchange mechanism (as part of the low level branch of the Hadley circulation), that determines to a large extent the precipitation distribution in the region, i.e., that associated with the Inter-Tropical Convergence Zone (ITCZ). Monsoonal circulations also play an important role in determining the warm season precipitation distribution over the eastern tropical Pacific through a large variety of air–sea–land interaction mechanisms. Westward traveling waves, tropical cyclones, low latitude cold air intrusions, and other synoptic and mesoscale perturbations associated with the ITCZ are also important elements that modulate the annual rainfall cycle. The low-level jets of the Gulf of California, the Intra-Americas Sea (Gulf of Mexico and Caribbean Sea) and Chocó, Colombia are prominent features of the eastern tropical Pacific low-level circulations related to sub-regional and regional scale precipitation patterns. Observations show that the Intra-Americas Low-Level Jet intensity varies with El Niño/Southern Oscillation (ENSO) phases, however its origin and role in the westward propagation and development of disturbances that may hit the eastern tropical Pacific, such as easterly waves and tropical cyclones, are still unclear. Changes in the intensity of the trade winds in the Caribbean Sea and the Gulf of Mexico (associated with eastern tropical Pacific wind jets) exert an important control on precipitation by means of wind–topography interactions. Gaps in the mountains of southern Mexico and Central America allow strong wind jets to pass over the continent imprinting a unique signal in sea surface temperatures and ocean dynamics of the eastern tropical Pacific.The warm pools of the Americas constitute an important source of moisture for the North American Monsoon System. The northeastern tropical Pacific is a region of intense cyclogenetic activity, just west of the coast of Mesoamerica. Over the oceanic regions, large-scale properties of key variables such as precipitation, moisture, surface energy fluxes and wind stress curl are still uncertain, which inhibits a more comprehensive view of the region and stresses the importance of regional field experiments. Progress has been substantial in the understanding of the ocean and atmospheric dynamics of the eastern tropical Pacific, however, recent observational evidence such as that of a shallow meridional circulation cell in that region, in contrast to the classic concept of the Hadley-type deep meridional circulation, suggests that more in situ observations to validate theories are still necessary.This paper is part of a comprehensive review of the oceanography of the eastern tropical Pacific Ocean.