EDEMpy入门：用Python接口实现离散元仿真自动化

简介面向颗粒动力学模拟与EDEM二次开发的学习者这份EDEMpy技术资料包系统整理了从离散元基础到Python接口应用的完整内容可帮助用户解决颗粒系统建模、模拟控制与数据后处理中的实际问题。全部资料共14个文件压缩包整体仅5.16MB以示例工程压缩包、打包源码、安装包、PDF文档和Python脚本为主涵盖发布说明与入门指南、几何输出脚本以及粘结颗粒、分离指数、停留时间、压缩、螺旋输送、犁铲等多类典型模拟案例。已有855人学习下载适合希望借助Python定制EDEM模拟流程、实现批量自动化分析的工程师与研究人员。资源内附的示例工程和源码可配合PDF文档系统学习接口调用、边界条件设置、结果提取与可视化展示同时通过输送系统、粉碎机、搅拌器等实际场景帮助读者把技术要点快速迁移到自身课题中。整体定位紧凑兼顾理论原理与工程实践是了解和使用EDEMpy时的实用参考资料。1. EDEMpy是什么把EDEM颗粒模拟从GUI拖进Python脚本EDEMpy是EDEM官方为颗粒离散元仿真提供的Python接口它让原本依赖鼠标点击的建模、运行和后处理工作变成一串可复现的命令。在批量测试参数、做CFD-DEM耦合或者搭建自动化优化流程时这个接口能省掉大量重复操作。这套技术资料里包含了0.1.1和0.1.3的源码包、Release Notes、Getting Started以及BondedParticle、SegregationIndex、ResidenceTime、Compression、ScrewAugerTorque、Plough等六七个示例项目适合刚接触EDEMpy的人从零搭环境也适合已经有EDEM基础的人直接抄作业。需要说明的是EDEMpy并不是把EDEM的全部功能都暴露给Python它的核心价值在于把“模型参数—求解控制—结果提取”写成脚本让颗粒动力学研究可以批量开展。2. EDEMpy环境搭建与API调用从tar.gz到第一个DEM仿真脚本2.1 版本选择与环境依赖拿到压缩包后先不要急着装。EDEMpy和EDEM主版本有对应关系0.1.1和0.1.3对应的EDEM版本范围不同。我一般会先打开EDEMpy-0.1.1_ReleaseNotes.pdf确认当前安装的EDEM主版本是否在兼容列表里。旧版EDEMpy对Python版本也比较挑剔0.1.x时代常需要在Python2.7或者较老的Python3环境里运行原因在于官方发布时绑定的rpy2和sip版本比较旧。安装方式很简单推荐直接用pip从本地源码包安装pip install edempy-0.1.3.tar.gz如果不想污染全局环境也可以用虚拟环境python -m venv edemenv edemenv\Scripts\activate pip install edempy-0.1.3.tar.gz安装完成后检查导入是否正常import edempy print(edempy.__version__)如果输出版本号说明底层依赖没问题。常见错误是提示找不到sip或PyQt5因为EDEMpy在部分版本里要借助Qt的库和EDEM主程序通信此时安装对应依赖后重试。资源里的EDEMpy_1.2.rar和edempy-0.1.x.tar.gz在API上并不完全兼容下文代码以0.1.x为主1.2版的调用方式要用dir(edempy)自行核对。下面是一个简单对照表文件内容用途edempy-0.1.1.tar.gz0.1.1 Python源码包兼容较低版本EDEMedempy-0.1.3.tar.gz0.1.3 Python源码包修复部分反馈优先使用EDEMpy_1.2.rar1.x版本分发包API有变化适合新版EDEM2.2 启动EDEM并连接EDEMpyEDEMpy的启动逻辑通常是先启动EDEM主程序再通过本地通信接口把Python进程和EDEM进程绑定。最简单的写法是直接用EDEMApplication指定可执行文件路径import edempy dm edempy.EDEMApplication(rD:\Program Files\EDEM\EDEM.exe) dm.Start()参数是EDEM安装目录下的可执行文件如果你已经手动打开了EDEM也可以不传路径直接连接到当前实例。Start()返回后Python进程就拿到了这个EDEM实例的控制句柄。这一步最常见的报错是找不到EDEM可执行文件或者端口被占用建议启动后等待几秒再执行后续指令。然后是加载模型和运行模拟。EDEM中建模完成的几何和颗粒工厂保存在.dem文件里可以先用EDEM GUI建好基础模型再用EDEMpy加载运行也可以完全用代码创建。下面是一个典型的加载并运行流程sim dm.CreateSimulation() sim.LoadModel(hopper.dem) sim.Reset() sim.RunSimulation(totalTime1.0, timeStep0.001) sim.Close()LoadModel的参数是模型文件路径Reset把颗粒状态恢复到初始时刻重复试验时很有用RunSimulation的totalTime是模拟总时长单位秒timeStep是固定时间步长实际取值要和EDEM的Rayleigh时间步对齐。时间步太大颗粒穿透严重太小计算量指数上升。2.3 核心API结构了解启动逻辑后再看EDEMpy的API结构。0.1.x版的接口基本分为三层Simulation层负责加载模型、运行、暂停、恢复Model层负责访问几何体、颗粒、接触模型和数据Data层负责从结果文件提取颗粒位置、速度、力、能量。下面用一个代码片段串起来model sim.GetModel() for particle in model.Particles: pos particle.GetPosition() if pos.z 0.1: v particle.GetVelocity() print(particle:, particle.GetId(), pos, v)GetModel()拿到当前模型的根对象Particles是已有颗粒的迭代器GetPosition()返回坐标为(x, y, z)的listGetVelocity()返回颗粒在三方向的速度。这里只打印高度小于0.1米的颗粒可用于快速过滤某个层位。EDEMpy的数据接口返回的是列表不是numpy数组要算统计量时建议自己转成np.array后续分析才顺手。3. 颗粒模型与边界条件用EDEMpy定义物料属性和相互作用3.1 颗粒形状与尺寸分布EDEM的颗粒模型分为球体、多球聚合和自定义多边形。EDEMpy可以读取并修改已有几何中的颗粒类型也可以新建颗粒工厂。比如给一个给料口设置正态粒径分布factory sim.AddParticleFactory(feed) factory.SetParticleShape(sphere) factory.SetRadius(0.005) factory.SetDistribution(normal, 0.005, 0.0005)SetDistribution第一个参数是分布类型EDEMpy内置了normal、uniform等第二个参数是平均值第三个是标准差。如果要用粘结颗粒比如资源包里的BondedParticle_Example.zip还需要在颗粒层面定义键数目和键面积EDEMpy中通过SetBondedContacts控制。多球聚合模型适合形状不规则的矿石计算成本比球体高不少不要一开始就用。3.2 材料属性与接触模型颗粒模拟的结果高度依赖材料本构参数EDEMpy的好处是这些参数都可以脚本化适合做参数标定。常用材料属性如下表参数对应EDEMpy方法单位说明密度SetDensitykg/m³影响重力与惯性剪切模量SetShearModulusPa越小越容易发生接触变形泊松比SetPoissonRatio无影响接触面积恢复系数SetRestitution无碰撞反弹程度静摩擦系数SetStaticFriction无影响堆积角滚动摩擦系数SetRollingFriction无影响动态休止角代码示例coal sim.AddMaterial(coal) coal.SetDensity(1400) coal.SetShearModulus(1e7) coal.SetPoissonRatio(0.25) steel sim.AddMaterial(steel) contact sim.AddContactModel(Hertz-Mindlin, coal, steel) contact.SetRestitution(0.3) contact.SetStaticFriction(0.5) contact.SetRollingFriction(0.01)AddMaterial返回材料对象再用它创建接触模型。第三个参数Hertz-Mindlin是EDEM默认的非线性接触模型适合无粘性颗粒如果要模拟湿颗粒换成JKR模型多一个表面能参数。SetStaticFriction对堆积角影响最大很多标定脚本会优先扫这个参数。EDEM自带了一些标准材料库但工程上最好用实验标定值替代默认值。3.3 几何边界与运动条件实际工程里边界往往不是静止壁面。螺旋输送机有转动叶片犁有直线运动。ScrewAugerTorque_Example.zip和Plough_Example.zip就是这类场景。在EDEMpy里给几何体设置运动很简单screw sim.AddGeometry(screw.stl) screw.SetMotion(rotation, speed1.5, axis[0, 0, 1]) plough sim.AddGeometry(plough.stl) plough.SetMotion(translation, velocity[1.0, 0.0, 0.0])第一段让螺旋几何绕Z轴以1.5 rad/s旋转第二段让犁沿X方向以1 m/s移动。SetMotion第一参数是运动类型speed是角速度axis是旋转轴translation用velocity指定线速度。注意几何体STL文件必须和EDEM默认的CAD坐标系一致否则会出现边界穿透。EDEMtoStl_EDEM输出stl几何文件.rar里的EDEMtoStl.py可以把EDEM内嵌几何导出成STL再做坐标变换。EDEMpy还支持虚拟边界和周期性边界。模拟螺旋输送机时用周期性边界能减少计算量设置方法是sim.SetPeriodicAxis(x, 0.2)含义是沿X方向每隔0.2米循环。这个参数不能乱设周期长度要大于最大颗粒直径的5倍否则颗粒会跨过边界自相交。Compression_Example.zip里对边界运动和时间步的配合有很细的示例压板速度过快时接触力会突变需要把时间步再缩小。4. 数据后处理与可视化EDEMpy提取结果并用matplotlib/ParaView展示4.1 读取结果文件中的颗粒数据模拟完成后EDEM会输出.dem结果文件里面按时间步保存颗粒和几何的状态。EDEMpy的Data层能直接读这些文件不需要重新打开GUI。常见的读取方式是先加载结果文件再指定时间点获取颗粒列表import edempy import numpy as np sim edempy.LoadResult(run_final.dem) for t in np.arange(0, 1.0, 0.1): particles sim.GetParticles(t) positions [p.GetPosition() for p in particles] velocities [p.GetVelocity() for p in particles]GetParticles(t)的参数是物理时间EDEMpy会自动找到离该时刻最近的结果帧并返回颗粒对象。如果发现取到的帧和请求的时间对不上需要检查EDEM输出频率设置通常用sim.SetSaveTime控制。上面的循环把每个时刻的颗粒坐标和速度存进列表后续计算分布和时序曲线都很方便。4.2 计算平均速度和偏析指数后处理不能只停留在“把数据导出”要能回答工程问题。分析颗粒流是否达到稳态时可以计算某段区域内的平均速度def avg_velocity(particles, x_min, x_max): v_sum np.zeros(3) count 0 for p in particles: pos p.GetPosition() if x_min pos[0] x_max: v_sum np.array(p.GetVelocity()) count 1 return v_sum / count if count else np.zeros(3)这个函数接受颗粒列表和空间范围返回区域平均速度向量。颗粒对象在读取后不会自行释放循环里不要保留过多历史对象的引用否则内存会涨得很快。ResidenceTime_Example.zip里其实就是一个类似逻辑按高度切片用“颗粒数量-时间曲线”拟合平均停留时间。资源里的SegregationIndex_Example.zip演示了偏析指数的计算。料堆划分为多层后统计每层小颗粒占比再计算标准差layer_frac [] segregation_index np.std(layer_frac)偏析指数越大说明颗粒分层越严重。EDEMpy在这里的价值不是帮你算而是帮你把所有层的数据按键批量取出来。常用后处理接口对比如下接口返回适用场景GetParticles(t)List[Particle]单帧查询GetParticleData(...)np.ndarray时间序列统计GetContacts(t)List[Contact]接触力分析4.3 几何导出与外部可视化有时需要在ParaView里和CFD结果叠加显示此时需要把EDEM几何导出为STL。EDEMtoStl.py在这个场景下很实用python EDEMtoStl.py -i run.dem -o geometry.stl-i指定EDEM结果文件-o指定输出的STL文件名。脚本会提取结果文件中的几何体和边界表面生成可导入模型。之后ParaView打开再叠加颗粒点云或者速度场。注意导出的STL只有三角网格没有颗粒数据颗粒仍然要用EDEMpy读取再写成CSV或VTK放入ParaView。我一般会把颗粒坐标直接导出成VTKdef write_particles_vtk(particles, path): with open(path, w) as f: f.write(# vtk DataFile Version 3.0\n) f.write(particles\nASCII\nDATASET POLYDATA\n) f.write(fPOINTS {len(particles)} double\n) for p in particles: x, y, z p.GetPosition() f.write(f{x} {y} {z}\n)这个函数把颗粒位置写成VTK点集供ParaView渲染。对于上千帧数据建议用二进制VTK或直接使用EDEM结果导入插件否则文件会非常大。EDEMpy也提供了GetParticleData接口能一次返回整个时间序列数组适合做统计而不是逐帧查询。5. 脚本自动化与案例排错从批量运行到EDEMpy常见坑5.1 批量参数扫描做敏感性分析时比如扫静摩擦系数EDEMpy的脚本化优势最明显for friction in [0.3, 0.4, 0.5]: sim.Reset() contact.SetStaticFriction(friction) sim.RunSimulation(1.0) sim.SaveResult(ffriction_{friction}.dem)这里要注意Reset()之后contact对象是否仍然有效。在部分EDEMpy版本中Reset会重建接触列表必须先重新获取接触对象再设置参数。稳妥做法是把参数设置写成独立函数每次Reset后调用。5.2 常见排错与资源说明实际使用EDEMpy时有几个错误我踩过多次错误现象可能原因检查方向启动EDEM失败可执行文件路径错误检查EDEMApplication路径导入edempy报错缺少Qt依赖安装PyQt5或sipReset后颗粒丢失颗粒工厂未绑定重新执行LoadModel运行时间过长时间步设置过大对照Rayleigh时间步导出STL空白结果帧没有几何数据确认EDEMtoStl.py参数最后说一个我常用的验证方法先用BondedParticle_Example.zip跑通标准示例再把实际模型替换进去。这个示例包含粘结颗粒的完整参数能跑通它说明EDEMpy环境和EDEM版本匹配正常。如果BondedParticle都失败优先看EDEMpy-0.1.1_ReleaseNotes.pdf中的兼容性矩阵不要盲目升级EDEMpy版本。另外跑批量时建议把EDEM的图形界面输出关掉通过sim.SetRenderMode(0)降低内存和GPU占用这样同样一台机器能同时跑两到三个算例。本文还有配套的精品资源点击获取