the present invention in one example corresponding to fig2 a volume of 16 u00d7 16 u00d7 8 m 3 is used a fault consisting of two planar surfaces is introduced to the volume the first planar surface dips 45 degrees w and strikes 0 degrees i e n s the second planar surface dips 45 degrees sw and strikes n45w in this example a constant restoration slip vector is applied to every point on the first planar surface of the upper fault block the restoration vectors have equal x and z components and zero y component and represents sliding of the upper fault block against the fault surface note that the y axis is parallel to the n s compass direction and the z axis is positive downward the x axis is perpendicular to the plane containing the directions of the x and z axes the restoration vectors are based on restoring a particular horizontal line on the fault surface a line that represents the intersection of a horizontal horizon and the fault surface many points on the fault surface but not those on the intersection may be subject to an inequality constraint and the repulsion method is used to enforce the constraint the boundary of the fault blocks other than the fault surfaces are under a given traction in order to maintain equilibrium of forces involved in this test zero traction is used to simulate the case of no resistance to restoration where a portion of the fault is a planar surface the deformation of the fault block consists only of space translation resulting in no strain in another example corresponding to fig3 5 the.