Andrew T. M. Phillips

An open source lower limb model: Hip joint validation

Musculoskeletal lower limb models have been shown to be able to predict hip contact forces (HCFs) that are comparable to in vivo measurements obtained from instrumented prostheses. However, the muscle recruitment predicted by these models does not necessarily compare well to measured electromyographic (EMG) signals. In order to verify if it is possible to accurately estimate HCFs from muscle force patterns consistent with EMG measurements, a lower limb model based on a published anatomical dataset (Klein Horsman et al., 2007. Clinical Biomechanics. 22, 239-247) has been implemented in the open source software OpenSim. A cycle-to-cycle hip joint validation was conducted against HCFs recorded during gait and stair climbing trials of four arthroplasty patients (Bergmann et al., 2001. Journal of Biomechanics. 34, 859-871). Hip joint muscle tensions were estimated by minimizing a polynomial function of the muscle forces. The resulting muscle activation patterns obtained by assessing multiple powers of the objective function were compared against EMG profiles from the literature. Calculated HCFs denoted a tendency to monotonically increase their magnitude when raising the power of the objective function; the best estimation obtained from muscle forces consistent with experimental EMG profiles was found when a quadratic objective function was minimized (average overestimation at experimental peak frame: 10.1% for walking, 7.8% for stair climbing). The lower limb model can produce appropriate balanced sets of muscle forces and joint contact forces that can be used in a range of applications requiring accurate quantification of both. The developed model is available at the website https://simtk.org/home/low_limb_london.

A quadratic assignment formulation of the molecular conformation problem

The molecular conformation problem is discussed, and a concave quadratic global minimization approach for solving it is described. This approach is based on a quadratic assignment formulation of a discrete approximation to the original problem.

Global optimization of fractional programs

Dinkelbachs global optimization approach for finding the global maximum of the fractional programming problem is discussed. Based on this idea, a modified algorithm is presented which provides both upper and lower bounds at each iteration. The convergence of the lower and upper bounds to the global maximum function value is shown to be superlinear. In addition, the special case of fractional programming when the ratio involves only linear or quadratic terms is considered. In this case, the algorithm is guaranteed to find the global maximum to within any specified tolerance, regardless of the definiteness of the quadratic form.

The femur as a musculo-skeletal construct: a free boundary condition modelling approach.

Previous finite element studies of the femur have made simplifications to varying extents with regard to the boundary conditions used during analysis. Fixed boundary conditions are generally applied to the distal femur when examining the proximal behaviour at the hip joint, while the same can be said for the proximal femur when examining the distal behaviour at the knee joint. While fixed boundary condition analyses have been validated against in vitro experiments it remains a matter of debate as to whether the numerical and experimental models are indicative of the in vivo situation. This study presents a finite element model in which the femur is treated as a complete musculo-skeletal construct, spanning between the hip and knee joints. Linear and non-linear implementations of a free boundary condition modelling approach are applied to the bone through the explicit inclusion of muscles and ligaments spanning both the hip joint and the knee joint. A non-linear force regulated, muscle strain based activation strategy was found to result in lower observed principal strains in the cortex of the femur, compared to a linear activation strategy. The non-linear implementation of the model in particular, was found to produce hip and knee joint reaction forces consistent with in vivo data from instrumented implants.

Hip abduction can prevent posterior edge loading of hip replacements

Edge loading causes clinical problems for hard‐on‐hard hip replacements, and edge loading wear scars are present on the majority of retrieved components. We asked the question: are the lines of action of hip joint muscles such that edge loading can occur in a well‐designed, well‐positioned acetabular cup? A musculoskeletal model, based on cadaveric lower limb geometry, was used to calculate for each muscle, in every position within the complete range of motion, whether its contraction would safely pull the femoral head into the cup or contribute to edge loading. The results show that all the muscles that insert into the distal femur, patella, or tibia could cause edge loading of a well‐positioned cup when the hip is in deep flexion. Patients frequently use distally inserting muscles for movements requiring deep hip flexion, such as sit‐to‐stand. Importantly, the results, which are supported by in vivo data and clinical findings, also show that risk of edge loading is dramatically reduced by combining deep hip flexion with hip abduction. Patients, including those with sub‐optimally positioned cups, may be able to reduce the prevalence of edge loading by rising from chairs or stooping with the hip abducted.

Application of a falsification strategy to a musculoskeletal model of the lower limb and accuracy of the predicted hip contact force vector

In the literature, lower limb musculoskeletal models validated against in vivo measured hip contact forces (HCFs) exhibit a tendency to overestimate the HCFs magnitude and predict inaccurate components of the HCF vector in the transverse plane. In order to investigate this issue, a musculoskeletal model was forced to produce HCFs identical to those measured and the resulting joint equilibrium equations were studied through both a general approach and a static optimization framework. In the former case, the existence of solutions to the equilibrium equations was investigated and the effect of varying the intersegmental moments and the muscle tetanic stress assessed: for a value of 100 N/cm(2) and moments calculated from an inverse dynamics analysis on average only 62% of analyzed frames were solvable for level walking and 70% for stair climbing. In the static optimization study, the model could reproduce the experimental HCFs but the recruited muscles were unable to simultaneously equilibrate the hip intersegmental moments without the contribution of reserve moment actuators. Without constraints imposed on the HCFs, the predicted HCF vectors presented maximum angle deviations up to 22° for level walking and 33° for stair climbing during the gait stance phase. The influence of the medio-lateral HCF component on the solvability of the equilibrium equations and the muscle recruitment alteration when the model was forced to produce the experimental HCFs suggest that a more accurate geometrical representation of the gluteal muscles is mandatory to improve predictions of the HCF vector yielded by the static optimization technique.

CGU: An Algorithm for Molecular Structure Prediction

A global optimization method is presented for predicting the minimum energy structure of small protein-like molecules. This method begins by collecting a large number of molecular conformations, each obtained by finding a local minimum of a potential energy function from a random starting point. The information from these conformera is then used to form a convex quadratic global underestimating function for the potential energy of all known conformers. This underestimator is an L1 approximation to all known local minima, and is obtained by a linear programming formulation and solution. The minimum of this underestimator is used to predict the global minimum for the function, allowing a localized conformer search to be performed based on the predicted minimum. The new set of conformers generated by the localized search serves as the basis for another quadratic underestimation step in an iterative algorithm. This algorithm has been used to predict the minimum energy structures of heteropolymers with as many as 48 residues, and can be applied to a variety of molecular models. The results obtained also show the dependence of the native conformation on the sequence of hydrophobic and polar residues.

Modeling the biomechanics of fetal movements

Fetal movements in the uterus are a natural part of development and are known to play an important role in normal musculoskeletal development. However, very little is known about the biomechanical stimuli that arise during movements in utero, despite these stimuli being crucial to normal bone and joint formation. Therefore, the objective of this study was to create a series of computational steps by which the forces generated during a kick in utero could be predicted from clinically observed fetal movements using novel cine-MRI data of three fetuses, aged 20–22 weeks. A custom tracking software was designed to characterize the movements of joints in utero, and average uterus deflection of

Evaluating a suitable level of model complexity for finite element analysis of the intact acetabulum

6.95 \pm 0.41