大尺寸分布颗粒离散元数值计算中接触检测算法的改进研究.pdfVIP

大尺寸分布颗粒离散元数值计算中接触检测算法的 改进研究 1 1 2 黄昕 ，刘义伦 ，李方 1.中南大学机电工程学院，湖南长沙（410083 ）； 2. 比亚迪股份有限公司第三事业部研究部，广东深圳（518116 ） 摘 要：针对大尺寸分布颗粒离散元数值计算中传统接触检测算法的缺点，提出了改进网 格法，分析了原算法和改进算法的原理，比较了新旧算法在效率和存储空间上的差异，并 通过应用于炭素糊料压实过程的模拟，验证了新算法的优越性。改进方法提高了大尺寸分 布颗粒离散元数值计算中接触检测算法的效率。 关键词：离散元法；接触检测；邻居搜索；网格法；大尺寸分布 1 引言 颗粒体在实际生产和试验中表现出的许多复杂行为，通常无法直接使用现有基本理论， 尤其是连续介质理论来解释。离散元法把所分析的系统当作离散颗粒的集合体，采用简单 方程即可对高度复杂系统的准静态和动态行为进行模拟分析，并能够直接获得系统内部的 复杂行为信息，从而进一步假定和分析离散物质的细观结构，为准确预测和分析现有连续 介质理论无法解释和分析的力学行为提供了理论依据和研究方法。在岩土[1]、农业[2]、化 学[3] [4]等工程行业得到广泛应用，并逐渐从二维扩展到三维。 离散元的计算过程主要包括三个方面：(1).单元间以及单元与边界的相对位置关系判断 及计算，通常称为接触检测；(2).单元运动判断及计算；(3).其它等效物理场如应力、应变 等的计算。其中第一步的计算量最大，耗时最多。接触检测最早采用直接搜索法，算法的 时间复杂度为O(N2) ，计算时间随颗粒单元N 的增大呈几何级数增长，在实际工程问题中 无法应用。为了减小计算量，一般将接触检测过程分两步来完成：(1).邻居搜索，即找出在 给定时间内与系统中每个单元可能发生接触的其它单元 (a) 或边界；(2).接触搜索，即根据颗粒的邻居数据确定该单 元与邻居单元或边界是否真的存在接触，如果存在接触 则进行相互作用计算。由于颗粒离散元通常将颗粒简化 为球形，因此接触搜索过程相对简单，邻居搜索过程计 算量较大。应用于离散元法数值计算的邻居搜索算法很 多，主要有邻居链表法[5] [6]、网格法(Cell Method)[6] [7] (b) [8]、边界盒法(Bounding Box)[9] [10]、四叉树(Quad trees) 和八叉树(Oct trees)[11] 、三角剖分法(Dynamic Triangulations)[12]等。其中邻居链表法的平均计算复杂 度比直接搜索法要略小，但仍为 O(N2)[3]，只适合小规 模的计算。边界盒法效率一般但所需存储空间较大。四 （八）叉树和三角剖分法实现起来比较复杂，实际应用 (a) 二维 (b) 三维 较少。网格法的效率较高且易于实现，因而应用最为广 图1 需要进行接触判断的邻居格子 作者简介：黄昕(1978-)，男,湖北大悟人，博士研究生，主要研究方向：颗粒物质压实及其数值计算。 -1- 泛。 2 网格法 网格法将分析区域划分为规则的网格，网格坐标固定不动，与颗粒单元的运动无关。 如果颗粒根据其当前位置临时被分配到指定的网格中，则只可能与该网格及其相邻网格中 的颗粒发生相互作用。对于2D/3D 系统只有9/27 个不同网格中的颗粒可能是邻居。由于颗 粒对只需检测一次，因此并非所有相邻网格中的颗粒都需要检测，而只需检测当前