Fusion of radiostereometric analysis data into computed tomography space: application to the elbow joint

Improvement of joint prostheses is dependent upon information concerning the biomechanical properties of the joint. Radiostereometric analysis (RSA) and electromagnetic techniques have been applied in previous cadaver and in vivo studies on the elbow joint to provide valuable information concerning joint motion axes. However, such information is limited to mathematically calculated positions of the axes according to an orthogonal coordinate system and is difficult to relate to individual skeletal anatomy. The aim of this study was to evaluate the in vivo application of a new fusion method to provide three-dimensional (3D) visualization of flexion axes according to bony landmarks. In vivo RSA data of the elbow joint's flexion axes was combined with data obtained by 3D computed tomography (CT). Results were obtained from five healthy subjects after one was excluded due to an instable RSA marker. The median error between imported and transformed RSA marker coordinates and those obtained in the CT volume was 0.22 mm. Median maximal rotation error after transformation of the rigid RSA body to the CT volume was 0.003 degrees . Points of interception with a plane calculated in the RSA orthogonal coordinate system were imported into the CT volume, facilitating the 3D visualization of the flexion axes. This study demonstrates a successful fusion of RSA and CT data, without significant loss of RSA accuracy. The method could be used for relating individual motion axes to a 3D representation of relevant joint anatomy, thus providing important information for clinical applications such as the development of joint prostheses.     

MRI vs CT-based 2D-3D auto-registration accuracy for quantifying shoulder motion using biplane video-radiography

Biplane 2D-3D registration approaches have been used for measuring 3D, in vivo glenohumeral (GH) joint kinematics. Computed tomography (CT) has become the gold standard for reconstructing 3D bone models, as it provides high geometric accuracy and similar tissue contrast to video-radiography. Alternatively, magnetic resonance imaging (MRI) would not expose subjects to radiation and provides the ability to add cartilage and other soft tissues to the models. However, the accuracy of MRI-based 2D-3D registration for quantifying glenohumeral kinematics is unknown. We developed an automatic 2D-3D registration program that works with both CT- and MRI-based image volumes for quantifying joint motions. The purpose of this study was to use the proposed 2D-3D auto-registration algorithm to describe the humerus and scapula tracking accuracy of CT- and MRI-based registration relative to radiostereometric analysis (RSA) during dynamic biplanar video-radiography. The GH kinematic accuracy (RMS error) was 0.6–1.0 mm and 0.6–2.2° for the CT-based registration and 1.4–2.2 mm and 1.2–2.6° for MRI-based registration. Higher kinematic accuracy of CT-based registration was expected as MRI provides lower spatial resolution and bone contrast as compared to CT and suffers from spatial distortions. However, the MRI-based registration is within an acceptable accuracy for many clinical research questions.     

Accuracy of a CT-based bone contour registration method to measure relative bone motions in the hindfoot

For measuring the in-vivo range of motion of the hindfoot, a CT-based bone contour registration method (CT-BCM) was developed to determine the three-dimensional position and orientation of bones. To validate this technique, we hypothesized that the range of motion in the hindfoot is equally, accurately measured by roentgen stereophotogrammetric analysis (RSA) as by the CT-BCM technique.Tantalum bone markers were placed in the distal tibia, talus and calcaneus of one cadaver specimen. With a fixed lower leg, the cadaveric foot was held in neutral and subsequently loaded in eight extreme positions. Immediately after acquiring a CT-scan with the foot in a position, RSA radiographs were made. Bone contour registration and RSA was performed. Helical axis parameters were calculated for talocrural and subtalar joint motion from neutral to extreme positions and between opposite extreme positions. Differences between CT-BCM and RSA were calculated.Compared with RSA, the CT-BCM data registered an overall root mean square difference (RMSd) of 0.21° for rotation about the helical axis, and 0.20 mm translation along the helical axis for the talocrural and subtalar joint and for all motions combined. The RMSd of the position and direction of the helical axes was 3.3 mm and 2.4°, respectively. The latter errors were larger with smaller helical rotations.The differences are similar to those reported for validated RSA and thus are not clinically relevant. Concluding, CT-BCM is an accurate and accessible alternative for studying joint motion, as it does not have the risk of infection and overlapping bone markers.     

A joint coordinate system proposal for the study of the trapeziometacarpal joint kinematics

The International Society of Biomechanics (ISB) has recommended a standardisation for the motion reporting of almost all human joints. This study proposes an adaptation for the trapeziometacarpal joint.

The definition of the segment coordinate system of both trapezium and first metacarpal is based on functional anatomy. The definition of the joint coordinate system (JCS) is guided by the two degrees of freedom of the joint, i.e. flexion–extension about a trapezium axis and abduction–adduction about a first metacarpal axis. The rotations obtained using three methods are compared on the same data: the fixed axes sequence proposed by Cooney et al., the mobile axes sequence proposed by the ISB and our alternative mobile axes sequence. The rotation amplitudes show a difference of 9° in flexion–extension, 2° in abduction–adduction and 13° in internal–external rotation.

This study emphasizes the importance of adapting the JCS to the functional anatomy of each particular joint.     

Evaluation of knee functional calibration with and without the effect of soft tissue artefact

Functional calibration methods were devised to improve repeatability and accuracy of the knee flexion–extension axis, which is used to define the medio-lateral axis of the femur coordinate system in gait analysis. Repeatability of functional calibration methods has been studied extensively in healthy individuals, but not accuracy in the absence of a benchmark knee axis. We captured bi-plane fluoroscopy data of the knee joint in 17 subjects with unilateral total knee arthroplasty during treadmill walking. The prosthesis provided a benchmark knee axis to evaluate the functional calibration methods. Stereo-photogrammetry data of thigh and shank marker clusters were captured simultaneously to investigate the effect of soft tissue artefact (STA). Three methods were tested, the Axis Transformation Technique (ATT) finds the best single fixed axis of rotation, 2DofKnee finds the axis that minimises knee varus–valgus and trajAJC finds the axis perpendicular to the trajectory, in the transverse plane of the femur, of a point located on the longitudinal axis of the tibia. Using fluoroscopy data, functional axes formed an angle of less than 2° in the transverse plane with the benchmark axis. True internal–external range of movement was correlated with decreased accuracy for ATT, while varus–valgus range of movement was correlated with decreased accuracy for 2DofKnee and trajAJC. STA had negative impact on accuracy and variability. Using stereo-photogrammetry data, the accuracy of 2DofKnee was 1.7°(SD: 5.1°), smaller than ATT 2.9°(SD: 5.1°) but not to trajAJC 1.7°(SD: 5.2°). Our results confirm that of previous studies, which utilised the femur condylar axis as reference.     

A novel functional calibration method for real-time elbow joint angles estimation with magnetic-inertial sensors

Magnetic-inertial measurement units (MIMUs) are often used to measure the joint angles between two body segments. To obtain anatomically meaningful joint angles, each MIMU must be computationally aligned (i.e., calibrated) with the anatomical rotation axes. In this paper, a novel four-step functional calibration method is presented for the elbow joint, which relies on a two-degrees-of-freedom elbow model. In each step, subjects are asked to perform a simple task involving either one-dimensional motions around some anatomical axes or a static posture. The proposed method was implemented on a fully portable wearable system, which, after calibration, was capable of estimating the elbow joint angles in real time. Fifteen subjects participated in a multi-session experiment that was designed to assess accuracy, repeatability and robustness of the proposed method. When compared against an optical motion capture system (OMCS), the proposed wearable system showed an accuracy of about 4° along each degree of freedom. The proposed calibration method was tested against different MIMU mountings, multiple repetitions and non-strict observance of the calibration protocol and proved to be robust against these factors. Compared to previous works, the proposed method does not require the wearer to maintain specific arm postures while performing the calibration motions, and therefore it is more robust and better suited for real-world applications.     

A new non-orthogonal decomposition method to determine effective torques for three-dimensional joint rotation

This paper describes a new non-orthogonal decomposition method to determine effective torques for three-dimensional (3D) joint rotation. A rotation about a joint coordinate axis (e.g. shoulder internal/external rotation) cannot be explained only by the torque about the joint coordinate axis because the joint coordinate axes usually deviate from the principal axes of inertia of the entire kinematic chain distal to the joint. Instead of decomposing torques into three orthogonal joint coordinate axes, our new method decomposes torques into three "non-orthogonal effective axes" that are determined in such a way that a torque about each effective axis produces a joint rotation only about one of the joint coordinate axes. To demonstrate the validity of this new method, a simple internal/external rotation of the upper arm with the elbow flexed at 90 degrees was analyzed by both orthogonal and non-orthogonal decomposition methods. The results showed that only the non-orthogonal decomposition method could explain the cause-effect mechanism whereby three angular accelerations at the shoulder joint are produced by the gravity torque, resultant joint torque, and interaction torque. The proposed method would be helpful for biomechanics and motor control researchers to investigate the manner in which the central nervous system coordinates the gravity torque, resultant joint torque, and interaction torque to control 3D joint rotations.     

Shoulder joint kinetics and dynamics during underwater forward arm elevation

Aquatic exercises are widely implemented into rehabilitation programs. However, both evaluating their mechanical demands on the musculoskeletal system and designing protocols to provide progressive loading are difficult tasks. This study reports for the first time shoulder joint kinetics and dynamics during underwater forward arm elevation performed at speeds ranging from 22.5 to 90°/s. Net joint moments projected onto anatomical axes of rotation, joint power, and joint work were calculated in 18 participants through a novel approach coupling numerical fluid flow simulations and inverse dynamics. Joint dynamics was revealed from the 3D angle between the joint moment and angular velocity vectors, identifying three main functions—propulsion, stabilization, and resistance. Speeds <30°/s necessitated little to no power at all, whereas peaks about 0.20 W⋅kg⁻¹ were seen at 90°/s. As speed increased, peak moments were up to 61 × higher at 90 than at 22.5°/s, (1.82 ± 0.12%BW⋅AL vs 0.03 ± 0.01%BW⋅AL, P < 0.038). This was done at the expense of a substantial decrease in the joint moment contribution to joint stability though, which goes against the intuition that greater stabilization is required to protect the shoulder from increasing loads. Slow arm elevations (<30°/s) are advantageous for joint mobility gain at low mechanical solicitation, whereas the intensity at 90°/s is high enough to stimulate muscular endurance improvements. Simple predictive equations of shoulder mechanical loading are provided. They allow for easy design of progressive protocols, either for the postoperative shoulder or the conditioning of athlete targeting very specific intensity regions.     

Elbow joint biomechanics for preclinical evaluation of total elbow prostheses

Total elbow arthroplasty is a clinically successful procedure, yet long-term implant survival rates have historically lagged behind those reported for total hips and knees. Clinical complications associated with implant wear, osteolysis, stem loosening and device fracture have been implicated as reasons for limited long-term survivorship. Unfortunately, there is little published information on the biomechanics and method(s) for preclinical evaluation of total elbow prostheses that could provide insight into the mechanisms of failure. Additionally, there are no consensus testing standards or summaries of loading profiles of the humero-ulnar joint associated with a range of activities of daily living. Such data would facilitate the standardized preclinical assessment of total elbow devices such is commonplace for other large joints. The objective of the work here is therefore to provide a comprehensive review of elbow joint biomechanics as it relates to preclinical evaluation of total elbow implants. This summary includes a review of elbow joint forces, kinematics, the types and frequency of humero-ulnar joint motions associated with activities of daily living and clinical outcomes, as well as proposing a methodology for deriving humero-ulnar joint reaction force magnitudes and vector orientations as a function of a known mass/force at the hand. From these data, a scalable, bi-axial loading profile is proposed as a foundation for the development of clinically relevant, laboratory simulations for assessment of total elbow prostheses performance.     

Influence of muscle anatomical cross-sectional area on the moment arm length of the triceps brachii muscle at the elbow joint

The purpose of this study was to test the hypothesis that the musculotendon moment arm length is affected by the muscle anatomical cross-sectional area. The moment arm length of the triceps brachii (TB) muscle at 30°, 50°, 70°, 90°, 110° elbow flexion positions was measured in sagittal magnetic resonance images (MRI) of 18 subjects as the perpendicular distance between the center of the pulley of the humerus to the line through the center of the TB tendon. The moment arm increased as the elbow flexion angle decreased, from 1.74±0.13 cm at 110° to 2.39±0.14 cm at 30°. The maximal anatomical cross-sectional area of the TB muscle was significantly correlated with the moment arms at all joint positions (r=0.545–0.803, p<0.05). Furthermore, the circumference of the upper arm was also significantly correlated with the moment arms at all joint positions, except for 70° (r=0.504–0.702, p<0.05). These results indicate that the moment arm length of the TB muscle is affected by the muscle anatomical cross-sectional area.     

In-vivo range of motion of the subtalar joint using computed tomography

Understanding in vivo subtalar joint kinematics is important for evaluation of subtalar joint instability, the design of a subtalar prosthesis and for analysing surgical procedures of the ankle and hindfoot. No accurate data are available on the normal range of subtalar joint motion. The purpose of this study was to introduce a method that enables the quantification of the extremes of the range of motion of the subtalar joint in a loaded state using multidetector computed tomography (CT) imaging. In 20 subjects, an external load was applied to a footplate and forced the otherwise unconstrained foot in eight extreme positions. These extreme positions were foot dorsiflexion, plantarflexion, eversion, inversion and four extreme positions in between the before mentioned positions. CT images were acquired in a neutral foot position and each extreme position separately. After bone segmentation and contour matching of the CT data sets, the helical axes were determined for the motion of the calcaneus relative to the talus between four pairs of opposite extreme foot positions. The helical axis was represented in a coordinate system based on the geometric principal axes of the subjects’ talus. The greatest relative motion between the calcaneus and the talus was calculated for foot motion from extreme eversion to extreme inversion (mean rotation about the helical axis of 37.3±5.9°, mean translation of 2.3±1.1 mm). A consistent pattern of range of subtalar joint motion was found for motion of the foot with a considerable eversion and inversion component.     

Shoulder,elbow, and wrist joint angle excursions vary by gesture during touchscreen interaction

Repeated gesturing on touchscreen computing devices has become part of professional, personal, or school use by persons of all ages. Few studies have compared kinematics among joint motions and gestures during touchscreen interaction. We aimed to quantify the relative contributions of the shoulder, elbow and wrist to completion of several gestures to aid understanding of touchscreen ergonomics. Joint angles of the shoulder, elbow, and wrist were recorded for 22 seated participants while they interacted with a 10.1″ tablet computer held on an easel. Joint excursions at the shoulder, elbow, and wrist were all on average ≤20° during touchscreen interaction. The greatest excursion measured was shoulder rotation for swipe right with a mean of 15.5(±6.0)°. Index finger tap on a touchscreen was completed by participants with less than 5° of mean joint excursion at the shoulder, elbow and wrist. Tap, pinch and stretch gestures demonstrated significantly more wrist flexion/extension (p < 0.05) than shoulder flexion/extension, ab/adduction and rotation. Also, swipe left, right and up involved more shoulder rotation (p < 0.05) than wrist flexion/extension. These results suggest that when gestures are repeated frequently, the relative risk of overuse injury at the shoulder, elbow, or wrist may depend on the gesture being repeated.     

Prescribing joint co-ordinates during model preparation to improve inverse kinematic estimates of elbow joint angles

To appropriately use inverse kinematic (IK) modelling for the assessment of human motion, a musculoskeletal model must be prepared 1) to match participant segment lengths (scaling) and 2) to align the model׳s virtual markers positions with known, experimentally derived kinematic marker positions (marker registration). The purpose of this study was to investigate whether prescribing joint co-ordinates during the marker registration process (within the modelling framework OpenSim) will improve IK derived elbow kinematics during an overhead sporting task. To test this, the upper limb kinematics of eight cricket bowlers were recorded during two testing sessions, with a different tester each session. The bowling trials were IK modelled twice: once with an upper limb musculoskeletal model prepared with prescribed participant specific co-ordinates during marker registration – MR_PC – and once with the same model prepared without prescribed co-ordinates – MR; and by an established direct kinematic (DK) upper limb model. Whilst both skeletal model preparations had strong inter-tester repeatability (MR: Statistical Parametric Mapping (SPM1D)=0% different; MR_PC: SPM1D=0% different), when compared with DK model elbow FE waveform estimates, IK estimates using the MR_PC model (RMSD=5.2±2.0°, SPM1D=68% different) were in closer agreement than the estimates from the MR model (RMSD=44.5±18.5°, SPM1D=100% different). Results show that prescribing participant specific joint co-ordinates during the marker registration phase of model preparation increases the accuracy and repeatability of IK solutions when modelling overhead sporting tasks in OpenSim.     

An image-based kinematic model of the tibiotalar and subtalar joints and its application to gait analysis in children with Juvenile Idiopathic Arthritis

In vivo estimates of tibiotalar and the subtalar joint kinematics can unveil unique information about gait biomechanics, especially in the presence of musculoskeletal disorders affecting the foot and ankle complex. Previous literature investigated the ankle kinematics on ex vivo data sets, but little has been reported for natural walking, and even less for pathological and juvenile populations. This paper proposes an MRI-based morphological fitting methodology for the personalised definition of the tibiotalar and the subtalar joint axes during gait, and investigated its application to characterise the ankle kinematics in twenty patients affected by Juvenile Idiopathic Arthritis (JIA). The estimated joint axes were in line with in vivo and ex vivo literature data and joint kinematics variation subsequent to inter-operator variability was in the order of 1°. The model allowed to investigate, for the first time in patients with JIA, the functional response to joint impairment. The joint kinematics highlighted changes over time that were consistent with changes in the patient’s clinical pattern and notably varied from patient to patient. The heterogeneous and patient-specific nature of the effects of JIA was confirmed by the absence of a correlation between a semi-quantitative MRI-based impairment score and a variety of investigated joint kinematics indexes. In conclusion, this study showed the feasibility of using MRI and morphological fitting to identify the tibiotalar and subtalar joint axes in a non-invasive patient-specific manner. The proposed methodology represents an innovative and reliable approach to the analysis of the ankle joint kinematics in pathological juvenile populations.     

Design of intramedullary prostheses to prevent bone loss: predictions based on damage-stimulated remodelling

The adaptation of bone around intramedullary fixated prostheses, such as the femoral component of the hip joint or the radial component of the elbow joint, is well documented in follow-up studies. Bone adaptation takes the form of proximal bone atrophy accompanied, in some cases, by distal bone hypertrophy. A mechanistic model has been formulated to predict bone adaptation based on the concept that the continuous processes of damage and repair regulate bone adaptation. We apply the model to investigate the significance of two features of intramedullary prosthesis design on bone adaptation: prosthesis Young's modulus and the presence of a prosthesis collar. Results, as well as indicating some characteristics of accumulative-damage stimulated bone adaptation, predict that a low Young's modulus stem will very much reduce the extent of bone loss whereas the presence of a collar will have no significant effect. The results predict that a collarless low stiffness prosthesis is one possible approach for improving the secondary stability of intramedullary-fixated orthopaedic implants.     

Feasibility of using combined EMG and kinematic signals for prosthesis control: A simulation study using a virtual reality environment

Transhumeral amputation has a significant effect on a person’s independence and quality of life. Myoelectric prostheses have the potential to restore upper limb function, however their use is currently limited due to lack of intuitive and natural control of multiple degrees of freedom. The goal of this study was to evaluate a novel transhumeral prosthesis controller that uses a combination of kinematic and electromyographic (EMG) signals recorded from the person’s proximal humerus. Specifically, we trained a time-delayed artificial neural network to predict elbow flexion/extension and forearm pronation/supination from six proximal EMG signals, and humeral angular velocity and linear acceleration. We evaluated this scheme with ten able-bodied subjects offline, as well as in a target-reaching task presented in an immersive virtual reality environment. The offline training had a target of 4° for flexion/extension and 8° for pronation/supination, which it easily exceeded (2.7° and 5.5° respectively). During online testing, all subjects completed the target-reaching task with path efficiency of 78% and minimal overshoot (1.5%). Thus, combining kinematic and muscle activity signals from the proximal humerus can provide adequate prosthesis control, and testing in a virtual reality environment can provide meaningful data on controller performance.     

Biomechanical causes of trapeziometacarpal arthroplasty failure

Trapeziometacarpal joint prosthesis revision has been widely reported, mainly due to loosening of the trapezium cup. Our hypothesis is that current prostheses do not sufficiently respect the kinematics of this joint. CT scan acquisitions enabled us to determine the position of the first metacarpal relative to the trapezium in three different characteristic postures, in subjects in different stages of arthrosis. A CAD model of a current prosthesis was inserted into the numerical 3D model of the joint under the different postures. In the numerical model, we observe penetration of the cup by the head of the prosthesis. This virtual penetration could, in vivo, amount to overstressing the prosthetic elements, which would lead to loosening of the cup or of the metacarpal stem and luxation of the prosthesis.     

Patient-Specific Design and Biomechanical Evaluation of a Novel Bipolar Femoral Hemi-Knee Prosthesis

While total knee replacement is successful, hemiarthroplasty is necessary for some young, obese and active patients who are especially not suitable for unicompartmental or total knee prostheses. Hemiarthroplasty also provides an opportunity for children with bone tumors. The design ofhemiarthroplasty should be patient-specific to reduce contact stress and friction as well as instability, compared to conventional hemi-knee prosthesis. A novel bipolar hemi-knee prosthesis with two flexion stages was developed according to a healthy male＇s knee morphological profile. The motion mode of the bipolar hemi-knee prosthesis was observed through roentgenoscopy in vitro experiment. The biomechanical properties in one gait cycle were evaluated though finite element simulation. The bipolar hemi-knee prosthesis was found to produce knee flexion at two stages through X-ray images. The first stage is the motion from upright posture to a specified 60~ flexion, followed by the second stage of motion subsequently to deep flexion. The finite element simulation results also show that the designed hemi-knee prosthesis has the ability to reduce stresses on the joint contact surfaces. Therefore, it is possible for the bipolar hemi-knee prosthesis to provide better biotribological performances because it can reduce stresses and potentially wear on the opposing contacting surface during a gait cycle, orovidin~ a t~romisin~ treatment strate~v in future Joint renair znd renlneement     

Repeatability,reproducibility, and agreement of three computational methods to approximate the functional flexion-extension axis of the tibiofemoral joint using 3D bone models of the femur

Abstract

Background: Closely approximating the functional flexion-extension (FE) axis of the tibiofemoral joint in 3D models of the femur is important when computing joint motions which are physiologic. The objectives were to 1) develop methods to approximate the functional FE axis based on fitting circles, a tapered cylinder, and spheres to the posterior condyles, 2) determine the repeatability and reproducibility of each method, and 3) determine limits of agreement between pairs of axes.

Methods: For each method, the respective axis was determined in forty 3D bone models of the distal femur. Varus-valgus angles and internal-external axial angles were computed relative to standard planes.

Results: Repeatability and reproducibility were comparable for the tapered cylinder-based and sphere-based methods and better than that for the circle-based method. Limits of agreement were tightest when comparing the sphere-based and tapered cylinder-based axes. However, limits of agreement for the internal-external axial angle were wide at +3.6° to ?3.9° whereas limits of agreement were tighter at +1.4° to ?0.7° for the varus-valgus angle.

Conclusion: The tapered cylinder-based and sphere-based methods offer advantages of better repeatability and reproducibility over the circle-based method. However, the sphere-based and tapered cylinder-based axes are not interchangeable owing to wide limits of agreement for the internal-external axial angle. The tapered cylinder-based axis is preferred intuitively over the sphere-based axis because the spheres require fitting in both the sagittal and coronal planes whereas the tapered cylinder requires fitting in the sagittal plane only which is the plane of motion in flexion-extension.     

Computational analysis of glenohumeral joint growth and morphology following a brachial plexus birth injury

Children affected with brachial plexus birth injury (BPBI) undergo muscle paralysis. About 33% of affected children experience permanent osseous deformities of the glenohumeral joint. Recent evidence suggests that some cases experience restricted muscle longitudinal growth in addition to paralysis and reduced range of motion at the shoulder and elbow. It is unknown whether altered loading due to paralysis, muscle growth restriction and contracture, or static loading due to disuse is the primary driver of joint deformity after BPBI. This study uses a computational framework integrating finite element analysis and musculoskeletal modeling to examine the mechanical factors contributing to changes in bone growth and morphometry following BPBI. Simulations of 8 weeks of glenohumeral growth in a rat model of BPBI predicted that static loading of the joint is primarily responsible for joint deformation consistent with experimental measures of bone morphology, whereas dynamic loads resulted in normal bone growth. Under dynamic loading, glenoid version angle (GVA), glenoid inclination angle (GIA), and glenoid radius of curvature (GRC) (−1.3°, 38.2°, 2.5 mm respectively) were similar to the baseline values (−1.8°, −38°, 2.1 mm respectively). In the static case with unrestricted muscle growth, these measures increased in magnitude (5.2°, −48°, 3.5 mm respectively). More severe joint deformations were observed in GIA and GRC when muscle growth was restricted (GVA: 3.6°, GIA: −55°, GRC: 4.0 mm). Predicted morphology was consistent with literature reports of in vivo glenoid morphology following postganglionic BPBI. This growth model provides a framework for understanding the most influential mechanical factors driving glenohumeral deformity following BPBI.