Anatomically based finite element geometries are becoming increasingly popular in physiological modelling, owing to the demand for modelling that links organ function to spatially distributed properties at the protein, cell and tissue level. We present a collection of anatomically based finite element geometries of the musculo-skeletal system and other organs suitable for use in continuum analysis. These meshes are derived from the widely used Visible Human (VH) dataset and constitute a contribution to the world wide International Union of Physiological Sciences (IUPS) Physiome Project (www.physiome.org.nz). The method of mesh generation and fitting of tricubic Hermite volume meshes to a given dataset is illustrated using a least-squares algorithm that is modified with smoothing (Sobolev) constraints via the penalty method to account for sparse and scattered data. A technique ("host mesh" fitting) based on "free-form" deformation (FFD) is used to customise the fitted (generic) geometry. Lung lobes, the rectus femoris muscle and the lower limb bones are used as examples to illustrate these methods. Geometries of the lower limb, knee joint, forearm and neck are also presented. Finally, the issues and limitations of the methods are discussed.

No abstract

Musculoskeletal modeling and optimization theory are often used to determine muscle forces in vivo. However, convincing quantitative evaluation of these predictions has been limited to date. The present study evaluated model predictions of knee muscle forces during walking using in vivo measurements of joint contact loading acquired from an instrumented implant. Joint motion, ground reaction force, and tibial contact force data were recorded simultaneously from a single subject walking at slow, normal, and fast speeds. The body was modeled as an 8-segment, 21-degree-of-freedom articulated linkage, actuated by 58 muscles. Joint moments obtained from inverse dynamics were decomposed into leg-muscle forces by solving an optimization problem that minimized the sum of the squares of the muscle activations. The predicted knee muscle forces were input into a 3D knee implant contact model to calculate tibial contact forces. Calculated and measured tibial contact forces were in good agreement for all three walking speeds. The average RMS errors for the medial, lateral, and total contact forces over the entire gait cycle and across all trials were 140 AE 40 N, 115 AE 32 N, and 183 AE 45 N, respectively. Muscle coordination predicted by the model was also consistent with EMG measurements reported for normal walking. The combined experimental and modeling approach used in this study provides a quantitative framework for evaluating model predictions of muscle forces in human movement. ß

This study presents an investigation of the changes in foot posture, joint kinematics, joint moments and joint contact forces in the lower extremity following a 5 k treadmill run. A relationship between knee and ankle joint loading and foot posture index (FPI) is developed. Twenty recreational male heel-strike runners participated in this study. All participants had a history of running exercise and were free from lower extremity injuries and foot deformities. Foot posture was assessed from a six-item FPI to quantitatively classify high supination to high pronation foot poses. The FPI is scored using a combination of observations and foot palpations. The three-dimensional marker trajectories, ground reaction force and surface electromyography (EMG) were recorded at pre and post-gait sessions conducted over-ground and 5 k running was conducted on a treadmill. Joint kinematics, joint moments and joint contact forces were computed in OpenSim. Simulated EMG activations were compared against experimental EMG to validate the model. A paired sample t -test was conducted using a 1D statistical parametric mapping method computed temporally. Hip joint moments and contact forces increased during initial foot contact following 5 k running. Knee abduction moment and superior-inferior knee contact force increased, whereas the knee extension moment decreased. Ankle plantarflexion moment and ankle contact forces increased during stance. FPI was found to be moderately correlated with peak knee and ankle moments. Recreational male runners presented increased static foot pronation after 5 k treadmill running. These findings suggest that following mid distance running foot pronation may be an early indicator of increased lower limb joint loading. Furthermore, the FPI may be used to quantify the changes in knee and ankle joint moments.

A 3D anatomically based patient-specific finite element (FE) model of patello-femoral (PF) articulation is presented to analyse the main features of patella biomechanics, namely, patella tracking (kinematics), quadriceps extensor forces, surface contact and internal patella stresses. The generic geometries are a subset from the model database of the International Union of Physiological Sciences (IUPS) (http://www.physiome.org.nz) Physiome Project with soft tissue derived from the widely used visible human dataset, and the bones digitised from an anatomically accurate physical model with muscle attachment information. The models are customised to patient magnetic resonance images using a variant of free-form deformation, called 'host-mesh' fitting. The continuum was solved using the governing equation of finite elasticity, with the multibody problem coupled through contact mechanics. Additional constraints such as tissue incompressibility are also imposed. Passive material properties are taken from the literature and implemented for deformable tissue with a non-linear micro-structurally based constitutive law. Bone and cartilage are implemented using a 'St-Venant Kirchoff' model suitable for rigid body rotations. The surface fibre directions have been estimated from anatomy images of cadaver muscle dissections and active muscle contraction was based on a steady-state calcium-tension relation. The 3D continuum model of muscle, tendon and bone is compared with experimental results from the literature, and surgical simulations performed to illustrate its clinical assessment capabilities (a Maquet procedure for reducing patella stresses and a vastus lateralis release for a bipartite patella). Finally, the model limitations, issues and future improvements are discussed.

Physique changes during pregnancy lead to gait characteristic variations. This study aimed to analyse gait of pregnant individuals throughout pregnancy and post-partum. Sixteen healthy pregnant women volunteered as participants and had their lower limb kinematics analysed through a VICON three-dimensional motion system and plantar pressure measured with a Novel EMED force plate. Significant changes were observed in pelvic anterior motion, hip and ankle joint kinematics. Mean pressure distribution and COP trajectory deviation altered accordingly with increased pregnancy time, compared with post-partum. This longitudinal study of pregnant gait biomechanics in T2, T3 and PP reveals lower extremity kinematic and foot pressure alterations to adapt to pregnancy related changes, and the COP trajectory highlights a falling risk during pregnancy, particularly in T3.

Foot morphology and function has received increasing attention from both biomechanics researchers and footwear manufacturers. In this study, 168 habitually unshod runners (90 males whose age, weight & height were 23±2.4years, 66±7.1kg & 1.68±0.13m and 78 females whose age, weight & height were 22±1.8years, 55±4.7kg & 1.6±0.11m) (Indians) and 196 shod runners (130 males whose age, weight & height were 24±2.6years, 66±8.2kg & 1.72±0.18m and 66 females whose age, weight & height were 23±1.5years, 54±5.6kg & 1.62±0.15m)(Chinese) participated in a foot scanning test using the easy-foot-scan (a three-dimensional foot scanning system) to obtain 3D foot surface data and 2D footprint imaging. Foot length, foot width, hallux angle and minimal distance from hallux to second toe were calculated to analyze foot morphological differences. This study found that significant differences exist between groups (shod Chinese and unshod Indians) for foot length (female p = 0.001), width (female p = 0.001), hallux angle (male and female p = 0.001) and the minimal distance (male and female p = 0.001) from hallux to second toe. This study suggests that significant differences in morphology between different ethnicities could be considered for future investigation of locomotion biomechanics characteristics between ethnicities and inform last shape and design so as to reduce injury risks and poor performance from mal-fit shoes.

This study presents the kinematics and plantar pressure characteristics of eight elite national-level badminton athletes and eight recreational college-level badminton players while performing a right-forward lunge movement in a laboratory-simulated badminton court. The hypothesis was that recreational players would be significantly different from elite players in kinematics and plantar pressure measures. Vicon motion capture and Novel insole plantar pressure measurement were simultaneously taken to record the lower extremity kinematics and foot loading during stance. Recreational players showed significantly higher peak pressure in the lateral forefoot (P = 0.002) and force time integral in the lateral forefoot (P = 0.013) and other toes (P = 0.005). Elite athletes showed higher peak pressure in the medial forefoot (P = 0.003), hallux (P = 0.037) and force time integral in the medial forefoot (P = 0.009). The difference in landing techniques for the lunge step between elite athletes and recreational players was observed with peak ankle eversion (-38.2°±2.4° for athletes and -11.1°±3.9° for players, P = 0.015); smaller knee range of motion in the coronal and transverse planes, with differences in peak knee adduction (28.9°±6.8° for athletes and 15.7°±6.2° for players, P = 0.031); peak knee internal rotation (20.3°±1.3° for athletes and 11.8°±3.2° for players, P = 0.029) and peak hip flexion (77.3°±4.1° for athletes and 91.3°±9.3° for players, P = 0.037).