Spinal loads during dynamic full flexion and return to standing posture in different age and sex groups: A musculoskeletal model study.

During forward flexion, spine motion varies due to age and sex differences. Previous studies showed that lumbar/pelvis range of flexion (RoF) and lumbo-pelvic ratio (L/P) are age/sex dependent. How variation of these parameters affects lumbar loading in a normal population requires further assessment. We aimed to estimate lumbar loads during dynamic flexion-return cycle and the differences in peak loads (compression) and corresponding trunk inclinations due to variation in lumbar/pelvis RoF and L/P. Based on in vivo L/P (0.11-3.44), temporal phases of flexion (early, middle, and later), the lumbar (45-55°) and hip (60-79°) RoF; full flexion-return cycles of six seconds were reconstructed for three age groups (20-35, 36-50 and 50+ yrs.) in both sexes. Six inverse dynamic analyses were performed with a 50th percentile model, and differences in peak loads and corresponding trunk inclinations were calculated. Peak loads at L4-L5 were 179 N higher in younger males versus females, but 228 N and 210 N lower in middle-aged and older males, respectively, compared to females. Females exhibited higher trunk inclinations (6°-20°) than males across all age groups. Age related differences in L4-L5 peak loads and corresponding trunk inclinations were found up to 415 N and 19° in males and 152 N and 13° in females. With aging, peak loads were reduced in males but were found non-monotonic in females, whereas trunk inclinations at peak loads were reduced in both sexes from young to middle/old age groups. In conclusion, lumbar loading and corresponding trunk inclinations varied notably due to age/sex differences. Such data may help distinguishing normal or pathological condition of the lumbar spine.

Subject-specific trunk segmental masses prediction for musculoskeletal models using artificial neural networks.

Accurate determination of body segment parameters is crucial for studying human movement and joint forces using musculoskeletal (MSK) models. However, existing methods for predicting segment mass have limited generalizability and sensitivity to body shapes. With recent advancements in machine learning, this study proposed a novel artificial neural network-based method for computing subject-specific trunk segment mass and center of mass (CoM) using only anthropometric measurements. We first developed, trained, and validated two artificial neural networks that used anthropometric measurements as input to predict body shape (ANN1) and tissue mass (ANN2). Then, we calculated trunk segmental mass for two volunteers using the predicted body shape and tissue mass. The body shape model (ANN1) was tested on 279 subjects, and maximum deviation between the predicted body shape and the original was 28 mm. The tissue mass model (ANN2) was evaluated on 223 subjects, which when compared to ground truth data, had a mean error of less than 0.51% in the head, trunk, legs, and arms. We also compared the two volunteer's trunk segment mass with experimental data and found similar trend and magnitude. Our findings suggested that the proposed method could serve as an effective and convenient tool for predicting trunk mass.

BlazePose-Seq2Seq: Leveraging Regular RGB Cameras for Robust Gait Assessment.

Evaluation of human gait through smartphone-based pose estimation algorithms provides an attractive alternative to costly lab-bound instrumented assessment and offers a paradigm shift with real time gait capture for clinical assessment. Systems based on smart phones, such as OpenPose and BlazePose have demonstrated potential for virtual motion assessment but still lack the accuracy and repeatability standards required for clinical viability. Seq2seq architecture offers an alternative solution to conventional deep learning techniques for predicting joint kinematics during gait. This study introduces a novel enhancement to the low-powered BlazePose algorithm by incorporating a Seq2seq autoencoder deep learning model. To ensure data accuracy and reliability, synchronized motion capture involving an RGB camera and ten Vicon cameras were employed across three distinct self-selected walking speeds. This investigation presents a groundbreaking avenue for remote gait assessment, harnessing the potential of Seq2seq architectures inspired by natural language processing (NLP) to enhance pose estimation accuracy. When comparing BlazePose alone to the combination of BlazePose and 1D convolution Long Short-term Memory Network (1D-LSTM), Gated Recurrent Unit (GRU) and Long Short-Term Memory (LSTM), the average mean absolute errors decreased from 13.4° to 5.3° for fast gait, from 16.3° to 7.5° for normal gait, and from 15.5° to 7.5° for slow gait at the left ankle joint angle respectively. The strategic utilization of synchronized data and rigorous testing methodologies further bolsters the robustness and credibility of these findings.

Muscle-inspired bi-planar cable routing: a novel framework for designing cable driven lower limb rehabilitation exoskeletons (C-LREX).

There is a growing interest in the research and development of Cable Driven Rehabilitation Devices (CDRDs) due to multiple inherent features attractive to clinical applications, including low inertia, lightweight, high payload-to-weight ratio, large workspace, and modular design. However, previous CDRDs have mainly focused on modifying motor impairment in the sagittal plane, despite the fact that neurological disorders, such as stroke, often involve postural control and gait impairment in multiple planes. To address this gap, this work introduces a novel framework for designing a cable-driven lower limb rehabilitation exoskeleton which can assist with bi-planar impaired posture and gait. The framework used a lower limb model to analyze different cable routings inspired by human muscle architecture and attachment schemes to identify optimal routing and associated parameters. The selected cable routings were safeguarded for non-interference with the human body while generating bi-directional joint moments. The subsequent optimal cable routing model was then implemented in simulations of tracking reference healthy trajectory with bi-planar impaired gait (both in the sagittal and frontal planes). The results showed that controlling joints independently via cables yielded better performance compared to dependent control. Routing long cables through intermediate hinges reduced the peak tensions in the cables, however, at a cost of induced additional joint forces. Overall, this study provides a systematic and quantitative in silico approach, featured with accessible graphical user interface (GUI), for designing subject-specific, safe, and effective lower limb cable-driven exoskeletons for rehabilitation with options for multi-planar personalized impairment-specific features.

Development and validation of a novel ankle joint musculoskeletal model.

An improved understanding of contact mechanics in the ankle joint is paramount for implant design and ankle disorder treatment. However, existing models generally simplify the ankle joint as a revolute joint that cannot predict contact characteristics. The current study aimed to develop a novel musculoskeletal ankle joint model that can predict contact in the ankle joint, together with muscle and joint reaction forces. We modelled the ankle joint as a multi-axial joint and simulated contact mechanics between the tibia, fibula and talus bones in OpenSim. The developed model was validated with results from experimental studies through passive stiffness and contact. Through this, we found a similar ankle moment-rotation relationship and contact pattern between our study and experimental studies. Next, the musculoskeletal ankle joint model was incorporated into a lower body model to simulate gait. The ankle joint contact characteristics, kinematics, and muscle forces were predicted and compared to the literature. Our results revealed a comparable peak contact force and the same muscle activation patterns in four major muscles. Good agreement was also found in ankle dorsi/plantar-flexion and inversion/eversion. Thus, the developed model was able to accurately model the ankle joint and can be used to predict contact characteristics in gait.

Design of a lightweight universal talus implant using topology optimization.

Total talus replacement is a promising alternative treatment for talus fractures complicated by avascular necrosis and collapse. This surgical option replaces the human talus bone with a customized talus implant and can maintain ankle joint functionality compared to traditional treatment (e.g., ankle fusion). However, the customized implant is costly and time-consuming due to its customized nature. To circumvent these drawbacks, universal talus implants were proposed. While they showed clinically satisfactory results, existing talus implants are heavier than biological talus bones as they are solid inside. This can lead to unequal weight between the implant and biological talus bone, and therefore leading to other complications. The reduction of the implants' weight without compromising its performance and congruency with surrounding bones is a potential solution. Therefore, this study aims to design a lightweight universal talus implant using topology optimization. This is done through establishing the loading and boundary conditions for three common foot postures: neutral, dorsi- and plantar-flexion. The optimized implant performance in terms of mass, contact characteristics with surrounding joint cartilage and stress distributions is studied using a 3D Finite Element (FE) model of the ankle joint. The mass of the optimized implant is reduced by approximately 66.6% and its maximum stresses do not exceed 70 MPa, resulting in a safety factor of 15.7. Moreover, the optimized and solid implants show similar contact characteristics. Both implants produced peak contact pressures that were approximately 19.0%-196% higher than those produced by the biological talus. While further mechanical testing under in-vivo loading conditions is required to determine clinical feasibility, preliminarily, the use of a lightweight universal implant is expected to provide the patient with a more natural feel, and a reduced waiting period until surgery.

Sensitivity of subject-specific upper body musculoskeletal model predictions to mass scaling methods.

Accurate predictions of spinal loads in subject-specific musculoskeletal models require precise body segment parameters, including segment mass and center of mass (CoM) locations. Existing upper body models often assume a constant percentage of total body mass to calculate segmental masses, disregarding inter-individual variability and limiting their predictive capacity. This study evaluated the sensitivity of subject-specific upper body musculoskeletal model predictions to body mass scaling methods. The upper body segmental masses and corresponding CoM of six male subjects with varying body mass indices were computed using two mass scaling methods: the constant-percentage-based (CPB) scaling method, commonly used in AnyBody software; and our recently developed body-shape-based (BSB) method. Subsequently, these values were used by a validated musculoskeletal model to predict the muscle and disc forces in upright and flexed postures. The discrepancies between the results of the two scaling methods were compared across subjects and postures. Maximum deviations in thorax masses reached up to 7.5% of total body weight (TBW) in overweight subjects, while maximum CoM location differences of up to 35 mm were observed in normal weight subjects. The root mean squared errors (RMSE) of the CPB results, calculated with the BSB results as baseline, showed that the muscle and shear forces of the two scaling methods were quite similar (<4.5% of TBW). Though, there were small to moderate differences in compressive forces (6.5-16.0% of TBW). Thus, the compressive forces predicted with CPB method should be used with caution, particularly for overweight and obese subjects.

Bi-Planar Trajectory Tracking with a Novel 3DOF Cable Driven Lower Limb Rehabilitation Exoskeleton (C-LREX).

Although Cable-driven rehabilitation devices (CDRDs) have several advantages over traditional link-driven devices, including their light weight, ease of reconfiguration, and remote actuation, the majority of existing lower-limb CDRDs are limited to rehabilitation in the sagittal plane. In this work, we proposed a novel three degrees of freedom (DOF) lower limb model which accommodates hip abduction/adduction motion in the frontal plane, as well as knee and hip flexion/extension in the sagittal plane. The proposed model was employed to investigate the feasibility of using bi-planar cable routing to track a bi-planar reference healthy trajectory. Various possible routings of four cable configurations were selected and studied with the 3DOF model. The optimal locations of the hip cuffs were determined using optimization. When compared with the five-cable routing configuration, the four-cable routing produced higher joint forces, which motivated the future study of other potential cable routing configurations and their ability to track bi-planar motion.

A Framework for Determining the Performance and Requirements of Cable-Driven Mobile Lower Limb Rehabilitation Exoskeletons.

The global increase in the number of stroke patients and limited accessibility to rehabilitation has promoted an increase in the design and development of mobile exoskeletons. Robot-assisted mobile rehabilitation is rapidly emerging as a viable tool as it could provide intensive repetitive movement training and timely standardized delivery of therapy as compared to conventional manual therapy. However, the majority of existing lower limb exoskeletons continue to be heavy and induce unnecessary inertia and inertial vibration on the limb. Cable-driven exoskeletons can overcome these issues with the provision of remote actuation. However, the number of cables and routing can be selected in various ways posing a challenge to designers regarding the optimal design configuration. In this work, a simulation-based generalized framework for modelling and assessment of cable-driven mobile exoskeleton is proposed. The framework can be implemented to identify a 'suitable' configuration from several potential ones or to identify the optimal routing parameters for a given configuration. For a proof of concept, four conceptual configurations of cable-driven exoskeletons (one with a spring) were developed in a manner where both positive and negative moments could be generated for each joint (antagonistic configuration). The models were analyzed using the proposed framework and a decision metric table has been developed based on the models' performance and requirements. The weight of the metrics can be adjusted depending on the preferences and specified constraints. The maximum score is assigned to the configuration with minimum requirement or error, maximum performance, and vice versa. The metric table indicated that the 4-cable configuration is a promising design option for a lower limb rehabilitation exoskeleton based on tracking performance, model requirements, and component forces exerted on the limb.

Sensitivity of the Cervical Disc Loads, Translations, Intradiscal Pressure, and Muscle Activity Due to Segmental Mass, Disc Stiffness, and Muscle Strength in an Upright Neutral Posture.

Musculoskeletal disorders of the cervical spine have increased considerably in recent times. To understand the effects of various biomechanical factors, quantifying the differences in disc loads, motion, and muscle force/activity is necessary. The kinematic, kinetic, or muscle response may vary in a neutral posture due to interindividual differences in segmental mass, cervical disc stiffness, and muscle strength. Therefore, our study aimed to develop an inverse dynamic model of the cervical spine, estimate the differences in disc loads, translations, intradiscal pressure, and muscle force/activity in a neutral posture and compare these results with data available in the literature. A head-neck complex with nine segments (head, C1-T1) was developed with joints having three rotational and three translational degrees of freedom, 517 nonlinear ligament fibers, and 258 muscle fascicles. A sensitivity analysis was performed to calculate the effect of segmental mass (5th to 95th percentile), translational disc stiffness (0.5-1.5), and muscle strength (0.5-1.5) on the cervical disc loads (C2-C3 to C7-T1), disc translations, intradiscal pressure, and muscle force/activity in a neutral posture. In addition, two axial external load conditions (0 and 40 N) were also considered on the head. The estimated intradiscal pressures (0.2-0.56 MPa) at 0 N axial load were comparable to in vivo measurements found in the literature, whereas at 40 N, the values were 0.39-0.93 MPa. With increased segmental mass (5th to 95th), the disc loads, translations, and muscle forces/activities increased to 69% at 0 N and 34% at 40 N axial load. With increased disc stiffness (0.5-1.5), the maximum differences in axial (<1%) and shear loads (4%) were trivial; however, the translations were reduced by 67%, whereas the differences in individual muscle group forces/activities varied largely. With increased muscle strength (0.5-1.5), the muscle activity decreased by 200%. For 40 vs. 0 N, the differences in disc loads, translations, and muscle forces/activities were in the range of 52-129%. Significant differences were estimated in disc loads, translations, and muscle force/activity in the normal population, which could help distinguish between normal and pathological cervical spine conditions.

The evaluation of artificial talus implant on ankle joint contact characteristics: a finite element study based on four subjects.

Customized talus implants have been regarded as a better treatment alternative to talus avascular necrosis than traditional surgical fusion because of its ability to maintain joint mobility while ameliorating pain. Despite the use of ankle hemiarthroplasty clinically, the cartilage contact characteristics of adjacent bones remain unclear. This study aims to use finite element modeling to evaluate the contact characteristics of three types of cobalt-chrome talus implants in three postures, in four subjects. This study also compared the contact area, contact pressure, and peak contact pressure of the implant models with a reference biological model. Among the various biological and implant models, our results showed that the biological models generally had the largest contact areas and smallest peak contact pressures, whereas the implant-type models had smaller contact areas and relatively larger peak contact pressure. Moreover, among the three implant types, customized-scale models showed a larger total contact area than that of the SSM-scale and universal-scale models, but their variation was relatively limited. The results from this study can have significance in future endeavors into ankle joint modeling, as well as being able to improve implant design to enhance recovery outcomes for patients who may benefit from talar replacement.

Polycarbonate-urethane coating can significantly improve talus implant contact characteristics.

Talus implants can be utilized in cases of talus avascular necrosis and has been regarded as a promising treatment method. However, existing implants are made of stiff materials that directly oppose natural cartilage. The risk of long-term cartilage wear and bone fracture from the interaction between the cartilage and stiff implant surfaces has been documented in post-hemiarthroplasty of the hip, knee and ankle joints. The aim is to explore the effects of adding a layer of compliant material (polycarbonate-urethane; PCU) over a stiff material (cobalt chromium) in talus implants. To do so, we obtained initial ankle geometry from four cadaveric subjects in neutral standing to create the finite element models. We simulated seven models for each subject: three different types of talus implants, each coated with and without PCU, and a biological model. In total, we constructed 28 finite element models. By comparing the contact characteristics of the implant models with their respective biological model counterparts, our results showed that PCU coated implants have comparable contact area and contact pressure to the biological models, whereas stiff material implants without the PCU coating all have relatively higher contact pressure and smaller contact areas. These results confirmed that adding a layer of compliant material coating reduces the contact pressure and increases the contact area which in turn reduces the risk of cartilage wear and bone fracture. The results also suggest that there can be clinical benefits of adding a layer of compliant material coating on existing stiff material implants, and can provide valuable information towards the design of more biofidelic implants in the future.

Present and future of gait assessment in clinical practice: Towards the application of novel trends and technologies.

Background: Despite being available for more than three decades, quantitative gait analysis remains largely associated with research institutions and not well leveraged in clinical settings. This is mostly due to the high cost/cumbersome equipment and complex protocols and data management/analysis associated with traditional gait labs, as well as the diverse training/experience and preference of clinical teams. Observational gait and qualitative scales continue to be predominantly used in clinics despite evidence of less efficacy of quantifying gait. Research objective: This study provides a scoping review of the status of clinical gait assessment, including shedding light on common gait pathologies, clinical parameters, indices, and scales. We also highlight novel state-of-the-art gait characterization and analysis approaches and the integration of commercially available wearable tools and technology and AI-driven computational platforms. Methods: A comprehensive literature search was conducted within PubMed, Web of Science, Medline, and ScienceDirect for all articles published until December 2021 using a set of keywords, including normal and pathological gait, gait parameters, gait assessment, gait analysis, wearable systems, inertial measurement units, accelerometer, gyroscope, magnetometer, insole sensors, electromyography sensors. Original articles that met the selection criteria were included. Results and significance: Clinical gait analysis remains highly observational and is hence subjective and largely influenced by the observer's background and experience. Quantitative Instrumented gait analysis (IGA) has the capability of providing clinicians with accurate and reliable gait data for diagnosis and monitoring but is limited in clinical applicability mainly due to logistics. Rapidly emerging smart wearable technology, multi-modality, and sensor fusion approaches, as well as AI-driven computational platforms are increasingly commanding greater attention in gait assessment. These tools promise a paradigm shift in the quantification of gait in the clinic and beyond. On the other hand, standardization of clinical protocols and ensuring their feasibility to map the complex features of human gait and represent them meaningfully remain critical challenges.

Computational modelling of hip resurfacing arthroplasty investigating the effect of femoral version on hip biomechanics.

AIM: How reduced femoral neck anteversion alters the distribution of pressure and contact area in Hip Resurfacing Arthroplasty (HRA) remains unclear. The purpose of this study was to quantitatively describe the biomechanical implication of different femoral neck version angles on HRA using a finite element analysis. MATERIALS AND METHODS: A total of sixty models were constructed to assess the effect of different femoral neck version angles on three different functional loads: 0°of hip flexion, 45°of hip flexion, and 90° of hip flexion. Femoral version was varied between 30° of anteversion to 30° of retroversion. All models were tested with the acetabular cup in four different positions: (1) 40°/15° (inclination/version), (2) 40°/25°, (3) 50°/15°, and (4) 50°/25°. Differences in range of motion due to presence of impingement, joint contact pressure, and joint contact area with different femoral versions and acetabular cup positions were calculated. RESULTS: Impingement was found to be most significant with the femur in 30° of retroversion, regardless of acetabular cup position. Anterior hip impingement occurred earlier during hip flexion as the femur was progressively retroverted. Impingement was reduced in all models by increasing acetabular cup inclination and anteversion, yet this consequentially led to higher contact pressures. At 90° of hip flexion, contact pressures and contact areas were inversely related and showed most notable change with 30° of femoral retroversion. In this model, the contact area migrated towards the anterior implant-bone interface along the femoral neck. CONCLUSION: Femoral retroversion in HRA influences impingement and increases joint contact pressure most when the hip is loaded in flexion. Increasing acetabular inclination decreases the area of impingement but doing so causes a reciprocal increase in joint contact pressure. It may be advisable to study femoral neck version pre-operatively to better choose hip resurfacing arthroplasty candidates.

Prediction of human male trunk mass distribution using anthropometric measurements: A feasibility study.

The current study proposes a new method to predict the body shape and mass distribution of the trunk (Tl-L5) of a human male using 15 anthropometric measurements acquired at various locations of the body. Trunk cross-sectional images adopted from the Visible Human male project database were segmented into fat, bone, and lean tissue. Assuming that all male subjects have similar cross-sectional composition at a given body height percentile, areas of the segmented cross-sectional images of the Visible Human male along the trunk were scaled to match those of the predicted body shape. The trunk mass distribution of the target subject can then be computed using the density values of fat, bone, and lean tissue. Comparison of the predicted body shape circumference with ground truth values measured using digital and actual measurements yielded maximum mean error of 13.3 mm and 30.3 mm, respectively. The accuracy of the image segmentation was evaluated, and the results showed a high Jaccard index (>0.95). The proposed method was able to predict the trunk mass distribution of two volunteers with a maximum deviation of 384 g at T4 level and a minimum deviation of 12 g at L4 level and the corresponding centers of mass fell within the experimental data at most levels. Thus, our method can be considered as a feasible option to calculate subject-specific trunk mass distribution.

The importance of intervertebral disc material model on the prediction of mechanical function of the cervical spine.

BACKGROUND: Linear elastic, hyperelastic, and multiphasic material constitutive models are frequently used for spinal intervertebral disc simulations. While the characteristics of each model are known, their effect on spine mechanical response requires a careful investigation. The use of advanced material models may not be applicable when material constants are not available, model convergence is unlikely, and computational time is a concern. On the other hand, poor estimations of tissue's mechanical response are likely if the spine model is oversimplified. In this study, discrepancies in load response introduced by material models will be investigated. METHODS: Three fiber-reinforced C2-C3 disc models were developed with linear elastic, hyperelastic, and biphasic behaviors. Three different loading modes were investigated: compression, flexion and extension in quasi-static and dynamic conditions. The deformed disc height, disc fluid pressure, range of motion, and stresses were compared. RESULTS: Results indicated that the intervertebral disc material model has a strong effect on load-sharing and disc height change when compression and flexion were applied. The predicted mechanical response of three models under extension had less discrepancy than its counterparts under flexion and compression. The fluid-solid interaction showed more relevance in dynamic than quasi-static loading conditions. The fiber-reinforced linear elastic and hyperelastic material models underestimated the load-sharing of the intervertebral disc annular collagen fibers. CONCLUSION: This study confirmed the central role of the disc fluid pressure in spinal load-sharing and highlighted loading conditions where linear elastic and hyperelastic models predicted energy distribution different than that of the biphasic model.

Analysis of four methods of measuring three-dimensional pelvic tilt in the lateral decubitus position.

Proper alignment of acetabular implantation is necessary for good patient outcomes and preventing complications or additional surgeries in total hip arthroplasty (THA). Rotation of the pelvis in lateral decubitus (LD) is typically not accounted for as surgeons use the surgery table as a reference plane to align implants. This study compared four techniques to measure 3D pelvic rotations in LD position using computer models. CT scans of 19 subjects in LD position were used to create 3D pelvis models. Pelvic rotations were measured by three users using four techniques: landmark (LM), defined plane (DP), anterior pelvic plane (APP), and plane of best symmetry (POBS) methods. Measurements were analysed for intra-user reliability and relationships between methods were assessed using intraclass correlation coefficient, standard error of measurement, and coefficient of determination. The POBS method exhibited the highest inter-observer reliability and is recommended as a new measurement technique to measure pelvic rotations before THA surgery. The LM method exhibited low reliability but may be suitable for expert users familiar with pelvic landmarks. The APP method should only be used when the patients' APP planes are known in neutral standing and side-lying postures. The DP method is not recommended due to substantial individual variability. Graphical abstract.

Investigation of the Average Shape and Principal Variations of the Human Talus Bone Using Statistic Shape Model.

Due to the complexity of articular interconnections and tenuous blood supply to the talus, talus fractures are often associated with complications (e.g., avascular necrosis). Currently, surgically fusing the talus to adjacent bones is widely used as treatment to talus fractures, but this procedure can greatly reduce mobility in the ankle and hindfoot. Alternatively, customized talus implants have shown an overall satisfactory patient feedback but with the limitation of high expenses and time-consuming manufacturing process. In order to circumvent these disadvantages, universal talus implants have been proposed as a potential solution. In our study, we aimed to develop a methodology using Statistical Shape Model (SSM) to simulate the talus, and then evaluate the feasibility of the model to obtain the mean shape needed for universal implant design. In order to achieve this, we registered 98 tali (41 females and 57 males) and used the registered dataset to train our SSM. We used the mean shape derived from the SSM as the basis for our talus implant template, and compared our template with that of previous works. We found that our SSM mean shape talus implant was geometrically similar to implants from other works, which used a different method for the mean shape. This suggests the feasibility of SSM as a method of finding mean shape information for the development of universal implants. A second aim of our study was to investigate if one scalable talus implant can accommodate all patients. In our study, we focused on addressing this from a geometric perspective as there are multiple factors impacting this (e.g., articular surface contact characteristics, implant material properties). Our initial findings are that the first two principal components should be afforded consideration for the geometrical accuracy of talus implant design. Additional factors would need to be further evaluated for their role in informing universal talus implant design.

Effects of nucleus pulposus location on spinal loads and joint centers of rotation and reaction during forward flexion: A combined finite element and Musculoskeletal study.

The mechanical environment of the intervertebral disc and spinal loads are intimately associated with low back pain (LBP) caused by mechanical load. The cause and effect relationship between postures and pain has been therapeutically investigated and widely used to guide patient care. Shift of the nucleus pulposus (NP) inside the intervertebral disc when the spine changes posture in the sagittal plane, known also as NP migration, has been observed and quantified in in-vivo studies. However, the effects of this phenomenon on biomechanics of the disc and the surrounding spinal structures has not been investigated. This study aimed to quantify the effects of the NP location in a moderate flexion posture on spinal load and load-sharing using a Finite Element (FE) model of the L4-L5 functional spinal unit (FSU) driven by muscle forces, reaction forces, ligament forces and disc moment predicted by a Musculoskeletal (MSK) model of the trunk. The results showed that intradiscal pressure (IDP) and compressive force are sensitive to the nucleus location while the effects on the center of rotation (CoR), center of reaction (CTR) and moment rotation curves were negligible. Also, our findings revealed that FE models should consider the effects of NP location during bending to predict more realistic results as the nucleus displacement caused by disc bulge predicted by these models is much smaller than the real shift observed in in-vivo. In addition, this study confirmed that position of the rigid joint in MSK models that fix this latter to the CoR, must be modelled carefully for more accurate muscle forces and spinal loads prediction.