共查询到20条相似文献，搜索用时 31 毫秒
1.
The distinct lattice spring model (DLSM) is a newly developed numerical tool for modeling rock dynamics problems, i.e. dynamic failure and wave propagation. In this paper, parallelization of DLSM is presented. With the development of parallel computing technologies in both hardware and software, parallelization of a code is becoming easier than before. There are many available choices now. In this paper, Open Multi‐Processing (OpenMP) with multicore personal computer (PC) and message passing interface (MPI) with cluster are selected as the environments to parallelize DLSM. Performances of these parallel DLSM codes are tested on different computers. It is found that the parallel DLSM code with OpenMP can reach a maximum speed‐up of 4.68× on a quad‐core PC. The parallel DLSM code with MPI can achieve a speed‐up of 40.886× when 256 CPUs are used on a cluster. At the end of this paper, a high‐resolution model with four million particles, which is too big to handle by the serial code, is simulated by using the parallel DLSM code on a cluster. It is concluded that the parallelization of DLSM is successful. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
2.
A full rate-dependent cohesive law is implemented in the distinct lattice spring method (DLSM) to investigate the dynamic fracturing behavior of brittle materials. Both the spring ultimate deformation and spring strength are dependent on the spring deformation rate. From the simulation results, it is found that the dynamic crack propagation velocity can be well predicted by the DLSM through the implemented full rate-dependent cohesive law. Furthermore, a numerical investigation on dynamic branching is also conducted by using the DLSM code. 相似文献
3.
A validation study of the distinct lattice spring model (DLSM) for wave propagation problems is performed. DLSM is a microstructure-based numerical model, which is meshless and has advantages in modelling dynamic problems where stress wave propagation is important. To verify the applicability of DLSM to modelling wave propagation through a discontinuous medium, the virtual wave source (VWS) method is used to obtain analytical solutions for wave propagation across a jointed rock mass. Numerical modelling results of the commercial code UDEC are selected as the reference. The effects of particle size and lattice rotation angle on wave propagation are first studied. Then, the results of wave transmission across a single joint with a different joint stiffness and across multiple parallel joints with different joint spacings are derived with DLSM, UDEC and VWS. These results are in good agreement with each other. Therefore, the capability of DLSM to model P-wave propagation across jointed rock mass is verified, which provides confidence for the further application of DLSM to modelling more complex problems. 相似文献
4.
A 3D distinct lattice spring model (DLSM) is proposed where matter is discretized into individual particles linked by springs. The presented model is different from the conventional lattice spring models where a shear spring is introduced to model the multibody force by evaluating the spring deformation from the local strain rather than the particle displacement. By doing this, the proposed model can represent the diversity of Poisson's ratio without violating the rotational invariance. The local strain of the spring is calculated through a least square method which makes the model possessing meshless properties. Because of this and explicitly representing the microstructure, DLSM is able to model dynamic fracturing problems and can be used to study the microstructure influences. The material parameters inputted in the model is the conventional material parameters, e.g. the elastic modules and the Poisson's ratio. Relationships between microscopic spring parameters and macroscopic material constants are derived based on the Cauchy–Born rules and the hyperelastic theory. Numerical examples are presented to show the abilities and properties of DLSM in modeling elastic and dynamic failure problems. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
5.
This study reports a Graphics Processing Unit (GPU)-based parallelization of the Distinct Lattice Spring Model (DLSM) for geomechanics simulation. The DLSM is a newly developed numerical model for rock dynamics problems, i.e., dynamic failure and wave propagation. Despite its applicability, one of the drawbacks of this model is the high computational load for practical simulations. To tackle this problem, a GPU with a Compute Unified Device Architecture (CUDA) is adopted to parallelize the DLSM code. The performance of the GPU DLSM code is tested on two computers equipped with modern GPU cards. The results show that significant performance improvements are gained from GPU parallelization of the DLSM code (the maximum speed up achieved was 23×). 相似文献
6.
Numerical modelling of concrete cracking requires robust models able to describe opening and propagation of cracks. Structural concrete codes provide practical relations to describe crack openings. However, these empirical methods were developed for specific structures and cannot be used for general applications. Here, a continuous modelling approach based on damage mechanics is used to compute crack openings in a tie‐beam concrete structure. We propose a post‐processing method to extract crack openings from a continuum damage finite element computation. This method can be applied to all continuum damage/plasticity models. The tie‐beam concrete is characterized by a weak stress gradient; this aspect complicates predictions of crack positions and number. A stochastic method is used to take into account the spatial variability in concrete properties and create a spatially correlated random property field. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
7.
Formulation and algorithmic treatment of a rate‐dependent plastic–damage model modified to capture large tensile cracking in cyclic‐loaded concrete structures are presented in detail for a three‐dimensional implementation. The plastic–damage model proposed by Lee and Fenves in 1998 was founded based on isotropic damaged elasticity in combination with isotropic multi‐hardening plasticity to simulate cracking and crushing of concrete under cyclic or dynamic loadings. In order that the model can capture large crack opening displacements, which are inevitable in plain concrete structures, the excessive increase in plastic strain causing unrealistic results in cyclic behaviors is prevented when the tensile plastic–damage variable controlling the evolution of tensile damage is larger than a critical value. In such a condition, the crack opening/closing mechanism becomes similar to discrete cracking. The consistent tangent operator required to accelerate convergence rate is also formulated for the large cracking state including viscoplasticity. The validation and performance of the modified algorithm implemented in a special finite element program is exemplified through several single‐element tests as well as three structural applications. The last example examines the model in the seismic fracture analysis of Koyna dam as a benchmark problem and the resulting crack profile is compared with the available experiment. The numerical experimentations well demonstrate that the developed model whose modification is necessary to properly simulate the cyclic behavior of plain concrete subjected to large tensile strains is robust and reasonably accurate. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
8.
9.
Crack damage evolution of shale is crucial to the hydraulic fracturing treatment and engineering stability. Although many effects have been done on the macroscopic characteristics of shale, yet the microscopic failure mechanism is not well understood. A uniaxial compressive test on black shale was conducted under topographic monitoring using in situ X-ray micro-tomography (µCT). A series of high-resolution reconstruction images were obtained by carrying out CT scans at six key points throughout the test to obtain the internal structure of shale sample. In addition, the CT values for the purpose of crack damage evolution in shale were identified. Clear 2D/3D CT images, CT value analysis and image segmentation analysis reveal that the sample experiences compression, damage, cracking, crack propagation, and collapse stages. Crack geometry and distribution in the shale sample is visualized by rendered CT images, and a combined tension and shear failure mode is observed from the fracture rose diagram. This work suggests that formation and propagation of fractures are influenced by the stratified structure and weak cementation medium between layers. 相似文献
10.
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. 相似文献
11.
裂纹监测对岩石损伤演化的认识至关重要。为研究岩石裂纹扩展及损伤变形特性,开展了含不同倾角（0°～90°）预制裂隙的标准细黄砂岩样的单轴压缩试验。利用三维数字图像相关技术（3D-DIC）获取岩样三维空间坐标下的应变分布,并结合声发射从光学与声学的角度监测了裂纹的扩展演化。由此提出了一种裂纹主应变的计算方法,定量表征岩石劣化的损伤变量D值。最后,探讨了由声发射与损伤变量D值确定岩样特征强度的影响因素。结论如下：（1）裂纹主应变反映了岩样受荷过程中同源裂纹在时间上的变化速率与空间上的扩展趋势,能较好地表征岩石的开裂行为;（2）声发射适用于确定岩样的起裂应力,不适用于损伤应力的确定,损伤变量D值所确定的起裂应力滞后于声发射,但适用于损伤应力特征值的确定;（3）结合声发射与DIC技术确定的归一化起裂应力范围为0.63～0.94、归一化损伤应力的范围为0.83～0.99;（4）预有裂隙会影响岩石的材料力学性能。随着倾角的增加,岩石的起裂应力、损伤应力及峰值应力呈增长的趋势,由于难以形成局部应变场聚集,裂纹的萌生与起裂更加困难。结果表明,3D-DIC技术的利用可以提高对岩石开裂行为的理解,对岩石的损伤监测与判识更有重要的意义。 相似文献
12.
The present paper focuses on selected plasticity and damage‐plasticity models for describing the 3D material behavior of concrete. In particular, a plasticity model and a damage‐plasticity model are reviewed and evaluated. Based on the results of the evaluation, enhancements are proposed, aiming at improving the correspondence between predicted and observed material behavior and aiming at implementing a robust and efficient stress update algorithm in a finite element program for performing large‐scale 3D numerical simulations of concrete structures. The capabilities of the concrete models are demonstrated by 3D numerical simulations of benchmark tests with combined bending and torsional loading and combined compression and shear loading and by a large‐scale 3D finite element analysis of a model test of a concrete arch dam. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
13.
Wen Chen Christian La Borderie Olivier Maurel Gilles Pijaudier‐Cabot Franck Rey‐Bethbeder 《国际地质力学数值与分析法杂志》2014,38(5):457-474
The results reported in this paper deal with the simulation of damage in cohesive geomaterials such as rocks or concrete subjected to dynamic loads. The practical objective is to stimulate the production of tight gas reservoirs with a technique that is an alternative to hydraulic fracturing. The principle is that when subjected to dynamic loads, cohesive materials such as concrete, rocks or ceramics exhibit distributed micro‐cracking as opposed to localised cracking observed under static loads. Hence, a low permeability rock containing gas trapped into occluded pores can be fragmented with the help of dynamic loads, and gas can be extracted in a much more efficient way compared with hydraulic fracturing, where only large macro cracks are formed with very few connections between occluded pores. At the stage of laboratory development of this technique, compressive underwater shock waves have been used to increase the intrinsic permeability of concrete specimens. In a previous study, pressure waves generated by pulsed arc electrohydraulic discharges in water were used in order to induce micro‐cracking and an increase of average permeability of concrete hollow cylinders subjected to confinement stresses (equivalent to geostatic stresses). We discuss here a 3‐D anisotropic constitutive model aimed at describing the dynamic response of these specimens. It is based on rate‐dependent continuum damage constitutive relations. Crack closure effects and damage‐induced anisotropy are included in the model. The directional growth of damage is related to the directional growth of material intrinsic permeability. Numerical simulations of damage induced by shock waves show good agreement with the experiments for various confinement levels of the specimens. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
14.
为了探索煤层顶板中水平井向目标层穿层压裂的裂缝扩展规律,以华北石炭-二叠纪煤田为例,运用断裂力学、损伤力学以及流体力学等经典理论并结合现场实测资料,开展了压裂缝延伸距离与压裂时间时空演化规律的建模与验证。首先,基于原生裂缝特性、渗透特性以及压裂射孔段附加应力等因素,提出了顶板水平井垂向造缝的起裂压力计算公式;其次,在考虑裂缝性煤岩体损伤效应的基础上,引入Dougill损伤因子,将该计算模型拓展为延伸压力计算模型;最后,基于改进的经典PKN裂缝模型和压裂液滤失理论,建立了连续穿层工况下压裂缝延伸距离与压裂施工时间的函数关系。实践验证表明,根据理论模型合理调配时间参数,可以控制穿层裂缝的延伸距离。 相似文献
15.
在振动载荷作用下,岩石内部细观缺陷的发展将导致其宏观力学性能的劣化。为了提高超声波振动碎岩效率,采用数值模拟与实验研究相结合的方法,对超声波振动下岩石内部的损伤过程进行分析;采用有限元与统计技术相结合的方法,建立岩石二维模型,引入损伤增量因子的概念,分析岩石非均匀性对岩石裂纹扩展过程的影响。研究结果表明:1)岩石材料在超声波振动下的损伤过程分为萌生、扩展、贯通3个阶段,随着非均匀系数的增加,岩石原始裂纹的扩展效率逐渐加快。2)损伤因子存在一个临界值,约为0.005 4:当某一时刻岩石的损伤因子低于0.005 4时,损伤增量因子与非均匀系数的响应关系不明显;当损伤因子高于0.005 4时,损伤增量因子随着非均匀系数的增大而急剧增大,此时岩石非均匀性对损伤因子的影响不容忽视。 相似文献
16.
Three‐dimensional simulations of tensile cracks in geomaterials by coupling meshless and finite element method 
下载免费PDF全文
下载免费PDF全文 Failure in geotechnical engineering is often related to tension‐induced cracking in geomaterials. In this paper, a coupled meshless method and FEM is developed to analyze the problem of three‐dimensional cracking. The radial point interpolation method (RPIM) is used to model cracks in the smeared crack framework with an isotropic damage model. The identification of the meshless region is based on the stress state computed by FEM, and the adaptive coupling of RPIM and FEM is achieved by a direct algorithm. Mesh‐bias dependency, which poses difficulties in FEM‐based cracking simulations, is circumvented by a crack tracking algorithm. The performance of our scheme is demonstrated by two numerical examples, that is, the four‐point bending test on concrete beam and the surface cracks caused by tunnel excavation. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
17.
对预制交叉裂隙岩石试样进行裂纹扩展试验,研究了裂纹萌生、扩展、聚合过程,分析了主裂隙和轴向载荷夹角及主次裂隙夹角对裂纹起裂应力和聚合应力的影响,并用混合有限元-离散元程序,即图形处理器并行化的3D Y-HFDEM代码对试验进行了仿真计算,实现了岩石破坏从连续介质向非连续介质的过渡,对裂纹类及损伤破坏模式进行了识别,捕捉到了试验中难以发现的现象。研究表明：随主裂隙与轴向载荷夹角增加,裂纹聚合区的拉伸裂纹数量增加;裂纹起裂和聚合应力与主裂隙与轴向载荷夹角成正比;主次裂隙夹角增加,岩石的破坏模式由拉伸破坏转为剪切破坏,交叉裂隙加剧岩石破碎程度;主裂隙尖端萌生扩展的拉伸-剪切混合裂缝引起的破坏在岩石破坏中占主导地位,是导致岩体失去承载能力的主控裂纹;混合有限元-离散元仿真软件GPGPU并行化的3D Y-HFDEM IDE在岩石裂纹扩展研究中具有优势,可以捕捉实验室难以发现的损伤断裂类型,可以作为岩石裂纹扩展研究的有力工具。 相似文献
18.
A numerical study of the influence from pre‐existing cracks on granite rock fragmentation at percussive drilling 下载免费PDF全文
下载免费PDF全文 Mahdi Saadati Pascal Forquin Kenneth Weddfelt Per‐Lennart Larsson Francois Hild 《国际地质力学数值与分析法杂志》2015,39(5):558-570
The aim of this study is to investigate the effect of pre‐existing, or structural, cracks on dynamic fragmentation of granite. Because of the complex behavior of rock materials, a continuum approach is employed relying upon a plasticity model with yield surface locus as a quadratic function of the mean pressure in the principal stress space coupled with an anisotropic damage model. In particular, Bohus granite rock is investigated, and the material parameters are chosen based on previous experiments. The equation of motion is discretized using a finite element approach, and the explicit time integration method is employed. The pre‐existing cracks are introduced in the model by considering sets of elements with negligible tensile strength that leads to their immediate failure when loaded in tension even though they still carry compressive loads as crack closure occurs because of compressive stresses. Previously performed edge‐on impact tests are reconsidered here to validate the numerical model. Percussive drilling is simulated, and the influence of the presence of pre‐existing cracks is studied. The results from the analysis with different crack lengths and orientations are compared in terms of penetration stiffness and fracture pattern. It is shown that pre‐existing cracks in all investigated cases facilitate the drilling process. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
19.
Drying of deformable porous media results in their shrinkage, and it may cause cracking provided that shrinkage deformations are hindered by kinematic constraints. Herein, we focus on slow drying of an initially water‐saturated sample of a microheterogeneous poroelastic material damaged by parallel mesocracks. The cracking risk is analyzed by means of the thermodynamics‐based microporoelasticity model described in the companion paper (Part I), which is extended toward consideration of the hierarchical organization of cracked argillite. Drying of a material sample is studied in a framework where macrodisplacements in direction of the crack normal are blocked, while elsewise macrostress‐free boundary conditions prevail. The model implies that the opening/closure behavior of the cracks is governed by an effective pressure, in which the average crack (under)pressure, making the crack opening smaller, competes with the average micropore (under)pressure that makes the crack opening larger. The driving force for crack propagation is a quadratic function of this effective pressure. The model proposes that if drying shrinkage deformations are hindered by kinematic constraints, onset of cracking becomes possible once air entry into the cracks is observed. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
20.
Numerical modeling of dynamic rock fracture with a combined 3D continuum viscodamage‐embedded discontinuity model 下载免费PDF全文
下载免费PDF全文 This paper deals with numerical modeling of dynamic failure phenomena in rate‐sensitive quasi‐brittle materials, such as rocks, with initial microcrack populations. To this end, a continuum viscodamage‐embedded discontinuity model is developed and tested in full 3D setting. The model describes the pre‐peak nonlinear and rate‐sensitive hardening response of the material behavior, representing the fracture‐process zone creation, by a rate‐dependent continuum damage model. The post‐peak response, involving the macrocrack creation accompanied by exponential softening, is formulated by using an embedded displacement discontinuity model. The finite element implementation of this model relies upon the linear tetrahedral element, which seems appropriate for explicit dynamic analyses involving stress wave propagation. The problems of crack locking and spreading typical of embedded discontinuity models are addressed in this paper. A combination of two remedies, the inclusion of viscosity in the spirit of Wang's viscoplastic consistency approach and introduction of isotropic damaging into the embedded discontinuity model, is shown to be effective in the present explicit dynamics setting. The model performance is illustrated by several numerical simulations. In particular, the dynamic Brazilian disc test and the Kalthoff–Winkler experiment show that the present model provides realistic predictions with the correct failure modes and rate‐dependent tensile strengths of rock at different loading rates. The ability of initial embedded discontinuity populations to model the initial microcrack populations in rocks is also successfully tested. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献