Multibody dynamics-based musculoskeletal modeling for gait analysis: a systematic review.

Beyond qualitative assessment, gait analysis involves the quantitative evaluation of various parameters such as joint kinematics, spatiotemporal metrics, external forces, and muscle activation patterns and forces. Utilizing multibody dynamics-based musculoskeletal (MSK) modeling provides a time and cost-effective non-invasive tool for the prediction of internal joint and muscle forces. Recent advancements in the development of biofidelic MSK models have facilitated their integration into clinical decision-making processes, including quantitative diagnostics, functional assessment of prosthesis and implants, and devising data-driven gait rehabilitation protocols. Through an extensive search and meta-analysis of over 116 studies, this PRISMA-based systematic review provides a comprehensive overview of different existing multibody MSK modeling platforms, including generic templates, methods for personalization to individual subjects, and the solutions used to address statically indeterminate problems. Additionally, it summarizes post-processing techniques and the practical applications of MSK modeling tools. In the field of biomechanics, MSK modeling provides an indispensable tool for simulating and understanding human movement dynamics. However, limitations which remain elusive include the absence of MSK modeling templates based on female anatomy underscores the need for further advancements in this area.

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.

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.

Hip Joint Contact Pressure Distribution During Pavlik Harness Treatment of an Infant Hip: A Patient-Specific Finite Element Model.

Developmental dysplasia of the hip (DDH) in infants under 6 months of age is typically treated by the Pavlik harness (PH). During successful PH treatment, a subluxed/dislocated hip is spontaneously reduced into the acetabulum, and DDH undergoes self-correction. PH treatment may fail due to avascular necrosis (AVN) of the femoral head. An improved understanding of mechanical factors accounting for the success/failure of PH treatment may arise from investigating articular cartilage contact pressure (CCP) within a hip during treatment. In this study, CCP in a cartilaginous infant hip was investigated through patient-specific finite element (FE) modeling. We simulated CCP of the hip equilibrated at 90 deg flexion at abduction angles of 40 deg, 60 deg, and 80 deg. We found that CCP was predominantly distributed on the anterior and posterior acetabulum, leaving the superior acetabulum (mainly superolateral) unloaded. From a mechanobiological perspective, hypothesizing that excessive pressure inhibits growth, our results qualitatively predicted increased obliquity and deepening of the acetabulum under such CCP distribution. This is the desired and observed therapeutic effect in successful PH treatment. The results also demonstrated increase in CCP as abduction increased. In particular, the simulation predicted large magnitude and concentrated CCP on the posterior wall of the acetabulum and the adjacent lateral femoral head at extreme abduction (80 deg). This CCP on lateral femoral head may reduce blood flow in femoral head vessels and contribute to AVN. Hence, this study provides insight into biomechanical factors potentially responsible for PH treatment success and complications.

Geometric analysis of the talus and development of a generic talar prosthetic.

BACKGROUND: Trauma to the talus can result in fracture, avascular necrosis and structural collapse. Treatment has been limited to surgical fusion and total ankle arthroplasty. Total ankle arthroplasty may not be an appropriate treatment for avascular necrosis while surgical fusion of the joint limits mobility. Custom-made implants have recently been used to address these limitations but have lengthy delays between injury and surgery and higher associated costs. A generic talar prosthesis available in various sizes may serve as a suitable alternative. METHODS: The geometric variation between shapes of individual tali was determined using 3D geometric models of 91 tali created from CT-scan data. Comparisons were done to determine if tali are one shape. The best shape was determined for each sex, and was compared to determine if a unisex implant would be possible. A geometric template for the implant in multiple sizes was created and compared to the models. RESULTS: The average of the average deviation between tali after volume scaling was found to be less than 1mm on the main articulating surfaces. One shape group was found for the talus. The female and male tali were found to be similar and a unisex implant template was created. CONCLUSIONS: Ten generic talar implant sizes were determined to be sufficient to match the size and shape of the 91 tali examined in this study.

Leg dominance may not be a predictor of asymmetry in peak joint moments and ground reaction forces during sit-to-stand movements.

Sit-to-stand transfer is a common prerequisite for many daily tasks. Literature often assumes symmetric behavior across the left and right side. Although this assumption of bilateral symmetry is prominent, few studies have validated this supposition. This pilot study uniquely quantifies peak joint moments and ground reaction forces (GRFs), using a Euclidian norm approach, to evaluate bilateral symmetry and its relation to lower limb motor-dominance during sit to stand in ten healthy males. Peak joint moments and GRFs were determined using a motion capture system and inverse dynamics. This analysis included joint moment contributions from all three body planes (sagittal, coronal, and axial) as well as vertical and shearing GRFs. A paired, one-tailed t test was used, suggesting asymmetrical joint moment development in all three lower extremity joints as well as GRFs (P < .05). Furthermore, using an unpaired two-tailed t test, asymmetry developed during these movements does not appear to be predictable by participants' lower limb motor-dominance (P < .025). Consequently, when evaluating sit-to-stand it is suggested the effects of asymmetry be considered in the interpretation of data. The absence of a relationship between dominance and asymmetry prevents the suggestion that one side can be tested to infer behavior of the contralateral.

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.

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.

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.

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.

An in vitro study on the dimensional stability of a vinyl polyether silicone impression material over a prolonged storage period.

STATEMENT OF PROBLEM: A new elastomeric impression material, a vinyl polyether silicone, has been commercially introduced. Its dimensional stability at different pouring times has not been reported. PURPOSE: The purpose of this study was to evaluate the dimensional stability of vinyl polyether silicone (VPES) impressions as a function of delayed-pouring time for up to 2 weeks after performing a recommended clinical disinfection procedure. MATERIAL AND METHODS: The 3 medium-body impression materials tested were EXA'lence 370, Imprint 3, and Impregum Penta soft. Impressions of a cylindrical metal model, which served as a control, were made and poured in a Type V stone after being disinfected in 2.5% buffered glutaraldehyde solution. Die diameter and anteroposterior and cross-arch measurements on each cast were made and compared to direct measurements of the control with a digital micrometer. The linear dimensional changes were compared and analyzed by using the Mann-Whitney U test (α=.05). RESULTS: Considering the absolute values for the mean percentage dimensional change, VPES cast measurements were below 1.0% (P<.001), with the majority being minimal (≤0.34%) at all pour times. Measurements of the casts were larger than those of the control in all but one specimen. The percentage dimensional changes of the individual die diameters were higher than those of the anteroposterior and cross-arch linear measurements. Casts produced from VPES impressions had similar dimensional changes to those of VPS at 1 week of storage and similar changes to those of PE at 2 weeks of storage. The behavior of VPES impressions varied from being the most accurate at immediate-pour to gradually exhibiting higher dimensional changes as time progressed. CONCLUSIONS: Casts produced from a disinfected regular set VPES (EXA'lence 370 monophase) demonstrated excellent dimensional stability at different pour times and were comparable to the tested VPS and PE impression materials.

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.

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 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.

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.

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.

Calibration of hyperelastic material properties of the human lumbar intervertebral disc under fast dynamic compressive loads.

Under fast dynamic loading conditions (e.g. high-energy impact), the load rate dependency of the intervertebral disc (IVD) material properties may play a crucial role in the biomechanics of spinal trauma. However, most finite element models (FEM) of dynamic spinal trauma uses material properties derived from quasi-static experiments, thus neglecting this load rate dependency. The aim of this study was to identify hyperelastic material properties that ensure a more biofidelic simulation of the IVD under a fast dynamic compressive load. A hyperelastic material law based on a first-order Mooney-Rivlin formulation was implemented in a detailed FEM of a L2-L3 functional spinal unit (FSU) to represent the mechanical behavior of the IVD. Bony structures were modeled using an elasto-plastic Johnson-Cook material law that simulates bone fracture while ligaments were governed by a viscoelastic material law. To mimic experimental studies performed in fast dynamic compression, a compressive loading velocity of 1 m/s was applied to the superior half of L2, while the inferior half of L3 was fixed. An exploratory technique was used to simulate dynamic compression of the FSU using 34 sets of hyperelastic material constants randomly selected using an optimal Latin hypercube algorithm and a set of material constants derived from quasi-static experiments. Selection or rejection of the sets of material constants was based on compressive stiffness and failure parameters criteria measured experimentally. The two simulations performed with calibrated hyperelastic constants resulted in nonlinear load-displacement curves with compressive stiffness (7335 and 7079 N/mm), load (12,488 and 12,473 N), displacement (1.95 and 2.09 mm) and energy at failure (13.5 and 14.7 J) in agreement with experimental results (6551 ± 2017 N/mm, 12,411 ± 829 N, 2.1 ± 0.2 mm and 13.0 ± 1.5 J respectively). The fracture pattern and location also agreed with experimental results. The simulation performed with constants derived from quasi-static experiments showed a failure energy (13.2 J) and a fracture pattern and location in agreement with experimental results, but a compressive stiffness (1580 N/mm), a failure load (5976 N) and a displacement to failure (4.8 mm) outside the experimental corridors. The proposed method offers an innovative way to calibrate the hyperelastic material properties of the IVD and to offer a more realistic simulation of the FSU in fast dynamic compression.