Numerical manifold method for dynamic nonlinear analysis of saturated porous media     

H. W. Zhang  L. Zhou 《国际地质力学数值与分析法杂志》2006,30(9):927-951

It is well known that the Babuska–Brezzi stability criterion or the Zienkiewicz–Taylor patch test precludes the use of the finite elements with the same low order of interpolation for displacement and pore pressure in the nearly incompressible and undrained cases, unless some stabilization techniques are introduced for dynamic analysis of saturated porous medium where coupling occurs between the displacement of solid skeleton and pore pressure. The numerical manifold method (NMM), where the interpolation of displacement and pressure can be determined independently in an element for the solution of u–p formulation, is derived based on triangular mesh for the requirement of high accurate calculations from practical applications in the dynamic analysis of saturated porous materials. The matrices of equilibrium equations for the second‐order displacement and the first‐order pressure manifold method are given in detail for program coding. By close comparison with widely used finite element method, the NMM presents good stability for the coupling problems, particularly in the nearly incompressible and undrained cases. Numerical examples are given to illustrate the validity and stability of the manifold element developed. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

Restricted Three-Body Problem: An Approximation of its General Solution – Part One – the Manifold of Symmetric Periodic Solutions     

Kostantinos Papadakis  Constantine Goudas 《Astrophysics and Space Science》2006,305(2):99-124

This paper gives the results of a programme attempting to exploit ‘la seule bréche’ (Poincaré, 1892, p. 82) of non-integrable systems, namely to develop an approximate general solution for the three out of its four component-solutions of the planar restricted three-body problem. This is accomplished by computing a large number of families of ‘solutions précieuses’ (periodic solutions) covering densely the space of initial conditions of this problem. More specifically, we calculated numerically and only for μ = 0.4, all families of symmetric periodic solutions (1st component of the general solution) existing in the domain D:(x ₀ ∊ [−2,2],C ∊ [−2,5]) of the (x ₀, C) space and consisting of symmetric solutions re-entering after 1 up to 50 revolutions (see graph in Fig. 4). Then we tested the parts of the domain D that is void of such families and established that they belong to the category of escape motions (2nd component of the general solution). The approximation of the 3rd component (asymmetric solutions) we shall present in a future publication. The 4th component of the general solution of the problem, namely the one consisting of the bounded non-periodic solutions, is considered as approximated by those of the 1st or the 2nd component on account of the `Last Geometric Theorem of Poincaré' (Birkhoff, 1913). The results obtained provoked interest to repeat the same work inside the larger closed domain D:(x ₀ ∊ [−6,2], C ∊ [−5,5]) and the results are presented in Fig. 15. A test run of the programme developed led to reproduction of the results presented by Hénon (1965) with better accuracy and many additional families not included in the sited paper. Pointer directions construed from the main body of results led to the definition of useful concepts of the basic family of order n, n = 1, 2,… and the completeness criterion of the solution inside a compact sub-domain of the (x ₀, C) space. The same results inspired the ‘partition theorem’, which conjectures the possibility of partitioning an initial conditions domain D into a finite set of sub-domains D i that fulfill the completeness criterion and allow complete approximation of the general solution of this problem by computing a relatively small number of family curves. The numerical results of this project include a large number of families that were computed in detail covering their natural termination, the morphology, and stability of their member solutions. Zooming into sub-domains of D permitted clear presentation of the families of symmetric solutions contained in them. Such zooming was made for various values of the parameter N, which defines the re-entrance revolutions number, which was selected to be from 50 to 500. The areas generating escape solutions have being investigated. In Appendix A we present families of symmetric solutions terminating at asymptotic solutions, and in Appendix B the morphology of large period symmetric solutions though examples of orbits that re-enter after from 8 to 500 revolutions. The paper concludes that approximations of the general solution of the planar restricted problem is possible and presents such approximations, only for some sub-domains that fulfill the completeness criterion, on the basis of sufficiently large number of families.  相似文献   

用遗传算法解算病态方程   总被引：7，自引：1，他引：6  

曾群意  欧吉坤 《大地测量与地球动力学》2003,23(3):93-97

对应用遗传算法解决病态方程问题进行了探讨。利用拟合法而不是通过法方程求解参数，从而避免了法方程系数求逆，使病态方程的解答有了较好的结果。通过模拟计算并和其他方法进行比较，证明该方法是可行的和有效的。  相似文献   

Fast iterative solution of large undrained soil‐structure interaction problems     

Kok‐Kwang Phoon  Swee‐Huat Chan  Kim‐Chuan Toh  Fook‐Hou Lee 《国际地质力学数值与分析法杂志》2003,27(3):159-181

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.  相似文献   

On the Assessment of Failure in Slope Stability Analysis by the Finite Element Method   总被引：3，自引：0，他引：3  

H. Zheng  D. F. Liu  C. G. Li 《Rock Mechanics and Rock Engineering》2008,41(4):629-639

  相似文献   

多传感器观测下带乘性噪声系统的逆向滤波与反褶积融合算法     

褚东升  李玉全  梁猛  张玲 《中国海洋大学学报(自然科学版)》2002,32(5):837-842

针对多传感器观测环境下带乘性噪声系统的逆向最优滤波与反褶积融合估计问题 ,本文提出了 1种基于极大似然准则的最优融合算法。该算法中各单传感器间并行计算 ,并且融合中心与单传感器处理中心间无反向通讯 ,因而执行效率较高。仿真表明 ,该融合算法产生的逆向滤波与反褶积比单传感器处理结果有较明显提高  相似文献   

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

杨玉霞  刘新国 《中国海洋大学学报(自然科学版)》2003,33(2):317-323

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

Investigating spatio-temporal change in spawning activity by Atlantic mackerel between 1977 and 1998 using generalized additive models   总被引：1，自引：3，他引：1  

Beare  D. J.; Reid  D. G. 《ICES Journal of Marine Science》2002,59(4):711-724

  相似文献   

一类广义Sine-Gordon方程的概周期解   总被引：4，自引：0，他引：4  

朴大雄  邱汶华 《中国海洋大学学报(自然科学版)》2006,36(6):892-894

用最大值原理和Banach压缩映像原理研究一类广义Sine-Gordon方程的概周期解问题,证明概周期解的存在‖u‖L∞<π2中的唯一性。  相似文献   

关于GAOR、GSSOR的特征值关系式     

李瑞明 《中国海洋大学学报(自然科学版)》1997,(1)

对两种广义迭代方法ＧＡＯＲ和ＧＳＳＯＲ分别导出了其迭代矩阵和Ｊａｃｏｂｉ迭代矩阵的特征值之间的关系式，这些结果推广了已有的结果  相似文献