A better understanding of how the shape and density of the shoulder vary among members of a population can help design more effective population-based orthopedic implants. The main objective of this study was to develop statistical shape models (SSMs) and statistical density models (SDMs) of the shoulder to describe the main modes of variability in the shape and density distributions of shoulder bones within a population in terms of principal components (PCs). These PC scores were analyzed, and significant correlations were observed between the shape and density distributions of the shoulder and demographics of the population, such as sex and age. Our results demonstrated that when the overall body sizes of male and female donors were matched, males still had, on average, larger scapulae and thicker humeral cortical bones. Moreover, we concluded that age has a weak but significant inverse effect on the density within the entire shoulder. Weak and moderate, but significant, correlations were also found between many modes of shape and density variations in the shoulder. Our results suggested that donors with bigger humeri have bigger scapulae and higher bone density of humeri corresponds with higher bone density in the scapulae. Finally, asymmetry, to some extent, was noted in the shape and density distributions of the contralateral bones of the shoulder. These results can be used to help guide the designs of population-based prosthesis components and pre-operative surgical planning.

Soft tissues observed in clinical medical images are often prestrained in tension by internal pressure or tissue hydration. For a native disc, nucleus swelling occurs in equilibrium with osmotic pressure induced by the high concentration of proteoglycan in the nucleus. The objective of this computational study was to investigate the effects of nucleus swelling on disc geometry, fiber orientation, and mechanical behavior by comparing those of prestrained and zero-pressure (unswelled) discs. Thermoelastic analysis techniques were repurposed in order to determine the zero-pressure disc geometry which, when pressurized, matches the prestrained disc geometry observed in clinical images. The zero-pressure geometry was then used in simulations to approximately represent a degenerated disc, which loses the ability of nucleus swelling but has not undergone distinct soft tissue remodeling/disruption. Our simulation results demonstrated that the loss of nucleus swelling caused a slight change in the disc geometry and fiber orientation, but a distinct deterioration in the resistance to intervertebral rotations including sagittal bending, lateral bending, and axial torsion. Different from rotational loading, in compression (with a displacement of 0.45 mm applied), a much larger stiffness (3.02 KN/mm) and a greater intradiscal pressure (IDP) (0.61 MPa) were measured in the zero-pressure disc, compared to the prestrained disc (1.41 KN/mm and 0.52 MPa). This computational study could be useful to understand mechanisms of disc degeneration, and guide the future design of disc tissue engineering material and biomimic disc implants.

Total disk arthroplasty (TDA) using an artificial disk (AD) is an attractive surgical technique for the treatment of spinal disorders, since it can maintain or restore spinal motion (unlike interbody fusion). However, adverse surgical outcomes of contemporary lumbar TDAs have been reported. We previously proposed a new mobile-bearing AD design concept featuring a biconcave ultrahigh-molecular-weight polyethylene (UHMWPE) mobile core. The objective of this study was to develop an artificial neural network (NN) based multiobjective optimization framework to refine the biconcave-core AD design considering multiple TDA performance metrics, simultaneously. We hypothesized that there is a tradeoff relationship between the performance metrics in terms of range of motion (ROM), facet joint force (FJF), and polyethylene contact pressure (PCP). By searching the resulting three-dimensional (3D) Pareto frontier after multiobjective optimization, it was found that there was a “best-tradeoff” AD design, which could balance all the three metrics, without excessively sacrificing each metric. However, for each single-objective optimum AD design, only one metric was optimal, and distinct sacrifices were observed in the other two metrics. For a commercially available biconvex-core AD design, the metrics were even worse than the poorest outcomes of the single-objective optimum AD designs. Therefore, multiobjective design optimization could be useful for achieving native lumbar segment biomechanics and minimal PCPs, as well as for improving the existing lumbar motion-preserving surgical treatments.

Much of our understanding of the role of elbow ligaments to overall joint biomechanics has been developed through in vitro cadaver studies using joint motion simulators. The principle of superposition can be used to indirectly compute the force contributions of ligaments during prescribed motions. Previous studies have analyzed the contribution of different soft tissue structures to the stability of human elbow joints, but have limitations in evaluating the loads sustained by those tissues. This paper introduces a unique, hybrid experimental-computational technique for measuring and simulating the biomechanical contributions of ligaments to elbow joint kinematics and stability. in vitro testing of cadaveric joints is enhanced by the incorporation of fully parametric virtual ligaments, which are used in place of the native joint stabilizers to characterize the contribution of elbow ligaments during simple flexion–extension (FE) motions using the principle of superposition. Our results support previously reported findings that the anterior medial collateral ligament (AMCL) and the radial collateral ligament (RCL) are the primary soft tissue stabilizers for the elbow joint. Tuned virtual ligaments employed in this study were able to restore the kinematics and laxity of elbows to within 2 deg of native joint behavior. The hybrid framework presented in this study demonstrates promising capabilities in measuring the biomechanical contribution of ligamentous structures to joint stability.

In the present study a direct comparison was made between in vitro total hip wear testing and a computational analysis considering the effects of time and a nonlinear stress-strain relationship for ultrahigh molecular weight polyethylene (UHMWPE) at 37 ° C . The computational simulation was made correct through calibration to experimental volumetric wear results, and the predicted damage layout on the acetabular liner surface was compared with results estimated from laser scanning of the actual worn specimens. The wear rates for the testing specimens were found to be 17.14 ± 1.23 mg / 10 6 cycles and 19.39 ± 0.79 mg / 10 6 cycles , and the cumulative volumetric wear values after 3 × 10 6 cycles were 63.70 mm 3 and 64.02 mm 3 for specimens 1 and 2, respectively. The value of the calibrated wear coefficient was found to be 5.32 ( 10 − 10 ) mm 3 / N mm for both specimens. The major difference between the computational and experimental wear results was the existence of two damage vectors in the experimental case. The actual location of damage was virtually the same in both cases, and the maximum damage depth of the computational model agreed well with the experiment. The existence of multiple wear vectors may indicate the need for computational approaches to account for multidirectional sliding or strain hardening of UHMWPE. Despite the limitation in terms of describing the overall damage layout, the present computational model shows that simulation can mimic some of the behavior of in vitro wear.