I'm trying to model impact simulation of an upper beam rotates and hit another beam below it. I have used multibody dynamic physics along with the contact formulation using the penalty factor method. The forces are due gravity only. I'm using comsol 5.0 and as I have read in the comsol documents (structural mechanics module) that the contact formulation does not work with impact problems and we need always to check conservation of energy and momentum to check the results. On the other hand, I have solved the problem analytically using conservation of energy and momentum. The shear in comsol, was higher than the analytical solution, but the moment and deflection were smaller than the analytical. In addition, I'm having large contact force due to impact leading to increase shear and reduces the velocity after impact which leads to decrease moment and deflection in comsol as compared with analytical solution. I just want to match the results because I'm getting twice the shear from comsol and 1/6 moment and deflection from comsol as compared with analytical solution.

Any suggestions would be very helpful, thanks in advance.