Application of extended distinct element method with lattice model to collapse analysis of RC bridges     

Limin Sun  Cen Zhou  Dong Qin  Lichu Fan 《地震工程与结构动力学》2003,32(8):1217-1236

This paper proposes a simple lattice model for collapse analysis of RC bridges subjected to earthquakes by using the extended distinct element method (EDEM). In the model, a concrete element consists of lumped masses connected to one another by springs, and a reinforcement bar is represented by a discrete model or an integrated model. The proposed lattice model is simple but its parameters are reasonably defined. It has fewer element nodes and connecting springs, which will be of benefit by shortening the CPU time. The processes to determine the initial stiffness of concrete and steel springs, the parameters of the constitutive model and the fracture criteria for springs are described. A re‐contact spring model is also proposed to simulate the re‐contact of the concrete after fracture of springs; and a general grid searching method is used to decrease the CPU time for judging re‐contact after fracture. The lattice model is assessed by numerical simulations and experiments. As an application, a damaged single‐column pier subjected to the Kobe Earthquake in 1995 is analysed by EDEM with the proposed model. The simulation results indicate that the proposed model predicts well qualitatively the collapse process of RC bridges. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

孟湾槽在百色市异常降雪天气中的作用     

邓海光  唐毓勇 《云南地理环境研究》2004,(Z1)

通过对1983年12月28日、2002年12月26日百色市两次异常降雪天气过程的对比分析,发现冬季孟加拉湾槽的发展东移在其中起到了关键性的作用。由于孟加拉湾槽前暖湿气流强盛,这种强盛的西南暖湿气流在低空冷空气堆(冷锋)上发生强烈的系统性上升运动,从而导致百色市及其邻近地区出现异常的雨雪天气。  相似文献   

Polygonal grain‐based distinct element modeling for mechanical behavior of brittle rock          下载免费PDF全文

J.‐W. Park  C. Park  J.‐W. Song  E.‐S. Park  J.‐J. Song 《国际地质力学数值与分析法杂志》2017,41(6):880-898

The microstructure of rock was numerically reproduced by a polygonal grain‐based model, and its mechanical behavior was examined by performing the uniaxial compression test and Brazilian tests via the Universal Distinct Element Code. The numerical results of the model demonstrated good agreement with the experimental results obtained with rock specimens in terms of the stress–strain behavior, strength characteristics, and brittle fracture phenomenon. An encouraging result is that the grain‐based model‐Universal Distinct Element Code model can reproduce a low ratio of tensile to compressive strength of 1/20 to 1/10 without the need for an additional process. This finding is ascribed to the fact that the geometrical features of polygons can effectively capture the effects of angularity, finite rotation, and interlocking of grains that exist in reality. A numerical methodology to monitor the evolution of micro‐cracks was developed, which enabled us to examine the progressive process of the failure and distinguish the contribution of tensile cracking to the process from that of shear cracking. From the observations of the micro‐cracking process in reference to the stress–strain relation, crack initiation stress, and crack damage stress, it can be concluded that the failure process of the model closely resembles the microscopic observations of rock. We also carried out a parametric study to examine the relationships between the microscopic properties and the macroscopic behavior of the model. Depending on the micro‐properties, the model exhibited a variety of responses to the external load in terms of the strength and deformation characteristics, the evolution of micro‐cracks, and the post‐peak behavior. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

Implementation of configuration dependent stiffness proportional damping for the dynamics of rigid multi-block systems     

Yun?Byeong?Chae  Jae?Kwan?KimEmail author 《地震工程与工程振动(英文版)》2003,2(1):87-97

The distinct element method (DEM) has been used successfully for the dynamic analysis of rigid block systems. One of many difficulties associated with DEM is modeling of damping. In this paper, new procedures are proposed for the damping modeling and its numerical implementation in distinct element analysis of rigid multi-block systems. The stiffness proportional damping is constructed for the prescribed damping ratio, based on the non-zero fundamental frequency effective during the time interval while the boundary conditions remain essentially constant. At this time interval, the fundamental frequency can be estimated without complete eigenvalue analysis. The damping coefficients will vary while the damping ratio remains the same throughout the entire analysis. A new numerical procedure is developed to prevent unnecessary energy loss that can occur during the separation phases. These procedures were implemented in the development of the distinct element method for the dynamic analyses of piled multi-block systems. The analysis results for the single-block and two-block systems were in a good agreement with the analytic predictions. Applications to the seismic analyses of piled fourblock systems revealed that the new procedures can make a significant difference and may lead to much-improved results.  相似文献   

DEM analyses of crack propagation in flawed rock mass under uniaxial compressionEI北大核心CSCD     

Jiang Ming-jing  Zhang Ning  Shen Zhi-fu  Chen He 《岩土力学》2015,(11):3293-+

A contact model for rock is established and imbedded into a DEM software by summarizing the bond granule tests. DEM simulation of uniaxial compression test on the pre-cracked Lac du Bonnet granite is performed, and then stress distributions are further analyzed and compared with the theoretical results. Different fracture criteria are employed to predict the crack initiation angles that are compared with theoretical ones. The results show that the failure modes obtained from DEM simulation are similar to experimental results, and stress distributions in DEM simulation are qualitatively similar to theoretical values. When the angle of pre-crack is small, the lateral stresses are compressive and tensile. The compressive strains concentrate at two edges, resulting in the tensile strains in the up-and downward cracks. When the angle of the pre-crack is large enough, the stress concentration is unobvious, leading to a discrepancy between the DEM and theoretical results. The crack extension angle resulting from uniaxial compression measured from DEM tests are in good agreement with those acquired from experimental tests. These angles are consistent with theoretical predictions by the maximum circumferential stress criterion and the maximum energy release rate criterion.  相似文献   

Modelling failure mechanisms to explain rock slope change along the Isle of Purbeck coast,UK     

Robert J. Allison  Owen G. Kimber 《地球表面变化过程与地形》1998,23(8):731-750

Results are presented of distinct element computer modelling used to examine rates and mechanisms of change in rock slopes and cliffs, where material intact properties determine process and form but the most significant controls are the joint pattern and cross-joint properties. The modelling approach does not appear to have been used before in a geomorphological context and provides an alternative approach for examining cliff development. Field and laboratory data have been collected for the Portland Limestone outcrop of the Isle of Purbeck, central southern England. The Portland Limestome is a hard, shelly, crystalline sediment of the Upper Jurassic. It has a regular discontinuity pattern throughout the outcrop in Purbeck. While joint orientation remains relatively constant, bedding changes from horizontal to vertical, a consequence of the Purbeck Monocline. There are resulting implications for spatial variations in rock slope evolution. The modelling exercise enhances previous knowledge on rock failure mechanisms and slope development along the Purbeck coast and demonstrates its potential in research where landforms are developed in lithified, jointed rock masses. © 1998 John Wiley & Sons, Ltd.  相似文献   

不同重力场下静力触探试验离散元数值分析     

蒋明镜  王新新 《岩土力学》2013,34(3):863-873

采用离散元商业软件PFC2D对7种重力场下静力触探试验进行了对比分析试验,以研究高、低重力场下静力触探试验的异同点。数值试验结果表明：不同重力场下土体受静力触探贯入主要影响区域不同,低重力场下主要是上部土体受到影响,高重力场下主要是下部土体受到影响;归一化贯入阻力的最大值和稳定值、归一化球应力和偏应力的最大值均与重力加速度的倒数呈线性关系,均随重力加速度倒数的增加而增大;土体经历了明显的加载、卸载过程,重力加速度越小,加载、卸载现象越明显,同一重力加速度下上部土体加载、卸载程度大于下部土体。  相似文献   

基于三维流-固耦合模型的油井出砂细观机制研究   总被引：1，自引：0，他引：1  

刘先珊  许明 《岩土力学》2013,34(8):2363-2370

油井出砂机制的研究是提高油藏产能和石油开采成本减小的关键课题,而常规的宏观力学理论和方法不能全面反映油藏开采过程中油井出砂的发生和发展。鉴于砂岩储层的物理性质和射孔试验特征,从岩土力学的角度建立基于柱坐标系的三维颗粒流数值模型,与理论分析成果进行比较,以说明该细观数值模型可行性,有效地模拟出砂过程中的渗流及流-固耦合效应。在该基础上,综合考虑流体压力梯度力和拖曳力,基于PFC3D模型模拟流体不同运动时的砂岩性态。数值分析得到的模型宏观应力图形说明流体运动对砂岩力学特性的影响不可忽略,且在相同条件下,流量越大,砂岩的塑性区越大,形成砂岩破坏出砂的几率也越大。同时,不同工况的砂岩黏结分布和颗粒转动图形表明,相同条件下流量越大,颗粒间平行黏结破坏越多,颗粒转动越大,失去黏结约束的颗粒也越多,出砂量就越大,可见两种细观特征图形与宏观应力图形变化规律一致,该模型可用于油井出砂机制的研究,可为出砂量预测及出砂控制提供新的研究思路。  相似文献   

基于微观胶结厚度模型的深海能源土宏观力学特性离散元分析   总被引：1，自引：0，他引：1  

蒋明镜  贺洁  周雅萍 《岩土力学》2013,34(9):2672-2681

首先,引入蒋明镜等提出的考虑水合物胶结厚度的深海能源土粒间微观胶结模型,用以反映能源土颗粒之间水合物微观胶结接触力学特性;其次,采用C++语言将模型程序化,并将其引入离散单元法中;然后,对选定的水合物饱和度经过实际二维离散元模拟调算,得出相应的水合物胶结尺寸,以修正水合物临界胶结厚度、最小胶结厚度及胶结宽度,进而确定水合物微观胶结参数;最后,根据所确定的胶结参数,针对不同水合物饱和度试样进行能源土宏观力学特性离散元双轴试验模拟,并从应力-应变、体变、剪胀角等方面与Masui等所进行的能源土室内三轴试验进行对比分析。结果表明：采用考虑粒间胶结厚度的水合物微观胶结模型,能够定性反映深海能源土的宏观力学特性,能源土试样的峰值强度、黏聚力和剪胀角均随水合物饱和度的增加而增加,但水合物饱和度对内摩擦角的影响规律不明朗;能源土试样的峰值强度、残余强度及体积剪缩量随着有效围压的增大而增大;剪胀角随有效围压的增大而减小。  相似文献   

基于离散元的强度折减法分析岩质边坡稳定性   总被引：26，自引：0，他引：26  

雷远见  王水林 《岩土力学》2006,27(10):1693-1698

将通用离散元UDEC与强度折减法结合,对含多结构面的岩质边坡的稳定性进行了分析。通过对节理岩质边坡的UDEC模型中的可变形块体和节理单元的强度参数进行折减,使模型不能再达到平衡状态,此时的折减系数就是边坡的安全系数,另外,由对应的边坡块体的速度矢量可以确定滑动面和边坡的破坏形态。通过与传统的条分法的结果比较,表明基于UDEC的强度折减法是一种可靠、有效的方法,为复杂节理岩质边坡的滑动面确定与安全系数计算开辟了新的途径。  相似文献