Determination of Poisson Ratio of Bovine Extraocular Muscle by Computed X-Ray Tomography

The Poisson ratio (PR) is a fundamental mechanical parameter that approximates the ratio of relative change in cross sectional area to tensile elongation. However, the PR of extraocular muscle (EOM) is almost never measured because of experimental constraints. The problem was overcome by determining changes in EOM dimensions using computed X-ray tomography (CT) at microscopic resolution during tensile elongation to determine transverse strain indicated by the change in cross-section. Fresh bovine EOM specimens were prepared. Specimens were clamped in a tensile fixture within a CT scanner (SkyScan, Belgium) with temperature and humidity control and stretched up to 35% of initial length. Sets of 500–800 contiguous CT images were obtained at 10-micron resolution before and after tensile loading. Digital 3D models were then built and discretized into 6–8-micron-thick elements. Changes in longitudinal thickness of each microscopic element were determined to calculate strain. Green's theorem was used to calculate areal strain in transverse directions orthogonal to the stretching direction. The mean PR from discretized 3D models for every microscopic element in 14 EOM specimens averaged 0.457 ± 0.004 (SD). The measured PR of bovine EOM is thus near the limit of incompressibility.


Introduction
Since extraocular muscles (EOMs) are manipulated mechanically during strabismus surgery to correct binocular misalignment, the mechanical properties of the EOMs should be understood in order to optimize surgical results. With increasing demand for accuracy in simulation of orbital mechanics, �nite element analysis (FEA) is becoming increasingly attractive. However, casual estimation of mechanical parameters in FEA can lead to serious errors in simulation. In order to accurately determine the biomechanical properties of orbital tissues, scientists in the �eld have employed a variety of experimental techniques. While conventional tensile elongation tests have been performed to investigate the uniaxial force and length relationship for EOMs [1][2][3][4][5], micro/nano indentation has permitted measurement of the compressive modulus of other orbital tissues [6,7]. Despite such efforts [1,[5][6][7][8][9][10][11], many material parameters of orbital tissues have yet to be de�ned. Experimental technique and the theoretical constitutive framework should be appropriate for each tested orbital tissue. For instance, triborheometry, which treats a solid from rheometric perspective, was employed for characterizing amorphous specimens such as orbital connective and fatty tissue [10]. A variety of other biomechanical methods [12] have been employed to characterize more comprehensive constitutive models for EOMs that capture their timedependent relationships between stress and strain [1,5,8].
e Poisson ratio (PR) is a critical mechanical parameter required to de�ne comprehensively the elastic behavior of a material. e PR is the ratio of the transverse contraction strain to the axial extension strain, which can be obtained during simple tensile elongation. Normally, the PR of a material ranges between 0 and 0.5, depending on the material's compressibility. However, the PR can be higher than 0.5 or lower than 0 for materials having complex matrices and inner structures [13,14]. Most so tissues are considered to be elastomeric materials with high bulk modulus relative to Young's modulus, so the PR is expected to approximate 0.5 [15].
In general, the PR can be measured by static or dynamic methods. Static methods, such as classical tensile or compressive testing, are most widely used in solid mechanics [15,16]. In static determinations, the PR is calculated from transverse and axial deformations due to uniaxial stress. For dynamic determination, the PR is determined from the natural frequency of the transverse and axial waves in the material [17][18][19], most commonly elicited by ultrasound perturbation.
e PR has been typically assumed to be between 0.35 and 0.49 for so tissues [20][21][22]. However, it has been estimated that a 20% error in the PR would result in errors of 3.8% and 4.4% in the biaxial �exural strength and the indentation modulus of a material, respectively [15]. Such errors could propagate and compound during the iterative computations in FEA. Clearly, there is a need to minimize errors by accurate experimental determination of the PR. By employing a novel X-ray computed tomographic (CT) imaging method for precise determination of strain, we aimed to extract accurate static PR for EOMs.

Specimen Preparation.
Fresh heads of adult cows were obtained from a nearby abattoir. In the laboratory, orbits were carefully dissected for extraction of EOM and connective tissue. Transport time from abattoir to laboratory was approximately 30 min; the additional time elapsed to dissect the EOMs averaged 3 hrs. Aer extraction, EOMs were maintained in lactated Ringer's solution at 37 ∘ C. To minimize axial damage to EOM �bers, each specimen was initially prepared in the shape of an approximately 7 mm long prism with a 4 mm × 2 mm cross-section. For consistency, samples were prepared from the transverse center of each EOM. Given that clamping of both ends was necessary, the actual tested length was the 10 mm of middle portion of each specimen and in every case avoided the terminal tendon. Specimen preparation time was approximately 45 min.

Experiment.
A high-resolution micro-CT scanner (Model 1172, SkyScan, Belgium) incorporating a tensile loading �xture was used to image deformed and undeformed states of 14 freshly prepared bovine EOM specimens. Scanning was accomplished by revolving and longitudinally translating the specimen and tensile �xture between a �xed, collimated X-ray source and a �xed detector (Figure 1). Once loaded in the tensile �xture, the undeformed specimen was �rst imaged at 10 micron spatial resolution. Aer the specimen was elongated 30%-35%, which is well within the linear elastic region [1], the deformed specimen was again imaged. In order to prevent dehydration, corn oil was applied on each specimen before placement into the tensile �xture. Figure 1 shows schematics of undeformed and deformed states of the specimen in the tensile �xture.

3D Reconstruction of EOMs.
Aer EOM specimens were scanned, 500 cross-sectional area (CSA) images for undeformed and 800 CSA images for deformed states were used to create 3D reconstruction of the entire length of each EOM specimen using Matlab (Version R2010a, e MathWorks, Inc., Massachusetts) image processing tools and SolidWorks CAD soware (version 2011, Dassault Systèmes SolidWorks Corp., Massachusetts). More deformed than undeformed image planes were required since the deformed specimen was elongated. e CSA in each image plane was connected to that in the adjacent planes in order to generate 3D reconstruction by using the lo feature in Solidworks. (1), the PR , for materials undergoing deformations exceeding 1% [23], is expressed as the negative ratio of transverse to axial true strain [24,25]:

Poisson Ratio. As shown in
where and are the transverse and axial engineering strains, respectively. Equation (1) can be rearranged and can be expressed as shown in (2): Recognizing the fact that an in�nitesimal element in the CSA undergoes a plane deformation when the EOM is subjected to loading, the transverse strain in (2) can be shown as ln( o ) and (2) can be rearranged as (3) with o and being the length of side of the element in undeformed and deformed con�guration: Equation (3) is valid under the assumption of isotropy or transverse isotropy, which is appropriate for an EOM [25]. Hence the CSA of the square in the deformed con�guration can be expressed as (4) where 0 is the initial CSA: As reported by Vergari et al. [25], the instantaneous CSA can be expressed as (5) by summation of all the elements in the EOM CSA: where and o are the instantaneous and initial EOM CSA values, respectively. e CT scanner employed was speci�cally designed to image CSAs suitable for 3D reconstruction. Aer 3D reconstruction for both deformed and undeformed states, each model was then uniformly discretized into Deformed specimen X-ray X-ray x z y F 1: Both undeformed and deformed states of EOM specimens were imaged. As shown above, the X-ray direction was orthogonal to the revolving specimen �xture. 8 micron thick elements. e CSA for each element was then computed using Green's theorem as shown in (6): Finally, the PR for each discretized element was calculated using (7): . ln o ln 1 .
For more precise evaluation of PR, the CAD surface reconstruction excluded regions near clamping plates that are in�uenced by the clamping forces. Aer the PR was computed for all the elements within each specimen using (7), the average PR was calculated for each of the 14 specimens tested. Figure 2 shows the 3D reconstruction from 500 CSAs (Figure 2(a)) used to build the undeformed model of an EOM specimen (Figure 2(b)). Once 3D reconstruction was completed for each specimen for both undeformed and deformed states, the models were discretized into elements with uniform thickness of 6 microns and 8.1 microns, respectively. During preliminary experiments it was veri�ed that all 6 anatomical EOMs exhibited similar PR values. us PR values for 6 EOMs were not differentiated by anatomical EOM. Figure 3 contrasts undeformed and deformed states of the same EOM specimen.

3D Reconstruction.
Assuming uniform stretch throughout each specimen, the PR for each element within each EOM specimen was F 3: Finished 3D models for both undeformed and deformed states of an EOM specimen. e original length of the specimen was 7 mm, and the �nal length of the specimen aer the 35% deformation was 9.45 mm. computed from the CSAs and actual elongated length. Figure  4 shows the average PR values for each of the 14 specimens.

Discussion
Micro-CT imaging of EOM specimens effectively characterized the PR for bovine EOM during tensile elongation, which is a critical mechanical parameter. e present investigation is the �rst to evaluate bovine EOM PR directly from CSA measurements by noncontact imaging during large deformation. Prior studies were based on measurement of small, linear transverse deformations [26,27]. In the present investigation, mean PR of 14 bovine EOM was computed to be 0.457 ± 0.004 (SD). While isotropic media cannot have PRs exceeding 0.5, orthotropic or transversely isotropic materials (such as tendons) sometimes have PRs exceeding 1 [25]. A PR value approximating 0.5 has interesting implications for tendon behavior: assuming conservation of mass, a PR value smaller than 0.5 implies an increase in volume during loading leading to a decrease in density.
On the other hand, a PR exceeding 0.5 implies a volume reduction and an increase in density. Still, the mean PR for all EOM specimens tested in the present study was 0.457, so EOM under axial loading can be considered to a good approximation to behave as an incompressible material. e small volume variations measured might be caused by water loss, as suggested by Lynch et al. [27]. However, since in the present experiment corn oil was employed to coat the specimen to avoid dehydration, measured small volume variations are probably due to internal rearrangements of the �ber structures [25]. e precise PR value for EOM reported in the present study should facilitate quantitative modeling of ocular motor biomechanics and is represented in a theoretical framework practical for graphical simulation of quasistatic ocular motility using FEM. e present paper introduces micro-CT imaging as a noncontacting approach to compute the PR from specimen geometry. As it has been presented, micro-CT imaging technology coupled with 3D reconstruction based on the volumetric specimen changes allows more accurate evaluation of Poisson's ratio for bovine EOM specimens.

Conclusion
e current investigation describes a method to determine the PR for bovine EOM based upon CSAs from discretized deformed and undeformed specimens obtained from micro-CT imaging during quasistatic loading. e study also demonstrated that 3D reconstruction of micro CT imaging can be successfully performed for both deformed and undeformed EOM specimens. In the present study, the PR, a critical parameter for quantitative mechanical characterization of so tissues that is necessary for modeling and simulation, was determined for bovine EOM specimens to be near the criterion for incompressibility.

Con�ict of �nterests
e authors declare that none of them has any �nancial or personal relationships with people of organization that can inappropriately in�uence his work or the conclusions drawn from this investigation.