最近在用近场动力学Peridynamics, PD做一些冲击破坏的问题。这次跑的是薄板中心冲击的算例计算了 300 步能观察到裂纹已从冲击中心向四周充分扩展。当然用 PD 处理这类强动态断裂问题最突出的优势就是无需预设裂纹路径也不用做网格重构从冲头接触板件、应力波扩散传播到损伤萌生、多裂纹径向延伸的全过程都能自然模拟很适合研究冲击载荷下的断裂演化机理。我目前用的这套 PD 计算框架从数值计算到后处理的链路很完整不用自己从零搭建后处理脚本就能拿到很系统的量化结果对做参数分析和机理研究省了不少重复工作。咋说呢现在用的这套代码基础场量输出很完备三维位移模量云图可以直观呈现冲击中心的变形集中效应与应力波扩散形态全场位移分布、最大位移量都能直接读取损伤结果既支持三维全域展示也可输出厚度方向投影的损伤云图还能自定义损伤阈值本次取 0.5自动统计损伤面积与占比该时刻损伤面积比约 35.95%定量评估破坏程度非常方便很实用的一点是自带裂纹骨架提取功能可自动识别出 7 条主裂纹分别标注每条主裂纹的长度与延伸角度同时统计总裂纹长度把离散的损伤场直接转化为可定量分析的断裂几何参数省去了自己用图像处理提取裂纹的步骤动态历程数据也很完整损伤面积、裂纹长度的时间演化规律冲头的位移与速度变化以及接触力、系统动能的变化曲线都能自动输出不管是校验能量守恒性还是对比不同工况下的断裂响应都有充足的量化依据额外附带位移、速度、加速度、PD 体力等物理量的分布直方图做结果收敛性校验、场量统计特征分析的时候很省心。我自己也是边学边摸索从影响域半径、时间步长、键强度的参数标定到后处理结果的工程解读踩了不少坑也攒了一些实操经验。整体感觉这套框架对于想入门 PD 冲击断裂研究、不想重复造轮子的同学来说上手门槛友好分析功能也足够扎实。等我后面建个群吧感兴趣的 pd 小伙伴可以进群大家可以一起讨论。