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.     

Numerical modelling of the impact crater depth-diameter dependence in an acoustically fluidized target

2D numerical modelling of impact cratering has been utilized to quantify an important depth-diameter relationship for different crater morphologies, simple and complex. It is generally accepted that the final crater shape is the result of a gravity-driven collapse of the transient crater, which is formed immediately after the impact. Numerical models allow a quantification of the formation of simple craters, which are bowl-shaped depressions with a lens of rock debris inside, and complex craters, which are characterized by a structural uplift. The computation of the cratering process starts with the first contact of the impactor and the planetary surface and ends with the morphology of the final crater. Using different rheological models for the sub-crater rocks, we quantify the influence on crater mechanics. To explain the formation of complex craters in accordance to the threshold diameter between simple and complex craters, we utilize the Acoustic Fluidization model. We carried out a series of simulations over a broad parameter range with the goal to fit the observed depth/diameter relationships as well as the observed threshold diameters on the Moon, Earth and Venus.     

基于统一强度理论的地基临界荷载公式

基于统一强度理论,按照临界荷载公式的推导方法,得到地基临界荷载的统一解公式。Mohr-Coulomb强度理论计算的结果为其特例。把计算结果与《建筑地基基础设计规范》中的结果进行了对比分析,认为当不同程度地考虑中间主应力σ2的影响时,可以有效地发挥地基的强度潜能。并通过实例计算了统一强度理论参数b取不同值时地基承载力特征值,并进行了对比分析。     

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.     

Comments on ‘A comparative study of ANN and neuro-fuzzy for the prediction of dynamic constant of rockmass’

Singh et al (2005) examined the potential of the ANN and neuro-fuzzy systems application for the prediction of dynamic constant of rockmass. However, the model proposed by them has some drawbacks according to fuzzy logic principles. This discussion will focus on the main fuzzy logic principles which authors and potential readers should take into consideration.     

Discussion on “Experimental study on fracture toughness and tensile strength of a clay” [Engineering Geology 94 (2007) 64–75]

The paper [Wang, J.-J., Zhu, J.-G., Chiu, C.F., Zhang, H., 2007. Experimental study on fracture toughness and tensile strength of a clay. Engineering Geology 94, 64–75.] focuses on two important fracture parameters of soils: tensile strength and fracture toughness. These parameters control the behaviour of soils in a wide range of situations, from the design of a simple footing to much complicated fracture behaviour of clay liners or covers. The authors have done extensive laboratory work to determine these two parameters and their laborious and complicated experimental program needs praise. However, some of the points raised in their conclusions, based on the analysis and comparison with the data from the literature, need to be discussed.     

Shear strength characterization of municipal solid waste at the Suzhou landfill, China

The current practice of slope stability analysis for a municipal solid waste (MSW) landfill usually overlooks the dependence of waste properties on the fill age or embedment depth. Changes in shear strength of MSW as a function of fill age were investigated by performing field and laboratory studies on the Suzhou landfill in China. The field study included sampling from five boreholes advanced to the bottom of the landfill, cone penetration tests and monitoring of pore fluid pressures. Twenty-six borehole samples representative of different fill ages (0 to 13 years) were used to perform drained triaxial compression tests. The field and laboratory study showed that the waste body in the landfill can be sub-divided into several strata corresponding to different ranges of fill age. Each of the waste strata has individual composition and shear strength characteristics. The triaxial test results showed that the MSW samples exhibited a strain-hardening and contractive behavior. As the fill age of the waste increased from 1.7 years to 11 years, the cohesion mobilized at a strain level of 10% was found to decrease from 23.3 kPa to 0 kPa, and the mobilized friction angle at the same strain level increasing from 9.9° to 26°. For a confinement stress level greater than 50 kPa, the shear strength of the recently-placed MSW seemed to be lower than that of the older MSW. This behavior was consistent with the cone penetration test results. The field measurement of pore pressures revealed a perched leachate mound above an intermediate cover of soils and a substantial leachate mound near the bottom of the landfill. The measurements of shear strength properties and pore pressures were utilized to assess the slope stability of the Suzhou landfill.     

Weibull statistical analysis of granite bending strength

Summary  This paper describes and discusses the adequacy of Weibull statistical analysis to analyse the bending strength of granite. The experimental results show that strength variability is related with a specific origin of failure. This conclusion is based on analysing the influence of the surface condition (extrinsic defects) on the bending strength results treated by the Weibull statistics. The conclusions drawn from this study have been validated by analysing the results of the critical flaw dimension estimated by applying the linear elastic fracture mechanics (LEFM) formulae. Results obtained from fractographic examination also have been used to describe the location of the origin of the fracture and understand the distribution of defects; i.e., there is a unimodal distribution of defects (intrinsic defects), despite the fact that some outlier values are normally observed in the fractured surfaces. Correspndence: P. M. Amaral, Department of Materials Engineering, Instituto Superior Técnico, Technical University of Lisbon, Av. Rovisco Pais, 1049-001 Lisboa, Portugal