Latest Articles Include:
- Editorial Board and publication information
- J Biomech 43(5):IFC (2010)
- Directing bone marrow-derived stromal cell function with mechanics
Potier E Noailly J Ito K - J Biomech 43(5):807-817 (2010)
Because bone marrow-derived stromal cells (BMSCs) are able to generate many cell types, they are envisioned as source of regenerative cells to repair numerous tissues, including bone, cartilage, and ligaments. Success of BMSC-based therapies, however, relies on a number of methodological improvements, among which better understanding and control of the BMSC differentiation pathways. Since many years, the biochemical environment is known to govern BMSC differentiation, but more recent evidences show that the biomechanical environment is also directing cell functions. Using in vitro systems that aim to reproduce selected components of the in vivo mechanical environment, it was demonstrated that mechanical loadings can affect BMSC proliferation and improve the osteogenic, chondrogenic, or myogenic phenotype of BMSCs. These effects, however, seem to be modulated by parameters other than mechanics, such as substrate nature or soluble biochemical environment. This paper revie ws and discusses recent experimental data showing that despite some knowledge limitation, mechanical stimulation already constitutes an additional and efficient tool to drive BMSC differentiation. - The effect of positioning on the biomechanical performance of soft shell hip protectors
Choi WJ Hoffer JA Robinovitch SN - J Biomech 43(5):818-825 (2010)
Wearable hip protectors represent a promising strategy for reducing risk for hip fracture from a sideways fall. However, small changes in pad positioning may influence their protective benefit. Using a mechanical hip impact simulator, we investigated how three marketed soft shell hip protectors attenuate and redistribute the impact force applied to the hip, and how this depends on displacement from their intended position by 2.5 or 5 cm superiorly, posteriorly, inferiorly or anteriorly. For centrally-placed protectors, peak pressure was reduced 93% below the unpadded value by a 16 mm horseshoe-shaped protector, 93% by a 14 mm horseshoe protector, and 94% by a 16 mm continuous protector. In unpadded trials, 83% of the total force was applied to the skin overlying the proximal femur (danger zone). This was lowered to 19% by the centrally placed 16 mm horseshoe protector, to 34% by the 14 mm horseshoe, and to 40% by the 16 mm continuous protector. Corresponding reductions in peak force delivered to the femoral neck (relative to unpadded) were 45%, 38%, and 20%, respectively. The protective benefit of all three protectors decreased with pad displacement. For example, displacement of protectors by 5 cm anteriorly caused peak femoral neck force to increase 60% above centrally-placed values, and approach unpadded values. These results indicate that soft shell hip protectors provide substantial protective benefits, but decline in performance with small displacements from their intended position. Our findings confirm the need for correct and stable positioning of hip protectors in garment design. - Structural behaviour and strain distribution of the long bones of the human lower limbs
Cristofolini L Conti G Juszczyk M Cremonini S Sint Jan SV Viceconti M - J Biomech 43(5):826-835 (2010)
Although stiffness and strength of lower limb bones have been investigated in the past, information is not complete. While the femur has been extensively investigated, little information is available about the strain distribution in the tibia, and the fibula has not been tested in vitro. This study aimed at improving the understanding of the biomechanics of lower limb bones by: (i) measuring the stiffness and strain distributions of the different low limb bones; (ii) assessing the effect of viscoelasticity in whole bones within a physiological range of strain-rates; (iii) assessing the difference in the behaviour in relation to opposite directions of bending and torsion. The structural stiffness and strain distribution of paired femurs, tibias and fibulas from two donors were measured. Each region investigated of each bone was instrumented with 8–16 triaxial strain gauges (over 600 grids in total). Each bone was subjected to 6–12 different loading configurations. Te sts were replicated at two different loading speeds covering the physiological range of strain-rates. Viscoelasticity did not have any pronounced effect on the structural stiffness and strain distribution, in the physiological range of loading rates explored in this study. The stiffness and strain distribution varied greatly between bone segments, but also between directions of loading. Different stiffness and strain distributions were observed when opposite directions of torque or opposite directions of bending (in the same plane) were applied. To our knowledge, this study represents the most extensive collection of whole-bone biomechanical properties of lower limb bones. - Biomechanical model of human cornea based on stromal microstructure
Studer H Larrea X Riedwyl H Büchler P - J Biomech 43(5):836-842 (2010)
The optical characteristics of the human cornea depends on the mechanical balance between the intra-ocular pressure and intrinsic tissue stiffness. A wide range of ophthalmic surgical procedures alter corneal biomechanics to induce local or global curvature changes for the correction of visual acuity. Due to the large number of surgical interventions performed every day, a deeper understanding of corneal biomechanics is needed to improve the safety of these procedures and medical devices. The aim of this study is to propose a biomechanical model of the human cornea, based on stromal microstructure. The constitutive mechanical law includes collagen fiber distribution based on X-ray scattering analysis, collagen cross-linking, and fiber uncrimping. Our results showed that the proposed model reproduced inflation and extensiometry experimental data [Elsheikh et al., Curr. Eye Res., 2007; Elsheikh et al., Exp. Eye Res., 2008] successfully. The mechanical properties obtained for different age groups demonstrated an increase in collagen cross-linking for older specimens. In future work such a model could be used to simulate non-symmetric interventions, and provide better surgical planning. - Optimal handgrip height of four-wheeled walker on various road conditions to reduce muscular load for elderly users with steady walking
Takanokura M - J Biomech 43(5):843-848 (2010)
A four-wheeled walker is a valuable tool for assisting elderly persons with walking. The handgrip height is one of the most important factor determining the usefulness of the walker. However, the optimal handgrip height for elderly users has not been considered from a biomechanical viewpoint. In this study, the handgrip height was optimized by a two-dimensional mechanical model to reduce muscular loads in the lower body as well as in the upper body with various road conditions during steady walking. A critical height of the handgrip existed at 48% of the body height for the user regardless of gender and body dimension. A lower handgrip relieved muscular load for stooping users with a lower standing height. The stooping user pushed the handgrip strongly in the perpendicular direction by leaning the upper body on the walker. However, upright users with a higher standing height should use a four-wheeled walker with a higher handgrip for maintaining his or her upright postu re. For downhill movement, the optimal handgrip height depended on the slope angle and the friction coefficient between the road and the wheels of the walker. On a low-friction downhill such as asphalt with a steeper slope angle, the user was required to maintain an erect trunk with a higher handgrip and to press on the handgrip strongly in the perpendicular direction. Movement on a low-friction road was easier for users on a flat road and an uphill road, but it compelled distinct effort from users when moving downhill. - Contact mechanics and elastohydrodynamic lubrication in a novel metal-on-metal hip implant with an aspherical bearing surface
Meng Q Gao L Liu F Yang P Fisher J Jin Z - J Biomech 43(5):849-857 (2010)
Diameter and diametral clearance of the bearing surfaces of metal-on-metal hip implants and structural supports have been recognised as key factors to reduce the dry contact and hydrodynamic pressures and improve lubrication performance. On the other hand, application of aspherical bearing surfaces can also significantly affect the contact mechanics and lubrication performance by changing the radius of the curvature of a bearing surface and consequently improving the conformity between the head and the cup. In this study, a novel metal-on-metal hip implant employing a specific aspherical bearing surface, Alpharabola, as the acetabular surface was investigated for both contact mechanics and elastohydrodynamic lubrication under steady-state conditions. When compared with conventional spherical bearing surfaces, a more uniform pressure distribution and a thicker lubricant film thickness within the loaded conjunction were predicted for this novel Alpharabola hip implant. Th e effects of the geometric parameters of this novel acetabular surface on the pressure distribution and lubricant thickness were investigated. A significant increase in the predicted lubricant film thickness and a significant decrease in the dry contact and hydrodynamic pressures were found with appropriate combinations of these geometric parameters, compared with the spherical bearing surface. - Effect of acetabular labral tears, repair and resection on hip cartilage strain: A 7 T MR study
Greaves LL Gilbart MK Yung AC Kozlowski P Wilson DR - J Biomech 43(5):858-863 (2010)
Background Tears of the acetabular labrum are frequently present in patients with groin pain. While it is clear that the labrum contributes to the surface area articulating with the femoral head, it is not clear whether labral repair yields different load distribution in the hip compared to labral resection. Purpose Determine whether labral repair reduces cartilage strain more effectively than labral resection. Methods Six human cadaveric hips (mean age 37 years) were loaded in a simulated single-leg stance within the bore of a 7 T MR scanner. After cartilage had reached a steady-state thickness distribution, a scan of the cartilage was acquired with a voxel size of 0.1×0.1×0.3 mm. This method was repeated for each of six specimens when the labrum was intact, after a surgically simulated labral tear, after an arthroscopic labral repair and after labral resection. Cartilage thickness and strain in an anterosuperior region of interest were measured from the MR scans. A paired t-test was used to compare mean and maximum cartilage strain when the labrum was intact vs. torn, torn vs. repaired and repaired vs. resected. Three-dimensional patterns of cartilage strain distribution were qualitatively compared for the different labral conditions. Results For the number of specimens tested we found no change in mean and maximum cartilage strain, and little obvious change in the pattern of cartilage strain distribution after a simulated labral tear. Labral repair caused a 2% decrease in mean cartilage strain compared to a torn labrum (p=0.014). Labral resection caused a 4% and 6% increase in mean and maximum cartilage strain, respectively, compared to labral repair (p=0.02), and the cartilage strain distribution was elevated throughout the region of interest. Conclusion Based on our ex vivo findings of increased cartilage strain after labral resection when compared to labral repair, we have demonstrated the associated consequences to the mechanical environment of the cartilage following surgical treatment of the labrum. - No need for particle tracing: From accumulating fluid properties to novel blood coagulation model in the lattice Boltzmann method
Moiseyev G Bar-Yoseph PZ - J Biomech 43(5):864-870 (2010)
The accumulation of a fluid property from the standpoint of a particle moving with non-steady fluid flow (i.e., platelet/blood-cell damage index in pulsating blood flow) is a challenged computational problem due to the current need for particle-tracing methods. The method we developed (dubbed VPI) enables the approximation of the Lagrangian integral in real-time for any point in space and time for the entire domain and which is easily integrated into the the lattice Boltzmann method. As an illustrative numerical example we applied our method to a blood coagulation model which was shown to accurately capture the coagulation characteristics observed in experiments, and therefore opening a door for more detailed study of systems which are currently hard to study using particle tracing methods. - Assessment of amputee socket–stump–residual bone kinematics during strenuous activities using Dynamic Roentgen Stereogrammetric Analysis
Papaioannou G Mitrogiannis C Nianios G Fiedler G - J Biomech 43(5):871-878 (2010)
The design, construction, and fitting of artificial limbs remain to this day an art, dependent on the accumulated expertise of the practitioner/prosthetist. Socket fitting is cost ineffective, time consuming, and a source of inconvenience for the amputee. Stump–skin slippage within the socket can cause discomfort, internal limb pain, and eventually skin ulcers as a result of excessive pressure and shear within the socket. This study presents a new method of assessment of three-dimensional (3D) socket–stump kinematics/slippage of strenuous activities using Biplane Dynamic Roentgen Stereogrammetric Analysis instrumentation. Ten below knee amputees participated in the study. A more holistic representation of the downward slippage trend of all proximal side skin markers with respect to the socket, and an even more characteristic and of higher magnitude downward-and anterioposterior slippage (maximum slippage: 151 mm for the fast-stop task and 19 mm for the step-down tas k) between the distal markers after impact, was possible for both tasks for all amputees. Displacement between skin-to-skin marker pairs reached maximum values of approximately 10 mm for the step-down trials and up to 24 mm for the fast stop trials. Maximum skin strain was dependent on the position of the skin markers. Distally positioned skin marker pairs demonstrated mainly anterioposterior displacement between each other (maximum relative strain: 13–14%). Maximum relative strain for the proximal markers was 8–10%. This highly accurate, in-vivo, patient-specific, unobtrusive dynamic information, presented using 3D visualization tools that were up to now unavailable to the clinician-prosthetist, can significantly impact the iterative cycle of socket fitting and evaluation. - Finite element analysis of barbed sutures in skin and tendon tissues
Ingle NP King MW Zikry MA - J Biomech 43(5):879-886 (2010)
Barbed surgical sutures are a new type of knotless suture that are currently being used clinically in cosmetic and plastic surgery procedures for faster healing and better cosmesis. Clinical studies are also underway to evaluate their performance in other deep tissue applications. However, little is known about their intrinsic mechanical behavior and their interactions with surrounding tissues. The primary objective of the current study was to analyze the mechanical behavior of barbed sutures using a finite element analysis approach. First, the effect of applying a point-pressure load to the tip of the barb and measuring its effect on barb displacement was studied. Second, the effect of an applied displacement to a barb anchored either in skin or tendon material for both the suture and the surrounding tissue. The results indicate that the flexibility of the barb can be increased or decreased by changing the barb geometry. It was concluded that the barb geometry and desi gn need to be modified for use with different types of tissue. For example, in order to achieve the best mechanical anchoring with skin tissue the barb should be more flexible compared to the one designed to work with tendon tissue. The uniqueness of this study is that it is the first to establish a virtual prototyping and designing method for barbed sutures. For example, a new and improved virtual design of barb geometry is proposed and validated. It also provides the first report on how to develop a virtual bench top suture/tissue pullout testing environment. - A cell based modelling framework for skeletal tissue engineering applications
Geris L Van Liedekerke P Smeets B Tijskens E Ramon H - J Biomech 43(5):887-892 (2010)
In this study, a cell based lattice free modelling framework is proposed to study cell aggregate behaviour in bone tissue engineering applications. The model encompasses cell-to-cell and cell–environment interactions such as adhesion, repulsion and drag forces. Oxygen, nutrients, waste products, growth factors and inhibitors are explicitly represented in the model influencing cellular behaviour. Furthermore, a model for cell metabolism is incorporated representing the basic enzymic reactions of glycolysis and the Krebs cycle. Various types of cell death such as necrosis, apoptosis and anoikis are implemented. Finally, an explicit model of the cell cycle controls the proliferation process, taking into account the presence or absence of various metabolites, sufficient space and mechanical stress. Several examples are presented demonstrating the potential of the modelling framework. The behaviour of a synchronised cell aggregate under ideal circumstances is simulated, cl early showing the different stages of the cell cycle and the resulting growth of the aggregate. Also the difference in aggregate development under ideal (normoxic) and hypoxic conditions is simulated, showing hypoxia induced necrosis mainly in the centre of the aggregate grown under hypoxic conditions. The next step in this research will be the application of this modelling framework to specific experimental set-ups for bone tissue engineering applications. - Joint moment contributions to swing knee extension acceleration during gait in children with spastic hemiplegic cerebral palsy
Goldberg EJ Requejo PS Fowler EG - J Biomech 43(5):893-899 (2010)
Inadequate peak knee extension during the swing phase of gait is a major deficit in individuals with spastic cerebral palsy (CP). The biomechanical mechanisms responsible for knee extension have not been thoroughly examined in CP. The purpose of this study was to assess the contributions of joint moments and gravity to knee extension acceleration during swing in children with spastic hemiplegic CP. Six children with spastic hemiplegic CP were recruited (age=13.4±4.8 years). Gait data were collected using an eight-camera system. Induced acceleration analysis was performed for each limb during swing. Average joint moment and gravity contributions to swing knee extension acceleration were calculated. Total swing and stance joint moment contributions were compared between the hemiplegic and non-hemiplegic limbs using paired t-tests (p<0.05). Swing limb joint moment contributions from the hemiplegic limb decelerated swing knee extension significantly more than those of the non-hemiplegic limb and resulted in significantly reduced knee extension acceleration. Total stance limb joint moment contributions were not statistically different. Swing limb joint moment contributions that decelerated knee extension appeared to be the primary cause of inadequate knee extension acceleration during swing. Stance limb muscle strength did not appear to be the limiting factor in achieving adequate knee extension in children with CP. Recent research has shown that the ability to extend the knee during swing is dependent on the selective voluntary motor control of the limb. Data from individual participants support this concept. - Cells in 3D matrices under interstitial flow: Effects of extracellular matrix alignment on cell shear stress and drag forces
Pedersen JA Lichter S Swartz MA - J Biomech 43(5):900-905 (2010)
Interstitial flow is an important regulator of various cell behaviors both in vitro and in vivo, yet the forces that fluid flow imposes on cells embedded in a 3D extracellular matrix (ECM), and the effects of matrix architecture on those forces, are not well understood. Here, we demonstrate how fiber alignment can affect the shear and pressure forces on the cell and ECM. Using computational fluid dynamics simulations, we show that while the solutions of the Brinkman equation accurately estimate the average fluid shear stress and the drag forces on a cell within a 3D fibrous medium, the distribution of shear stress on the cellular surface as well as the peak shear stresses remain intimately related to the pericellular fiber architecture and cannot be estimated using bulk-averaged properties. We demonstrate that perpendicular fiber alignment of the ECM yields lower shear stress and pressure forces on the cells and higher stresses on the ECM, leading to decreased permeabil ity, while parallel fiber alignment leads to higher stresses on cells and increased permeability, as compared to a cubic lattice arrangement. The Spielman–Goren permeability relationships for fibrous media agreed well with CFD simulations of flow with explicitly considered fibers. These results suggest that the experimentally observed active remodeling of ECM fibers by fibroblasts under interstitial flow to a perpendicular alignment could serve to decrease the shear and drag forces on the cell. - Mathematical modeling and simulation of seated stability
Tanaka ML Ross SD Nussbaum MA - J Biomech 43(5):906-912 (2010)
Various methods have been used to quantify the kinematic variability or stability of the human spine. However, each of these methods evaluates dynamic behavior within the stable region of state space. In contrast, our goal was to determine the extent of the stable region. A 2D mathematical model was developed for a human sitting on an unstable seat apparatus (i.e., the "wobble chair"). Forward dynamic simulations were used to compute trajectories based on the initial state. From these trajectories, a scalar field of trajectory divergence was calculated, specifically a finite time Lyapunov exponent (FTLE) field. Theoretically, ridges of local maxima within this field are expected to partition the state space into regions of qualitatively different behavior. We found that ridges formed at the boundary between regions of stability and failure (i.e., falling). The location of the basin of stability found using the FTLE field matched well with the basin of stability dete rmined by an alternative method. In addition, an equilibrium manifold was found, which describes a set of equilibrium configurations that act as a low dimensional attractor in the controlled system. These simulations are a first step in developing a method to locate state space boundaries for torso stability. Identifying these boundaries may provide a framework for assessing factors that contribute to health risks associated with spinal injury and poor balance recovery (e.g., age, fatigue, load/weight, and distribution). Furthermore, an approach is presented that can be adapted to find state space boundaries in other biomechanical applications. - Muscle fatigue does not lead to increased instability of upper extremity repetitive movements
Gates DH Dingwell JB - J Biomech 43(5):913-919 (2010)
Muscle fatigue alters neuromuscular responses. This may lead to increased sensitivity to perturbations and possibly to subsequent injury risk. We studied the effects of muscle fatigue on movement stability during a repetitive upper extremity task. Twenty healthy young subjects performed a repetitive work task, similar to sawing, synchronized with a metronome before and after performing each of two fatiguing tasks. The first fatigue task (LIFT) primarily fatigued the shoulder flexor muscles, while the second fatigue task (SAW) fatigued all of the muscles of the arm. Subjects performed each task in random order on two different days at least seven days apart. Instantaneous mean EMG frequencies (IMNF) decreased over both fatiguing tasks indicating that subjects did experience significant muscle fatigue. The slopes of the IMNF over time and the decreases in maximum force measurements demonstrated that the LIFT fatigue task successfully fatigued the shoulder flexors to a gre ater extent than any other muscle. On average, subjects exhibited more locally stable shoulder movements after the LIFT fatigue task (p=0.035). They also exhibited more orbitally stable shoulder (p=0.021) and elbow (p=0.013) movements after the SAW fatigue task. Subjects also had decreased cocontraction at the wrist post-fatigue for both tasks (p=0.001) and at the shoulder (p<0.001) for the LIFT fatigue task. Therefore, increased dynamic stability of these repeated movements cannot be explained by increased muscle cocontraction. Possible alternative mechanisms are discussed. - Functional properties of cartilaginous tissues engineered from infrapatellar fat pad-derived mesenchymal stem cells
Buckley CT Vinardell T Thorpe SD Haugh MG Jones E McGonagle D Kelly DJ - J Biomech 43(5):920-926 (2010)
Articular cartilage has a poor intrinsic capacity for self-repair. The advent of autologous chondrocyte implantation has provided a feasible method to treat cartilage defects. However, the associated drawbacks with the isolation and expansion of chondrocytes from autologous tissue has prompted research into alternative cell sources such as mesenchymal stem cells (MSCs) which have been found to exist in the bone marrow as well as other joint tissues such as the infrapatellar fat pad (IFP), synovium and within the synovial fluid itself. In this work we assessed the chondrogenic potential of IFP-derived porcine cells over a 6 week period in agarose hydrogel culture in terms of mechanical properties, biochemical content and histology. It was found that IFP cells underwent robust chondrogenesis as assessed by glycosaminoglycan (1.47±0.22% w/w) and collagen (1.44±0.22% w/w) accumulation after 42 days of culture. The 1 Hz dynamic modulus of the engineered tissue at this time point was 272.8 kPa (±46.8). The removal of TGF-β3 from culture after 21 days was shown to have a significant effect on both the mechanical properties and biochemical content of IFP constructs after 42 days, with minimal increases occurring from day 21 to day 42 without continued supplementation of TGF-β3. These findings further strengthen the case that the IFP may be a promising cell source for putative cartilage repair strategies. - A strain-hardening bi-power law for the nonlinear behaviour of biological soft tissues
Nicolle S Vezin P Palierne JF - J Biomech 43(5):927-932 (2010)
Biological soft tissues exhibit a strongly nonlinear viscoelastic behaviour. Among parenchymous tissues, kidney and liver remain less studied than brain, and a first goal of this study is to report additional material properties of kidney and liver tissues in oscillatory shear and constant shear rate tests. Results show that the liver tissue is more compliant but more strain hardening than kidney. A wealth of multi-parameter mathematical models has been proposed for describing the mechanical behaviour of soft tissues. A second purpose of this work is to develop a new constitutive law capable of predicting our experimental data in the both linear and nonlinear viscoelastic regime with as few parameters as possible. We propose a nonlinear strain-hardening fractional derivative model in which six parameters allow fitting the viscoelastic behaviour of kidney and liver tissues for strains ranging from 0.01 to 1 and strain rates from 0.0151 s−1 to 0.7 s−1. - Mechanical interaction between cells and fluid for bone tissue engineering scaffold: Modulation of the interfacial shear stress
Blecha LD Rakotomanana L Razafimahery F Terrier A Pioletti DP - J Biomech 43(5):933-937 (2010)
An analytical model of the fluid/cell mechanical interaction was developed. The interfacial shear stress, due to the coupling between the fluid and the cell deformation, was characterized by a new dimensionless number Nfs. For Nfs above a critical value, the fluid/cell interaction had a damping effect on the interfacial shear stress. Conversely, for Nfs below this critical value, interfacial shear stress was amplified. As illustration, the role of the dynamic fluid/cell mechanical coupling was studied in a specific biological situation involving cells seeded in a bone scaffold. For the particular bone scaffold chosen, the dimensionless number Nfs was higher than the critical value. In this case, the dynamic shear stress at the fluid/cell interface is damped for increasing excitation frequency. Interestingly, this damping effect is correlated to the pore diameter of the scaffold, furnishing thus target values in the design of the scaffold. Correspondingly, an efficient c ell stimulation might be achieved with a scaffold of pore size larger than 300 μm as no dynamic damping effect is likely to take place. The analytical model proposed in this study, while being a simplification of a fluid/cell mechanical interaction, brings complementary insights to numerical studies by analyzing the effect of different physical parameters. - The 3D path of body centre of mass during adult human walking on force treadmill
Tesio L Rota V Chessa C Perucca L - J Biomech 43(5):938-944 (2010)
Three-dimensional (3D) path of the body centre of mass (CM) over an entire stride was computed from ground reaction forces during walking at constant average speed on a treadmill mounted on 3D force sensors. Data were obtained from 18 healthy adults at speeds ranging from 0.30 to 1.40 m s−1, in 0.1 m s−1 increments. Six subsequent strides were analyzed for each subject and speed (total strides=1296). The test session lasted about 30 min (10 min for walking). The CM path had an upward concave figure-of-eight shape that was highly consistent within and across subjects. Vertical displacement of the CM increased monotonically as a function of walking speed. The forward and particularly lateral displacements of the CM showed a U-shaped relationship to speed. The same held for the total 3D displacement (25.6–16.0 cm, depending on the speed). The results provide normative benchmarks and suggest hypotheses for further physiologic and clinical research. The familiar invert ed pendulum model might be expanded to gyroscopic, "spin-and-turn" models. Abnormalities of the 3D path might flag motor impairments and recovery. - Simultaneous prediction of muscle and contact forces in the knee during gait
Lin YC Walter JP Banks SA Pandy MG Fregly BJ - J Biomech 43(5):945-952 (2010)
Musculoskeletal models are currently the primary means for estimating in vivo muscle and contact forces in the knee during gait. These models typically couple a dynamic skeletal model with individual muscle models but rarely include articular contact models due to their high computational cost. This study evaluates a novel method for predicting muscle and contact forces simultaneously in the knee during gait. The method utilizes a 12 degree-of-freedom knee model (femur, tibia, and patella) combining muscle, articular contact, and dynamic skeletal models. Eight static optimization problems were formulated using two cost functions (one based on muscle activations and one based on contact forces) and four constraints sets (each composed of different combinations of inverse dynamic loads). The estimated muscle and contact forces were evaluated using in vivo tibial contact force data collected from a patient with a force-measuring knee implant. When the eight optimization pr oblems were solved with added constraints to match the in vivo contact force measurements, root-mean-square errors in predicted contact forces were less than 10 N. Furthermore, muscle and patellar contact forces predicted by the two cost functions became more similar as more inverse dynamic loads were used as constraints. When the contact force constraints were removed, estimated medial contact forces were similar and lateral contact forces lower in magnitude compared to measured contact forces, with estimated muscle forces being sensitive and estimated patellar contact forces relatively insensitive to the choice of cost function and constraint set. These results suggest that optimization problem formulation coupled with knee model complexity can significantly affect predicted muscle and contact forces in the knee during gait. Further research using a complete lower limb model is needed to assess the importance of this finding to the muscle and contact force estimation proce ss. - Compressive axial mechanical properties of rat bone as functions of bone volume fraction, apparent density and micro-ct based mineral density
Cory E Nazarian A Entezari V Vartanians V Müller R Snyder BD - J Biomech 43(5):953-960 (2010)
Mechanical testing has been regarded as the gold standard to investigate the effects of pathologies on the structure–function properties of the skeleton. With recent advances in computing power of personal computers, virtual alternatives to mechanical testing are gaining acceptance and use. We have previously introduced such a technique called structural rigidity analysis to assess mechanical strength of skeletal tissue with defects. The application of this technique is predicated upon the use of relationships defining the strength of bone as a function of its density for a given loading mode. We are to apply this technique in rat models to assess their compressive skeletal response subjected to a host of biological and pharmaceutical stimulations. Therefore, the aim of this study is to derive a relationship expressing axial compressive mechanical properties of rat cortical and cancellous bone as a function of equivalent bone mineral density, bone volume fraction or a pparent density over a range of normal and pathologic bones. We used bones from normal, ovariectomized and partially nephrectomized animals. All specimens underwent micro-computed tomographic imaging to assess bone morphometric and densitometric indices and uniaxial compression to failure. We obtained univariate relationships describing 71–78% of the mechanical properties of rat cortical and cancellous bone based on equivalent mineral density, bone volume fraction or apparent density over a wide range of density and common skeletal pathologies. The relationships reported in this study can be used in the structural rigidity analysis introduced by the authors to provide a non-invasive method to assess the compressive strength of bones affected by pathology and/or treatment options. - Effect of cell seeding and mechanical loading on vascularization and tissue formation inside a scaffold: A mechano-biological model using a lattice approach to simulate cell activity
Checa S Prendergast PJ - J Biomech 43(5):961-968 (2010)
Achieving successful vascularization remains one of the main problems in bone tissue engineering. After scaffold implantation, the growth of capillaries into the porous construct may be too slow to provide adequate nutrients to the cells in the scaffold interior and this inhibits tissue formation in the scaffold core. Often, prior to implantation, a controlled cell culture environment is used to stimulate cell proliferation and, once in place, the mechanical environment acting on the tissue construct is determined by the loading conditions at the implantation site. To what extent do cell seeding conditions and the construct loading environment have an effect on scaffold vascularization and tissue growth? In this study, a mechano-biological model for tissue differentiation and blood vessel growth was used to determine the influence of cell seeding on vascular network development and tissue growth inside a regular-structured bone scaffold under different loading condition s. It is predicted that increasing the number of cells seeded homogeneously reduces the rate of vascularization and the maximum penetration of the vascular network, which in turn reduces bone tissue formation. The seeding of cells in the periphery of the scaffold was predicted to be beneficial for vascularization and therefore for bone growth; however, tissue formation occurred more slowly during the first weeks after implantation compared to homogeneous seeding. Low levels of mechanical loading stimulated bone formation while high levels of loading inhibited bone formation and capillary growth. This study demonstrates the feasibility of computational design approaches for bone tissue engineering. - The influence of sagittal center of pressure offset on gait kinematics and kinetics
Haim A Rozen N Wolf A - J Biomech 43(5):969-977 (2010)
Objectives Kinetic patterns of the lower extremity joints have been shown to be influenced by modification of the location of the center of pressure (CoP) of the foot. The accepted theory is that a shifted location of the CoP alters the distance between the ground reaction force and the center of the joint, thereby modifying torques during gait. Various footwear designs have been reported to significantly alter the magnitude of sagittal joint torques during gait. However, the relationship between the CoP and the kinetic patterns in the sagittal plane has not been examined. The aim of this study was to evaluate the association between the sagittal location of the CoP and gait patterns during gait in healthy men. Methods A foot-worn biomechanical device which allows controlled manipulation of the CoP location was utilized. Fourteen healthy men underwent successive gait analysis with the device set to convey three different sagittal locations of the CoP: neutral, anterior offset and posterior offset. Results CoP translation in the sagittal plane (i.e., from posterior to anterior) significantly related with an ankle dorsiflexion torque and a knee extension torque shift throughout the stance phase. Likewise, an anterior translation of the CoP significantly reduced the extension torque at the hip during pre-swing. Conclusions The study results confirm a direct correlation between sagittal offset of the CoP and the magnitude of joint torques throughout the lower extremity. - Polymerisation stress modelling in acrylic bone cement
Briscoe A New A - J Biomech 43(5):978-983 (2010)
Fatigue failure of the cement mantle has been proposed as one of the failure processes contributing to aseptic loosening of cemented joint replacements. It has also been suggested that fatigue failure is dramatically accelerated by residual stress generated during the cement polymerisation process. Previous computational models of the polymerisation process have investigated only the latter part of polymerisation by assuming both instantaneous hardening of the material (a stress locking point) and that all residual stress results from thermal shrinkage after this stress locking point. In this study, finite element models which use the local degree of polymerisation to calculate material properties and shrinkage have been used to predict residual stresses in two models of total hip replacement cement mantles. Results indicate that the final value of cement mantle stress may not be the highest stresses that the cement is subjected to during the polymerisation process. Two models are presented, a 2-dimensional model, which was adapted from a similar model in the literature (Lennon and Prendergast, 2002) and a 3-dimensional concentric-cylinders model. In both cases a chemical kinetics model was used to predict the progress of the polymerisation reaction and a second linear model used to predict cement mechanical properties and density, and so stress generation and volume change, over time. There was good agreement of the results of the 2D model with its counterpart in the literature. For the 3D model, the final residual stress magnitudes and patterns showed good agreement with similar physical and computational models in the literature. - A constitutive formulation of vascular tissue mechanics including viscoelasticity and softening behaviour
Peña E Alastrué V Laborda A Martínez MA Doblaré M - J Biomech 43(5):984-989 (2010)
Nearly all soft tissues, among which the vascular tissue is included, present a certain degree of viscoelastic response. This behaviour may be attributed in part to fluid transport within the solid matrix, and to the friction between its fluid and solid constituents. After being preconditioned, the tissue displays highly repetitive behaviour, so that it can be considered pseudo-elastic, that is, elastic but behaving differently in loading and unloading. Because of this reason, very few constitutive laws accounting for the viscoelastic behaviour of the tissue have been developed. Nevertheless, the consideration of this inelastic effect is of crucial importance in surgeries—like vascular angioplasty—where the mentioned preconditioning cannot be considered since non-physiological deformation is applied on the vessel which, in addition, can cause damage to the tissue. A new constitutive formulation considering the particular features of the vascular tissue, such as anis otropy, together with these two inelastic phenomena is presented here and used to fit experimental stress–stretch curves from simple tension loading–unloading tests and relaxation test on porcine and ovine vascular samples. - Stress relaxation of human ankles is only minimally affected by knee and ankle angle
Tian M Hoang PD Gandevia SC Bilston LE Herbert RD - J Biomech 43(5):990-993 (2010)
Comprehensive characterization of stress relaxation in musculotendinous structures is needed to create robust models of viscoelastic behavior. The commonly used quasi-linear viscoelastic (QLV) theory requires that the relaxation response be independent of tissue strain (length). This study aims to characterize stress relaxation in the musculotendinous and ligamentous structures crossing the human ankle (ankle-only structures and the gastrocnemius muscle–tendon unit, which crosses the ankle and knee), and to determine whether stress relaxation is independent of the length of these structures. Two experiments were conducted on 8 healthy subjects. The first experiment compared stress relaxation over 10 min at different gastrocnemius muscle–tendon unit lengths keeping the length of ankle-joint only structures fixed. The second experiment compared stress relaxation at different lengths of ankle-joint only structures keeping gastrocnemius muscle–tendon unit length fixed . Stress relaxation data were fitted with a two-term exponential function (T=G0+G1e−λ1t+G2e−λ2t). The first experiment demonstrated a significant effect of gastrocnemius muscle–tendon unit length on G1, and the second experiment demonstrated an effect of the length of ankle-joint only structures on G2, λ1 and λ2 (p<0.05). Nonetheless, the size of effects on stress relaxation was small (ΔG/G<10%), similar to experimental variability. We conclude that stress relaxation in the relaxed human ankle is minimally affected by changing gastrocnemius muscle–tendon unit length or by changing the lengths of ankle-joint only structures. Consequently quasi-linear viscoelastic models of the relaxed human ankle can use a common stress relaxation modulus at different knee and ankle angles with minimal error. - Using relative velocity vectors to reveal axial rotation about the medial and lateral compartment of the knee
Anderst WJ Tashman S - J Biomech 43(5):994-997 (2010)
A new technique is presented that utilizes relative velocity vectors between articulating surfaces to characterize internal/external rotation of the tibio-femoral joint during dynamic loading. Precise tibio-femoral motion was determined by tracking the movement of implanted tantalum beads in high-speed biplane X-rays. Three-dimensional, subject-specific CT reconstructions of the femur and tibia, consisting of triangular mesh elements, were positioned in each analyzed frame. The minimum distance between subchondral bone surfaces was recorded for each mesh element comprising each bone surface, and the relative velocity between these opposing closest surface elements was determined in each frame. Internal/external rotation was visualized by superimposing tangential relative velocity vectors onto bone surfaces at each instant. Rotation about medial and lateral compartments was quantified by calculating the angle between these tangential relative vectors within each compartm ent. Results acquired from 68 test sessions involving 23 dogs indicated a consistent pattern of sequential rotation about the lateral condyle (approximately 60 ms after paw strike) followed by rotation about the medial condyle (approximately 100 ms after paw strike). These results imply that axial knee rotation follows a repeatable pattern within and among subjects. This pattern involves rotation about both the lateral and medial compartments. The technique described can be easily applied to study human knee internal/external rotation during a variety of activities. This information may be useful to define normal and pathologic conditions, to confirm post-surgical restoration of knee mechanics, and to design more realistic prosthetic devices. Furthermore, analysis of joint arthrokinematics, such as those described, may identify changes in joint mechanics associated with joint degeneration. - Changes in force and tendinous tissue elongation during the early phase of tetanic summation in in vivo human tibialis anterior muscle
Ohta Y Shima N Yabe K - J Biomech 43(5):998-1001 (2010)
The purpose of this study was to determine the changes that occur in tendinous tissue properties during the early phase of tetanic summation in the in vivo human tibialis anterior muscle (TA). The torque response and tendinous tissue elongation following single stimuli, two-pulse trains, and three-pulse trains were recorded in the TA during isometric contractions. The elongation, compliance, and lengthening velocity of tendinous tissue were determined by real-time ultrasonography. The contribution of the response to the second stimulation (C2) was obtained by subtracting the response to the single stimulation (C1) from the response of doublet. The third contribution (C3) was obtained by subtracting the response to the doublet from that of the triplet. C2 (7.8±0.5 Nm) and C3 (7.3±0.6 Nm) had torque responses significantly higher than C1 (3.6±0.7 Nm). In contrast, the elongations of tendinous tissue for C2 (2.8±0.4 mm) and C3 (1.7±0.2 mm) were significantly lower tha n for C1 (4.9±0.3 mm), indicating that the summation pattern of tendinous tissue elongation is different from the summation pattern of torque response. In addition, this showed considerable difference both between C1 (0.12±0.01 mm/N; 83±4.6 mm/s) and C2 (0.03±0.005 mm/N; 50±6.3 mm/s) and between C1 and C3 (0.02±0.002 mm/N; 39±6.4 mm/s) in the compliance and lengthening velocity of tendinous tissue. These results suggest that changes in tendinous tissue properties between first and second contraction are related to different summation patterns of force and tendinous tissue elongation during early phase of tetanic summation. - An innovative application of a small-scale motion analysis technique to quantify human skin deformation in vivo
Mahmud J Holt CA Evans SL - J Biomech 43(5):1002-1006 (2010)
This study highlights a new experimental method developed to measure full-field deformation of human skin in vivo. The technique uses a small-scale Qualisys (Sweden) 3D motion capture system and an array of reflective markers placed on the forearm of five healthy volunteers. A load of up to 1.5 N was applied to induce skin deformation by pulling a fine wire attached to the centre of the marker configuration. Loading and marker displacements were recorded simultaneously. 3D marker trajectory data was generated for three different load directions. Tests were repeated to investigate accuracy and repeatability. Calibration results indicate the accuracy of the motion capture system with an average residual of 0.05 mm. The procedure was found to be repeatable and accurate for five repeated tests of measured displacements with a maximum variance of 5%. Experimental data are presented to demonstrate robustness and the ability to produce significant outputs. For all five subject s, at 1 N load, the mean and standard deviations of skin axial and lateral displacements were found to be 11.7±1.6 mm and 12.3±3.3 mm, respectively. The axial displacements ratio (u90/u0) ranges from 0.63 to 1.45 with mean±standard deviation of 0.982±0.34 and 0.982±0.32 for left and right arms, respectively. The experiments generated useful and accurate data that can be used to study the viscoelastic, hyperelastic or anisotropic behaviour of human skin. The measured displacements will be analysed further to determine the mechanical properties of skin using inverse Finite Element Analysis and Ogden model. - Quantitative detection of cartilage surfaces and ligament geometry of the wrist using an imaging cryomicrotome system
Dvinskikh NA Blankevoort L Foumani M Spaan JA Streekstra GJ - J Biomech 43(5):1007-1010 (2010)
Biomechanical models may aid in improving diagnosis and treatment of wrist joint disorders. As input, geometrical information is required for model development. Previous studies acquired some elements of the average wrist joint geometry. However, there is a close geometric functional match between articulating surfaces and ligament geometry. Therefore, biomechanical models need to be fed with the geometric data of individual joints. This study is aimed at acquiring geometric data of cartilage surfaces and ligaments from individual wrist joints by using a cryomicrotome imaging system and the evaluation of inter- and intra-observer variability of the data. The 3D geometry of 30 cartilage surfaces and 15 ligaments in three cadaver wrists was manually detected and quantitatively reconstructed. The inter- and intra-observer variability of the cartilage surface detection was 0.14 and 0.19 mm, respectively. For the position of the radius attachment of the dorsal radiocarpal ligament (DRC), the observer variations were 0.12 and 0.65 mm, for intra-/inter-observer, respectively. For the DRC attachment on the triquetrum, the observer variations were 0.22 and 1.19 mm. Anatomic reconstruction from 3D cryomicrotome images offer a method to obtain unique geometry data of the entire wrist joint for modeling purposes. - A new approach for assigning bone material properties from CT images into finite element models
Chen G Schmutz B Epari D Rathnayaka K Ibrahim S Schuetz MA Pearcy MJ - J Biomech 43(5):1011-1015 (2010)
Generation of subject-specific finite element (FE) models from computed tomography (CT) datasets is of significance for application of the FE analysis to bone structures. A great challenge that remains is the automatic assignment of bone material properties from CT Hounsfield Units into finite element models. This paper proposes a new assignment approach, in which material properties are directly assigned to each integration point. Instead of modifying the dataset of FE models, the proposed approach divides the assignment procedure into two steps: generating the data file of the image intensity of a bone in a MATLAB program and reading the file into ABAQUS via user subroutines. Its accuracy has been validated by assigning the density of a bone phantom into a FE model. The proposed approach has been applied to the FE model of a sheep tibia and its applicability tested on a variety of element types. The proposed assignment approach is simple and illustrative. It can be ea sily modified to fit users' situations.
No comments:
Post a Comment