 Research article
 Open Access
 Published:
Novel implant design of the proximal interphalangeal joint using an optimized rolling contact joint mechanism
Journal of Orthopaedic Surgery and Researchvolume 14, Article number: 212 (2019)
Abstract
Background
The aims of this study were to propose a novel implant design for the proximal interphalangeal joint (PIPJ) of the hand using a rolling contact joint (RCJ) mechanism and to derive an optimal implant design based on human PIPJ kinematics.
Methods
In total, 10 participants with normal PIPJs were enrolled in this study. True lateral finger radiographs were obtained in 10° increments from 0º (full extension) to 120° flexion of PIPJ. Radiographs were used to determine the average center of rotation, which formed the basis of a mathematical expression of the PIPJ kinematics. The variations in extensor tendon excursions in relation to the range of motion of PIPJ were determined using results from previous cadaveric studies. As the next step, a PIPJ implant design using an RCJ mechanism that was most consistent with the mathematically expressed PIPJ kinematics and tendon excursions was determined using a constrained optimization algorithm.
Results
The final proposed PIPJ implant had a relatively constant center of rotation over the entire PIPJ range of motion among the participants. In addition, the extensor tendon excursions of the proposed implant as applied to the phalangeal bones were similar to those of the human tendon. The proposed PIPJ implant achieved an acceptable position of the RCJ surface on the proximal and middle phalanges, which was derived from the constrained optimization algorithm.
Conclusions
A novel PIPJ implant design using an RCJ mechanism demonstrated acceptable outcomes in terms of PIPJ human kinematics and tendon excursions.
Introduction
Joint replacement arthroplasty of the proximal interphalangeal joint (PIPJ) can be performed as the salvage procedures when the joint is destroyed due to reasons such as degenerative arthritis, inflammatory arthritis, and posttraumatic arthritis [1]. In particular, posttraumatic arthritis is a wellknown complication of articular fracture associated with PIPJ dislocation, and total PIPJ replacement arthroplasty might be one of the best options for reducing the pain and rangeofmotion limitation in patients with posttraumatic arthritis due to the PIPJ dislocation [2, 3]. Several studies have been conducted to develop implants used in PIPJ replacement arthroplasty.
Various implant designs for the PIPJ have been proposed in clinical practice. In particular, the following two types of designs have been suggested: the constrained design, that has a hinge structure, and the unconstrained design, that uses surface replacement [4]. An optimal combination of designs and materials to provide durability and biocompatibility to the PIPJ implant has not yet been established [5]. Moreover, the PIPJ implants that have been developed thus far are less durable and are able to achieve less improvement in the functional range of motion than hip or knee joint prostheses [6].
The field of robotic arm applications has recently contributed to the development of implant surfaces using a rolling contact joint (RCJ) mechanism [7, 8]. The RCJ mechanism is defined as a type of restraint system that allows two circular components, connected with flexible straps, to rolling in relation to other contact surface without slip [9]. An implant using the RCJ mechanism offers the advantage of a greater range of motion for the joint, and a lower frictional force that reduces the risk of a chronic inflammatory response to wear debris compared with other implant types [10, 11]. In addition, the RCJ mechanism can help achieve greater stability because of its constrained type design compared to implants using unconstrained designs [12].
Design optimization involves choosing parameters for optimal design from among many alternatives. It requires a cost function that measures the goodness of the design. The cost function is a function of the design parameters and is designed to be minimized when the design is optimal. As per the design optimizations, constraints about the parameters can be needed. We focused on the nonlinear constrained optimization in the present study, because the design optimization used in this study is nonlinear with some constraints. There are many methods to achieve nonlinear constrained optimization, such as Lagrange multiplier, interiorpoint method, and sequential quadratic programming [13, 14]. The Lagrange multiplier is a method for finding the local minima of a function subject to equality constraints. Both the interiorpoint method and sequential quadratic programming use the Lagrange multiplier to solve a nonlinear constrained optimization.
The implant design for the upper extremity joints using the RCJ mechanism was first introduced by Bora [15]. However, this design was not tailored to the PIPJ and did not consider the size and kinematics of this joint. Therefore, the present study aimed to propose a novel implant design of the hand PIPJ using an RCJ mechanism and derive an optimal implant design based on the human PIPJ kinematics determined using plain radiographs.
Materials and methods
Study participants and acquisition of plain radiographic data
In total, 10 patients (8 female and 2 male) who underwent open reduction and internal fixation for unilateral distal radius fractures, were recruited from December 2018 to March 2019. The inclusion criteria included patients who were aged between 20 and 59 years, who had no history of hand trauma, and who had no history of medication use that could affect bone metabolism. Patients with hand osteoarthritis, including that affecting the PIPJ, or with other upper extremity musculoskeletal conditions were excluded from this study.
We used a true lateral plain radiograph to evaluate the kinematics of the proximal and middle phalanges in relation to their ranges of motion. All images were obtained during surgery after fixation of the distal radius fracture with PIPJ of the second finger positioned at the center of the image. For each participant, 10–12 images were obtained across the PIPJ range of motion from 0 (full extension) to 120° flexion, at intervals of approximately 10°. True lateral plain radiographs of the finger were defined as images in which the two condyles of the proximal phalangeal head overlapped.
Analysis and optimization of data on human PIPJ kinematics and PIPJ implant design using the RCJ mechanism were performed using Matlab® software (version R2018a, Mathworks, Inc., Natick, MA, USA). The NX computeraided design (CAD) program (version 8.5, Siemens, Munich, Germany) was used to formulate the positions of the proximal and middle phalanges in a twodimensional (2D) plane based on digitized plain radiographic data.
The study protocol is compliant with the Declaration of Helsinki and was reviewed and approved by the Institutional Review Board of the University Hospital in November 2018 (B1810/497002). Written informed consent was obtained from all participants.
Configuration and parameters of the proposed PIPJ implant
Hillberry’s RCJ implant designs, which has been applied to knee joint prostheses, were used [16, 17]. This design allows motion around 1 degree of freedom (1DOF), comprising PIPJ flexion or extension in the sagittal plane and includes two components: the middle phalangeal component (MPC) and proximal phalangeal component (PPC). Each component has one stem and one head with a circular joint surface; the components are linked by three flexible straps of equal width (Fig. 1a). Among them, two straps are symmetrically located in relation to the third strap, which is located at the center. On the basis of a previous cadaveric study of PIPJ morphology, the width of the radioulnar heads of both components was set a 10 mm and of the straps at 2.5 mm [18].
To determine the mounting positions of PPC and MPC, 2D coordinate frames Σ_{p} and Σ_{m} are placed at the proximal and middle phalangeal bones, respectively. The radii of the rolling contact surfaces (RCSs) of the PPC and MPC heads are denoted by r_{p} and r_{m}, respectively. The positions of the centers of the two circles represented in Σ_{p} and Σ_{m} are denoted by \( {}^p{p}_p \) and \( {}^m{p}_m \), respectively. The two RCSs can be assembled with a nonzero offset angle ϕ (Fig. 1b). The shape of the dorsal part of each implant head is important for determining the magnitude of extensor tendon excursion. These two circular surfaces can be defined by their radii r_{pe} and r_{me} and the center positions \( {}^p{p}_{pe} \) and \( {}^m{p}_{me} \). The proposed PIPJ implant has 13 design parameters including those related to rolling contact motion (r_{p}, x_{p}, y_{p}, r_{m}, x_{m}, y_{m}, ϕ) and extensor excursion (r_{pe}, x_{pe}, y_{pe}, r_{me}, x_{me}, y_{me}), where (x_{p}, y_{p}), (x_{m}, y_{m}), (x_{pe}, y_{pe}), and (x_{me}, y_{me}) represent the components of the position vectors of \( {}^p{p}_p \), \( {}^p{p}_m \), \( {}^p{p}_{pe} \), and \( {}^p{p}_{me} \), respectively.
The relationship between two coordinate frames can be expressed by a rotation matrix and a position vector of the origin of a target frame with respect to a base frame. The rotation matrix \( {}^p{R}_m \) and the position vector \( {}^p{o}_m \) of Σ_{m} with respect to Σ_{p} are as follows:
To simplify reporting, the flexion angle θ was defined in the clockwise direction and the minus signs at θ were applied as shown in Eq. (1). Using this equation, an arbitrary position vector in Σ_{m} can be transformed to a position vector in Σ_{p} as follows:
where the 2 × 1 vectors \( {}^pp={\left[{x}_p,{y}_p\right]}^T \) and \( {}^mp={\left[{x}_m,{y}_m\right]}^T \) indicate the same position with respect to the different coordinate frames and each vector consists of two components. For a compact representation, Eq. (2) can be simplified to the homogeneous transformation (HT) matrix as follows:
where \( {}^p{T}_m \) denotes the 3 × 3 HT matrix of Σ_{p} with respect to Σ_{m}. The notations \( \overset{\sim }{{}^pp}={\left[{x}_p,{y}_p,1\right]}^T \) and \( \overset{\sim }{{}^mp}={\left[{x}_m,{y}_m,1\right]}^T \) represent 3 × 1 vectors with the same components as \( {}^pp \) and \( {}^mp \), respectively, and the third unit component.
Using the design parameters, the HT matrix from Σ_{p} to Σ_{m} can be calculated. As illustrated in the Fig. 1b, given the flexion angle θ, the angle θ_{r} from \( {}^p{p}_p \) to \( {}^p{p}_m \) as follows:
Thus, using θ and θ_{r}, the position vector \( {}^p{o}_m \) of Σ_{m} with respect to Σ_{p} can be rewritten as follows:
Therefore, the HT matrix from Σ_{p} to Σ_{m} is as follows:
This HT matrix indicates the relative motion between the proximal and middle phalangeal bones according to the PIPJ flexion angle when the proposed implant is applied. This relative motion is compared with the human kinematics data to derive the optimal design parameters.
Determination of the average center of rotation
We hypothesized that the human PIPJ is a hinge joint with a 1DOF of motion around the center of rotation in the sagittal plane, and we proposed an estimation method that gives the center of rotation of human PIPJ, i.e., average center of rotation (ACR), using plain radiographic data to explain the ideal human PIPJ kinematics.
Using the total number (n) of radiographic images, the corresponding n number of HT matrices relating the proximal and middle phalanges are calculated. The contours of the proximal and middle phalangeal bones are selected using the CAD program, and the coordinate frames are attached at the centers of the contours (Fig. 2a). The long axis of the proximal and middle phalangeal bones are defined as x axes, with the axis perpendicular to the x axis defined as y. Σ_{b}, Σ_{pi}, and Σ_{mi} indicate the base frame fixed at an absolute position in the image, the frame of the proximal phalanx of ith radiographic image, and the frame of the middle phalanx of ith radiographic image, respectively (Fig. 2b); each of these three are set in the same orientation as that of the first image of PIPJ in full extension (Fig. 2a). HTs for the ith image, \( {}^b{T}_{\mathrm{p}i} \) and \( {}^b{T}_{mi} \), are obtained using the CAD program (Fig. 2b). For each image, \( {}^p{T}_{mi} \), which represents the relative position and orientation between the proximal and middle phalanges, is calculated using a simple matrix operation as follows:
The ACR position is calculated using the acquired n numbers of HT matrices assuming that \( {}^p{c}_p=\left({x}_{pc},{y}_{pc}\right) \) and \( {}^m{c}_m=\left({x}_{mc},{y}_{mc}\right) \) are the ACR positions representing the true center of rotation on the proximal and middle phalanges, respectively. \( {}^p{c}_p \) is fixed in Σ_{p}; however, \( {}^p{c}_{mi} \), which is the position of \( {}^m{c}_m \) with respect to Σ_{p}, varies according to the images. In order for these positions to represent the optimal position of the center of rotation, the following least square distance error should be minimized (Fig. 2c):
The optimal solution for Eq. (8) is achieved using a pseudoinverse matrix. The error around the ith image is expressed using the following HT:
This is arranged to the following desired equation:
where A_{i} and B_{i} represent a 2 × 4 matrix and a 2 × 1 vector, respectively, and X = [x_{mc} y_{mc} x_{pc} y_{pc} ]^{T} is a vector containing the ACR positions for optimizing. By combining Eq. (10) for all n of HTs, a resultant equation, for which the error should be minimized, is achieved as follows:
where the sizes of A and B are 2n × 4 and 2n × 1 matrices, respectively. Using a pseudoinverse matrix, the solution that is used to minimize the error of Eq. (8) is obtained as follows:
where A^{+} denotes the pseudoinverse matrix of A and, if A is fullrank, A^{+} is equivalent to (A^{T}A)^{−1}A^{T} and the optimal solution \( \hat{X} \) contains the ACRs, that is \( \hat{X}={\left[{\hat{x}}_{mc}{\hat{y}}_{mc}{\hat{x}}_{pc}{\hat{y}}_{pc}\right]}^T \). The position \( \left({\hat{x}}_{pc},{\hat{y}}_{pc}\right) \) best represents the center of rotation with a minimum error when the proximal phalangeal bone is considered as a fixed object; similarly, \( \left({\hat{x}}_{mc},{\hat{y}}_{mc}\right) \) is the corresponding position when the middle phalangeal bone is considered as a fixed object. Since the human PIPJ was hypothesized to be a hinge joint, its HT could be obtained based on the ACR derived above, as follows:
Magnitude of tendon excursions in accordance with the PIPJ flexion angles
The magnitude of extensor tendon excursions was determined in accordance with the PIPJ flexion angles suggested in previous cadaveric studies [19]. The effects of PIPJ motion on the flexor and extensor tendon excursions were evaluated under condition in which the motion of other upper extremity joints was restricted. The excursion of the extensor tendon based on the extensor digitorum communis increased by 0.08 mm with each 1° increase in PIPJ forward flexion [19]. Therefore, the magnitude of human PIPJ excursion can be modeled as follows:
The magnitude of extensor tendon excursion according to the flexion angles of the proposed PIPJ implant were determined using the design parameters. As mentioned above, the two circular surfaces on the dorsal part of each component of the implant head affect the magnitude of extensor tendon excursion and are illustrated in Fig. 3. Thus, the length of the extensor tendon with respect to θ can be denoted as follows:
with \( {\phi}_e=\mathrm{atan}2\left({L}_y,\kern0.5em {L}_x\right)\mathrm{asin}\left(\frac{r_{pe}{r}_{me}}{L}\right) \),
where L is the distance between \( {}^p{p}_{pe} \) and \( {}^p{p}_{me} \); L_{x} and L_{y} are the lengths of L projected to the x and y axes, respectively; ϕ_{e} is the clockwise angle between the x axis and extensor tendon tangential to the two circles; and L_{RCJ}(θ) is the length of the extensor tendon of the proposed PIPJ implant. We defined ∆L_{RCJ}(θ) as the difference between the length of the proposed extensor tendon at the flexion angle θ and the length at 0° (full extension), calculated as follows:
This implies the magnitude of extensor tendon excursion of the proposed PIPJ implant. Therefore, it was essential to properly determine the design parameters so that the function of Eq. (16) could become as consistent as possible with Eq. (14).
Determination of RCJ design parameters using optimization
We adopted a constrained optimization algorithm using sequential quadratic programming [20] to obtain the optimal RCJ design parameters. Radiographic images were obtained during PIPJ motion from 0 (full extension) to 120° flexion, which reflects the normal range of motion in humans. The proposed PIPJ implant was designed within the same range of motion. The optimization method used in this study comprised one cost function with two components and eight constraints.
The first component of the cost function was designed to allow RCJ to have a similar motion as the human PIPJ. The motions of the human PIPJ and the proposed PIPJ implant were expressed using the HT matrices obtained in Eqs. (6) and (13). The distance error between the human PIPJ and the proposed PIPJ implant motion can be calculated using the position difference between the origin of the two HT matrices. The mean squared error of the distance between the human PIPJ and the proposed PIPJ implant motion is obtained as follows:
where the first and second terms in the summation represent the origins of Σ_{m} in the human PIPJ and the proposed PIPJ implant, respectively. In addition, ‖∙‖^{2} denotes the square of the Euclidean norm of the vector, where the third component of the vector is ignored. The total n of errors in the distance from 0 to 120° of the flexion angle θ were averaged. In this optimization, the distance error at every 1° was considered; therefore, n was set to 121.
The second component of the cost function was designed to minimize the difference between the extensor tendon excursions of the human PIPJ and the proposed PIPJ implant. Using Eqs. (14) and (16), the mean squared error was calculated as follows:
Therefore, the overall cost function was defined as
where k_{p} and k_{E} are the weighted coefficients of C_{p} and C_{E}, respectively; these determine the priority and weight between the joint motion and excursion length optimizations. In our optimization, k_{p} and k_{E} were given equivalent weightings and set at 1.
In order for all parameters to represent the optimal PIPJ implant using the RCJ mechanism, C should be minimized as follows:
To obtain a feasible RCJ implant design, eight constraints were applied, while minimizing Eq. (20) as follows:
r_{p}≥ 5, r_{m}≥ 5 (23)
r_{pe}≥ 2, r_{me}≥ 2 (24)
Figure 4 illustrates the positions and angles related to these constraints. The notations (x_{de}, y_{de}) and (x_{ve}, y_{ve}) represent the positions of the dorsal and volar end points, respectively, on the RCS of PPC; θ_{d} and θ_{v} represent the angles of the tangential lines of these dorsal and volar end points in relation to the x axis, respectively; and (x_{I}, y_{I}) and (x_{V}, y_{V}) represent the positions of the dorsal infliction point and volar vertex of the proximal phalangeal head, respectively. Equation (21) and (22) place the RCS of PPC at a suitable position for surgery and indicate that the tangential lines of the dorsal and volar end points of the RCS are located distal to the dorsal infliction point (x_{I}, y_{I}) and volar vertex (x_{V}, y_{V}), respectively. On the basis of Eq. (23), the radii of each RCS were set at ≥ 5 mm. These constraints allowed the 2.5mm flexible strap to endure sufficient contact stress (Hz) [21]. Equation (24) restricts the radii of each circle at the dorsal side of the proposed PIPJ implant hand to ≥ 2 mm in order to prevent the extensor tendon from folding and kinking during flexion. Equations (25) and (26) represent the tangential positions of the dorsal circles of the proposed PIPJ implant head in relation to the dorsal end point (x_{de}, y_{de}).
Results
Position of the ACR
To represent the standard position of the ACR for each participant, a virtual circle was defined, which circumscribed the triangle connecting the dorsal infliction point, volar vertex, and distal vertex of the proximal phalangeal head. The center of the circle was denoted as C_{PH} (Fig. 5), and the relative positions of ACRs were expressed on the x and y axes (Fig. 6). The proximal ACRs were located at a mean of 0.04 ± 0.32 mm (range − 0.34 to 0.59 mm) and − 0.21 ± 0.46 mm (range − 1.37 to 0.43 mm) from C_{PH} on the x and y axes, respectively. The position of \( {}^m{c}_m \) on the x and y axes, which was centered on \( {}^p{c}_p \) in accordance with the PIPJ range of motion, did not vary with the range of motion and was located within 1.27 mm of \( {}^p{c}_p \) (Fig. 7). In addition, the proximal ACRs for each participant tended to be positioned inferior to C_{PH} (Fig. 8).
Determination of the RCJ design parameters
The mean values of r_{p}, r_{m}, r_{pe}, and r_{me} determined using the constrained optimization algorithm were 11.49 ± 1.53 mm (range 8.25–13.27 mm), 5.00 ± 0.00 mm, 2.15 ± 0.28 mm (range 2.00–2.87 mm), and 2.03 ± 0.09 mm (range 2.00–2.30 mm), respectively. Table 1 shows the long–axis length of the proximal and middle phalanges, anteroposterior diameters of the proximal phalangeal head and middle phalangeal base, and values of parameters determined for each participant.
The mean value of the position errors between the centers of rotation of the PIPJ implant derived from the optimized algorithm and ACRs of the human PIPJ determined using plain radiographs was 1.25 ± 0.32 mm (range 0.63–1.81 mm), and the mean difference between the extensor tendon excursions of the proposed PIPJ implant and those of the human PIPJ was 0.16 ± 0.04 mm (range 0.1–0.24 mm). Both mean values were obtained for each of the 10 participants and were calculated as the mean of values measured from 0 (full extension) to 120° flexion in 1° increments. Therefore, the proposed PIPJ implant exhibited an acceptable kinematic range compared with that of the human PIPJ.
Using the determined parameters, the proposed PIPJ implant was designed using Matlab® (Fig. 9); a concordant threedimensional (3D) model was created using the CAD program (Fig. 10).
Discussion
The most essential criterion for developing a PIPJ implant using the RCJ mechanism is to achieve a joint design that precisely simulates human biomechanical properties. The RCJ motion should be similar to the human PIPJ kinematics, and the extensor and flexor tendon excursion should be within acceptable magnitudes as compared to those of the human PIPJ [22]. Accurate motion analysis is essential for reducing the human PIPJ kinematics. The initial step of our analysis involved determination of the PIPJ center of rotation.
Various methods to determine the center of rotation have been developed, which are based on calculating the finite displacement of a rigid body [23, 24]; however, this finite center of rotation becomes inaccurate if the displacement approaches zero. Thus, iterative numerical optimizations can be used to achieve the true center of rotation [25]. ACR proposed in this study provides a simple and efficient closedform solution of ACR in the sagittal plane using data from multiple radiographic images to minimize the distance error from the optimal center of rotation.
In addition to the flexion–extension in the sagittal plane, PIPJ kinematics includes abduction–adduction, axial rotation, and sagittal rotation [26]. However, the proposed PIPJ implant was based on Hillberry’s RCJ design, which reproduced only 1DOF of motion, namely flexion–extension, in sagittal plane. Abduction–adduction or rotational motion represents a small proportion of the total PIPJ motion [25], and flexion–extension in the sagittal plane better reflects the motions performed in daily life [27]. Among the PIPJ implant designs reported to date, a silicone implant with a volar approach has been found to exhibit the best durability and least complications [28]; this implant was a 1DOF hinge type joint allowing flexion–extension in the sagittal plane, consistent with our proposed PIPJ implant using an RCJ mechanism. Recently, surface replacement arthroplasties with ≥ 2DOF have shown poorer clinical outcomes than silicone implants [28]. Thus, our PIPJ implant with 1DOF has a suitable design and can reproduce human PIPJ kinematics.
The present study proposed a constrained type of PIPJ implant with 1DOF characterized by a hinge joint. Previous cadaveric motion analysis study has reported that the displacement of the center of rotation in the sagittal plane through the PIPJ range of motion was ≤ 0.30 mm [25]. ACRs in this study showed position errors within 1.27 mm, which were not correlated with the PIPJ flexion angle at a relatively constant position (Fig. 7 ). Therefore, the proposed PIPJ implant with a 1DOF hinge joint positioned on ACR appears to be appropriate.
A constrained optimization algorithm was applied to determine the design parameters of the proposed PIPJ implant. The proximal position of the RCJ surface, in relation to the neck of the proximal phalanx, was detected without including the condition for the RCJ surface position in the optimization algorithm. A joint rotation occurs on the RCJ surface in implants using an RCJ mechanism [17]; therefore, the RCJ surface of the proposed PIPJ implant was optimized to a position close to the center of rotation derived from human PIPJ kinematics. However, removal of a substantial amount of proximal phalanx bone stock is necessary during surgery to achieve the proximal position of the RCJ surface proximal to the neck of the proximal phalanx. In addition, PIPJ could become unstable as a result of injury to the collateral ligament during implantation [26]. Therefore, we included a condition in the constrained optimization algorithm requiring the tangential line extending from the dorsal and volar sides of the RCJ surface to be positioned distal to the dorsal infliction point and volar vertex of the proximal phalangeal head. Consequently, the optimized RCJ surface was located inside the head of proximal phalanx, and it could preserve a considerable amount of bone stock from the proximal phalangeal head and avoid damage to the collateral ligament during surgery, thereby contributing to improved postoperative PIPJ stability.
The interiorpoint method is an effect and rapid method for a large number of parameters and constraints. However, its solution can be relatively less accurate than those of other algorithms because the internal process of the algorithm keeps iterates away from inequality constraint boundaries. The sequential quadratic programming is appropriate for small or mediumscale optimization; however, it is relatively slow during the optimization of a large number of parameters and constraints. In this study, there were 13 parameters and eight constraints, and the sequential quadratic programming was adopted to solve the optimization.
The optimized results demonstrated that the diameter of the RCS of MPC was 5 mm and that of PPC was 8–13 mm; these findings can be explained by two factors. First, because the human PIPJ is hinged around the center of rotation, the optimized RCS diameter should be minimized in order to mimic human kinematics. However, to adequately endure the contact stress, the RCS diameter should be set at ≥ 5 mm; therefore, the optimized diameters of MPC and PPC should approach this value. Second, ACR determined using plain radiographic data was located near the neck of the proximal phalanx. As the RCJ center of rotation was located on RCS, a larger PPC diameter was required to enable RCS to be located near the neck of the proximal phalanx, thus, resembling human PIPJ kinematics.
Decreased length of the phalangeal bone may contribute to relative extensor tendon elongation. A previous cadaveric study has suggested that a decrease of 1 mm could lead to a 12° extension lag [29], that could cause functional limitation of the hand. In this study, equal weightings were given to the cost function components of the extensor tendon and the center of rotation to reduce the potential for bone–tendon length discrepancy following PIPJ implant insertion; thus, optimization was performed to maintain the long–axis phalangeal bone length. The final implant design demonstrated that the proposed PIPJ implant was on average 1.1 mm shorter than the sum of the lengths of the normal proximal and middle phalanges. However, given the compensation provided by the extensor tendon excursion reserve effect and the intrinsic muscle function [30], the proposed PIPJ implant design could show a few degrees of extension lag on PIPJ that does not significantly affect the hand function.
This study has certain limitations. First, because we used 2D plain radiographs, 3D PIPJ motions such as abduction–adduction or rotation could not be evaluated. Second, although the size and position of the proposed PIPJ implant were designed to reproduce human PIPJ kinematics and anatomy, the various technical problems that might occur during total PIPJ replacement arthroplasty cannot be ruled out. Third, because the 3D structures were projected in a 2D plane, the original size and contour of the PIPJ structures could be distorted. Fourth, because the PIPJ range of motion and magnitude of extensor excursion may be dependent to each other, it may cause the error in the cost function. However, the errors calculated from the two cost functions were 1.25 ± 0.32 mm and 0.16 ± 0.04 mm, respectively; these were small enough to be clinically significant. Therefore, the feasibility of this study is maintained.
Conclusions
To our knowledge, this is the first study to design a PIPJ implant using the RCJ mechanism to fit the human PIPJ structure. Further, the concept of ACR was newly introduced to determine the center of rotation during the human PIPJ range of motion. A novel PIPJ implant design using the RCJ mechanism demonstrated acceptable features in terms of tendon excursions and human PIPJ kinematics. Future studies on the biomechanical properties of implant materials and strap applications may help accelerate the clinical application of the RCJ mechanism implants in total PIPJ arthroplasty.
Availability of data and materials
The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.
Abbreviations
 1DOF:

1 degree of freedom
 2D:

Twodimensional
 2DOF:

2 degrees of freedom
 3D:

Threedimensional
 ACR:

Average center of rotation
 CAD:

Computeraided design
 HT:

Homogenous transform
 MPC:

Middle phalangeal component
 PIPJ:

Proximal interphalangeal joint
 PPC:

Proximal phalangeal component
 RCJ:

Rolling contact joint
 RCS:

Rolling contact surface
References
 1.
Rizzo M, Beckenbaugh RD. Proximal interphalangeal joint arthroplasty. J Am Acad Orthop Surg. 2007;15(3):189–97.
 2.
Mangelson JJ, Stern PJ, Abzug JM, Chang J, Osterman AL. Complications following dislocations of the proximal interphalangeal joint. Instr Course Lect 2014;63:123–30.
 3.
Kamnerdnakta S, Huetteman HE, Chung KC. Complications of proximal interphalangeal joint injuries. prevention and treatment. Hand Clin. 2018;34(2):267–88.
 4.
Zhu AF, Rahgozar P, Chung KC. Advances in proximal interphalangeal joint arthroplasty: biomechanics and biomaterials. Hand Clin. 2018;34(2):185–94.
 5.
DeDeugd CM, Rizzo M. Surgical exposure of the proximal interphalangeal joint. Hand Clin. 2018;34(2):127–38.
 6.
Aversano FJ, Calfee RP. Salvaging a failed proximal interphalangeal joint implant. Hand Clin. 2018;34(2):217–27.
 7.
Jeanneau Al, Herder J, Laliberté T, Gosselin Cm. A compliant rolling contact joint and its application in a 3DOF planar parallel mechanism with kinematic analysis. ASME 2004 International Design Engineering Technical Conferences and Computers and Information in Engineering Conference; 28 Sep  2 Oct; Salt Lake City: ASME; 2004. p. 689698.
 8.
Catalano MG, Grioli G, Farnioli E, Serio A, Piazza C, Bicchi A. Adaptive synergies for the design and control of the Pisa/IIT SoftHand. Int J Rob Res. 2014;33(5):768–82.
 9.
Nelson TG, Herder JL. Developable compliantaided rollingcontact mechanisms. Mech Mach Theo. 2018;126:225–42.
 10.
Kim S, In H, Song J, Cho K. Force characteristics of rolling contact joint for compact structure. 2016 6th IEEE International Conference on Biomedical Robotics and Biomechatronics (BioRob); 2629 June. Singapore: IEEE; 2016. p. 1207–12.
 11.
Bitar D, Parvizi J. Biological response to prosthetic debris. World J Orthop. 2015;6(2):172–89.
 12.
Slocum AH Jr, Cervantes TM, Seldin EB, Varanasi KK. Analysis and design of rollingcontact joints for evaluating bone plate performance. Med Eng Phys. 2012;34(7):1009–18.
 13.
Gill PE, Murray W, Saunders MA, Wright MH. Recent developments in constrained optimization. J Comput Appl Math. 1988;22(2):257–70.
 14.
Schmidt M. An interiorpoint method for nonlinear optimization problems with locatable and separable nonsmoothness. EURO J Comput Opt. 2015;3(4):309–48.
 15.
Bora FW, Jr, inventor; Bora FW, Jr., assignee. Prosthetic joint. US patent 4267608A. 1981.
 16.
Hillberry BM, Hall AS, Jr, inventors; Purdue Research Foundation, assignee. Rolling contact prosthetic knee joint. .US patent 3945053A. 1976.
 17.
Slocum AH Jr. Rolling contact orthopedic joint design [PhD]. Boston: Massachusettes Institue of Technology; 2013.
 18.
Lawrence T, Trail IA, Noble J. Morphological measurements of the proximal interphalangeal joint. J Hand Surg Br. 2004;29(3):244–9.
 19.
Elliot D, McGrouther DA. The excursions of the long extensor tendons of the hand. J Hand Surg Br. 1986;11(1):77–80.
 20.
Zhao W, Wang L, Yin Y, Wang B, Tang Y. Sequential quadratic programming enhanced backtracking search algorithm. Frontiers of Computer Science. 2018;12(2):316–30.
 21.
Sanders AP, Brannon RM. Assessment of the applicability of the Hertzian contact theory to edgeloaded prosthetic hip bearings. J Biomech. 2011;44(16):2802–8.
 22.
Elloy MA, Wright JTM, Cavendish ME. The basic requirements and design criteria for total joint prostheses. Acta Orthop Scand. 1976;47(2):193–202.
 23.
Challis JH. Estimation of the finite center of rotation in planar movements. Med Eng Phys. 2001;23(3):227–33.
 24.
McCane B, Abbott JH, King T. On calculating the finite centre of rotation for rigid planar motion. Med Eng Phys. 2005;27(1):75–9.
 25.
Hess F, Furnstahl P, Gallo LM, Schweizer A. 3D analysis of the proximal interphalangeal joint kinematics during flexion. Comput Math Methods Med. 2013;2013:138063.
 26.
Pang EQ, Yao J. Anatomy and biomechanics of the finger proximal interphalangeal joint. Hand Clin. 2018;34(2):121–6.
 27.
Bain GI, Polites N, Higgs BG, Heptinstall RJ, McGrath AM. The functional range of motion of the finger joints. J Hand Surg Eur Vol. 2015;40(4):406–11.
 28.
Yamamoto M, Chung KC. Implant arthroplasty: selection of exposure and implant. Hand Clin. 2018;34(2):195–205.
 29.
Markiewitz AD. Complications of hand fractures and their prevention. Hand Clin. 2013;29(4):601–20.
 30.
Vahey J, Wegner DA, Hastings H 3rd. Effect of proximal phalangeal fracture deformity on extensor tendon function. J Hand Surg Am. 1998;23(4):673–81.
Acknowledgements
None declared
Funding
This research was supported by the Seoul National University Bundang Hospital Research Fund (grant no. 142017001) and in part by the convergence technology development program for bionic arm through the National Research Foundation of Korea (NRF) funded by the Ministry of Science & ICT (no. 2015M3C1B2052817). The funders had no role in the study design, data collection, analysis and interpretation, and writing the manuscript.
Author information
Affiliations
Contributions
SWH, HSG, and YJK designed the study. SWH and JSY contributed to the acquisition of data. SWH, JSY, and YJK contributed to the analysis and interpretation of data. SWH, YJK, and HSG prepared and edited the manuscript. All authors read and approved the final manuscript.
Corresponding author
Correspondence to Hyun Sik Gong.
Ethics declarations
Ethics approval and consent to participate
This study protocol was reviewed and approved by the Institutional Review Board of the University Hospital in November 2018(B1810/497002) and written informed consent was obtained from all participants.
Consent for publication
Not applicable
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Proximal interphalangeal joint implant
 Rolling contact joint mechanism
 Average center of rotation
 Tendon excursion
 Constrained optimization