INVERSE DDA OF A SYSTEM OF BLOCKS IN SPHERE COORDINATES

1　IntroductionDiscontinuousDeformationAnalysis(DDA)developedbyShiin 1 980’s[1 ,2 ] maybeusedtocalculatethedeformationanddisplacementsinamulti blocksystemandtoanalyzetherelationbe tweenforcesanddisplacementsinthesystem .Ifthesurveydisplacementofeverypointordirectionisknown,alldisplacementsanddeformationsofblocks,eventheglide ,theclosureandthestretchoftheboundary planeofblocks ,canbecalculated .TheresultisoptimallyfittedaccordingtotheLeastSquarePrinciple.DDAismainlyusedinrockblockssystems…     

Analysis of seabed instability using element free Galerkin method

Wave-induced seabed instability, either momentary liquefaction or shear failure, is an important topic in ocean and coastal engineering. Many factors, such as seabed properties and wave parameters, affect the seabed instability. A non-dimensional parameter is proposed in this paper to evaluate the occurrence of momentary liquefaction. This parameter includes the properties of the soil and the wave. The determination of the wave-induced liquefaction depth is also suggested based on this non-dimensional parameter. As an example, a two-dimensional seabed with finite thickness is numerically treated with the EFGM meshless method developed early for wave-induced seabed responses. Parametric study is carried out to investigate the effect of wavelength, compressibility of pore fluid, permeability and stiffness of porous media, and variable stiffness with depth on the seabed response with three criteria for liquefaction. It is found that this non-dimensional parameter is a good index for identifying the momentary liquefaction qualitatively, and the criterion of liquefaction with seepage force can be used to predict the deepest liquefaction depth.     

Estimation of Land Subsidence Based on Groundwater Flow Model

This article presents an approach for estimating land subsidence due to withdrawal of groundwater. The proposed method calculates the groundwater seepage in 3-D-condition and calculates the land subsidence one-dimensionally. The governing equation on groundwater seepage is based on the three-dimensional mass conservation law and the principle of effective stress. The land subsidence calculation method is derived based on the following assumptions: (1) displacements occur only in the vertical direction, and (2) in vertical direction the total stresses do not change. The governing equation is solved by numerical method, i.e., finite element method (FEM) in spatial discretization and finite difference method (FDM) in time series discretization. In FEM Galerkin method is adopted and in FDM, lumped matrix method is employed. The proposed method is calibrated via analyzing 1-D consolidation problem and the results are compared with those from Terzaghi's one-dimensional consolidation theory and oedometer test. The proposed method is employed to analyze the consolidation of a soft layer due to withdrawal of groundwater from an aquifer under it. Moreover, this method is also applied to a field case of land subsidence due to groundwater pumping in a gas production field in Japan. The analytical results are compared with the field observed data. The results show that this approach simulates the field case well.     

基于决策树的我国农业数据挖掘分析

决策树是用于分类的常用建模方法。首先对分类的概念和决策树方法分别进行了总体介绍,在此基础上对我国30个省市自治区的乡村劳动力、耕地面积以及农业总产值信息数据进行了挖掘分析,在运用决策树对数据进行分类过程中对连续数据采用聚类分析的方法进行离散化处理,从而避免了原始经验分类方法的主观性。最后,通过上述决策树分类方法,生成我国农业情况的决策树,获得相关空间分类规则,并对其进行分析说明。结果表明,决策树分类方法适合我国农业情况。     

A continuous/discontinuous Galerkin framework for modeling coupled subsurface and surface water flow

We consider conjunctive surface-subsurface flow modeling, where surface water flow is described by the shallow water equations and ground water flow by Richards’ equation for the vadose zone. Coupling between the models is based on the continuity of flux and water pressure. Numerical approximation of the coupled model using the framework of discontinuous Galerkin (DG) methods is formulated. In the subsurface, the local discontinuous Galerkin (LDG) method is used to approximate ground water velocity and hydraulic head; a DG method is also used to approximate surface water velocity and elevation. This approach allows for a weak coupling of the models and the use of different approximating spaces and/or meshes within each regime. A simplified LDG method based on continuous approximations to water head is also described. Numerical results that investigate physical and numerical aspects of surface–subsurface flow modeling are presented. This work was supported by National Science Foundation grant DMS-0411413.     

A coupling of mixed and discontinuous Galerkin finite-element methods for poroelasticity

In this paper, we formulate a finite-element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation model of poroelasticity. We approximate the flow variables by a mixed finite-element space and the displacement by a family of discontinuous Galerkin methods. Theoretical convergence error estimates are derived and, in particular, are shown to be independent of the constrained specific storage coefficient, c _o . This suggests that our proposed algorithm is a potentially effective way to combat locking, or the nonphysical pressure oscillations, which sometimes arise in numerical algorithms for poroelasticity.     

基于DEM的分布式水文模型构建方法

基于 DEM的分布式水文模型是现代水文学同高科技 (如计算机技术、 3S技术等 )相结合的产物 ,是研究变化环境下水文循环与水资源演化规律的理想工具 ,代表了水文模型的最新发展方向。从 DEM的特性出发 ,本文探讨并总结了分布式水文模型的特点、建模思路和模型基本结构框图。在流域离散化方面 ,重点介绍了分布式水文模型常用的三种单元划分方法 ;最后 ,针对分布式水文模型构建问题 ,从“输入模块”、“单元水文模型”、“河网汇流模型”三方面 ,阐述了分布式水文模型微结构构建方法。     

1966年以来华北地区一系列七级大震破裂过程的数值模拟

用非连续变形分析方法(DDA+FEM)数值模拟，在华北地区各地块相互制约的块体系统环境中，地块边界断层上发生1966年邢台地震、 1969年渤海地震、 1975年海城地震、 1976年唐山地震等不同类型破裂模式大震的破裂过程. 数值模拟结果给出每个大震释放的主应力场，最大剪应力变化等值线图，地震前后位移变化矢量图，发震断层滑移随时间变化以及走滑错距和应力降等震源参数. 这些结果与地震的震源机制，用地震波资料研究得到的震源参数,宏观等震线,地表观测的水平位移矢量图基本一致. 其中1969年渤海地震正交破裂模式的结果与宏观等震线及小震分布图像更接近. 1976年唐山地震复杂震源模式与该震早期余震分布图像更相符. 表明用DDA+FEM方法进行数值模拟研究，能较好地模拟地震破裂过程.     

ON THE VERTICAL DISCRETIZATION OF ATMOSPHERIC MODELS-PART Ⅰ:EXISTING PROBLEMS,DATA PROCESSING AND CHOICE OF VERTICAL COORDINATE AND BASIS FUNCTIONS

Discussions are carried out on the vertical discretization of current atmospheric models.It is pointed out that there exist problems in the integration of the hydrostatic equation and the computation of vertical advection,vertical diffusion and so on.Then some possible ways for solving or alleviating them are suggested.Finally,the choice of vertical coordinate and basis functions is discussed.     

The Present Space-Time Motion and Deformation Features of the Northeastern Margin of the Qinghai-Xizang（Tibet） Block and Its Adjacent Area

On the basis of Discontinuous Deformation Analysis (DDA), and considering the moderate intrusion of specific block boundaries to different extents, the first-order block motion model is established for the northeastern margin of Qinghai-Xizang(Tibet) block and the kinematical model for depicting deformation of small regions as well by using GPS observations of three periods (1991, 1999 and 2001). By simulating, we obtained the motion features of the firstorder blocks between the large WWN faults on the sides of the studied region, the distribution features of the principal strain rate field and the inhomogeneous motion features with spacetime of the faults in the northern boundary of the Qinghai-Xizang (Tibet) block.