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.     

Long-Term Effect of TMD on Vibration Control of An MDOF Offshore Fixed Platform

A three-dimensional fixed offshore platform in deep water modeled by the finite element method is studied in this paper. Analysis of the dynamic response of the MDOF structure is realized taking the non-linearity of the wave drag force and the wave-structure interaction into account. The structural response statistics, which have Gaussian distributions, are used to evaluate the vibration effect of the structure without TMD and with TMD. And an optimal method to design TMD controlling the first mode of the multi-mode structure is proposed. Moreover, the probabilities of occurrence of sea states at the platform site are considered for prediction of the long-term effect of a TMD. Simulation results demonstrate that the long-term effect of a well-designed TMD is good and the practical use is possible due to the good stability of its optimal parameters under different sea states.     

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.     

二次特征值问题的最佳向后扰动分析

二次特征值问题 (QEP)的主要的求解方法之一是转化为广义特征值问题 (GEP) ,然后用求解广义特征值的方法 (比如 QZ方法 )求解。本文研究由此获得的计算解的范数意义下的最佳向后扰动分析 ,所得结果是 Tisseur最近所得结果的加强。     

受扰线性系统的前馈-反馈最优控制

针对具有外部持续扰动的线性系统，研究前馈-反馈最优控制律的设计问题。给出了最优控制律的存在唯一性条件。并提出了最优控制律的设计算法。利用滤波器解决了前馈控制的物理不可实现问题。仿真结果表明，此算法易于实现，与传统的反馈最优控制相比对抑制外部扰动具有较强的鲁棒性。     

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

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

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

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

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.     

GIS中直线元内插点精度及对误差带的影响

基于误差传播定律，考虑参数r误差影响，推导了线元内插点的精度计算公式，讨论内插点精度对线元误差带的影响，并对影响的结果进行了分析，得到了一些有益的结论。     

地学自由曲面径向基函数网络重建的详细算法

径向基函数(RadialBasisFunction,简称RBF)神经网络是一种理想的地学离散数据网格化工具,能够适应各种不同分布形式和边界条件的数据,收敛速度较快,可以逼近任何复杂曲面。这里详细介绍了RBF神经网络的算法。适当径向基函数的形式和偏差系数是使用RBF神经网络作地学曲面重建的关键。大量的实际数据验证结果表明,当选用Gauss型径向基函数时,一般可获得比较理想的网格化效果,同时具备"曲面平滑"和"拟合度高"的特点。