Algorithm research on surface-to-surface contact for multibody contact problem

被引：0
|
作者：
Li G. [1 ]
Liu C. [1 ]
Xiao R. [1 ]
Ma H. [2 ,3 ]
机构：
[1] Shandong Taishan Pumped Storage Power Station Co., Ltd., Tai'an
[2] State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin, China Institute of Water Resources and Hydropower Research, Beijing
[3] Earthquake Engineering Research Center, China Institute of Water Resources and Hydropower Research, Beijing
来源：
Ma, Huaifa (mahf@iwhr.com) | 1600年 / China Water Power Press卷 / 51期
关键词：
Finite element; Global search; Local search; Multibody contact; Surface-to-surface contact model;
D O I：
10.13243/j.cnki.slxb.20190926
中图分类号：
学科分类号：
摘要：
The surface-to-surface contact algorithm is focused on in this paper. (1) Combining the master-slave surface algorithm with the bit-code algorithm, a shared real constant array is defined to identify the master-surface element and the slave-surface element of the "contact-face couples". Refering to the coding ideas of the bit-code method, the potential contact surfaces in three-dimensional space are sorted as a one-dimensional array to realize global search of contact surface. (2) The combination of the point-surface algorithm and the inside-outside contact search algorithms is used to judge whether the projection point of the master node along the average normal direction is internal or external in the target surface element, and determines the local coordinates of the contact point by the area coordinates of the element, and then obtains the contact association matrix to complete local contact search, thus avoiding solving nonlinear equations of local contact point coordinates. (3) The algorithm preposed treats each contact body as a sub-region, and each sub-region can be independently divided by finite elements. The quasi-Gauss iterative method is used to implicitly solve the displacement increment and contact force. The FORTRAN source code program is compiled in accordance with the proposed contact algorithm strategy. (4) Calculation results obtained from the numerical examples are almost completely matched the theoretical solution of the classic contact problem. At the same time, the contact algorithm presented has been verified by using similar contact calculation functions provided by ANSYS commercial software because their calculation results are almost identical. Another advantage of the algorithm proposed in this paper that the contact area decomposition has been taken into account during the pre-processing of contact search, which can conveniently perform parallel computional processing for high concrete dams and high slope stability of rock masses. © 2020, China Water Power Press. All right reserved.
引用
收藏
页码：597 / 605
页数：8
相关论文
共 15 条
  • [1] HALLQUIST J O, GOUDREAU G L, BENSON D J., Sliding interfaces with contact-impact in large scale lagrangian computation, Computer Methods in Applied Mechanics and Engineering, 5, 1, pp. 107-137, (1985)
  • [2] ZHONG Z H, NILSSON L., A unified contact algorithm based on territory concept, Computer Methods in Applied Mechanics and Engineering, 130, pp. l-16, (1996)
  • [3] ZHONG Z H, NILSSON L., Automatic contact searching algorithm of dynamic finite element analysis, Computers and Structures, 52, 2, pp. 187-197, (1994)
  • [4] BENSON D J, HALLQUIST J O., A single surface contact algorithm for the post-buckling analysis of shell structures, Computer Methods in Applied Mechanics and Engineering, 78, 2, pp. 141-163, (1990)
  • [5] MATS OLDENBURG, MLSSON LARSGUNNAR, The position code algorithm for contact searching, Int. J. for Numerical Methods in Engineering, 37, pp. 359-386, (1994)
  • [6] HU CHEN, ZHOU LEI, MENGYAN ZANG, LC-Grid: a linear global contact search algorithm for finite element analysis, Comput. Mech, 54, pp. 1285-1301, (2014)
  • [7] TED BELYTSCHKO, NEAL MARK O, Contact-impact by the pinball algorithm with penalty and Lagrangian methods, Int. J. for Numerical Methods in Engineering, 31, pp. 547-572, (1991)
  • [8] HALLQUIST J O, WAINSCOTT B, SCHWEIZERHOF K., Improved simulation of thin-sheet metal forming using LS-DYNA3D on parallel computers, Journal of Materials Processing Technology, 51, pp. 144-157, (1995)
  • [9] ULAGA S, ULBIN M, FLASKER J., Contact problems of gears using Overhauser splines, Int. J. of Mechanical Sciences, 41, pp. 385-395, (1999)
  • [10] WANG PING, EIIJI NAKAMACHI, The inside-outside contact search algorithm for finite element analysis, Int. J. for Numerical Methods in Engineering, 40, pp. 3665-3685, (1997)
12