Abstract:To solve the contact/impact problem of flexible bodies in arbitrary shape, the finite element method is widely used to discretize the contact bodies. In the finite element method, two contact models are mainly used to compute the contact force, i.e., penalty function method and additional constraint method, which are different in constraint imposition strategy. The penalty function method regards the contact effect as a force function of local penetration at the contact point and its rate. This method has gained significant popularity because it does not bring extra dimensions to the dynamic equations and does not need to solve constraint equations either. However, as the non-penetration constraint is not precisely satisfied in the contact process when using the penalty function method, the accuracy of the numerical simulation depends on the choice of the penalty parameter. On the other hand, the additional constraint method can strictly satisfy the contact constraint condition by introducing the Lagrange multipliers into the dynamic equations, but the method poses some numerical difficulties due to the additional effort required to solve the multipliers. Considering the advantages and disadvantages of the two contact methods, the interactive mode method is proposed. This method divides the whole model into local static module and main dynamics module. The static module establishes a local finite element model of the contact region to compute the contact force, and the main dynamics module is used to obtain the kinematic variables of the whole body. In the simulation, the two modules are coupled by exchanging displacements and forces in each time step. In the current integration step, the main dynamics module provides the displacements of the boundaries of the local contact region at first, the values are transferred to the local finite model to compute the contact force next, and then the contact force is fed back to the dynamics module for the calculation of the next step. The proposed method combines the advantages of both the additional constraint method and the penalty function method, in which not only the artificial selection of penalty parameter is avoided, but also the non-penetration constraint of local contact region is satisfied and the numerical solution is convenient. The validity of the proposed method is verified by the comparison between simulation results and experimental results of a rod-plate impact case. Furthermore, a multi-point impact problem of a slider sliding in the gap chute is presented to validate the proposed method of dealing with the general impact problem.

An improved point-to-surface contact algorithm with penalty method for peridynamics

Analysis of Short-range Contact Forces in Peridynamics Endowed with an Improved Nonlocal Contact Model

Penalty-based surface-to-surface contact algorithm for SPH method

A penetration correction contact model involving friction based on Peridynamics

A Contact Method for B-spline Material Point Method with Application in Impact and Penetration Problems

A Novel Method of Designing Random Rough Surface Considering Its Contact Force–Deformation Characteristic Requirement

Study of discrete element model for non-sphere particles base on super-quadrics

A generalized three-dimensional hybrid contact method for smoothed particle hydrodynamics

A Fluid/cloth Coupling Method for High Velocity Collision Simulation

INFLUENCE OF SURFACE EFFECT CORRECTION ON PERIDYNAMIC SIMULATION OF DYNAMIC FRACTURES IN BRITTLE MATERIALS

“Touch‐aware” contact model for peridynamics modeling of granular systems

An efficient and robust GPGPU-parallelized contact algorithm for the combined finite-discrete element method

Physics-based numerical implementation framework towards multi-scale contact problem

Dynamic Modeling Method of Flexible Bodies with Contact/impact Based on Interactive Mode

Essential insights into particle interfaces: From arbitrary switching FEM-PD coupling to effective mitigating surface effects

A novel peridynamics refinement method with dual-horizon peridynamics

A Deterministic FE Contact Analysis of 3D Rough Surfaces with Textures and Comparison with Classic Statistical Contact Models

A New Deterministic Model for Mixed Lubricated Point Contact with High Accuracy

Contact Algorithms for the Material Point Method in Impact and Penetration Simulation

Numerical Modeling on Crack Propagation Based on a Multi-Grid Bond-Based Dual-Horizon Peridynamics

A Peridynamic Model for Contact Problems Involving Fracture