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.     

重力插值方法研究

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

现代黄河三角洲土地利用变化分析

人类对海岸带的开发利用影响着海岸的变迁。土地利用和土地覆被变化可以说是自然与人文过程交叉最为密切的问题，是人地关系的具体体现。通过对黄河三角洲土地利用变化的研究，可了解人为因素对三角洲海岸带的开发利用状况及其对海岸侵蚀的影响。根据1987年和2001年两期卫星遥感图像和2001年野外调查数据，运用地理信息系统方法及数理统计方法，分析了现代黄河三角洲地区1987年至2001年土地利用／覆被的数量变化和空间变化特征；揭示了该地区各类土地利用数量变化的幅度、速度和土地利用空间变化的主要类型；阐明了该区各类土地利用变化的区域特点。该地区土地利用变化速度比整个环渤海地区速度快，草地、滩地和不明地类这3种类型大量减少，其他类型增多。从变化率来看，养殖区、库塘、油田用地增长最快。库塘、海滩养殖区等类型的增加，滨海草地的减少，使海岸侵蚀加剧。     

北太平洋亚热带环流区表层水中226Ra和228Ra的分布与变化

对1999年9～10月采自北太平洋亚热带环流区的19份表层海水样品的Ra同位素分析表明。研究海域表层水中的^226Ra、^228Ra放射性比度分别介于0．67～0．92、0．08～0．30Bq／m^3之间，平均值分别为0．74、0．11Bq/m^3．^226Ra/^228Ra)A．R．活度比的变化范围为0．11～0．44，平均值为0．19．上述数值明显低于近岸海域水体的相应值，表现为典型的开阔大洋水的特征．从空间分布的特征看，研究海域Ra同位素含量与^226Ra/^228Ra)A.R.值均呈均匀分布态势．将本研究结果与历史数据进行对比后发现，本研究获得的^226Ra、^228Ra放射性比度比20世纪60～80年代得到的数据来得低，可能与水体层化作用加强导致的Ra补充量的减少以及生物生产力升高导致的Ra迁出量的增加有关．北太平洋亚热带环流区表层水中Ra同位素的时间变化与文献报道的该海域叶绿素a、硅酸盐、磷酸盐含量与初级生产力的历史变化趋势相吻合．     

遥感与GIS支持下的海洋渔业空间分布研究——以东海为例

基于海洋渔业多年来生产实践所表明的海洋渔场与海洋水文要素密切相关性,将研究区的遥感融合信息与生产数据进行GIS空间配准,同时建立了诸多要素数据的空间聚类模式,利用该模式提取了水文信息和中心渔场信息相关联的空间分布规律的隐伏信息,由此阐述了实现海洋渔业现代化,应用空间遥感融合信息和GIS支持下指导海洋渔业生产的渔情预报有着重要意义.     

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.     

Separation of Waves by Interpolation Method

1 .ＩｎｔｒｏｄｕｃｔｉｏｎＴｏｍｅａｓｕｒｅｔｈｅｒｅｆｌｅｃｔｉｏｎｏｆｉｎｃｉｄｅｎｔｗａｖｅｓｐｒｏｄｕｃｅｄｂｙａｐｈｙｓｉｃａｌｍｏｄｅｌｉｎａｗａｖｅｆｌｕｍｅｉｓａｃｏｍｍｏｎｐｒｏｂｌｅｍｗｉｔｈｐｈｙｓｉｃａｌｔｅｓｔｓｏｆｗａｖｅａｃｔｉｏｎｏｎｃｏａｓｔａｌｓｔｒｕｃｔｕｒｅｓ .Ｗａｖｅｓｇｅｎｅｒａｔｅｄｂｙｔｈｅｗａｖｅ ｍａｋｅｒｐｒｏｐａｇａｔｅｆｏｒｗａｒｄｉｎｔｈｅｗａｖｅｆｌｕｍｅａｎｄａｒｅｒｅｆｌｅｃｔｅｄｂｙｔｈｅｐｈｙｓｉｃａｌｍｏｄｅｌ,ａｎｄｔｈｅ…