基于R树的分布式并行空间索引机制研究

为提高分布式并行计算环境下海量空间数据管理与并行化处理的效率,基于并行空间索引机制的研究,设计一种多层并行R树空间索引结构。该索引结构以高效率的并行空间数据划分策略为基础,以经典的并行计算方法论为依据,使其结构设计在保证能够获得较好的负载平衡性能的前提下,更适合于海量空间数据的并行化处理。以空间范围查询并行处理的系统响应时间为性能评估指标,通过实验证明并行空间索引结构具有设计合理、性能高效的特点。     

A 3D discrete FEM iterative algorithm for solving the water pipe cooling problems of massive concrete structures

Water pipe cooling has been widely used for the temperature control and crack prevention of massive concrete structures such as high dams. Because both under‐cooling and over‐cooling may reduce the efficiency of crack prevention, or even lead to great harm to structures, we need an accurate and robust numerical tool for the prediction of cooling effect. Here, a 3D discrete FEM Iterative Algorithm is introduced, which can simulate the concrete temperature gradient near the pipes, as well as the water temperature rising along the pipes. On the basis of the heat balance between water and concrete, the whole temperature field of the problem can be computed exactly within a few iteration steps. Providing the pipe meshing tool for building the FE model, this algorithm can take account of the water pipe distribution, the variation of water flow, water temperature, and other factors, while the traditional equivalent algorithm based on semi‐theoretical solutions can only solve problems with constant water flow and water temperature. The validation and convergence are proved by comparing the simulated results and analytical solutions of two standard second‐stage cooling problems. Then, a practical concrete block with different cooling schemes is analyzed and the influences of cooling factors are investigated. In the end, detailed guidance for pipe system optimization is provided. Copyright © 2015 John Wiley & Sons, Ltd.     

Selection of sand models and identification of parameters using an enhanced genetic algorithm

Numerous constitutive models of granular soils have been developed during the last few decades. As a consequence, how to select an appropriate model with the necessary features based on conventional tests and with an easy way of identifying parameters for geotechnical applications has become a major issue. This paper aims to discuss the selection of sand models and parameters identification by using genetic algorithm. A real‐coded genetic algorithm is enhanced for the optimization with high efficiency. Models with gradually varying features (elastic‐perfectly plastic modelling, nonlinear stress–strain hardening, critical state concept and two‐surface concept) are selected from numerous sand models as examples for optimization. Conventional triaxial tests on Hostun sand are selected as the objectives in the optimization. Four key points are then discussed in turn: (i) which features are necessary to be accounted for in constitutive modelling of sand; (ii) which type of tests (drained and/or undrained) should be selected for an optimal identification of parameters; (iii) what is the minimum number of tests that should be selected for parameter identification; and (iv) what is the suitable and least strain level of objective tests to obtain reliable and reasonable parameters. Finally, a useful guide, based on all comparisons, is provided at the end of the discussion. Copyright © 2015 John Wiley & Sons, Ltd.     

Modelling combined open channel flow by artificial neural networks

Combined open channel flow is encountered in many hydraulic engineering structures and processes, such as irrigation ditches and wastewater treatment facilities. Extensive experimental studies have conducted to investigate combined flow characteristics. Nevertheless, there is no simple relationship that can fully describe the velocity profiles in a turbulent flow. The artificial neural network (ANN) has great computational capability for solving various complex problems, such as function approximation. The main objective of this study is to evaluate the applicability of the ANN for simulating velocity profiles, velocity contours and estimating the discharges accordingly. The velocity profiles measured by an acoustic doppler velocimeter in the open channel of the Chihtan purification plant, Taipei, with different discharges at fixed measuring section and different depths are presented. The total number of data sets is 640 and the data sets are split into two subsets, i.e. training and validation sets. The backpropagation algorithm is used to construct the neural network. The results demonstrate that the velocity profiles can be modelled by the ANN, and the ANN constructed can nicely fit the velocity profiles and can precisely predict the discharges for the conditions investigated. Copyright © 2005 John Wiley & Sons, Ltd.     

位场DFT算法研究

本文推广了经典的抽样定理,并据此导出了函数有限离散傅里叶变换误差方程（简称DFT误差方程,下同）。该方程把有限离散傅里叶变换中固有的离散效应和有限效应表示为确切的数学形式。离散效应被表示为一个含整参变量（参变量取0,1,…,N-1）的复无穷级数;有限效应被表示为一个含整参变量（参变量取0,1,…,N-1）的复无穷级数的DFT。 基于DFT误差方程和位函数特点,作者提出了两种位场数值傅里叶变换新算法--移样法和等效源续尾叠样法。移样法可近百倍地提高位场数值傅里叶反变换的精度,等效源续尾叠样法可数十倍地提高正变换精度。两种算法都不需要增加资料长度和取样密度,因而基本不需要增加计算机时间和内存。文中给出了算例。     

快速序贯算法

本文讨论了序贯算法的基本原理，与ｋａｌｍａｎ滤波的关系，和序贯算法的渐近性态，得到无先验信息的序贯算法不受初值的影响。实例分析表明，此序贯算法可适于方程系数阵的求逆、求解和病态方程的快速解算。     

神经网络方法在地震预报研究中的初步应用

本文把神经网络方法引进地震预报研究当中。使用地震频次，最大震级，平均震级，等价地震次数等多项地震活动性指标作为神经网络的输入，未来时段内的最大地震震级作为其输出，可以对某一固定地区的最大地震震级作出中近期预报。选用的神经网络模型为含两个中间层的前向模型，并采用ＢＰ算法。所得结果表明，用神经网络方法可以在一定精度范围内使震级预报的内检符合率达到１００％，在本文的例子中，外推预报准确率达到６０％以上。     

遗传算法及其在GOMS模型反演中的应用效果分析

几何光学交互遮蔽模型（GOMS）是一种重要的遥感前向模型,它较好解释了“热点”现象,具有较强的前向模拟能力。但由于其固有的非线性性,给反演带来困难。本文尝试采用近年来兴起的并行随机全局寻优算法－遗传算法对GOMS进行反演,并针对传统遗传算法的不足进行了改进。在使用相同先验知识的条件下,将该算法与目前最有效的约束非线性最优化确定性搜索算法－逐步二次规划法对GOMS模型的反演效果进行了比较,结果表明,逐步二次规划法搜索效率较高,但结果受初值的影响很大,初值选择不当,易收敛于局部最优解,而遗传算法具有全局最优的收敛效果,但局部搜索效率较差。在某些对精度要求不高,而对搜索效率要求较高的场合,可以采用遗传算法与确定性搜索算法相结合的混合遗传算法,以提高算法的搜索效率,获得较为满意的效果。     

多维地震动各方向相关性对拱坝动应力的影响

本文将随机振动的虚拟激励法与拱坝-地基动力相互作用FE-BE-IBE时域模型结合，发展了一个可以考虑多维随机地震动作用下的拱坝动力响应计算模型，并用Monte Garlo方法对模型进行了验证，计算结果表明，地震动分量的相关性对结构的动力响应存在一定影响，合理考虑地震动各方向分量的相关性可以更好地计算实际地震作用下的拱坝动力响应。     

A Parametric Approach for Dealing with Compositional Rounded Zeros

In this work, a parametric approach for replacing data below the detection limit, also known as rounded zeros, in compositional data sets is proposed. Compositional rounded zeros correspond to small proportions of some whole that cannot be reliably detected by the analytical instruments under given operating conditions. This kind of zeros appear frequently in the data collection process in geosciences. They must be treated in an adequate way before some multivariate analysis can be applied. Our procedure results from a modification of the Expectation-Maximization (EM) algorithm and is based on the additive log-ratio transformation. Its coherence with the nature of compositional data and with basic operations in the simplex sample space is checked. Using real data sets, we find that this approach improves other parametric and non-parametric techniques for compositional rounded zeros.