Generalized plasticity state parameter‐based model for saturated and unsaturated soils. Part 1: Saturated state     

Diego Manzanal  José Antonio Fernández Merodo  Manuel Pastor 《国际地质力学数值与分析法杂志》2011,35(12):1347-1362

The behavior of granular materials is known to depend on its loose or dense nature, which in turns depends both on density and confining pressure. Many models developed in the past require the use of different sets of constitutive parameters for the same material under different confining pressures. The purpose of this paper is to extend a basic generalized plasticity model for sands proposed by Pastor, Zienkiewicz and Chan by modifying the main ingredients of the model flow—rule, loading–unloading discriminating direction and plastic modulus—to include a dependency on the state parameter. The proposed model is tested against the available experimental data on three different sands, using for each of them a single set of material parameters, finding a reasonably good agreement between experiments and predictions. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

Stability analysis of tunnels driven in jointed rocks by means of a homogenized limit analysis approach     

Samir Maghous  José E. B. Winiawer  Vanessa F. Pasa Dutra 《国际地质力学数值与分析法杂志》2014,38(18):2009-2032

Conceived as a potential alternative to the classical design methods employed for analyzing the stability of underground structures driven in jointed rocks, the homogenization approach stems from the heuristic idea that, from a macroscopic point of view, a rock mass cut by a network of joints may be perceived as a homogenized continuum. The strength properties of the latter can be theoretically obtained from the failure conditions of its individual constituents: rock matrix and joint interfaces. At the material level, the limit analysis reasoning is used in the context of homogenization to formulate the homogenized strength criterion of a jointed rock mass in the particular situation of a single set of parallel joints. As it could be expected, the obtained closed‐form expressions show the strength anisotropy induced by joint preferential orientation. The support functions (π functions) associated with the homogenized strength criterion are also determined in both plane strain and three‐dimensional cases. This criterion is then applied to the investigation of stability analysis of a tunnel excavated in a jointed rock mass. Upper bounds estimated of the stability factor are derived from the implementation of the kinematic approach directly on the homogenized underground structure. Finally, the approach is applied to analyze and discuss the collapse of the Pinheiros subway station (São Paulo, Brazil). Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

Three‐dimensional shakedown solutions for anisotropic cohesive‐frictional materials under moving surface loads     

Juan Wang  Hai‐Sui Yu 《国际地质力学数值与分析法杂志》2014,38(4):331-348

Previous work on three‐dimensional shakedown analysis of cohesive‐frictional materials under moving surface loads has been entirely for isotropic materials. As a result, the effects of anisotropy, both elastic and plastic, of soil and pavement materials are ignored. This paper will, for the first time, develop three‐dimensional shakedown solutions to allow for the variation of elastic and plastic material properties with direction. Melan's lower‐bound shakedown theorem is used to derive shakedown solutions. In particular, a generalised, anisotropic Mohr–Coulomb yield criterion and cross‐anisotropic elastic stress fields are utilised to develop anisotropic shakedown solutions. It is found that shakedown solutions for anisotropic materials are dominated by Young's modulus ratio for the cases of subsurface failure and by shear modulus ratio for the cases of surface failure. Plastic anisotropy is mainly controlled by material cohesion ratio, the rise of which increases the shakedown limit until a maximum value is reached. The anisotropic shakedown limit varies with frictional coefficient, and the peak value may not occur for the case of normal loading only. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

Generalized plasticity state parameter‐based model for saturated and unsaturated soils. Part II: Unsaturated soil modeling     

Diego Manzanal  Manuel Pastor  José Antonio Fernández Merodo 《国际地质力学数值与分析法杂志》2011,35(18):1899-1917

The aim of this paper is to extend the generalized plasticity state parameter‐based model presented in part 1 to reproduce the hydro‐mechanical behavior of unsaturated soils. The proposed model is based on two pairs of stress–strain variables and a suitable hardening law taking into account the bonding—debonding effect of suction and degree of saturation. A generalized state parameter for unsaturated state is proposed to reproduce soil behavior using a single set of material parameters. Generalized plasticity gives a suitable framework to reproduce not only monotonic stress path but also cyclic behavior. The hydraulic hysteresis during a drying—wetting cycle and the void ratio effect on the hydraulic behavior is introduced. Comparison between model simulations and a series of experimental data available, both cohesive and granular, are given to illustrate the accuracy of the enhanced generalized plasticity equation. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

A nonstationary extreme value distribution for analysing the cessation of karst spring discharge          下载免费PDF全文

Yan Liu  Yonghong Hao  Yonghui Fan  Tongke Wang  Xueli Huo  Youcun Liu  Tian‐Chyi J. Yeh 《水文研究》2014,28(20):5251-5258

The effects of climate change and population growth in recent decades are leading us to consider their combined and potentially extreme consequences, particularly regarding hydrological processes, which can be modeled using a generalized extreme value (GEV) distribution. Most of the GEV models were based on a stationary assumption for hydrological processes, in contrast to the nonstationary reality due to climate change and human activities. In this paper, we present the nonstationary generalized extreme value (NSGEV) distribution and use it to investigate the risk of Niangziguan Springs discharge decreasing to zero. Rather than assuming the location, scale, and shape parameters to be constant as one might do for a stationary GEV distribution analysis, the NSGEV approach can reflect the dynamic processes by defining the GEV parameters as functions of time. Because most of the GEV model is designed to evaluate maxima (e.g. flooding, represented by positive numbers), and spring discharge cessation is a ?minima’, we deduced an NSGEV model for minima by applying opposite numbers, i.e. negative instead of positive numbers. The results of the model application to Niangziguan Springs showed that the probability of zero discharge at Niangziguan Springs will be 1/80 in 2025, and 1/10 in 2030. After 2025, the rate of decrease in spring discharge will accelerate, and the probability that Niangziguan Springs will cease flowing will dramatically increase. The NSGEV model is a robust method for analysing karst spring discharge. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

Improved algorithms for generalized finite element simulations of three-dimensional hydraulic fracture propagation     

Nathan Shauer  Carlos Armando Duarte 《国际地质力学数值与分析法杂志》2019,43(18):2707-2742

This paper reports improvements to algorithms for the simulation of 3-D hydraulic fracturing with the Generalized Finite Element Method (GFEM). Three optimizations are presented and analyzed. First, an improved initial guess based on solving a 3-D elastic problem with the pressure from the previous step is shown to decrease the number of Newton iterations and increase robustness. Second, an improved methodology to find the time step that leads to fracture propagation is proposed and shown to decrease significantly the number of iterations. Third, reduced computational cost is observed by properly recycling the linear part of the coupled stiffness matrix. Two representative examples are used to analyze these improvements. Additionally, a methodology to include the leak-off term is presented and verified against asymptotic analytical solutions. Conservation of mass is shown to be well satisfied in all examples.  相似文献   

Cracking risk of partially saturated porous media—Part I: Microporoelasticity model     

Bernhard Pichler  Luc Dormieux 《国际地质力学数值与分析法杂志》2010,34(2):135-157

Drying of deformable porous media results in their shrinkage, and it may cause cracking provided that shrinkage deformations are hindered by kinematic constraints. This is the motivation to develop a thermodynamics‐based microporoelasticity model for the assessment of cracking risk in partially saturated porous geomaterials. The study refers to 3D representative volume elements of porous media, including a two‐scale double‐porosity material with a pore network comprising (at the mesoscale) 3D mesocracks in the form of oblate spheroids, and (at the microscale) spherical micropores of different sizes. Surface tensions prevailing in all interfaces between solid, liquid, and gaseous matters are taken into account. To establish a thermodynamics‐based crack propagation criterion for a two‐scale double‐porosity material, the potential energy of the solid is derived, accounting—in particular—for mesocrack geometry changes (main original contribution) and for effective micropore pressures, which depend (due to surface tensions) on the pore radius. Differentiating the potential energy with respect to crack density parameter yields the thermodynamical driving force for crack propagation, which is shown to be governed by an effective macrostrain. It is found that drying‐related stresses in partially saturated mesocracks reduce the cracking risk. The drying‐related effective underpressures in spherical micropores, in turn, result in a tensile eigenstress of the matrix in which the mesocracks are embedded. This way, micropores increase the mesocracking risk. Model application to the assessment of cracking risk during drying of argillite is the topic of the companion paper (Part II). Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

利用Abel-Poisson径向基函数模型化局部重力场     

马志伟  陆洋  涂弋  朱传东  郗慧 《测绘学报》2016,45(9):1019-1027

多种类型高分辨率重力场数据的不断增加,使得在局部范围内精化重力场模型成为了可能。本文采用Abel-Poisson核将重力场量表示成有限个径向基函数线性求和的形式,对局部区域的多种重力场数据进行联合建模。为了提高运算速度,运用了基于自适应精化格网算法的最小均方根误差准则(RMS)来求解径向基函数平均带宽。以南海核心地区为例,联合两种不同类型、不同分辨率的重力场资料(大地水准面起伏6'×6'、重力异常2'×2'),构建了局部区域高分辨率的重力场模型。所建模型表示的重力场参量达到了2'×2'的分辨率,对原始的重力异常数据(2'×2')拟合的符合程度达到±0.8×10-5m/s2。结果表明,利用径向基函数方法进行局部重力场建模,避免了球谐函数建模收敛慢的问题,有效提高了模型表示重力场的分辨率。  相似文献   

南沙群岛海域长腹剑水蚤(Oithona spp.)的种类组成、数量分布及其与环境因子的关系   总被引：1，自引：0，他引：1  

杜飞雁  王亮根  王雪辉  宁加佳  谷阳光  李亚芳 《海洋与湖沼》2016,47(6):1176-1184

中小型桡足类是海洋食物网中重要的中间环节,在碳循环里细菌碳和藻类碳向高营养级转化的过程中扮演重要角色。本文通过分析2013年3—4月、6—7月、9—10月和11—12月南沙群岛海域(111°30′E—112°30′E,4°30′N—11°30′N)的生态环境调查数据,对长腹剑水蚤属(Oithona)的种类组成、数量分布及其与环境因子的关系进行了研究。结果表明:(1)南沙群岛海域共出现长腹剑水蚤15种(包括1个未定种),其中线长腹水蚤(Oithona linearis)、粗长腹剑水蚤(O.robusta)和长刺长腹剑水蚤(O.longispina)为南沙群岛海域首次记录,长腹剑水蚤属种类组成沿水深变化的现象较为明显,该属虽然在南沙群岛海域广泛分布,但单种的出现频率较低,其种类出现的季节性较强;(2)南沙群岛海域的长腹剑水蚤可以划分为常见广布组群Ⅰ和低频组群Ⅱ;(3)长腹剑水蚤年均栖息密度为30ind./m3,季节变化呈单峰型,秋季最高、冬季最低,高数量区主要出现在受南沙西部沿岸流和东部沿岸流影响的近岸海域;(4)在0—750m水深范围内,长腹剑水蚤数量沿水深梯度的垂直变化明显,呈现由浅至深数量明显减少的趋势,长腹剑水蚤主要集中分布在0—75m的水层中,75—500m数量急剧降低至最低后,在500—750m基本维持不变;(5)GAM分析表明,纬度、海水温度、盐度和叶绿素a浓度对长腹剑水蚤数量均有影响,以纬度的影响最为显著,其中长腹剑水蚤属最适温度为28.6—29.2℃,最适盐度为32.6—33.2,最适叶绿素a浓度为5—10μg/L;(6)长腹剑水蚤属更适栖息于低温和低盐,且受沿岸流影响的环境中。因此,其数量分布的季节变化与季风驱动的沿岸流势力强弱密切相关。  相似文献   

雅鲁藏布江流域极端降水模拟及预估     

高佳佳  杜军 《冰川冻土》2021,43(2):580-588

为了解雅鲁藏布江流域汛期极端降水的变化规律,推算一定重现期的极端降水量分位数,通过百分位法、Hill图法、年交叉率法选取阈值,借助广义帕累托分布函数（GPD）对流域极端降水频率进行了分析。结果表明：99百分位时的阈值为流域内各站点的最佳阈值,且各站点超阈值序列通过了M-K的平稳性检验,无明显突变。拟合效果通过K-S检验,各站点拟合的极端降水理论频数和实测频数基本相符。尺度参数的大值区位于流域下游,表明该地区的极值波动大;形状参数正值区位于流域中上游地区,说明发生破纪录降水事件的概率较大,拟合结果与实际观测一致。从5年一遇和10年一遇的极端降水值来看,雅江流域除拉孜站外,其他地区降水极值均超过30 mm,日喀则地区的降水极值达50 mm;各地区20年一遇和30年一遇的降水极值增长的非常缓慢。通过与实际极端降水值对比分析得出,GPD拟合计算出的重现期水平基本符合实际,即具有一定的合理性。  相似文献