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.     

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.     

Solution of an unstable axisymmetric Cauchy–Poisson problem of dispersive water waves for a spectrum with compact support

 It has been known that the axisymmetric Cauchy–Poisson problem for dispersive water waves is well posed in the sense of stability. Thereby time evolution solutions of wave propagation depend continuously on initial conditions. However, in this paper, it is demonstrated that the axisymmetric Cauchy–Poisson problem is ill posed in the sense of stability for a certain class of initial conditions, so that the propagating solutions do not depend continuously on the initial conditions. In order to overcome the difficulty of the discontinuity, Landweber–Fridman's regularization, famous and well known in applied mathematics, are introduced and investigated to learn whether it is applicable to the present axisymmetric wave propagation problem. From the numerical experiments, it is shown that stable and accurate solutions are realized by the regularization, so that it can be applicable to the determination of the ill-posed Cauchy–Poisson problem.     

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.     

Optimizing interpolation of shoot density data from a Posidonia oceanica seagrass bed

A case study on the optimization of Posidonia oceanica density interpolation, using a data set from a large meadow at Porto Conte Bay (NW Sardinia, Italy), is presented. Ordinary point kriging, cokriging and a weighted average based on inverse square distance were used to interpolate density data measured in 36 sampling stations. The results obtained from different methods were then compared by means of a leave‐one‐out cross‐validation procedure. The scale at which interpolation was carried out was defined on the basis of the Hausdorff dimension of the variogram. Optimizing spatial scale and data points search strategy allowed obtaining more accurate density estimates independently of the interpolation method.     

重力异常垂向一阶导数的一种简便算法

利用拉格朗日插值导出了重力异常垂向一阶导数的计算公式，给出了地面以及地面之上不同高度的求导系数。该公式可以计算地面上的重力垂向一阶导数，还可以直接计算地面之上任意高度上的垂向一阶导数。鉴于该公式除系数不同外，与上延公式完全相同，因此，程序设计尤为简单。使用本公式对模型数据和实际资料进行了处理，证明了本算法的实用性。     

基于钻孔数据的三维地质建模与可视化研究

三维地质模型能够完整准确地表达复杂地质现象的边界条件及地质体内含的各种地质构造,直观地再现地质单元的空间展布及其相互关系,最大限度地提高地质分析的直观性和准确性,因此三维地质建模已引起地质、采矿、岩土工程等诸多领域的日益重视。笔者在综合考察多种建模方法的基础上,结合工程实践,提出了一种实用的工程地质三维建模实施方法。该方法以钻孔资料和地质剖面数据为三维地质建模的源数据,允许用户手工修正模型剖面地层分界线,从而实现对建模过程的干预和控制,克服了单纯依靠钻孔数据建模结果不精确且难以修正的问题。最后,通过一个研究实例展示了该方法的实际建模效果。     

随机Dirichlet问题解的调和展开

针对重力学随机Dirichlet问题,通过适当地对边界检验函数的分解,并在随机边界样本空间中提取确定性部分的对偶基,本文将随机Dirichlet问题的一般解展开为一随机系数的调和级数形式。     

An integrable case of a rotational motion analogous to that of Lagrange and Poisson for a gyrostat in a Newtonian force field

The scope of the present paper is to provide analytic solutions to the problem of the attitude evolution of a symmetric gyrostat about a fixed point in a central Newtonian force field when the potential function isV ⁽²⁾.We assume that the center of mass and the gyrostatic moment are on the axis of symmetry and that the initial conditions are the following: (t ₀)=₀, (t ₀)=₀, (t ₀)=(t ₀)=₀, ₁(t ₀)=0, ₂(t ₀)=0 and ₃(t ₀)= ₃ ⁰ .The problem is integrated when the third component of the total angular momentum is different from zero (B ₁ 0). There now appear equilibrium solutions that did not exist in the caseB ₁=0, which can be determined in function of the value ofl ₃ ^r (the third component of the gyrostatic momentum).The possible types of solutions (elliptic, trigonometric, stationary) depend upon the nature of the roots of the functiong(u). The solutions for Euler angles are given in terms of functions of the timet. If we cancel the third component of the gyrostatic momentum (l ₃ ^r =0), the obtained solutions are valid for rigid bodies.     

Interpolation by regularized spline with tension: I. Theory and implementation

Bivariate and trivariate functions for interpolation from scattered data are derived. They are constructed by explicit minimization of a general smoothness functional, and they include a tension parameter that controls the character of the interpolation function (e.g., for bivariate case the surface can be tuned from a membrane to a thin steel plate), Tension can be applied also in a chosen direction, for modeling of phenomena with a simple type of anisotropy. The functions have regular derivatives of all orders everywhere. This makes them suitable for analysis of surface geometry and for direct application in models where derivatives are necessary. For processing of large datasets (thousands of data points), which are now common in geosciences, a segmentation algorithm with a flexible size of overlapping neighborhood is presented. Simple examples demonstrating flexibility and accuracy of the functions are presented.On leave from the Department of Physical Geography and Cartography, Comenius University, Mlynská dolina, Bratislava, Czechoslovakia.On leave from the Institute of Physics, Dúbravská cesta 9, Bratislava, Czechoslovakia.