Ephemeral penalty functions for contact-impact dynamics

The use of penalty functions to treat a class of structural contact-impact problems is investigated, with emphasis on ones in which the impact phenomena are primarily nondestructive in nature and in which only the gross characterization of the response is required. The dynamic equations of motion are integrated by the difference method. The penalty is represented as an ephemeral fictitious nonlinear spring that is inserted on anticipation of contact. The magnitude and variation of the penalty force is determined through energy balancing considerations. The 'bell shape' of the penalty force function for positive gap was found to be satisfactory, as it depends on only two parameters that can be directly assigned the physical meaning of force and distance. The determination of force law parameters by energy balance worked well. The incorporation of restitution coefficients by the area balancing method yielded excellent results, and no substantial modifications are anticipated. Extensional penalty springs are obviously sufficient for the simple examples treated.