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.     

交通流量分析与实时最优控制

通过将车流量的增大或减小转化为路长权重的变化。将交通流量的动态问题转化为静态问题，用解决最短路问题的Dijkstra方法，给出交通流量实时最优控制的可行性模型及其有效算法。     

模糊划分矩阵在岩土参数概率分布中的应用

讨论如何在小样本条件下用已有的过程经验与试验资料确定岩土参数概率分布,用模糊划分矩阵与BAYES方法相结合,给出由小样本试验数据确定岩土参数的概率分布。     

量子化学方法研究氯代苯酚类化合物对水生物毒性的构效关系

对20个氯代苯酚类化合物进行了CNDO/2量子化学计算,讨论了化合物电子结构与其对水生物翻车鱼(Bluegill)、花鳉鱼(Guppy)、和虾(Shrimp)半致死量之间的关系,分别获得了表示其构效关系的三个线性方程,它们的显著性均远高于α=0.01水平。结果表明:氯代苯酚化合物苯环上碳原子的兀电荷之和(∑Qπ_R)越大苯环上碳原子的净电荷之和(∑Q_R)越大,或LUMO轨道能(E_(LUMO))越低,均使化合物对Bluegill和Guppy水生物的毒性增大;而化合物苯环上羟基邻位的碳原子的兀电荷布居(P_3)越大,同时HOMO轨道能E_(HOMO)越高,则化合物对Shrimp毒性也越大。据此,可预测氯代苯酚类化合物对上述三种水生物的毒性。     

山西省石炭系与二叠系界线

１８８２年李希霍芬在太原西山建立石炭系的太原系含煤建造以来，石炭系与二叠系的分界意见，归纳起来有三：①以Ｐｓｅｕｄｏｓｃｈｗａｇｅｒｉｎａ（ｓ．１）带之底为界，下面为石炭系，上面为二叠系；②以晋东南的灰白色“佳祥砂岩”与太原西山黄绿色“骆驼脖子砂岩”对比，砂岩底面为石炭系与二叠系的分界线，下面是石炭系山西组，上面是二叠系石盒子组，山西组的地层时代定为晚石炭世；③以晋东南灰白色第三砂岩（陵川砂岩）与太原西山“北岔沟砂岩”对比，砂岩底面为石炭系与二叠系的分界线，下面是石炭系太原组，上面是二叠系山西组，山西组的地层时代为早二叠世。笔者认为石炭系与二叠系的划分，以第③方案较适宜。     

A MESO-β SCALE SIMULATION OF THE EFFECTS OF BOREAL FOREST ECOSYSTEM ON THE LOWER ATMOSPHERE

Based on the Intensive Field Campaign(IFC-1)data of Boreal Ecosystem-Atmosphere Study(BOREAS).a three-dimensional meso-β scale model is used to simulate the effect of boreal forests onthe lower atmosphere.A fine horizontal resolution of 2 km×2 km is used in order to distinguish thevegetative heterogeneity in the boreal region.A total of 20×25 grid points cover the entire sub-modeling area in BOREAS' South Study Area(SSA).The ecosystem types and their coverage ineach grid square are extracted from the North American Land Cover Characteristics Data Base(NALCCD)generated by the U.S.Geographical Survey(USGS)and the University of Nebraska-Lincoln(UNL).The topography of the study area is taken from the Digital Elevation Map(DEM)of USGS.The model outputs include the components of the energy balance budget within the canopyand at the ground.the turbulence parameters in the atmospheric boundary layer and the wind.temperature and humidity profiles extending up to a height of 1500 m.In addition to the fine timeand spatial step,the unique feature of the present model is the incorporation of both dynamic andbiological effects of the Boreal forest into the model parameterization scheme.The model resultscompare favorably with BOREAS' IFC-1 data in 1994 when the forest was in the luxuriant growingperiod.     

Numerical calculations in the orbital determination of an artificial satellite for a long arc

Numerical methods are usually used for the computation of ephemerides with perturbations for the precise orbital determination of an artificial satellite. But their numerical stability will be encountered in a long arc. In this case the use the improved Encke special perturbation methods has been suggested. The results of this paper show that Encke's method does indeed have a certain effectiveness, but cannot yet completely resolve the numerical stability, and the more efficient method is to use the energy integral or its variational relation to control the growth of the along-track error in general numerical calculations so that the aim of stabilization can be achieved.     

A stepwise method for determining the number of component distributions in a mixture

This paper considers the problem of estimatingm, the number of components in a finite mixture of distributions from a parametric family. A step-up procedure using the bootstrap method is proposed. Some properties of the procedure are illustrated with simulation studies. An example of the method, applied to orientation of beach clasts, is given.     

Notes on Hori Method for Canonical Systems

In this paper a slightly different approach is proposed for the process of determining the functions S _m and H _m ^* of the algorithm of the canonical version of Hori method. This process will be referred to as integration theory of the mth order equation of the method. It will be shown that the ordinary differential equation with an auxiliary parameter t ^* as independent variable, introduced through Hori auxiliary system, can be replaced by a partial differential equation in the time t. In this way, the mth order equation of the algorithm assumes a form very similar to the one of other perturbation methods. In virtue of this new approach of the integration theory for Hori method, Lagrange's variational equations introduced by Sessin are revised. As an example, the Duffing equation is solved through this new approach.     

电测深ρ＿s差值法的应用效果

用电测深研究地电断面、寻找隐伏多金属矿体时，以往采用ρ＿ｓ（或η＿ｓ）等值线断面图进行地质推断解释，因受浅源信号的干扰，效果差。本文提出了采用ρ＿ｓ（或η＿ｓ）差值法，能排除浅源干扰信号，反映出地电断面的相对变化，突出了深部微弱信息，取得了较好的地质解释效果。