Solution of quasi‐static large‐strain problems by the material point method     

L. Beuth  P. A. Vermeer 《国际地质力学数值与分析法杂志》2011,35(13):1451-1465

A procedure for solving quasi‐static large‐strain problems by the material point method is presented. Owing to the Lagrangian–Eulerian features of the method, problems associated with excessive mesh distortions that develop in the Lagrangian formulations of the finite element method are avoided. Three‐dimensional problems are solved utilizing 15‐noded prismatic and 10‐noded tetrahedral elements with quadratic interpolation functions as well as an implicit integration scheme. An algorithm for exploiting the numerical integration procedure on the computational mesh is proposed. Several numerical examples are shown. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

辐射方案对GRAPES全球模式的影响   总被引：10，自引：4，他引：6  

杨学胜  沈元芳  徐国强 《大气科学》2009,33(3):593-599

辐射过程是中期数值天气预报模式和气候模式中一个重要的物理过程。本文针对GRAPES (Global/Regional Assimilation and PrEdiction System) 全球模式非静力半隐式半拉格朗日的特点, 进行了不同辐射方案的敏感性比较试验, 以期找到一个适合于GRAPES全球模式的辐射方案。用到的方案包括WRF模式中的RRTM长波辐射方案、 NASA的Goddard短波辐射方案以及ECMWF的长短波辐射方案。试验结果表明, ECMWF的辐射方案更适合于GRAPES模式, 30天平均的北半球500 hPa高度场距平相关系数较WRF的辐射方案提高了近2%。另外, 针对实际运行中存在的辐射计算不稳定问题, 对辐射方案进行了气压订正和温度订正, 并对2005年7月影响我国的海棠台风进行了模拟试验, 气压订正后的120小时模式预报在对副热带高压和台风的描述上更接近于分析, 可用预报时效也有了明显提高。  相似文献   

多尺度大气数值预报的技术进展   总被引：2，自引：1，他引：1       下载免费PDF全文

彭新东  李兴良 《应用气象学报》2010,21(2):129-138

随着计算机和计算理论的发展,数值模式正在向全球化、精细化发展,以适应多尺度、多目的应用的要求,从而模糊了大气环流模式和中尺度数值模式的界限,其主要手段就是改进模式动力框架、离散化手段、计算方法、物理过程的普适性。在实际应用中如何提高模式的多尺度计算性能则是问题关键。该文从模式球面坐标系、网格构造、离散化方法的动力特性、垂直坐标和地形处理以及对物理过程的要求出发,探讨分析多尺度大气数值模式的特点:全球/区域可选、非静力近似、具有良好的频散关系和详细的物理过程,垂直高度坐标和"剪切"地形对多尺度通用模式的改良十分重要。除上述特点外,模式所采用的计算方法也应该最大限度地描述大气动力过程特性,采用高性能计算方案有利于多尺度预报。结合当前多尺度预报的国际研究热点和开发前沿,探讨我国新一代多尺度数值预报系统GRAPES的进一步发展及改进方向。  相似文献   

B‐bar based algorithm applied to meshfree numerical schemes to solve unconfined seepage problems through porous media          下载免费PDF全文

Pedro Navas  Susana Lpez‐Querol  Rena C. Yu  Bo Li 《国际地质力学数值与分析法杂志》2016,40(6):962-984

The object of this work is to establish a meshfree framework for solving coupled, steady and transient problems for unconfined seepage through porous media. The Biot's equations are formulated in displacements (or u − w) assuming an elastic solid skeleton. The free surface location and its evolution in time are obtained by interpolation of pore water pressures throughout the domain. Shape functions based on the principle of local maximum entropy are chosen for the meshfree approximation schemes. In order to avoid the locking involved in the fluid phase of the porous media, a B‐bar based algorithm is devised to compute the average volumetric strain in a patch composed of various integration points. The efficiency of such an implementation for one phase problems is shown through the Benchmark problem, Cook's membrane loaded by a distributive shear load. The proposed methodology is firstly applied to various classical examples in unconfined steady seepage problems through earth dams, then to the dynamic consolidation of a soil column. The results obtained for both problems are quite satisfactory and demonstrate the feasibility of the proposed method in solving coupled problems in porous media. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

Incompressible or nearly incompressible soil dynamic behaviour—a new staggered algorithm to circumvent restrictions of mixed formulation     

Maosong Huang  Shiming Wu  O. C. Zienkiewicz 《Soil Dynamics and Earthquake Engineering》2001,21(2)

The restriction of the Babuska–Brezzi stability criteria for the mixed formulation requires non-uniform interpolation of displacement and pore-pressure variables, which complicates the coding in the finite element analysis. This will certainly preclude the use of many otherwise useful elements. A possible solution to overcome this problem introduced earlier by the authors requires a semi-explicit form. However, the effectiveness of the semi-explicit procedure is dependent on the time step length used, which presents difficulties for transient analysis. In this paper, an unconditionally stable staggered implicit-implicit algorithm is developed, which is generalised from the semi-explicit form. For non-uniform mesh configurations, the results indicate a significant improvement compared with those obtained by using the semi-explicit procedure.  相似文献   

Modelling Two-Phase Incompressible Flow in Porous Media Using Mixed Hybrid and Discontinuous Finite Elements     

D. Nayagum  G. Schäfer  R. Mosé 《Computational Geosciences》2004,8(1):49-73

In this paper, we present a numerical model for simulating two-phase (oil–water and air–water) incompressible and immiscible flow in porous media. The mathematical model which is based on a fractional flow formulation is formed of two nonlinear partial differential equations: a mean pressure equation and a water saturation equation. These two equations can be solved in a sequential manner. Two numerical methods are used to discretize the equations of the two-phase flow model: mixed hybrid finite elements are used to treat the pressure equation, h-based Richards' equation and the diffusion term in the saturation equation, the advection term in the saturation equation is treated with the discontinuous finite elements. We propose a better way to calculate the nonlinear coefficients contained in our equations on each element of the discretized domain. In heterogeneous porous media, the saturation becomes discontinuous at the interface between two porous media. We show in this paper how to use the capillary pressure–saturation relationship in order to handle the saturation jump in the mixed hybrid finite element method. The two-phase flow simulator is verified against analytical solutions for some flow problems treated by other authors.  相似文献   

Formulation of discontinuous deformation analysis (DDA) — an implicit discrete element model for block systems     

Lanru Jing 《Engineering Geology》1998,49(3-4):371-381

Continuum-based numerical methods have played a leading role in the numerical solution of problems in rock mechanics and engineering geology. However, for fractured rocks, a continuum assumption often leads to difficult parameters to define and over-simplified geometry to be realistic. In such case, discrete representations of fractures and individual blocks must be adopted. In this paper, a newly emerged member in the family of discrete element methods (DEM), the discontinuous deformation analysis (DDA), is presented, including its variational principle, governing equations, solution techniques and contact representation and detection algorithms. Its relative advantages and shortcomings are compared with the explicit distinct element method and the finite element method. An example of the analysis of tunnel stability is provided to demonstrate the capability of this new method.  相似文献   

Stretch-free generalized normal moveout correction     

Jorge H. Faccipieri  Tiago A. Coimbra  Rodrigo Bloot 《Geophysical Prospecting》2019,67(1):52-68

The effective application of normal moveout correction processes mainly depends on four factors: the chosen traveltime approximation, the stretching associated with the given traveltime, crossing events and phase changes, the last two being inherent to the seismic data. In this context, we conduct a quantitative analysis on stretching considering a general traveltime expression depending on half-offset and midpoint coordinates. Through this analysis, we propose a mathematically proven procedure to eliminate stretching, which can be applied to any traveltime approximation. The proposed method is applied to synthetic and real data sets, considering different traveltime approximations and achieved complete elimination of stretching.  相似文献   

基于Comsol的饱和多孔介质动力方程的数值模拟及应用   总被引：1，自引：0，他引：1       下载免费PDF全文

杨贝贝  李伟华  赵成刚  李梦姿 《地震工程学报》2017,39(2):321-328

两相饱和多孔介质的动力响应问题在地震工程领域具有重要的研究意义,由于涉及到固相和液相的动力耦合,使得该问题的求解尤为复杂。本文利用Comsol在求解多场耦合问题上的优点,针对Biot饱和多孔介质u-U耦合形式下的波动方程特征,经过一系列微分算子运算和矩阵变换得到导数形式下的波动方程,基于Comsol Multiphysics提供的广义偏微分方程模式对变形后的波动方程进行求解,并把改进后的无限元边界应用到无限域动力问题的模拟中。通过与饱和多孔介质动力响应的解析解进行对比,验证模型求解技术的可行性和正确性,并在此基础上讨论饱和土地基中空沟隔振效果与饱和土体参数孔隙率、泊松比的关系。通过研究分析,可以为饱和土地基中空沟隔振设计提供一些有价值的参考。  相似文献   

沿程非均匀分布植物区对海啸波消减的数值研究     

蒋昌波  徐进  邓斌  陈杰  屈科 《海洋工程》2020,38(3):1-11

为研究近岸植物的非均匀分布对海啸波的消减作用,采用非静压单相流模型计算了5种不同密度分布植物与孤立波的相互作用,分析不同密度分布植物对海啸的消波特性。结果表明在植物总株数相同时,具有不同分布密度植物的消波效果相近,透射波波高差小于3%、最大流速差小于5.4%,但是不同植物密度分布情况下植物区前的反射系数有所不同,反射系数最大差为80.2%。此外,所有分布的植物区对波高和流速的透射系数随着波高和植物区长度变化的趋势是一样的,在不同入射波高和植物区长度情况下5种分布均呈现出相似的沿程衰减机制。  相似文献