Upper bound limit analysis of uplift failure in pressurized sealed rock tunnels

The rock around tunnels used for gas storage is subject to high pressures, reaching 30 MPa in the case of compressed air energy storage. Uplift failure of the overlaying rock mass up to the surface represents the main hazard scenario in such cases. The present paper investigates this problem by using the upper bound theorem of limit analysis assuming a continuum rock mass model obeying the Mohr Coulomb failure criterion with tension cut‐off. Tools of the calculus of variations are used to assess the geometry of the failure surface. The effects of geometrical and geotechnical parameters on uplift pressure are analyzed systematically. Charts are then provided, which enable a quick estimation of the upper bound of uplift pressure across a wide range of geotechnical and geometrical conditions.     

采用误差迭代补偿的高精度相位解缠方法

最小二乘相位解缠算法是一种全局算法,快速稳定简单易实现,但该算法未绕过相位不连续区易造成误差的全局扩散,导致结果存在较大误差.针对最小二乘相位解缠算法的缺陷,提出一种高精度相位解缠方法.该方法在基于快速傅里叶变换的最小二乘解缠算法基础上,采用误差迭代补偿技术,补偿解缠误差.仿真和实测相位数据的实验结果表明,该方法可有效提高解缠精度,且在仿真实验的定量分析中得出,在强相位噪声条件下,该方法的解缠精度比最小二乘算法提高一个数量级.     

地铁隧道三维激光扫描数据配准方法

针对传统的迭代最近点算法(ICP)用于多测站点云数据配准时计算效率低的问题,该文提出了一种基于特征点的ICP改进算法,该方法利用相邻两测站数据进行配准的实现。首先采用体素化格网方法对两点云数据集进行精简处理,并计算精简处理后每一点的法向量;然后利用kd-tree最近邻查询搜索特征点之间的对应关系;并通过估计出的最优变换矩阵更新至全局变换,以提高配准精度。实验结果表明,改进的ICP算法在地铁隧道点云数据配准中的效率高于其他的配准方法,为隧道变形监测工作的进行提供保证。     

利用区域地基GPS的电离层层析成像研究

电离层是近地空间的重要组成部分,如何对电离层的异常扰动进行合理监测与预报一直是空间物理领域的研究课题。本文将计算机层析成像技术引入到电离层扰动监测中,利用大量的区域地基GPS观测数据,借助代数层析迭代算法反演得到三维电离层电子密度;并将层析结果与国际电离层参考模型IRI2007进行对比分析,结果表明:地基GPS层析所得的电离层电子密度与IRI2007基本一致,但层析结果精度略高于IRI2007模型。     

重磁异常界面反演中的计算方法

迭代法和子空间法为重力异常界面反演最优化求解计算的两种方法。传统的最优化计算立意虽然简单明确，但迭代计算过程中必须的简化使得求解不稳定，加进一些限制条件之后，使计算量明显增大，以致使求解变得困难。通过引入子空间计算在矩阵求解中减少了工作量，使得加约束的最优化计算变得容易实现。通过理论分析，表明利用子空间方法的解的计算结果准确且效率较高。     

Combining Sample and Census Data in Small Area Estimates: Iterative Proportional Fitting with Standard Software*

The combination of detailed sample data with less detailed but fully enumerated marginal subtotals is the focus of a wide range of research. In this article we advocate careful modeling of sample data, followed by Iterative Proportional Fitting (IPF). The modeling aims to estimate accurately the interaction or odds ratios of complex tables, which is information not contained in the marginal subtotals. IPF ensures consistency with the subtotals. We advance this work in three practical ways. First, we show that detailed small‐area estimates of both counts and proportional distributions usually gain accuracy by combining data for larger areas containing the small areas, and we illustrate the multilevel framework to achieve these estimates. Second, we find that a general classification or socioeconomic typology of the small areas is even more associated with the within‐area interactions than is membership of the larger area. Third, we show how the Statistical Package for the Social Sciences (SPSS) can be used for IPF in any number of dimensions and with any structure of constraining marginal subtotals. Throughout, we use an example taken from the 1991 U.K. Census. These data allow us to evaluate various methods combining 100 percent tabulations and the Samples of Anonymised Records.     

神农架松柏镇水文地球化学研究

近年来，神农架松柏镇居民的结石病时有发生．并有增长的势头。为了解地质环境尤其是饮用水对病患的影响，笔者现场对患病情况以及各种水源(包括河水、泉水、自来水)的来源、化学组成和放射性进行了调查，运用水质常规方法和原子吸收光谱分别测定了离子和微量元素。同时，与地质环境类似的木鱼镇进行了对比调查。研究结果显示，两地地表水、地下水的各种指标均未超出“生活饮用水水源水质卫生标准”。(CJ3020-93)，且松柏镇水体中的放射性(γ射线)低于木鱼镇，而居民的患病率却远高于木鱼镇，可见松柏镇的水环境与居民的结石病无直接关系。     

Non‐iterative time integration scheme for non‐linear dynamic FEM analysis

This paper proposes a non‐iterative time integration (NITI) scheme for non‐linear dynamic FEM analysis. The NITI scheme is constructed by combining explicit and implicit schemes, taking advantage of their merits, and enables stable computation without an iteration process for convergence even when used for non‐linear dynamic problems. Formulation of the NITI scheme is presented and its stability is studied. Although the NITI scheme is not unconditionally stable when applied to non‐linear problems, it is stable in most cases unless stiffness hardening occurs or the problem has a large velocity‐dependent term. The NITI scheme is applied to dynamic analysis of the non‐linear soil–structure system and computation results are compared with those by the central difference method (CDM). Comparison shows that the stability of the NITI scheme is superior to that of the CDM. Accuracy of the NITI scheme is verified because its results are identical with those by the CDM in which the time step is set as 1/10 of that for the NITI scheme. The application of the NITI scheme to the mesh‐partitioned FEM is also proposed. It is applied to dynamic analysis of the linear soil–structure system. It yields the same results as a conventional single‐domain FEM analysis using the Newmark β method. This result verifies the usability of mesh‐partitioned FEM analysis using the NITI scheme. Copyright © 2003 John Wiley& Sons, Ltd.     

Asteroid Mean Elements: Higher Order and Iterative Theories

Mean orbital elements are obtained from osculating ones by removing the short periodic perturbations. Large catalogues of asteroid mean elements need to be computed, as a first step in the computation of proper elements, used to study asteroid families. The algorithms for this purpose available so far are only accurate to first order in the masses of the perturbing planet; the mean elements have satisfactory accuracy for most of the asteroid belt, but degraded accuracy in the neighbourhoods of the main mean motion resonances, especially the 2:1. We investigate a number of algorithms capable of improving this approximation; they belong to the two classes of Breiter-type methods and iterative methods. The former are obtained by applying some higher order numerical integration scheme, such as Runge–Kutta, to the differential equation whose solution is a transformation removing the fast angular variables from the equations; they can be used to compute a full second order theory, however, only if the full second order determining function is explicitly computed, and this is computationally too cumbersome for a complicated problem such as the N-body. The latter are fixed point iterative schemes, with the first order theory as an iteration step, used to compute the inverse map from mean to osculating elements; formally the method is first order, but because they implement a fixed frequency perturbation theory, they are more accurate than conventional single iteration methods; a similar method is already in use in our computation of proper from mean elements. Many of these methods are tested on a sample of asteroid orbits taken from the Themis family, up to the edge of the 2:1 resonance, and the dispersion of the values of the computed mean semimajor axis over 100 000 years is used as quality control. The results of these tests indicate that the iterative methods are superior, in this specific application, to the Breiter methods, in accuracy and reliability. This is understood as the result of the cancellations occurring between second order perturbation terms: the incomplete second order theory, resulting from the use of a Breiter method with the first order determining function only, can be less accurate than complete, fixed frequency theories of the first order. We have therefore computed new catalogues of asteroid mean and proper elements, incorporating an iterative algorithm in both steps (osculating to mean and mean to proper elements). This new data set, significantly more reliable even in the previously degraded regions of Themis and Cybele, is in the public domain. This revised version was published online in July 2006 with corrections to the Cover Date.     

On the colocation iterative solution model and algorithm for gravity anomaly

ＯｎｔｈｅｃｏｌｏｃａｔｉｏｎｉｔｅｒａｔｉｖｅｓｏｌｕｔｉｏｎｍｏｄｅｌａｎｄａｌｇｏｒｉｔｈｍｆｏｒｇｒａｖｉｔｙａｎｏｍａｌｙＹＵＡＮ－ＸＩＹＡＮＧ（杨元喜）ａｎｄＣＨＡＮＧ－ＪＩＡＮＬＩＵ（刘长建）ＺｈｅｎｇｚｈｏｕＩｎｓｔｉｔｕｔｅｏｆＳ...