Latest Articles Include:
- Editorial Board and Publication Information
- J Biomech 43(13):IFC (2010)
- Computational analysis of flow structure and particle deposition in a single asthmatic human airway bifurcation
Zhang H Papadakis G - J Biomech 43(13):2453-2459 (2010)
This paper aims to improve current understanding of flow structure and particle deposition in asthmatic human airways. A single, symmetric airway bifurcation, corresponding to generations 10–11 of Weibel's model, is investigated through validated numerical simulations. The parent airway segment is modelled as a smooth circular tube. The child segments are considered asthmatic and their cross-section is modelled as a constricted tube with sinusoidal folds uniformly distributed along the circumference. The flow structure and particle deposition pattern for normal (i.e., healthy) and asthmatic airway bifurcations are compared and discussed. The numerical results reveal that the secondary flow in the asthmatic airway bifurcation is much stronger than in the healthy one, resulting in higher particle deposition. The effects of size of the lumen area and number of folds on particle deposition and pressure drop are also investigated. It is found that particle deposition ef! ficiency is significantly affected by lumen area of the asthmatic segment (the smaller the lumen area, the higher the particle deposition efficiency). The effect of number of folds is small. Particle deposition efficiency also increases with Reynolds number. The pressure drop in the asthmatic airway bifurcation depends mainly on size of the lumen area. The effect of number of folds becomes important for strongly collapsed airways. - Type and orientation of yielded trabeculae during overloading of trabecular bone along orthogonal directions
Shi X Sherry Liu X Wang X Edward Guo X Niebur GL - J Biomech 43(13):2460-2466 (2010)
Trabecular architecture plays a major role in bone mechanics. Osteoporosis leads to a transition from a plate-like to a more rod-like trabecular morphology, which may contribute to fracture risk beyond that predicted by changes in density. In this study, microstructural finite element analysis results were analyzed using individual trabeculae segmentation (ITS) to identify the type and orientation of trabeculae where tissue yielded during compressive overloads in two orthogonal directions. For both apparent loading conditions, most of the yielded tissue was found in longitudinally oriented plates. However, the primary loading mode of yielded trabeculae was axial compression with superposed bending for on-axis loading in contrast to bending for transverse loading. For either loading direction, most plate-like trabeculae yielded in the same loading mode, regardless of their orientation. In contrast, rods oriented parallel to the loading axis yielded in compression, while! rods oblique or perpendicular to the loading axis yielded in combined bending and tension. The predominance of tissue yielding in plates during both on-axis and transverse overloading explains why on-axis overloading is detrimental to the off-axis mechanical properties. At the same time, a large fraction of the tissue in rod-like trabeculae parallel to the loading direction yielded in both on-axis and transverse loading. Hence, rods may be more likely to be damaged and potentially resorbed by damage mediated remodeling. - Effect of boot shaft stiffness on stability joint energy and muscular co-contraction during walking on uneven surface
Böhm H Hösl M - J Biomech 43(13):2467-2472 (2010)
Increased boot shaft stiffness may have a noticeable impact on the range of motion of the ankle joint. Therefore, the ability of the ankle joint to generate power for propulsion might be impaired. This might result in compensatory changes at the knee and hip joint. Besides, adaptability of the subtalar joint to uneven surface might be reduced, which could in turn affect stability. The aim of the study was therefore to investigate the influence of boot shaft stiffness on biomechanical gait parameters. Fifteen healthy young adults walked over coarse gravel wearing two different hiking boots that differed by 50% in passive shaft stiffness. Leg kinematics, kinetics and electromyography were measured. Gait velocity and indicators for stability were not different when walking with the hard and soft boot shaft over the gravel surface. However, the hard boot shaft decreased the ankle range of motion as well as the eccentric energy absorbed at the ankle joint. As a consequence, compensatory changes at the knee joint were observed. Co-contraction was increased, and greater eccentric energy was absorbed. Therefore, the efficiency of gait with hard boots might be decreased and joint loading at the knee might be increased, which might cause early fatigue of knee muscles during walking or hiking. The results of this study suggest that stiffness and blocking of joint motion at the ankle should not be equated with safety. A trade-off between lateral stiffness and free natural motion of ! the ankle joint complex might be preferable. - Internal forces and moments in the femur of the rat during gait
Wehner T Wolfram U Henzler T Niemeyer F Claes L Simon U - J Biomech 43(13):2473-2479 (2010)
The rat is of increasing importance for experimental studies on fracture healing. The healing outcome of long bone fractures is strongly influenced by mechanical factors, such as the interfragmentary movement. This movement depends on the stability of the fracture fixation and the musculoskeletal loads. However, little is known about these loads in rats. The musculoskeletal loads during gait were estimated using an inverse-dynamic musculoskeletal model of the right hindlimb of the rat. This model was based on a micro-CT scan of the lower extremities and an anatomical study using 15 rat cadavers. Kinematics were reconstructed from X-ray movies, taken simultaneously from two perpendicular directions during a gait cycle. The ground reaction forces were taken from the literature. The muscle forces were calculated using an optimization procedure. The internal forces and moments varied over the gait cycle and along the femoral axis. The greatest internal force (up to 7 times bodyweight) acted in the longitudinal direction. The greatest internal moment (up to 13.8 bodyweight times millimeter) acted in the sagittal plane of the femur. The validity of the model was corroborated by comparing the estimated strains caused by the calculated loads on the surface of the femoral mid-shaft with those from the literature. Knowledge of the internal loads in the femur of the rat allows adjustment of the biomechanical properties of fixation devices in fracture healing studies to the desired interfragmentary movement. - Scaffold microarchitecture determines internal bone directional growth structure: A numerical study
Sanz-Herrera JA Doblaré M García-Aznar JM - J Biomech 43(13):2480-2486 (2010)
A number of successful results have been reported in bone tissue engineering, although the routine clinical practice has not been reached so far. One of the reasons is the poor understanding of the role of each scaffold design parameter in its functional performance, which yields an uncertain outcome of each clinical application. Specifically, the role of internal scaffold microarchitectural shape on the regeneration rate and distribution of newly formed bone is still unknown. This work is focused on the in-silico determination of the role of scaffold microstructural anisotropy in bone tissue regeneration. A multiscale approach of the problem is established distinguishing between macroscopic region domain (bone organ and scaffold) and microscopic domain (scaffold microstructure). Results show that, once the scaffold microstructure is properly interconnected and the porosity is sufficiently high, similar rates of bone regeneration are found. However, the main conclusion! of the work is that initial scaffold microstructural anisotropy has important consequences since it determines the spatial distribution of the newly formed tissue. - Shoulder model validation and joint contact forces during wheelchair activities
- J Biomech 43(13):2487-2492 (2010)
Chronic shoulder impingement is a common problem for manual wheelchair users. The loading associated with performing manual wheelchair activities of daily living is substantial and often at a high frequency. Musculoskeletal modeling and optimization techniques can be used to estimate the joint contact forces occurring at the shoulder to assess the soft tissue loading during an activity and to possibly identify activities and strategies that place manual wheelchair users at risk for shoulder injuries. The purpose of this study was to validate an upper extremity musculoskeletal model and apply the model to wheelchair activities for analysis of the estimated joint contact forces. Upper extremity kinematics and handrim wheelchair kinetics were measured over three conditions: level propulsion, ramp propulsion, and a weight relief lift. The experimental data were used as input to a subject-specific musculoskeletal model utilizing optimization to predict joint contact forces ! of the shoulder during all conditions. The model was validated using a mean absolute error calculation. Model results confirmed that ramp propulsion and weight relief lifts place the shoulder under significantly higher joint contact loading than level propulsion. In addition, they exhibit large superior contact forces that could contribute to impingement. This study highlights the potential impingement risk associated with both the ramp and weight relief lift activities. Level propulsion was shown to have a low relative risk of causing injury, but with consideration of the frequency with which propulsion is performed, this observation is not conclusive. - Glenohumeral contact forces in reversed anatomy shoulder replacement
Masjedi M Johnson GR - J Biomech 43(13):2493-2500 (2010)
A major requirement to design an implant is to develop our understanding of the applied internal forces during everyday activities. In the absence of any basic apparatus for measuring forces directly, it is essential to rely on modelling. The major aim of this study was therefore to understand the biomechanical function of subjects with the reversed anatomy Bayley−Walker prosthesis, using an inverse dynamic shoulder model. In this context, the muscle and joint forces of 12 Bayley–Walker subjects were compared to those of 12 normal subjects during 12 activities of daily living. Maximum glenohumeral contact forces for normal and Bayley–Walker subjects were found to be 77% (±15) and 137% (±21) body weight for lifting a 2 kg shopping bag, and the least forces 29% (±4) and 67% (±8) body weight for reaching to opposite axilla, respectively. For normal subjects, middle deltoid, supraspinatus and infraspinatus were found to be the most active muscles across the subjects and tasks. On the other hand, for implanted subjects with a lack of rotator cuff muscles, the middle deltoid and coracobrachialis muscles were found to be the most active. The biomechanical model can therefore be used in order to gain knowledge about the pathology as well as possible post surgical rehab for subjects with reversed shoulder replacement. - Differential regulation of immature articular cartilage compressive moduli and Poisson's ratios by in vitro stimulation with IGF-1 and TGF-β1
Williams GM Dills KJ Flores CR Stender ME Stewart KM Nelson LM Chen AC Masuda K Hazelwood SJ Klisch SM Sah RL - J Biomech 43(13):2501-2507 (2010)
Mechanisms of articular cartilage growth and maturation have been elucidated by studying composition-function dynamics during in vivo development and in vitro culture with stimuli such as insulin-like growth factor-1 (IGF-1) and transforming growth factor-beta 1 (TGF-β1). This study tested the hypothesis that IGF-1 and TGF-β1 regulate immature cartilage compressive moduli and Poisson's ratios in a manner consistent with known effects on tensile properties. Bovine calf articular cartilage from superficial-articular (S) and middle-growth (M) regions were analyzed fresh or following culture in medium with IGF-1 or TGF-β1. Mechanical properties in confined (CC) and unconfined (UCC) compression, cartilage matrix composition, and explant size were assessed. Culture with IGF-1 resulted in softening in CC and UCC, increased Poisson's ratios, substantially increased tissue volume, and accumulation of glycosaminoglycan (GAG) and collagen (COL). Culture with TGF-β1 promot! ed maturational changes in the S layer, including stiffening in CC and UCC and increased concentrations of GAG, COL, and pyridinoline crosslinks (PYR), but little growth. Culture of M layer explants with TGF-β1 was nearly homeostatic. Across treatment groups, compressive moduli in CC and UCC were positively related to GAG, COL, and PYR concentrations, while Poisson's ratios were negatively related to concentrations of these matrix components. Thus, IGF-1 and TGF-β1 differentially regulate the compressive mechanical properties and size of immature articular cartilage in vitro. Prescribing tissue growth, maturation, or homeostasis by controlling the in vitro biochemical environment with such growth factors may have applications in cartilage repair and tissue engineering. - Upper limb joint kinetics during manual wheelchair propulsion in patients with different levels of spinal cord injury
Gil-Agudo A Del Ama-Espinosa A Pérez-Rizo E Pérez-Nombela S Pablo Rodríguez-Rodríguez L - J Biomech 43(13):2508-2515 (2010)
The purpose of this study was to compare the forces and moments of the whole upper limb, analyzing forces and moments at the shoulder, elbow and wrist joints simultaneously during manual wheelchair propulsion of persons with different levels of spinal cord injury (SCI) on a treadmill. Fifty-one people participated in this study and were grouped by their level of SCI: C6 tetraplegia (G1), C7 tetraplegia (G2), high paraplegia (G3), and low paraplegia (G4). An inverse dynamic model was defined to compute net joint forces and moments from segment kinematics, the forces acting on the pushrim, and subject anthropometrics. Right side, upper limb kinematic data were collected with four camcorders (Kinescan–IBV). Kinetic data were recorded by replacing the wheels with SmartWheels (Three Rivers Holdings, LLC). All participants propelled the wheelchair at 3 km/h for 1 min. The most noteworthy findings in both our tetraplegic groups in relation to paraplegic groups were increase! d superior joint forces in the shoulder (G1 and G2 vs G3 p<0.001; G1 and G2 vs G4 p<0.01), elbow (G1 vs G3 p<0.001; G1 vs G4 p<0.05) and wrist (G1 vs G4 p<0.001), an increased adduction moment in the shoulder (G1 vs G3 p<0.001; G1 vs G4 p<0.01; G2 vs G3 and G4 p<0.05) and the constancy of the moments of force of the wrist the fact that they reached their lowest values in the tetraplegic groups. This pattern may increase the risk of developing upper limb overuse injuries in tetraplegic subjects. - Low oxygen tension is a more potent promoter of chondrogenic differentiation than dynamic compression
Meyer EG Buckley CT Thorpe SD Kelly DJ - J Biomech 43(13):2516-2523 (2010)
During fracture healing and microfracture treatment of cartilage defects mesenchymal stem cells (MSCs) infiltrate the wound site, proliferate extensively and differentiate along a cartilaginous or an osteogenic lineage in response to local environmental cues. MSCs may be able to directly sense their mechanical environment or alternatively, the mechanical environment could act indirectly to regulate MSC differentiation by inhibiting angiogenesis and diminishing the supply of oxygen and other regulatory factors. Dynamic compression has been shown to regulate chondrogenesis of MSCs. In addition, previous studies have shown that a low oxygen environment promotes in vitro chondrogenesis of MSCs. The hypothesis of this study is that a low oxygen environment is a more potent promoter of chondrogenic differentiation of MSCs embedded in agarose hydrogels compared to dynamic compression. In MSC-seeded constructs supplemented with TGF-β3, GAG and collagen accumulation was higher! in low oxygen conditions compared to normoxia. For normoxic and low oxygen culture GAG accumulation within the agarose hydrogel was inhomogeneous, with low levels of GAG measured in the annulus of constructs maintained in normoxic conditions. Dynamic compression did not significantly increase GAG or collagen accumulation in normoxia. However under low oxygen conditions, dynamic compression reduced GAG accumulation compared to free-swelling controls, but remained higher than comparable constructs maintained in normoxic conditions. This study demonstrates that continuous exposure to low oxygen tension is a more potent pro-chondrogenic stimulus than 1 h/day of dynamic compression for porcine MSCs embedded in agarose hydrogels. - Stabilizing to disruptive transition of focal adhesion response to mechanical forces
Kong D Ji B Dai L - J Biomech 43(13):2524-2529 (2010)
Strong mechanical forces can, obviously, disrupt cell–cell and cell–matrix adhesions, e.g., cyclic uniaxial stretch induces instability of cell adhesion, which then causes the reorientation of cells away from the stretching direction. However, recent experiments also demonstrated the existence of force dependent adhesion growth (rather than dissociation). To provide a quantitative explanation for the two seemingly contradictory phenomena, a microscopic model that includes both integrin–integrin interaction and integrin–ligand interaction is developed at molecular level by treating the focal adhesion as an adhesion cluster. The integrin clustering dynamics and integrin–ligand binding dynamics are then simulated within one unified theoretical frame with Monte Carlo simulation. We find that the focal adhesion will grow when the traction force is higher than a relative small threshold value, and the growth is dominated by the reduction of local chemical potential! energy by the traction force. In contrast, the focal adhesion will rupture when the traction force exceeds a second threshold value, and the rupture is dominated by the breaking of integrin–ligand bonds. Consistent with the experiments, these results suggest a force map for various responses of cell adhesion to different scales of mechanical force. - Advanced human carotid plaque progression correlates positively with flow shear stress using follow-up scan data: An in vivo MRI multi-patient 3D FSI study
Yang C Canton G Yuan C Ferguson M Hatsukami TS Tang D - J Biomech 43(13):2530-2538 (2010)
Although it has been well-accepted that atherosclerosis initiation and early progression correlate negatively with flow wall shear stresses (FSS), increasing evidence suggests mechanisms governing advanced plaque progression are not well understood. Fourteen patients were scanned 2–4 times at 18 month intervals using a histologically validated multi-contrast magnetic resonance imaging (MRI) protocol to acquire carotid plaque progression data. Thirty-two scan pairs (baseline and follow-up scans) were formed with slices matched for model construction and analysis. 3D fluid–structure interaction (FSI) models were constructed and plaque wall stress (PWS) and flow shear stress (FSS) were obtained from all matching lumen data points (400–1000 per plaque; 100 points per matched slice) to quantify correlations with plaque progression measured by vessel wall thickness increase (WTI). Using FSS and PWS data from follow-up scan, 21 out of 32 scan pairs showed a significant ! positive correlation between WTI and FSS (positive/negative/no significance ratio=21/8/3), and 26 out of 32 scan pairs showed a significant negative correlation between WTI and PWS (positive/negative/no significance ratio=2/26/4). The mean FSS value of lipid core nodes (n=5294) from all 47 plaque models was 63.5 dyn/cm2, which was 45% higher than that from all normal vessel nodes (n=27553, p<0.00001). The results from this intensive FSI study indicate that flow shear stress from follow-up scan correlates positively with advanced plaque progression which is different from what has been observed in plaque initiation and early-stage progression. It should be noted that the correlation results do not automatically lead to any causality conclusions. - In vivo estimation of the short-range stiffness of cross-bridges from joint rotation
van Eesbeek S de Groot JH van der Helm FC de Vlugt E - J Biomech 43(13):2539-2547 (2010)
Short-range stiffness (SRS) is a mechanical property of muscles that is characterized by a disproportionally high stiffness within a short length range during both lengthening and shortening movements. SRS is attributed to the cross-bridges and is beneficial for stabilizing a joint during, e.g., postural conditions. Thus far, SRS has been estimated mainly on isolated mammalian muscles. In this study we presented a method to estimate SRS in vivo in the human wrist joint. SRS was estimated at joint level in the angular domain (N m/rad) for both flexion and extension rotations of the human wrist in nine healthy subjects. Wrist rotations of 0.15 rad at 3 rad/s were imposed at eight levels of voluntary contraction ranging from 0 to 2.1 N m by means of a single axis manipulator. Flexion and extension SRS of the wrist joint was estimated consistently and accurately using a dynamic nonlinear model that was fitted onto the recorded wrist torque. SRS increased monotonically with torque in a way consistent with previous studies on isolated muscles. It is concluded that in vivo measurement of joint SRS represents the population of coupled cross-bridges in wrist flexor and extensor muscles. In its current form, the presented technique can be used for clinical applications in many neurological and muscular diseases where altered joint torque and (dissociated) joint stiffness are important clinical parameters. - Control of center of mass motion state through cuing and decoupling of spontaneous gait parameters in level walking
Espy DD Yang F Pai YC - J Biomech 43(13):2548-2553 (2010)
Can the center of mass (COM) motion state, i.e., its position and velocity relative to the base of support (BOS), which dictate gait stability, be predictably controlled by the global gait parameters of step length and gait speed, or by extension, cadence? The precise relationships among step length and gait speed, and the COM motion state are unknown, partially due to the interdependence between step length and gait speed and the difficulty in independent control of both parameters during spontaneous level walking. The purposes of this study were to utilize simultaneous audio-visual cuing to independently manipulate step length and gait speed, and to determine the extent to which the COM position and velocity can be subsequently controlled. Fifty-six young adults were trained at one of the three gait patterns in which both the step length and gait speed were targeted simultaneously. The results showed that the cuing could successfully "decouple" gait speed from st! ep length. Although this approach did yield reliable control of the COM velocity through manipulation of gait speed (R2=0.97), the manipulation of step length yielded less precise control of COM position (R2=0.60). This latter control appears to require manipulation of an additional degree-of-freedom at the local segment level, such that the inclusion of trunk inclination with step length improved the prediction of COM position (R2=0.80). - Variability in kinematic coupling assessed by vector coding and continuous relative phase
Miller RH Chang R Baird JL Van Emmerik RE Hamill J - J Biomech 43(13):2554-2560 (2010)
Variability in the spatio-temporal coordination of human movement kinematics is often assessed by vector coding and continuous relative phase (CRP). To facilitate appropriate comparisons between the findings of studies that have used different techniques to assess variability, the purposes of this study were: (1) to determine if both vector coding and CRP behave according to dynamical systems theories on variability and state space transitions; and (2) to determine if trends in coordination variability during movement are consistent when using either vector coding or CRP. We present both a theoretical case (the Lorenz Attractor) and two experimental cases (rearfoot–forefoot coupling during overground walking for 22 subjects; the effect of treadmill speed on thigh-leg coupling for five subjects). In the theoretical case, variability quantified by CRP agreed with dynamical systems theory on state space transitions more so than variability quantified by vector coding. I! n experimental cases, this distinction was less clear, although CRP appeared to be a more conservative metric for variability. The magnitudes (all p<0.001) and timings (all p<0.04) of peaks in variability during the stance phase of overground walking depended on whether vector coding or CRP was used for two couplings. Similar distinctions were observed for peaks during the stride cycle of treadmill locomotion (all effect sizes >2.8). However, changes in the average variability during the stride cycle as speed increased were consistent for both methods (all effect sizes <0.2). The results suggest that comparisons between the findings of studies that have quantified variability using CRP and those that have used vector coding should be made with caution. - In vivo cartilage contact strains in patients with lateral ankle instability
Bischof JE Spritzer CE Caputo AM Easley ME Deorio JK Nunley JA Defrate LE - J Biomech 43(13):2561-2566 (2010)
Damage to the anterior talofibular ligament (ATFL) and cacaneofibular ligament (CFL) during an ankle sprain may be linked to the development of osteoarthritis. Although altered tibiotalar kinematics have been demonstrated, the effects of lateral ankle instability (LAI) on in vivo cartilage strains have not been described. We hypothesized that peak cartilage strains increase, and the location is shifted in patients with ATFL injuries. We used 3-D MRI models and biplanar fluoroscopy to evaluate in vivo cartilage contact strains in seven patients with unilateral LAI. Subjects had chronic unilateral ATFL injury or combined ATFL and CFL injury, and were evaluated with increasing load while stepping onto a force plate. Peak cartilage strain and the location of the peak strain were measured using the contralateral normal ankle as a control. Ankles with LAI demonstrated significantly increased peak strain when compared with ATFL-intact controls. For example, at 100% body weigh! t, peak strain was 29±8% on the injured side compared to 21±5% on the intact side. The position of peak strain on the injured ankle also showed significant anterior translation and medial translation. At 100% body weight, the location of peak strain in the injured ankle translated anteriorly by 15.5±7.1 mm and medially by 12.9±4.3 mm relative to the intact ankle. These changes correspond to the region of clinically observed osteoarthritis. Chronic LAI, therefore, may contribute to the development of tibiotalar cartilage degeneration due to altered cartilage strains. - Effect of sulfated glycosaminoglycan digestion on the transverse permeability of medial collateral ligament
Henninger HB Underwood CJ Ateshian GA Weiss JA - J Biomech 43(13):2567-2573 (2010)
Dermatan and chondroitin sulfate glycosaminoglycans (GAGs) comprise over 90% of the GAG content in ligament. Studies of their mechanical contribution to soft tissues have reported conflicting results. Measuring the transient compressive response and biphasic material parameters of the tissue may elucidate the contributions of GAGs to the viscoelastic response to deformation. The hypotheses of the current study were that digestion of sulfated GAGs would decrease compressive stress and aggregate modulus while increasing the permeability of porcine medial collateral ligament (MCL). Confined compression stress relaxation experiments were carried out on porcine MCL and tissue treated with chondroitinase ABC (ChABC). Results were fit to a biphasic constitutive model to derive permeability and aggregate modulus. Bovine articular cartilage was used as a benchmark tissue to verify that the apparatus provided reliable results. GAG digestion removed up to 88% of sulfated GAGs fro! m the ligament. Removal of sulfated GAGs increased the permeability of porcine MCL nearly 6-fold versus control tissues. Peak stress decreased significantly. Bovine articular cartilage exhibited the typical reduction of GAG content and resultant decreases in stress and modulus and increases in permeability with ChABC digestion. Given the relatively small amount of GAG in ligament (<1% of tissue dry weight) and the significant change in peak stress and permeability upon removal of GAGs, sulfated GAGs may play a significant role in maintaining the apposition of collagen fibrils in the transverse direction, thus supporting dynamic compressive loads experienced by the ligament during complex joint motion. - The effects of aponeurosis geometry on strain injury susceptibility explored with a 3D muscle model
Rehorn MR Blemker SS - J Biomech 43(13):2574-2581 (2010)
In the musculoskeletal system, some muscles are injured more frequently than others. For example, the biceps femoris longhead (BFLH) is the most commonly injured hamstring muscle. It is thought that acute injuries result from large strains within the muscle tissue, but the mechanism behind this type of strain injury is still poorly understood. The purpose of this study was to build computational models to analyze the stretch distributions within the BFLH muscle and to explore the effects of aponeurosis geometry on the magnitude and location of peak stretches within the model. We created a three-dimensional finite element (FE) model of the BFLH based on magnetic resonance (MR) images. We also created a series of simplified models with a similar geometry to the MR-based model. We analyzed the stretches predicted by the MR-based model during lengthening contractions to determine the region of peak local fiber stretch. The peak along-fiber stretch was 1.64 and was located ! adjacent to the proximal myotendinous junction (MTJ). In contrast, the average along-fiber stretch across all the muscle tissue was 0.95. By analyzing the simple models, we found that varying the dimensions of the aponeuroses (width, length, and thickness) had a substantial impact on the location and magnitude of peak stretches within the muscle. Specifically, the difference in widths between the proximal and distal aponeurosis in the BFLH contributed most to the location and magnitude of peak stretch, as decreasing the proximal aponeurosis width by 80% increased peak average stretches along the proximal MTJ by greater than 60% while slightly decreasing stretches along the distal MTJ. These results suggest that the aponeurosis morphology of the BFLH plays a significant role in determining stretch distributions throughout the muscle. Furthermore, this study introduces the new hypothesis that aponeurosis widths may be important in determining muscle injury susceptibility. - Tissue and cellular morphological changes in growth plate explants under compression
Amini S Veilleux D Villemure I - J Biomech 43(13):2582-2588 (2010)
The mechanisms by which mechanical loading may alter bone development within growth plates are still poorly understood. However, several growth plate cell or tissue morphological parameters are associated with both normal and mechanically modulated bone growth rates. The aim of this study was to quantify in situ the three-dimensional morphology of growth plate explants under compression at both cell and tissue levels. Growth plates were dissected from ulnae of immature swine and tested under 15% compressive strain. Confocal microscopy was used to image fluorescently labeled chondrocytes in the three growth plate zones before and after compression. Quantitative morphological analyses at both cell (volume, surface area, sphericity, minor/major radii) and tissue (cell/matrix volume ratio) levels were performed. Greater chondrocyte bulk strains (volume decrease normalized to the initial cell volume) were found in the proliferative (35.4%) and hypertrophic (41.7%) zones, wi! th lower chondrocyte bulk strains (24.7%) in the reserve zone. Following compression, the cell/matrix volume ratio decreased in the reserve and hypertrophic zones by 24.3% and 22.6%, respectively, whereas it increased by 35.9% in the proliferative zone. The 15% strain applied on growth plate explants revealed zone-dependent deformational states at both tissue and cell levels. Variations in the mechanical response of the chondrocytes from different zones could be related to significant inhomogeneities in growth plate zonal mechanical properties. The ability to obtain in situ cell morphometry and monitor the changes under compression will contribute to a better understanding of mechanisms through which abnormal growth can be triggered. - Linear vs. non-linear mapping of peak power using surface EMG features during dynamic fatiguing contractions
Gonzalez-Izal M Malanda A Rodríguez-Carreño I Navarro-Amézqueta I Gorostiaga EM Farina D Falla D Izquierdo M - J Biomech 43(13):2589-2594 (2010)
This study compares a non-linear (neural network) and a linear (linear regression) power mapping using a set of features of the surface electromyogram recorded from the vastus medialis and lateralis muscles. Fifteen healthy participants performed 5 sets of 10 repetitions leg press using the individual maximum load corresponding what they could perform 10 times (10RM) with 120 s of rest between them. The following sEMG variables were computed from each extension contraction and used as inputs to both approaches: mean average value (MAV), median frequency (Fmed), the spectral parameter proposed by Dimitrov (FInsm5), average (over the observation interval) of the instantaneous mean frequency obtained from a Choi–Williams distribution (MFM), and wavelet indices ratio between moments at different scales (WIRM1551, WIRM1M51, WIRM1522, WIRE51, and WIRW51). The non-linear mapping (neural network) provided higher correlation coefficients and signal-to-noise ratios values (alt! hough not significantly different) between the actual and the estimated changes of power compared to linear mapping (linear regression) using the sEMG variables alone and a combination of WIRW51 and MFM (obtained by a stepwise multiple linear regression). In conclusion, non-linear mapping of force loss during dynamic knee extension exercise showed higher signal-to-noise ratio and correlation coefficients between the actual and estimated power output compared to linear mapping. However, since no significant differences were observed between linear and non-linear approaches, both were equally valid to estimate changes in peak power during fatiguing repetitive leg extension exercise. - Knee contact force in subjects with symmetrical OA grades: Differences between OA severities
Richards C Higginson JS - J Biomech 43(13):2595-2600 (2010)
In using musculoskeletal models, researchers can calculate muscle forces, and subsequently joint contact forces, providing insight into joint loading and the progression of such diseases as osteoarthritis (OA). The purpose of this study was to estimate the knee contact force (KCF) in patients with varying degrees of OA severity using muscle forces and joint reaction forces derived from OpenSim. Walking data was obtained from healthy individuals (n=14) and those with moderate (n=10) and severe knee OA (n=2). For each subject, we generated 3D, muscle-actuated, forward dynamic simulations of the walking trials. Muscle forces that reproduced each subject's gait were calculated. KCFs were then calculated using the vector sum of the muscle forces and joint reaction forces along the longitudinal axis of the femur. Moderate OA subjects exhibited a similar KCF pattern to healthy subjects, with lower second peaks (p=0.021). Although subjects with severe OA had similar initial ! peak KCF to healthy and moderate OA subjects (more than 4 times BW), the pattern of the KCF was very different between groups. After an initial peak, subjects with severe OA continually unloaded the joint, whereas healthy and moderate OA subjects reloaded the knee during late stance. In subjects with symmetric OA grades, there appears to be differences in loading between OA severities. Similar initial peaks of KCF imply that reduction of peak KCF may not be a compensatory strategy for OA patients; however, reducing duration of high magnitude loads may be employed. - The convex wrapping algorithm: A method for identifying muscle paths using the underlying bone mesh
Desailly E Sardain P Khouri N Yepremian D Lacouture P - J Biomech 43(13):2601-2607 (2010)
Associating musculoskeletal models to motion analysis data enables the determination of the muscular lengths, lengthening rates and moment arms of the muscles during the studied movement. Therefore, those models must be anatomically personalized and able to identify realistic muscular paths. Different kinds of algorithms exist to achieve this last issue, such as the wired models and the finite elements ones. After having studied the advantages and drawbacks of each one, we present the convex wrapping algorithm. Its purpose is to identify the shortest path from the origin to the insertion of a muscle wrapping over the underlying skeleton mesh while respecting possible non-sliding constraints. After the presentation of the algorithm, the results obtained are compared to a classically used wrapping surface algorithm (obstacle set method) by measuring the length and moment arm of the semitendinosus muscle during an asymptomatic gait. The convex wrapping algorithm gives an ! efficient and realistic way of identifying the muscular paths with respect to the underlying bones mesh without the need to define simplified geometric forms. It also enables the identification of the centroid path of the muscles if their thickness evolution function is known. All this presents a particular interest when studying populations presenting noticeable bone deformations, such as those observed in cerebral palsy or rheumatic pathologies. - Symmetry indices based on accelerometric data in trotting horses
Halling Thomsen M Tolver Jensen A Sørensen H Lindegaard C Haubro Andersen P - J Biomech 43(13):2608-2612 (2010)
Detection and quantification of lameness in horses consists primarily of a subjective assessment, whereby both intra- and inter-observer disagreements exist, especially with low grade lameness. Therefore, clinically applicable methods are needed for reliable, objective assessments. The aim of this study was to describe three symmetry indices derived from a simple accelerometric method and investigate these in sound trotting horses. The indices describe the overall symmetry of the gait, the symmetry of loads placed on the limbs and the symmetry in timing between left and right steps. These symmetry indices were able to quantify the high degree of symmetry of the trot in sound horses that has been described in earlier studies using other gait analysis methods. Also, we have analysed the variances and have found high repeatability for all three indices. This provides a basis for future investigations of the symmetry indices and their potential for objective detection and quantification of lameness in horses. - Solutions for determining equibiaxial substrate strain for dynamic cell culture
Chiang MY Cheng T Pakstis L Dunkers J - J Biomech 43(13):2613-2617 (2010)
In this work, empirical and analytical solutions of equibiaxial strain on a flexible substrate are derived for a dynamic cell culture system. The empirical formula, which fulfills the mechanistic conditions of the culture system, is based on a regression analysis from finite element analyses for a substrate undergoing large strains (<15%). The analytical (closed-form) solution is derived from the superposition of two elastic responses induced in the equibiaxial strain culture system after applying pressure to a substrate undergoing small strains (microstrains). There is good agreement between the strain predicted from the solutions and from the direct measurement. Using material and geometric properties of the culture system, the solutions developed here are straightforward and can be used to circumvent experimental measurements or finite element analysis to establish substrate pressure–strain relationships. - A probabilistic finger biodynamic model better depicts the roles of the flexors during unloaded flexion
Li K Zhang X - J Biomech 43(13):2618-2624 (2010)
Previous deterministic finger biomechanical models predicted that the flexor digitorum superficialis (FDS) was silent and the flexor digitorum profundus (FDP) was the only active flexor during finger flexion. Experimental studies in vivo, however, recorded activities of both flexors. In this study, in an attempt to elucidate the roles of the flexors, a probabilistic biodynamic model of the index finger was constructed to estimate the muscle–tendon forces during an experimentally measured index finger flexion movement. A Monte-Carlo simulation was performed with four model parameters, including moment arms, physiological cross sectional areas (PCSA), passive torques, and anthropometric measures as independent random variables. The muscle-tendon forces at each time point were determined using a nonlinear optimization technique. The model predicted that both FDS and FDP contributed to sustaining the movement and the FDS was not necessarily silent. The two distinct force patterns observed in vivo in experimental studies were also corroborated by the simulation. These findings, contrary to previous deterministic models' predictions but in agreement with experimental measurements, explained the observed coactivation of FDS and FDP, and resolved the controversy regarding the roles of the flexors in finger movement dynamics. - A non-optimized follower load path may cause considerable intervertebral rotations
Dreischarf M Zander T Bergmann G Rohlmann A - J Biomech 43(13):2625-2628 (2010)
Osseoligamentous spinal specimens buckle under even a small vertical compressive force. To allow higher axial forces, a compressive follower load (FL) was suggested previously that approximates the curvature of the spine without inducing intervertebral rotation in both the frontal and the sagittal planes. In in vitro experiments and finite element analyses, the location of the FL path is subjected to estimation by the investigator. Such non-optimized FLs may induce bending and so far it is still unknown how this affects the results of the study and their comparability. A symmetrical finite element model of the lumbar spine was employed to simulate upright standing while applying a follower load. In analogy to in vitro experiments, the path of this FL was estimated seven times by different members of our institute's spine group. Additionally, an optimized FL path was determined and additional moments of ±7.5 Nm were applied to simulate flexion and extension. Application of the optimized 500 N compressive FL causes only a marginal alteration of the curvature (cardan angle L1–S1 in sagittal plane <0.25°). An individual estimation of the FL path, however, results in flexions of up to 10.0° or extensions of up to 12.3°. The resulting angles for the different non-optimized FL paths depend on the magnitude of the bending moment applied and whether a differential or an absolute measurement is taken. A preceding optimization of the location of the FL path would increase the comparability of different studies. - A visco-hyperelastic model for skeletal muscle tissue under high strain rates
Lu YT Zhu HX Richmond S Middleton J - J Biomech 43(13):2629-2632 (2010)
In this paper, a visco-hyperelastic skeletal muscle model is developed. The constitutive relation is based on the definition of a Helmholtz free energy function. It is assumed that the Helmholtz energy can be decomposed into volumetric and isochoric parts; furthermore, the isochoric energy can be decoupled into hyperelastic and viscous parts. The model developed involves 14 material parameters and its performance is evaluated by comparing the finite element simulation results with the published experimental studies on the New Zealand white rabbit tibialis anterior muscle. Results show that this model is able to describe the visco-hyperelastic behaviour of both passive and active skeletal muscle tissues under high strain rates (10/s and 25/s). - Assessment of a greater trochanter-based method of locating the hip joint center
Weinhandl JT O'Connor KM - J Biomech 43(13):2633-2636 (2010)
An alternative, yet unverified, predictive method that places the hip joint center (HJC) at one-quarter of the distance from the ipsolateral to the contralateral greater trochanter (GT method) is currently widely used in the biomechanics community. Therefore, the objective of this study was to confirm that this method is a viable option for estimating HJC coordinates. To accomplish this, HJC coordinates in the pelvic anatomical coordinate system were estimated via the GT method, a functional method, and the regression equations proposed by Bell et al. (1990). The HJC coordinated estimated by the functional method served as a baseline measurement. The results of this study demonstrate that all three methods evaluated offer repeatable estimates of HJC location. In comparison to the functional method, the GT method yielded a HJC estimate that was 7.6 mm medial, 12.2 mm posterior, and 4.8 mm proximal. On the other hand, the Bell regression equations estimated the HJC to be! 2.6 mm medial, 7.2 mm posterior, and 21.7 mm proximal relative to the functional method. Additionally, the total 3D difference between the GT and functional methods was 23.5 mm compared to the 30.8 mm difference between the Bell and functional methods. These results suggest that the GT method is a viable option for estimating HJC coordinates. - Full field strain measurements of collagenous tissue by tracking fiber alignment through vector correlation
- J Biomech 43(13):2637-2640 (2010)
Full field strain measurements of biological tissue during loading are often limited to the quantification of fiduciary marker displacements on the tissue surface. These marker measurements can lack the necessary spatial resolution to characterize non-uniform deformation and may not represent the deformation of the load-bearing collagen microstructure. To overcome these potential limitations, a method was developed to track the deformation of the collagen fiber microstructure in ligament tissue. Using quantitative polarized light imaging, fiber alignment maps incorporating both direction and alignment strength at each pixel were generated during facet capsular ligament loading. A grid of virtual markers was superimposed over the tissue in the alignment maps, and the maximization of a vector correlation calculation between fiber alignment maps was used to track marker displacement. Tracking error was quantified through comparisons to the displacements of excised ligamen! t tissue (n=3); separate studies applied uniaxial tension to isolated facet capsular ligament tissue (n=4) to evaluate tracking capabilities during large tissue deformations. The average difference between virtual marker and tissue displacements was 0.07±0.06 pixels. This error in marker location produced principal strain measurements of 1.2±1.6% when markers were spaced 4 pixels apart. During tensile tissue loading, substantial inhomogeneity was detected in the strain field using vector correlation tracking, and the location of maximum strain differed from that produced by standard tracking techniques using coarser meshes. These findings provide a method to directly measure fiber network strains using quantitative fiber alignment data, enabling a better understanding of structure–function relationships in tissues at different length scales. - Automated image analysis method for quantifying damage accumulation in tendon
Sereysky JB Andarawis-Puri N Ros SJ Jepsen KJ Flatow EL - J Biomech 43(13):2641-2644 (2010)
Tendon pathology is frequently sub-clinical prior to frank rupture, denoting the need for non-destructive methods of assessing disease presence and progression. Despite the lack of clinical presentation, previous studies have observed that distinct changes to the tendon microstructure are present, and that such qualitative changes have a dose–response relationship with the level of damage accumulated. These initial findings suggest that there is value in investigating the physical nature of damage within tendon, not only to better understand the pathological process, but also to gain insight into reparative processes and develop treatments. However, a necessary first step towards carrying out these avenues of research is to develop diagnostic tools for quantitatively assessing the level of damage present. In this study, we established a dose–response relationship between a quantitative measure of structural damage and the level of global damage induced. Furthermore, we developed and validated an automated technique for quantifying matrix disorganization (damage), which correlates and co-localizes strongly with manual quantification. In combination, these findings allow us to measure the amount of damage accumulation of a region of tendon on a clinical scale, rapidly and accurately. - Beyond the Virtual Intracranial Stenting Challenge 2007: Non-Newtonian and flow pulsatility effects
Cavazzuti M Atherton M Collins M Barozzi G - J Biomech 43(13):2645-2647 (2010)
The Virtual Intracranial Stenting Challenge 2007 (VISC'07) is becoming a standard test case in computational minimally invasive cerebrovascular intervention. Following views expressed in the literature and consistent with the recommendations of a report, the effects of non-Newtonian viscosity and pulsatile flow are reported. Three models of stented cerebral aneurysms, originating from VISC'07 are meshed and the flow characteristics simulated using commercial computational fluid dynamics (CFD) software. We conclude that non-Newtonian and pulsatile effects are important to include in order to discriminate more effectively between stent designs. - All joint moments significantly contribute to trunk angular acceleration
Nott CR Zajac FE Neptune RR Kautz SA - J Biomech 43(13):2648-2652 (2010)
Computationally advanced biomechanical analyses of gait demonstrate the often counter-intuitive roles of joint moments on various aspects of gait such as propulsion, swing initiation, and balance. Each joint moment can produce linear and angular acceleration of all body segments (including those on which the moment does not directly act) due to the dynamic coupling inherent in the interconnected musculoskeletal system. This study presents quantitative relationships between individual joint moments and trunk control with respect to balance during gait to show that the ankle, knee, and hip joint moments all affect the angular acceleration of the trunk. We show that trunk angular acceleration is affected by all joints in the leg with varying degrees of dependence during the gait cycle. Furthermore, it is shown that inter-planar coupling exists and a two-dimensional analysis of trunk balance neglects important out-of-plane joint moments that affect trunk angular accelerati! on. - Letter to the Editor regarding "All joint moments significantly contribute to trunk angular acceleration"
Winter DA - J Biomech 43(13):2653 (2010)
- Assumption of a 'gravity only region' for gravity correction of passive joint moment data may be problematic
Clarke EC Martin JH Herbert R - J Biomech 43(13):2653-2655 (2010)
- Response to "Assumption of a 'gravity only region' for gravity correction of passive joint moment data may be problematic"
Anderson DE Madigan ML Nussbaum MA - J Biomech 43(13):2655-2656 (2010)
- Authors' response to "Comments on 'Validation of a musculoskeletal model of wheelchair propulsion and its application to minimizing shoulder joint forces'"
Dubowsky SR - J Biomech 43(13):2656 (2010)
No comments:
Post a Comment