Finite element modelling of the glenohumeral capsule can help assess the tested region during a clinical exam

The objective of this research was to examine the efficacy of evaluating the region of the glenohumeral capsule being tested by clinical exams for shoulder instability using finite element (FE) models of the glenohumeral joint. Specifically, the regions of high capsule strain produced by glenohumeral joint positions commonly used during a clinical exam were identified. Kinematics that simulated a simple translation test with an anterior load at three external rotation angles were applied to a validated, subject-specific FE model of the glenohumeral joint at 60° of abduction. Maximum principal strains on the glenoid side of the inferior glenohumeral ligament (IGHL) were significantly higher than the maximum principal strains on the humeral side, for all three regions of the IGHL at 30° and 60° of external rotation. These regions of localised strain indicate that these joint positions might be used to test the glenoid side of the IGHL during this clinical exam, but are not useful for assessing the humeral side of the IGHL. The use of FE models will facilitate the search for additional joint positions that isolate high strains to other IGHL regions, including the humeral side of the IGHL.     

Knee joint passive stiffness and moment in sagittal and frontal planes markedly increase with compression

Knee joints are subject to large compression forces in daily activities. Due to artefact moments and instability under large compression loads, biomechanical studies impose additional constraints to circumvent the compression position–dependency in response. To quantify the effect of compression on passive knee moment resistance and stiffness, two validated finite element models of the tibiofemoral (TF) joint, one refined with depth-dependent fibril-reinforced cartilage and the other less refined with homogeneous isotropic cartilage, are used. The unconstrained TF joint response in sagittal and frontal planes is investigated at different flexion angles (0°, 15°, 30° and 45°) up to 1800 N compression preloads. The compression is applied at a novel joint mechanical balance point (MBP) identified as a point at which the compression does not cause any coupled rotations in sagittal and frontal planes. The MBP of the unconstrained joint is located at the lateral plateau in small compressions and shifts medially towards the inter-compartmental area at larger compression forces. The compression force substantially increases the joint moment-bearing capacities and instantaneous angular rigidities in both frontal and sagittal planes. The varus–valgus laxities diminish with compression preloads despite concomitant substantial reductions in collateral ligament forces. While the angular rigidity would enhance the joint stability, the augmented passive moment resistance under compression preloads plays a role in supporting external moments and should as such be considered in the knee joint musculoskeletal models.     

Finding consistent strain distributions in the glenohumeral capsule between two subjects: implications for development of physical examinations

The anterior-inferior glenohumeral capsule is the primary passive stabilizer to the glenohumeral joint during anterior dislocation. Physical examinations following dislocation are crucial for proper diagnosis of capsule pathology; however, they are not standardized for joint position which may lead to misdiagnoses and poor outcomes. To suggest joint positions for physical examinations where the stability provided by the capsule may be consistent among patients, the objective of this study was to evaluate the distribution of maximum principal strain on the anterior-inferior capsule using two validated subject-specific finite element models of the glenohumeral joint at clinically relevant joint positions. The joint positions with 25 N anterior load applied at 60° of glenohumeral abduction and 10°, 20°, 30° and 40° of external rotation resulted in distributions of strain that were similar between shoulders (r2 ≥ 0.7). Furthermore, those positions with 20-40° of external rotation resulted in capsule strains on the glenoid side of the anterior band of the inferior glenohumeral ligament that were significantly greater than in all other capsule regions. These findings suggest that anterior stability provided by the anterior-inferior capsule may be consistent among subjects at joint positions with 60° of glenohumeral abduction and a mid-range (20-40°) of external rotation, and that the glenoid side has the greatest contribution to stability at these joint positions. Therefore, it may be possible to establish standard joint positions for physical examinations that clinicians can use to effectively diagnose pathology in the anterior-inferior capsule following dislocation and lead to improved outcomes.     

A combined passive and active musculoskeletal model study to estimate L4-L5 load sharing

A number of geometrically-detailed passive finite element (FE) models of the lumbar spine have been developed and validated under in vitro loading conditions. These models are devoid of muscles and thus cannot be directly used to simulate in vivo loading conditions acting on the lumbar joint structures or spinal implants. Gravity loads and muscle forces estimated by a trunk musculoskeletal (MS) model under twelve static activities were applied to a passive FE model of the L4-L5 segment to estimate load sharing among the joint structures (disc, ligaments, and facets) under simulated in vivo loading conditions. An equivalent follower (FL), that generates IDP equal to that generated by muscle forces, was computed in each task. Results indicated that under in vivo loading conditions, the passive FE model predicted intradiscal pressures (IDPs) that closely matched those measured under the simulated tasks (R² = 0.98 and root-mean-squared-error, RMSE = 0.18 MPa). The calculated equivalent FL compared well with the resultant force of all muscle forces and gravity loads acting on the L4-L5 segment (R² = 0.99 and RMSE = 58 N). Therefore, as an alternative approach to represent in vivo loading conditions in passive FE model studies, this FL can be estimated by available in-house or commercial MS models. In clinical applications and design of implants, commonly considered in vitro loading conditions on the passive FE models do not adequately represent the in vivo loading conditions under muscle exertions. Therefore, more realistic in vivo loading conditions should instead be used.     

Collision tumor with inflammatory breast carcinoma and malignant phyllodes tumor: a case report and literature review

Background

Pelvic reconstruction after hemipelvectomy can greatly improve the weight-bearing stability of the supporting skeleton and improve patients’ quality of life. Although an autograft can be used to reconstruct pelvic defects, the most suitable choice of autograft, i.e., the use of either femur or tibia, has not been determined. We aimed to analyze the mechanical stresses of a pelvic ring reconstructed using femur or tibia after hemipelvectomy using finite element (FE) analysis.

Methods

FE models of normal and reconstructed pelvis were established based on computed tomography images, and the stress distributions were analyzed under physiological loading from 0 to 500 N in both intact and restored pelvic models using femur or tibia.

Results

The vertical displacement of the intact pelvis was less than that of reconstructed pelvis, but there was no significant difference between the two reconstructed models. In FE analysis, the stress distribution of the intact pelvic model was bilaterally symmetric and the maximum stresses were located at the sacroiliac joint, arcuate line, ischiatic ramus, and ischial tuberosity. The maximum stress in each part of the reconstructed pelvis greatly exceeded that of the intact model. The maximum von Mises stress of the femur was 13.9 MPa, and that of the tibia was 6.41 MPa. However, the stress distribution was different in the two types of reconstructed pelvises. The tibial reconstruction model induced concentrated stress on the tibia shaft making it more vulnerable to fracture. The maximum stress on the femur was concentrated on the connections between the femur and the screws.

Conclusions

From a biomechanical point of view, the reconstruction of hemipelvic defects with femur is a better choice.     

Biomechanical comparison of conventional and anatomical calcaneal plates for the treatment of intraarticular calcaneal fractures – a finite element study

Initial stability is essential for open reduction internal fixation of intraarticular calcaneal fractures. Geometrical feature of a calcaneal plate is influential to its endurance under physiological load. It is unclear if conventional and pre-contoured anatomical calcaneal plates may exhibit differently in biomechanical perspective. A Sanders’ Type II-B intraarticular calcaneal fracture model was reconstructed to evaluate the effectiveness of calcaneal plates using finite element methods. Incremental vertical joint loads up to 450 N were exerted on the subtalar joint to evaluate the stability and safety of the calcaneal plates and bony structure. Results revealed that the anatomical calcaneal plate model had greater average structural stiffness (585.7 N/mm) and lower von Mises stress on the plate (774.5 MPa) compared to those observed in the conventional calcaneal plate model (stiffness: 430.9 N/mm; stress on plate: 867.1 MPa). Although both maximal compressive and maximal tensile stress and strain were lower in the anatomical calcaneal plate group, greater loads on fixation screws were found (average 172.7 MPa compared to 82.18 MPa in the conventional calcaneal plate). It was noted that high magnitude stress concentrations would occur where the bone plate bridges the fracture line on the lateral side of the calcaneus bone. Sufficient fixation strength at the posterolateral calcaneus bone is important for maintaining subtalar joint load after reduction and fixation of a Sanders’ Type II-B calcaneal fracture. In addition, geometrical design of a calcaneal plate should worth considering for the mechanical safety in practical usage.     

Interactive graph-cut segmentation for fast creation of finite element models from clinical ct data for hip fracture prediction

In this study, we propose interactive graph cut image segmentation for fast creation of femur finite element (FE) models from clinical computed tomography scans for hip fracture prediction. Using a sample of N = 48 bone scans representing normal, osteopenic and osteoporotic subjects, the proximal femur was segmented using manual (gold standard) and graph cut segmentation. Segmentations were subsequently used to generate FE models to calculate overall stiffness and peak force in a sideways fall simulations. Results show that, comparable FE results can be obtained with the graph cut method, with a reduction from 20 to 2–5 min interaction time. Average differences between segmentation methods of 0.22 mm were not significantly correlated with differences in FE derived stiffness (R² = 0.08, p = 0.05) and weakly correlated to differences in FE derived peak force (R² = 0.16, p = 0.01). We further found that changes in automatically assigned boundary conditions as a consequence of small segmentation differences were significantly correlated with FE derived results. The proposed interactive graph cut segmentation software MITK-GEM is freely available online at https://simtk.org/home/mitk-gem.     

Modulation of shoulder muscle and joint function using a powered upper-limb exoskeleton

Robotic-assistive exoskeletons can enable frequent repetitive movements without the presence of a full-time therapist; however, human-machine interaction and the capacity of powered exoskeletons to attenuate shoulder muscle and joint loading is poorly understood. This study aimed to quantify shoulder muscle and joint force during assisted activities of daily living using a powered robotic upper limb exoskeleton (ArmeoPower, Hocoma). Six healthy male subjects performed abduction, flexion, horizontal flexion, reaching and nose touching activities. These tasks were repeated under two conditions: (i) the exoskeleton compensating only for its own weight, and (ii) the exoskeleton providing full upper limb gravity compensation (i.e., weightlessness). Muscle EMG, joint kinematics and joint torques were simultaneously recorded, and shoulder muscle and joint forces calculated using personalized musculoskeletal models of each subject’s upper limb. The exoskeleton reduced peak joint torques, muscle forces and joint loading by up to 74.8% (0.113 Nm/kg), 88.8% (5.8%BW) and 68.4% (75.6%BW), respectively, with the degree of load attenuation strongly task dependent. The peak compressive, anterior and superior glenohumeral joint force during assisted nose touching was 36.4% (24.6%BW), 72.4% (13.1%BW) and 85.0% (17.2%BW) lower than that during unassisted nose touching, respectively. The present study showed that upper limb weight compensation using an assistive exoskeleton may increase glenohumeral joint stability, since deltoid muscle force, which is the primary contributor to superior glenohumeral joint shear, is attenuated; however, prominent exoskeleton interaction moments are required to position and control the upper limb in space, even under full gravity compensation conditions. The modeling framework and results may be useful in planning targeted upper limb robotic rehabilitation tasks.     

The use of preserved tissue in finite element modelling of fresh ovine vertebral behaviour

The aim of this study was to investigate whether the predicted finite element (FE) stiffness of vertebral bone is altered when using images of preserved rather than fresh tissue to generate specimen-specific FE models. Fresh ovine vertebrae were used to represent embalmed (n = 3) and macerated dry-bone (n = 3) specimens and treated accordingly. Specimens were scanned pre- and post-treatment using micro-computed tomography. A constant threshold level derived from these images was used to calculate the respective bone volume fraction (BV/TV) from which the conversion factor validated for fresh tissue was used to determine material properties that were assigned to corresponding FE models. Results showed a definite change in the BV/TV between the fresh and the preserved bone. However, the changes in the predicted FE stiffness were not generally greater than the variations expected from assignment of loading and boundary conditions. In conclusion, images of preserved tissue can be used to generate FE models that are representative of fresh tissue with a tolerable level of error.     

Comparison of an EMG-based and a stress-based method to predict shoulder muscle forces

The estimation of muscle forces in musculoskeletal shoulder models is still controversial. Two different methods are widely used to solve the indeterminacy of the system: electromyography (EMG)-based methods and stress-based methods. The goal of this work was to evaluate the influence of these two methods on the prediction of muscle forces, glenohumeral load and joint stability after total shoulder arthroplasty. An EMG-based and a stress-based method were implemented into the same musculoskeletal shoulder model. The model replicated the glenohumeral joint after total shoulder arthroplasty. It contained the scapula, the humerus, the joint prosthesis, the rotator cuff muscles supraspinatus, subscapularis and infraspinatus and the middle, anterior and posterior deltoid muscles. A movement of abduction was simulated in the plane of the scapula. The EMG-based method replicated muscular activity of experimentally measured EMG. The stress-based method minimised a cost function based on muscle stresses. We compared muscle forces, joint reaction force, articular contact pressure and translation of the humeral head. The stress-based method predicted a lower force of the rotator cuff muscles. This was partly counter-balanced by a higher force of the middle part of the deltoid muscle. As a consequence, the stress-based method predicted a lower joint load (16% reduced) and a higher superior–inferior translation of the humeral head (increased by 1.2 mm). The EMG-based method has the advantage of replicating the observed cocontraction of stabilising muscles of the rotator cuff. This method is, however, limited to available EMG measurements. The stress-based method has thus an advantage of flexibility, but may overestimate glenohumeral subluxation.     

Stress analysis of irradiated human tooth enamel using finite element methods

The objectives of this project were to use finite element methods to determine how changes in the elastic modulus due to oral cancer therapeutic radiation alter the distribution of mechanical stresses in teeth and to determine if observed failures in irradiated teeth correlate with changes in mechanical stresses. A thin slice section finite element (FE) model was constructed from micro CT sections of a molar tooth using MIMICS and 3-Matic software. This model divides the tooth into three enamel regions, the dentin-enamel junction (DEJ) and dentin. The enamel elastic modulus was determined in each region using nano indentation for three experimental groups namely – control (non-radiated), in vitro irradiated (simulated radiotherapy following tooth extraction) and in vivo irradiated (extracted subsequent to oral cancer patient radiotherapy) teeth. Physiological loads were applied to the tooth models at the buccal and lingual cusp regions for all three groups (control, in vitro and in vivo). The principal tensile stress and the maximum shear stress were used to compare the results from different groups since it has been observed in previous studies that delamination of enamel from the underlying dentin was one of the major reasons for the failure of teeth following therapeutic radiation. From the FE data, we observed an increase in the principal tensile stress within the inner enamel region of in vivo irradiated teeth (9.97 ± 1.32 MPa) as compared to control/non-irradiated teeth (8.44 ± 1.57 MPa). Our model predicts that failure occurs at the inner enamel/DEJ interface due to extremely high tensile and maximum shear stresses in in vivo irradiated teeth which could be a cause of enamel delamination due to radiotherapy.     

Can a partial volume edge effect reduction algorithm improve the repeatability of subject-specific finite element models of femurs obtained from CT data?

The reliability of patient-specific finite element (FE) modelling is dependent on the ability to provide repeatable analyses. Differences of inter-operator generated grids can produce variability in strain and stress readings at a desired location, which are magnified at the surface of the model as a result of the partial volume edge effects (PVEEs). In this study, a new approach is introduced based on an in-house developed algorithm which adjusts the location of the model's surface nodes to a consistent predefined threshold Hounsfield unit value. Three cadaveric human femora specimens were CT scanned, and surface models were created after a semi-automatic segmentation by three different experienced operators. A FE analysis was conducted for each model, with and without applying the surface-adjustment algorithm (a total of 18 models), implementing identical boundary conditions. Maximum principal strain and stress and spatial coordinates were probed at six equivalent surface nodes from the six generated models for each of the three specimens at locations commonly utilised for experimental strain guage measurement validation. A Wilcoxon signed-ranks test was conducted to determine inter-operator variability and the impact of the PVEE-adjustment algorithm. The average inter-operator difference in stress values was significantly reduced after applying the adjustment algorithm (before: 3.32 ± 4.35 MPa, after: 1.47 ± 1.77 MPa, p = 0.025). Strain values were found to be less sensitive to inter-operative variability (p = 0.286). In summary, the new approach as presented in this study may provide a means to improve the repeatability of subject-specific FE models of bone obtained from CT data.     

Mechanical properties of the avian acrocoracohumeral ligament and its role in shoulder stabilization in flight

Control of movement in the avian shoulder joint is fundamental to understanding the avian wingstroke. The acrocoracohumeral ligament (AHL) is thought to play a key role in stabilizing the glenoid and balancing the pectoralis in gliding flight. If the AHL has to be taut to balance the pectoralis, then it must constrain glenohumeral motion during flapping flight as well. However, birds vary wing kinematics depending on flight speed and behavior. How can a passive ligament accommodate such varying joint movements? Herein, mechanical testing and 3-D modeling are used to link the mechanical properties and morphology of the AHL to its functional role during flapping flight. The bone-ligament-bone complex of the pigeon (Columba livia) fails at a tensile loading of 141 ± 18 N (± s .D., n = 10) or 39 times body weight, which corresponds to a failure stress of 51 MPa, well above expected loads during flight. Simulated AHL length changes, comparisons to glenohumeral kinematics from the literature, and manipulations of partially dissected pigeon specimens all support the hypothesis that the AHL remains taut through downstroke and most of upstroke while becoming slack during the downstroke/upstroke transition. The digital AHL model provides a mechanism for explaining how the AHL can stabilize the shoulder joint under a broad array of humeral paths by constraining the coordination of glenohumeral degrees of freedom.     

A musculoskeletal shoulder simulation of moment arms and joint reaction forces after medialization of the supraspinatus footprint in rotator cuff repair

A non-anatomical reinsertion of the supraspinatus medially to the original footprint to avoid over-tensioning of the tendon in large and retracted tears is one surgical option in rotator cuff (RC) repair. The purpose of the study was to determine the biomechanical effects on the glenohumeral joint with regard to this surgical technique. A modified musculoskeletal computational shoulder model was used to evaluate the change in moment arms and muscle forces of the RC and the co-contracting muscles and the alteration of the joint reaction forces (compressive and shear forces) after reinsertion of the supraspinatus 5?mm, 10?mm, 15?mm and 20?mm medially to the original footprint. A medialization of the supraspinatus reduces its moment arm in glenohumeral abduction. In case of a medialization of the attachment of 15?mm and 20?mm, the supraspinatus restricts glenohumeral abduction at 54° and 68°. In glenohumeral forward flexion and in lower degrees of internal rotation the moment arm of the supraspinatus increases for a medialized tendon attachment and decreases in external rotation in relation to the anatomical condition. A medialization of the supraspinatus insertion point yields in an increase in muscle force for abduction, internal and external rotation. In the present model a medially non-anatomic reinsertion reduces significantly the compressive glenohumeral joint reaction and the glenohumeral stability. Moreover, the results show that a medialization of the supraspinatus leads to a reduction of the supraspinatus moment arm especially in abduction. This leads to an increase of a compensatory supraspinatus load for stabilization the humerus in space, which may potentially cause a postoperative overload of the tendon-bone-complex.     

Steeper posterior tibial slope markedly increases ACL force in both active gait and passive knee joint under compression

The role of the posterior tibial slope (PTS) in anterior cruciate ligament (ACL) risk of injury has been supported by many imaging studies but refuted by some in vitro works. The current investigation was carried out to compute the effect of ±5^o change in PTS on knee joint biomechanics in general and ACL force/strain in particular. Two validated finite element (FE) models of the knee joint were employed; one active lower extremity musculoskeletal model including a complex FE model of the knee joint driven by in vivo kinematics/kinetics collected in gait of asymptomatic subjects, and the other its isolated unconstrained passive tibiofemoral (TF) joint considered under 1400 N compression at four different knee flexion angles (0°–45°). In the TF model, the compression force was applied at the joint mechanical balance point causing no rotations in sagittal and frontal planes.     

In vivo measurement of shoulder joint loads during activities of daily living

Until recently the contact loads acting in the glenohumeral joint have been calculated using musculoskeletal models or measured in vitro. Now, contact forces and moments are measured in vivo using telemeterized shoulder implants. Mean total contact forces from four patients during eight activities of daily living are reported here.Lifting a coffee pot (1.5 kg) with straight arm caused an average force of 105.0%BW (%body weight) (range: 90–124.6%BW), while setting down the coffee pot in the same position led to higher forces of 122.9%BW on the average (105.3–153.4%BW). The highest joint contact forces were measured when the straight arm was abducted or elevated by 90° or more, with a weight in the hand. Lifting up 2 kg from a board up to head height caused a contact force of 98.3%BW (93–103.6%BW); again, setting it down on the board led to higher forces of 131.5%BW (118.8–144.1%BW). In contrast to previously calculated high loads, the contact force during passive holding of a 10 kg weight laterally was only 12.3%BW (9.2–17.9%BW), but when lifting it up to belt height it increased to 91.5%BW (87–95%BW).The moments transferred inside the joint at our patients varied much more than did the forces both inter and intra-individually.Our data suggest that patients with shoulder problems or during the first post-operative weeks after shoulder fractures or joint replacements should avoid certain activities encountered during daily living e.g. lifting or holding a weight with an outstretched arm. Some energy-related optimization criteria used in the literature for analytical musculoskeletal shoulder models must now be reconsidered.     

Finite element modeling of arachnid slit sensilla: II. Actual lyriform organs and the face deformations of the individual slits

Arachnid slit sensilla respond to minute strains in the exoskeleton. After having applied Finite Element (FE) analysis to simplified arrays of five straight slits (Hößl et al. J Comp Physiol A 193:445–459, 2007) we now present a computational study of the effects of more subtle natural variations in geometry, number and arrangement of slits on the slit face deformations. Our simulations show that even minor variations in these parameters can substantially influence a slit’s directional response. Using white-light interferometric measurements of the surface deformations of a lyriform organ, it is shown that planar FE models are capable of predicting the principal characteristics of the mechanical responses. The magnitudes of the measured and calculated slit face deformations are in good agreement. At threshold, they measure between 1.7 and 43 nm. In a lyriform organ and a closely positioned loose group of slits, the detectable range of loads increases to approximately 3.5 times the range of the lyriform organ alone. Stress concentration factors (up to ca. 29) found in the vicinity of the slits were evaluated from the models. They are mitigated due to local thickening of the exocuticle and the arrangement of the chitinous microfibers that prevents the formation of cracks under physiological loading conditions.     

The effect of muscle loading on the kinematics of in vitro glenohumeral abduction

This in vitro study evaluated the effects of four different muscle-loading ratios on active glenohumeral joint abduction. Eight cadaveric shoulders were tested using a shoulder simulator designed to reproduce unconstrained abduction of the humerus via computer-controlled pneumatic actuation. Forces were applied to cables that were sutured to tendons or fixed to bone, to simulate loading of the supraspinatus, subscapularis, infraspinatus/teres minor, and anterior, middle, and posterior deltoid muscles. Four sets of muscle-loading ratios were employed, based on: (1) equal loads, (2) average physiological cross-sectional areas (pCSAs), (3) constant values of the product of electromyographic (EMG) data and pCSAs, and (4) variable ratios of the EMG and pCSA data which changed as a function of abduction angle. The investigator generated passive motions with no muscle loads simulated. Repeatability was quantified by five successive trials of the passive and simulated active motions. There was improved repeatability in the simulated active motions versus passive motions, significant for abduction angles less than 40 degrees (p=0.02). No difference was found in the repeatability of the four different muscle-loading ratios for simulated active motions (p0.067 for all angles). The improved repeatability of active over passive motion suggests simulated active motion should be employed for in vitro simulations of shoulder motion.     

Computational modelling of mobile bearing TKA anterior–posterior dislocation

Anterior–posterior stability in an unconstrained mobile-bearing total knee arthroplasty (TKA) and one with rotational constraints is compared in a computational model based on an ASTM test. Both TKA designs dislocate at loads greater than reported maximum in vivo forces. The posterior drawer forces (mean: 3027 N vs. 1817 N) needed to induce subluxation increase with a greater anterior jump distance (12 mm vs. 7 mm; refers to the vertical height of the anterior or posterior border of the tibial insert's articulating surface). The posterior jump distance for both tested TKA differed by 1.5 mm and had minimal effect on the magnitude of the anterior drawer forces at dislocation in mid-flexion (unconstrained vs. constrained: 445 N vs. 412 N). The unconstrained insert dislocated by means of spin-out whereas in the constrained TKA the femur dislocated from the bearing during posterior drawer and the bearing from the baseplate during anterior drawer. MCL function is an important consideration during ligament balancing since a ± 10% variation in MCL tension affects dislocation forces by ± 20%. The simulation platform provided the means to investigate TKA designs in terms of anterior–posterior stability as a function of knee flexion, collateral ligament function and mechanical morphology.     

Subchondral bone microarchitecture and failure mechanism under compression: A finite element study

Subchondral bone (SCB) microdamage is commonly observed in traumatic joint injuries and has been strongly associated with post-traumatic osteoarthritis (PTOA). Knowledge of the three-dimensional stress and strain distribution within the SCB tissue helps to understand the mechanism of SCB failure, and may lead to an improved understanding of mechanisms of PTOA initiation, prevention and treatment. In this study, we used high-resolution micro-computed tomography (µCT)-based finite element (FE) modelling of cartilage-bone to evaluate the failure mechanism and the locations of SCB tissue at high-risk of initial failure under compression. The µCT images of five cartilage-bone specimens with an average SCB thickness of 1.23 ± 0.20 mm were used to develop five µCT-based FE models. The FE models were analysed under axial compressions of approximately 30 MPa applied to the cartilage surface while the bone edges were constrained. Strain and stress-based failure criteria were then applied to evaluate the failure mechanism of the SCB tissue under excessive compression through articular cartilage. µCT-based FE models predicted two locations in the SCB at high-risk of initial failure: (1) the interface of the calcified-uncalcified cartilage due to excessive tension, and (2) the trabecular bone beneath the subchondral plate due to excessive compression. µCT-based FE models of cartilage-bone enabled us to quantify the distribution of the applied compression which was transferred through the articular cartilage to its underlying SCB, and to investigate the mechanism and the mode of SCB tissue failure. Ultimately, the results will help to understand the mechanism of injury formation in relation to PTOA.