# The Influence of Pin Deviation on the Fracture Correction and the Fixator Adjustment with Sensitivity and Kinematic Analysis

Authors: Xia Zhao; Jianfeng Li
Journal: BioMed Research International
 (2018)

---

## Abstract

Fracture correction is important orthopedics operation and can be performed by unilateral external fixator. Due to pin deviations being inevitable during the operation, searching the impacts of pin deviation on the fracture correction and fixator joint adjustment are important. This study puts forward evaluation index with sensitivity analysis for investigating the impact of three orientation deviations and three position deviations on fracture correction. Meanwhile, based on a clinical case from a 28-year-old female, the influence of pin deviation on the adjustment of fixator joints was analyzed by inverse kinematic method. Different pin deviations have different sensitivity; the sensitivity of orientation deviation is relatively larger than position deviation. The existence of pin deviation will result in the change of the adjustment value of fixator joints. In addition, the experiments of seven situations were established to verify the analytical results. This study revealed the sensitivity of different pin deviations which can also be used to predict the adjustment value of fixator joints and the accuracy of fracture correction. This research also helps to reduce operation time and decrease damage to soft tissue by reducing the frequency of inserting pins.


---

## Body

## 1. Introduction

Bone fracture is becoming a very prominent research topic mainly due to aging of population and traffic accident [1]. Fracture reduction is usually performed in orthopedics doctors to ensure that fractured bones were restored to their original location. It has been shown that fracture site using mechanical means through external fixators may have beneficial effects on bone healing and remodeling [2–4]. If fracture has not been accurately corrected, fracture mal-union or other fracture dislocation problems can lead to complications, such as early degenerative disease caused by abnormal joint contact pressures [5–7]. Complex types of malalignment often involve pediatric orthopaedic patients and they are often treated by osteotomy [8]. Regulating the rotational and translational joints of the fixator is usually necessary to reduce fracture and correct complex deformities [9, 10]. But the pins with different position deviations and orientation deviation will affect the accuracy of fracture correction. Clinically, the procedure of inserting pins into bone is usually performed artificially before an external fixator is applied for reduction. Under ideal situations, pins are inserted perpendicular to the bone and the fixator is applied in neutral or straight position [11]. However, in the actual situation, the existence of pin deviations is inevitable and it is difficult to insert the pins into the bone accurately and appropriately as ideal situation. Because the fixator has its own limitation, removing the attached fixator, redrilling new pin-holes, and then installing fixator on the new drilled pin-holes are sometimes necessary, but this kind of behavior will be harmful to bone tissue and extend the operation time [12]. If orthopedic surgeon knows how the pin deviation affects fracture correction and the fixator joint regulation rule, then they will be very skilled using external fixator. Meanwhile, if surgeons had more freedom to better arrange the pin-bone relationship and were not distracted to the subsequent reduction again, the quality of operation would be greatly improved.To the author’s knowledge, although many prior studies discussed the sensitivities of two biomechanical methods for fracture healing [13], the adjustability of two unilateral external fixators (Othofixs Dynamic Axial fixator and Dynafixs Standard Tibia fixator) for fracture reduction [14, 15] and the possibility of their nonaxial dynamization after fracture reduction [16], however, where pin deviation has the greatest impact on the accuracy of fracture reduction and the influence of pin deviations on the adjustment values of fixator joint, have been rarely investigated due to lack of a theoretical analysis.So the aim of this study was, based on a clinical case, to discuss the effect of different pin deviations on fracture correction through sensitivity analysis method, through this study where pin deviation that has the greatest impact on fracture dislocation can be found. Secondly, how to adjust the fixator joints could help correct fracture completely when different pin deviations exist, which was also investigated through inverse kinematic analysis method where joints values can be directly obtained through analytic formulas. The establishment of the experiment also proves the correctness of the analytical results. This is of practical importance since from a physician’s standpoint knowing how to accurately correct fracture dislocation is by understanding the relationship between the pin deviations and the adjustment of fixator joints, which can also help to provide theoretical guidance and improve the design of fixator.


## 2. Methods

### 2.1. Kinematic Model

A 6-degree-of-freedom (6-DOF) unilateral external fixator (Dyna-extor, Korean) was used in this study. In this study, only one bone pin at each bone segment is considered for rigid-body in kinematic analysis. This usage method has been reported in a previous article [17].Taking proximal bone (B) as a reference, supposing distal bone has threeZ1, Z2, Z3 axes, having three orientation deviations (φ1  φ2  φ3) and three position deviations (s1  s2  s3), which may exist in the process of inserting pins into bone, shown in Figure 1. Detailed joint parameters of unilateral external fixator are shown in Table 1. The parameters of pin deviation for unilateral external fixator are shown in Table 2.Table 1
The joint parameters of unilateral external fixator.
Joint variable Kinematic pair limit value of joints Rotational axis a 5 Distal revolute pair 0~360 (degree) Z 5 a 6 Distal revolute pair 0~120 (degree) Z 6 d Moving pair 0~40(mm) Z 7 a 8 Proximal revolute pair 0~360 (degree) Z 8 a 9 Proximal revolute pair 0~120 (degree) Z 9 a 10 Proximal revolute pair 0~360 (degree) Z 10Table 2
The parameters of pin deviation for unilateral external fixator.
Symbol of pin deviations pin deviations Error range of angle deviation Angle deviation axis φ 1 Distal angel deviation 0~20 (degree) Z 1 φ 2 Distal angel deviation 0~20 (degree) Z 2 φ 3 Distal angel deviation 0~20 (degree) Z 3 s 1 Distal position deviation 0~10mm Z 1 s 2 Distal position deviation 0~10mm Z 2 s 3 Distal position deviation 0~10mm Z 3Figure 1
(a) Kinematic model of bone and fixator system, “A”—distal bone; “1”—distal pin; “2”—distal pin clamp; “3”—the telescoping mechanism; “4”—proximal pin clamp; “5”—proximal pin; “B”—proximal bone. (b): The physical model of fixator and bone system,L0=L5=75mm; L2=L3=30mm; L1=L4=65mm; D=90mm; “a5, a6”— two rotational variables near distal bone, “a8  a9, a10”—three rotational variables near proximal bone; “d”—translational variable; “φ1, φ2, φ3”—angle deviations of pin; “s1, s2, s3”— position deviations of pin.
(a) (b)


### 2.2. Establishment of Sensitivity Analysis Model

The precision of fracture reduction is affected by the position deviation and orientation deviation of pin. The effects of different pin deviations on the precision of fracture reduction are different. In order to find out which pin deviation has the greatest impact on the precision of fracture correction, the sensitivity analysis of pin deviation is needed, which describes the effect of degree of the small variation of pin deviations on the fracture correction [18].The error calculation model of the fixator device can be established, which can be expressed explicitly.(1)B=M4AM54M65M76M87M98M109M1210MB12=Ms1,s2,s3,φ1,φ2,φ3B is the continuous differentiable function ofs1,s2,s3,φ1,φ2,φ3; the first-order Taylor series expansion of (1) can be expressed as follows:(2)Ms1+Δs1,s2+Δs2,s3+Δs3,φ1+Δφ1,φ2+Δφ2,φ3+Δφ3=Ms1,s2,s3,φ1,φ2,φ3+∂M∂s1Δs1+∂M∂s2Δs2+∂M∂s3Δs3+∂M∂φ1Δφ1+∂M∂φ2Δφ2+∂M∂φ3Δφ3+OΔs12+Δs22+Δs32+Δφ12+Δφ22+Δφ32In (2), Δs1, Δs2, Δs3, Δφ1, Δφ2, Δφ3 are the small amount of the s1,s2,s3,φ1,φ2,φ3; removing the high order terms, the following formula can be obtained.(3)Ms1+Δs1,s2+Δs2,s3+Δs3,φ1+Δφ1,φ2+Δφ2,φ3+Δφ3=Ms1,s2,s3,φ1,φ2,φ3+∂M∂s1Δs1+∂M∂s2Δs2+∂M∂s3Δs3+∂M∂φ1Δφ1+∂M∂φ2Δφ2+∂M∂φ3Δφ3=B+ΔBAnd then the sensitive matrix of pin deviations can be obtained as follows; in (4), J is error sensitivity matrix or error transfer matrix.(4)ΔM=ΔB=∂M∂s1Δs1+∂M∂s2Δs2+∂M∂s3Δs3+∂M∂φ1Δφ1+∂M∂φ2Δφ2+∂M∂φ3Δφ3=JΔs1,Δs2,Δs3,Δφ1,Δφ2,Δφ3Inverse kinematic analysis method was adopted to solveJ. In the mathematical analysis of the kinematic chain, MBA is the transformation matrix from the proximal bone segment to the distal bone segment which could be expressed as follows:(5)MBA=M4AM54M65M87M98M109M1210MB12=nxoxaxpxnyoyaypynzozazpz0001M 4   A refers to homogeneous transformation matrix from distal bone to distal pin.M54 and M65 are the rotational matrix of distal revolute joint, M76 denotes the translational matrix of the telescoping mechanism, and M87M98M109 are the rotational matrix of proximal ball joint. M1210 expresses matrix from proximal pin to proximal bone and MB12 represents displacement transformation matrix from the proximal pin to fracture site.Inverse kinematic analysis method was adopted to solve J. Firstly, multiply [MB12]−1 [M1210]−1 [M109]−1 on both right sides of (5) and obtain (6). [MB12]−1 [M1210]−1 [M109]−1 are the inverse matrix of [MB12] [M1210] [M109], respectively.(6)MBAMB12-1M1210-1M109-1=M4AM54M65M76M87M98(7)H=M4AM54M65M76M87M98=H11H12H13H14H21H22H23H24H31H32H33H34H41H42H43H44(8)M=MBAMB12-1M1210-1M109-1=M11M12M13M14M21M22M23M24M31M32M33M34M41M42M43M44Equation (7) denotes the transformation matrix from the coordinate system ZA to Z9 in a counterclockwise direction; (8) refers to transformation matrix from the coordinate system of ZA to Z9 in a clockwise direction in Figure 1(a). Given that the equations on both sides of (6) are equal, H and M are equal and first derivatives of H and M are equal, so (9)-(10) is derived.(9)H14=M14H24=M24H34=M34According to (9), a5, a6, d, can be solved. Detailed process can be shown in appendix.(10)H14′=M14′H24′=M24′H34′=M34′According to (10), (11) can be obtained as follows:(11)A1φ1′+A2φ2′+A3φ3′+A4s3′-A5s2′=Bs1′+A6Bs3′+A7Bφ2′-A8Bφ3′B1φ1′+B2φ2′+B3φ3′+B4s2′+B5s3′=B6Bφ1′+B7Bφ2′+B8Bφ3′+B9Bs2′-B10Bs3′s1′-C1φ2′+C2φ3′-C3s3′=C4Bφ′-C5Bφ2′+C6Bφ3′+C7Bs2′+C8Bs3′Multiplying both sides of (5) by [MB14]−1[M1410]−1[M109]−1[M98]−1 [M87]−1[M76]−1, (12) is obtained as follows.(12)M4AM54M65=MBAMB14-1M1410-1M109-1M98-1M87-1M76-1(13)P=M4AM54M65=P11P12P13P14P21P22P23P24P31P32P33P34P41P42P43P44(14)K=MBAMB14-1M1410-1M109-1M98-1M87-1M76-1=K11K12K13K14K21K22K23K24K31K32K33K34K41K42K43K44P denotes the transformation matrix fromZA to Z6 in a counterclockwise direction; K refers to the transformation matrix from ZA to Z6 in a clockwise direction, as shown in Figure 1(a). P and K are equal; thus, first derivatives of K and P are also equal, so (14) is derived.(15)K14=P14K24=P24K34=P34K13=P13According to (15), a8, a9, a10 can be solved. Detailed process can be expressed in appendix.(16)K14′=P14′K24′=P24′K34′=P34′According to (16), (17) can be obtained as follows:(17)D1Bφ2′+D2Bφ3′+Bs1′+D3Bs3′=D4φ1′+D5φ2′+D6φ3′-D7s2′+D8s3′E1Bφ1′+E2Bφ2′+E3Bφ3′+E4Bs2′-E5Bs3′=E6φ1′-E7φ2′+E8φ3′+E9s2′+E10s3′F1Bφ1′-F2Bφ2′+F3Bφ3′+F4Bs2′+F5Bs3′=F6φ2′+F7φ3′+s1′-F8s3′According to (11) and (17), (18) can be obtained as follows:(18)A1A2A30-A5A4B1B2B30B4B50-C1C210-C3D4D5D60-D7D8E6-E7E80E9E100F6F710-F8×φ1φ2φ3s1s2s3=0A7-A810A6B6B7B80B9-B10C4-C5C60C7C80D1D210D3E1E2E30E4E5F1-F2F30F4F5×Bφ1Bφ2Bφ3Bs1Bs2Bs3According to (18), (19) can be obtained as follows:(19)0A7-A810A6B6B7B80B9-B10C4-C5C60C7C80D1D210D3E1E2E30E4E5F1-F2F30F4F5-1A1A2A30-A5A4B1B2B30B4B50-C1C210-C3D4D5D60-D7D8E6-E7E80E9E100F6F710-F8×φ1φ2φ3s1s2s3=Bφ1Bφ2Bφ3Bs1Bs2Bs3Error transfer matrix can be obtained as follows:(20)J=0A7-A810A6B6B7B80B9-B10C4-C5C60C7C80D1D210D3E1E2E30E4E5F1-F2F30F4F5-1A1A2A30-A5A4B1B2B30B4B50-C1C210-C3D4D5D60-D7D8E6-E7E80E9E100F6F710-F8The detailed solution process and concrete expressions of A1 ~A8, B1 ~B10, C1 ~C8, D1 ~D8, E1 ~E10, F1 ~F8 in (11) and (17)~(20) can be explained in the appendix.


### 2.3. Sensitivity Evaluation Index of Pin Deviation

The influence degree of pin deviation on the fracture correction can be expressed by the error sensitivity coefficient. The corresponding amount of various types of deviation of pins is shown in Table3.Table 3
The corresponding amount of various types deviation of pins.
Types of pin deviation Fracture dislocation Deviation transfer matrix Error model in fracture site Δ s 1 Δ B s 1 J s 1 Δ B s 1 = J s 1 Δ s 1 Δ s 2 Δ B s 2 J s 2 Δ B s 2 = J s 2 Δ s 2 Δ s 3 Δ B s 3 J s 3 Δ B s 3 = J s 3 Δ s 3 Δ φ 1 Δ B φ 1 J φ 1 Δ B φ 1 = J φ 1 Δ φ 1 Δ φ 2 Δ B φ 2 J φ 2 Δ B φ 2 = J φ 2 Δ φ 2 Δ φ 3 Δ B φ 3 J φ 3 Δ B φ 3 = J φ 3 Δ φ 3In (21), qs is displacement deviation at the fracture site; qφ is angular misalignment at the fracture site shown as Figure 2. Js is error transfer matrix of position deviation of pin; Jφ is error transfer matrix of orientation deviation of pin.(21)qs=Jsδsqφ=JφδφFigure 2
A clinical case. (a) Preoperative X-ray film for a case of tibia deformity. (b) Postoperative anteroposterior X-ray film of limb. (c) Photo showing the limb after application of the unilateral external fixator.
(a) (b) (c)According to (21), using Lagrange operator to establish Lagrange equation,(22)Ls=δsTJsTJsδs-λsδsTδs-1Lφ=δφTJφTJφδφ-λφδφTδφ-1In (22), λs and λφ are Lagrange multipliers. According to (22), construct the extreme value of Lagrange conditional equation:(23)∂Ls∂δs=0:JsTJsδs-λsδs=0∂Lφ∂δφ=0:JφTJφδφ-λφδφ=0In (23), Lagrange multipliers of λs and λφ are the eigenvalues of JsTJs and JφTJφ, respectively. The square of displacement deviation qs and angular misalignment qφ are shown as follows:(24)qs2=ΔqsTΔqs=δsTJsTJsδs=δsTλsδs=λsqφ2=ΔqφTΔqφ=δφTJφTJφδφ=δφTλφδφ=λφThe extreme values ofΔBs and ΔBφ are the root of λs and λφ, respectively:(25)qs=λs=JsTJsqφ=λφ=JφTJφΔ B s and ΔBφ have three extreme components, respectively, which can be written as follows:(26)fs1=λs1=Js1TJs1fs2=λs2=Js2TJs2fs3=λs3=Js3TJs3In (26), λs1, λs2, and λs3 are the sensitivity coefficient of displacement dislocation ΔBs at the fracture site in XA, YA, ZA, three directions, respectively, caused by position deviation of pin.(27)fφ1=λφ1=Jφ1TJφ1fφ2=λφ2=Jφ2TJφ2fφ3=λφ3=Jφ3TJφ3In (25), λφ1, λφ2, λφ3 are the sensitivity coefficient of angular misalignment ΔBφ at the fracture site in three directions of XA, YA, ZA, respectively, which caused orientation deviation of pin.Takeλs and λφ as sensitivity evaluation index to assess the impact of position deviation of pin and orientation deviation of pin on the displacement dislocation and angular dislocation, respectively.


### 2.4. Clinical Application

In order to better study the effect of pin deviation on fixator joint adjustment and verify the sensitivity of pin deviation, a case was used.In clinical practice, 28-year-old female was injured in a road traffic accident and treated with a unilateral external fixator. A 28-year-old female motorcyclist was brought into hospital after having been hit by a car from her right side. She had no medical history, took no regular medication, and had no allergies. The patient provided written informed consent for publication. The case has AP radiograph showing tibia deformity of 20° internal rotation and 20° external rotation and has 15° angular deformity and 5mm displacement deviations in three planes. By adjusting the fixator joints ofa5, a6, d, a10, a9, and a8 with -6°, 20.8°, 9.3mm, -10.6°, -3.3°, and -12° respectively, the fracture deformity can be completely corrected.Follow-up continued for a mean of 18 months (10–37 months). On the second day after surgery all patients had radiographs (anteroposterior, lateral, and mortise views) (Figure2(a)). Partial weight bearing was allowed six weeks postoperatively and was gradually increased according to clinical and radiological evidence of union up to full weight bearing after complete union (Figure 2(b)) at a mean of 12 weeks (10–16 weeks). X-ray examination was repeated every month for a period of six months to evaluate bone union and fracture consolidation. Removal of the unilateral fixator was done for all cases after complete union at a mean period of 14 weeks (12–17 weeks). Removal of the fixator was followed by muscle strengthening exercises and physiotherapy.


## 2.1. Kinematic Model

A 6-degree-of-freedom (6-DOF) unilateral external fixator (Dyna-extor, Korean) was used in this study. In this study, only one bone pin at each bone segment is considered for rigid-body in kinematic analysis. This usage method has been reported in a previous article [17].Taking proximal bone (B) as a reference, supposing distal bone has threeZ1, Z2, Z3 axes, having three orientation deviations (φ1  φ2  φ3) and three position deviations (s1  s2  s3), which may exist in the process of inserting pins into bone, shown in Figure 1. Detailed joint parameters of unilateral external fixator are shown in Table 1. The parameters of pin deviation for unilateral external fixator are shown in Table 2.Table 1
The joint parameters of unilateral external fixator.
Joint variable Kinematic pair limit value of joints Rotational axis a 5 Distal revolute pair 0~360 (degree) Z 5 a 6 Distal revolute pair 0~120 (degree) Z 6 d Moving pair 0~40(mm) Z 7 a 8 Proximal revolute pair 0~360 (degree) Z 8 a 9 Proximal revolute pair 0~120 (degree) Z 9 a 10 Proximal revolute pair 0~360 (degree) Z 10Table 2
The parameters of pin deviation for unilateral external fixator.
Symbol of pin deviations pin deviations Error range of angle deviation Angle deviation axis φ 1 Distal angel deviation 0~20 (degree) Z 1 φ 2 Distal angel deviation 0~20 (degree) Z 2 φ 3 Distal angel deviation 0~20 (degree) Z 3 s 1 Distal position deviation 0~10mm Z 1 s 2 Distal position deviation 0~10mm Z 2 s 3 Distal position deviation 0~10mm Z 3Figure 1
(a) Kinematic model of bone and fixator system, “A”—distal bone; “1”—distal pin; “2”—distal pin clamp; “3”—the telescoping mechanism; “4”—proximal pin clamp; “5”—proximal pin; “B”—proximal bone. (b): The physical model of fixator and bone system,L0=L5=75mm; L2=L3=30mm; L1=L4=65mm; D=90mm; “a5, a6”— two rotational variables near distal bone, “a8  a9, a10”—three rotational variables near proximal bone; “d”—translational variable; “φ1, φ2, φ3”—angle deviations of pin; “s1, s2, s3”— position deviations of pin.
(a) (b)


## 2.2. Establishment of Sensitivity Analysis Model

The precision of fracture reduction is affected by the position deviation and orientation deviation of pin. The effects of different pin deviations on the precision of fracture reduction are different. In order to find out which pin deviation has the greatest impact on the precision of fracture correction, the sensitivity analysis of pin deviation is needed, which describes the effect of degree of the small variation of pin deviations on the fracture correction [18].The error calculation model of the fixator device can be established, which can be expressed explicitly.(1)B=M4AM54M65M76M87M98M109M1210MB12=Ms1,s2,s3,φ1,φ2,φ3B is the continuous differentiable function ofs1,s2,s3,φ1,φ2,φ3; the first-order Taylor series expansion of (1) can be expressed as follows:(2)Ms1+Δs1,s2+Δs2,s3+Δs3,φ1+Δφ1,φ2+Δφ2,φ3+Δφ3=Ms1,s2,s3,φ1,φ2,φ3+∂M∂s1Δs1+∂M∂s2Δs2+∂M∂s3Δs3+∂M∂φ1Δφ1+∂M∂φ2Δφ2+∂M∂φ3Δφ3+OΔs12+Δs22+Δs32+Δφ12+Δφ22+Δφ32In (2), Δs1, Δs2, Δs3, Δφ1, Δφ2, Δφ3 are the small amount of the s1,s2,s3,φ1,φ2,φ3; removing the high order terms, the following formula can be obtained.(3)Ms1+Δs1,s2+Δs2,s3+Δs3,φ1+Δφ1,φ2+Δφ2,φ3+Δφ3=Ms1,s2,s3,φ1,φ2,φ3+∂M∂s1Δs1+∂M∂s2Δs2+∂M∂s3Δs3+∂M∂φ1Δφ1+∂M∂φ2Δφ2+∂M∂φ3Δφ3=B+ΔBAnd then the sensitive matrix of pin deviations can be obtained as follows; in (4), J is error sensitivity matrix or error transfer matrix.(4)ΔM=ΔB=∂M∂s1Δs1+∂M∂s2Δs2+∂M∂s3Δs3+∂M∂φ1Δφ1+∂M∂φ2Δφ2+∂M∂φ3Δφ3=JΔs1,Δs2,Δs3,Δφ1,Δφ2,Δφ3Inverse kinematic analysis method was adopted to solveJ. In the mathematical analysis of the kinematic chain, MBA is the transformation matrix from the proximal bone segment to the distal bone segment which could be expressed as follows:(5)MBA=M4AM54M65M87M98M109M1210MB12=nxoxaxpxnyoyaypynzozazpz0001M 4   A refers to homogeneous transformation matrix from distal bone to distal pin.M54 and M65 are the rotational matrix of distal revolute joint, M76 denotes the translational matrix of the telescoping mechanism, and M87M98M109 are the rotational matrix of proximal ball joint. M1210 expresses matrix from proximal pin to proximal bone and MB12 represents displacement transformation matrix from the proximal pin to fracture site.Inverse kinematic analysis method was adopted to solve J. Firstly, multiply [MB12]−1 [M1210]−1 [M109]−1 on both right sides of (5) and obtain (6). [MB12]−1 [M1210]−1 [M109]−1 are the inverse matrix of [MB12] [M1210] [M109], respectively.(6)MBAMB12-1M1210-1M109-1=M4AM54M65M76M87M98(7)H=M4AM54M65M76M87M98=H11H12H13H14H21H22H23H24H31H32H33H34H41H42H43H44(8)M=MBAMB12-1M1210-1M109-1=M11M12M13M14M21M22M23M24M31M32M33M34M41M42M43M44Equation (7) denotes the transformation matrix from the coordinate system ZA to Z9 in a counterclockwise direction; (8) refers to transformation matrix from the coordinate system of ZA to Z9 in a clockwise direction in Figure 1(a). Given that the equations on both sides of (6) are equal, H and M are equal and first derivatives of H and M are equal, so (9)-(10) is derived.(9)H14=M14H24=M24H34=M34According to (9), a5, a6, d, can be solved. Detailed process can be shown in appendix.(10)H14′=M14′H24′=M24′H34′=M34′According to (10), (11) can be obtained as follows:(11)A1φ1′+A2φ2′+A3φ3′+A4s3′-A5s2′=Bs1′+A6Bs3′+A7Bφ2′-A8Bφ3′B1φ1′+B2φ2′+B3φ3′+B4s2′+B5s3′=B6Bφ1′+B7Bφ2′+B8Bφ3′+B9Bs2′-B10Bs3′s1′-C1φ2′+C2φ3′-C3s3′=C4Bφ′-C5Bφ2′+C6Bφ3′+C7Bs2′+C8Bs3′Multiplying both sides of (5) by [MB14]−1[M1410]−1[M109]−1[M98]−1 [M87]−1[M76]−1, (12) is obtained as follows.(12)M4AM54M65=MBAMB14-1M1410-1M109-1M98-1M87-1M76-1(13)P=M4AM54M65=P11P12P13P14P21P22P23P24P31P32P33P34P41P42P43P44(14)K=MBAMB14-1M1410-1M109-1M98-1M87-1M76-1=K11K12K13K14K21K22K23K24K31K32K33K34K41K42K43K44P denotes the transformation matrix fromZA to Z6 in a counterclockwise direction; K refers to the transformation matrix from ZA to Z6 in a clockwise direction, as shown in Figure 1(a). P and K are equal; thus, first derivatives of K and P are also equal, so (14) is derived.(15)K14=P14K24=P24K34=P34K13=P13According to (15), a8, a9, a10 can be solved. Detailed process can be expressed in appendix.(16)K14′=P14′K24′=P24′K34′=P34′According to (16), (17) can be obtained as follows:(17)D1Bφ2′+D2Bφ3′+Bs1′+D3Bs3′=D4φ1′+D5φ2′+D6φ3′-D7s2′+D8s3′E1Bφ1′+E2Bφ2′+E3Bφ3′+E4Bs2′-E5Bs3′=E6φ1′-E7φ2′+E8φ3′+E9s2′+E10s3′F1Bφ1′-F2Bφ2′+F3Bφ3′+F4Bs2′+F5Bs3′=F6φ2′+F7φ3′+s1′-F8s3′According to (11) and (17), (18) can be obtained as follows:(18)A1A2A30-A5A4B1B2B30B4B50-C1C210-C3D4D5D60-D7D8E6-E7E80E9E100F6F710-F8×φ1φ2φ3s1s2s3=0A7-A810A6B6B7B80B9-B10C4-C5C60C7C80D1D210D3E1E2E30E4E5F1-F2F30F4F5×Bφ1Bφ2Bφ3Bs1Bs2Bs3According to (18), (19) can be obtained as follows:(19)0A7-A810A6B6B7B80B9-B10C4-C5C60C7C80D1D210D3E1E2E30E4E5F1-F2F30F4F5-1A1A2A30-A5A4B1B2B30B4B50-C1C210-C3D4D5D60-D7D8E6-E7E80E9E100F6F710-F8×φ1φ2φ3s1s2s3=Bφ1Bφ2Bφ3Bs1Bs2Bs3Error transfer matrix can be obtained as follows:(20)J=0A7-A810A6B6B7B80B9-B10C4-C5C60C7C80D1D210D3E1E2E30E4E5F1-F2F30F4F5-1A1A2A30-A5A4B1B2B30B4B50-C1C210-C3D4D5D60-D7D8E6-E7E80E9E100F6F710-F8The detailed solution process and concrete expressions of A1 ~A8, B1 ~B10, C1 ~C8, D1 ~D8, E1 ~E10, F1 ~F8 in (11) and (17)~(20) can be explained in the appendix.


## 2.3. Sensitivity Evaluation Index of Pin Deviation

The influence degree of pin deviation on the fracture correction can be expressed by the error sensitivity coefficient. The corresponding amount of various types of deviation of pins is shown in Table3.Table 3
The corresponding amount of various types deviation of pins.
Types of pin deviation Fracture dislocation Deviation transfer matrix Error model in fracture site Δ s 1 Δ B s 1 J s 1 Δ B s 1 = J s 1 Δ s 1 Δ s 2 Δ B s 2 J s 2 Δ B s 2 = J s 2 Δ s 2 Δ s 3 Δ B s 3 J s 3 Δ B s 3 = J s 3 Δ s 3 Δ φ 1 Δ B φ 1 J φ 1 Δ B φ 1 = J φ 1 Δ φ 1 Δ φ 2 Δ B φ 2 J φ 2 Δ B φ 2 = J φ 2 Δ φ 2 Δ φ 3 Δ B φ 3 J φ 3 Δ B φ 3 = J φ 3 Δ φ 3In (21), qs is displacement deviation at the fracture site; qφ is angular misalignment at the fracture site shown as Figure 2. Js is error transfer matrix of position deviation of pin; Jφ is error transfer matrix of orientation deviation of pin.(21)qs=Jsδsqφ=JφδφFigure 2
A clinical case. (a) Preoperative X-ray film for a case of tibia deformity. (b) Postoperative anteroposterior X-ray film of limb. (c) Photo showing the limb after application of the unilateral external fixator.
(a) (b) (c)According to (21), using Lagrange operator to establish Lagrange equation,(22)Ls=δsTJsTJsδs-λsδsTδs-1Lφ=δφTJφTJφδφ-λφδφTδφ-1In (22), λs and λφ are Lagrange multipliers. According to (22), construct the extreme value of Lagrange conditional equation:(23)∂Ls∂δs=0:JsTJsδs-λsδs=0∂Lφ∂δφ=0:JφTJφδφ-λφδφ=0In (23), Lagrange multipliers of λs and λφ are the eigenvalues of JsTJs and JφTJφ, respectively. The square of displacement deviation qs and angular misalignment qφ are shown as follows:(24)qs2=ΔqsTΔqs=δsTJsTJsδs=δsTλsδs=λsqφ2=ΔqφTΔqφ=δφTJφTJφδφ=δφTλφδφ=λφThe extreme values ofΔBs and ΔBφ are the root of λs and λφ, respectively:(25)qs=λs=JsTJsqφ=λφ=JφTJφΔ B s and ΔBφ have three extreme components, respectively, which can be written as follows:(26)fs1=λs1=Js1TJs1fs2=λs2=Js2TJs2fs3=λs3=Js3TJs3In (26), λs1, λs2, and λs3 are the sensitivity coefficient of displacement dislocation ΔBs at the fracture site in XA, YA, ZA, three directions, respectively, caused by position deviation of pin.(27)fφ1=λφ1=Jφ1TJφ1fφ2=λφ2=Jφ2TJφ2fφ3=λφ3=Jφ3TJφ3In (25), λφ1, λφ2, λφ3 are the sensitivity coefficient of angular misalignment ΔBφ at the fracture site in three directions of XA, YA, ZA, respectively, which caused orientation deviation of pin.Takeλs and λφ as sensitivity evaluation index to assess the impact of position deviation of pin and orientation deviation of pin on the displacement dislocation and angular dislocation, respectively.


## 2.4. Clinical Application

In order to better study the effect of pin deviation on fixator joint adjustment and verify the sensitivity of pin deviation, a case was used.In clinical practice, 28-year-old female was injured in a road traffic accident and treated with a unilateral external fixator. A 28-year-old female motorcyclist was brought into hospital after having been hit by a car from her right side. She had no medical history, took no regular medication, and had no allergies. The patient provided written informed consent for publication. The case has AP radiograph showing tibia deformity of 20° internal rotation and 20° external rotation and has 15° angular deformity and 5mm displacement deviations in three planes. By adjusting the fixator joints ofa5, a6, d, a10, a9, and a8 with -6°, 20.8°, 9.3mm, -10.6°, -3.3°, and -12° respectively, the fracture deformity can be completely corrected.Follow-up continued for a mean of 18 months (10–37 months). On the second day after surgery all patients had radiographs (anteroposterior, lateral, and mortise views) (Figure2(a)). Partial weight bearing was allowed six weeks postoperatively and was gradually increased according to clinical and radiological evidence of union up to full weight bearing after complete union (Figure 2(b)) at a mean of 12 weeks (10–16 weeks). X-ray examination was repeated every month for a period of six months to evaluate bone union and fracture consolidation. Removal of the unilateral fixator was done for all cases after complete union at a mean period of 14 weeks (12–17 weeks). Removal of the fixator was followed by muscle strengthening exercises and physiotherapy.


## 3. Results

In this section, we discuss the sensitivities of different pin deviations. Meanwhile, based on a clinical case, the effect of pin deviation on fracture correction is also analyzed, which can verify the sensitivity of pin deviation from an analytical point of view. Moreover, the experiment was also established to verify the analytic results.


### 3.1. Simulation Results of Pin Deviations Sensitivity

#### 3.1.1. Results of Sensitivity Analysis for Single Pin Deviation

We first investigate the sensitivity of single pin deviation in the appropriate range. The bigger the sensitivity is, the greater the effect of the pin deviation on the precision of the fracture reduction is.Figure3(a) shows sensitivity curves of three single orientation deviations; it can be noted that φ2 has the largest sensitivity when compared with φ3 and φ1. Besides φ2 and φ3 in the entire range of pin deviations, the sensitivity was relatively larger than φ1. This indicates that residual fracture deformity of Bφ2 is most likely to appear. The possibility of fracture deformity of Bφ1 caused by φ1 is minimal. Meanwhile, with the increase of the pin deviation, the sensitivity of all pin deviations present an increasing trend; this indicates that the increase in pin deviation may lead to increased sensitivity for fracture dislocations.Figure 3
The sensitivity of single pin deviation. (a) The sensitivity of orientation deviation of single pin. (b) The sensitivity of positional deviation of single pin.
(a) (b)From Figure3(b), it can be seen that the sensitivity of s3 is higher than that of s1 and s2.With the increase of the position deviation, the sensitivity also displays an increasing trend.The sensitivity ofs1 is minimal, meaning that it has the least impact on fracture displacement dislocation of Bs1, and Bs3 has the biggest influence in displacement dislocation due to s3 having the biggest sensitivity compared with s2 and s1.Comparing Figures3(a) and 3(b), the maximum sensitivity is φ2. Orientation deviation of φ3, φ2, and φ1 is relatively larger than position deviation of s3, s2, and s1, respectively. These sensitivity results show that angular dislocation is more likely to occur than displacement dislocation caused by pin deviation, which imply that the orientation deviation has a greater impact on the accuracy of fracture correction than the position deviation and which may be needed to adjust more joints with more amount to ensure fracture correction completely.In order to investigate whether the sensitivity of double pin deviations has the same regularity, the following sensitivity simulation was performed.


#### 3.1.2. Results of Sensitivity Analysis for Double Deviations

The sensitivity results from Figure4 help operation doctor to seek the place of maximum sensitivity when there are two kinds of pin deviations simultaneously, as well as predict the change of the fixator’s regulation. In Figure 4(a), the sensitivity peak is 0.6986 at φ1=20°, φ2=20°; from Figures 4(b) and 4(c), the maximum sensitivity is 0.7618 and 0.6103; they appear in (φ3=20°  φ2=20°) and (φ3=20°  φ1=20°), respectively. By comparing these three values, it is revealed that the sensitivity of φ2 is greater than φ3 and φ1, because due to existing pin deviation of φ2, the sensitivity values of 0.6986 and 0.7618 are bigger than 0.6103. Meanwhile, under the condition of having φ2, the value of 0.7618 is bigger than 0.6986, indicating the sensitivity of φ3 is higher than that of φ1. Besides, from Figure 4(a), the higher sensitivity of 0.612 also appears in the (φ2=20°  φ1=0°); this illustrates again that φ2 has a greater influence on fracture correction than φ1, because the larger deviation sensitivity tends to occur to the side of φ2. For Figure 4(b), the bigger sensitivity is at (φ3=0°  φ2=20°), which is 0.7618, and for Figure 4(c) the bigger sensitivity that is 0.561 appears at (φ3=20°  φ1=0°). These results also illustrated that φ2 has larger impact than φ3 and φ1 for fracture dislocation. φ3 has a greater influence on fracture deformity correction than φ1.Figure 4
The sensitivity of position and orientation deviation of double deviations of pin. (a), (b), and (c) The sensitivity of position deviation of double deviations of pin. (d), (e), and (f) The sensitivity of orientation deviation of double deviations of pin.
(a) (b) (c) (d) (e) (f)For Figures4(d), 4(e), and 4(f), the maximum sensitivity values exist in s1=10mm, s2=10mm; s3=10mm, s2=10mm; s1=10mm, s3=10mm; and they are 0.45437, 0.6371, and 0.59342, respectively. By comparing these three values, the conclusions same as Figure 3(b) can be obtained as follows: the sensitivity of s3, s2, and s1 shows a gradual increasing trend. In Figure 4(d), the higher sensitive value appears at the s2=0mm, s1=10mm and s1=0mm, s2=10mm; they are 0.38 and 0.4218, respectively. Obviously, the impact of s2 on fracture displacement correction is greater than s1. For Figures 4(e) and 4(f), the higher sensitivity values exist at (s3=10mm s2=0mm) and (s1=0mm s3=10mm); they all incline the side of s3; this indicates that the sensitivity of s3 is bigger than s1 and s2.So we can put forward a hypothesis; as long as we know the maximum sensitivity of single deviation, the maximum sensitivity of multiple deviations can be found. In order to verify this hypothesis, the following section analyzes the sensitivity of the multiple deviations of pin.


#### 3.1.3. Results of Sensitivity Analysis for Multiple Deviations of Pin

For Figure5, the volume of blue ball and red ball indicates the sensitivity of orientation deviation and position sensitivity, which have three different orientation deviations (φ1  φ2  φ3) and three different position deviations (s1  s2  s3) at the same time, respectively. With the increase of deviation, the volume of the ball increases gradually, so the sensitivity is also increased. For Figure 5(a), the growth rate of the spherical volume in the φ2 direction was significantly larger than that in φ1 and φ3 directions. The change rate in the direction of φ1 is the smallest. For Figure 5(b), the growth rate of the spherical volume in s3 direction was relatively bigger than that in s2 and s1 directions. The change rate in the direction of s1 is also the smallest as the same results in Figures 4(d)–4(f).Figure 5
The sensitivity of multiple deviations of pin. (a) The sensitivity of multiple position deviations of pin. (b) The sensitivity of multiple orientation deviations of pin.
(a) (b)Comparing Figure5(a) with Figure 5(b), the maximum position sensitivity appears in φ1=20°, φ2=20°, φ3=20° and s1=10mm, s2=10mm, s3=10mm, respectively; the results from this figure are consistent with previous assumptions since the maximum sensitivity also appears in the point which is also the maximum sensitivity in the single deviation.Because the accuracy of fracture correction is guaranteed by adjusting fixator joints, it is necessary to investigate how to adjust fixator that can compensate fracture dislocations caused by pin deviation.


### 3.2. Results of the Effect of Pin Deviations on the Fixator Adjustment

Based on a clinical case in 2.4 (20° internal rotation, 20° external rotation and has 15° angular deformity, and 5mm displacement deviations on three planes). The effect of pin deviations on the fixator adjustment was analyzed.


#### 3.2.1. The Effect of Orientation Pin Deviation on Fixator Adjustment

Under the condition of ideal situation—having no pin deviation,a5  a6 d a10  a9  a8 can be obtained by (9) and (15), they are -5.99°, 20.8024°, 9.3537mm, -10.5515°, -3.2705°, and -12.1301°, respectively. In the same way, when we have pin deviations, fixator joints also can be solved. Next we investigate the effect of pin deviation on fixator’s joints as follows.Figure6 displays the bar charts of the effect of different orientation deviation on the fixator joint adjustment (a5  a6 d a8  a9  a10). According to the results listed in Figures 6(a)–6(f), the results can be observed that the increase of orientation deviation of φ1 will lead to rapid decrease of fixator joints of a8. Since the regulating range of a8 is only 36°, the existence of φ1 can help a8 avoid reaching the maximum value of regulation and expand the solution domain of the fixator.Figure 6
The influence of orientation deviation on the fixator joint adjustment. (a) The influence of orientation deviation on the adjustment of the fixator joint -a5. (b) The influence of orientation deviation on the adjustment of the fixator joint a6. (c) The influence of orientation deviation on the adjustment of the fixator joint d. (d) The influence of orientation deviation on the adjustment of the fixator joint a8. (e) The influence of orientation deviation on the adjustment of the fixator joint a9. (f) The influence of orientation deviation on the adjustment of the fixator joint a10.
(a) (b) (c) (d) (e) (f)It can be shown that the adjustment joint ofa5 is mainly affected by φ2 and a5 and a10 need to adjust greater quantity to ensure fracture correction accurately when compared with the condition of having no deviation, but this sharp increase in a5 and a10 will result in fracture correction incompletely, because the adjustment range of the fixator is limited. On the contrary, the existence of φ2 can make d decrease a lot; this will expand solution domain of d.Similarly, with the increase ofφ3, to ensure the fracture was corrected completely, the adjustment amount of a6 needs to decrease greatly and a9 will reduce a certain value, while φ3 did not cause a lot of change to other joints compared with ideal situation.The greater the sensitivity of pin deviation, the greater the adjustment amount of the fixator. According to the above results thatφ2 affects more joints with greater value than φ3 and φ1, it can be also proven that the sensitivity of φ2 is larger than φ3 and φ1; nevertheless the sensitivity of φ1 is minimal.In addition to this, an important property was found in Figures6(a)–6(f); when there are more deviations at the same time, the adjustment value of the joint is equal to the sum of the adjustment value of each deviation.


#### 3.2.2. The Effect of Position Deviation of Pin on the Fixator Joint Adjustment

For Figures7(a)–7(f), the effects of different position deviation on fixator adjustment are presented. As the sensitivity analysis above, there are more joints affected by s3; with the increase of s3, the joint regulation of a5 increased significantly, a10 declined with a rapid rate, but a8 and a9 show a slightly downward trend. The appearance of s3 leads a5 to reach the limit value in advance and also cause a great residual fracture dislocation and then will delay fracture healing.Figure 7
The influence of position deviation on the fixator joint adjustment. (a) The influence of position deviation on the adjustment of fixator jointa5. (b) The influence of position deviation on the adjustment of fixator joint a6. (c) The influence of position deviation on the adjustment of fixator joint d. (d) The influence of position deviation on the adjustment of fixator joint a8. (e) The influence of position deviation on the adjustment of fixator joint a9. (f) The influence of position deviation on the adjustment of fixator joint a10.
(a) (b) (c) (d) (e) (f)When the deviation ofs2 increases, a9 generated a higher regulatory value, a6 and d had lower value compared with ideal situation; the existence of s2 also makes the fixator joint a9 not capable enough to correct deformity which need greater regulation value.The presence of the position deviation ofs1 makes d increase significantly and also bring little increase in a9 and little decrease in a5. Similar to previous position sensitivity analysis, it is conceivable that s3 has impact on more fixator joints with high value; the effect of s1 on fixator joint is minimal. It also can be found that the influence of pin deviation on the fixator's regulation has superposition characteristic, meaning that just knowing the influence of each pin deviation on adjustment amount of fixator joints and accordingly how to adjust fixator can also be predicted when there are many deviations; the adjustment value of fixator is the sum of adjustment value of every pin deviation.


### 3.3. Experimental Verification

The experiment of fixator-bone system for deformity correction at ideal situation shown in Figure8 was established, in order to better verify the analytical results. In experiment, the bone deformation has 20° internal rotation, 20° external rotation, 15° angular deformity, and 5mm displacement deviations on three planes, which is the same as the clinical case. The adjustment values of fixator joints which need to correct deformity are also in accordance with clinical practice, meaning the experiment is reasonable. Moreover, seven situations were also established to verify the analytical results; experimental results are shown in Table 4. The parameters of fixator-bone system in experiment are consistent with analytical model in Figure 1. In the process of experiment, one side of bone was fixed; by adjusting fixator joints, the deformity can be completely corrected.Table 4
The adjustment values of fixator joints at seven situations for experimental results.
Four situations a 5(degree) a 6(degree) d(mm) a 8(degree) a 9(degree) a 10(degree) Ideal situation -5.5 20 9 -10 -3 -12 φ 1=20° 0.5 8.5 3 4 -9.5 -9 φ 2=20° 16 21 -8 -16.5 -6.5 -20.5 φ 3=20° 5.5 -5 7.5 -14 2.5 -9.3 s 1=10mm -5 18.5 18.5 -11.5 -1.2 -10.5 s 2=10mm -5.6 15 6.0 -12 2 -8.5 s 3=10mm -12 20 10.5 -6.5 -4.5 -5Figure 8
The experiment of fixator-bone system for bone deformity at ideal situation. (a) The experiment of fixator-bone system before deformity correction. (b) The experiment of fixator-bone system after deformity correction.
(a) (b)Table4 shows the adjustment values of fixator joints at seven situations in experiment. Take the ideal situation as reference, when compared with pin deviations of φ2 and φ3, the adjustment value of the fixator joint of a8 has the biggest discrepancy between the ideal situation and the existence of pin deviation φ1. The presence of φ2 causes a greater change than φ1 and φ3 in the amount of a5, a8, and a10. Compared with φ1 and φ2, φ3 has the biggest effect on a6 and a9. In summary, when compared with ideal situation, the existence of pin deviation will affect the regulating amount of fixator joint, and φ2 will affect more joints and a greater amount of adjustment is needed; φ1 has only bigger influences on one joint of a8 and φ3 causes bigger effect on two joints of a6 and a9. That also indicates that φ2 has the bigger sensitivity than φ1 and φ3, meaning φ2 has the biggest effect on the fixator joint adjustment. In addition, φ1 has the minimum sensitivity compared to φ2 and φ3. These experimental results are consistent with analytical results for angular deviations.For positional deviations, take the ideal situation with no pin deviations as a reference; compared withs2 and s3, the existence of s1 makes a bigger change in fixator joint d. Meanwhile, s2 has a great effect on a9 and a6. The existence of s3 makes a great contribution to the bigger changes on a5, a8, and a10 than s1 and s2. These results indicate that s3 has the biggest sensitivity on the fixator joint adjustment and s1 has minimum sensitivity, which can also verify the correctness of the analytical results.Table5 shows the discrepancy between experiment results and analytical results for fixator joint adjustment value. The errors between experiment results and analytic results are very small and within twenty percent, meaning that the experimental results can verify the correctness of the analysis results.Table 5
The errors between experimental and analytical results for fixator joint adjustment values.
Four situations a 5(degree) a 6(degree) d(mm) a 8(degree) a 9(degree) a 10(degree) Ideal situation 8% 3.8% 3.78% 5.2% 8.27% 1.05% φ 1=20° 16% 2.96% 13.4% 7.69% 1.56% 1.26% φ 2=20° 4.7% 0.213% 3.12% 1.17% 1% 1.7% φ 3=20° 1.7% 1.82% 5.4% 3.05% 0.56% 19.19% s 1=10mm 7.4% 2.32% 1.14% 5.41% 9.7% 0.52% s 2=10mm 5% 1.51% 4.11% 9.33% 3.7% 5.357% s3=10mm 0.35% 2.43% 3% 19.7% 3.17% 0.85%


## 3.1. Simulation Results of Pin Deviations Sensitivity

### 3.1.1. Results of Sensitivity Analysis for Single Pin Deviation

We first investigate the sensitivity of single pin deviation in the appropriate range. The bigger the sensitivity is, the greater the effect of the pin deviation on the precision of the fracture reduction is.Figure3(a) shows sensitivity curves of three single orientation deviations; it can be noted that φ2 has the largest sensitivity when compared with φ3 and φ1. Besides φ2 and φ3 in the entire range of pin deviations, the sensitivity was relatively larger than φ1. This indicates that residual fracture deformity of Bφ2 is most likely to appear. The possibility of fracture deformity of Bφ1 caused by φ1 is minimal. Meanwhile, with the increase of the pin deviation, the sensitivity of all pin deviations present an increasing trend; this indicates that the increase in pin deviation may lead to increased sensitivity for fracture dislocations.Figure 3
The sensitivity of single pin deviation. (a) The sensitivity of orientation deviation of single pin. (b) The sensitivity of positional deviation of single pin.
(a) (b)From Figure3(b), it can be seen that the sensitivity of s3 is higher than that of s1 and s2.With the increase of the position deviation, the sensitivity also displays an increasing trend.The sensitivity ofs1 is minimal, meaning that it has the least impact on fracture displacement dislocation of Bs1, and Bs3 has the biggest influence in displacement dislocation due to s3 having the biggest sensitivity compared with s2 and s1.Comparing Figures3(a) and 3(b), the maximum sensitivity is φ2. Orientation deviation of φ3, φ2, and φ1 is relatively larger than position deviation of s3, s2, and s1, respectively. These sensitivity results show that angular dislocation is more likely to occur than displacement dislocation caused by pin deviation, which imply that the orientation deviation has a greater impact on the accuracy of fracture correction than the position deviation and which may be needed to adjust more joints with more amount to ensure fracture correction completely.In order to investigate whether the sensitivity of double pin deviations has the same regularity, the following sensitivity simulation was performed.


### 3.1.2. Results of Sensitivity Analysis for Double Deviations

The sensitivity results from Figure4 help operation doctor to seek the place of maximum sensitivity when there are two kinds of pin deviations simultaneously, as well as predict the change of the fixator’s regulation. In Figure 4(a), the sensitivity peak is 0.6986 at φ1=20°, φ2=20°; from Figures 4(b) and 4(c), the maximum sensitivity is 0.7618 and 0.6103; they appear in (φ3=20°  φ2=20°) and (φ3=20°  φ1=20°), respectively. By comparing these three values, it is revealed that the sensitivity of φ2 is greater than φ3 and φ1, because due to existing pin deviation of φ2, the sensitivity values of 0.6986 and 0.7618 are bigger than 0.6103. Meanwhile, under the condition of having φ2, the value of 0.7618 is bigger than 0.6986, indicating the sensitivity of φ3 is higher than that of φ1. Besides, from Figure 4(a), the higher sensitivity of 0.612 also appears in the (φ2=20°  φ1=0°); this illustrates again that φ2 has a greater influence on fracture correction than φ1, because the larger deviation sensitivity tends to occur to the side of φ2. For Figure 4(b), the bigger sensitivity is at (φ3=0°  φ2=20°), which is 0.7618, and for Figure 4(c) the bigger sensitivity that is 0.561 appears at (φ3=20°  φ1=0°). These results also illustrated that φ2 has larger impact than φ3 and φ1 for fracture dislocation. φ3 has a greater influence on fracture deformity correction than φ1.Figure 4
The sensitivity of position and orientation deviation of double deviations of pin. (a), (b), and (c) The sensitivity of position deviation of double deviations of pin. (d), (e), and (f) The sensitivity of orientation deviation of double deviations of pin.
(a) (b) (c) (d) (e) (f)For Figures4(d), 4(e), and 4(f), the maximum sensitivity values exist in s1=10mm, s2=10mm; s3=10mm, s2=10mm; s1=10mm, s3=10mm; and they are 0.45437, 0.6371, and 0.59342, respectively. By comparing these three values, the conclusions same as Figure 3(b) can be obtained as follows: the sensitivity of s3, s2, and s1 shows a gradual increasing trend. In Figure 4(d), the higher sensitive value appears at the s2=0mm, s1=10mm and s1=0mm, s2=10mm; they are 0.38 and 0.4218, respectively. Obviously, the impact of s2 on fracture displacement correction is greater than s1. For Figures 4(e) and 4(f), the higher sensitivity values exist at (s3=10mm s2=0mm) and (s1=0mm s3=10mm); they all incline the side of s3; this indicates that the sensitivity of s3 is bigger than s1 and s2.So we can put forward a hypothesis; as long as we know the maximum sensitivity of single deviation, the maximum sensitivity of multiple deviations can be found. In order to verify this hypothesis, the following section analyzes the sensitivity of the multiple deviations of pin.


### 3.1.3. Results of Sensitivity Analysis for Multiple Deviations of Pin

For Figure5, the volume of blue ball and red ball indicates the sensitivity of orientation deviation and position sensitivity, which have three different orientation deviations (φ1  φ2  φ3) and three different position deviations (s1  s2  s3) at the same time, respectively. With the increase of deviation, the volume of the ball increases gradually, so the sensitivity is also increased. For Figure 5(a), the growth rate of the spherical volume in the φ2 direction was significantly larger than that in φ1 and φ3 directions. The change rate in the direction of φ1 is the smallest. For Figure 5(b), the growth rate of the spherical volume in s3 direction was relatively bigger than that in s2 and s1 directions. The change rate in the direction of s1 is also the smallest as the same results in Figures 4(d)–4(f).Figure 5
The sensitivity of multiple deviations of pin. (a) The sensitivity of multiple position deviations of pin. (b) The sensitivity of multiple orientation deviations of pin.
(a) (b)Comparing Figure5(a) with Figure 5(b), the maximum position sensitivity appears in φ1=20°, φ2=20°, φ3=20° and s1=10mm, s2=10mm, s3=10mm, respectively; the results from this figure are consistent with previous assumptions since the maximum sensitivity also appears in the point which is also the maximum sensitivity in the single deviation.Because the accuracy of fracture correction is guaranteed by adjusting fixator joints, it is necessary to investigate how to adjust fixator that can compensate fracture dislocations caused by pin deviation.


## 3.1.1. Results of Sensitivity Analysis for Single Pin Deviation

We first investigate the sensitivity of single pin deviation in the appropriate range. The bigger the sensitivity is, the greater the effect of the pin deviation on the precision of the fracture reduction is.Figure3(a) shows sensitivity curves of three single orientation deviations; it can be noted that φ2 has the largest sensitivity when compared with φ3 and φ1. Besides φ2 and φ3 in the entire range of pin deviations, the sensitivity was relatively larger than φ1. This indicates that residual fracture deformity of Bφ2 is most likely to appear. The possibility of fracture deformity of Bφ1 caused by φ1 is minimal. Meanwhile, with the increase of the pin deviation, the sensitivity of all pin deviations present an increasing trend; this indicates that the increase in pin deviation may lead to increased sensitivity for fracture dislocations.Figure 3
The sensitivity of single pin deviation. (a) The sensitivity of orientation deviation of single pin. (b) The sensitivity of positional deviation of single pin.
(a) (b)From Figure3(b), it can be seen that the sensitivity of s3 is higher than that of s1 and s2.With the increase of the position deviation, the sensitivity also displays an increasing trend.The sensitivity ofs1 is minimal, meaning that it has the least impact on fracture displacement dislocation of Bs1, and Bs3 has the biggest influence in displacement dislocation due to s3 having the biggest sensitivity compared with s2 and s1.Comparing Figures3(a) and 3(b), the maximum sensitivity is φ2. Orientation deviation of φ3, φ2, and φ1 is relatively larger than position deviation of s3, s2, and s1, respectively. These sensitivity results show that angular dislocation is more likely to occur than displacement dislocation caused by pin deviation, which imply that the orientation deviation has a greater impact on the accuracy of fracture correction than the position deviation and which may be needed to adjust more joints with more amount to ensure fracture correction completely.In order to investigate whether the sensitivity of double pin deviations has the same regularity, the following sensitivity simulation was performed.


## 3.1.2. Results of Sensitivity Analysis for Double Deviations

The sensitivity results from Figure4 help operation doctor to seek the place of maximum sensitivity when there are two kinds of pin deviations simultaneously, as well as predict the change of the fixator’s regulation. In Figure 4(a), the sensitivity peak is 0.6986 at φ1=20°, φ2=20°; from Figures 4(b) and 4(c), the maximum sensitivity is 0.7618 and 0.6103; they appear in (φ3=20°  φ2=20°) and (φ3=20°  φ1=20°), respectively. By comparing these three values, it is revealed that the sensitivity of φ2 is greater than φ3 and φ1, because due to existing pin deviation of φ2, the sensitivity values of 0.6986 and 0.7618 are bigger than 0.6103. Meanwhile, under the condition of having φ2, the value of 0.7618 is bigger than 0.6986, indicating the sensitivity of φ3 is higher than that of φ1. Besides, from Figure 4(a), the higher sensitivity of 0.612 also appears in the (φ2=20°  φ1=0°); this illustrates again that φ2 has a greater influence on fracture correction than φ1, because the larger deviation sensitivity tends to occur to the side of φ2. For Figure 4(b), the bigger sensitivity is at (φ3=0°  φ2=20°), which is 0.7618, and for Figure 4(c) the bigger sensitivity that is 0.561 appears at (φ3=20°  φ1=0°). These results also illustrated that φ2 has larger impact than φ3 and φ1 for fracture dislocation. φ3 has a greater influence on fracture deformity correction than φ1.Figure 4
The sensitivity of position and orientation deviation of double deviations of pin. (a), (b), and (c) The sensitivity of position deviation of double deviations of pin. (d), (e), and (f) The sensitivity of orientation deviation of double deviations of pin.
(a) (b) (c) (d) (e) (f)For Figures4(d), 4(e), and 4(f), the maximum sensitivity values exist in s1=10mm, s2=10mm; s3=10mm, s2=10mm; s1=10mm, s3=10mm; and they are 0.45437, 0.6371, and 0.59342, respectively. By comparing these three values, the conclusions same as Figure 3(b) can be obtained as follows: the sensitivity of s3, s2, and s1 shows a gradual increasing trend. In Figure 4(d), the higher sensitive value appears at the s2=0mm, s1=10mm and s1=0mm, s2=10mm; they are 0.38 and 0.4218, respectively. Obviously, the impact of s2 on fracture displacement correction is greater than s1. For Figures 4(e) and 4(f), the higher sensitivity values exist at (s3=10mm s2=0mm) and (s1=0mm s3=10mm); they all incline the side of s3; this indicates that the sensitivity of s3 is bigger than s1 and s2.So we can put forward a hypothesis; as long as we know the maximum sensitivity of single deviation, the maximum sensitivity of multiple deviations can be found. In order to verify this hypothesis, the following section analyzes the sensitivity of the multiple deviations of pin.


## 3.1.3. Results of Sensitivity Analysis for Multiple Deviations of Pin

For Figure5, the volume of blue ball and red ball indicates the sensitivity of orientation deviation and position sensitivity, which have three different orientation deviations (φ1  φ2  φ3) and three different position deviations (s1  s2  s3) at the same time, respectively. With the increase of deviation, the volume of the ball increases gradually, so the sensitivity is also increased. For Figure 5(a), the growth rate of the spherical volume in the φ2 direction was significantly larger than that in φ1 and φ3 directions. The change rate in the direction of φ1 is the smallest. For Figure 5(b), the growth rate of the spherical volume in s3 direction was relatively bigger than that in s2 and s1 directions. The change rate in the direction of s1 is also the smallest as the same results in Figures 4(d)–4(f).Figure 5
The sensitivity of multiple deviations of pin. (a) The sensitivity of multiple position deviations of pin. (b) The sensitivity of multiple orientation deviations of pin.
(a) (b)Comparing Figure5(a) with Figure 5(b), the maximum position sensitivity appears in φ1=20°, φ2=20°, φ3=20° and s1=10mm, s2=10mm, s3=10mm, respectively; the results from this figure are consistent with previous assumptions since the maximum sensitivity also appears in the point which is also the maximum sensitivity in the single deviation.Because the accuracy of fracture correction is guaranteed by adjusting fixator joints, it is necessary to investigate how to adjust fixator that can compensate fracture dislocations caused by pin deviation.


## 3.2. Results of the Effect of Pin Deviations on the Fixator Adjustment

Based on a clinical case in 2.4 (20° internal rotation, 20° external rotation and has 15° angular deformity, and 5mm displacement deviations on three planes). The effect of pin deviations on the fixator adjustment was analyzed.


### 3.2.1. The Effect of Orientation Pin Deviation on Fixator Adjustment

Under the condition of ideal situation—having no pin deviation,a5  a6 d a10  a9  a8 can be obtained by (9) and (15), they are -5.99°, 20.8024°, 9.3537mm, -10.5515°, -3.2705°, and -12.1301°, respectively. In the same way, when we have pin deviations, fixator joints also can be solved. Next we investigate the effect of pin deviation on fixator’s joints as follows.Figure6 displays the bar charts of the effect of different orientation deviation on the fixator joint adjustment (a5  a6 d a8  a9  a10). According to the results listed in Figures 6(a)–6(f), the results can be observed that the increase of orientation deviation of φ1 will lead to rapid decrease of fixator joints of a8. Since the regulating range of a8 is only 36°, the existence of φ1 can help a8 avoid reaching the maximum value of regulation and expand the solution domain of the fixator.Figure 6
The influence of orientation deviation on the fixator joint adjustment. (a) The influence of orientation deviation on the adjustment of the fixator joint -a5. (b) The influence of orientation deviation on the adjustment of the fixator joint a6. (c) The influence of orientation deviation on the adjustment of the fixator joint d. (d) The influence of orientation deviation on the adjustment of the fixator joint a8. (e) The influence of orientation deviation on the adjustment of the fixator joint a9. (f) The influence of orientation deviation on the adjustment of the fixator joint a10.
(a) (b) (c) (d) (e) (f)It can be shown that the adjustment joint ofa5 is mainly affected by φ2 and a5 and a10 need to adjust greater quantity to ensure fracture correction accurately when compared with the condition of having no deviation, but this sharp increase in a5 and a10 will result in fracture correction incompletely, because the adjustment range of the fixator is limited. On the contrary, the existence of φ2 can make d decrease a lot; this will expand solution domain of d.Similarly, with the increase ofφ3, to ensure the fracture was corrected completely, the adjustment amount of a6 needs to decrease greatly and a9 will reduce a certain value, while φ3 did not cause a lot of change to other joints compared with ideal situation.The greater the sensitivity of pin deviation, the greater the adjustment amount of the fixator. According to the above results thatφ2 affects more joints with greater value than φ3 and φ1, it can be also proven that the sensitivity of φ2 is larger than φ3 and φ1; nevertheless the sensitivity of φ1 is minimal.In addition to this, an important property was found in Figures6(a)–6(f); when there are more deviations at the same time, the adjustment value of the joint is equal to the sum of the adjustment value of each deviation.


### 3.2.2. The Effect of Position Deviation of Pin on the Fixator Joint Adjustment

For Figures7(a)–7(f), the effects of different position deviation on fixator adjustment are presented. As the sensitivity analysis above, there are more joints affected by s3; with the increase of s3, the joint regulation of a5 increased significantly, a10 declined with a rapid rate, but a8 and a9 show a slightly downward trend. The appearance of s3 leads a5 to reach the limit value in advance and also cause a great residual fracture dislocation and then will delay fracture healing.Figure 7
The influence of position deviation on the fixator joint adjustment. (a) The influence of position deviation on the adjustment of fixator jointa5. (b) The influence of position deviation on the adjustment of fixator joint a6. (c) The influence of position deviation on the adjustment of fixator joint d. (d) The influence of position deviation on the adjustment of fixator joint a8. (e) The influence of position deviation on the adjustment of fixator joint a9. (f) The influence of position deviation on the adjustment of fixator joint a10.
(a) (b) (c) (d) (e) (f)When the deviation ofs2 increases, a9 generated a higher regulatory value, a6 and d had lower value compared with ideal situation; the existence of s2 also makes the fixator joint a9 not capable enough to correct deformity which need greater regulation value.The presence of the position deviation ofs1 makes d increase significantly and also bring little increase in a9 and little decrease in a5. Similar to previous position sensitivity analysis, it is conceivable that s3 has impact on more fixator joints with high value; the effect of s1 on fixator joint is minimal. It also can be found that the influence of pin deviation on the fixator's regulation has superposition characteristic, meaning that just knowing the influence of each pin deviation on adjustment amount of fixator joints and accordingly how to adjust fixator can also be predicted when there are many deviations; the adjustment value of fixator is the sum of adjustment value of every pin deviation.


## 3.2.1. The Effect of Orientation Pin Deviation on Fixator Adjustment

Under the condition of ideal situation—having no pin deviation,a5  a6 d a10  a9  a8 can be obtained by (9) and (15), they are -5.99°, 20.8024°, 9.3537mm, -10.5515°, -3.2705°, and -12.1301°, respectively. In the same way, when we have pin deviations, fixator joints also can be solved. Next we investigate the effect of pin deviation on fixator’s joints as follows.Figure6 displays the bar charts of the effect of different orientation deviation on the fixator joint adjustment (a5  a6 d a8  a9  a10). According to the results listed in Figures 6(a)–6(f), the results can be observed that the increase of orientation deviation of φ1 will lead to rapid decrease of fixator joints of a8. Since the regulating range of a8 is only 36°, the existence of φ1 can help a8 avoid reaching the maximum value of regulation and expand the solution domain of the fixator.Figure 6
The influence of orientation deviation on the fixator joint adjustment. (a) The influence of orientation deviation on the adjustment of the fixator joint -a5. (b) The influence of orientation deviation on the adjustment of the fixator joint a6. (c) The influence of orientation deviation on the adjustment of the fixator joint d. (d) The influence of orientation deviation on the adjustment of the fixator joint a8. (e) The influence of orientation deviation on the adjustment of the fixator joint a9. (f) The influence of orientation deviation on the adjustment of the fixator joint a10.
(a) (b) (c) (d) (e) (f)It can be shown that the adjustment joint ofa5 is mainly affected by φ2 and a5 and a10 need to adjust greater quantity to ensure fracture correction accurately when compared with the condition of having no deviation, but this sharp increase in a5 and a10 will result in fracture correction incompletely, because the adjustment range of the fixator is limited. On the contrary, the existence of φ2 can make d decrease a lot; this will expand solution domain of d.Similarly, with the increase ofφ3, to ensure the fracture was corrected completely, the adjustment amount of a6 needs to decrease greatly and a9 will reduce a certain value, while φ3 did not cause a lot of change to other joints compared with ideal situation.The greater the sensitivity of pin deviation, the greater the adjustment amount of the fixator. According to the above results thatφ2 affects more joints with greater value than φ3 and φ1, it can be also proven that the sensitivity of φ2 is larger than φ3 and φ1; nevertheless the sensitivity of φ1 is minimal.In addition to this, an important property was found in Figures6(a)–6(f); when there are more deviations at the same time, the adjustment value of the joint is equal to the sum of the adjustment value of each deviation.


## 3.2.2. The Effect of Position Deviation of Pin on the Fixator Joint Adjustment

For Figures7(a)–7(f), the effects of different position deviation on fixator adjustment are presented. As the sensitivity analysis above, there are more joints affected by s3; with the increase of s3, the joint regulation of a5 increased significantly, a10 declined with a rapid rate, but a8 and a9 show a slightly downward trend. The appearance of s3 leads a5 to reach the limit value in advance and also cause a great residual fracture dislocation and then will delay fracture healing.Figure 7
The influence of position deviation on the fixator joint adjustment. (a) The influence of position deviation on the adjustment of fixator jointa5. (b) The influence of position deviation on the adjustment of fixator joint a6. (c) The influence of position deviation on the adjustment of fixator joint d. (d) The influence of position deviation on the adjustment of fixator joint a8. (e) The influence of position deviation on the adjustment of fixator joint a9. (f) The influence of position deviation on the adjustment of fixator joint a10.
(a) (b) (c) (d) (e) (f)When the deviation ofs2 increases, a9 generated a higher regulatory value, a6 and d had lower value compared with ideal situation; the existence of s2 also makes the fixator joint a9 not capable enough to correct deformity which need greater regulation value.The presence of the position deviation ofs1 makes d increase significantly and also bring little increase in a9 and little decrease in a5. Similar to previous position sensitivity analysis, it is conceivable that s3 has impact on more fixator joints with high value; the effect of s1 on fixator joint is minimal. It also can be found that the influence of pin deviation on the fixator's regulation has superposition characteristic, meaning that just knowing the influence of each pin deviation on adjustment amount of fixator joints and accordingly how to adjust fixator can also be predicted when there are many deviations; the adjustment value of fixator is the sum of adjustment value of every pin deviation.


## 3.3. Experimental Verification

The experiment of fixator-bone system for deformity correction at ideal situation shown in Figure8 was established, in order to better verify the analytical results. In experiment, the bone deformation has 20° internal rotation, 20° external rotation, 15° angular deformity, and 5mm displacement deviations on three planes, which is the same as the clinical case. The adjustment values of fixator joints which need to correct deformity are also in accordance with clinical practice, meaning the experiment is reasonable. Moreover, seven situations were also established to verify the analytical results; experimental results are shown in Table 4. The parameters of fixator-bone system in experiment are consistent with analytical model in Figure 1. In the process of experiment, one side of bone was fixed; by adjusting fixator joints, the deformity can be completely corrected.Table 4
The adjustment values of fixator joints at seven situations for experimental results.
Four situations a 5(degree) a 6(degree) d(mm) a 8(degree) a 9(degree) a 10(degree) Ideal situation -5.5 20 9 -10 -3 -12 φ 1=20° 0.5 8.5 3 4 -9.5 -9 φ 2=20° 16 21 -8 -16.5 -6.5 -20.5 φ 3=20° 5.5 -5 7.5 -14 2.5 -9.3 s 1=10mm -5 18.5 18.5 -11.5 -1.2 -10.5 s 2=10mm -5.6 15 6.0 -12 2 -8.5 s 3=10mm -12 20 10.5 -6.5 -4.5 -5Figure 8
The experiment of fixator-bone system for bone deformity at ideal situation. (a) The experiment of fixator-bone system before deformity correction. (b) The experiment of fixator-bone system after deformity correction.
(a) (b)Table4 shows the adjustment values of fixator joints at seven situations in experiment. Take the ideal situation as reference, when compared with pin deviations of φ2 and φ3, the adjustment value of the fixator joint of a8 has the biggest discrepancy between the ideal situation and the existence of pin deviation φ1. The presence of φ2 causes a greater change than φ1 and φ3 in the amount of a5, a8, and a10. Compared with φ1 and φ2, φ3 has the biggest effect on a6 and a9. In summary, when compared with ideal situation, the existence of pin deviation will affect the regulating amount of fixator joint, and φ2 will affect more joints and a greater amount of adjustment is needed; φ1 has only bigger influences on one joint of a8 and φ3 causes bigger effect on two joints of a6 and a9. That also indicates that φ2 has the bigger sensitivity than φ1 and φ3, meaning φ2 has the biggest effect on the fixator joint adjustment. In addition, φ1 has the minimum sensitivity compared to φ2 and φ3. These experimental results are consistent with analytical results for angular deviations.For positional deviations, take the ideal situation with no pin deviations as a reference; compared withs2 and s3, the existence of s1 makes a bigger change in fixator joint d. Meanwhile, s2 has a great effect on a9 and a6. The existence of s3 makes a great contribution to the bigger changes on a5, a8, and a10 than s1 and s2. These results indicate that s3 has the biggest sensitivity on the fixator joint adjustment and s1 has minimum sensitivity, which can also verify the correctness of the analytical results.Table5 shows the discrepancy between experiment results and analytical results for fixator joint adjustment value. The errors between experiment results and analytic results are very small and within twenty percent, meaning that the experimental results can verify the correctness of the analysis results.Table 5
The errors between experimental and analytical results for fixator joint adjustment values.
Four situations a 5(degree) a 6(degree) d(mm) a 8(degree) a 9(degree) a 10(degree) Ideal situation 8% 3.8% 3.78% 5.2% 8.27% 1.05% φ 1=20° 16% 2.96% 13.4% 7.69% 1.56% 1.26% φ 2=20° 4.7% 0.213% 3.12% 1.17% 1% 1.7% φ 3=20° 1.7% 1.82% 5.4% 3.05% 0.56% 19.19% s 1=10mm 7.4% 2.32% 1.14% 5.41% 9.7% 0.52% s 2=10mm 5% 1.51% 4.11% 9.33% 3.7% 5.357% s3=10mm 0.35% 2.43% 3% 19.7% 3.17% 0.85%


## 4. Discussion

In this paper, the sensitivity analyses of pin deviation on fracture dislocation and the influence of pin deviation on fixator joint adjustment were presented; these researches provide important information and rules to doctor and researchers.For sensitivity analysis, firstly, when there is only one kind of pin deviation, pin deviation having the maximum sensitivity can be found. Althoughφ1, φ2, φ3 are all the orientation deviation and s1, s2, s3 are all the position deviation, the sensitivity are not the same. According to the value of the maximum sensitivity, for orientation deviation, it can be judged that the sensitivity of φ2 is greater than φ3 and φ1; the sensitivity of φ3 is larger than that of φ1. For position deviation, s2 has the highest sensitivity and s1 has the least sensitivity compared with s2. Comparing orientation deviation with position deviation, the sensitivity of orientation deviation is relatively large. This phenomenon explained in previous study that angle deformation was difficult to correct which need to regulate more fixator joints [19]. Secondly, when there were many kinds of deviations, the location of the maximum sensitivity can be found by seeking the maximum sensitivity in single deviation. This discovery will reduce a lot time of computation, and as long as we know the position of the maximum sensitivity of the single deviation, the location of maximum sensitivity can also be found for the multiple deviations.Based on clinical case, the effects of pin deviation on fixator joint are also investigated. Firstly, for orientation deviation,φ1 have a great influence on fixator joint for a8. φ2 have a great influence on a5, d, and a10. φ3 have a great influence on a6 and a9. Through the comparison, the impact of φ1, φ2, and φ3 on fixator joints, due to φ2, affects more joints with more adjustment value and from this point of view, in which it can be drawn that φ2 has the biggest influence on the adjustment of fixator; the effect of φ1 is minimal. Meanwhile, for position deviation, s3 has a great influence on a5, a8, and a10. Due to the existence of s2, the regulatory value of a9 and a6 changed a lot. s1 made a great change in d. Comparing s2, s1, and s3, s3 can affect more joints, meaning when the deviation of s3 occurs, it is required to adjust more joints with greater amount of regulation in order to correct fracture completely; s3 has greater influence on the adjustment of fixator joint than s1 and s2. The impact of s1 is minimal. The results of the effects of pin deviation on fixator joints verify the sensitivity analysis from a point of view. The greater the sensitivity of pin deviation, the more the joints affected. Another interesting finding is that the impact of pin deviation on fixator joints has the property of superposition; that is, if any deviations exist at the same time, that will have an impact on more fixator joints, which are equal to the sum of the effect of each deviation on fixator joint adjustment. Meanwhile, based on the clinical case, seven situations of pin deviation in experimental were set up to verify the analytic results. Moreover, in clinical practice, at the ideal situation, the adjustment values of fixator joints are consistent with analytical results and experimental resultsThrough these analytic results from these findings, the orthopedist according to the actual situations reasonably increases or decreases the adjustment value of fixator joints, which could help correct the fracture residual dislocation caused by pin deviation; it also helps in reducing anesthesia time, lowering patient morbidity, and potentially lowering the risk of infection and radiation exposure, which also helps medical personnel to flexibly correct the fracture. The decrease of joint regulation due to the existence of pin deviation will increase the joint solution domain so as to expand the application scope of various bone deformities and avoid reaching the limit value of joints. Without knowing these information, complex reduction pattern could only be determined through trial and error, and pin needs to be inserted into bone more than once, which could prolong operation time and duplicate drilling attempts which may leave more bone tunnels that could potentially affect ligamentous reconstructions or complex fracture reduction [20].The analysis technique and results reported herewith can provide useful information for the medical personnel who use external fixators in the management of long bone fractures and for device manufacturers to enhance their products’ efficacy in adjustability, modify their current fixators, or develop new devices. These results should also be considered for possible manual or motorized adjustment in the next generation of external fixator.


## 5. Conclusions

Different pin deviation has different sensitivity; the greater the sensitivity of pin deviation, the greater the effect on the accuracy of fracture correction and the adjustment of fixator. The regulatory analysis and sensitivity analysis of pin deviation are validated by each other, based on a clinical case. The greater sensitivity of pin deviation means that it may be more likely to cause fracture dislocation and then need to adjust more joints. Moreover, the experiments are established to verify the analytic results. This research could provide the important theoretical foundation for orthopedics and these are crucial steps toward its successful applications.



---

*Source: 9267570-2018-10-22.xml*