Liakou, AnnaDetournay, EmmanuelDenoël, Vincent2016-06-022016-06-022016-06-02https://hdl.handle.net/11299/180885These MATLAB codes demonstrate the method proposed by the paper of the same name. HeMo evaluates the dynamic response of a cantilever beam against two symmetrical stops ((end)-point and/or horizontal wall). The impact motion is simulated by two methods: a) Moreau's midpoint method and b) Moreau-Jean's theta method. The second code, Rayleigh, measures the response of a Rayleigh beam for unit impulse load with beta=0.1 (Fourier series - 6 modes).A computationally efficient technique to simulate the dynamic response of a beam colliding with rigid obstacles is described in this paper. The proposed method merges three key concepts. First, a low order discretization scheme that maximizes the number of nodes of the discrete model, where impacts are detected, at the expense of the degree of continuity of the constructed displacement field. Second, the constrained problem is transformed into an unconstrained one by formulating the impact by means of a complementarity Signorini's law involving the impulse generated by the collision and the pre- and post- impact velocity linked via a coefficient of restitution. Third, Moreau's midpoint time-stepping scheme developed within the context of colliding rigid bodies, is used to advance the solution.beam impact obstaclesDatasethttp://dx.doi.org/10.13020/D6F30G