离散元

摘要： 针对矿石破碎模型与实际矿石性质差异造成的破碎仿真不准确及磨矿工艺参数难以确定的问题,采用离散元法建立矿石颗粒黏结模型,重点研究基础颗粒的数量组合对颗粒黏结模型精度的影响。依据落重试验破碎的63~53 mm粒级矿石,构建了5组不同粒径颗粒的数量组合及对应的颗粒黏结模型。通过对颗粒黏结模型的破碎仿真及实际矿石的破碎试验,得到破碎后粒度分布结果,对比分析了仿真与试验的破碎后粒度分布差异。研究结果表明:当大颗粒体积不超过总体积的50%,且中间颗粒与大颗粒体积比为0.21~0.22时,破碎仿真粒度的分布与实际矿石破碎试验结果较为一致。在模型体积填充率不变的情况下,不同基础颗粒的数量组合对破碎后颗粒团簇粒度有较大影响,选择合适的基础颗粒数量组合能够提高颗粒黏结模型的精度。

摘要： 秸秆还田机还田弯刀的力学状态,直接影响还田质量及刀具可靠性。在含水率为30%,坚实度为700 kPa的黑土土壤中进行耕作参数优化试验,当耕深为200 mm、刀轴旋转速度为240 r/min、单位秸秆量为3.5 kg/m^(2)时,刀轴扭矩平均值最小。依据优化的耕作参数及测定的仿真接触参数和材料参数,采用离散元法建立还田弯刀—土壤和还田弯刀—土壤—秸秆的仿真互作模型,对比分析有、无秸秆两种状态下,弯刀所受扭矩的变化,并进行试验验证。仿真试验表明:有、无秸秆扭矩最大值分别为25.8 N·m、21.7 N·m,平均值分别为19.5 N·m、17.4 N·m,刀辊在有秸秆耕作时平均所消耗功率大于无秸秆状态的12.1%。验证试验表明:有、无秸秆状态扭矩最大值分别为32.3 N·m、26.7 N·m,平均值分别为23.7 N·m、21.3 N·m,有秸秆状态平均扭矩要高于无秸秆状态扭矩的11.3%左右。对比分析可知,有秸秆的试验值和仿真值,二者变化趋势基本相同,但试验值高于仿真值,且数据波动多、变化快。

摘要： 为准确表征沥青混合料的特性,文章提出混合料数字试件生成算法。首先通过随机切割算法,生成反映粗集料形状、级配、棱角性的数字集料;其次利用数字图像处理技术获取旋转压实试件内部空隙分布规律,并按照真实试件的分布规律,通过删除小粒径单元球使得数字试件中内部空隙与真实试件的一致;最后通过动态模量虚拟试验与室内试验对该数字试件生成算法进行对比验证。结果表明,虚拟试验仿真结果与室内试验结果较为相符,该数字试件生成算法切实可行,为揭示沥青混合料细观行为提供了新的手段。

摘要： 为提高绿豆精密排种过程离散元仿真模拟试验所用仿真参数的准确度,进一步优化排种部件,基于绿豆种子的本征参数,采用Hertz Mindlin with bonding粘结模型建立种子仿真模型,分别采用自由落体碰撞法、斜面滑动法、斜面滚动法对绿豆种子与接触材料(有机玻璃、Somos8000树脂)间仿真参数进行标定,结果表明:绿豆与有机玻璃碰撞恢复系数、静摩擦因数、滚动摩擦因数分别为0.445、0.458、0.036,与Somos8000树脂碰撞恢复系数、静摩擦因数、滚动摩擦因数分别为0.434、0.556、0.049。以种间接触参数为因素,以实测堆积角与仿真堆积角相对误差为指标,进行了最陡爬坡试验、三因素五水平旋转组合设计试验,以最小相对误差为优化目标,对试验数据寻优分析得到:绿豆种间碰撞恢复系数、静摩擦因数、滚动摩擦因数分别为0.3、0.23、0.03。对标定结果进行排种验证试验,结果表明:仿真试验漏吸率与台架试验漏播率最大相对误差为4.71%、重吸率与重播率最大相对误差为4.94%、单粒率与合格率最大相对误差为0.98%,证明标定结果可靠。该研究结果可为绿豆精密排种装置的设计与仿真优化提供理论参考。

摘要： 为寻求丸化包衣过程中紫花苜蓿种子进行离散元仿真模拟试验时最优接触参数组合,采用物理试验与仿真试验相结合的方法对离散元仿真参数进行标定。首先,采用物理试验测得紫花苜蓿种子的基本物性参数(密度、千粒质量、含水率、外形尺寸、泊松比和弹性模量)以及接触参数(碰撞恢复系数、静摩擦系数和滚动摩擦系数),以物理试验结果为根据,确定仿真试验参数范围,开展Plackett-Burman试验对仿真参数进行显著性筛选,得到对仿真休止角影响显著的3个参数分别为紫花苜蓿种子-紫花苜蓿种子碰撞恢复系数、紫花苜蓿种子-紫花苜蓿种子静摩擦系数以及紫花苜蓿种子-钢板静摩擦系数;进一步以物理试验休止角与仿真试验休止角相对误差值为指标,进行最陡爬坡试验,优化显著性参数的最优取值范围;之后开展Box-Behnken试验,并以物理试验休止角为目标值对二阶回归方程进行寻优,获得紫花苜蓿种子最优接触参数组合:紫花苜蓿种子-紫花苜蓿种子碰撞恢复系数为0.188,紫花苜蓿种子-紫花苜蓿种子静摩擦系数为0.684,紫花苜蓿种子-钢板静摩擦系数为0.371。最后以优化解进行仿真试验得到休止角为32.88°,与物理试验休止角相对误差仅为0.422%,且进行t检验得到P>0.05,说明仿真试验结果与物理试验无显著差异。结果表明,所得最优参数可用于紫花苜蓿种子离散元仿真,同时标定的参数可为紫花苜蓿种子的其他仿真试验提供参考。

摘要： 为研究行车荷载作用下沥青砂浆的迁移规律,借助离散元软件PFC2D建立了沥青混合料二维离散元模型,通过改变沥青砂浆Burgers模型的宏观参数,模拟了40、50和60°C温度下的沥青混合料黏弹状态,并利用时温等效原理,进行虚拟汉堡车辙试验,分析了沥青砂浆的迁移规律和温度对沥青砂浆迁移的影响。研究结果表明:在行车荷载作用下,随着试件深度的增加,沥青砂浆迁移量逐渐减小,但4 cm深度以下的沥青砂浆迁移量可忽略不计;温度越高,沥青砂浆Burgers模型宏观参数越小,沥青砂浆流动性越大,位移量越大。

摘要： 为解决黑水虻在畜禽粪便处理后续分离环节中应用筛分机械作业难、缺乏准确离散元仿真模型等问题,以黑水虻为研究对象,基于EDEM仿真软件模拟,选取系统中“Hertz-Mindlin with JKR”接触模型,以堆积角为评价指标,对黑水虻离散元仿真模型参数标定进行研究。首先,通过Plackett-Burman试验筛选对黑水虻堆积有显著影响的3个参数项(黑水虻剪切模量、黑水虻间静摩擦系数、黑水虻间滚动摩擦系数);然后,结合抽板法堆积物理试验,利用最陡爬坡试验确定显著性参数的最优区域;最后,进一步以堆积角为响应值,基于二次回归正交旋转组合试验得到堆积角与显著性参数的二阶回归模型,并以实际堆积角为目标,针对显著性参数进行寻优,确定EDEM仿真试验中黑水虻的最佳接触参数,得到最佳组合:黑水虻剪切模量8.67 MPa、黑水虻间静摩擦系数0.43、黑水虻间滚动摩擦系数0.32。运用最佳参数组合进行仿真分析,得到堆积角的均值为35.84°,与物理试验测得堆积角34.66°的相对误差为3.40%。研究结果表明,使用该方法对黑水虻离散元仿真参数的标定具有可行性,为基于EDEM软件设计与研究黑水虻筛分机械提供有效物料特性数据基础。

摘要： 为了提高桩土荷载分担比,一些建筑工程桩筏基础下的复合地基采用了加筋垫层,然而加筋垫层复合地基的变形协调与荷载传递的影响还缺乏深入研究。采用自制的二维多活动门试验装置与椭圆钢棒相似土填料,开展了无筋褥垫层随桩间土下沉过程的模型试验,获得了不同垫层厚度下无筋垫层的变形与应力协调特性。在模型试验基础上,采用颗粒离散元DEM数值模拟建立了无筋与加筋垫层数值模型,揭示了有、无加筋条件下的褥垫层工作特性区别,并探讨了褥垫层厚度、桩间净距以及加筋位置等对褥垫层工作特性的影响。结果表明:当加筋高度距桩顶超过100 mm时,桩土应力比曲线和褥垫层变形均接近未加筋的情况;加筋体设置高度小于100 mm时,下方颗粒垫层的变形协调会受到影响,格栅与下部垫层之间产生脱空现象,从而导致桩土应力比随着桩间土下沉量的增加而出现陡升;陡升段起点对应的桩间土下沉量随着加筋位置的提高而增加,工程中可能会导致复合地基桩体超过承载力而引起复合地基发生破坏。

摘要： 从数学领域的随机填充法入手,建立不同沙粒尺寸自组织轨道结构离散元分析模型,并借助“胞分法”保证沙粒粒径改变前后颗粒体系质量守恒。开展风沙区铁路轨道结构动测试验,标定接触模型的微观参数,验证离散元模型的正确性。在此基础上,利用不同沙粒粒径的数值试验,进行系统的道床宏细观力学特性分析。结果表明:沙粒内禀尺度愈小,道砟平均接触力愈小;道床内部颗粒之间的接触力概率密度随着接触力的增大而逐渐减小,沙粒粒径愈小,较大接触力出现的概率越低;峰值荷载作用下接触力玫瑰图表明,沙粒胞元愈多,道床中颗粒数目愈多,道砟颗粒接触力方向各向异性愈强,而道床中强接触力愈少;随着沙粒内禀尺度的逐渐减小,道床中高应力区域逐渐减小,低应力区域占比越来越高,道床应力均匀化程度增强。风沙道床研究中沙粒粒径合理尺度应介于0.51~1.02 mm之间。

摘要： 为研究以废弃钢渣为骨料制备的透水混凝土力学性能,选用两种不同粒径的钢渣骨料制备透水混凝土,通过轴压及劈裂拉伸加载对其基本力学性能进行试验;并结合PFC 3D软件针对钢渣透水混凝土的破坏形态及破坏过程进行了数值模拟。结果表明:骨料粒径大小对钢渣透水混凝土的基本力学性能存在较大影响,骨料粒径越大其静态强度越高;同时通过模拟发现,加载工况下大粒径骨料断裂的位置分布较为均匀,加载过程中所消耗的能量较大,而小粒径试件断裂位置集中在某一条裂缝发育路径上,最终造成试件沿着主导裂缝发生破坏,并且钢渣透水混凝土的破坏过程存在明显的阶段性特征,在宏观裂缝出现前断裂数稳步增长,以荷载达到峰值点为临界点,断裂数开始急剧增加。

摘要： 西南某露天铁矿开采形成了高陡边坡,边坡坡体中有不同构造的断层相互切割,节理较为发育,同时受软弱结构面影响比较明显.采用离散元法对高陡边坡稳定性进行计算,分析边坡在天然状态及开挖后的稳定性情况.模拟结果表明,由于受优势节理面控制,边坡坡脚处容易向临空面产生较大位移,可能发生浅表层局部失稳,整体稳定性较好.同时,也体现了离散元法在解决节理岩质边坡稳定性问题上具有较强适用性和优越性.

摘要： 用3DEC离散元数值模拟软件为深部硬岩矿床建立矿体模型,结合工程实际提出了在矿体中开挖巷道诱导矿岩致裂以实现矿体非爆连续开采.对诱导致裂巷道竖向、水平两种开挖形式进行了开挖模拟,模拟结果显示,竖向开挖比水平开挖利用地压诱导致裂矿岩的效果更好;掘进机上下分层开挖巷道时,先开挖下部的位移量较大,为7～9mm,塑性区主要集中在巷道两帮,主要以剪切、拉伸破坏为主,更容易利用地压破岩,研究成果成功地在开阳磷矿进行了深部高应力硬岩矿床非爆连续开采.

摘要： 砂土的密实状态对其力学性质影响很大,离散元软件在模拟砂土的力学性质方面已经有了许多进展,但是对砂土相对密度研究还比较少.为了研究离散元软件在模拟砂土相对密度方面的适用性,应用二维离散元软件,采用重力沉积法生成数值试样,进行了一系列数值试验,并与室内试验所得结果进行了对比,得出以下结论:对于二维数值试样,当试样与颗粒的尺寸比L/R＞70时,试样尺寸对最大、最小孔隙比没有影响;颗粒为圆形时,数值试样的最大、最小孔隙比都比室内试验的结果小,数值试验的结果不能反映实际砂土相对密度的变化规律;颗粒形状对相对密度的影响很大,当颗粒形状为非圆形时,数值试样孔隙比的变化趋势以及孔隙比变化范围随d50的变化趋势都与室内试验的结果基本一致,但是孔隙比变化范围随dmax/dmin的变化趋势与室内试验的结果相反.

摘要： 对于膜结构褶皱效应的模拟近些年来受到广泛关注,恰当的褶皱理论对于准确高效地模拟膜面的力学行为至关重要.本文提出一种基于离散元(DEM)的数值模拟方法,用于分析膜面在褶皱状态下的力学特性.本文将该方法用于一个矩形受剪膜面的静力学响应模拟,此算例被选作是验证模拟方法有效性的基准模型.在模拟中,整个膜面被离散成许多相互不接触的质量元,它们之间由远程作用力相连,以此来模拟面内刚度.并定义了褶皱判断准则,额外添加了由褶皱效应引起的作用力.本文将得到的数值模拟结果与参考文献中的有限元模拟结果进行了对比,得到了较好的一致性.

摘要： 采煤机螺旋滚筒承担着割煤和装煤的双重作用,滚筒的装载是一个很复杂的工作过程,很难用数理模型来评价和计算滚筒的装载效率.基于离散元理论,对滚筒的装载机理进行了分析,并建立了滚筒装载的三维模型,模拟了滚筒的装载,为对滚筒的研究提供了一种新的思路.

摘要： 土拱效应是桩承式路堤中影响荷载传递的关键因素.基于前人的室内模型试验,建立了桩承式路堤离散元(DEM)数值分析模型.基于路堤中应力偏转规律对土拱随桩土相对位移的形成规律进行了分析.在该基础上提出了合理拱轴线土拱模型,并引入荷载传递系数α对土拱与土拱下方路堤填料间的荷载传递进行量化分析.模拟结果表明,桩土相对位移引起路堤中应力主方向发生偏转并形成虚拟土拱;土拱形态及高度随桩土相对位移的变化而变化,最大拱高约08倍桩净间距;α随土拱高度的增加呈对数关系减小.

摘要： 盾构在掘进的过程中,设定合理的支护力是确保隧道开挖面前方土体稳定的关键.利用二维离散元方法,本文研究了动荷载作用下在干砂地层中盾构开挖的合理支护力.首先,通过动荷载作用下隧道开挖面卸荷模型研究卸荷位移S、地层扰动大小与开挖面支护力的关系;其次,通过动荷载作用下盾构掘进模型研究开挖面支护力与地层损失的定量关系,最终从强度控制、地层扰动控制、开挖面地层损失控制三种不同角度对盾构掘进中开挖面所需支护力的大小进行了分析比较,得出基于开挖面地层损失控制的支护压力比最大,为0.99～1.07.并将结果与静载下进行对比,验证动荷载下的盾构施工控制要求更加严格.

摘要： 提出了一种用于分析铁路道床动力特性的离散元-有限元耦合模型.道碴材料具有非连续、非均匀和各项异性等非线性性质,采用离散单元法对其进行分析.由于道碴石子几何形状的不规则性质,道碴颗粒采用镶嵌单元构造.相对于道碴材料,更加密实的路堤和地基采用有限元分析,并在离散元与有限元的接触面处提出一种新的接触算法.利用此模型分析了铁路道床在正弦循环载荷作用下的力学行为.结果表明,枕木的沉降量主要来自于道碴石子的重新排列以及路堤和地基的沉降.此外,道碴中的力链分布于枕木下方,并与其下部结构的应力分布相吻合.同时,道碴材料能够将交通载荷分散并传递至下部结构,降低路堤应力水平,从而提高整体结构的承载能力.

摘要： 煤层气勘探开发过程中井壁稳定性问题突出,传统力学模型难以评价其井眼失效机理.为此,本文在分析破碎性煤层物理性质的基础上,以沁水盆地南部沁端区块3号煤层为例,将虚拟煤岩体定义为基质和裂隙两部分,采用Monte Carlo方法生成了随机的裂隙网络,同时分别赋予相应的岩石力学参数.通过通用离散元程序UDEC进行了破碎性煤层段井壁稳定数值模拟分析,研究了井筒液柱压力对井壁稳定性的影响.研究结果表明,微欠平衡和微过平衡钻进均可以满足煤层段井壁稳定性的需要.研究成果对类似井段的井壁稳定性研究有借鉴意义.

摘要： 引入GPU并行技术,研究了FEM-DEM(有限元-离散元)动态耦合并行计算的数值方法,通过埋没颗粒介质中受迫振动梁的动力学行为来验证方法的有效性.采用面向对象程序设计方法,考虑滚动接触的颗粒接触细观模型并建立颗粒-振动梁接触判断准则,提出了基于GPU(图形处理器)的DEM-FEM耦合运算模型并将 CUDA(统一计算设备架构)超算技术引入并行耦合模拟计算.计算结果表明,该模型能够模拟颗粒物质与结构的动力学耦合作用过程,颗粒介质与简谐激励下振动梁的强耦合作用使原本的线弹性梁呈现出明显的非线性动力学行为,CUDA超算技术的引入可以显著提高大规模耦合计算效率.

摘要： 本发明公开了一种基于离散‑连续耦合的离散元虚拟压缩试验数值仿真方法，首先是建立离散元虚拟数值试验模型，然后利用离散元应力伺服原理，将数值试样伺服到初始固结应力状态，其次利用有限元单元替代侧向刚性墙体，建立离散‑连续耦合数值试验模型，最后施加试验围压与轴向加载速度，开展离散‑连续耦合虚拟数值压缩试验。与现有技术相比，本发明的虚拟数值压缩试验方法可以在试验中实现侧向围压的柔性施加，能够真实地模拟室内压缩试验。通过虚拟数值试验仿真实例，验证了本发明的基于离散‑连续耦合的离散元虚拟压缩试验仿真方法的优势性。

摘要： 本发明涉及一种离散元组合颗粒及其离散元堆积试验模拟方法，包括各种球形度的Clump颗粒，所述的Clump颗粒为由单元球A、单元球B、单元球C和单元球D四个单元球构成，所述的单元球A、单元球B、单元球C的半径相等，所述的单元球A、单元球B、单元球C以圆周分布方式排列于单元球D的四周，且单元球A、单元球B、单元球C均与单元球D相切或相交，所述Clump颗粒的球形度n的取值为0.33≤n≤1。本发明采用离散元方法开发Clump颗粒并进行了堆积试验模拟，可根据颗粒球形度统计结果直接采用实际颗粒的摩擦系数在离散元中模拟砂土颗粒的力学特征，避免了传统圆球颗粒使用严重偏大摩擦系数的缺陷，更逼近实际砂土的力学特性。

摘要： 本发明涉及一种涡轮榫槽喷丸离散元‑有限元耦合多尺度模拟方法，步骤为：(1)根据榫槽结构对称性对欲研究的榫槽结构进行切割，赋予高应变率本构模型；(2)根据实际工艺情况设置喷口与榫槽的相对位置，定义喷口尺寸、喷口角度，完成榫槽与喷口建模；(3)基于高应变率下应力应变数据获取位错演化模型参数，建立榫槽材料的位错演化模型，关联晶粒尺寸、位错密度与宏观应力应变；(4)利用DEM模块定义粒子生成器*Particle，设置弹丸参数，将榫槽与喷口进行组合，设置弹丸与榫槽之间接触特性；(5)进行数值仿真，获取榫槽结构强化之后表面完整性参数。

摘要： 本发明提供了基于有限元与离散元耦合的有砟轨道建模分析方法及系统，方法包括：根据有砟轨道的实际设计施工情况，将有砟轨道简化；获取各结构的真实几何参数建立有限元模型；根据各结构实际采用的材料和物理性能，对各有限元部件赋予材料属性和单元属性，并划分网格；根据各部件的实际约束状态设置相应的边界条件和自由面；施加结构荷载；将模型数据输出为INP文件；将需要被替换为离散颗粒道砟层的部分单独作为一个道砟层模型部件导出，并提取其几何参数，导入至离散元软件中；生成颗粒数据并进行处理；在INP文件中修改数据行，导入处理好的颗粒数据并创建节点集、单元集，定义颗粒属性，并创建接触属性与接触对，得到有砟轨道模型文件。

摘要： 本发明公开了一种有限元‑离散元耦合数值模拟程序(FDEM)输入参数快速标定方法，步骤如下：(1)采用单轴压缩、巴西劈裂和三轴压缩室内试验获得岩石宏观物理力学参数；(2)采用单轴压缩模拟试验峰前阶段曲线标定节理罚值Pf；(3)基于算法改进和单个节理单元压剪破坏标定法向接触刚度Pn；(4)采用单轴压缩模拟试验全过程曲线标定II型断裂能GII；(5)采用直接拉伸试验标定I型断裂能GI；(6)再次采用单轴压缩模拟试验标定切向接触刚度Pt；(7)采用三轴压缩、巴西劈裂和直接剪切模拟试验验算所标定参数的可靠性。本发明提出的标定方法效率高，消除了参数取值对网格尺寸的依赖性，能够适用于工程尺度的模拟。

摘要： 本发明提供了高铁无砟轨道基床表层渗蚀有限元‑离散元耦合评估方法，方法包括：获取车辆‑无砟轨道‑饱和路基系统动力有限元模型的初始水力梯度；计算第n‑1步孔隙率以及基床表层离散元模型的形变；调整车辆‑无砟轨道‑饱和路基系统动力有限元模型的底座板和基床表层接触参数、渗透系数，并计算第n步水力梯度；YADE软件计算n+1孔隙率以及基床表层离散元模型的形变；根据有限元模型的轨道路基等宏观响应参数判别孔隙率变化是否超出阈值并获取基床表层细颗粒损失动态全过程。本发明准确反映底座板与基床表层接触状况，获取基床翻浆冒泥病害细观过程，并将宏观与细观进行耦合评估，有助于高铁运营预警，保障行车安全。

摘要： 本发明公开了一种用于离散元计算即时估算离散介质密实度的方法及存储介质，其方法包括：选择N个代表性抛石，对代表性抛石进行三维轮廓扫描作为抛石体颗粒生成的模板；设置抛石体的投放区域，利用自重下落法随机生成抛石体颗粒，进行分层自重沉积，构成自然堆积的抛石体三维基床模型；采用锤夯/振动模拟方法对抛石体三维基床模型进行密实，在间隔时间或间隔时步内更新密实区域地表面貌的变化，估算出密实区域的局部密实度变化与整体区域的密实度变化。本发明通过设计工程的需求，在PFC3D计算平台中模拟荷载施加后对抛石体密实度的影响，利用密实度与荷载施加的模拟工艺参数的关系，来指导实际工程参数优化，提高工程实施工艺的合理性和准确性。

摘要： 本发明公开了一种基于离散‑连续耦合的离散元虚拟压缩试验数值仿真方法，首先是建立离散元虚拟数值试验模型，然后利用离散元应力伺服原理，将数值试样伺服到初始固结应力状态，其次利用有限元单元替代侧向刚性墙体，建立离散‑连续耦合数值试验模型，最后施加试验围压与轴向加载速度，开展离散‑连续耦合虚拟数值压缩试验。与现有技术相比，本发明的虚拟数值压缩试验方法可以在试验中实现侧向围压的柔性施加，能够真实地模拟室内压缩试验。通过虚拟数值试验仿真实例，验证了本发明的基于离散‑连续耦合的离散元虚拟压缩试验仿真方法的优势性。

摘要： 本发明涉及一种离散元组合颗粒及其离散元堆积试验模拟方法，包括各种球形度的Clump颗粒，所述的Clump颗粒为由单元球A、单元球B、单元球C和单元球D四个单元球构成，所述的单元球A、单元球B、单元球C的半径相等，所述的单元球A、单元球B、单元球C以圆周分布方式排列于单元球D的四周，且单元球A、单元球B、单元球C均与单元球D相切或相交，所述Clump颗粒的球形度n的取值为0.33≤n≤1。本发明采用离散元方法开发Clump颗粒并进行了堆积试验模拟，可根据颗粒球形度统计结果直接采用实际颗粒的摩擦系数在离散元中模拟砂土颗粒的力学特征，避免了传统圆球颗粒使用严重偏大摩擦系数的缺陷，更逼近实际砂土的力学特性。

摘要： 本发明公开了一种基于有限元‑离散元耦合的岩石破裂声发射模拟分析方法，通过获取岩石破裂时当前时间步的声发射数据，根据所述声发射数据获得声发射事件信息，所述声发射事件信息包括声发射事件的状态和声发射事件包含的微裂纹对应的微裂纹信息；基于声发射事件的微裂纹信息，获得所述状态为非激活态的声发射事件的破裂区域，计算所述破裂区域内的所有有限单元的节点力变化分量；基于所有的所述节点力变化分量，获得力矩张量；基于所述力矩张量，获得宏观裂纹的分析结果。与现有技术相比，能够区分裂纹的类型并确定裂纹的变形运动特征。