基于三角剖分局部搜索的三维滑移面算法研究.pdfVIP

第四届全国计算爆炸力学会议论文集 * 基于三角剖分局部搜索的三维滑移面算法 姜玉曦 周海兵 张树道 （北京应用物理与计算数学研究所，北京 100088 ） 摘 要：武器动力学过程的拉氏模拟要求能够处理复杂接触碰撞界面的高效的接触算法，而确定接 触点的局部接触搜索算法是接触算法的一个重要部分。CHAP3D 是应用相容拉氏方法开发的流体动 力学程序，其主要特点是在离散方法中引入形心，将单元剖分成若干个四面体角单元，然后应用离 散相容微分算子保证流体运动微分控制方程的总能守恒。因此提出了三角剖分局部搜索的算法，通 过接触面单元四个节点的位置平均引入形心，将一个四边形的接触曲面单元剖分为四个法向确定的 三角形接触面子单元，由此求得接触点。将基于此局部搜索方法的滑移面算法应用到 CHAP3D 中， 数值结果显示此方法与相容拉氏方法的离散特点相符并且能够克服传统方法由于网格扭曲变形带来 的接触点求解困难。 关键词：局部接触搜索；接触算法；滑移面；拉氏算法 A Three-dimensional Contact Algorithm for Sliding Surface Based on Triangular Subdivision Local Search Jiang Yuxi Zhou Haibing Zhang shudao ( Beijing Institute of Applied Physics and Computational Mathematics, Beijing, 100088 ) Abstract: The numerical simulation for the weapon dynamics requests an efficient contact algorithm which can deal with the complex contact-impact interface. One of the most important parts of the contact algorithm is the local contact search. CHAP3D is a hydrodynamic code developed by using compatible Lagrangian numerical method. This method divides the discretized cell into some tetrahedrons by using the shape hearts of the cell and the faces of the cell, and preserves total energy conservation by compatible differential operator discretization. So a triangular subdivision local contact search algorithm has been developed. A shape heart is introduced by averaging the positions of the four contact segment nodes in this method. Then a four-node quadrilateral contact surface segment may be subdivided into four triangular sub-segments which normal vectors are confirmed with the shape heart point. Therefore the contact point is achieved by the geometry method. The contact algorithm