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