Least Squares Approach to Inverse Problems in Musculoskeletal Biomechanics

Inverse simulations of musculoskeletal models computes the internal forces such as muscle and joint reaction forces, which are hard to measure, using the more easily measured motion and external forces as input data. Because of the difficulties of measuring muscle forces and joint reactions, simulations are hard to validate. One way of reducing errors for the simulations is to ensure that the mathematical problem is well-posed. This paper presents a study of regularity aspects for an inverse simulation method, often called forward dynamics or dynamical optimization, that takes into account both measurement errors and muscle dynamics. The simulation method is explained in detail. Regularity is examined for a test problem around the optimum using the approximated quadratic problem. The results shows improved rank by including a regularization term in the objective that handles the mechanical over-determinancy. Using the 3-element Hill muscle model the chosen regularization term is the norm of the activation. To make the problem full-rank only the excitation bounds should be included in the constraints. However, this results in small negative values of the activation which indicates that muscles are pushing and not pulling. Despite this unrealistic behavior the error maybe small enough to be accepted for specific applications. These results is a starting point start for achieving better results of inverse musculoskeletal simulations from a numerical point of view.

Improving Patients’ Experience and Outcome of Total Joint Replacement: the RESTORE Programme

BACKGROUND: Total hip replacements (THRs) and total knee replacements (TKRs) are common elective procedures. In the REsearch STudies into the ORthopaedic Experience (RESTORE) programme, we explored the care and experiences of patients with osteoarthritis after being listed for THR and TKR up to the time when an optimal outcome should be expected. OBJECTIVE: To undertake a programme of research studies to work towards improving patient outcomes after THR and TKR. METHODS: We used methodologies appropriate to research questions: systematic reviews, qualitative studies, randomised controlled trials (RCTs), feasibility studies, cohort studies and a survey. Research was supported by patient and public involvement. RESULTS: Systematic review of longitudinal studies showed that moderate to severe long-term pain affects about 7–23% of patients after THR and 10–34% after TKR. In our cohort study, 10% of patients with hip replacement and 30% with knee replacement showed no clinically or statistically significant functional improvement. In our review of pain assessment few research studies used measures to capture the incidence, character and impact of long-term pain. Qualitative studies highlighted the importance of support by health and social professionals for patients at different stages of the joint replacement pathway. Our review of longitudinal studies suggested that patients with poorer psychological health, physical function or pain before surgery had poorer long-term outcomes and may benefit from pre-surgical interventions. However, uptake of a pre-operative pain management intervention was low. Although evidence relating to patient outcomes was limited, comorbidities are common and may lead to an increased risk of adverse events, suggesting the possible value of optimising pre-operative management. The evidence base on clinical effectiveness of pre-surgical interventions, occupational therapy and physiotherapy-based rehabilitation relied on small RCTs but suggested short-term benefit. Our feasibility studies showed that definitive trials of occupational therapy before surgery and post-discharge group-based physiotherapy exercise are feasible and acceptable to patients. Randomised trial results and systematic review suggest that patients with THR should receive local anaesthetic infiltration for the management of long-term pain, but in patients receiving TKR it may not provide additional benefit to femoral nerve block. From a NHS and Personal Social Services perspective, local anaesthetic infiltration was a cost-effective treatment in primary THR. In qualitative interviews, patients and health-care professionals recognised the importance of participating in the RCTs. To support future interventions and their evaluation, we conducted a study comparing outcome measures and analysed the RCTs as cohort studies. Analyses highlighted the importance of different methods in treating and assessing hip and knee osteoarthritis. There was an inverse association between radiographic severity of osteoarthritis and pain and function in patients waiting for TKR but no association in THR. Different pain characteristics predicted long-term pain in THR and TKR. Outcomes after joint replacement should be assessed with a patient-reported outcome and a functional test. CONCLUSIONS: The RESTORE programme provides important information to guide the development of interventions to improve long-term outcomes for patients with osteoarthritis receiving THR and TKR. Issues relating to their evaluation and the assessment of patient outcomes are highlighted. Potential interventions at key times in the patient pathway were identified and deserve further study, ultimately in the context of a complex intervention.

A predictive computational fluid dynamics study of carotid artery

Abstract not available

Angiogenesis in the pathogenesis of inflammatory joint and lung diseases

This paper reviews hypotheses about roles of angiogenesis in the pathogenesis of inflammatory disease in two organs, the synovial joint and the lung. Neovascularisation is a fundamental process for growth and tissue repair after injury. Nevertheless, it may contribute to a variety of chronic inflammatory diseases, including rheumatoid arthritis, osteoarthritis, asthma, and pulmonary fibrosis. Inflammation can promote angiogenesis, and new vessels may enhance tissue inflammation. Angiogenesis in inflammatory disease may also contribute to tissue growth, disordered tissue perfusion, abnormal ossification, and enhanced responses to normal or pathological stimuli. Angiogenesis inhibitors may reduce inflammation and may also help to restore appropriate tissue structure and function

Quantification of variable palmar ligaments around the triquetrum-hamate joint determined by lunate type

The ligaments of the wrist are highly variable and poorly described, which is more obvious on the ulnar side of the wrist. Previous studies highlighted the potential differences within the ligaments of the wrist but no consensus has been reached. Poor tissue description and inconsistent use of terminology hindered the reproducibility of the results. Improved understanding of the morphological variations between carpal bones may facilitate improved understanding of the ligamentous structure within the wrist. This study aims to identify the potential variations between carpal bones that could be used to separate palmar ligamentous patterns around the triquetrum-hamate joint into subgroups within the sample population. Investigations were performed following a detailed nomenclature and a clear definition of ligamentous structures to facilitate detailed description and reproducible results. Quantitative analyses were conducted using 3D modelling technique. Histological sections were then analysed to identify the structure of each ligamentous attachment. Variable patterns of ligamentous attachments were identified. Differences were not only obvious between samples but also between the right and left hands of the same person. These identifications suggested that the palmar ligamentous patterns around the triquetrum-hamate joint are best described as a spectrum with a higher affinity of the triquetrum-hamate-capitate ligament and the lunate-triquetrum ligament to be associated with type I lunate wrists on one extreme and type II lunate wrists with the palmar triquetrum-hamate ligament, triquetrum-hamate-capitate ligament and palmar radius-lunate-triquetrum ligament attachments at the other extreme. Histological analyses confirmed pervious established work regarding the mechanical role of ligaments in wrist joint biomechanics. Also, there were no significant differences between the quantitative data obtained from the Genelyn-embalmed and unembalmed specimens (p>0.05). The current study demonstrated variable ligamentous patterns that suggest different bone restraints and two different patterns of motion. These findings support previous suggestions regarding separating the midcarpal joint into two distinct functional types. Type I wrists were identified with ligamentous attachments that are suggestive of rotating/translating hamate whilst type II wrists identified with ligamentous attachments that are suggestive of flexing/extending hamate motion based upon the patterns of the ligamentous attachments in relation to the morphological features of the underlying lunate type of the wrist. This opens the horizon for particular consideration and/or modification of surgical procedures, which may enhance the clinical management of wrist dysfunction.

Numerical simulation of strain-adaptive bone remodelling in the ankle joint

Background: The use of artificial endoprostheses has become a routine procedure for knee and hip joints while ankle arthritis has traditionally been treated by means of arthrodesis. Due to its advantages, the implantation of endoprostheses is constantly increasing. While finite element analyses (FEA) of strain-adaptive bone remodelling have been carried out for the hip joint in previous studies, to our knowledge there are no investigations that have considered remodelling processes of the ankle joint. In order to evaluate and optimise new generation implants of the ankle joint, as well as to gain additional knowledge regarding the biomechanics, strain-adaptive bone remodelling has been calculated separately for the tibia and the talus after providing them with an implant. Methods: FE models of the bone-implant assembly for both the tibia and the talus have been developed. Bone characteristics such as the density distribution have been applied corresponding to CT scans. A force of 5,200 N, which corresponds to the compression force during normal walking of a person with a weight of 100 kg according to Stauffer et al., has been used in the simulation. The bone adaptation law, previously developed by our research team, has been used for the calculation of the remodelling processes. Results: A total bone mass loss of 2% in the tibia and 13% in the talus was calculated. The greater decline of density in the talus is due to its smaller size compared to the relatively large implant dimensions causing remodelling processes in the whole bone tissue. In the tibia, bone remodelling processes are only calculated in areas adjacent to the implant. Thus, a smaller bone mass loss than in the talus can be expected. There is a high agreement between the simulation results in the distal tibia and the literature regarding. Conclusions: In this study, strain-adaptive bone remodelling processes are simulated using the FE method. The results contribute to a better understanding of the biomechanical behaviour of the ankle joint and hence are useful for the optimisation of the implant geometry in the future.

Gait analysis of the hind limb in Labrador Retrievers with and without cranial cruciate ligament disease

Cranial cruciate ligament (CCL) deficiency is the leading cause of lameness affecting the stifle joints of large breed dogs, especially Labrador Retrievers. Although CCL disease has been studied extensively, its exact pathogenesis and the primary cause leading to CCL rupture remain controversial. However, weakening secondary to repetitive microtrauma is currently believed to cause the majority of CCL instabilities diagnosed in dogs. Techniques of gait analysis have become the most productive tools to investigate normal and pathological gait in human and veterinary subjects. The inverse dynamics analysis approach models the limb as a series of connected linkages and integrates morphometric data to yield information about the net joint moment, patterns of muscle power and joint reaction forces. The results of these studies have greatly advanced our understanding of the pathogenesis of joint diseases in humans. A muscular imbalance between the hamstring and quadriceps muscles has been suggested as a cause for anterior cruciate ligament rupture in female athletes. Based on these findings, neuromuscular training programs leading to a relative risk reduction of up to 80% has been designed. In spite of the cost and morbidity associated with CCL disease and its management, very few studies have focused on the inverse dynamics gait analysis of this condition in dogs. The general goals of this research were (1) to further define gait mechanism in Labrador Retrievers with and without CCL-deficiency, (2) to identify individual dogs that are susceptible to CCL disease, and (3) to characterize their gait. The mass, location of the center of mass (COM), and mass moment of inertia of hind limb segments were calculated using a noninvasive method based on computerized tomography of normal and CCL-deficient Labrador Retrievers. Regression models were developed to determine predictive equations to estimate body segment parameters on the basis of simple morphometric measurements, providing a basis for nonterminal studies of inverse dynamics of the hind limbs in Labrador Retrievers. Kinematic, ground reaction forces (GRF) and morphometric data were combined in an inverse dynamics approach to compute hock, stifle and hip net moments, powers and joint reaction forces (JRF) while trotting in normal, CCL-deficient or sound contralateral limbs. Reductions in joint moment, power, and loads observed in CCL-deficient limbs were interpreted as modifications adopted to reduce or avoid painful mobilization of the injured stifle joint. Lameness resulting from CCL disease affected predominantly reaction forces during the braking phase and the extension during push-off. Kinetics also identified a greater joint moment and power of the contralateral limbs compared with normal, particularly of the stifle extensor muscles group, which may correlate with the lameness observed, but also with the predisposition of contralateral limbs to CCL deficiency in dogs. For the first time, surface EMG patterns of major hind limb muscles during trotting gait of healthy Labrador Retrievers were characterized and compared with kinetic and kinematic data of the stifle joint. The use of surface EMG highlighted the co-contraction patterns of the muscles around the stifle joint, which were documented during transition periods between flexion and extension of the joint, but also during the flexion observed in the weight bearing phase. Identification of possible differences in EMG activation characteristics between healthy patients and dogs with or predisposed to orthopedic and neurological disease may help understanding the neuromuscular abnormality and gait mechanics of such disorders in the future. Conformation parameters, obtained from femoral and tibial radiographs, hind limb CT images, and dual-energy X-ray absorptiometry, of hind limbs predisposed to CCL deficiency were compared with the conformation parameters from hind limbs at low risk. A combination of tibial plateau angle and femoral anteversion angle measured on radiographs was determined optimal for discriminating predisposed and non-predisposed limbs for CCL disease in Labrador Retrievers using a receiver operating characteristic curve analysis method. In the future, the tibial plateau angle (TPA) and femoral anteversion angle (FAA) may be used to screen dogs suspected of being susceptible to CCL disease. Last, kinematics and kinetics across the hock, stifle and hip joints in Labrador Retrievers presumed to be at low risk based on their radiographic TPA and FAA were compared to gait data from dogs presumed to be predisposed to CCL disease for overground and treadmill trotting gait. For overground trials, extensor moment at the hock and energy generated around the hock and stifle joints were increased in predisposed limbs compared to non predisposed limbs. For treadmill trials, dogs qualified as predisposed to CCL disease held their stifle at a greater degree of flexion, extended their hock less, and generated more energy around the stifle joints while trotting on a treadmill compared with dogs at low risk. This characterization of the gait mechanics of Labrador Retrievers at low risk or predisposed to CCL disease may help developing and monitoring preventive exercise programs to decrease gastrocnemius dominance and strengthened the hamstring muscle group.

Multi-body simulation of a canine hind limb: model development, experimental validation and calculation of ground reaction forces

Background: Among other causes the long-term result of hip prostheses in dogs is determined by aseptic loosening. A prevention of prosthesis complications can be achieved by an optimization of the tribological system which finally results in improved implant duration. In this context a computerized model for the calculation of hip joint loadings during different motions would be of benefit. In a first step in the development of such an inverse dynamic multi-body simulation (MBS-) model we here present the setup of a canine hind limb model applicable for the calculation of ground reaction forces. Methods: The anatomical geometries of the MBS-model have been established using computer tomography- (CT-) and magnetic resonance imaging- (MRI-) data. The CT-data were collected from the pelvis, femora, tibiae and pads of a mixed-breed adult dog. Geometric information about 22 muscles of the pelvic extremity of 4 mixed-breed adult dogs was determined using MRI. Kinematic and kinetic data obtained by motion analysis of a clinically healthy dog during a gait cycle (1 m/s) on an instrumented treadmill were used to drive the model in the multi-body simulation. Results and Discussion: As a result the vertical ground reaction forces (z-direction) calculated by the MBS-system show a maximum deviation of 1.75%BW for the left and 4.65%BW for the right hind limb from the treadmill measurements. The calculated peak ground reaction forces in z- and y-direction were found to be comparable to the treadmill measurements, whereas the curve characteristics of the forces in y-direction were not in complete alignment. Conclusion: In conclusion, it could be demonstrated that the developed MBS-model is suitable for simulating ground reaction forces of dogs during walking. In forthcoming investigations the model will be developed further for the calculation of forces and moments acting on the hip joint during different movements, which can be of help in context with the in silico development and testing of hip prostheses.

Relationship between internal derangement of temporomandibular joint and changes in body posture

Background: The presence of body posture changes among patients with temporomandibular disorders (TMD) has been a controversial issue in the literature, in which it supporters point out the muscular origin as the main etiological factors, mainly associated with postural changes in head. Due to this controversy, it is pertinent to check whether this relationship exists on the most common etiology of TMD, the disk displacement, which translates a biomechanical internal disorder of the temporomandibular joint (TMJ). Objectives: Assess body posture changes in subjects with internal derangement of the TMJ when compared to subjects without this biomechanical dysfunction, characterize the patterns of the jaw movements and assess to the muscle activation during jaw movements. Methods: 21 subjects with TMJ disc displacement (DD) (test group) and 21 subjects without any TMD (control group) was assessed for body posture changes through evaluation of several body segments by posturography and also was evaluated the postural balance reactions through the center of mass during jaw movements using a balance platform. For the characterization of the jaw movement patterns it was done a kinematic analysis during jaw movements (active ROM and path of the jaw). For the muscle activation during jaw movements it was evaluated the masseter, sternocleidomastoid and spinae erector muscles by surface electromyography (EMG). Results Discussion: Both groups show forward head posture and extension of the cervical spine, not noticing any other significant body posture changes in subjects with DD, and if we had to see in detail, in general, subjects without TMD shows more body posture changes than subjects with DD. The pattern of jaw movements is similar in both groups, but in subjects with DD the closing movements are more instable than the opening movements, related to a less effective movement control to counteract the force of gravity and the disk displacement. The bilateral muscle activation during jaw movements is higher in subjects with DD, likely related to a less stable pattern of movement which leads in a higher muscle activation to guide the movement and ensure the best as possible articular stability. Conclusion: The disk displacement with reduction should be viewed as part of a set of signs and symptoms that require an accurate musculoskeletal and psychosocial assessment towards an earlier diagnosis for reduction and control of the functional limiting factors. In this direction, it seems that the relevant set of limiting signs and symptoms deserve a particular attention by health care practitioners involved in the assessment and treatment of TMD, in order to define effective therapeutic options.

Elevated hepatocyte growth factor levels in osteoarthritis osteoblasts contribute to their altered response to bone morphogenetic protein-2 and reduced mineralization capacity

International audience

Upper Limb Posture Estimation in Robotic and Virtual Reality-based Rehabilitation.

New motor rehabilitation therapies include virtual reality (VR) and robotic technologies. In limb rehabilitation, limb posture is required to (1) provide a limb realistic representation in VR games and (2) assess the patient improvement. When exoskeleton devices are used in the therapy, the measurements of their joint angles cannot be directly used to represent the posture of the patient limb, since the human and exoskeleton kinematic models differ. In response to this shortcoming, we propose a method to estimate the posture of the human limb attached to the exoskeleton. We use the exoskeleton joint angles measurements and the constraints of the exoskeleton on the limb to estimate the human limb joints angles. This paper presents (a) the mathematical formulation and solution to the problem, (b) the implementation of the proposed solution on a commercial exoskeleton system for the upper limb rehabilitation, (c) its integration into a rehabilitation VR game platform, and (d) the quantitative assessment of the method during elbow and wrist analytic training. Results show that this method properly estimates the limb posture to (i) animate avatars that represent the patient in VR games and (ii) obtain kinematic data for the patient assessment during elbow and wrist analytic rehabilitation.

Biomechanical and morphological adaptations in women during pregnancy and the postpartum period in walking

A gravidez é uma fase especial da vida, com diversas alterações nos sistemas hormonais, anatómicos, e na composição corporal da mulher. No entanto, não é claro que alterações biomecânicas tridimensionais ocorrerem. Através do acompanhamento da mulher na gravidez e pós-parto, os objetivos da presente tese foram: 1) determinar os parâmetros temporais e espaciais do ciclo da marcha; 2) descrever a cinemática angular do membro inferior; 3) calcular os momentos e potências articulares do tornozelo, joelho e coxofemoral, utilizando o cálculo por dinâmica inversa; 4) descrever as magnitudes dos picos dos momentos e potências articulares dos membros inferiores; 5) identificar possíveis diferenças entre as fases de recolha relativamente aos parâmetros biomecânicos; 6) descrever longitudinalmente a composição corporal as alterações morfológicas; 7) analisar a influência das alterações antropométricas na cinética articular. Os resultados mostram que as mulheres mantêm os parâmetros temporais e espaciais da marcha. A cinemática angular do membro inferior tem o mesmo padrão, no entanto, a magnitude de alguns picos, especialmente na bacia e coxofemoral durante a fase terminal do apoio, pré-balanço e de balanço, apresentam alterações significativas. A coxofemoral é a articulação com mais alterações na cinética articular, com um aumento da carga interna associada aos momentos articulares da coxofemoral no plano transversal. No entanto, diversos momentos e potências articulares revelam uma diminuição significativa para o final da gravidez e/ou um aumento entre alguns trimestres da gravidez e o pós-parto. Como esperado, a maioria das variáveis associadas à composição corporal e às dimensões corporais tem um aumento significativo durante a gravidez e uma diminuição no pós-parto. Os modelos desenvolvidos para prever a carga interna aplicada ao membro inferior da grávida através de variáveis antropométricas, incluem quatro modelos com variáveis associadas à quantidade de gordura, quatro modelos com variáveis associadas à massa corporal global, três modelos que incluem a massa livre de gordura, e um modelo que inclui a forma do tronco. Os altos valores do R2 ajustado, mostram que as alterações na composição corporal e morfologia, determinam em grande parte a cinética articular da mulher nesta fase particular da vida.

Three-dimensional kinetic adaptations of gait throughout pregnancy and postpartum

Biomechanical adaptations that occur during pregnancy can lead to changes on gait pattern. Nevertheless, these adaptations of gait are still not fully understood. The purpose was to determine the effect of pregnancy on the biomechanical pattern of walking, regarding the kinetic parameters. A three-dimensional analysis was performed in eleven participants. The kinetic parameters in the joints of the lower limb during gait were compared at the end of the first, second, and third trimesters of pregnancy and in the postpartum period, in healthy pregnant women. The main results showed a reduction in the normalized vertical reaction forces, throughout pregnancy, particularly the third peak. Pregnant women showed, during most of the stance phase, medial reaction forces as a motor response to promote the body stability. Bilateral changes were observed in hip joint, with a decrease in the participation of the hip extensors and in the eccentric contraction of hip flexors. In ankle joint a decrease in the participation of ankle plantar flexors was found. In conclusion, the overall results point to biomechanical adjustments that showed a decrease of the mechanical load of women throughout pregnancy, with exception for few unilateral changes of hip joint moments.

Three-dimensional kinetic adaptations of gait throughout pregnancy and postpartum

Biomechanical adaptations that occur during pregnancy can lead to changes on gait pattern. Nevertheless, these adaptations of gait are still not fully understood. The purpose was to determine the effect of pregnancy on the biomechanical pattern of walking, regarding the kinetic parameters. A three-dimensional analysis was performed in eleven participants. The kinetic parameters in the joints of the lower limb during gait were compared at the end of the first, second, and third trimesters of pregnancy and in the postpartum period, in healthy pregnant women. The main results showed a reduction in the normalized vertical reaction forces, throughout pregnancy, particularly the third peak. Pregnant women showed, during most of the stance phase, medial reaction forces as a motor response to promote the body stability. Bilateral changes were observed in hip joint, with a decrease in the participation of the hip extensors and in the eccentric contraction of hip flexors. In ankle joint a decrease in the participation of ankle plantar flexors was found. In conclusion, the overall results point to biomechanical adjustments that showed a decrease of the mechanical load of women throughout pregnancy, with exception for few unilateral changes of hip joint moments.