Reaction Constraint Calculation
For each mass point, for each surface:
Transform into local (undeformed) coordinates
Perform inside/outside test
If inside, calculate gradient vector
Transform vector to global coordinates
Replace perpendicular force component with RC