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

A Multi-time-step Discrete Element Method for Bar Structures

作　　者：李佟 王倩[1,2] 金先龙 LI Tong;WANG Qian;JIN Xianlong(School of Mechanical Engineering,Shanghai Jiao Tong University,Shanghai 200240,China;State Key Laboratory of Mechanical System and Vibration,Shanghai Jiao Tong University,Shanghai 200240,China)

机构地区：[1]上海交通大学机械与动力工程学院,上海200240 [2]上海交通大学机械系统与振动国家重点实验室,上海200240

出　　处：《计算物理》2022年第4期395-402,共8页Chinese Journal of Computational Physics

基　　金：国家重点研发计划(2016YFB0201800);国家自然科学基金(11772192)资助项目

摘　　要：为了提升离散元法处理连续介质问题的计算效率,提出一种基于重叠颗粒的离散元分区异步长计算方法来处理连续介质问题。该方法采用颗粒分割将连续介质划分为若干子分区,各分区采用velocity-Verlet积分格式求解运动方程。相邻分区通过重叠颗粒构成局部耦合区域,边界数据传递过程中不涉及插值和截断过程。数值算例表明,该方法在保证高精度的同时有效地降低了计算时间。在实际应用中可以有针对性地将连续介质划分为不同尺度颗粒的分区,根据问题规模及分区颗粒尺度特性采用不同时间步长,节省存储空间且大幅提升计算效率。To improve computational efficiency of discrete element method(DEM),referring to domain decomposition method and sub-cycle method in finite element method,a discrete element partition asynchronous time step computing method based on overlapping particles was proposed to deal with continuous medium problems.In the method,the continuum is divided into sub-domains.The velocity-Verlet integral scheme is used to solve the motion equation in sub-domains.Local action regions are formed with overlapping particles in adjacent sub-domains.Interpolation and truncation are not involved in transfer process of boundary data.Numerical examples show that accuracy of the method is high and the computation time is reduced effectively.In practical application,continuous media could be divided into several sub-domains with different particles size,and different time steps are adopted according to particle size characteristics of the sub-domains,with which storage space could be saved and computing efficiency improved greatly.

关 键 词：离散元方法 连接键 重叠颗粒 异步长 连续介质 

分 类 号：O341[理学—固体力学]