输出

PD3D元素没有可用的元素输出。节点输出包括Abaqus/Explicit分析中通?？捎玫乃惺涑霰淞?见Abaqus/Explicit Output Variable ldentifier)。

局限性

离散元法分析会受到以下限制：

l应力、应变和其他类似连续体单元输出的体积平均输出不适用于DEM分析。

lPD3D元素只支持球形形状。

lPD3D元素不能成为刚体定义的一部分。

l带有使用PD3D元素和BEAM类型多点约束的刚性簇的模型不能在域并行式下运行。

l不可能指定PD3D元素之间或PD3D元素和其他元素之间的粘结或热接触。

lPD3D单元之间或PD3D单元与其他单元之间的接触忽略滚动摩擦。

lPD3D元素之间的接触不支持用户定义的表面交互。

l虽然在Abaqus/iewer中支持此功能，但在Abaqus/CAE中不支持此功能。您可以使用Abaqus/CAE中的现有功能来生成质量元素，编写输入文件，然后手动编辑输入文件以将质量元素转换为粒子?；蛘?，您可以使用C3D8R元素创建网格，编写输入文件，然后使用脚本将这些元素转换为粒子。将实体元素转换为SPH粒子的Pvthon脚本示例描述在达索系统知识库的“从实体网格生成SPH粒子元素”中。http：/support.3ds.com/knowledge-base/.

DEM计算分布在并行域，除非多个离散部分定义不同的a阻尼参数(这将降低并行可扩展性)。如果使用多个CPU，DEM分析会受到以下限制：

lDEM辅助节点不支持触点输出。

l不支持除整个模型之外的能量历史输出。

l无法激活动态负载平衡。

l如果任何DEM粒子参与一般接触，所有DEM粒子必须包含在一般接触定义中。

l建议每个域至少10，000个DEM粒子，以达到良好的可扩展性。

l如果使用大量的CPU，可能需要大幅增加内存使用量。

输入文件模板

下面的例子说明了一个离散元法分析：

*HEADING

…

*ELEMENT, TYPE=PD3D, ELSET=name

Element number, node number

…

*DISCRETE SECTION, ELSET=name

Particle radius

**

*INITIAL CONDITIONS, TYPE=VELOCITY

Data lines to define velocity initial conditions

*NSET, NSET=name, ELSET=name

*SURFACE, NAME=name

,

**

*SURFACE INTERACTION

**

*STEP

*DYNAMIC, EXPLICIT

*DLOAD

Data lines to define gravity load

**

*CONTACT

*CONTACT INCLUSIONS

*CONTACT PROPERTY ASSIGNMENT

**

*CONTACT CONTROLS ASSIGNMENT

*OUTPUT, FIELD

*END STEP

参考资料

l“颗粒组合体的离散单元法”，地质工程，第29卷，第47-65页，1979年。

lJohnson,K.L.，K.Kendall，和A.D.Roberts,“表面能与弹性固体的接触”，“伦敦皇家学会议事录”，第324卷，第301-313页，1971年。

l“相似尺寸物体的接触检测算法”，国立成功大学工程数值方法研究所硕士论文，第43卷，第131-149页，1998。

l0'Sullivan，C.和J.D.Bray,“为使用中央差分时间积分方案的离散单元模拟选择适当的时间步长”，“工程计算”，第21卷(2/3/4)，第278至303页，2004年。

l“微粒系统离散微粒模拟之研究进展”，国立成功大学化学工程研究所硕士论文，2006。。朱海鹏，

l朱海鹏，“微粒系统离散微粒模拟之研究””，国立成功大学化学工程研究所硕士论文，2007。

温馨提示：

此文档为达索官方英文文档翻译，尽管我们已经尽力确保准确性，但在翻译过程中可能会有一些错误或细微差别。如果想要了解官方原版，可联系客服进行索取。