摘要
大步长静力隐式有限元方法具有快速的计算速度和较高的计算精度,在金属板料的成形模拟领域中得到了广泛的应用.在接触搜寻阶段,由大步长引起的大量历史接触信息改变的特点使接触搜寻问题变得更困难,因此有必要针对该问题设计一种高效稳健的搜寻算法.基于主从面法的思想,提出了一种结合位置码法和ADT(Alternating Ditigal Tree,交替数字二叉树)法的改进接触搜寻算法,并集成到自主研发的QuickForm软件中;通过数值算例对该方法的计算性能进行了测试.结果表明:改进算法具有良好的搜寻性能,平均搜寻效率提升了20%左右.
The large step implicit FEMs have been widely used in the field of sheet metal forming due to its fast calculation speed and highly accuracy.Many historical contact information may be changed for the current time step in implicit FEMs which makes the contact searching problem more difficult,and a robust contact searching algorithm is necessary.Combining location code algorithm and Alternating Ditigal Tree(ADT) structure,an algorithm based on master-slave method was proposed for large step implicit FEMs,and it was integrated into the QuickForm software.The experimental results show that the algorithm has good search performance.For many complex forming parts,the average search efficiency is up to about 20%.
引文
[1] 谭继锦.汽车有限元法 [M].2版.北京:人民交通出版社,2012.
[2] 陈吉清,刘朝阳,兰凤崇.中国人体上肢碰撞损伤有限元模型的开发 [J].华南理工大学学报 (自然科学版),2017,45(8):21- 27.CHEN Ji-qing,LIU Chao-yang,LAN Feng-chong.Develop-ment of finite element model of impact injury to upper limb of Chinese human body [J].Journal of South China University of Technology(Natural Science Edition),2017,45(8):21- 27.
[3] 马晓春,沈卫兵.有限元数值模拟技术在汽车冲压件成形中的应用 [J].浙江工业大学学报,2007,35(1):100- 104.MA Xiao-chun,SHEN Wei-bing.Application of numerical simulation of finite element in molding process of automobile pressing parts [J].Journal of Zhejiang University of Technology,2007,35(1):100- 104.
[4] HALLQUIST J O.LS-DYNAtheoretical manual [EB/OL].(2019- 02- 20) [2019- 04- 18]http://ftp.lstc.com/anonymous/outgoing/jday/manuals/DRAFT_Theory.pdf
[5] BENSON D J,HALLQUIST J O.A single surface contact algorithm for the post-buckling analysis of shell structures [J].Computer Methods in Applied Mechanics & Engineering,1990,78(2):141- 163.
[6] CHEN H,LEI Z,ZANG M.LC-Grid:a linear global contact search algorithm for finite element analysis [J].Computational Mechanics,2014,54(5):1285- 1301.
[7] OLDENBURG M,NILSSON L.The position code algorithm for contact searching [J].International Journal for Numerical Methods in Engineering,1994,37(3):359- 386.
[8] 张鹏.板料冲压成形仿真接触搜索算法改进及程序实现 [D].大连:大连理工大学,2014.
[9] 刘泉声,何军,吴月秀.数值流形方法接触检索算法的改进 [J].岩石力学与工程学报,2016,35(1):40- 49.LIU Quan-sheng,HE Jun,WU Yue-xiu.An improvement of contact detection algorithm of numerical manifold me-thod [J].Chinese Journal of Rock Mechanics and Engineering,2016,35(1):40- 49.
[10] BRUNEEL H C J,DE-RYCKE I.QuickTrace:a fast algorithm to detect contact [J].International Journal for Numerical Methods in Engineering,2010,54(2):299- 316.
[11] 陈成军,柳明,陈小伟,等.基于八叉树的全局接触搜索算法研究 [J].计算力学学报,2017,34(3):322- 329.CHEN Cheng-jun,LIU Ming,CHEN Xiao-wei,et al.Development of global contact searching method based on octree algorithm [J].Chinese Journal of Computational Mechanics,2017,34(3):322- 329.
[12] 陈成军,柳明,陈小伟,等.一种基于面-面局部搜索的接触算法 [J].计算力学学报,2018,35(1):105- 110.CHEN Cheng-jun,LIU Ming,CHEN Xiao-wei,et al.A contact-impact algorithm based on the segment-to-segmant local searching method [J].Chinese Journal of Computational Mechanics,2018,35(1):105- 110.
[13] SANTASUSANA M,IRAZABAL J,ONATE E,et al.The double hierarchy method.A parallel 3D contact method for the interaction of spherical particles with rigid FE boundaries using the DEM [J].Computational Particle Mechanics,2016,3(3):407- 428.
[14] 钟阳,钟志华,李光耀,等.机械系统接触碰撞界面显式计算的算法综述 [J].机械工程学报,2011,47(13):44- 58.ZHONG Yang,ZHONG Zhi-hua,LI Guang-yao,et al.Review on contact algorithms calculating the contact-impact interface in mechanical system with the explicit FEM [J].Chinese Journal of Mechanical Engineering,2011,47(13):44- 58.
[15] BONET J,PERAIRE J.An alternating digital tree (ADT) algorithm for 3D geometric searching and intersection problems [J].International Journal for Numerical Me-thods in Engineering,1991,1(1):1 - 17.
[16] KUBLI W,REISSNER J,Optimization of sheet-metal for-ming processes using the special purpose program AUTOFORM [J].Journal of Materials Processing Technolo-gy,1995,50(1/2/3/4):292- 305.
[17] 李鹏,高振勋,蒋崇文,等.重叠网格装配中的一种改进ADT搜索方法 [J].北京航空航天大学学报,2017,43(6):1182- 1190.LI Peng,GAO Zhen-xun,JIANG Chong-wen,et al.Improved ADT searching method in overlapping grid assembly [J].Journal of Beijing University of Aeronautics and Astronautics,2017,43(6):1182- 1190.