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.     

概率赋范空间的S-凸和K-凸

文献[1][2]认为M－PN空间是局部凸的。通过定理说明这个结论对t-模T（a,b)≥min(a,b)时成立。又通过反例表明，当t-模型T（a,b)=max(a b-1,0）时无论邻域Nθ（ε，λ）概率有界还是概率半有界结论均不成立。还讨论了Takahashi凸性并对S－凸和K－凸作了比较。     

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.     

Seismic Hazard Analysis — Quo vadis?

The paper is dedicated to the review of methods of seismic hazard analysis currently in use, analyzing the strengths and weaknesses of different approaches. The review is performed from the perspective of a user of the results of seismic hazard analysis for different applications such as the design of critical and general (non-critical) civil infrastructures, technical and financial risk analysis. A set of criteria is developed for and applied to an objective assessment of the capabilities of different analysis methods. It is demonstrated that traditional probabilistic seismic hazard analysis (PSHA) methods have significant deficiencies, thus limiting their practical applications. These deficiencies have their roots in the use of inadequate probabilistic models and insufficient understanding of modern concepts of risk analysis, as have been revealed in some recent large scale studies. These deficiencies result in the lack of ability of a correct treatment of dependencies between physical parameters and finally, in an incorrect treatment of uncertainties. As a consequence, results of PSHA studies have been found to be unrealistic in comparison with empirical information from the real world. The attempt to compensate these problems by a systematic use of expert elicitation has, so far, not resulted in any improvement of the situation. It is also shown that scenario-earthquakes developed by disaggregation from the results of a traditional PSHA may not be conservative with respect to energy conservation and should not be used for the design of critical infrastructures without validation. Because the assessment of technical as well as of financial risks associated with potential damages of earthquakes need a risk analysis, current method is based on a probabilistic approach with its unsolved deficiencies.

Traditional deterministic or scenario-based seismic hazard analysis methods provide a reliable and in general robust design basis for applications such as the design of critical infrastructures, especially with systematic sensitivity analyses based on validated phenomenological models. Deterministic seismic hazard analysis incorporates uncertainties in the safety factors. These factors are derived from experience as well as from expert judgment. Deterministic methods associated with high safety factors may lead to too conservative results, especially if applied for generally short-lived civil structures. Scenarios used in deterministic seismic hazard analysis have a clear physical basis. They are related to seismic sources discovered by geological, geomorphologic, geodetic and seismological investigations or derived from historical references. Scenario-based methods can be expanded for risk analysis applications with an extended data analysis providing the frequency of seismic events. Such an extension provides a better informed risk model that is suitable for risk-informed decision making.     

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

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

长期非平稳荷载调谐质量减振阻尼器优化设计

研究外荷载为长期非平稳随机过程。考虑长期荷载的特性 ,采用 1个概率谱密度函数来反映长期非平稳随机荷载及其特征 ;概率谱密度函数是基于大量的一般谱密度函数的统计特性获得。以延长结构的抗疲劳使用寿命为目标函数 ,提出了调谐质量减振阻尼器的优化设计方法 ,这在实际工程中有着极为广阔的应用前景。本文旨在从理论上发展长期非平稳随机荷载作用下调谐质量减振阻尼器的优化设计方法 ;文中采用长期波浪实测数据 ,给出了 1个数值算例说明整个设计过程。     

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

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

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

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

中国海相探矿权区块定量评价方法

在研究了我国海相油气探矿权区块基本油气地质条件、勘探程度、技术经济条件等要素和特点的基础上,给出了48个评价参数和归一化取值参照表,建立了探矿权区块评价工作流程、评分参数及计算公式。提出对我国海相勘探层探矿权区块评价应该注意:与国外海相碳酸盐岩评价的差别;与国内陆相碎屑岩评价的差别;成藏组合的划分;评价内容的有效性、不确定性和评价参数的可信度;成败经验的总结和勘探程度的研究。