计算物理

• •    下一篇

面向杆系结构的离散元异步长计算方法

李佟, 王倩, 金先龙 #br#
  

  1. 1.  上海交通大学机械与动力工程学院,上海200240;2.上海交通大学机械系统与振动国家重点实验室,上海200240
  • 出版日期:2021-10-13 发布日期:2021-10-13
  • 通讯作者: 金先龙 E-mail:jxlong@sjtu.edu.cn
  • 作者简介:李佟(1994-),男,辽宁铁岭,博士,主要从事离散元数值模拟和并行计算研究 . E-mail:tong.li@sjtu.edu.cn

Multi-time-step DEM Computational Method for Bar Structures

 LI Tong, WANG Qian, JIN Xianlong #br#   

  1. 1. School of Mechanical Engineering, Shanghai Jiao Tong University, Shanghai 200240, China; 

    2. State Key Laboratory of Mechanical System and Vibration, Shanghai Jiao Tong University,Shanghai 200240, China;

  • Online:2021-10-13 Published:2021-10-13

摘要: 离散元法在处理连续介质问题的应用中越来越广泛,其计算时间与颗粒尺寸紧密相关.面对大规模应用问题,全区域采用单一的时间步长会造成计算资源的极大浪费.为了提升计算效率,提出一种基于重叠颗粒的离散元分区异步长计算方法来处理连续介质问题.该方法采用颗粒分割将连续介质划分为若干子分区,各分区采用velocity-Verlet积分格式求解运动方程.相邻分区通过重叠颗粒构成局部耦合区域,边界数据传递过程中不涉及插值和截断过程.数值算例表明该方法在保证高精度的同时有效地降低了计算时间.在实际应用中可以有针对性地将连续介质划分为不同尺度颗粒的分区,根据问题规模及分区颗粒尺度特性采用不同时间步长,节省存储空间并且大幅提升计算效率.

关键词: font-family:宋体, ">离散元方法;连接键;重叠颗粒;异步长font-family:宋体, ">;连续介质

Abstract: Discrete element method (DEM) is more and more widely used to deal with continuum problems. The calculation time is closely related to the particle size, and a single time step will cause a great waste of computing resources. In order to improve the computational efficiency, referring to the domain decomposition method and sub-cycle method in the finite element method, a discrete element partition asynchronous time step computing method based on overlapping particles was proposed to deal with the continuous medium problem. In this method, the continuum was divided into several sub-domains. The sub-domains used the velocity-Verlet integral scheme to solve the motion equation. The local action region was formed by overlapping particles in adjacent sub-domains. Interpolation and truncation were not involved in transfer process of boundary data. Numerical examples showed that the accuracy of the proposed method was high and the computation time was further reduced effectively. In practical application, continuous media could be divided into several sub-domains with different particles size, and different time steps were adopted according to the particles size characteristics of the sub-domains, by which storage space could be saved and computing efficiency greatly improved

Key words: discrete element method, linked bar, overlapping particles, asynchronous step, continuous medium