›› 2014, Vol. 35 ›› Issue (9): 2717-2724.

• 数值分析 • 上一篇    下一篇

基于OpenMP的二维有限元-离散元并行分析方法

严成增,郑 宏,孙冠华,葛修润   

  1. 中国科学院武汉岩土力学研究所 岩土力学与工程国家重点实验室,武汉 430071
  • 收稿日期:2013-05-25 出版日期:2014-09-10 发布日期:2014-09-16
  • 作者简介:严成增,男,1986年生,博士研究生,主要从事离散元有限元耦合、离散元、颗粒物质物理与力学、多物理场耦合、DDA、大规模高性能并行计算以及软件工程等方面的研究工作。
  • 基金资助:

    国家自然科学基金项目(No. 11202223, No. 11172313);国家重点基础研究发展计划(973)项目(No. 2011CB013505)。

Parallel analysis of two-dimensional finite-discrete element method based on OpenMP

YAN Cheng-zeng, ZHENG Hong, SUN Guan-hua, GE Xiu-run   

  1. State Key Laboratory of Geomechanics and Geotechnical Engineering, Institute of Rock and Soil Mechanics, Chinese Academy of Sciences,Wuhan 430071, China
  • Received:2013-05-25 Online:2014-09-10 Published:2014-09-16

摘要: Munjiza提出的有限元-离散元耦合分析方法(FEDM)是分析岩石破裂过程的一种十分有效的方法。然而,为了克服网格依赖性,需要将岩体剖分成非常细小的三角形单元,且三角形单元之间不共用节点,导致问题的变量数目巨大,计算非常耗时。为了提高计算效率,基于OpenMP(open multi-processing)多核并行技术实现了有限元-离散元法的并行化,克服了并行化过程中存在的数据竞争,实现了并行程序的负载平衡。提出了一套将串行程序并行化的策略,即首先确定串行程序的热点区域,然后尽可能地将热点区域并行化,尽量使用私有变量来规避数据竞争;若各线程间仍存在数据竞争,可采用动态链表数据结构,先将数据存于动态链表中,最后在并行区域外,将存于各个动态链表中的数据进行合并,这样可以规避数据竞争,同时避免了使用临界区或锁,从而提高了程序的并行化效果。开发了并行版本FDEM程序,将所提出的方法用于大规模工程问题的求解,最后通过陡崖塌落的算例说明该方法的有效性。

关键词: 有限元-离散元法, 并行计算, 负载均衡

Abstract: The coupling method of finite- discrete element method(FDEM) by Munjiza is an effective procedure for simulating rock failure. In order to overcome the mesh dependency, it is necessary to partition the rock mass into very small triangular elements; and the elements do not share the nodes, creating a vast number of variables, therefore, the calculation is very time consuming. The method overcomes competition among the existing data in the parallelization and balances loads in a parallel program. A set of strategies for parallelizing serial programs are proposed, determining the hot zones of the serial program, then parallelizing the hot zones as far as possible and using as many private variable as possible to circumvent competition. If there are still data competition among the treads, dynamic lists can be used to store these competitive data. Finally, the data stored in the various dynamic lists are merged outside the parallel zones. Data competition is accordingly avoided with no use of critical section or lock; and the effect of the program parallelization is hence increased. By the proposed method, the parallel version of the FDEM program is developed, which can be used to solve large-scale problem. A numerical example about cliffs collapse shows the effectiveness of the proposed method.

Key words: finite-discrete element method(FDEM), parallel computing, load balancing

中图分类号: 

  • O 241
[1] 张友良,谭 飞,张礼仁,施明明. 岩土工程亿级单元有限元模型可扩展并行计算[J]. , 2016, 37(11): 3309-3316.
[2] 严成增 ,郑 宏 ,孙冠华 ,葛修润,. 基于FDEM-Flow研究地应力对水力压裂的影响[J]. , 2016, 37(1): 237-246.
[3] 严成增,郑 宏,葛修润. 基于统一标定的势接触力计算[J]. , 2015, 36(1): 249-256.
[4] 付晓东,盛 谦,张勇慧. 基于OpenMP的非连续变形分析并行计算方法[J]. , 2014, 35(8): 2401-2407.
[5] 杨 勋 ,王欢欢 ,余克勤 ,金先龙 , . 行波激励下防波堤地震动力响应分析[J]. , 2014, 35(6): 1775-1781.
[6] 刘耀儒,杨 强,薛利军,周维垣. 基于三维非线性有限元的边坡稳定分析方法[J]. , 2007, 28(9): 1894-1898.
[7] 付朝江 ,张 武 . 基于粗细网格的有限元并行分析方法[J]. , 2006, 27(5): 807-810.
[8] 李邵军; 冯夏庭; 安红刚;. 大型洞室群软岩置换方案优化的并行实现[J]. , 2003, 24(3): 335-338.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 徐金明,羌培,张鹏飞. 粉质黏土图像的纹理特征分析[J]. , 2009, 30(10): 2903 -2907 .
[2] 梁桂兰,徐卫亚,谈小龙. 基于熵权的可拓理论在岩体质量评价中的应用[J]. , 2010, 31(2): 535 -540 .
[3] 马文涛. 基于灰色最小二乘支持向量机的边坡位移预测[J]. , 2010, 31(5): 1670 -1674 .
[4] 于琳琳,徐学燕,邱明国,闫自利,李鹏飞. 冻融作用对饱和粉质黏土抗剪性能的影响[J]. , 2010, 31(8): 2448 -2452 .
[5] 王 伟,刘必灯,周正华,王玉石,赵纪生. 刚度和阻尼频率相关的等效线性化方法[J]. , 2010, 31(12): 3928 -3933 .
[6] 王海波,徐 明,宋二祥. 基于硬化土模型的小应变本构模型研究[J]. , 2011, 32(1): 39 -43 .
[7] 曹光栩,宋二祥,徐 明. 山区机场高填方地基工后沉降变形简化算法[J]. , 2011, 32(S1): 1 -5 .
[8] 刘华丽 ,朱大勇 ,钱七虎 ,李宏伟. 边坡三维端部效应分析[J]. , 2011, 32(6): 1905 -1909 .
[9] 刘年平 ,王宏图 ,袁志刚 ,刘竟成. 砂土液化预测的Fisher判别模型及应用[J]. , 2012, 33(2): 554 -557 .
[10] 王卫东 ,李永辉 ,吴江斌 . 超长灌注桩桩-土界面剪切模型及其有限元模拟[J]. , 2012, 33(12): 3818 -3824 .