三维离散元法及其在脆性材料破坏问题中的应用.pdfVIP

第十届全国振动理论及应用学术会议论文集(2011) 三维离散元法及其在脆性材料破坏问题中的应用 成名 （上海飞机设计研究院，上海200232 ） 摘要：离散元法的三维模型较于其二维模型，具有更强的适用性和更广阔的工程应用前景。本文通过能量 等效准则推导出二十七球模型的弹簧系数，建立三维离散元法，并将其应用于脆性材料破坏问题。通过模 拟数值算例，并分别与相应实验的结果进行比较，验证了本文的三维离散元法在解决脆性材料破坏问题中 的适用性和准确性。 关键词： 离散元法；三维问题；脆性材料；破坏 中图分类号：O389 Three-Dimensional Discrete Element Method and its Application in Modeling Severe Damage Problems of Brittle Materials Cheng Ming (Shanghai Aircraft Design and Research Institute, Shanghai 200232, P.R. China) Abstract: In this paper, a scheme of the discrete element method (DEM) for numerical simulation of three-dimensional (3-D) severe damage problems of Brittle Materials is proposed. Using the principle of energy equivalence, the three-dimensional DEM models is developed. The model is also validated by numerical examples, and the numerical results have good agreement with the results obtained by the corresponding experiments. Key words: Discrete element method; Three-dimensional problem; Brittle Materials; Severe Damage 1 引言 脆性材料的破坏问题是一个复杂的力学过程，其实质是计算模型由连续介质过渡到非连 续介质的动态转化过程。目前的计算方法对于这样复杂的过程尚不能提供很好的解决方案。 譬如有限差分法很难追踪移动边界和交界面，并难以处理复杂的、不规则的几何形状；而有 限元法在处理大变形问题时会出现严重的网格畸变，并在处理破坏问题时存在单元侵蚀、质 量丢失的现象。理想的解决方案是使某种数值方法能够应用在包括连续介质、非连续介质和 两者相转化过程的统一计算中。具有散体计算传统优势的离散元法，通过建立其在连续介质 力学中的理论及算法，就能在上述问题中展露出优越性和良好的发展前景。 离散元法在岩土工程和粉体工程中已经得到了广泛的应用[1] 。但是传统的离散元法计算 连续介质问题时精度远远不如有限元法，连续介质离散方案缺乏严密的理论基础，这些都严 重制约了离散元法的应用。Yiwai 、Morii 提出了弹簧圆盘离散元模型[2][3]，刘凯欣等在此基 [4][5][6] [7] 础上建立了弹性圆盘离散元模型 和弹塑性圆盘离散元模型 的理论，并验证了它们在 计算连续介质力学问题方面的计算精度及效率。 本文应用三维离散元模型并引入 Mohr-Coulomb 型破坏准则和与应变率有关的材料破 坏系数来模拟脆性材料的破坏过程，通过将计算结果与文献进行对比，证明了该算法的有效 性和可靠性。 2 三维离散元法的基本方程 1158 第十六部分 大飞机中的动力学问题