共查询到20条相似文献，搜索用时 406 毫秒
1.
We present a multiscale investigation on the initiation and development of compaction bands in high-porosity sandstones based on an innovative hierarchical multiscale approach. This approach couples the finite element method and the discrete element method (DEM) to offer direct, rigorous linking of the microscopic origins and mechanisms with complex macroscopic phenomena observed in granular rocks such as strain localization and failure. To simulate compaction band in granular cementitious sandstone, we adopt a bonded contact model with normal and tangential interparticle cohesions in the DEM and propose a dual-porosity structure consisting of macro-pores and interstitial voids for the representative volume element to mimic the typical meso-structure of high-porosity sandstones. In the absence of particle crushing, our multiscale analyses identify debonding and pore collapses as two major contributors to the formation of compaction bands. The critical pressures predicted by our simulations, corresponding to surges of debonding and pore collapse events, agree well with the estimations from field data. The occurrence patterns of compaction band are found closely related to specimen heterogeneity, porosity and confining pressure. Other deformation band patterns, including shear-enhanced compaction bands and compactive shear bands, were also observed under relatively low confining pressure conditions with a rough threshold at \(0.55P^{*}\) (\(P^{*}\) is the critical pressure) on the failure envelop. Key microscopic characteristics attributable to the occurrence of these various deformation patterns, including fabric anisotropy, particle rotation, debonding and pore collapse, are examined. Shear-enhanced compaction bands and pure compaction bands bear many similarities in terms of these microscopic characteristics, whereas both differ substantially from compactive shear bands. 相似文献
2.
Modelling shear band is an important problem in analysing failure of earth structures in soil mechanics. Shear banding is the result of localization of deformation in soil masses. Most finite element schemes are unable to model discrete shear band formation and propagation due to the difficulties in modelling strain and displacement discontinuities. In this paper, a framework to generate shear band elements automatically and continuously is developed. The propagating shear band is modelled using discrete shear band elements by splitting the original finite element mesh. The location or orientation of the shear band is not predetermined in the original finite element mesh. Based on the elasto‐perfect plasticity with an associated flow rule, empirical bifurcation and location criteria are proposed which make band propagation as realistic as possible. Using the Mohr–Coulomb material model, various results from numerical simulations of biaxial tests and passive earth pressure problems have shown that the proposed framework is able to display actual patterns of shear banding in geomaterials. In the numerical examples, the occurrence of multiple shear bands in biaxial test and in the passive earth pressure problem is confirmed by field and laboratory observations. The effects of mesh density and mesh alignment on the shear band patterns and limit loads are also investigated. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
3.
为了建立能够表征煤体真实孔隙结构的数字模型,进而进行数值模拟研究,通过μCT225kVFCB高精度CT系统扫描得到了大柳塔煤矿长焰煤CT数据,可观测到的最小孔隙直径为1.94 μm。使用基于Matlab语言的三维重建程序结合逆向工程技术,提出了一种将煤体CT三维数据转换为CAD数字模型的方法。以Ansys建立的煤体有限元模型为例,对煤体瓦斯渗流进行了模拟,分析了煤体孔隙内的速度及压力分布规律,并计算了沿X、Y、Z方向的渗透系数。计算结果表明：在微观尺度下（< 100 μm）煤体渗透系数呈现各向异性,其受煤体结构的影响较明显。提出的基于CT三维重建结合逆向工程技术构建的煤体CAD数字模型不仅可进行有限元分析,同时也可以被EDEM等离散元分析软件所使用,拓宽了煤体CT三维数据的应用范围,丰富了煤体在微观尺度上的研究方法。 相似文献
4.
不同应力路径下超固结黏土试样变形局部化分析 总被引：2，自引：1，他引：1
基于改进伏斯列夫面超固结黏土三维本构模型,利用有限元软件ABAQUS材料子程序接口,采用回映应力更新算法,实现了该模型在有限元分析中的应用。通过该模型与比奥固结理论的耦合,对超固结比为8的超固结黏土在三轴压缩、三轴伸长及平面应变应力条件下的变形局部化问题,进行了水-土耦合弹塑性有限元分析。分析结果表明：剪切带带内、带外点经历不同应力路径;剪切带带外单元经历了体缩、剪胀及被吸水体缩过程,而剪切带带内单元一直保持剪胀趋势;剪切带的形成伴随着剪胀,剪切带内、外出现了负的孔压,且孔压的分布也具有局部化特性。关于剪切带带内、外的孔隙水压及体变变化趋势与剪切速率有关,而平面应变介于三轴压缩与三轴伸长之间,但平面应变较早出现剪切带。孔隙水的迁移速度影响剪切带带内单元的剪胀,进而影响剪切带的形成及发展;而围压和弱单元位置也对剪切带的形成也有影响。 相似文献
5.
Stress-induced evolution of anisotropic thermal conductivity of dry granular materials 总被引：1，自引：1，他引：0
Jinhyun Choo Young Jin Kim Jung Hwoon Lee Tae Sup Yun Jangguen Lee Young Seok Kim 《Acta Geotechnica》2013,8(1):91-106
Various factors, such as the volumetric fraction of constituents, mineralogy, and pore fluids, affect heat flow in granular materials. Although the stress applied on granular materials controls the formation of major pathways for heat flow, few studies have focused on a detailed investigation of its significance with regard to the thermal conductivity and anisotropy of the materials. This paper presents a numerical investigation of the stress-induced evolution of anisotropic thermal conductivity of dry granular materials with supplementary experimental results. Granular materials under a variety of stress conditions in element testing are analyzed by the three-dimensional discrete element method, and quantitative variations in their anisotropic effective thermal conductivity are calculated via the network model and conductivity tensor measurements. Results show that the directional development of contact area and fabric under anisotropic stress conditions leads to the evolution of anisotropy in thermal conductivity. The anisotropy induced in thermal conductivity by shear stress is higher than that induced by compressive stress because shear stress causes more significant changes in microstructural configurations and boundary conditions. The shear-stress-induced evolution of anisotropy between principal thermal conductivities depends on dilatancy as well as shearing mode, and the shear-driven discontinuity localizes the conductivity. Factors involved in the stress-induced evolution and their implications on the thermal conductivity characterization are discussed. 相似文献
6.
考虑强度各向异性的边坡稳定有限元分析 总被引：1，自引：0，他引：1
天然沉积的土层总是表现出一定程度的强度各向异性,但现有的边坡稳定有限元分析极少考虑各向异性的影响。对大型有限元软件ABAQUS进行二次开发,使其能够考虑土体黏聚力随大主应力方向的变化,动态更新增量迭代求解过程中边坡不同位置处的抗剪强度,进而提出具备安全系数自动搜索功能的各向异性边坡稳定分析方法。计算结果表明,均匀边坡的有限元解与极限分析上限解相差很小。如果采用土体固结方向的黏聚力并按各向同性评价缓坡的稳定性,可能严重高估安全系数,尤其是在黏聚力较高的情况下。与极限分析不同,建立的强度更新有限元模型能够分析成层边坡的稳定性。 相似文献
7.
A new method is proposed for the development of a class of elastoplastic thermomicromechanical constitutive laws for granular
materials. The method engenders physical transparency in the constitutive formulation of multiscale phenomena from the particle
to bulk. We demonstrate this approach for dense, cohesionless granular media under quasi-static loading conditions. The resulting
constitutive law—expressed solely in terms of particle scale properties—is the first of its kind. Micromechanical relations
for the internal variables, tied to nonaffine deformation, and their evolution laws, are derived from a structural mechanical
analysis of a particular mesoscopic event: confined, elastoplastic buckling of a force chain. It is shown that the constitutive
law can reproduce the defining behavior of strain-softening under dilatation in both the mesoscopic and macroscopic scales,
and reliably predict the formation and evolution of shear bands. The thickness and angle of the shear band, the distribution
of particle rotation and the evolution of the normal contact force anisotropy inside the band, are consistent with those observed in discrete element simulations and physical experiments. 相似文献
8.
The quantitative analysis of the pore characteristics of granular materials has been often challenging due to arbitrarily shaped geometry of pores despite its significant implications. In this study, we investigate the size distribution and orientation of pores in dilative and contractive assemblies in the direct shear test by performing 3D discrete element simulations in conjunction with image processing of pore geometry. We quantitatively define unit pores by the Delaunay Tessellation followed by pore mergence and fitting them with ellipsoids. It is observed that the evolution of pore size distribution depends on the dilatancy of assemblies. Results also show that the direction of principal stresses governs the orientations of pores during shearing, with respect to the size of pores. This study highlights that the dominant factors of the pore characteristics upon shearing are stress anisotropy and particle mobilization to make the internal structure stable. 相似文献
9.
This paper investigates two aspect of the direct simple shear (DSS) using three-dimensional finite element analysis. Firstly, the different total stress paths followed by DDS devices that impose constant cross-sectional area using a stiff external boundary, and those that use a constant total stress lateral boundary condition are explored. This is done by conducting finite element analysis of a single cubic element. It is shown that this element may be subjected to perfect simple shear using four different boundary condition types. Each boundary condition type results in the same effective stress path, but different total stress paths and excess pore pressures. The boundary condition types are related to DSS devices in use. The second aspect investigated is the fact that no DSS device can impose true simple shear conditions, as they are unable to generate the required complementary shear stress on the vertical boundaries. Full three-dimensional finite element analysis of a UWA/Berkeley type simple shear device, which has a constant total stress lateral boundary, is then presented. The results are compared to the ideal solution and effects of stress non-homogeneity on the friction angle and undrained shear strength, deduced from a standard interpretation, are established. 相似文献
10.
不同应力路径下剪切带的数值模拟 总被引：4，自引：1，他引：3
采用回映应力更新算法,编写了基于伏斯列夫面的超固结黏土本构关系模型子程序,嵌入非线性有限元软件ABAQUS。通过对单元试验进行三轴压缩、三轴伸长及平面应变等问题的模型预测,再现了超固结黏土在不同初始超固结比和应力路径时的变形和强度特性,从而验证了子程序的正确性。借助该本构模型,对三轴压缩、三轴伸长及平面应变应力路径下超固结黏土体变形局部化问题,进行了三维数值模拟。分析结果表明：超固结黏土在三轴压缩及伸长状态时,土体变形局部化在应力-应变关系软化时出现,而平面应变状态时,在应力-应变关系硬化阶段出现,其超固结黏土的剪胀特性在剪切带的形成过程中起重要作用。 相似文献
11.
This paper describes the mechanical behavior of artificially cemented sands with strong, intermediate, and weak bond strengths, using experimentation and 3D discrete element method (DEM) simulation. The focus is on the features of bond breakage and the associated influences on the stress–strain responses. Under triaxial shearing, the acoustic emission rate captured in the experiment and the bond breakage rate recorded in the simulations show resemblance to the stress–strain response, especially for strongly and intermediately cemented samples, where a strain softening response is observed. The simulations further reveal the shear band formation coincides with the development of bond breakage locations due to the local weakness caused by the bond breakages. Strain softening and volumetric dilation are observed inside the shear band, while the region outside the shear band undergoes elastic unloading. The weakly cemented sample exhibits a strain hardening response instead; bond breakages and the associated local weaknesses are always randomly formed such that no persistent shear band is observed. Note that in the DEM simulation, the flexible membrane boundary is established by a network of bonded membrane particles; the membrane particle network is further partitioned into finite triangular elements. The associated algorithm can accurately distribute the applied confining pressure onto the membrane particles and determine the sample volume. 相似文献
12.
为分析砂土在复杂应力条件下的剪切力学特性,采用商业离散元软件PFC3D对单粒组中密砂的空心扭剪试验进行了仿真模拟,分析了数值试样的应力-应变关系,研究了不同剪切方向下离散介质的强度、体积应变特性以及中主应力比对它们的影响,再现了力链在加载过程中的演化,并对剪切带的倾角做了深入分析。同时,从细观上看,以颗粒接触数和纯转动率变量为中心,观察了试样内部颗粒的运动状态,对比了不同剪切方向下剪切带内外颗粒接触数与纯转动位移的变化。最后,将数值试验结果与已有的室内试验结果进行了对比。此研究实现了复杂应力条件下空心扭剪试样的三维离散元模拟,加深了对空心扭剪试验过程和结果的理解和解释。 相似文献
13.
The local pore spaces in granular materials tend to be aligned parallel to the major principal stress direction upon particle mobilization. Manifestation of this response has been numerically validated in our previous studies with the aid of discrete element method modeling and image processing techniques during creep and shearing. We now extend the modeling of pore geometry, constructed with spherical particles, to assemblies of particle clumps. Two-dimensional simulations are performed for both loose and dense assemblies of spherical particles and particle clumps. Each particle packing is bound by rigid or flexible walls and subjected to biaxial compression and the particle mobilization effect on the evolution of pore orientation is explored. Randomly shaped pores surrounded by adjacent particles are geometrically quantified by Delaunay tessellation and fitted with ellipses. Results show that localization is apparent in dense assemblies, in particular for clumped particle packing, while loose assemblies exhibit diffusive failure. Small pores within well-defined shear bands tend to align either parallel to the direction of the shear band or perpendicular to the major principal stress. On the other hand, small pores within the blocks and large pores have a tendency to become elongate towards the major principal stress direction. This study reveals for the first time that pore orientation is dependent upon particle shape, pore size, and assembly conditions on the pore and global scales. 相似文献
14.
为研究岩石颗粒胶结方式对储层岩石弹性和渗流性质的影响,采用过程模拟法构建了三维数字岩心,在此基础上,分别利用有限元方法和格子玻尔兹曼方法研究了胶结物均匀生长、沿孔隙生长和沿喉道生长3种胶结方式对岩石弹性和渗流特性的影响规律.结果表明:岩石颗粒胶结方式会影响岩石刚性和孔隙连通性,引起岩石弹性模量和渗透率的变化.在相同孔隙度下,胶结物沿喉道生长形成的岩石抗压性最强,渗透率最小;沿孔隙生长形成的岩石抗压性最弱,渗透率最大.3种胶结方式下,岩石弹性模量随着胶结物含量增加而增大,变化率近似相等;岩石渗透率随着胶结物含量的增加而减小,岩石渗透性对颗粒胶结方式的变化更敏感. 相似文献
15.
软黏土层一维有限应变固结的超静孔压消散研究 总被引：1，自引：0，他引：1
根据土力学固结理论计算分析软黏土层固结过程的超静孔隙水压力值,确定软黏土体固结过程的强度增长,对排水固结法处理软土地基至关重要。软黏土层固结过程中土体变形较大时,有限应变固结理论和小应变固结理论计算分析软黏土固结所得结果差异较大。利用非线性有限元法及程序,通过对软黏土层固结工程算例的计算结果分析,研究了有限应变固结理论和小应变固结理论计算分析软黏土层一维固结超静孔压值消散的差异;探讨了软黏土体一维固结过程中,几何非线性、土体渗透性变化和压缩性变化对超静孔隙水压力消散的影响。研究结果表明,当土体的变形较大时,有限应变固结理论计算出的超静孔压要比小应变固结理论得到的值消散的更快。考虑土体固结过程中渗透性的变化时,超静孔压消散变慢;可用软黏土渗透性变化指数ck 反映渗透性变化对超静孔压消散的影响,渗透性变化指数ck值越小、超静孔压消散越慢。固结过程中软黏土压缩性的大小及变化也影响超静孔压的消散,可用软黏土的压缩指数cc反映固结过程中压缩性的大小及变化对超静孔压消散的影响,软黏土的压缩指数cc越小,固结过程软黏土层中的超静孔压消散越快。 相似文献
16.
饱和砂土局部变形带模拟的有限元数值实现 总被引：1，自引：0，他引：1
基于有限变形理论,推导了Newton-Raphson 迭代算法在k+1步增量表达的矩阵形式,实现了饱和砂土变形局部化的有限元数值计算,得到了饱和砂土发生局部化变形的准则。基于Galerkin 方法,得到了位移场和应力场的空间离散化矩阵方程;由土体局部变形带的连续性条件,引入第1切线算子,推导出了砂土等颗粒状媒介发生局部化变形的必要条件。基于此核心算法,编制了有限元计算程序,模拟了饱和砂土在不排水条件下平面压缩过程中剪切带的形成与发展;通过比较分析,研究了有限元网格粗细对于土体局部变形带的影响,结果表明,网格粗细的病态依赖只是微小的,它只与变形条带的宽度有关,对于土体所表现出来的其他力学特性没有影响。 相似文献
17.
Undrained deformation of dilative sand generates negative excess pore pressure. It enhances the strength, which is called dilative hardening. This increased suction is not permanent. The heterogeneity at the grain scale triggers localisations causing local volume changes. The negative hydraulic gradient drives fluid into dilating shear zones. It loosens the soil and diminishes the shear strength. It is essential to understand the mechanism behind this internal drainage and to capture it numerically. The purpose of this paper is to develop a macroscopic constitutive relationship for the undrained deformation of saturated dense sand in the presence of a locally fully or partially drained shear band. Separate constitutive relations are generated for the band and intact material. Both time and scale dependence during pore fluid diffusion in saturated sand are captured, eliminating the mesh dependency for finite element implementations. The model is applied to the Gauss points that satisfy the bifurcation criterion. The proposed method is calibrated to recreate the undrained macroscopic response bestowed by an extra-small mesh. The microscopic behaviours inside and outside shear band predicted by this model are qualitatively in good agreement with individual material point behaviours inside and outside the shear band in the extra-small mesh. Depending on the loading rate and the shear band thickness, the response inside the band can be fully or partially drained, which governs the ultimate global strength. The calibrated model is exploited to simulate an upscaled biaxial compression test with semipermeable boundaries. 相似文献
18.
T型触探仪近年来在国外的海洋软土地基原位测试中得到了广泛应用,但贯入速率、地基土强度的各向异性与渐进软化等因素对于贯入阻力的影响没有得到系统的研究。在大型有限元软件ABAQUS平台上进行二次开发,针对基于Tresca屈服准则的理想弹塑性模型进行了相应改进,使之可以比较合理地反映上述因素的影响,进而对软土中T型触探仪的贯入机制进行了比较系统的数值分析。计算结果表明,贯入速率、地基土强度的各向异性与强度软化效应对于T型触探仪的贯入阻力系数影响较大。通过与有关极限分析上限解的对比分析,在一定程度上验证了有限元分析结果的合理性。 相似文献
19.
A homogenization strategy for granular materials is presented and applied to a three-dimensional discrete element method (DEM), that uses superellipsoids as particles. Macroscopic quantities are derived from the microscopic quantities resulting from a DEM simulation by averaging over representative volume elements (RVEs). The implementation of an RVE is described in detail regarding the definition and discretization of the RVE boundary. The homogenization strategy is validated by DEM simulations of compression and shear tests of cohesionless granular assemblies. Finally, an elasto-plastic material model is fitted to the resulting stress–strain curves. 相似文献
20.
在波浪循环荷载作用下饱和软黏土地基出现孔隙水压力升高,并导致不排水强度弱化,严重影响防波堤的承载性能。考虑静偏应力影响,基于最大孔隙水压力发展模型和正常固结软黏土不排水强度公式,推导出软黏土不排水强度随循环荷载作用次数和应力水平变化的动态折减规律。结合软黏土不排水强度动态折减规律和M-C屈服准则,在有限元软件ABAQUS上实现软黏土不排水强度循环弱化分析的数值开发和动力运算过程。运用该动力有限元方法对天津港防波堤地基软黏土的动、静三轴试验进行数值模拟运算。结果表明,最大孔隙水压力发展曲线以及循环荷载作用后不排水强度的数值预测结果与动三轴试验结果吻合良好。另外,动力有限元方法（DFEM）能够表示土体强度在循环荷载作用下的具体弱化过程。 相似文献