An Efficient Modelling-Simulation-Analysis Workflow to Investigate Stump-Socket Interaction Using Patient-Specific, Three-Dimensional, Continuum-Mechanical, Finite Element Residual Limb Models

Determining a musculoskeletal system's pre-stretched state using continuum-mechanical forward modelling and joint range optimization

The subject-specific range of motion (RoM) of a musculoskeletal joint system is balanced by pre-tension levels of individual muscles, which affects their contraction capability. Such an inherent pre-tension or pre-stretch of muscles is not measureable with in vivo experiments. Using a 3D continuum mechanical forward simulation approach for motion analysis of the musculoskeletal system of the forearm with 3 flexor and 2 extensor muscles, we developed an optimization process to determine the muscle fibre pre-stretches for an initial arm position, which is given human dataset. We used RoM values of a healthy person to balance the motion in extension and flexion. The performed sensitivity study shows that the fibre pre-stretches of the m. brachialis, m. biceps brachii and m. triceps brachii with 91 % dominate the objective flexion ratio, while m. brachiradialis and m. anconeus amount 7.8 % and 1.2 % . Within the multi-dimensional space of the surrogate model, 3D sub-spaces of primary variables, namely the dominant muscles and the global objective, flexion ratio, exhibit a path of optimal solutions. Within this optimal path, the muscle fibre pre-stretch of two flexors demonstrate a negative correlation, while, in contrast, the primary extensor, m. triceps brachii correlates positively to each of the flexors. Comparing the global optimum with four other designs along the optimal path, we saw large deviations, e.g., up to 15∘ in motion and up to 40% in muscle force. This underlines the importance of accurate determination of fibre pre-stretch in muscles, especially, their role in pathological muscular disorders and surgical applications such as free muscle or tendon transfer.

Finite Element Model of the Shoulder with Active Rotator Cuff Muscles: Application to Wheelchair Propulsion

The rotator cuff is prone to injury, remarkably so for manual wheelchair users. To understand its pathomechanisms, finite element models incorporating three-dimensional activated muscles are needed to predict soft tissue strains during given tasks. This study aimed to develop such a model to understand pathomechanisms associated with wheelchair propulsion. We developed an active muscle model associating a passive fiber-reinforced isotropic matrix with an activation law linking calcium ion concentration to tissue tension. This model was first evaluated against known physiological muscle behavior; then used to activate the rotator cuff during a wheelchair propulsion cycle. Here, experimental kinematics and electromyography data was used to drive a shoulder finite element model. Finally, we evaluated the importance of muscle activation by comparing the results of activated and non-activated rotator cuff muscles during both propulsion and isometric contractions. Qualitatively, the muscle constitutive law reasonably reproduced the classical Hill model force-length curve and the behavior of a transversally loaded muscle. During wheelchair propulsion, the deformation and fiber stretch of the supraspinatus muscle-tendon unit pointed towards the possibility for this tendon to develop tendinosis due to the multiaxial loading imposed by the kinematics of propulsion. Finally, differences in local stretch and positions of the lines of action between activated and non-activated models were only observed at activation levels higher than 30%. Our novel finite element model with active muscles is a promising tool for understanding the pathomechanisms of the rotator cuff for various dynamic tasks, especially those with high muscle activation levels.

Linking cortex and contraction-Integrating models along the corticomuscular pathway

Computational models of the neuromusculoskeletal system provide a deterministic approach to investigate input-output relationships in the human motor system. Neuromusculoskeletal models are typically used to estimate muscle activations and forces that are consistent with observed motion under healthy and pathological conditions. However, many movement pathologies originate in the brain, including stroke, cerebral palsy, and Parkinson's disease, while most neuromusculoskeletal models deal exclusively with the peripheral nervous system and do not incorporate models of the motor cortex, cerebellum, or spinal cord. An integrated understanding of motor control is necessary to reveal underlying neural-input and motor-output relationships. To facilitate the development of integrated corticomuscular motor pathway models, we provide an overview of the neuromusculoskeletal modelling landscape with a focus on integrating computational models of the motor cortex, spinal cord circuitry, α-motoneurons  and skeletal muscle in regard to their role in generating voluntary muscle contraction. Further, we highlight the challenges and opportunities associated with an integrated corticomuscular pathway model, such as challenges in defining neuron connectivities, modelling standardisation, and opportunities in applying models to study emergent behaviour. Integrated corticomuscular pathway models have applications in brain-machine-interaction, education, and our understanding of neurological disease.

Next-generation devices to diagnose residuum health of individuals suffering from limb loss: A narrative review of trends, opportunities, and challenges

OBJECTIVES

There is a need for diagnostic devices that can assist prosthetic care providers to better assess and maintain residuum health of individuals suffering from neuromusculoskeletal dysfunctions associated with limb loss. This paper outlines the trends, opportunities, and challenges that will facilitate the development of next-generation diagnostic devices.

DESIGN

Narrative literature review.

METHODS

Information about technologies suitable for integration into next-generation diagnostic devices was extracted from 41 references. We considered the invasiveness, comprehensiveness, and practicality of each technology subjectively.

RESULTS

This review highlighted a trend toward future diagnostic devices of neuromusculoskeletal dysfunctions of the residuum capable to support evidence-based patient-specific prosthetic care, patient empowerment, and the development of bionic solutions. This device should positively disrupt the organization healthcare by enabling cost-utility analyses (e.g., fee-for-device business models) and addressing healthcare gaps due to labor shortages. There are opportunities to develop wireless, wearable and noninvasive diagnostic devices integrating wireless biosensors to measure change in mechanical constraints and topography of residuum tissues during real-life conditions as well as computational modeling using medical imaging and finite element analysis (e.g., digital twin). Developing the next-generation diagnostic devices will require to overcome critical barriers associated with the design (e.g., gaps between technology readiness levels of essential parts), clinical roll-out (e.g., identification of primary users), and commercialization (e.g., limited interest from investors).

CONCLUSIONS

We anticipate that next-generation diagnostic devices will contribute to prosthetic care innovations that will safely increase mobility, thereby improving the quality of life of the growing global population of individuals suffering from limb loss.

A biophysically guided constitutive law of the musculotendon-complex: modelling and numerical implementation in Abaqus

BACKGROUND AND OBJECTIVE

Many biomedical, clinical, and industrial applications may benefit from musculoskeletal simulations. Three-dimensional macroscopic muscle models (3D models) can more accurately represent muscle architecture than their 1D (line-segment) counterparts. Nevertheless, 3D models remain underutilised in academic, clinical, and commercial environments. Among the reasons for this is a lack of modelling and simulation standardisation, verification, and validation. Here, we strive towards a solution by providing an open-access, characterised, constitutive relation (CR) for 3D musculotendon models.

METHODS

The musculotendon complex is modelled following the state-of-the-art active stress approach and is treated as hyperelastic, transversely isotropic, and nearly incompressible. Furthermore, force-length and -velocity relationships are incorporated, and muscle activation is derived from motor-unit information. The CR was implemented within the commercial finite-element software package Abaqus as a user-subroutine. A masticatory system model with left and right masseters was used to demonstrate active and passive movement.

RESULTS

The CR was characterised by various experimental data sets and was able to capture a wide variety of passive and active behaviours. Furthermore, the masticatory simulations revealed that joint movement was sensitive to the muscle's in-fibre passive response.

CONCLUSIONS

This user-material provides a "plug and play" template for 3D neuro-musculoskeletal finite element modelling. We hope that this reduces modelling effort, fosters exchange, and contributes to the standardisation of such models.

Société de Biomécanique Young Investigator Award 2021: Numerical investigation of the time-dependent stress-strain mechanical behaviour of skeletal muscle tissue in the context of pressure ulcer prevention

BACKGROUND

Pressure-induced tissue strain is one major pathway for Pressure Ulcer development and, especially, Deep Tissue Injury. Biomechanical investigation of the time-dependent stress-strain mechanical behaviour of skeletal muscle tissue is therefore essential. In the literature, a viscoelastic formulation is generally assumed for the experimental characterization of skeletal muscles, with the limitation that the underlying physical mechanisms that give rise to the time dependent stress-strain behaviour are not known. The objective of this study is to explore the capability of poroelasticity to reproduce the apparent viscoelastic behaviour of passive muscle tissue under confined compression.

METHODS

Experimental stress-relaxation response of 31 cylindrical porcine samples tested under fast and slow confined compression by Vaidya and collaborators were used. An axisymmetric Finite Element model was developed in ABAQUS and, for each sample a one-to-one inverse analysis was performed to calibrate the specimen-specific constitutive parameters, namely, the drained Young's modulus, the void ratio, hydraulic permeability, the Poisson's ratio, the solid grain's and fluid's bulk moduli.

FINDINGS

The peak stress and consolidation were recovered for most of the samples (N=25) by the poroelastic model (normalised root-mean-square error ≤0.03 for fast and slow confined compression conditions).

INTERPRETATION

The strength of the proposed model is its fewer number of variables (N=6 for the proposed poroelastic model versus N=18 for the viscohyperelastic model proposed by Vaidya and collaborators). The incorporation of poroelasticity to clinical models of Pessure Ulcer formation could lead to more precise and mechanistic explorations of soft tissue injury risk factors.

Analyzing Dynamic Operational Conditions of Limb Prosthetic Sockets with a Mechatronics-Twin Framework

Lower limb prostheses offer a solution to restore the ambulation and self-esteem of amputees. One key component is the prosthetic socket that serves as the interface between prosthetic device and amputee stump and thereby has a wide range of impacts on efficient fitting, appropriate load transmission, operational stability, and control. For the design and optimization of a prosthetic socket, an understanding of the actual intra-socket operational conditions becomes therefore necessary. This is however a difficult task due to the inherent complexity and restricted observability of socket operation. In this study, an innovative mechatronics-twin framework that integrates advanced biomechanical models and simulations with physical prototyping and dynamic operation testing for effective exploration of operational behaviors of prosthetic sockets with amputees is proposed. Within this framework, a specific Stewart manipulator is developed to enable dynamic operation testing, in particular for a well-managed generation of dynamic intra-socket loads and behaviors that are otherwise difficult to observe or realize with the real amputees. A combination of deep learning and Bayesian Inference algorithms is then employed for analyzing the intra-socket load conditions and revealing possible anomalous.

A Physiology-Guided Classification of Active-Stress and Active-Strain Approaches for Continuum-Mechanical Modeling of Skeletal Muscle Tissue

The well-established sliding filament and cross-bridge theory explain the major biophysical mechanism responsible for a skeletal muscle's active behavior on a cellular level. However, the biomechanical function of skeletal muscles on the tissue scale, which is caused by the complex interplay of muscle fibers and extracellular connective tissue, is much less understood. Mathematical models provide one possibility to investigate physiological hypotheses. Continuum-mechanical models have hereby proven themselves to be very suitable to study the biomechanical behavior of whole muscles or entire limbs. Existing continuum-mechanical skeletal muscle models use either an active-stress or an active-strain approach to phenomenologically describe the mechanical behavior of active contractions. While any macroscopic constitutive model can be judged by it's ability to accurately replicate experimental data, the evaluation of muscle-specific material descriptions is difficult as suitable data is, unfortunately, currently not available. Thus, the discussions become more philosophical rather than following rigid methodological criteria. Within this work, we provide a extensive discussion on the underlying modeling assumptions of both the active-stress and the active-strain approach in the context of existing hypotheses of skeletal muscle physiology. We conclude that the active-stress approach resolves an idealized tissue transmitting active stresses through an independent pathway. In contrast, the active-strain approach reflects an idealized tissue employing an indirect, coupled pathway for active stress transmission. Finally the physiological hypothesis that skeletal muscles exhibit redundant pathways of intramuscular stress transmission represents the basis for considering a mixed-active-stress-active-strain constitutive framework.

A Sensor-Based Decision Support System for Transfemoral Socket Rectification

A decision support system (DSS) was developed that outputs suggestions for socket-rectification actions to the prosthetist, aiming at improving the fitness of transfemoral prosthetic socket design and reducing the time needed for the final socket design. For this purpose, the DSS employs a fuzzy-logic inference engine (IE) which combines a set of rectification rules with pressure measurements generated by sensors embedded in the socket, for deciding the rectification actions. The latter is then processed by an algorithm that receives, manipulates and modifies a 3D digital socket model as a triangle mesh formatted inside an STL file. The DSS results were validated and tested in an FEA simulation environment, by simulating and comparing the donning process among a good-fitting socket, a loose socket (poor-fit) and several rectified sockets produced by the proposed DSS. The simulation results indicate that volume reduction improves the pressure distribution over the stump. However, as the intensity of socket rectification increases, i.e., as volume reduction increases, high pressures appear in other parts of the socket which generate discomfort. Therefore, a trade-off is required between the amount of rectification and the balance of the pressure distributions experienced at the stump.

Key considerations for finite element modelling of the residuum-prosthetic socket interface

BACKGROUND

Finite element modelling has long been proposed to support prosthetic socket design. However, there is minimal detail in the literature to inform practice in developing and interpreting these complex, highly nonlinear models.

OBJECTIVES

To identify best practice recommendations for finite element modelling of lower limb prosthetics, considering key modelling approaches and inputs.

STUDY DESIGN

Computational modelling.

METHODS

This study developed a parametric finite element model using magnetic resonance imaging data from a person with transtibial amputation. Comparative analyses were performed considering socket loading methods, socket-residuum interface parameters and soft tissue material models from the literature, to quantify their effect on the residuum's biomechanical response to a range of parameterised socket designs.

RESULTS

These variables had a marked impact on the finite element model's predictions for limb-socket interface pressure and soft tissue shear distribution.

CONCLUSIONS

All modelling decisions should be justified biomechanically and clinically. In order to represent the prosthetic loading scenario in silico, researchers should (1) consider the effects of donning and interface friction to capture the generated soft tissue shear stresses, (2) use representative stiffness hyperelastic material models for soft tissues when using strain to predict injury and (3) interrogate models comparatively, against a clinically-used control.

Tumor ablation due to inhomogeneous anisotropic diffusion in generic three-dimensional topologies

In recent decades computer-aided technologies have become prevalent in medicine, however, cancer drugs are often only tested on in vitro cell lines from biopsies. We derive a full three-dimensional model of inhomogeneous -anisotropic diffusion in a tumor region coupled to a binary population model, which simulates in vivo scenarios faster than traditional cell-line tests. The diffusion tensors are acquired using diffusion tensor magnetic resonance imaging from a patient diagnosed with glioblastoma multiform. Then we numerically simulate the full model with finite element methods and produce drug concentration heat maps, apoptosis hotspots, and dose-response curves. Finally, predictions are made about optimal injection locations and volumes, which are presented in a form that can be employed by doctors and oncologists.

Exploring the role of transtibial prosthetic use in deep tissue injury development: a scoping review

Background

The soft tissue of the residual limb in transtibial prosthetic users encounters unique biomechanical challenges. Although not intended to tolerate high loads and deformation, it becomes a weight-bearing structure within the residuum-prosthesis-complex. Consequently, deep soft tissue layers may be damaged, resulting in Deep Tissue Injury (DTI). Whilst considerable effort has gone into DTI research on immobilised individuals, only little is known about the aetiology and population-specific risk factors in amputees. This scoping review maps out and critically appraises existing research on DTI in lower-limb prosthetic users according to (1) the population-specific aetiology, (2) risk factors, and (3) methodologies to investigate both.

Results

A systematic search within the databases Pubmed, Ovid Excerpta Medica, and Scopus identified 16 English-language studies. The results indicate that prosthetic users may be at risk for DTI during various loading scenarios. This is influenced by individual surgical, morphological, and physiological determinants, as well as the choice of prosthetic componentry. However, methodological limitations, high inter-patient variability, and small sample sizes complicate the interpretation of outcome measures. Additionally, fundamental research on cell and tissue reactions to dynamic loading and on prosthesis-induced alterations of the vascular and lymphatic supply is missing.

Conclusion

We therefore recommend increased interdisciplinary research endeavours with a focus on prosthesis-related experimental design to widen our understanding of DTI. The results have the potential to initiate much-needed clinical advances in surgical and prosthetic practice and inform future pressure ulcer classifications and guidelines.

Developing an Analogue Residual Limb for Comparative DVC Analysis of Transtibial Prosthetic Socket Designs

Personalised prosthetic sockets are fabricated by expert clinicians in a skill- and experience-based process, with research providing tools to support evidence-based practice. We propose that digital volume correlation (DVC) may offer a deeper understanding of load transfer from prosthetic sockets into the residual limb, and tissue injury risk. This study’s aim was to develop a transtibial amputated limb analogue for volumetric strain estimation using DVC, evaluating its ability to distinguish between socket designs. A soft tissue analogue material was developed, comprising silicone elastomer and sand particles as fiducial markers for image correlation. The material was cast to form an analogue residual limb informed by an MRI scan of a person with transtibial amputation, for whom two polymer check sockets were produced by an expert prosthetist. The model was micro-CT scanned according to (i) an unloaded noise study protocol and (ii) a case study comparison between the two socket designs, loaded to represent two-legged stance. The scans were reconstructed to give 108 µm voxels. The DVC noise study indicated a 64 vx subvolume and 50% overlap, giving better than 0.32% strain sensitivity, and ~3.5 mm spatial resolution of strain. Strain fields induced by the loaded sockets indicated tensile, compressive and shear strain magnitudes in the order of 10%, with a high signal:noise ratio enabling distinction between the two socket designs. DVC may not be applicable for socket design in the clinical setting, but does offer critical 3D strain information from which existing in vitro and in silico tools can be compared and validated to support the design and manufacture of prosthetic sockets, and enhance the biomechanical understanding of the load transfer between the limb and the prosthesis.

Predictive prosthetic socket design: part 1-population-based evaluation of transtibial prosthetic sockets by FEA-driven surrogate modelling

It has been proposed that finite element analysis can complement clinical decision making for the appropriate design and manufacture of prosthetic sockets for amputees. However, clinical translation has not been achieved, in part due to lengthy solver times and the complexity involved in model development. In this study, a parametric model was created, informed by variation in (i) population-driven residuum shape morphology, (ii) soft tissue compliance and (iii) prosthetic socket design. A Kriging surrogate model was fitted to the response of the analyses across the design space enabling prediction for new residual limb morphologies and socket designs. It was predicted that morphological variability and prosthetic socket design had a substantial effect on socket-limb interfacial pressure and shear conditions as well as sub-dermal soft tissue strains. These relationships were investigated with a higher resolution of anatomical, surgical and design variability than previously reported, with a reduction in computational expense of six orders of magnitude. This enabled real-time predictions (1.6 ms) with error vs the analytical solutions of < 4 kPa in pressure at residuum tip, and < 3% in soft tissue strain. As such, this framework represents a substantial step towards implementation of finite element analysis in the prosthetics clinic.

Influence of Gait Cycle Loads on Stress Distribution at The Residual Limb/Socket Interface of Transfemoral Amputees: A Finite Element Analysis

A Finite Element Analysis (FEA) was performed to evaluate the interaction between residual limb and socket when considering the dynamic loads of the gait cycle. Fourteen transfemoral amputees participated in this study, where their residual limbs (i.e., soft tissues and bone), and their sockets were reconstructed. The socket and the femur were defined as elastic materials, while the bulk soft tissues were defined as a hyperelastic material. Each model included the donning, standing, and gait cycle phase, with load and boundary conditions applied accordingly. The influence of adding the dynamic loads related to the gait cycle were compared against the modelling of the static load equivalent to the standing position resulting in changes of 23% ± 19% in the maximum values and in an increase in the size of the regions where they were located. Additionally, the possible correspondence between comfort and the location of peak loadbearing at the residual-limb/socket interface was explored. Consequently, the comfort perceived by the patient could be estimated based on the locations of the maximum stresses (i.e., if they coincide with the pressure tolerant or sensitive regions of the residual limb).

Analysis of Pressure Distribution in Transfemoral Prosthetic Socket for Prefabrication Evaluation via the Finite Element Method

In this study, we estimated and validated the pressure distribution profile between the residuum and two types of prosthetic sockets for transfemoral amputees by utilizing a finite element analysis. Correct shaping of the socket for an appropriate load distribution is a critical process in the design of lower-limb prosthesis sockets. The pressure distribution profile provides an understanding of the relationship between the socket design and the level of subject comfortability. Estimating the pressure profile is important, as it helps improve the prosthesis through an evaluation of the socket design before it undergoes the fabrication process. This study focused on utilizing a magnetic resonance imaging (MRI)-based three-dimensional (3D) model inside a predetermined finite element simulation. The simulation was predetermined by mimicking the actual socket-fitting environment. The results showed that the potential MRI-based 3D model simulation could be used as an estimation tool for a pressure distribution profile due to the high correlation coefficient value (R² > 0.8) calculated when the pressure profiles were compared to the experiment data. The simulation also showed that the pressure distribution in the proximal area was higher (~30%) than in the distal area of the prosthetic socket for every subject. The results of this study will be of tremendous interest for fabricators through the use of a finite element model as an alternative method for the prefabrication and evaluation of prosthetic sockets. In future prosthetic socket fabrications, less intervention will be required in the development of a socket, and the participation of the subject in the socket-fitting session will not be necessary. The results suggest that this study will contribute to expanding the development of an overall prefabrication evaluation system to allow healthcare providers and engineers to simulate the fit and comfort of transfemoral prosthetics.

Modelling the electrical activity of skeletal muscle tissue using a multi-domain approach

Electromyography (EMG) can be used to study the behaviour of the motor neurons and thus provides insights into the physiology of the central nervous system. However, due to the high complexity of neuromuscular control, EMG signals are challenging to interpret. While the exact knowledge of the excitation patterns of a specific muscle within an in vivo experimental setting remains elusive, simulations allow to systematically investigate EMG signals in a controlled environment. Within this context, simulations can provide virtual EMG data, which, for example, can be used to validate and optimise signal analysis methods that aim to estimate the relationship between EMG signals and the output of motor neuron pools. However, since existing methods, which are employed to compute EMG signals, exhibit deficiencies with respect to the physical model itself as well as with respect to numerical aspects, we propose a novel homogenised continuum model that closely resolves the electro-physiological behaviour of skeletal muscle tissue. The proposed model is based on an extension of the well-established bidomain model and includes a biophysically detailed description of the electrical activity within the tissue, which is due to the depolarisation of the muscle fibre membranes. In contrast to all other published EMG models, which assume that the electrical potential field for each muscle fibre can be calculated independently, the proposed model assumes that the electrical potential in the muscle fibres is coupled to the electrical potential in the extracellular space. We show that the newly proposed model is able to simulate realistic EMG signals and demonstrate the potential to employ the predicted virtual EMG signal in order to evaluate the goodness of automated decomposition algorithms.

Multiscale modeling of the neuromuscular system: Coupling neurophysiology and skeletal muscle mechanics

Mathematical models and computer simulations have the great potential to substantially increase our understanding of the biophysical behavior of the neuromuscular system. This, however, requires detailed multiscale, and multiphysics models. Once validated, such models allow systematic in silico investigations that are not necessarily feasible within experiments and, therefore, have the ability to provide valuable insights into the complex interrelations within the healthy system and for pathological conditions. Most of the existing models focus on individual parts of the neuromuscular system and do not consider the neuromuscular system as an integrated physiological system. Hence, the aim of this advanced review is to facilitate the prospective development of detailed biophysical models of the entire neuromuscular system. For this purpose, this review is subdivided into three parts. The first part introduces the key anatomical and physiological aspects of the healthy neuromuscular system necessary for modeling the neuromuscular system. The second part provides an overview on state-of-the-art modeling approaches representing all major components of the neuromuscular system on different time and length scales. Within the last part, a specific multiscale neuromuscular system model is introduced. The integrated system model combines existing models of the motor neuron pool, of the sensory system and of a multiscale model describing the mechanical behavior of skeletal muscles. Since many sub-models are based on strictly biophysical modeling approaches, it closely represents the underlying physiological system and thus could be employed as starting point for further improvements and future developments. This article is categorized under: Physiology > Mammalian Physiology in Health and Disease Analytical and Computational Methods > Computational Methods Models of Systems Properties and Processes > Organ, Tissue, and Physiological Models.