How can I simulate nonlinear stiffness (not bilinear) of a restraint if I have data?