离散元模拟

摘要： 页岩作为页岩气储层,在沉积过程中形成部分弱面,在力学特性上表现出各向异性特征.所以,使用离散元软件从微细观层面探讨深部页岩力学各向异性特征具有重要实践意义.基于页岩室内常规三轴压缩试验结果,采用离元程序PFC2D对常规三轴压缩下不同层理倾角页岩进行了颗粒流模拟研究,分析了层理倾角及围压对页岩力学特性的影响规律.结果表明:(1)页岩峰值强度与黏聚力随层理倾角的增加整体呈"U"形变化,但峰值强度在不同围压下的变化趋势有所区别;而内摩擦角随层理倾角的增大呈非线性变化.(2)层理倾角对页岩周围颗粒的位移方向及大小的影响随着层理面与轴向应力的夹角的增大而减小.(3)同一层理倾角试样最终破坏时的微裂纹总数随着围压的升高有所增加;同一围压下,试样最终破坏时的微裂纹数目,随着层理倾角的增加呈现先减少后增多的趋势.(4)同一层理倾角页岩的脆性随围压的增长整体呈下降趋势;低围压情况下,页岩脆性随层理倾角的增加呈两端大中间小的变化规律.

摘要： 反转构造在我国含油气盆地中广泛发育,对于油气勘探具有重要意义。本文采用离散元颗粒流方法研究了先存正断层的陡缓程度、倾向组合关系和距离挤压端的远近等因素对断层的反转构造发育及其反转量的影响。模拟结果显示,先存断层的产状及其距离挤压端的远近会影响先存断层的反转次序、反转量,进而影响盖层内反转构造的发育。当先存正断层倾向指向挤压端时,若缓倾正断层靠近挤压端,则缓倾断层发生反转,而远离挤压端的陡倾断层未发生反转;若陡倾断层靠近挤压端,则缓倾断层和陡倾断层同时发生反转。当先存正断层倾向背离挤压端时,无论靠近还是远离挤压端,缓倾断层和陡倾断层均发生反转。当先存正断层倾向相对(地堑)时,若缓倾正断层更靠近挤压端,则缓倾断层发生明显反转,而陡倾断层未发生反转;若陡倾正断层靠近挤压端,缓倾断层和陡倾断层同时发生反转。在先存断层均发生反转的模型中,无论先存断层倾角大小,靠近挤压端的先存断层的反转量最终均大于远离挤压端的先存断层的反转量。最终,离散元模拟结果还与西湖凹陷垂直于反转构造带走向的地震剖面进行对比,有很好的吻合性。

摘要： 为深入研究土工格栅的加筋机理,基于前人的研究成果,阐述拉拔试验中影响土工格栅与筋土界面特性的因素,介绍土工格栅和土体离散元建模以及拉拔荷载下筋土相互作用的可视化和量化分析结果,并提出有待深入研究的方向。研究表明:在试验研究方面,研究人员多采用传统试验方法,从试验边界条件、土工合成材料类型、土料类型及参数、加载条件等方面较为系统地研究各因素对筋土相互作用的影响;随着测试技术的发展,新型试验仪器可有效弥补传统测试方法的不足,并实现筋土界面作用分析的可视化。此外,离散元数值模拟是试验研究的重要补充,大量模拟结果可从细观层面更为深入地阐述筋土界面的作用机理。

摘要： 为了研究充填支架夯实机构的夯实角对矸石散体夯实效果的影响,以ZC1800/14/25型充填液压支架为原型,采用离散元数值模拟的方法构建矸石充填工作面空间模型,分析了推压板和矸石散体在夯实过程中的力链演化,研究了夯实角与夯实板推压力、矸石破碎程度、欠接顶量和矸石质心位置之间的关系。研究结果表明:随着夯实角的减小,夯实作用会影响到较大范围的矸石并压碎部分矸石,夯实板压力较大;随着夯实角的增大,夯实板压力降低,且推动的矸石量较少。最终确定夯实角为10°时,夯实机构对矸石散体的充填效果最好,研究结果可为充填工作面现场操作提供理论依据。

摘要： 滚石的运动特征研究是滚石灾害防治的重要前提,而坡表覆盖层和滚石自身特性是影响其运动特征的重要因素。以斜坡覆盖层材质为变量,开展物理模型试验,研究了混凝土坡面和土质坡面下滚石试块的运动特征;并在此基础上,以在建合壁津高速公路边坡为例,开展滚石运动特征数值模拟研究,以运动速度、转动速度、运动能量等为指标分析了滚石形状和尺寸对其运动过程的影响。结果表明:形状和尺寸对滚石运动特征影响较大,滚石越接近球形、尺寸越小,其运动速度越大、转动速度越大;滚石在下落过程中有30%以上的势能通过坡面碰撞损耗,在滚石防治工程中应考虑这部分损耗,而滚石转动能量只占总动能的1/30~1/10,因此一般不单独考虑滚石的转动过程。研究成果可为滚石防治工程的设计及优化提供参考。

摘要： 我国西南高山峡谷地区,常发生高位岩质斜坡失稳,形成碎屑流,造成巨大人员伤亡与经济财产损失。目前存在着许多关于岩质滑坡的研究理论和物理模型,但这些理论与模型都忽略了岩体结构这一属性,然而对岩质滑坡而言,岩体的结构特征不仅控制着斜坡变形失稳模式,同样影响失稳后滑体碎屑化运动过程,甚至最终的危害范围。因此,文章使用离散元数值模拟方法,通过模拟不同结构面强度、结构面密度、结构面方向和岩块强度下岩质滑坡碎屑化运动过程,研究了不同源区岩体结构对岩质滑坡破碎特征、运动形式和运动距离的影响规律及影响机制。研究结果如下:(1)质心运动距离随结构面强度增大而增大,结构面抗拉强度提升10倍,质心运动距离增加3%;(2)结构面密度的增加使得块体破碎率随之增加,但运动距离和分布面积都呈下降趋势;(3)水平向结构面岩体的前缘运动距离和质心运动距离与其他方向结构面工况相比均降低10%左右,分布面积则缩减接近30%,破碎率最大;(4)随着岩块强度增加,破碎率降低,分布面积最终缩减40%,完整块体工况运动距离最终增大15%左右。结果有助于进一步理解岩质滑坡动力学过程,指导山区防灾减灾工作。

摘要： 为探究各向异性对颗粒堆积体力学特性的影响,采用MatDEM离散元软件模拟了不同围压条件下的三轴排水试验,探究了砂土宏细观力学行为及其临界状态特性,利用组构张量对试样的各向异性进行量化分析,验证了应力-力-组构的相关关系。利用固结、剪切破坏后具有不同初始各向异性的试样,采用边界去除法进行堆积试验,分析了堆积体的休止角、基底压力分布、力链、各向异性等宏细观特性。结果表明:试样的初始各向异性程度越大,堆积体的休止角越大;剪切试样的初始应力拱及组构各向异性导致其堆积体的基底中心压力下降;堆积体的法向接触力各向异性对堆积体宏观特性起着主导作用。

摘要： 受结构面的控制,横观各向同性的现象在岩体工程中十分常见,其对于岩体的强度性质以及破坏行为有较为重要的影响。目前,对于岩体的强度估算,一般采用Hoek-Brown强度准则,但传统的Hoek-Brown准则缺乏对岩体横观各向同性性质的较为清晰的描述。为此,该文通过引入岩石体积节理比J_(v),节理断续度J_(c),岩体摩擦比系数J_(f),节理平直度J_(s)和节理弱化乘子γ5项指标,对Hoek-Brown强度准则进行了修正,并通过数值模拟验证,能较好地预测横观各向同性岩体的强度,为实际工程设计提供了理论依据。

摘要： 混凝土的粗骨料与砂浆之间存在界面过渡区(ITZ),过渡区的力学性质受骨料的形状、大小、含量、空间分布及初始裂隙等条件影响。为研究过渡区强度对混凝土材料性能的影响,采用离散元软件PFC2D建立了随机骨料分布的混凝土模型,并对其进行了单轴压缩试验。试验发现:当过渡区强度增大时,混凝土的峰值应力和初始裂隙产生时的应变逐渐增大,总裂隙数减少,破坏类型由受剪破坏转变为受拉破坏,并且过渡区强度增强会使混凝土的起裂区域由过渡区转移至砂浆聚集区。模拟结果可为混凝土界面过渡区的理论研究提供一定参考。

摘要： 为给能源桩的换热特性和承载能力特性的研究提供新思路,文中使用南京大学开发的离散元软件MatDEM,通过粒径优化方法逐层堆积离散元颗粒,建立单U型换热管的能源桩的三维传热模型,使用PHC能源桩现场实验对模型进行修正,模拟分析的桩体温度、桩周土体温度、桩身应力应变、土体位移,与现场实验结果对比结果一致,体现了单U型换热管的能源桩三维传热模型的合理性,为能源桩的离散元研究提供基础。

摘要： 沉积结构面的存在显著影响页岩等岩石沿不同方向的物理和力学特性.为研究层理对岩石水力压裂响应的影响,本文开展了基于微观结构的各向异性岩石水力压裂半定量离散元模拟.通过在由平行黏结模型组成的各向同性岩石中插入一系列光滑节理显式表征各向异性岩石的微观层理结构.采用颗粒尺度上的管域模型描述岩石的水力耦合行为.分别开展了各向同性和各向异性岩石在无围压和有围压条件下的水力压裂过程模拟.结果表明,各向同性岩石在无围压状态下,水力压裂裂缝扩展表现出一定的随机性,裂缝的起始和扩展主要沿着最弱黏结处进行.边界应力即围压的存在控制了各向同性岩石的裂缝沿着最大主应力方向发展.该结果与断裂力学预测结果一致.在无围压条件下,各向异性岩石的水力压裂裂缝沿着结构面的方向阶梯状发展.有围压条件下,各向异性岩石的水力压裂起始裂缝沿层理方向扩展至一定距离后转向最大主应力方向.本文证实了基于微观结构的离散元模型模拟页岩水力压裂的可行性,为后续建立定量的多尺度模型提供了基础.

摘要： 厚煤层(煤层厚度≥3.5m)产量占我国煤炭总产量的45％左右,其中,综合机械化放顶煤开采(综放开采)作为开采此类煤层的主要方法之一,一直面临着顶煤采出率有待提高的问题,其根本原因在于对顶煤破碎机理和放出规律缺乏深层次掌握,顶煤放出规律属于非理想松散介质的流动问题,目前以采用物理模拟试验与数值模拟计算为主.本文介绍了采用颗粒流程序PFC2D和PFC3D进行综放开采煤岩散体流动特性模拟的一般方法和步骤,并基于BBR研究体系,阐述了顶煤放出规律的主要研究内容,分析了影响煤岩散体流动特性的主要因素;最后对比了采用二维和三维离散元模拟研究煤岩散体流动特性的合理性和优缺点,提出了2.5维离散元快速模拟综放开采散体煤岩放出过程的方法.

摘要： 路面堵塞模型建立，常规排水路面堵塞分析，双层排水路面抗堵塞分析，随着渗水系数的增加，双层排水路面上层路面结构更快的达到最终堵塞状态，但总体的堵塞体积不变，下层路面结构更快的达到最终堵塞状态，但总体的堵塞体积减小。

摘要： 在山区,地震诱发滑坡造成的人员伤亡往往占到地震总伤亡人数的一半以上,尤其是在黄土地区,由于地震造成的黄土滑坡具有运动距离远等特征,其灾难性更严重.针对地震诱发黄土斜坡的失稳峰值加速度、临界位移和运动距离等问题,本文利用离散元(PFC)方法,通过对室内三轴实验的应力-应变曲线进行标定和高精度航拍三维地形数据进行转化,以宁夏海口村黄土斜坡为研究对象,开展三维地震作用下斜坡失稳破坏和运动过程的数值模拟研究.通过监测不同位置在地震作用下的应力分量,计算监测点的p、q值并与室内三轴实验所得到p-q破坏线进行比较,获得斜坡破坏过程中的应力路径.结合颗粒的监测位移,获得斜坡失稳破坏的临界位移,得出该黄土斜坡失稳峰值加速度为0.135g,临界位移为50cm.同时根据不同地面环境条件,预测了在不同摩擦系数下斜坡失稳后的危害范围,为黄土地区边坡工程的抗震设计及防震减灾工作提供一种新的可视化方法.

摘要： 堆石料广泛应用于土石坝和山区机场地基的填筑工程中,其湿化变形可能会影响工程的正常使用甚至威胁其安全性.颗粒破碎是导致湿化变形的重要原因,然而目前对于堆石料湿化过程中颗粒的破碎规律尚不明晰.本文采用离散单元法,分别模拟了侧限条件下单线法和双线法湿化变形试验,比较了单、双线法湿化变形的差异.采用Hardin破碎率表征颗粒破碎量,分析了湿化过程中颗粒破碎量与填筑体应变及颗粒间摩擦耗能等的关系.此外,还讨论了湿化路径、颗粒形状、材料软化系数对湿化变形的影响.研究结果表明:侧限条件下湿化变形随应力水平增加而增加,且单线法测得的湿化变形大于双线法.湿化路径对总变形量有一定影响,先湿化后加载可有效降低总变形量.加载过程中,颗粒破碎以棱角破碎为主,且多棱角颗粒的试样湿化变形大于浑圆颗粒试样.Hardin破碎率与试样摩擦耗能的关系受湿化路径的影响,与应变的关系不受湿化路径的影响,仅与试样初始状态和材料特性有关.

摘要： 散粒体的级配显著影响其变形特性.通过建立二维离散单元模型,选取不同围压和初始密度,对不同级配的散粒体在循环荷载作用下的累积变形进行了研究.研究发现:在相同初始孔隙率或者初始相对密度条件下,随着级配不均匀系数的增加,散粒体变形均呈现增长趋势;增加土体的初始密度和围压可以显著增加其内部配位数,从而有效抑制后期变形.在低围压下,所有试样均发生不同程度的剪胀,变形随着加载次数的增加而快速增长;当围压增加时,试样由剪胀转为剪缩,变形减小、趋于稳定.

摘要： 自密实混凝土在重力坝施工过程中受施工工艺、施工质量及施工环境等因素的影响,砂土颗粒间的胶结强度可能各不相同.为研究不同胶结强度的胶结砂土的力学特性,借助PFC3D自带的微观胶结模型,将胶结砂土中的胶结强度的分布定义为正态分布,并系统地分析了在不同方差和胶结含量下胶结砂土的力学特性.研究表明,胶结试样在数值试验中表现出应力应变软化特性及体变剪胀性,且随着胶结含量的增大,软化及剪胀程度提高;当偏应力达到峰值时,随着轴向应变的增大,体变由剪缩变化为剪胀.在假定的胶结砂样的胶结强度服从正态分布的情况下,峰值应力随着胶结含量成线性增大,但正态分布方差的变化对胶结砂样强度的影响不大.

摘要： 为了研究颗粒棱角对颗粒材料力学行为的影响,建立了具有不同棱角度的对称多面体颗粒,采用了一种简单并适合任意颗粒形状的接触本构模型,对三维离散元开源程序YADE进行了修改,研究了颗粒棱角度在模拟直剪试验中的影响以及接触力各向异性在剪切过程中的演化规律.研究结果表明,颗粒棱角度越小,颗粒间相互咬合自锁的作用越小,颗粒受剪更易转动,致使颗粒体系的剪切强度和剪胀性下降;竖向加载力越大,颗粒棱角度的影响越明显;法向接触力的各向异性在剪切过程中表现为先增后减最后趋向稳定的趋势;法向接触力的各向异性变化程度随颗粒棱角度的增大而增大.

摘要： 天然气水合物以胶结及孔隙填充等形式存在于深海能源土中,开采时因其分解会劣化地层力学特性进而引发海底事故,使得人们对能源土开采过程进行中力学特性的变化愈发重视.在前期室内试验的基础上,将一个温度-水压-力学二维微观胶结模型引入离散元商业软件PFC2D中,通过对排气、排水性较好的土体进行升温及降压法开采进行数值模拟,并将模拟结果与相同条件下的室内试验结果对比,验证了该胶结模型的适用性.进一步分析了颗粒接触分布与颗粒平均纯转动率(averaged pure rotation rate,APR)在水合物分解时的变化情况.升温分解过程中随温度升高,颗粒总接触分布各向异性程度增大;胶结接触逐渐减少并始终保持主方向为水平方向,无胶结接触增多并始终保持主方向为竖直方向;APR值逐渐增大且正负值分布逐渐趋于集中.降压分解过程中随反(水)压降低,颗粒总接触由各向同性分布逐渐发展为主方向为竖直方向的各向异性,APR值较小且分布均匀;恢复反压后,试样进一步破坏,颗粒总接触各向异性更加明显,APR值增大且正负值呈集中分布.

摘要： 以汶川地震触发的宝成铁路109隧道滑坡为例,在野外现场调查基础上,采用离散元数值模拟技术,输入实际记录地震波,对该边坡在地震作用下的动态响应及塑性区发展过程进行数值模拟,以研究地震作用下节理岩质边坡的失稳机理.并结合强度折减法,求出了在地震作用下该边坡的安全系数.分析结果表明:在地震作用初期,边坡表面出现拉破坏,并沿着坡面由边坡的下部向上部延伸,随着地震作用进行,拉破坏区逐渐向坡体内部发展,边坡下部的拉破坏区域逐渐扩大且向坡体上部发展,在地震作用后期,边坡表面存在着强烈的拉裂作用,当塑性拉破坏区发展到形成一个贯通的区域时,整个贯通区域上部的坡体便会发生滑塌破坏.

摘要： 本发明涉及一种离散元组合颗粒及其离散元堆积试验模拟方法，包括各种球形度的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)进行数值仿真，获取榫槽结构强化之后表面完整性参数。

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

摘要： 本发明公开了一种基于有限元‑离散元CT重构煤体的渗流模拟方法，通过工业CT扫描得到煤体的三维数据体，而后进行三维重构并去除“孤岛块体”，以得到可被仿真软件所调用的通用三维几何模型。通过采用离散元仿真方法，编写煤体断裂判别算法，确定煤体是否在水压应力作用下发生变形甚至断裂，同时基于有限元模拟方法，运用移动网格方法，实现断裂后的煤体碎块在煤体内随水分流动而运移，同时进行水分渗流压力场、渗流速度场及煤体水分增量的数值模拟，最大程度的展现多孔介质煤层在注水过程中煤体的变形与断裂现象及水分的运移规律。

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

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

摘要： 本发明公开了可用于施工过程模拟的有限元离散元多尺度耦合计算方法，属于岩土体工程技术领域，本发明革新了现有的有限元离散元多尺度耦合计算框架，构建了可变的离散元子结构序列Sdem，使之与施工过程中宏观有限元网格的减少或者增加对应起来，从而实现施工过程的模拟。与现有技术相比，本发明扩展了有限元离散元多尺度耦合计算方法的使用范围，可以更好地预测复杂岩土体材料在开挖、堆载等过程中的力学行为，为相关岩土石力学研究服务和工程安全性评价提供技术支撑。

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

摘要： 本发明公开了一种有限元‑离散元耦合(FDEM)恒阻大变形锚杆隧道围岩加固数值模拟方法，步骤如下：①通过设置功能型四边形节理单元实现恒阻装置在FDEM中的表征；②除恒阻单元外的其余单元采用线弹性本构方程，提出了等效弹性模量和等效剪切模量计算方法；③采用浇筑方法，将锚固剂浇筑成直接拉伸和变角剪切试样，获取锚固剂的拉伸与剪切本构方程；④采用锚杆拉拔试验和锚固剂‑岩石界面拉伸与剪切试验获得锚固剂‑锚杆和锚固剂‑岩石界面本构方程；⑤采用破裂节理单元粘结修复算法，实现锚杆滞后安装的模拟。本发明公开的恒阻大变形锚杆隧道围岩加固FDEM数值模拟方法，为恒阻大变形锚杆隧道围岩加固效果预测和加固机理研究提供了数值模拟工具。

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