Rotational moment arms of the medial hamstrings and adductors vary with femoral geometry and limb position: implications for the treatment of internally rotated gait

Persons with cerebral palsy frequently walk with a crouched, internally rotated gait. Spastic medial hamstrings or adductors are presumed to contribute to excessive hip internal rotation in some patients; however, the capacity of these muscles to produce internal rotation has not been adequately investigated. The purpose of this study was to determine the hip rotation moment arms of the medial hamstrings and adductors in persons with femoral anteversion deformities who walk with a crouched, internally rotated gait. A musculoskeletal model with a "deformable" femur was developed. This model was used, in conjunction with kinematic data obtained from gait analysis, to calculate the muscle moment arms for combinations of joint angles and anteversion deformities exhibited by 21 subjects with cerebral palsy and excessive hip internal rotation. We found that the semimembranosus, semitendinosus, and gracilis muscles in our model had negligible or external rotation moment arms when the hip was internally rotated or the knee was flexed -- the body positions assumed by the subjects during walking. When the femur was excessively anteverted, the rotational moment arms of the adductor brevis, adductor longus, pectineus, and proximal compartments of the adductor magnus in our model shifted toward external rotation. These results suggest that neither the medial hamstrings nor the adductors are likely to contribute substantially to excessive internal rotation of the hip and that other causes of internal rotation should be considered when planning treatments for these patients.     

Influence of altered gait patterns on the hip joint contact forces

Children who exhibit gait deviations often present a range of bone deformities, particularly at the proximal femur. Altered gait may affect bone growth and lead to deformities by exerting abnormal stresses on the developing bones. The objective of this study was to calculate variations in the hip joint contact forces with different gait patterns. Muscle and hip joint contact forces of four children with different walking characteristics were calculated using an inverse dynamic analysis and a static optimisation algorithm. Kinematic and kinetic analyses were based on a generic musculoskeletal model scaled down to accommodate the dimensions of each child. Results showed that for all the children with altered gaits both the orientation and magnitude of the hip joint contact force deviated from normal. The child with the most severe gait deviations had hip joint contact forces 30% greater than normal, most likely due to the increase in muscle forces required to sustain his crouched stance. Determining how altered gait affects joint loading may help in planning treatment strategies to preserve correct loading on the bone from a young age.     

Knee joint kinematics during walking influences the spatial cartilage thickness distribution in the knee

The regional adaptation of knee cartilage morphology to the kinematics of walking has been suggested as an important factor in the evaluation of the consequences of alteration in normal gait leading to osteoarthritis. The purpose of this study was to investigate the association of spatial cartilage thickness distributions of the femur and tibia in the knee to the knee kinematics during walking. Gait data and knee MR images were obtained from 17 healthy volunteers (age 33.2 ± 9.8 years). Cartilage thickness maps were created for the femoral and tibial cartilage. Locations of thickest cartilage in the medial and lateral compartments in the femur and tibia were identified using a numerical method. The flexion-extension (FE) angle associated with the cartilage contact regions on the femur, and the anterior-posterior (AP) translation and internal-external (IE) rotation associated with the cartilage contact regions on the tibia at the heel strike of walking were tested for correlation with the locations of thickest cartilage. The locations of the thickest cartilage had relatively large variation (SD, 8.9°) and was significantly associated with the FE angle at heel strike only in the medial femoral condyle (R(2)=0.41, p<0.01). The natural knee kinematics and contact surface shapes seem to affect the functional adaptation of knee articular cartilage morphology. The sensitivity of cartilage morphology to kinematics at the knee during walking suggests that regional cartilage thickness variations are influenced by both loading and the number of loading cycles. Thus walking is an important consideration in the analysis of the morphological variations of articular cartilage, since it is the dominant cyclic activity of daily living. The sensitivity of cartilage morphology to gait kinematics is also important in understanding the etiology and pathomechanics of osteoarthritis.     

Femoral neck anteversion: values, development, measurement, common problems

The femoral neck anteversion angle is an important factor for hip stability and normal walking. It is multifactoral result of evolution, heredity, fetal development, intrauterine position, and mechanical forces. Abnormal FNA sometimes can be associated with many clinical problems ranging from harmless intoeing gait in the early childhood, to disabling osteoarthritis of the hip and the knee in the adults. In most cases is associated with minor functional problems in children during growth, but cause a concern in parents for children future. The child must be examined carefully and an accurate diagnosis must be established. The most important part of care is observation of the children. If abnormal femoral neck anteversion produces severe functional disability, derotational osteotomy should be done, but delayed until late childhood.     

Non-rigid deformation to include subject-specific detail in musculoskeletal models of CP children with proximal femoral deformity and its effect on muscle and contact forces during gait

To account for proximal femoral deformities in children with cerebral palsy (CP), subject-specific musculoskeletal models are needed. Non-rigid deformation (NRD) deforms generic onto personalized bone geometry and thereby transforms the muscle points. The goal of this study was to determine to what extent the models and simulation outcomes in CP patients differ when including subject-specific detail using NRD or Magnetic Resonance Imaging (MRI)-based models. The NRD models slightly overestimated hip contact forces compared to MRI models and differences in muscle point positions and moment arm lengths (MALs) remained, although differences were smaller than for the generic model.     

Development of a surrogate model based on patient weight,bone mass and geometry to predict femoral neck strains and fracture loads

Osteoporosis and related bone fractures are an increasing global burden in our ageing society. Areal bone mineral density assessed through dual energy X-ray absorptiometry (DEXA), the clinically accepted and most used method, is not sufficient to assess fracture risk individually. Finite element (FE) modelling has shown improvements in prediction of fracture risk, better than aBMD from DEXA, but is not practical for widespread clinical use. The aim of this study was to develop an adaptive neural network (ANN)-based surrogate model to predict femoral neck strains and fracture loads obtained from a previously developed population-based FE model. The surrogate model performance was assessed in simulating two loading conditions: the stance phase of gait and a fall.The surrogate model successfully predicted strains estimated by FE (r² = 0.90–0.98 for level gait load case, r² = 0.92–0.96 for the fall load case). Moreover, an ANN model based on three measurements obtainable in clinics (femoral neck length (level gait) or maximum femoral neck diameter (fall), femoral neck bone mass, body weight) was able to give reasonable predictions (r² = 0.84–0.94) for all of the strain metrics and the estimated femoral neck fracture load. Overall, the surrogate model has potential for clinical applications as they are based on simple measures of geometry and bone mass which can be derived from DEXA images, accurately predicting FE model outcomes, with advantages over FE models as they are quicker and easier to perform.     

Prediction of femoral fracture load using finite element models: an examination of stress- and strain-based failure theories

Finite element (FE) models are often used to model bone failure. However, no failure theory for bone has been validated at this time. In this study, we examined the performance of nine stress- and strain-based failure theories, six of which could account for differences in tensile and compressive material strengths. The distortion energy, Hoffman and a strain-based Hoffman analog, maximum normal stress, maximum normal strain, maximum shear strain, maximum shear stress (tau(max)), Coulomb-Mohr, and modified Mohr failure theories were evaluated using automatically generated, computed tomographic scan-based FE models of the femur. Eighteen matched pairs of proximal femora were examined in two load configurations, one approximating joint loading during single-limb stance and one simulating impact from a fall. Mechanical testing was performed to assess model and failure theory performance in the context of predicting femoral fracture load. Measured and FE-computed fracture load were significantly correlated for both loading conditions and all failure criteria (p < or = 0.001). The distortion energy and tau(max) failure theories were the most robust of those examined, providing the most consistently strong FE model performance for two very different loading conditions. The more complex failure theories and the strain-based theories examined did not improve performance over the simpler distortion energy and tau(max) theories, and often degraded performance, even when differences between tensile and compressive failure properties were represented. The relatively strong performance of the distortion energy and tau(max) theories supports the hypothesis that shear/distortion is an important failure mode during femoral fracture.     

Real-time subject-specific monitoring of internal deformations and stresses in the soft tissues of the foot: a new approach in gait analysis

No technology is presently available to provide real-time information on internal deformations and stresses in plantar soft tissues of individuals during evaluation of the gait pattern. Because internal deformations and stresses in the plantar pad are critical factors in foot injuries such as diabetic foot ulceration, this severely limits evaluation of patients. To allow such real-time subject-specific analysis, we developed a hierarchal modeling system which integrates a two-dimensional gross structural model of the foot (high-order model) with local finite element (FE) models of the plantar tissue padding the calcaneus and medial metatarsal heads (low-order models). The high-order whole-foot model provides real-time analytical evaluations of the time-dependent plantar fascia tensile forces during the stance phase. These force evaluations are transferred, together with foot-shoe local reaction forces, also measured in real time (under the calcaneus, medial metatarsals and hallux), to the low-order FE models of the plantar pad, where they serve as boundary conditions for analyses of local deformations and stresses in the plantar pad. After careful verification of our custom-made FE solver and of our foot model system with respect to previous literature and against experimental results from a synthetic foot phantom, we conducted human studies in which plantar tissue loading was evaluated in real time during treadmill gait in healthy individuals (N = 4). We concluded that internal deformations and stresses in the plantar pad during gait cannot be predicted from merely measuring the foot-shoe force reactions. Internal loading of the plantar pad is constituted by a complex interaction between the anatomical structure and mechanical behavior of the foot skeleton and soft tissues, the body characteristics, the gait pattern and footwear. Real-time FE monitoring of internal deformations and stresses in the plantar pad is therefore required to identify elevated deformation/stress exposures toward utilizing it in gait laboratories to protect feet that are susceptible to injury.     

Effect of boundary conditions,impact loading and hydraulic stiffening on femoral fracture strength

Patient specific quantitative CT (QCT) imaging data together with the finite element (FE) method may provide an accurate prediction of a patient's femoral strength and fracture risk. Although numerous FE models investigating femoral fracture strength have been published, there is little consent on the effect of boundary conditions, dynamic loading and hydraulic strengthening due to intra-medullary pressure on the predicted fracture strength. We developed a QCT-derived FE model of a proximal femur that included node-specific modulus assigned based on the local bone density. The effect of three commonly used boundary conditions published in literature were investigated by comparing the resulting strain field due to an applied fracture load. The models were also augmented with viscoelastic material properties and subject to a realistic impact load profile to determine the effect of dynamic loads on the strain field. Finally, the effect of hydraulic strengthening was investigated by including node specific permeability and performing a coupled pore diffusion and stress analysis of the FE model. Results showed that all boundary conditions yield the same strain field patterns, but peak strains were 22% lower and fracture load was 18% higher when loaded at the greater trochanter than when loaded at the femoral head. Comparison of the dynamic models showed that material viscoelasticity was important, but inertial effects (vibration and shock) were not. Finally, pore pressure changes did not cause significant hydraulic strengthening of bone under fall impact loading.     

Relation between subject-specific hip joint loading, stress distribution in the proximal femur and bone mineral density changes after total hip replacement

In the prediction of bone remodelling processes after total hip replacement (THR), modelling of the subject-specific geometry is now state-of-the-art. In this study, we demonstrate that inclusion of subject-specific loading conditions drastically influences the calculated stress distribution, and hence influences the correlation between calculated stress distributions and changes in bone mineral density (BMD) after THR.For two patients who received cementless THR, personalized finite element (FE) models of the proximal femur were generated representing the pre- and post-operative geometry. FE analyses were performed by imposing subject-specific three-dimensional hip joint contact forces as well as muscle forces calculated based on gait analysis data. Average values of the von Mises stress were calculated for relevant zones of the proximal femur. Subsequently, the load cases were interchanged and the effect on the stress distribution was evaluated. Finally, the subject-specific stress distribution was correlated to the changes in BMD at 3 and 6 months after THR.We found subject-specific differences in the stress distribution induced by specific loading conditions, as interchanging of the loading also interchanged the patterns of the stress distribution. The correlation between the calculated stress distribution and the changes in BMD were affected by the two-dimensional nature of the BMD measurement.Our results confirm the hypothesis that inclusion of subject-specific hip contact forces and muscle forces drastically influences the stress distribution in the proximal femur. In addition to patient-specific geometry, inclusion of patient-specific loading is, therefore, essential to obtain accurate input for the analysis of stress distribution after THR.     

Specimen-specific predictions of contact stress under physiological loading in the human hip: validation and sensitivity studies

Hip osteoarthritis may be initiated and advanced by abnormal cartilage contact mechanics, and finite element (FE) modeling provides an approach with the potential to allow the study of this process. Previous FE models of the human hip have been limited by single specimen validation and the use of quasi-linear or linear elastic constitutive models of articular cartilage. The effects of the latter assumptions on model predictions are unknown, partially because data for the instantaneous behavior of healthy human hip cartilage are unavailable. The aims of this study were to develop and validate a series of specimen-specific FE models, to characterize the regional instantaneous response of healthy human hip cartilage in compression, and to assess the effects of material nonlinearity, inhomogeneity and specimen-specific material coefficients on FE predictions of cartilage contact stress and contact area. Five cadaveric specimens underwent experimental loading, cartilage material characterization and specimen-specific FE modeling. Cartilage in the FE models was represented by average neo-Hookean, average Veronda Westmann and specimen- and region-specific Veronda Westmann hyperelastic constitutive models. Experimental measurements and FE predictions compared well for all three cartilage representations, which was reflected in average RMS errors in contact stress of less than 25 %. The instantaneous material behavior of healthy human hip cartilage varied spatially, with stiffer acetabular cartilage than femoral cartilage and stiffer cartilage in lateral regions than in medial regions. The Veronda Westmann constitutive model with average material coefficients accurately predicted peak contact stress, average contact stress, contact area and contact patterns. The use of subject- and region-specific material coefficients did not increase the accuracy of FE model predictions. The neo-Hookean constitutive model underpredicted peak contact stress in areas of high stress. The results of this study support the use of average cartilage material coefficients in predictions of cartilage contact stress and contact area in the normal hip. The regional characterization of cartilage material behavior provides the necessary inputs for future computational studies, to investigate other mechanical parameters that may be correlated with OA and cartilage damage in the human hip. In the future, the results of this study can be applied to subject-specific models to better understand how abnormal hip contact stress and contact area contribute to OA.     

Influence of muscle groups’ activation on proximal femoral growth tendency

Muscle and joint contact force influence stresses at the proximal growth plate of the femur and thus bone growth, affecting the neck shaft angle (NSA) and femoral anteversion (FA). This study aims to illustrate how different muscle groups’ activation during gait affects NSA and FA development in able-bodied children. Subject-specific femur models were developed for three able-bodied children (ages 6, 7, and 11 years) using magnetic resonance images. Contributions of different muscle groups—hip flexors, hip extensors, hip adductors, hip abductors, and knee extensors—to overall hip contact force were computed. Specific growth rate for the growth plate was computed, and the growth was simulated in the principal stress direction at each element in the growth front. The predicted growth indicated decreased NSA and FA (of about \(0.1 {^{\circ }}\) over a four-month period) for able-bodied children. Hip abductors contributed the most, and hip adductors, the least, to growth rate. All muscles groups contributed to a decrease in predicted NSA (\(\sim \)0.01\({^{\circ }}\)–0.04\({^{\circ }})\) and FA (\(\sim \)0.004\({^{\circ }}\)–\(0.2{^{\circ }}\)), except hip extensors and hip adductors, which showed a tendency to increase the FA (\(\sim \)0.004\({^{\circ }}\)–\(0.02{^{\circ }}\)). Understanding influences of different muscle groups on long bone growth tendency can help in treatment planning for growing children with affected gait.     

Principles of determination and verification of muscle forces in the human musculoskeletal system: Muscle forces to minimise bending stress

While there are a growing number of increasingly complex methodologies available to model geometry and material properties of bones, these models still cannot accurately describe physical behaviour of the skeletal system unless the boundary conditions, especially muscular loading, are correct. Available in vivo measurements of muscle forces are mostly highly invasive and offer no practical way to validate the outcome of any computational model that predicts muscle forces. However, muscle forces can be verified indirectly using the fundamental property of living tissue to functional adaptation and finite element (FE) analysis. Even though the mechanisms of the functional adaptation are not fully understood, its result is clearly seen in the shape and inner structure of bones. The FE method provides a precise tool for analysis of the stress/strain distribution in the bone under given loading conditions. The present work sets principles for the determination of the muscle forces on the basis of the widely accepted view that biological systems are optimized light-weight structures with minimised amount of unloaded/underloaded material and hence evenly distributed loading throughout the structure. Bending loading of bones is avoided/compensated in bones under physiological loading. Thus, bending minimisation provides the basis for the determination of the musculoskeletal system loading. As a result of our approach, the muscle forces for a human femur during normal gait and sitting down (peak hip joint force) are obtained such that the bone is loaded predominantly in compression and the stress distribution in proximal and diaphyseal femur corresponds to the material distribution in bone.     

Experimentally tested computer modeling of stress fractures in rats

The objective of this study was to develop a finite-element (FE) modeling methodology for studying the etiology of a stress fracture (SF). Several variants of three-dimensional FE models of a rat hindlimb, which differed in length or stiffness of tissues, enabling the analyses of mechanical strains and stress in the tibia, were created. We compared the occurrence of SFs in an animal model to validate locations of peak strains/stresses in the FE models. Four Sprague-Dawley male rats, age ~7 wk, were subjected to mechanical cyclic loads of 1.2 Hz and ~6 N, which were delivered to their hindlimb for 30 min, 3 times/wk, up to 12 wk, by using a specially designed apparatus. The results showed that 1) FE modeling predicted the maximal strains/stresses (~220,0 με and ~29 MPa, respectively) between the mid- and proximal thirds of the tibia; 2) in a longer shin, greater and more inhomogeneous tensile strains/stresses were evident, at the same location; 3) anatomical variants in shin length influenced the strain/stress distributions to a greater extent with respect to changes in mechanical properties of tissues; and 4) bone stiffness was more dominant than muscle stiffness in affecting the strain/stress distributions. In the animal study, 35,000 loading cycles were associated with the formation of a SF. The location of the identified SF in the rat limb verified the FE model. We find the suggested model a valuable tool in studying various aspects of SFs.     

A finite element model of the lower limb during stance phase of gait cycle including the muscle forces

Background

Results of finite element (FE) analyses can give insight into musculoskeletal diseases if physiological boundary conditions, which include the muscle forces during specific activities of daily life, are considered in the FE modelling. So far, many simplifications of the boundary conditions are currently made. This study presents an approach for FE modelling of the lower limb for which muscle forces were included.

Methods

The stance phase of normal gait was simulated. Muscle forces were calculated using a musculoskeletal rigid body (RB) model of the human body, and were subsequently applied to a FE model of the lower limb. It was shown that the inertial forces are negligible during the stance phase of normal gait. The contact surfaces between the parts within the knee were modelled as bonded. Weak springs were attached to the distal tibia for numerical reasons.

Results

Hip joint reaction forces from the RB model and those from the FE model were similar in magnitude with relative differences less than 16%. The forces of the weak spring were negligible compared to the applied muscle forces. The maximal strain was 0.23% in the proximal region of the femoral diaphysis and 1.7% in the contact zone between the tibia and the fibula.

Conclusions

The presented approach based on FE modelling by including muscle forces from inverse dynamic analysis of musculoskeletal RB model can be used to perform analyses of the lower limb with very realistic boundary conditions. In the present form, this model can be used to better understand the loading, stresses and strains of bones in the knee area and hence to analyse osteotomy fixation devices.

     

Developmental Plasticity in the Ossification of the Proximal Femur of Heterocephalus glaber (Bathyergidae,Rodentia)

The proximal femoral morphology in rodents of different body sizes, locomotor modes, and from the three main rodent lineages (Sciuromorpha, Myomorpha, and Hystricomorpha) exhibit a separated condition of the femoral head and greater trochanter. We assessed the femoral ossification of eight species of all six genera of a subterranean lineage of mammals, the African mole-rats (Bathyergidae), including the naked mole-rat (Heterocephalus glaber). Here we report a surprising level of intraspecific variation in the ossification of the proximal femur of H. glaber, which presents both separated and coalesced conditions, regardless of sex and reproductive status. The other bathyergids, including chisel-tooth and scratch-diggers exhibit a separated condition, similar to the typical rodent condition. Because the coalesced condition is uncommon among rodents, our data suggests that the presence of two femoral morphologies in H. glaber represent developmental plasticity in this species. Such a dual condition may result from a constricted femoral head and greater trochanter morphology and slow skeletal growth rates, which could be also influenced by differential loading histories, such as magnitude and orientation of forces acting on the limb during ontogeny. This is the best documented case of intraspecific variation for this trait amongst non-human vertebrates, and its investigation is important to understanding the mechanisms of skeletal development and phenotypic plasticity in mammals.

     

3D finite element model of the diabetic neuropathic foot: A gait analysis driven approach

Diabetic foot is an invalidating complication of diabetes that can lead to foot ulcers. Three-dimensional (3D) finite element analysis (FEA) allows characterizing the loads developed in the different anatomical structures of the foot in dynamic conditions. The aim of this study was to develop a subject specific 3D foot FE model (FEM) of a diabetic neuropathic (DNS) and a healthy (HS) subject, whose subject specificity can be found in term of foot geometry and boundary conditions. Kinematics, kinetics and plantar pressure (PP) data were extracted from the gait analysis trials of the two subjects with this purpose. The FEM were developed segmenting bones, cartilage and skin from MRI and drawing a horizontal plate as ground support. Materials properties were adopted from previous literature. FE simulations were run with the kinematics and kinetics data of four different phases of the stance phase of gait (heel strike, loading response, midstance and push off). FEMs were then driven by group gait data of 10 neuropathic and 10 healthy subjects. Model validation focused on agreement between FEM-simulated and experimental PP.     

Leg and Joint Stiffness in Children with Spastic Diplegic Cerebral Palsy during Level Walking

Individual joint deviations are often identified in the analysis of cerebral palsy (CP) gait. However, knowledge is limited as to how these deviations affect the control of the locomotor system as a whole when striving to meet the demands of walking. The current study aimed to bridge the gap by describing the control of the locomotor system in children with diplegic CP in terms of their leg stiffness, both skeletal and muscular components, and associated joint stiffness during gait. Twelve children with spastic diplegia CP and 12 healthy controls walked at a self-selected pace in a gait laboratory while their kinematic and forceplate data were measured and analyzed during loading response, mid-stance, terminal stance and pre-swing. For calculating the leg stiffness, each of the lower limbs was modeled as a non-linear spring, connecting the hip joint center and the corresponding center of pressure, with varying stiffness that was calculated as the slope (gradient) of the axial force vs. the deformation curve. The leg stiffness was further decomposed into skeletal and muscular components considering the alignment of the lower limb. The ankle, knee and hip of the limb were modeled as revolute joints with torsional springs whose stiffness was calculated as the slope of the moment vs. the angle curve of the joint. Independent t-tests were performed for between-group comparisons of all the variables. The CP group significantly decreased the leg stiffness but increased the joint stiffness during stance phase, except during terminal stance where the leg stiffness was increased. They appeared to rely more on muscular contributions to achieve the required leg stiffness, increasing the muscular demands in maintaining the body posture against collapse. Leg stiffness plays a critical role in modulating the kinematics and kinetics of the locomotor system during gait in the diplegic CP.     

Effect of varus/valgus malalignment on bone strains in the proximal tibia after TKR: an explicit finite element study

Malalignment is the main cause of tibial component loosening. Implants that migrate rapidly in the first two post-operative years are likely to present aseptic loosening. It has been suggested that cancellous bone stresses can be correlated with tibial component migration. A recent study has shown that patient-specific finite element (FE) models have the power to predict the short-term behavior of tibial trays. The stresses generated within the implanted tibia are dependent on the kinematics of the joint; however, previous studies have ignored the kinematics and only applied static loads. Using explicit FE, it is possible to simultaneously predict the kinematics and stresses during a gait cycle. The aim of this study was to examine the cancellous bone strains during the stance phase of the gait cycle, for varying degrees of varus/valgus eccentric loading using explicit FE. A patient-specific model of a proximal tibia was created from CT scan images, including heterogeneous bone properties. The proximal tibia was implanted with a commercial total knee replacement (TKR) model. The stance phase of gait was simulated and the applied loads and boundary conditions were based on those used for the Stanmore knee simulator. Eccentric loading was simulated. As well as examining the tibial bone strains (minimum and maximum principal strain), the kinematics of the bone-implant construct are also reported. The maximum anterior-posterior displacements and internal-external rotations were produced by the model with 20 mm offset. The peak minimum and maximum principal strain values increased as the load was shifted laterally, reaching a maximum magnitude for -20 mm offset. This suggests that when in varus, the load transferred to the bone is shifted medially, and as the bone supporting this load is stiffer, the resulting peak bone strains are lower than when the load is shifted laterally (valgus). For this particular patient, the TKR design analyzed produced the highest cancellous bone strains when in valgus. This study has provided an insight in the variations produced in bone strain distribution when the axial load is applied eccentrically. To the authors' knowledge, this is the first time that the bone strain distribution of a proximal implanted tibia has been examined, also accounting for the kinematics of the tibio-femoral joint as part of the simulation. This approach gives greater insight into the overall performance of TKR.     

Growth plate mechanics and mechanobiology. A survey of present understanding

The longitudinal growth of long bones occurs in growth plates where chondrocytes synthesize cartilage that is subsequently ossified. Altered growth and subsequent deformity resulting from abnormal mechanical loading is often referred to as mechanical modulation of bone growth. This phenomenon has key implications in the progression of infant and juvenile musculoskeletal deformities, such as adolescent idiopathic scoliosis, hyperkyphosis, genu varus/valgus and tibia vara/valga, as well as neuromuscular diseases. Clinical management of these deformities is often directed at modifying the mechanical environment of affected bones. However, there is limited quantitative and physiological understanding of how bone growth is regulated in response to mechanical loading. This review of published work addresses the state of knowledge concerning key questions about mechanisms underlying biomechanical modulation of bone growth. The longitudinal growth of bones is apparently controlled by modifying the numbers of growth plate chondrocytes in the proliferative zone, their rate of proliferation, the amount of chondrocytic hypertrophy and the controlled synthesis and degradation of matrix throughout the growth plate. These variables may be modulated to produce a change in growth rate in the presence of sustained or cyclic mechanical load. Tissue and cellular deformations involved in the transduction of mechanical stimuli depend on the growth plate tissue material properties that are highly anisotropic, time-dependent, and that differ in different zones of the growth plate and with developmental stages. There is little information about the effects of time-varying changes in volume, water content, osmolarity of matrix, etc. on differentiation, maturation and metabolic activity of chondrocytes. Also, the effects of shear forces and torsion on the growth plate are incompletely characterized. Future work on growth plate mechanobiology should distinguish between changes in the regulation of bone growth resulting from different processes, such as direct stimulation of the cell nuclei, physico-chemical stimuli, mechanical degradation of matrix or cellular components and possible alterations of local blood supply.