Articular human joint modelling

Copyright @ Cambridge University Press 2009.The work reported in this paper encapsulates the theories and algorithms developed to drive the core analysis modules of the software which has been developed to model a musculoskeletal structure of anatomic joints. Due to local bone surface and contact geometry based joint kinematics, newly developed algorithms make the proposed modeller different from currently available modellers. There are many modellers that are capable of modelling gross human body motion. Nevertheless, none of the available modellers offer complete elements of joint modelling. It appears that joint modelling is an extension of their core analysis capability, which, in every case, appears to be musculoskeletal motion dynamics. It is felt that an analysis framework that is focused on human joints would have significant benefit and potential to be used in many orthopaedic applications. The local mobility of joints has a significant influence in human motion analysis, in understanding of joint loading, tissue behaviour and contact forces. However, in order to develop a bone surface based joint modeller, there are a number of major problems, from tissue idealizations to surface geometry discretization and non-linear motion analysis. This paper presents the following: (a) The physical deformation of biological tissues as linear or non-linear viscoelastic deformation, based on spring-dashpot elements. (b) The linear dynamic multibody modelling, where the linear formulation is established for small motions and is particularly useful for calculating the equilibrium position of the joint. This model can also be used for finding small motion behaviour or loading under static conditions. It also has the potential of quantifying the joint laxity. (c) The non-linear dynamic multibody modelling, where a non-matrix and algorithmic formulation is presented. The approach allows handling complex material and geometrical nonlinearity easily. (d) Shortest path algorithms for calculating soft tissue line of action geometries. The developed algorithms are based on calculating minimum ‘surface mass’ and ‘surface covariance’. An improved version of the ‘surface covariance’ algorithm is described as ‘residual covariance’. The resulting path is used to establish the direction of forces and moments acting on joints. This information is needed for linear or non-linear treatment of the joint motion. (e) The final contribution of the paper is the treatment of the collision. In the virtual world, the difficulty in analysing bodies in motion arises due to body interpenetrations. The collision algorithm proposed in the paper involves finding the shortest projected ray from one body to the other. The projection of the body is determined by the resultant forces acting on it due to soft tissue connections under tension. This enables the calculation of collision condition of non-convex objects accurately. After the initial collision detection, the analysis involves attaching special springs (stiffness only normal to the surfaces) at the ‘potentially colliding points’ and motion of bodies is recalculated. The collision algorithm incorporates the rotation as well as translation. The algorithm continues until the joint equilibrium is achieved. Finally, the results obtained based on the software are compared with experimental results obtained using cadaveric joints

Soft tissue structure modelling for use in orthopaedic applications and musculoskeletal biomechanics

We present our methodology for the three-dimensional anatomical and geometrical description of soft tissues, relevant for orthopaedic surgical applications and musculoskeletal biomechanics. The technique involves the segmentation and geometrical description of muscles and neurovascular structures from high-resolution computer tomography scanning for the reconstruction of generic anatomical models. These models can be used for quantitative interpretation of anatomical and biomechanical aspects of different soft tissue structures. This approach should allow the use of these data in other application fields, such as musculoskeletal modelling, simulations for radiation therapy, and databases for use in minimally invasive, navigated and robotic surgery

Muscle Synergies Facilitate Computational Prediction of Subject-Specific Walking Motions.

Researchers have explored a variety of neurorehabilitation approaches to restore normal walking function following a stroke. However, there is currently no objective means for prescribing and implementing treatments that are likely to maximize recovery of walking function for any particular patient. As a first step toward optimizing neurorehabilitation effectiveness, this study develops and evaluates a patient-specific synergy-controlled neuromusculoskeletal simulation framework that can predict walking motions for an individual post-stroke. The main question we addressed was whether driving a subject-specific neuromusculoskeletal model with muscle synergy controls (5 per leg) facilitates generation of accurate walking predictions compared to a model driven by muscle activation controls (35 per leg) or joint torque controls (5 per leg). To explore this question, we developed a subject-specific neuromusculoskeletal model of a single high-functioning hemiparetic subject using instrumented treadmill walking data collected at the subject's self-selected speed of 0.5 m/s. The model included subject-specific representations of lower-body kinematic structure, foot-ground contact behavior, electromyography-driven muscle force generation, and neural control limitations and remaining capabilities. Using direct collocation optimal control and the subject-specific model, we evaluated the ability of the three control approaches to predict the subject's walking kinematics and kinetics at two speeds (0.5 and 0.8 m/s) for which experimental data were available from the subject. We also evaluated whether synergy controls could predict a physically realistic gait period at one speed (1.1 m/s) for which no experimental data were available. All three control approaches predicted the subject's walking kinematics and kinetics (including ground reaction forces) well for the model calibration speed of 0.5 m/s. However, only activation and synergy controls could predict the subject's walking kinematics and kinetics well for the faster non-calibration speed of 0.8 m/s, with synergy controls predicting the new gait period the most accurately. When used to predict how the subject would walk at 1.1 m/s, synergy controls predicted a gait period close to that estimated from the linear relationship between gait speed and stride length. These findings suggest that our neuromusculoskeletal simulation framework may be able to bridge the gap between patient-specific muscle synergy information and resulting functional capabilities and limitations

The effect of boundary constraints on finite element modelling of the human pelvis

The use of finite element analysis (FEA) to investigate the biomechanics of anatomical systems critically relies on the specification of physiologically representative boundary conditions. The biomechanics of the pelvis has been the specific focus of a number of FEA studies previously, but it is also a key aspect in other investigations of, for example, the hip joint or new design of hip prostheses. In those studies, the pelvis has been modelled in a number of ways with a variety of boundary conditions, ranging from a model of the whole pelvic girdle including soft tissue attachments to a model of an isolated hemi-pelvis. The current study constructed a series of FEA models of the same human pelvis to investigate the sensitivity of the predicted stress distributions to the type of boundary conditions applied, in particular to represent the sacro-iliac joint and pubic symphysis. Varying the method of modelling the sacro-iliac joint did not produce significant variations in the stress distribution, however changes to the modelling of the pubic symphysis were observed to have a greater effect on the results. Over-constraint of the symphysis prevented the bending of the pelvis about the greater sciatic notch, and underestimated high stresses within the ilium. However, permitting medio-lateral translation to mimic widening of the pelvis addressed this problem. These findings underline the importance of applying the appropriate boundary conditions to FEA models, and provide guidance on suitable methods of constraining the pelvis when, for example, scan data has not captured the full pelvic girdle. The results also suggest a valid method for performing hemi-pelvic modelling of cadaveric or archaeological remains which are either damaged or incomplete

Biomechanics

Biomechanics is a vast discipline within the field of Biomedical Engineering. It explores the underlying mechanics of how biological and physiological systems move. It encompasses important clinical applications to address questions related to medicine using engineering mechanics principles. Biomechanics includes interdisciplinary concepts from engineers, physicians, therapists, biologists, physicists, and mathematicians. Through their collaborative efforts, biomechanics research is ever changing and expanding, explaining new mechanisms and principles for dynamic human systems. Biomechanics is used to describe how the human body moves, walks, and breathes, in addition to how it responds to injury and rehabilitation. Advanced biomechanical modeling methods, such as inverse dynamics, finite element analysis, and musculoskeletal modeling are used to simulate and investigate human situations in regard to movement and injury. Biomechanical technologies are progressing to answer contemporary medical questions. The future of biomechanics is dependent on interdisciplinary research efforts and the education of tomorrow’s scientists

Cancellous bone and theropod dinosaur locomotion. Part II—a new approach to inferring posture and locomotor biomechanics in extinct tetrapod vertebrates

This paper is the second of a three-part series that investigates the architecture of cancellous bone in the main hindlimb bones of theropod dinosaurs, and uses cancellous bone architectural patterns to infer locomotor biomechanics in extinct non-avian species. Cancellous bone is widely known to be highly sensitive to its mechanical environment, and therefore has the potential to provide insight into locomotor biomechanics in extinct tetrapod vertebrates such as dinosaurs. Here in Part II, a new biomechanical modelling approach is outlined, one which mechanistically links cancellous bone architectural patterns with three-dimensional musculoskeletal and finite element modelling of the hindlimb. In particular, the architecture of cancellous bone is used to derive a single ‘characteristic posture’ for a given species—one in which bone continuum-level principal stresses best align with cancellous bone fabric—and thereby clarify hindlimb locomotor biomechanics. The quasi-static approach was validated for an extant theropod, the chicken, and is shown to provide a good estimate of limb posture at around mid-stance. It also provides reasonable predictions of bone loading mechanics, especially for the proximal hindlimb, and also provides a broadly accurate assessment of muscle recruitment insofar as limb stabilization is concerned. In addition to being useful for better understanding locomotor biomechanics in extant species, the approach hence provides a new avenue by which to analyse, test and refine palaeobiomechanical hypotheses, not just for extinct theropods, but potentially many other extinct tetrapod groups as well

Adaptive Surrogate Modeling for Efficient Coupling of Musculoskeletal Control and Tissue Deformation Models

Background Finite element (FE) modeling and multibody dynamics have traditionally been applied separately to the domains of tissue mechanics and musculoskeletal movements, respectively. Simultaneous simulation of both domains is needed when interactions between tissue and movement are of interest, but this has remained largely impractical due to high computational cost. Method of Approach Here we present a method for concurrent simulation of tissue and movement, in which state of the art methods are used in each domain, and communication occurs via a surrogate modeling system based on locally weighted regression. The surrogate model only performs FE simulations when regression from previous results is not within a user-specified tolerance. For proof of concept and to illustrate feasibility, the methods were demonstrated on an optimization of jumping movement using a planar musculoskeletal model coupled to a FE model of the foot. To test the relative accuracy of the surrogate model outputs against those of the FE model, a single forward dynamics simulation was performed with FE calls at every integration step and compared with a corresponding simulation with the surrogate model included. Neural excitations obtained from the jump height optimization were used for this purpose and root mean square (RMS) difference between surrogate and FE model outputs (ankle force and moment, peak contact pressure and peak von Mises stress) were calculated. Results Optimization of jump height required 1800 iterations of the movement simulation, each requiring thousands of time steps. The surrogate modeling system only used the FE model in 5% of time steps, i.e. a 95% reduction of computation time. Errors introduced by the surrogate model were less than 1 mm in jump height and RMS errors of less than 2 N in ground reaction force, 0.25 Nm in ankle moment, and 10 kPa in peak tissue stress. Conclusion Adaptive surrogate modeling based on local regression allows efficient concurrent simulations of tissue mechanics and musculoskeletal movement

A novel musculoskeletal joint modelling for orthopaedic applications

This thesis was submitted for the degree of Docter of Philosophy and awarded by Brunel University.The objective of the work carried out in this thesis was to develop analytical and computational tools to model and investigate musculoskeletal human joints. It was recognised that the FEA was used by many researchers in modelling human musculoskeletal motion, loading and stresses. However the continuum mechanics played only a minor role in determining the articular joint motion, and its value was questionable. This is firstly due to the computational cost and secondly due to its impracticality for this application. On the other hand, there isn’t any suitable software for precise articular joint motion analysis to deal with the local joint stresses or non standard joints. The main requirement in orthopaedics field is to develop a modeller software (and its associated theories) to model anatomic joint as it is, without any simplification with respect to joint surface morphology and material properties of surrounding tissues. So that the proposed modeller can be used for evaluating and diagnosing different joint abnormalities but furthermore form the basis for performing implant insertion and analysis of the artificial joints. The work which is presented in this thesis is a new frame work and has been developed for human anatomic joint analysis which describes the joint in terms of its surface geometry and surrounding musculoskeletal tissues. In achieving such a framework several contributions were made to the 6DOF linear and nonlinear joint modelling, the mathematical definition of joint stiffness, tissue path finding and wrapping and the contact with collision analysis. In 6DOF linear joint modelling, the contribution is the development of joint stiffness and damping matrices. This modelling approach is suitable for the linear range of tissue stiffness and damping properties. This is the first of its kind and it gives a firm analytical basis for investigating joints with surrounding tissue and the cartilage. The 6DOF nonlinear joint modelling is a new scheme which is described for modelling the motion of multi bodies joined by non-linear stiffness and contact elements. The proposed method requires no matrix assembly for the stiffness and damping elements or mass elements. The novelty in the nonlinear modelling, relates to the overall algorithmic approach and handling local non-linearity by procedural means. The mathematical definition of joint stiffness is also a new proposal which is based on the mathematical definition of stiffness between two bodies. Based on the joint stiffness matrix properties, number of joint stiffness invariants was obtained analytically such as the centre of stiffness, the principal translational stiffnesses, and the principal rotational stiffnesses. In corresponding to these principal stiffnesses, their principal axes have been also obtained. Altogether, a joint is assessed by six principal axes and six principal stiffnesses and its centre of stiffness. These formulations are new and show that a joint can be described in terms of inherent stiffness properties. It is expected that these will be better in characterising a joint in comparison to laxity based characterisation. The development of tissue path finding and wrapping algorithms are also introduced as new approaches. The musculoskeletal tissue wrapping involves calculating the shortest distance between two points on a meshed surface. A new heuristic algorithm was proposed. The heuristic is based on minimising the accumulative divergence from the straight line between two points on the surface and the direction of travel on the surface (i.e. bone). In contact and collision based development, the novel algorithm has been proposed that detects possible colliding points on the motion trajectory by redefining the distance as a two dimensional measure along the velocity approach vector and perpendicular to this vector. The perpendicular distance determines if there are potentially colliding points, and the distance along the velocity determines how close they are. The closest pair among the potentially colliding points gives the “time to collision”. The algorithm can eliminate the “fly pass” situation where very close points may not collide because of the direction of their relative velocity. All these developed algorithms and modelling theories, have been encompassed in the developed prototype software in order to simulate the anatomic joint articulations through modelling formulations developed. The software platform provides a capability for analysing joints as 6DOF joints based on anatomic joint surfaces. The software is highly interactive and driven by well structured database, designed to be highly flexible for the future developments. Particularly, two case studies are carried out in this thesis in order to generate results relating to all the proposed elements of the study. The results obtained from the case studies show good agreement with previously published results or model based results obtained from Lifemod software, whenever comparison was possible. In some cases the comparison was not possible because there were no equivalent results; the results were supported by other indicators. The modelling based results were also supported by experiments performed in the Brunel Orthopaedic Research and Learning Centre