matlab力学仿真_多体仿真模拟颗粒状产品流动仿真

多体仿真(MBS)适用于机械系统精确的运动学和动力学性能分析。如果外部作用力对整个系统性能有显著影响，仅考虑机械系统是不够的。MBS的目标是建立一个考虑所有相关力的机械系统模型。

农业机械公司格立莫(Grimme)研发并制造马铃薯和甜菜种植机械。在本案例中，马铃薯颗粒在机器内流动的性能和载荷是非常重要的。对这种问题的一种常见方法是使用商业“离散元方法(DEM)”软件和Simpack进行联合仿真。为了避免额外的软件许可证成本和接口开发，在本仿真场景的范围内，研究了仅用Simpack仿真整个模型的可能性。使用“多边形接触方法(PCM)”定义所有物体之间的接触特性。

1

前言

以马铃薯播种机为例进行了可行性研究。为了建立这个模型，有必要建立一个与马铃薯颗粒流相互作用的橡胶弹性皮带，且皮带由电动振动装置进行激励。带的振动使马铃薯种子被分离到各自的勺杯中。

使用Simpack这个功能强大的MBS(多体系统)软件可以很容易实现机械部件的建模。主要问题是，在Simpack中是否能实现使用接触力元能够充分地表示许多类似颗粒体之间的接触，其中一些还具有复杂的几何结构。

另一个挑战是，是否能通过使用适当的参数以及良好的求解器设置，将计算时间保持在可接受的范围内，尽管模型有很多的自由度。

图 1: 土豆播种机的多体仿真模型

2

常规模型设置

高弹性皮带的建模是比较困难的挑战之一。通过把皮带离散为多个刚体段来实现(见图2)。大量的离散段可以更好地表示皮带的特性，并减少由多边形效应引起的带轮的振动激励。不过必须在离散段数和计算时间之间找到适当的平衡。各离散段之间通过衬套力元耦合。为了后续定义皮带段和带轮之间以及皮带段和土豆颗粒之间的接触力，有必要使用长方体和圆柱体对这些离散段进行建模。

所有和土豆颗粒有接触的机械部件也必须在模型中建立。在本案例中，机械部件指的是拾取室和向下运动的通道。这些部件几何图形从CAD文件中获取并转换为obj格式文件。然后，机械部件和颗粒体之间的接触力可以使用PCM接触方法精确建模。

图 2: 皮带建模示意图

3

接触定义

本模型中使用了两种不同的接触类型。第一种是曲线-曲线接触，用于定义皮带段和带轮之间的接触。与其它几何接触类型相比，这种接触类型的计算时间能明显缩短。由于曲线-曲线接触只能在具有曲面的接触部件之间实现，因此每个皮带段必须含有上述的小圆柱体曲面。

此外，必须定义颗粒之间以及颗粒与机器部件之间的接触，这使用Simpack软件PCM方法实现。任何几何都可用于定义接触面，唯一的要求是接触面必须是三角形网格形式的文件格式，如obj等。PCM接触允许接触面之间有多个接触，这是计算颗粒接触的先决条件。这种接触方法的另一个优点是高效的碰撞检测，因为不接触的表面对计算时间的影响很小。

这个播种机模型中接触力元的总数取决于颗粒的数量和皮带的离散程度。该模型含有108个皮带段和24个土豆颗粒，共有4309个接触力元。

4

使用Matlab生成模型数据

由于所需的力元数量特别多，如果使用标准的Simpack GUI操作，生成模型将非常耗时。因此，为了自动生成Simpack的模型文件，开发了一个Matlab脚本。

所有的模型参数，如皮带速度和土豆数量，都是在Matlab中预定义的。随后，执行该脚本生成用于仿真必要的数据文件，这些文件之后可以直接运行。对模型的修改总是需要修改Matlab脚本。在创建新的仿真模型时，通常可以采用脚本的部分内容，这能大大减少了建模所需的时间。

5

仿真场景

SODASRT 2积分器适用于所有的仿真过程。事实证明，将最大积分器步长限制在10-4秒是有效的，这样可以可靠地检测所有接触事件，并避免对积分器步长的较大更改。

首先运行仿真用于调整模型中的接触参数。通过各种实际的跌落冲击试验，确定了马铃薯与各种部件材料之间的实际接触特性。这些试验确定冲击力过程和恢复系数，与马铃薯的质量、跌落高度和接触材料等因素有关。然后，仿真研究的主要目的是优化计算时间，以实现尽可能多的颗粒数量仿真。

第二步，通过分析机械参数和部件几何参数的变化，优化颗粒流运动过程，降低土豆颗粒的受力。

6

结果

在一个简单的跌落试验中，对比仿真和实测的冲击力或冲击加速度曲线，显示具有良好的一致性。图3显示了当一个100克的土豆从50厘米的高度落在钢表面上时的冲击变化曲线。

在后续的产品开发模型分析中，作用于土豆颗粒的冲击力对Grimme公司来说是非常重要的。所有新研发的一个关键目标是最大程度地保护土豆。

图 3: 试验(蓝色)和仿真(红色)的冲击加速度曲线对比

仿真模型的质量进一步评估标准是直接观察对比颗粒流的运动情况。在试验和仿真过程中，除了马铃薯外，还使用了高尔夫球，因为高尔夫球具有更好的可比性和均匀的颗粒特性。通过台架试验获得的高速图像作为土豆颗粒的参考。

图4显示了台架试验高速成像图，左侧显示了仿真的类似状态。通过比较这个和其它视频，可以确定仿真模型中颗粒非常真实的运动性能。值得注意的是，即使是接触参数的微小变化也会对产品运行产生显著影响。

图 4: 仿真和试验的可视化对比

7

总结

以播种机为例，本文论证了Simpack在颗粒模型建模和仿真中的适用性。但是在模型结构、颗粒数量和计算时间方面存在一定的限制。

仿真结果与试验测试值具有很高的一致性。同时考虑了作用在颗粒上冲击力随时间变化以及整个颗粒流的运动性能。对于一些仿真模型，特别是那些颗粒较少的模型，可以不用联合仿真，仅在Simpack中就能对整个系统进行建模和仿真。

仿真动画请点击下面的视频观看。

扩展阅读：

基于柔性梁的受电弓/接触网动力学分析

齿轮箱坎贝尔图分析

使用柔性齿轮检测齿轮损坏的根本原因

为什么说齿轮仿真需要Simpack

Simpack Automotive 专业的汽车动力学仿真平台