共查询到18条相似文献，搜索用时 140 毫秒
1.
2.
《岩土力学》2015,(Z1):577-584
天然沉积砂土力学特性受各向异性及结构性影响明显,实际工程中不能忽视。为探究二者的影响,首先在二维离散元程序NS2D中采用椭圆颗粒模拟了重力场中颗粒长轴主方向为水平的各向异性净砂样,随后基于结构性砂土胶结厚度分布规律及室内试验提出了一个新的微观胶结接触模型并将其引入各向异性净砂样以模拟天然各向异性结构性砂土,最后对该离散元试样进行了双轴试验模拟,将模拟结果与室内试验结果对比以验证该模型的适用性,并对其微观力学特性变化进行研究。分析结果表明,随着剪切进行,各向异性结构性砂土呈明显应变软化及剪胀现象;胶结接触逐渐减少,且主方向始终为竖向方向;胶结破坏速率及胶结破坏率变化情况与宏观力学响应较一致,且胶结物多为拉剪破坏;土颗粒排列主方向始终为水平向,且水平向排列颗粒所占比例略微增大。 相似文献
3.
人工胶结砂土力学特性的离散元模拟 总被引：1，自引：1，他引：0
采用离散单元法(DEM)对胶结砂土力学特性进行模拟。将基于室内试验测得的理想胶结颗粒接触力学响应引入到开发的二维离散元程序(NS2D)中,模拟胶结砂土颗粒间的胶结作用。对不同胶结强度和围压的胶结砂土进行平面应变双轴压缩试验模拟,并将模拟结果与Wang和Leung[1]提供的人工胶结砂土的试验结果进行比较。最后对数值模拟中胶结试样的微观力学响应(接触力链、胶结点破坏率和位移场)进行分析。结果表明,离散元数值模拟能够有效地反映胶结砂土的主要力学特性,相比同一初始孔隙比的无胶结松散砂土,胶结砂土将具有更高的强度,应力-应变关系呈应变软化,体变为先剪缩后剪胀,且两者的差异随胶结强度的增大和围压的减小而越趋显著。此外,胶结砂土宏观力学响应(应力-应变关系和剪胀性)与其微观力学响应密切相关。 相似文献
4.
采用离散元法（DEM）研究胶结岩土材料在不同加载条件下的结构破损规律。首先,基于微观力学理论,考虑胶结岩土材料颗粒间胶结特性,给出表征结构性损伤的破损参数式。该式具有微观物理意义,但不能直接用于建立宏观本构模型。其次,采用二维离散元源程序NS2D模拟等向压缩、等应力比压缩以及双轴压缩试验,分析破损参数在不同加载条件下随宏观力学变量（体积应变和剪应变）的演变规律。最后,根据模拟结果提出破损参数数学表达式,其为大主应变的函数。研究结果表明：胶结强度、应力比以及围压均在一定程度上影响了数值试样的结构破损规律。在等向压缩和等应力比压缩试验中,容易用函数式描述数值试样破损参数随体积应变或偏应变的演变规律;而在双轴压缩条件下,由于数值试样有剪胀特性,破损参数随体积应变的演变规律则不易描述。建议的破损参数数学表达式能够较好地描述数值试样在不同加载条件下结构破损规律。 相似文献
5.
基于微观胶结厚度模型的深海能源土宏观力学特性离散元分析 总被引：1，自引：0，他引：1
首先,引入蒋明镜等提出的考虑水合物胶结厚度的深海能源土粒间微观胶结模型,用以反映能源土颗粒之间水合物微观胶结接触力学特性;其次,采用C++语言将模型程序化,并将其引入离散单元法中;然后,对选定的水合物饱和度经过实际二维离散元模拟调算,得出相应的水合物胶结尺寸,以修正水合物临界胶结厚度、最小胶结厚度及胶结宽度,进而确定水合物微观胶结参数;最后,根据所确定的胶结参数,针对不同水合物饱和度试样进行能源土宏观力学特性离散元双轴试验模拟,并从应力-应变、体变、剪胀角等方面与Masui等所进行的能源土室内三轴试验进行对比分析。结果表明：采用考虑粒间胶结厚度的水合物微观胶结模型,能够定性反映深海能源土的宏观力学特性,能源土试样的峰值强度、黏聚力和剪胀角均随水合物饱和度的增加而增加,但水合物饱和度对内摩擦角的影响规律不明朗;能源土试样的峰值强度、残余强度及体积剪缩量随着有效围压的增大而增大;剪胀角随有效围压的增大而减小。 相似文献
6.
接触模型的宏?细观参数标定是成功使用离散元方法的关键。在离散元的接触模型中线性接触模型与抗转线性接触模型均可用于模拟砂性土的力学行为,其中抗转线性接触模型在模拟密砂的剪胀性方面具备优势。采用抗转线性接触模型对室内密实砂土三轴试验进行了离散元模拟,验证了抗转线性接触模型的可靠性;进而系统分析了颗粒间摩擦系数、刚度比和抗转动系数等细观参数与砂土峰值内摩擦角、残余内摩擦角、峰值剪胀角等宏观参数的相关关系并进行了验证;揭示了偏应力作用下,细观参数对密实砂土试样内部剪切带宽度与倾角变化的影响规律,提出了考虑剪胀角的剪切带倾角经验公式。通过研究建立了抗转线性接触模型宏?细观参数的量化关系并给出了标定参数的具体流程图,提出了快速标定宏观参数的方法并应用实例进行了验证,为采用抗转线性接触模型精准模拟密实砂土的力学特性提供依据。 相似文献
7.
软土预压工程中,初始和诱发各向异性对软土力学性质的影响十分显著,而现有研究缺乏对初始和诱发各向异性的统一研究方法。采用离散单元法,以颗粒长宽比作为定量评价指标,构建真实形态的颗粒模型,生成5组不同沉积角的初始各向异性试样,并进行竖直和水平两方向加载的双轴模拟实验,研究了初始各向异性和诱发各向异性对软土力学特性影响;在细观层面,以颗粒为对象研究了颗粒接触形式和转动角度的变化规律,以接触为对象研究了配位数和接触法向各向异性的发展趋势,在此基础上探究抗剪强度指标与各向异性关系。结果表明:初始和诱发各向异性共同影响试样力学性质,当加载方向和软土沉积方向垂直时,土体有最大的峰值强度。颗粒接触形式中面面接触的比例随加载的进行逐渐增大,并影响着试样初始模量和抗剪强度,配位数和接触法向各向异性受颗粒接触形式的影响有不同的演化规律,并在加载后期趋于稳定;同时,初始各向异性试样相较各向同性试样有更大的黏聚力,诱发各向异性主要影响试样内摩擦角,进而影响试样抗剪强度。 相似文献
8.
《岩土力学》2015,(Z2):639-647
天然气水合物以胶结及孔隙填充等形式存在于深海能源土中,开采时因其分解会劣化地层力学特性进而引发海底事故,使得人们对能源土开采过程进行中力学特性的变化愈发重视。在前期室内试验的基础上,将一个温度-水压-力学二维微观胶结模型引入离散元商业软件PFC2D中,通过对排气、排水性较好的土体进行升温及降压法开采进行数值模拟,并将模拟结果与相同条件下的室内试验结果对比,验证了该胶结模型的适用性。进一步分析了颗粒接触分布与颗粒平均纯转动率(averaged pure rotation rate,APR)在水合物分解时的变化情况。升温分解过程中随温度升高,颗粒总接触分布各向异性程度增大;胶结接触逐渐减少并始终保持主方向为水平方向,无胶结接触增多并始终保持主方向为竖直方向;APR值逐渐增大且正负值分布逐渐趋于集中。降压分解过程中随反(水)压降低,颗粒总接触由各向同性分布逐渐发展为主方向为竖直方向的各向异性,APR值较小且分布均匀;恢复反压后,试样进一步破坏,颗粒总接触各向异性更加明显,APR值增大且正负值呈集中分布。 相似文献
9.
10.
借助离散单元法（DEM）建立海底管道上浮模型,模拟分析了不同胶结强度下胶结试样的力学特性;通过微观力学响应（局部变形、胶结点破坏分布和位移场）对海底管道上浮破坏机制进行分析;对比分析了海床砂土胶结强度对管道上浮土体破坏的影响（上浮抗力和地表隆起量）。结果表明：伴随土体变形和胶结破坏,上覆土体的渐进破坏主要集中在上覆土体区域内;土表隆起量随着管道上浮位移的增长而增长;胶结强度较高时,管道周围土体无回填行为,上覆土体有拉裂带形成,仅管道附近土体产生胶结破坏,峰值上浮抗力随海床砂土胶结强度的增加有显著提高 相似文献
11.
A discrete element modelling of bonded granulates and investigation on the bond effect on their behaviour are very important to geomechanics. This paper presents a two‐dimensional (2‐D) discrete element theory for bonded granulates with bond rolling resistance and provides a numerical investigation into the effect of bond rolling resistance on the yielding of bonded granulates. The model consists of mechanical contact models and equations governing the motion of bonded particles. The key point of the theory is that the assumption in the original bond contact model previously proposed by the authors (55th CSCE‐ASCE Conference, Hamilton, Ont., Canada, 2002; 313–320; J. Eng. Mech. (ASCE) 2005; 131 (11):1209–1213) that bonded particles are in contact at discrete points, is here replaced by a more reliable assumption that bonded particles are in contact over a width. By making the idealization that the bond contact width is continuously distributed with the normal/tangential basic elements (BE) (each BE is composed of spring, dashpot, bond, slider or divider), we establish a bond rolling contact model together with bond normal/tangential contact models, and also relate the governing equations to local equilibrium. Only one physical parameter β needs to be introduced in the theory in comparison to the original bond discrete element model. The model has been implemented into a 2‐D distinct element method code, NS2D. Using the NS2D, a total of 86 1‐D, constant stress ratio, and biaxial compressions tests have been carried out on the bonded granular samples of different densities, bonding strengths and rolling resistances. The numerical results show that: (i) the new theory predicts a larger internal friction angle, a larger yielding stress, more brittle behaviour and larger final broken contact ratio than the original bond model; (ii) the yielding stress increases nonlinearly with the increasing value of β, and (iii) the first‐yield curve (initiation of bond breakage), which define a zone of none bond breakage and which shape and size are affected by the material density, is amplified by the bond rolling resistance in analogous to that predicted by the original bond model. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
12.
Parametric study of the smooth‐joint contact model on the mechanical behavior of jointed rock 
下载免费PDF全文
下载免费PDF全文 The smooth‐joint contact model based on distinct element method has been widely used to represent discontinuity in the simulation of fractured rock mass, but there is rare efficient guidance for the selection of proper parameters of smooth‐joint contact model, which is the basement for using this model properly. In this paper, the effect of smooth joint parameters on the macroscopic properties and failure mechanism of jointed rock under triaxial compression test is investigated. The numerical results reveal that the friction coefficient of smooth joint plays a dominant role in controlling mechanical behaviors. The stiffness of smooth joint has a relative small influence on the mechanical behaviors. Poisson ratio decreases with the reduction of normal stiffness but increases with the reduction of shear stiffness. The reduction of smooth joint strength, which is determined by normal strength, cohesion, and friction angle of smooth joint, contributes to the breakage of bonded smooth joint and ultimately decreases the strength of the specimen. We proposed a detailed calibration process for smooth‐joint contact model according to the relationship between smooth‐joint parameters and mechanical properties. By following this process, the numerical results are validated against corresponding experimental results and good agreement between them can be found in stress‐strain curves and failure modes of different joint orientations. Further analyses from the microperspective are performed by looking at transmission of contact force, the nature and distribution of microcracks, and the particle displacement to show the failure process and failure modes. 相似文献
13.
14.
A discrete element model is proposed to examine rock strength and failure. The model is implemented by UDEC which is developed for this purpose. The material is represented as a collection of irregular-sized deformable particles interacting at their cohesive boundaries. The interface between two adjacent particles is viewed as a flexible contact whose stress–displacement law is assumed to control the material fracture and fragmentation process. To reproduce rock anisotropy, an innovative orthotropic cohesive law is developed for contact which allows the interfacial shear and tensile behaviours to be different from each other. The model is applied to a crystallized igneous rock and the individual and interactional effects of the microstructural parameters on the material compressive and tensile failure response are examined. A new methodical calibration process is also established. It is shown that the model successfully reproduces the rock mechanical behaviour quantitatively and qualitatively. Ultimately, the model is used to understand how and under what circumstances micro-tensile and micro-shear cracking mechanisms control the material failure at different loading paths. 相似文献
15.
A simple three‐dimensional distinct element modeling of the mechanical behavior of bonded sands 下载免费PDF全文
下载免费PDF全文 This paper presents a simple three‐dimensional (3D) Distinct Element Method (DEM) for numerical simulation of the mechanical behavior of bonded sands. First, a series of micro‐mechanical tests on a pair of aluminum rods glued together by cement with different bond sizes were performed to obtain the contact mechanical responses of ideally bonded granular material. Second, a 3D bond contact model, which takes into account the influences of bond sizes, was established by extending the obtained 2D experimental results to 3D case. Then, a DEM incorporating the new contact model was employed to perform a set of drained triaxial compression tests on the DEM bonded specimens with different cement contents under different confining pressures. Finally, the mechanical behavior of the bonded specimens was compared with the available experimental results. The results show that the DEM incorporating the simple 3D bond contact model is able to capture the main mechanical behavior of bonded sands. The bonded specimen with higher cement content under lower confining pressure exhibits more pronounced strain softening and shear dilatancy. The peak and residual strengths, the apparent cohesion and peak/residual friction angles, and the position and slope of the critical state line increase with increase in cement content. Microscopically, bond breakage starts when the system starts to dilate and the maximum rate of bond breakage coincides with the maximum rate of dilation. Bond breakage is primarily due to tension‐shear failure and the percentage of such failures is independent of both confining pressure and cement content. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
16.
A coupled 3‐dimensional bonded discrete element and lattice Boltzmann method for fluid‐solid coupling in cohesive geomaterials 下载免费PDF全文
下载免费PDF全文 This paper presents a 3D bonded discrete element and lattice Boltzmann method for resolving the fluid‐solid interaction involving complicated fluid‐particle coupling in geomaterials. In the coupled technique, the solid material is treated as an assembly of bonded and/or granular particles. A bond model accounting for strain softening in normal contact is incorporated into the discrete element method to simulate the mechanical behaviour of geomaterials, whilst the fluid flow is solved by the lattice Boltzmann method based on kinetic theory and statistical mechanics. To provide a bridge between theory and application, a 3D algorithm of immersed moving boundary scheme was proposed for resolving fluid‐particle interaction. To demonstrate the applicability and accuracy of this coupled method, a benchmark called quicksand, in which particles become fluidised under the driving of upward fluid flow, is first carried out. The critical hydraulic gradient obtained from the numerical results matches the theoretical value. Then, numerical investigation of the performance of granular filters generated according to the well‐acknowledged design criteria is given. It is found that the proposed 3D technique is promising, and the instantaneous migration of the protected soils can be readily observed. Numerical results prove that the filters which comply with the design criteria can effectively alleviate or eliminate the appearance of particle erosion in dams. 相似文献
17.
将基于圆化多边形离散单元法与有限元方法结合,提出一种可变形圆化多边形离散单元法。此法对块体离散元进行圆化处理,可较好地表征不规则块体外形,又保留了颗粒离散元计算高效的优势。在求解接触力时,消除了角点处法向奇异等问题,同时增强计算的稳定性和简化接触判断。同时对切向接触力计算模型进行修正,使得接触力计算效率得到提高。此法突破了圆化多边形刚体假设的限制,可以精确计算任意形状不规则离散单元之间的相互作用,对单元的运动和变形进行模拟。通过超静定梁冲击试验、不规则块体单轴压缩试验和料斗流动“卡阻”试验3个数值模拟算例,论证此法可以有效地捕捉单元的碰撞、分离和变形等空间运动和自身特性以及其细观力学表征。 相似文献
18.
Huillca Yoshiro Silva Matías Ovalle Carlos Quezada Juan Carlos Carrasco Sergio Villavicencio Gabriel Eduardo 《Acta Geotechnica》2021,16(3):699-709
Empirical evidence has shown that particle breakage affects the mechanical behaviour of granular materials. The source of this mechanism takes place at the particle scale, and the main consequence on the macromechanical behaviour is increasing compressibility. Due to the inverse correlation between particle size and particle crushing strength, coarse rockfill materials are particularly vulnerable to mechanical degradation due to particle breakage. However, such coarse materials do not fit in standard laboratory devices, and the alternative of large sample testing is usually unavailable or too expensive. Alternatively, recent works have proposed multi-scale approaches using the discrete element method (DEM) to carry out numerical testing of coarse crushable materials, although few studies have focused on size effects. This article presents the application of a DEM bonded-cell model to study particle size-strength correlation on angular rock aggregates. Each particle is modelled by a cluster of perfectly rigid polyhedral cells with Mohr–Coulomb contact law. Constant cell density within particles implies that the presence of potential fragmentation planes increases with size. Therefore, particle strength decreases with size. A comprehensive sensitivity analysis was carried out through 1477 particle crushing simulations in a given particle size. Based on published experimental data on calcareous rock aggregates, part of the simulations were used for calibration, and 97 additional simulations of a coarser size fraction were performed for validation. The results show a good agreement with the empirical data in terms of size effect and data scatter through Weibull statistics.
相似文献