Preconditioned IDR(s) iterative solver for non‐symmetric linear system associated with FEM analysis of shallow foundation     

H.H.T. Tran  K.C. Toh  K.K. Phoon 《国际地质力学数值与分析法杂志》2013,37(17):2972-2986

Non‐associated flow rule is essential when the popular Mohr–Coulomb model is used to model nonlinear behavior of soil. The global tangent stiffness matrix in nonlinear finite element analysis becomes non‐symmetric when this non‐associated flow rule is applied. Efficient solution of this large‐scale non‐symmetric linear system is of practical importance. The standard Krylov solver for a non‐symmetric solver is Bi‐CGSTAB. The Induced Dimension Reduction [IDR(s)] solver was proposed in the scientific computing literature relatively recently. Numerical studies of a drained strip footing problem on homogenous soil layer show that IDR(s = 6) is more efficient than Bi‐CGSTAB when the preconditioner is the incomplete factorization with zero fill‐in of global stiffness matrix K_ep (ILU(0)‐K_ep). Iteration time is reduced by 40% by using IDR(s = 6) with ILU(0)‐K_ep. To further reduce computational cost, the global stiffness matrix K_ep is divided into two parts. The first part is the linear elastic stiffness matrix K_e, which is formed only once at the beginning of solution step. The second part is a low‐rank matrix Δ, which is re‐formed at each Newton–Raphson iteration. Numerical studies show that IDR(s = 6) with this ILU(0)‐K_e preconditioner is more time effective than IDR(s = 6) with ILU(0)‐K_ep when the percentage of yielded Gauss points in the mesh is less than 15%. The total computation time is reduced by 60% when all the recommended optimizing methods are used. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

Analytical layer‐element method for non‐axisymmetric consolidation of multilayered soils     

Zhi Yong Ai  Wen Ze Zeng 《国际地质力学数值与分析法杂志》2012,36(5):533-545

A numerically efficient and stable method is developed to analyze Biot's consolidation of multilayered soils subjected to non‐axisymmetric loading in arbitrary depth. By the application of a Laplace–Hankel transform and a Fourier expansion, the governing equations are solved analytically. Then, the analytical layer‐element (i.e. a symmetric stiffness matrix) describing the relationship between generalized displacements and stresses of a layer is exactly derived in the transformed domain. Considering the continuity conditions between adjacent layers, the global stiffness matrix of multilayered soils is obtained by assembling the inter‐related layer‐elements. Once the solution in the Laplace–Hankel transformed domain that satisfies the boundary conditions has been obtained, the actual solution can be derived by the inversion of the Laplace–Hankel transform. Finally, numerical examples are presented to verify the theory and to study the influence of the layered soil properties and time history on the consolidation behavior. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

锡石的晶体结构和形态的关系          下载免费PDF全文

柴钢  王文魁 《岩石矿物学杂志》1988,7(4):374-381

本文对不同成因类型的锡石晶体形态进行了研究,并从PBC(周期性键链)理论出发,讨论了锡石晶体形态与晶体结构的关系。笔者认为,锡石晶体具有三种F面,即s{111}、e{101}、m{110},与Hartman对同样结构类型的金红石晶面性质划分有所差异。据锡石晶面结构性质所推导的理想晶体形态图,与晶体测量所得到的天然晶体形态图总体上相吻合。  相似文献   

北京地区气温诊断分析   总被引：1，自引：0，他引：1  

谢庄  李慧 《大气科学》1994,18(Z1):820-825

本文讨论了北京地区百余年气温的变化特征，依据加法模型理论，提出了一种气候诊断和预测方法─—逐步回归多重因子方法，其优点在于将序列的周期因子、谱成分及外界气候振动因子同时引入回归模型中。这一方法，对单序列的气候诊断和预测具有明显效果。  相似文献   

A STUDY ON THE VARIATIONS OF CIRCULATION PATTERN IN THE NORTHERN HEMISPHERE AND CLIMATE CHANGE IN CHINA FOR THE RECENT 100 YEARS          下载免费PDF全文

Zhang Jiji  Gao Xuejie 《Acta Meteorologica Sinica》1993,7(2):194-198

The time series of annual frequency of the circulation patterns W,C,E in the Northern Hemisphere and annual mean temperature of Beijing (T_B),Shanghai (T_S) and the whole country (T_C) in the period 1891-1990 are expanded as a long-term tendency,a periodic function and a random function.The coherence between circulation and temperature is calculated and analyzed.The abrupt change of circulation and temperature in the last 100 years is revealed with the methods of the accumulated anomalous curves,the signal/noise ratio as well as the Mann-Kendall test.The correlations between anomalous circulation and climate are discussed with some explanations.  相似文献   

农田土壤水分的随机模拟和预报   总被引：3，自引：0，他引：3  

申双和  周英 《南京气象学院学报》1993,16(3):324-328

  相似文献   

区域化变量的空间聚类分析     

陈明 《吉林大学学报(地球科学版)》1995,(2)

传统的聚类分析方法利用数据对平均值的偏差，研究变量之间的相互关系。当它应用于区域化变量的分类时，因缺乏对区域化变量的空间结构的考虑，对合理解释区域化变量之间的相互关系不利。空间聚类分析方法比较好地解决了这个问题。本文介绍了该方法的基本原理及其算法，并说明了它在非条件模拟及在某矿区中的应用情况  相似文献   

新的太阳系起源假说——再论巨星依次分离说     

牛广林 《吉林大学学报(地球科学版)》1987,(3)

本文以太阳系星球的物质状态、形态面貌、体积、物质密度、运动速度、运动方向和角动量等作为根据,研究它们在宇宙空间上的分布特点、相互关系、内在联系以及规律的基础上,提出太阳系(乃至宇宙星球)是在46亿年前,宇宙中一个巨大灼热的星球,在高速度左旋旋转(公转和自转)运动的过程中,从尾部分离出来的物质,经过……恒星、行星和卫星依次波状旋转分离运动产生的。  相似文献   

Global Bifurcations of Periodic solutions of the Hill Lunar Problem*     

A.J. Maciejewski  S.M. Rybicki 《Celestial Mechanics and Dynamical Astronomy》2001,81(4):279-297

  相似文献   

Unified Approach to Photon-Counting Microlaser Rangers, Transponders, and Altimeters   总被引：1，自引：0，他引：1  

John J. Degnan 《Surveys in Geophysics》2001,22(5-6):431-447

Unlike current manned systems, NASA's next generation SLR2000 Satellite Laser Ranging (SLR) station is fully autonomous, eye-safe, relatively compact and inexpensive, and, during daytime tracking, operates at signal-to-noise ratios several orders of magnitude below unity. Tiny, passivelyQ-switched microlasers generate ultra-short pulses with output energies on the order of 100 J at few kHz rates to achieve mm-levelranging precision to satellite altitudesof 20,000 km. Special ranging receivers, combined with Poisson statistical analysis of the received photon distribution, enable the system to rapidly and reliably identify and extract the single photon laser echoes from the solar background. The enhanced rate of return, combined with a uniform signal strength, can actually drive down both systematic and random range errors. The new SLR2000 technology has already spawned exciting new applications. Compact microlaser altimeters, capable of mapping the surface of a planet or other celestial body at multikilohertz rates, is one such application, and a high altitude, airborne version is currently being developed under NASA's Instrument Incubator Program. Interplanetary microlaser transponders would be capable of performing decimeter ranging or subnanosecond time transfer to spacecraft throughout the inner Solar System, resulting in improved knowledge of planetary motions and librations and enhanced General Relativity experiments.  相似文献