Convergence Behavior of Several Solvers in GAMS and Matlab

In this post I wanted to document the difference between the cadoux model and the chrono model for a few simple test cases.

Test Problems

Ball sliding on plane

In this test problem there is a ball in perfect contact with a plane, the ball has an initial velocity of $-2 \frac{m}{s}$ in the $x$ direction. friction is .1, step size = .005, with a step size of .01, path returns a weird solution.

Box sliding on plane

In this test problem there is a box in perfect contact with a plane, the box has an initial velocity of $-2 \frac{m}{s}$ in the $x$ direction. The friction coefficient is .3;s

Box sliding on $65^\circ$ Incline

In this test problem there is a box in perfect contact with a plane, it starts at rests and then slides down for one second. friction is .2

Box sliding to a stop on $25^\circ$ Incline

In this test problem there is a box in perfect contact with a plane, it starts with velocity of $-1 \frac{m}{s}$ parallel to the incline. It sides and comes to a stop, friction is .7