欢迎光临《计算力学学报》官方网站!
李昊,蒋明镜,阿比尔的,王多银.离散元中反力型刚性边界的高效模拟方法[J].计算力学学报,2025,42(4):607~614
本文二维码信息
码上扫一扫!
离散元中反力型刚性边界的高效模拟方法
An efficient simulation method of force-reactive rigid boundary in DEM
投稿时间:2024-02-05  修订日期:2024-06-10
DOI:10.7511/jslx20240205002
中文关键词:  离散单元法(DEM)  反力型刚性边界  边界条件  计算方法
英文关键词:discrete element method(DEM)  force-reactive rigid boundary  boundary condition  calculation method
基金项目:国家自然科学基金(51890911;52331010);国家重点研发计划(2022YFC3003403);重庆交通大学研究生科研创新项目(2022S0029)资助.
作者单位E-mail
李昊 重庆交通大学 河海学院, 重庆 400074  
蒋明镜 重庆交通大学 河海学院, 重庆 400074
苏州科技大学 土木工程学院, 苏州 215011
同济大学 土木工程学院, 上海 200092 
mingjing.jiang@usts.edu.cn,mingjing.jiang@tongji.edu.cn 
阿比尔的 重庆交通大学 河海学院, 重庆 400074  
王多银 重庆交通大学 河海学院, 重庆 400074  
摘要点击次数: 783
全文下载次数: 818
中文摘要:
      使用离散单元法(DEM)进行数值仿真分析时,对颗粒体施加精确的边界条件是保证模拟结果准确性的必要前提之一。在DEM中使用刚性边界对颗粒体施加应力边界条件时,若边界处的应力条件随边界位置的变化而变化(称该类边界为反力型边界),则采用现有的伺服算法对反力型刚性边界进行模拟时会存在难以收敛的问题。为避免这一不足,本文提出了一种可在DEM仿真中精确高效地计算弹性范围内轴对称分布反力型刚性边界相对平衡位置的迭代计算方法,并将该算法植入到DEM商业软件PFC3D v6.0中,验证模拟算法的可行性及程序编写的正确性;最后进一步将算法程序用于河海大学提出的一种新型静止侧压力系数的DEM仿真试验研究中。结果表明,平衡位置算法能够准确高效地模拟DEM仿真过程中任意时刻反力型刚性边界的相对平衡位置;使用该算法后,所得仿真试验结果的波动明显减弱并获得了实际分析价值,且计算过程中迭代失败概率很小,计算效率较高。因此,可认为本文提出的平衡位置算法可在DEM中对反力型刚性边界进行精确高效的模拟。
英文摘要:
      An accurate boundary condition should be used in the discrete element method(DEM) simulation to ensure the accuracy of results.When the rigid boundary is used to apply the stress boundary condition to the particle assembly,the stress condition of the boundary may be changed with the boundary position(this kind of boundary is called force-reactive boundary in this research).The servo method,which is used in current DEM simulation to control the boundary position,is more difficult to achieve the convergence condition when it is used to simulate the force-reactive rigid boundary.To avoid this problem,an efficient method was proposed to calculate the relative equilibrium position of the axisymmetric force-reactive rigid boundary accurately in elastic range with the iterative method.Then the method was developed into a commercial DEM software package called PFC3D v6.0,and was tested by a simple example to verify its feasibility and correctness.Finally,the verified program was used to the DEM simulation of a new apparatus for the earth pressure coefficient static test proposed by Hohai University.The results show that the efficient method can accurately calculate the relative equilibrium position of the force-reactive rigid boundary during the simulation.It is also found the calculation method proposed in this paper generates results of less fluctuation,and the iterative calculation still maintains high efficiency and low failure rate.Accordingly,the method proposed in this paper can accurately and efficiently simulate the force-reactive rigid boundary in the DEM simulation.
查看全文  查看/发表评论  下载PDF阅读器
您是第17744758位访问者
版权所有:《计算力学学报》编辑部
本系统由 北京勤云科技发展有限公司设计