△T球体异常的三元解释法

在△Ｔ异常解释课题中，一个新的角度参量─—视磁化特征角的拟定，最大限度地简化了球体磁场表达式。从而才有可能通过特殊解法，给出了球体磁场的极值与拐点公式。作者据此建立了四种不受正常场制约、仅需异常的局部持征曲线即可全面解释三度异常的解释方法。屏障分解、叠加判别、三元迭代的设计，为解释三度体叠加异常给出了简单而易行的途径。     

Fast iterative solution of large undrained soil‐structure interaction problems

In view of rapid developments in iterative solvers, it is timely to re‐examine the merits of using mixed formulation for incompressible problems. This paper presents extensive numerical studies to compare the accuracy of undrained solutions resulting from the standard displacement formulation with a penalty term and the two‐field mixed formulation. The standard displacement and two‐field mixed formulations are solved using both direct and iterative approaches to assess if it is cost‐effective to achieve more accurate solutions. Numerical studies of a simple footing problem show that the mixed formulation is able to solve the incompressible problem ‘exactly’, does not create pressure and stress instabilities, and obviate the need for an ad hoc penalty number. In addition, for large‐scale problems where it is not possible to perform direct solutions entirely within available random access memory, it turns out that the larger system of equations from mixed formulation also can be solved much more efficiently than the smaller system of equations arising from standard formulation by using the symmetric quasi‐minimal residual (SQMR) method with the generalized Jacobi (GJ) preconditioner. Iterative solution by SQMR with GJ preconditioning also is more elegant, faster, and more accurate than the popular Uzawa method. Copyright © 2003 John Wiley & Sons, Ltd.     

郑州国际会展中心展厅钢屋盖的动力特性分析

郑州国际会展中心展览厅钢屋盖采用了空间斜拉式索拱结构体系，为探讨该体系的动力特性，采用ANSYS有限元程序，分别建立钢屋盖与下部混凝土结构共同作用、钢屋盖单独作用以及无斜拉系统的钢屋盖单独作用时的3个空间计算模型，采用子空间迭代法计算其自振周期和振型，并加以比较。从而得出钢屋盖的动力特性，为进一步研究此类结构的受力特性奠定了基础。     

求解GM(1,1)模型新总体最小二乘算法

在GM(1,1)模型中系数矩阵和观测向量都是由原始序列组成的。系数矩阵中同样是有误差的,与观测向量中的误差一样,亦来源于原始序列,即它们误差同源。不同位置的相同元素应该有相同的改正数,采用传统总体最小二乘求解则不能达到此目的。针对这一缺陷,推导了一种新的总体最小二乘算法;并且通过算例验证了新方法的可行性和有效性。     

Helmert方差分量估计严密公式与简化公式等价性的证明

在分析了Helmert方差分量估计严密公式与简化公式迭代收敛结果一些统计性质的基础上,导出了关于收敛结果共同的求解方程,从理论上证明了两种结果的等价性,为使用简化公式提供了依据。     

由计算机产生的等距离区域分界线

本文介绍一种最近由作者在美国研制的,由计算机产生等距离区域分界线的程序及其数据结构,同时给出初步的试验结果。     

空间直角坐标至两类常用坐标的快速变换

由空间直角坐标至站心极坐标的间接变换方法存在着使用复杂和不够直接快速的问题,该文提出了空间直角坐标至大地坐标快速变换的二维迭代方法;通过分析空间直角坐标与站心极坐标的几何关系,提出了站心极坐标直接变换的几何解析法。计算结果表明,当迭代停止条件采用空间点位误差限值10-5 m时,大地坐标变换的二维迭代法一般2次收敛;站心极坐标几何解析法与常规间接变换结果完全相同,验证了文中所提方法的正确性和快速性。     

整体最小二乘的迭代解法

在引入整体最小二乘平差准则的基础上,推导了整体最小二乘的迭代解法;同时,引入多元函数隐函数求导的方法以确定未知参数对观测数据的线性信息,解决了整体最小二乘下的精度评定问题。给出了运用新的解法在拟合函数确定以及坐标转换参数确定等方面的应用实例,验证了新算法的可行性。     

迭代法建立的低频模型在储层反演中的应用

储层预测的精度直接制约着油气开发的经济性与有效性,而低频信息的补偿是改善储层反演效果的有利途径之一.常规的测井曲线内插低通滤波建模方法难以精细表征复杂地质背景下的低频模型.为降低研究区特殊地质体对储层反演带来的影响,采用迭代法建立低频模型,建模中分步考虑压实作用、特殊岩性等因素影响,并通过反演效果不断迭代更新低频模型,最终准确建立反映不同岩性(包含火山岩)的反演低频模型.研究表明:该方法能够保证储层反演结果的可靠性,为储层精细描述提供了资料;同时基于反演结果,对东海西湖凹陷W气田主力层(P1)有利储层进行重新刻画,优化了开发井水平段的位置,该井实施结果证实钻遇了优质储层.用迭代法建立的低频模型为复杂地质背景地区储层预测提供了一种更有效的建模方法.     

Accelerated initial stiffness schemes for elastoplasticity

Iterative methods for the solution of non‐linear finite element equations are generally based on variants of the Newton–Raphson method. When they are stable, full Newton–Raphson schemes usually converge rapidly but may be expensive for some types of problems (for example, when the tangent stiffness matrix is unsymmetric). Initial stiffness schemes, on the other hand, are extremely robust but may require large numbers of iterations for cases where the plastic zone is extensive. In most geomechanics applications it is generally preferable to use a tangent stiffness scheme, but there are situations in which initial stiffness schemes are very useful. These situations include problems where a nonassociated flow rule is used or where the zone of plastic yielding is highly localized. This paper surveys the performance of several single‐parameter techniques for accelerating the convergence of the initial stiffness scheme. Some simple but effective modifications to these procedures are also proposed. In particular, a modified version of Thomas' acceleration scheme is developed which has a good rate of convergence. Previously published results on the performance of various acceleration algorithms for initial stiffness iteration are rare and have been restricted to relatively simple yield criteria and simple problems. In this study, detailed numerical results are presented for the expansion of a thick cylinder, the collapse of a rigid strip footing, and the failure of a vertical cut. These analyses use the Mohr–Coulomb and Tresca yield criteria which are popular in soil mechanics. Copyright © 2000 John Wiley & Sons, Ltd.