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.     

饱和粉砂不稳定性的试验研究

通过对净砂和级配良好粉砂(含10 %粉土)进行一系列三轴固结不排水试验（CU）,研究了粉土、孔隙比和围压对饱和粉砂不稳定性的影响。试验结果表明,净砂与粉砂在不排水剪切条件下均会出现应变软化现象（即不稳定性）。同一围压下脆性指数(IB)随孔隙比增加,但不稳定线的应力比随孔隙比增加而减小。引用等效粒间孔隙比(ege)后,净砂和粉砂在ege-ln p?平面上拥有基本相同的临界状态线。在临界状态理论及等效粒间孔隙比的基础上,提出在同一修正状态参数(?ge)下净砂和级配良好粉砂有相似的不稳定性。     

SO4^2-改性铁铝复合层柱粘土的制备及其催化性能研究

由共聚法制得的铁铝复合交联剂交换钠基蒙脱石层间的水合钠离子,再经SO42-改性处理制备了铁铝复合层柱粘土固体超强酸催化剂.利用xRD、FT-IR和BET法对催化剂的结构进行了表征.以乙酸正丁酯的合成反应为探针反应对催化剂的活性进行了测试.通过正交实验确定了铁铝复合交联粘土SO42-改性最佳工艺条件.研究了催化剂对探针反应的选择性及重复使用情况.结果表明:在由实验确定的最佳改性工艺条件下,乙酸的转化率达92.2%,反应中无副产物产生,催化剂具有较好的重复使用性能.     

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

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

闽南：台湾浅滩渔场变性水团的初步分析

根据聚类方法，对闽南－台湾浅滩渔场变性水团进行划分，本海区共有６个变性水团。讨论各个变性水团在不同季节消长变化规律。     

聚类分析法在浅海水团分析中的应用及黄、东海变性水团的分析

本文引入变性水团概念,用聚类分析确定黄、东海浅海海域变性水团的边界及其变化。根据给定的标准进行聚类,确定了水团的个数、边界及混合区。 由聚类方法所得结果看出,在该海区中有八个变性水团。根据温、盐相对指标,将这些水团分为九种不同的特征。也可以把这些水团划分为三种盐度类型。在温—盐点聚图上,不同变性水团的温盐点,分布于一条曲线附近,它体现出逐级变性的特点。各变性水团的变化范围都很大,这说明该海区中水团变性强烈。在暖季中的增温降盐和冷季中的降温增盐,可认为是整个黄海和东海水体变性的特征。所有季节性水团都经历一个形成与消亡的过程。 文中讨论了变性水团和海洋环境的相互关系。水团变性是热力因素和动力因素共同作用而产生的。通过变性水团的分析,作者认为,在该海区水团的变性中,热力因素起着主要作用。变性水团界线的舌状分布与流向之间的关系是很明显的,可作为分析环流的旁证。最后描述了变性水团和底层中心渔场之间的关系。     

利用Jason-1卫星确定全球平均海面高

介绍了卫星测高数据处理的基本方法,分析研究了卫星测高的主要误差,提出了Jason-1的数据编辑准则,并用其测高数据确定了全球平均海面高,与CLS01模型进行了比较分析.     

基于权变理论的滑坡灾害监测预报思路探讨

文章简要评述了地质灾害基础理论与应用技术发展现状、滑坡灾害多种监测预报判据的利弊。利用综合信息处理决策方法,提出了基于权变理论的滑坡灾害监测预报新思路。分析了滑坡成灾的权变特征、环境因素和决策因素,建立了滑坡灾害预报决策概念模型。进一步探讨了在预报决策中应遵循的动态性及满意性原则,为提高地质灾害监测预报理论的科学性提供了新的理论依据与技术途径。