Strain estimation from single forms of distorted fossils-a computer graphics and MATLAB approach

Most of the existing methods of strain analysis can estimate strain in a single form of distorted brachiopod, or trilobite provided independent evidence, such as the association of the fossil with cleavage and/or stretching lineation is available for inferring the direction of maximum principal strain. This article proposes a simple computer graphics based method and its MATLAB code that determine the minimum amount of strain in a single distorted fossil form even if data for inferring the maximum principal strain direction are lacking. Our method is a rapid computer-graphics alternative to some of the existing analytical methods. In a distorted fossil form of original bilateral symmetry, the relative senses of angular shears along the hinge line and the median line are mutually opposite to each other. It follows, therefore, that the maximum principal strain direction lies within the acute angle between the hinge and the median lines in the plane of the fossil. Using this principle, our method performs several simulations such that each simulation retrodeforms the distorted fossil by assuming a particular orientation, lying within the acute angle between the hinge line and the median line, as the potential direction of the maximum principal strain. Each simulation of retrodeformation yields a potential strain ratio. The distribution of all the potential strain ratios, obtained by assuming different orientations as the potential directions of the maximum strain, is typically a parabola-like curve with a distinct vertex that corresponds to the minimum amount of strain in the distorted fossil. An entirely computer graphical approach is somewhat time-intensive because it involves a large number of retrodeformational simulations. We, therefore, give a MATLAB code, namely, the Minstrain, that rapidly retrodeforms the fossil and determines the minimum strain with precision