中国机械工程

• 先进材料加工工程 • 上一篇    下一篇

弹性力学中无网格和有限元耦合的元胞自动机算法

陈泽芸;袁卫锋   

  1. 西南科技大学制造科学与工程学院,绵阳,621010
  • 出版日期:2017-09-10 发布日期:2017-09-07
  • 基金资助:
    国家自然科学基金资助项目(11472232);
    四川省教育厅教改项目(14JGCX07);
    制造过程测试技术省部共建实验室基金资助项目(13zxzk05)
    National Natural Science Foundation of China (No. 11472232)

A Coupling CA Algorithm of Meshless and FEM in Elasticity

CHEN Zeyun;YUAN Weifeng   

  1. School of Manufacturing Science and Engineering,Southwest University of Science and Technology,Mianyang,Sichuan,621010
  • Online:2017-09-10 Published:2017-09-07
  • Supported by:
    National Natural Science Foundation of China (No. 11472232)

摘要: 结合有限元和无网格算法的优势,提出了一种元胞自动机算法用以求解二维弹性力学问题。该算法将二维模型离散成一系列节点,这些节点被分成有限元群和无网格群。有限元区域被定义在问题的边界附近,其中的任一节点和其周围相邻点的力学关系通过有限元单元建立;无网格区域定义在远离原理问题边界处,其中的节点之间的关系借用有限元中的位移插值概念建立。无论处于有限元区域还是无网格区域,任何一个节点都被置于元胞自动机的框架下进行处理,即节点的位移通过元胞自动机进行求解。与有限元方法相比,所提出的元胞自动机算法无需采用高斯消去法等传统系统求解器,而是通过元胞自动机的自动演化解决问题。依据该算法,有限元和无网格方法可以实现无缝连接。数值算例验证了该算法的新颖性和正确性。

关键词: 耦合算法, 元胞自动机, 有限元方法, 弹性力学

Abstract: A CA was developed to combine the advantages of FEM and meshless method which was used to analyze 2-dimensional elastic problems. By the CA, a 2-dimensional domain was discretized into a grid of nodes. These nodes were distributed randomly in the domain and they were divided into the FEM group and the meshless group. In the FEM zone which was mostly defined near the boundaries of the domain, FEM elements were used to establish the relationship between one node and other adjacent nodes. In the meshless zone which was normally away from the boundaries, the adjacent nodes were related by employing the concept of displacement interpolation used in FEM. However, each nodes, wherever they were in the FEM zone or meshless zone, were dealt with under the CA frame, that is, the unknown displacements of each nodes were evaluated by CA algorithm. Unlike FEM, the proposed CA solver worked out the results through automatic evolution, instead of Gaussian elimination or other conventional methods. Based on the current CA algorithm, FEM and meshless may be merged seamlessly. The novelty and the correctness of the current approach were proven by numerical examples.

Key words: coupling algorithm, cellular automaton(CA), finite element method(FEM), elasticity

中图分类号: