Reg. Penetrations in a Contact Analysis Code
I have developed a FORTRAN code to simulate the molding process of glass to make lenses. Hence, it involves a contact analysis between the glass and molds. The contact is modeled using the master and slave node concept and the mold is assumed to be rigid. For the application of this code, we are looking at highly precise dimensions to capture the curve of the mold on the glass after molding.
Now, I get very good results (curve shape) in a general sense (macro sense). But when I scalse down the position of nodes of the flexible body (glass) to microns, I see a lot of fluctuation/oscillation of the nodes and this happens due to penetrations of the glass into the rigid mold. These penetrations cause huge profile deviations as the accuracy we are looking here is of the order of microns.
I have tried changing the penalty parameter to a higher value, tried to offset penetration but always end up with the small penetrations. Can anyone suggest a way to get aroung this? How do I obtain absolutely no penetration? How do I obtain the exact curve of the mold surface on the glass?