Topologies of large deformation contact-aided compliant mechanisms (CCMs), with self and mutual contact, exemplified via path generation applications, are designed using the continuum synthesis approach. Design domain is parameterized using honeycomb tessellation. Assignment of material to each cell, and generation of rigid contact surfaces, are accomplished via suitably sizing and positioning negative circular masks using the stochastic hill-climber search. To facilitate contact analysis, boundary smoothing is implemented. Mean value coordinates are employed to compute shape functions, as many regular hexagonal cells get degenerated into irregular, concave polygons as a consequence of boundary smoothing. Both geometric and material nonlinearities are considered. The augmented Lagrange multiplier method with a formulated active set strategy is employed to incorporate both self and mutual contact. Synthesized contact-aided compliant continua trace paths with single, and importantly, multiple kinks and experience multiple contact interactions pertaining to both self and mutual contact modes. Copyright © 2019 by ASME.