Thanks Salomon,
As far as I know, the SHAKE algorithm can be used to constrain the bonds on LAMMPS either, but not for linear angles. I could not find right now, but I remember seeing a long discussion about this here at the mail list, maybe a year ago.
Regards,
Cassiano