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

Volume 45 Issue 3

Mar. 2025

Turn off MathJax

Article Contents

Abstract

References

Article Navigation > Explosion And Shock Waves > 2025 > 45(3): 034201

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, 2025, 45(3): 034201. 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, 2025, 45(3): 034201. 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, 2025, 45(3): 034201. doi: 10.11883/bzycj-2024-0073

PDF( 1714 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

Publish Date: 2025-03-05

Abstract

Abstract

To improve the accuracy and robustness of the explicit finite element algorithm based on penalty method for simulating large deformation impact-contact problem, a new large-deformation non-penetration contact algorithm based on forward incremental displacement central difference (FIDCD) is developed. On the one hand, according to FIDCD, the solving step of the dynamic equation is decomposed into an estimated step without considering contact and a correction step considering contact constraint. At the current moment, a contact force is applied through the penalty method to make the deformation of entities satisfy the non-penetration condition. The contact force is calculated by a soft constraint penalty stiffness, which helps to maintain stability of contact localization. It refines 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 is mapped to a known configuration for solution based on the arbitrary reference configurations (ARC) theory. It avoids 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 improves the robustness of the new explicit algorithm. Finally, the computational program written according to the new developed algorithm is 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 is 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 impact-contact algorithm based on the frog jump center difference scheme combining with penalty method.
- impact-contact,
- 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]. Mechanics of solids, 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 Ⅰ [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.

Relative Articles

[1]	REN Kai, ZHOU Hongjing, YANG Chen. A two-step iterative method for damage calculation of a ship hullsubjected to underwater close-up non-contact explosion[J]. Explosion And Shock Waves, 2023, 43(4): 044201. doi: 10.11883/bzycj-2022-0116
[2]	LI Ping, MA Tiechang, XU Xiangzhao, MA Tianbao. A GPU parallel staircase finite difference mesh generation algorithm based on the ray casting method[J]. Explosion And Shock Waves, 2020, 40(2): 024201. doi: 10.11883/bzycj-2019-0344
[3]	Li Dian, Zhu Xi, Hou Hailiang, Zhong Qiang. Finite element analysis of load characteristic of liquid-filled structure subjected to high velocity long-rod projectile penetration[J]. Explosion And Shock Waves, 2016, 36(1): 1-8. doi: 10.11883/1001-1455(2016)01-0001-08
[4]	Li Yan-yan, Zheng Zhi-jun, Yu Ji-lin, Wang Chang-feng. Finite element analysis on deformation modes of closed-cell metallic foam[J]. Explosion And Shock Waves, 2014, 34(4): 464-470. doi: 10.11883/1001-1455(2014)04-0464-07
[5]	CHEN Er-yun, ZHAO Gai-ping, DAI Ren, MA Da-wei. Adiscontinuousfiniteelementmethod forcompressiblemulti-componentflow[J]. Explosion And Shock Waves, 2010, 30(4): 419-423. doi: 10.11883/1001-1455(2010)04-0419-05
[6]	ZHANG Xiao-tian, JIA Guang-hui, HUANG Hai. Simulationofhypervelocity-impactdebrisclouds usingaLagrangeFEM withnodeseparation[J]. Explosion And Shock Waves, 2010, 30(5): 499-504. doi: 10.11883/1001-1455(2010)05-0499-06
[7]	ZHANG Zu-gen, LI Ying-lei, LI Ying-hua, CHEN Xi-meng. Influences of bar/specimen contact surfaces indentation on strain measurement in SHPB experiments[J]. Explosion And Shock Waves, 2009, 29(6): 573-578. doi: 10.11883/1001-1455(2009)06-0573-06
[8]	SUN Zhen-sheng, ZHANG Shi-ying. Discontinuous Galerkin method and computation of fluid-rigid interaction[J]. Explosion And Shock Waves, 2008, 28(1): 80-85. doi: 10.11883/1001-1455(2008)01-0080-06
[9]	ZHENG Bo, WANG An-wen. Finite element analysis for elastic-plastic dynamic postbuckling of bars subjected to axial impact[J]. Explosion And Shock Waves, 2007, 27(2): 126-130. doi: 10.11883/1001-1455(2007)02-0126-05
[10]	WANG Ji, WANG Xiao-jun, BIAN Liang. Linking of smoothed particle hydrodynamics method to standard finite element method and its application in impact dynamics[J]. Explosion And Shock Waves, 2007, 27(6): 522-528. doi: 10.11883/1001-1455(2007)06-0522-07
[11]	HU Jin-sheng, YANG Xiu-min, ZHOU Zao-sheng, ENG Guo-qiang, TANG De-gao. Experimental investigation on contact explosion damage effect to fiber reinforced concrete slab with soil bedding[J]. Explosion And Shock Waves, 2005, 25(2): 157-162. doi: 10.11883/1001-1455(2005)02-0157-06
[12]	REN Bo, WANG Cheng. Parallel impact/penetration finite element method simulation techniques[J]. Explosion And Shock Waves, 2005, 25(3): 260-264. doi: 10.11883/1001-1455(2005)03-0260-05