Robust explicit computational strategies based on penalty method for large-deformation impact problems

CHU Dongyang; RONG Yufei; ZHOU Zhangtao; WU Xingxing; WANG Jun; WANG Haikun

doi:10.11883/bzycj-2024-0073

Article Contents

Article Navigation > Explosion And Shock Waves > 2024 > Accepted Manuscript

CHU Dongyang, RONG Yufei, ZHOU Zhangtao, WU Xingxing, WANG Jun, WANG Haikun. Robust explicit computational strategies based on penalty method for large-deformation impact problems[J]. Explosion And Shock Waves. doi: 10.11883/bzycj-2024-0073

Citation:

CHU Dongyang, RONG Yufei, ZHOU Zhangtao, WU Xingxing, WANG Jun, WANG Haikun. Robust explicit computational strategies based on penalty method for large-deformation impact problems[J]. Explosion And Shock Waves. doi: 10.11883/bzycj-2024-0073

Citation:

CHU Dongyang, RONG Yufei, ZHOU Zhangtao, WU Xingxing, WANG Jun, WANG Haikun. Robust explicit computational strategies based on penalty method for large-deformation impact problems[J]. Explosion And Shock Waves. doi: 10.11883/bzycj-2024-0073

PDF( 1643 KB)

Robust explicit computational strategies based on penalty method for large-deformation impact problems

doi: 10.11883/bzycj-2024-0073

CHU Dongyang^{1, 2
,},
RONG Yufei^{1, 2},
ZHOU Zhangtao^{1, 2},
WU Xingxing^{1, 2},
WANG Jun^{1, 2},
WANG Haikun^{1, 2}

1.
China Ship Scientific Research Center, Wuxi 214082, Jiangsu, China
2.
TaiHu Laboratory on Deep Sea Technology and Science, Wuxi 214082, Jiangsu, China

Received Date: 2024-03-18
Rev Recd Date: 2024-08-09

Available Online: 2024-08-13

Abstract

Abstract

To improve the accuracy and robustness of the explicit FEM algorithm based on penalty method for simulating large deformation contact-impact problem, a new large-deformation non-penetration contact algorithm based on forward incremental displacement central difference (FIDCD) was developed. On the one hand, according to FIDCD, the solving step of the dynamic equation was decomposed into an estimated step without considering contact and a correction step considering contact constraint. At the current moment, a contact force was applied thorough the penalty method to make the deformation of entities satisfy the non-penetration condition. The contact force was calculated by a soft constraint penalty stiffness, which helped to maintain stability of contact localization. It enhanced the numerical accuracy of the explicit contact computation. On the other hand, to accurately calculate the large-deformation internal force of the next moment while only obtaining the displacement, the internal force term of the dynamic equation was mapped to a known configuration for solution based on the arbitrary reference configurations (ARC) theory. It avoided using the values of variables at intermediate configuration to approximate them, thereby improving the numerical accuracy of the large deformation computation. More rigorous contact algorithms and geometric nonlinear solution strategy can effectively suppress mesh distortion and non-physical penetration between entities during large-deformation impact simulation. This thus improved the robustness of the new explicit algorithm. Finally, the computational program written according to the new developed algorithm was applied to simulate several impact and penetration examples with different impact velocities. By comparing the simulation results with those obtained from commercial software, the correctness of the developed algorithm and computational program was verified. At the same time, it can also be proven that the algorithm proposed is more robust in simulating high-speed and large-deformation impact problems than the classical explicit contact-impact algorithm based on the frog jump center difference scheme combining with penalty method.
- contact-impact,
- explicit finite element method,
- contact algorithm,
- geometric large deformation,
- penalty method

FullText(HTML)

References(29)

References

[1]	BELYTSCHKO T, LIU W K, MORAN B, et al. Nonlinear finite elements for continua and structures [M]. 2nd ed. Chichester: John Wiley & Sons Inc. , 2014: 330–335.
[2]	钟阳, 钟志华, 李光耀, 等. 机械系统接触碰撞界面显式计算的算法综述 [J]. 机械工程学报, 2011, 47(13): 44–58. DOI: 10.3901/JME.2011.13.044. ZHONG Y, ZHONG Z H, LI G Y, et al. Review on Contact algorithms calculating the contact-impact interface in mechanical system with explicit FEM [J]. Journal of Mechanical Engineering, 2011, 47(13): 44–58. DOI: 10.3901/JME.2011.13.044.
[3]	ZHOU X, SHA D, TAMMA K K. A robust consistent configuration framework and formulation for 3D finite strain dynamic impact problems [J]. Computer Methods in Applied Mechanics and Engineering, 2008, 197(51/52): 4571–4590. DOI: 10.1016/j.cma.2008.06.001.
[4]	张雄, 王天舒, 刘岩. 计算动力学 [M]. 2版. 北京: 清华大学出版社, 2015: 249–265. ZHANG X, WANG T S, LIU Y. Computational dynamics [M]. 2nd ed. Beijing: Tsinghua University Press, 2015: 249–265.
[5]	HALLQUIST J O. LS-DYNA theory manual [M]. Livermore: Livermore Software Technology Corporation, 2006: 523–556.
[6]	陈成军, 陈小伟, 柳明. 接触-碰撞算法研究进展 [J]. 计算力学学报, 2018, 35(3): 261–274. DOI: 10.7511/jslx20160817002. CHEN C J, CHEN X W, LIU M. Review of research progress in contact-impact algorithms [J]. Chinese Journal of Computational Mechanics, 2018, 35(3): 261–274. DOI: 10.7511/jslx20160817002.
[7]	BOURAGO N G, KUKUDZHANOV V N. A review of contact algorithms [J]. Izv. RAN, MTT, 2005(1): 45–87.
[8]	Abaqus. Analysis user’s manual, version 6.14 [M]. Providence: Dassault Systemes Simulia Corporation, 2014: 151–155.
[9]	ZHONG Z H. Finite element procedures for contact-impact problems [M]. New York: Oxford University Press, 1993: 233–251.
[10]	WANG F J, WANG L P, CHENG J G, et al. Contact force algorithm in explicit transient analysis using finite-element method [J]. Finite Elements in Analysis and Design, 2007, 43(6/7): 580–587. DOI: 10.1016/j.finel.2006.12.010.
[11]	SHA D, TAMMA K K, LI M. Robust explicit computational developments and solution strategies for impact problems involving friction [J]. International Journal for Numerical Methods in Engineering, 1996, 39(5): 721–739. DOI: 10.1002/(SICI)1097-0207(19960315)39:5<721::AID-NME865>3.0.CO;2-J.
[12]	KOLMAN R, KOPAČKA J, GONZÁLEZ J A, et al. Bi-penalty stabilized technique with predictor–corrector time scheme for contact-impact problems of elastic bars [J]. Mathematics and Computers in Simulation, 2021, 189: 305–324. DOI: 10.1016/j.matcom.2021.03.023.
[13]	KOPAČKA J, TKACHUK A, GABRIEL D, et al. On stability and reflection-transmission analysis of the bipenalty method in contact-impact problems: a one-dimensional, homogeneous case study [J]. International Journal for Numerical Methods in Engineering, 2018, 113(10): 1607–1629. DOI: 10.1002/nme.5712.
[14]	SEWERIN F, PAPADOPOULOS P. On the finite element solution of frictionless contact problems using an exact penalty approach [J]. Computer Methods in Applied Mechanics and Engineering, 2020, 368: 113108. DOI: 10.1016/j.cma.2020.113108.
[15]	MOHERDAUI T F, NETO A G, WRIGGERS P. A second-order penalty-based node-to-segment contact using the Virtual Element Method [J]. Finite Elements in Analysis and Design, 2024, 237: 104183. DOI: 10.1016/j.finel.2024.104183.
[16]	王福军. 冲击接触问题有限元法并行计算及其工程应用 [D]. 北京: 清华大学, 2000: 72–82. WANG F J. Parallel computation of contact-impact problems with FEM and its engineering application [D]. Beijing: Tsinghua University, 2000: 72–82.
[17]	ZHOU X, SHA D, TAMMA K K. On a new concept and foundations of an arbitrary reference configuration (ARC) theory and formulation for computational finite deformation applications—Part Ⅰ: elasticity [J]. International Journal for Computational Methods in Engineering Science and Mechanics, 2006, 7(5): 331–351. DOI: 10.1080/15502280600790264.
[18]	ZHOU X, SHA D, TAMMA K K. On a new concept and foundations of an arbitrary reference configuration (ARC) theory and formulation for computational finite deformation applications—Part Ⅱ: elasto-plasticity [J]. International Journal for Computational Methods in Engineering Science and Mechanics, 2006, 7(5): 353–367. DOI: 10.1080/15502280600790314.
[19]	ZHOU X M, SHA D S, TAMMA K K, et al. A consistent configuration formulation involving continuum damage mechanics based Lagrangian hydrodynamic computational framework for 3D high- and hypervelocity impact/damage/penetration analysis [C]//47th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference. Newport, USA: AIAA, 2006. DOI: 10.2514/6.2006-1758.
[20]	SUN D W, LIU C, HU H Y. Dynamic computation of 2D segment-to-segment frictional contact for a flexible multibody system subject to large deformations [J]. Mechanism and Machine Theory, 2021, 158: 104197. DOI: 10.1016/j.mechmachtheory.2020.104197.
[21]	SUN D W, LIU C, HU H Y. Dynamic computation of 2D segment-to-segment frictionless contact for a flexible multibody system subject to large deformation [J]. Mechanism and Machine Theory, 2019, 140: 350–376. DOI: 10.1016/j.mechmachtheory.2019.06.011.
[22]	SERROUKH H K, MABSSOUT M, HERREROS M I. Updated Lagrangian Taylor-SPH method for large deformation in dynamic problems [J]. Applied Mathematical Modelling, 2020, 80: 238–256. DOI: 10.1016/j.apm.2019.11.046.
[23]	ZHENG X C, SEAID M, PISANÒ F, et al. A material point/finite volume method for coupled shallow water flows and large dynamic deformations in seabeds [J]. Computers and Geotechnics, 2023, 162: 105673. DOI: 10.1016/j.compgeo.2023.105673.
[24]	RÖTHLIN M, KLIPPEL H, WEGENER K. Meshless methods for large deformation elastodynamics [EB/OL]. arXiv: 1807.01117. (2018-07-05)[2024-03-18]. https://doi.org/ 10.48550/arXiv.1807.01117. DOI: 10.48550/arXiv.1807.01117.
[25]	WRIGGERS P. Computational contact mechanics [M]. 2nd ed. Berlin: Springer, 2006: 11–30. DOI: 10.1007/978-3-540-32609-0.
[26]	JOHNSON G R, COOK W H. Fracture characteristics of three metals subjected to various strains, strain rates, temperatures and pressures [J]. Engineering Fracture Mechanics, 1985, 21(1): 31–48. DOI: 10.1016/0013-7944(85)90052-9.
[27]	HEINSTEIN M W, MELLO F J, ATTAWAY S W, et al. Contact-impact modeling in explicit transient dynamics [J]. Computer Methods in Applied Mechanics and Engineering, 2000, 187(3/4): 621–640. DOI: 10.1016/S0045-7825(99)00342-4.
[28]	Livermore Software Technology Corporation. LS-DYNA^® keyword user’s manual, volume I [M]. Livermore: Livermore Software Technology Corporation, 2006: 1407–1170.
[29]	NOUR-OMID B, WRIGGERS P. A note on the optimum choice for penalty parameters [J]. Communications in Applied Numerical Methods, 1987, 3(6): 581–585. DOI: 10.1002/cnm.1630030620.