Pre-contoured dorsal and volar radial shaft fracture fixation plates were obtained from the manufacturer (Acumed™). In each plate, an array of 27, 1.5 mm holes were drilled through the plate across its surface (Fig. 1). These holes allowed physical measurements of the gap distribution (level of fit) between the surface of the plate and the surface of the bone. Three-dimensional CAD models of the plates were also obtained from the manufacturer, and an identical array of 1.5 mm holes was produced in each virtual CAD model (Fig. 1). These virtual models, together with virtual models of the bones, provided virtual measures of plate-to-bone fit and allowed a direct comparison to the physical measurements.
Six fresh frozen cadaver arms were thawed to room temperature and the radii were extracted and cleaned of soft tissue. Each radius was visually inspected to confirm that no prior fracture or deformities were present. They were coated with non-reflective paint and the midpoint between the proximal and distal ends of the bones were marked (Fig. 2). In addition, three small 1.5 mm diameter holes were drilled at each end of the bone and 1.5 mm metal pins were press-fitted into each hole protruding approximately 1.5 mm from the surface of the bone (Fig. 2). These pins were used as fiduciary markers to reproduce the physical alignment of the plate on the bone in the 3D printed bones and in the virtual models. The bones were scanned using a laser scanner (Artec Space Spider Scanner™; accuracy 0.05 mm; resolution 0.1 mm) and the scanned data were processed (Geomagic™) to obtain 3D models of the bones (Fig. 2). Three dimensional prints of each of the bones were produced (Connex Object500 by Stratasys, material VeroPureWhite RDG837, Resolution 30 microns, Accuracy 200 microns) and each of the 3D prints were scanned to obtain virtual models of the 3D printed bones.
The dorsal and volar plates were fixed to the dorsal and volar sides of each bone with the plate center aligned with the previously marked mid-distance between the distal and proximal ends of the bone (Fig. 3). While maintaining this central alignment, the plate could still be slightly adjusted circumferentially on the bone to improve fit as judged and accomplished by a fellowship-trained hand surgeon. Once this alignment was established, the plates were secured with one screw on each end. Care was taken not to overtighten the screws so as not to deform the plate (Fig. 3, PD and PV). The distances between the central and most distal holes, on each side of the plate, and the corresponding three fiduciary markers (d1 to d6 in Fig. 3) was measured with a caliper (Mitutoyo™, 0.01 mm accuracy). These distances were used to reproduce the position of the plates on each of the three-dimensional printed bones, as well as to reproduce the position of the virtual plates on the virtual bones (Fig. 3, VD and VV).
Using a digital depth-gage micrometer (Mitutoyo™, 0.01 mm accuracy) with a spindle diameter of 1.5 mm, the distance between the top surface of the plate (the one away from the bone) and the surface of the bone was measured at each of the 27, 1.5-mm drilled holes (Fig. 4, A). Each distance measurement was repeated three times and the average of the three measurements was used in the subsequent calculations. The gap distance between the bone and the plate at each of the 27 locations was then calculated based on the micrometer measures and the known plate dimensions. Using the central distal holes-to-fiduciary markers measurements (d1 to d6 in Fig. 3), the plate positions on the bones were reproduced on each of the 3D printed bones. Then, after fixation of the plate to the 3D printed bone, the physical plate-to-bone gap measurement at each of the 27 holes was repeated. Using the distances to the fiduciary markers (d1 to d6 in Fig. 3), the position of the plates on the virtual model of each bone was reproduced (Fig. 3, VD and VV). Then the virtual plate-to-bone distance at each of the 27 holes was measured by aligning a 1.5 mm diameter cylinder into each hole and measuring the distance from the underside of the plate at the specific hole to the point where the cylinder contacts the surface of the bone (Fig. 4, B). This allowed direct comparison of the plate-to-bone distance distribution between the physical and virtual conditions.
In each virtual plate (dorsal and volar), the limited bone contact region (the visible slightly elevated region in the physical models and the marked regions on each virtual plate in Fig. 1) was identified and marked (Geomagic™). This slightly elevated region on the plate was designed to be the only surface to contact the bone. For each of these regions, distance maps [22] describing the distance distribution between the bone limited contact region of the plate and the corresponding surface of the bone were obtained for both the dorsal side and the volar side of the bone (Fig. 5).
Using an optimal alignment algorithm based on least square error criterion in Geomagic™, for each specimen, the virtual model of the physical bone and the virtual model of the 3D printed bone were optimally aligned to each other. Distance maps describing the distance distribution between the two were generate (Fig. 6).
Data Processing
The 3D geometry of the optimally aligned virtual model of each physical bone and the corresponding virtual model of the 3D printed bone were compared to each other. The comparison was based on the distance distribution over the entire surfaces of the models described as distance maps (Fig. 6). From these, for each specimen, the mean distance and the standard deviation of the distance distribution were calculated. In addition, the volume of the virtual model of the physical bone and the 3D printed bone were calculated, and a paired t-test was performed to determine if there were significant differences in volume between the physical bone and the corresponding 3D printed bone.
The level of fit of the fracture fixation plate to the radius was quantified as the average and standard deviation of the plate-to-bone distance as measured across the 27 -hole array in the plates. For the virtual case, it was also quantified as the average distance and standard deviation across the distance between the bone limited contact region of the plate and the corresponding surface of the bone. Repeated Measure Analysis of Variance (ANOVA) was performed on the average and standard deviation distance between the plate and the bone for three different conditions. Physical bones, virtual bones, and 3D printed bones. The statistical significance level was set at p < 0.05. The analysis was performed separately for the volar side and for the dorsal side of the radius. A separate paired t-test was conducted for the virtual plate assessment to compare the level of fit as obtained from the 27-hole array measurements to that obtained from the distance map. This analysis provided a measure of how well the 27-hole array measurements approximated the full bone limited contact region.