This paper presents a methodology for contact detection between convex quadric surfaces using its implicit equations. With some small modifications in the equations, one can model superellipsoids, superhyperboloids of one or two sheets and supertoroids. This methodology is to be implemented on a multibody dynamics code, in order to simulate the interpenetration between mechanical systems, particularly, the simulation of collisions with motor vehicles and other road users, such as cars, motorcycles and pedestrians. The contact detection of two bodies is formulated as a convex nonlinear constrained optimization problem that is solved using two methods, an Interior Point method (IP) and a Sequential Quadratic Programming method (SQP), coded in MATLAB and FORTRAN environment, respectively. The objective function to be minimized is the distance between both surfaces. The design constraints are the implicit superquadrics surfaces equations and operations between its normal vectors and the distance itself. The contact points or the points that minimize the distance between the surfaces are the design variables. Computational efficiency can be improved by using Bounding Volumes in contact detection pre-steps. First one approximate the geometry using spheres, and then Oriented Bounding Boxes (OBB). Results show that the optimization technique suits for the accurate contact detection between objects modelled by implicit superquadric equations.