Numerical manifold method for dynamic nonlinear analysis of saturated porous media

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

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.     

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.     

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

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

On Szebehely's problem for holonomic systems involving generalized potential functions

We consider the problem of finding the generalized potential function V = U _i(q ¹, q ²,..., q ⁿ)q ⁱ + U(q ¹, q ²,...;q ⁿ) compatible with prescribed dynamical trajectories of a holonomic system. We obtain conditions necessary for the existence of solutions to the problem: these can be cast into a system of n – 1 first order nonlinear partial differential equations in the unknown functions U ₁, U ₂,...;, U _n, U. In particular we study dynamical systems with two degrees of freedom. Using adapted coordinates on the configuration manifold M ₂ we obtain, for potential function U(q ¹, q ²), a classic first kind of Abel ordinary differential equation. Moreover, we show that, in special cases of dynamical interest, such an equation can be solved by quadrature. In particular we establish, for ordinary potential functions, a classical formula obtained in different way by Joukowsky for a particle moving on a surface.Work performed with the support of the Gruppo Nazionale di Fisica Matematica (G.N.F.M.) of the Italian National Research Council.     

一类广义Hamilton系统的辛格式

通过引人泊松括号，分析了无限维Hamilton的性质，并将其推广到广义Hamilton系统，且从理论和实用角度讨论了这类广义Hamilton系统的辛格式构造问题，从而为辛几何算法在一般的时间发展方程的数值求解提供新的具体途径。     

A response‐based decoupling criterion for multiply‐supported secondary systems

It is often infeasible to carry out coupled analyses of multiply‐supported secondary systems for earthquake excitations. ‘Approximate’ decoupled analyses are then resorted to, unless the response errors due to those are significantly high. This study proposes a decoupling criterion to identify such cases where these errors are likely to be larger than an acceptable level. The proposed criterion is based on the errors in the primary system response due to decoupling and has been obtained by assuming (i) the input excitation to be an ideal white noise process, (ii) cross‐modal correlation to be negligible, and (iii) the combined system to be classically damped. It uses the modal properties of the undamped combined system, and therefore, a perturbation approach has been formulated to determine the combined system properties in case of light to moderately heavy secondary systems. A numerical study has been carried out to illustrate the accuracy achieved with the proposed perturbation formulation. The proposed decoupling criterion has been validated with the help of two example primary‐secondary systems and four example excitation processes. Copyright © 2002 John Wiley & Sons, Ltd.     

Nonlinear earthquake response analysis and energy calculation for seismic slit shear wall structures

Based on the concept of structural passive control, a new type of slit shear wall, with improved seismic performance when compared to an ordinary solid shear wall, was proposed by the authors in 1996. The idea has been verified by a series of pseudo-static and dynamic tests. In this paper a macro numerical model is developed for the wall element and the energy dissipation device. Then, nonlinear time history analysis is carried out for a 10-story slit shear wall model tested on a shaking table. Furthermore, the seismic input energy and the individual energy dissipated by the components are calculated by a method based on Newmark-β assumptions for this shear wall model, and the advantages of this shear wall are further demonstrated by the calculation results from the viewpoint of energy. Finally, according to the seismic damage criterion on the basis of plastic accumulative energy and maximum response, the optimal analysis is carried out to select design parameters for the energy dissipation device.     

小浪底工程岩体力学参数研究

本文系统研究了小浪底工程的岩体的强度特征,采用了室内试验,野外结构面统计,工程地质类比法、模型试验法和Hoek Brown的经验法,估算了工程岩体的力学指标,指出了上述各种方法中所存在的不足,认为在工程实践中,应尽可能使用多种估计方法,交叉校核所得参数,以提高参数的可靠性,保证工程的安全,节约投资。