ABAQUS边滚刀破岩控制仿真，关键参数怎么调

控制边滚刀（存在安装角度）公转且自转进行破岩，岩石材料DP，后续可加cohesive。

安装角度8°Example

Standard/Explicit 里 Job 提交完，status 里蹦出 “too many attempts made for this increment”，滚刀还没啃进岩石就停了，接触刚度、质量缩放、DP 屈服面三件事同时没调好才会这样。

边滚刀有安装角，刀刃法向不沿径向，公转靠参考点 RP 拖着刀座绕掌子面中心转，自转靠刀体与岩石切向接触反力自然带出来，别用刚体约束把自转也写死，写死之后破岩轨迹变成车螺纹，裂纹只沿一条线，岩屑不剥落，DP 材料塑性功算出来虚高。

参考点建在刀轴中心，用 Rigid Body 把刀体网格绑到 RP，RP 上施加旋转幅值曲线，公转角速度 ω，安装角 β 直接体现在刀座局部坐标系里，刀刃线在局部系画好再 rotate 到全局，别在全局里手挪节点，挪完质心偏了离心力项在 Explicit 里会把刀甩出接触面。

岩石 DP 参数，dilation angle 别上来填 15°，先 0° 跑通再慢慢加，摩擦 angle 按岩石三轴试验峰值减 2~3°，flow stress ratio 1.0 起，yield surface 用 *DRUCKER PRAGER 里默认 linear，Cohesion 先给单轴抗压的 0.3~0.5 倍，Young’s 模量别用实验室静态值，动态破岩用 0.6~0.8 倍，Poisson 0.25 凑合，密度按真实值，Explicit 里质量缩放只加在岩石最小单元，缩放因子控在 1e-7 秒级，再大惯性效应把裂纹方向带偏。

接触，刀-岩用 Surface-to-Surface Contact，法向 hard contact，切向 penalty，摩擦系数 0.15~0.3，边滚刀实际是滚动+滑移，μ 给 0.2 先跑，切向剪切里开 finite sliding，刀刃前角那一块网格 0.5~1 mm，岩石侵彻区三层扫掠，最外层 4~6 mm 过渡，单元 C3D8R 但关掉沙漏控制里的默认 stiffness，改 enhanced，hourglass 能量超 5% 就加密刀刃前方。

*AMPLITUDE, NAME=SPIN, 公转幅值给平滑阶跃，前 0.02 s 从 0 拉到目标 ω，别阶跃，阶跃一开始接触穿透 2 mm，接触对直接发散，自转不用单独幅值，接触反力矩 / 刀体绕轴惯量 自然出角加速度，刀轴处加个轴向小阻尼防止自转过冲，damping 系数 1e-4 量级，大了滚动变拖滑。

后续加 cohesive，先别用 COH3D8 铺满，用 VCCT 或者 bilinear traction 的 COH2D4/COH3D8 只在预切槽和节理面铺，cohesive 初始刚度 Kn 取岩石 E 的 10~50 倍，最大牵引力等于岩石抗拉强度，断裂能 Gf 用巴西劈裂算出来的 20~80 J/m²，边滚刀路过节理时若安装角使刀刃切向分力指向节理张开方向，裂纹直接沿弱面跑，不扩径，这种工况要把安装角往掌子面切深方向再倾 2~3°。

inp 里刀体那块：

*RIGID BODY, ANALYTICAL SURFACE=NO, REF NODE=RP-1

*ROTARY INERTIA, RP-1, I11, I22, I33

I11 绕刀轴，I22/I33 绕径向和切向，给错一个数量级自转启动慢半拍，破岩第一个齿痕浅一截，后面才稳态。

输出里 F-刀岩接触总力、RF-RP 反力、SDEG-cohesive 损伤、PE-DP 塑性应变、STATUS-单元删除状态全开，LE 应变别只看 U，边滚刀安装角导致刀体侧磨，U 里看不出，PE 在刀刃侧后角先起来才是装角对了的信号。

质量缩放和 DP 非线性一叠加，稳定步长掉到 1e-9，显存吃满，把岩石远场单元用 bulk viscosity 调 0.06/1.2 默认别动，近场才吃小步长，远场用 coarsen + tie，tie 界面别放在裂纹扩展主路径上，放上去裂纹到界面就拐弯，假象张应力。

安装角 β、公转半径 R、刀刃楔角 α 三个量，先固定 α=20°、R 按刀盘实际，只扫 β：0°（纯滚）、5°、8°、12°，β 超过 10° 岩石表面出现侧向犁沟，DP 里 shear failure 比例上去但体积破碎不下去，掌子面只起皱不起屑，这种 β 就得往回收，不是越大贯入越快。

接触对里把刀体设为主面，岩石为从面，从面单元不能比主面大，边滚刀刀刃 R 角那块主面曲率小，岩石单元再大就穿透，穿透超 0.1*t 接触力振荡，DP 应力回映在穿透区算出负静水压力，裂纹从没接触的地方冒出来，云图像开花。

Step 时间，破岩 Explicit 总时长按 公转一圈 / 3~5 个齿距 来给，别按真实掘进速度跑整圈，整圈显式算到地老天荒，看单个齿痕成形、裂纹网络、岩屑剥离就够了，稳态后再说延长。

DP 材料里 *DEPENDENCIES 先不挂温度和应变率，挂了要补 *RATE DEPENDENT 和 Johnson-Cook 类参数，边滚刀冲击区应变率 10²~10³，不补就低估强度，补了又和 DP 非关联流动打架，先稳态再上率效应。

刀体材料用线性弹或者简单弹塑，别给 H13 热处理后完整硬化曲线，滚刀本身不变形是前提，给硬了接触求解还稳，给软了刀体缩进去岩石没破刀先瘪，反力曲线全是刀体压缩，破岩力被吃掉一小半。

cohesive 加完第一次跑别看应力，先看 SDEG 从 0 到 1 的空间顺序，边滚刀先压密、再起裂、再沿节理或放射状主裂纹扩，这个顺序对了再调 Gf，顺序乱了就是参数权重歪了，摩擦角、胀裂角、cohesive 刚度三者里至少有一个在和物理作对。

免责声明：本文系网络转载或改编，未找到原创作者，版权归原作者所有。如涉及版权，请联系删