Search

Purpose: Muscles are the primary component responsible for the locomotion and change of posture of the human body. The physiologic basis of muscle force production and movement is determined by the muscle architecture (maximum muscle force, [Formula: see text], optimal muscle fiber length, [Formula: see text], tendon slack length, [Formula: see text], and pennation angle at optimal muscle fiber length, [Formula: see text]). The pennation angle is related to the maximum force production and to the range of motion. The aim of this study was to investigate a computational approach to calculate subject-specific pennation angle from magnetic resonance images (MRI)-based 3D anatomical model and to determine the impact of this approach on the motion analysis with personalized musculoskeletal models., Methods: A 3D method that calculates the pennation angle using MRI was developed. The fiber orientations were automatically computed, while the muscle line of action was determined using approaches based on anatomical landmarks and on centroids of image segmentation. Three healthy male volunteers were recruited for MRI scanning and motion capture acquisition. This work evaluates the effect of subject-specific pennation angle as musculoskeletal parameter in the lower limb, focusing on the quadriceps group. A comparison was made for assessing the contribution of personalized models on motion analysis. Gait and deep squat were analyzed using neuromuscular simulations (OpenSim)., Results: The results showed variation of the pennation angle between the generic and subject-specific models, demonstrating important interindividual differences, especially for the vastus intermedius and vastus medialis muscles. The pennation angle variation between personalized and generic musculoskeletal models generated significant variation in muscle moments and forces during dynamic motion analysis., Conclusions: A MRI-based approach to define subject-specific pennation angle was proposed and evaluated in motion analysis models. The significant differences obtained for the moments and muscle forces in quadriceps muscles indicate that a personalized approach in modeling the pennation angle can provide more individual details when investigating motion behaviors in specific subjects.

The growing numbers of short stem hip implants have redefined total hip arthroplasty with new stem geometries and possible functional differences. Several systematic reviews have reported good clinical results with this new class of stems, although kinematic alterations are still unclear in many aspects. The good clinical results obtained at the authors' institution led to the current study. The authors hypothesized that the geometric alignment of the prosthetic components may be closer to the anatomy of the healthy hip joint, thus leading to better function and clinical satisfaction. An examination via finite element analysis was chosen to model the hip joint and virtually implant a short and a standard straight stem. Findings indicated that anchoring of the short stem allowed favorable positioning in the proximal femur, with the femoral head already in the center of the cup. This positioning was not possible for the straight stem, which required further reduction of the femur by a significant translation into the cup, leading to abnormal soft-tissue balancing. The results from the simulation showed an absolute average deviation of ligamentous fiber strains of 6% for the short stem in 30° of flexion and extension versus 29% and 36% for the standard straight stem in 30° of flexion and extension, respectively. A femoral neck guided orientation of the short stem implant seems to allow a more anatomical reconstruction and thus a more balanced hip in terms of the modeled soft tissues. In contrast, the straight stem alters the head position and induces nonphysiological capsular strains., (Copyright 2015, SLACK Incorporated.)

Skeletal muscles are characterized by a large diversity in anatomical architecture and function. Muscle force and contraction are generated by contractile fiber cells grouped in fascicle bundles, which transmit the mechanical action between origin and insertion attachments of the muscle. Therefore, an adequate representation of fascicle arrangements in computational models of skeletal muscles is important, especially when investigating three-dimensional muscle deformations in finite element models. However, obtaining high resolution in vivo measurements of fascicle arrangements in skeletal muscles is currently still challenging. This motivated the development of methods in previous studies to generate numerical representations of fascicle trajectories using interpolation templates. Here, we present an alternative approach based on the hypothesis of a rotation and divergence free (Laplacian) vector field behavior which reflects observed physical characteristics of fascicle trajectories. To obtain this representation, the Laplace equation was solved in anatomical reconstructions of skeletal muscle shapes based on medical images using a uniform flux boundary condition on the attachment areas. Fascicle tracts were generated through a robust flux based tracing algorithm. The concept of this approach was demonstrated in two-dimensional synthetic examples of typical skeletal muscle architectures. A detailed evaluation was performed in an example of the anatomical human tibialis anterior muscle which showed an overall agreement with measurements from the literature. The utility and capability of the proposed method was further demonstrated in other anatomical examples of human skeletal muscles with a wide range of muscle shapes and attachment morphologies.

Ultrasound assessment of myocardial strain can provide valuable information on regional cardiac function. However, Doppler-based methods often used in practice for strain estimation suffer from angle dependency. In this study, a partial solution to that fundamental limitation is presented. We have previously reported using simulated data sets that spatial compounding of axial velocities obtained at three steering angles can theoretically outperform 2-D speckle tracking for 2-D strain estimation in the mouse heart. In this study, the feasibility of the method was analyzed in vivo using spatial compounding of Doppler velocities on six mice with myocardial infarction and five controls, and results were compared with those of tagged microscopic magnetic resonance imaging (μMRI). Circumferential estimates quantified by means of both ultrasound and μMRI could detect regional dysfunction. Between echocardiography and μMRI, a good regression coefficient was obtained for circumferential strain estimates (r = 0.69), whereas radial strain estimates correlated only moderately (r = 0.37). A second echocardiography was performed after μMRI to test the reproducibility of the compounding method. This yielded a higher correlation coefficient for the circumferential component than for the radial component (r = 0.74 circumferentially, r = 0.49 radially)., (2013 World Federation for Ultrasound in Medicine & Biology. Published by Elsevier Inc. All rights reserved)

Estimating myocardial strain in the mouse with clinical equipment remains difficult due to the high heart rate and the small size of the mouse heart. Measuring the strain component perpendicular to the ultrasound beam is especially challenging because of the lack of phase information in that direction and the large speckle width compared to the wall thickness. In this study, the performance of a Synthetic Lateral Phase (SLP) approach was contrasted to a standard and a regularized 2D Speckle Tracking (2D ST) algorithm using simulated data sets. SLP yielded higher rms errors for the lateral strain estimates than the regularized 2D ST (Lateral rms error: 0.087±0.012 vs. 0.052±0.010; p<0.05). No significant difference was found between the standard 2D ST and SLP. For the axial strain estimates, SLP produced higher rms errors than the standard 2D ST (Axial rms error: 0.063±0.012 vs. 0.040±0.008; p<0.05). 2D ST combined with geometric regularization showed thus to be the most accurate method., (Copyright © 2012 Elsevier B.V. All rights reserved.)

When imaging the heart, good temporal resolution is beneficial for capturing the information of short-lived cardiac phases (in particular, the isovolumetric phases). To increase the frame rate, parallel beamforming is a commonly used technique for fast cardiac imaging. Conventionally, a 4 multiple-line-acquisition (4MLA) system increases the frame rate by a factor of 4, making use of a broadened transmit beam to reduce block-like artifacts. As an alternative, it has been proposed to transmit an unfocused beam (i.e., plane wave or diverging wave) for which a large number of parallel receive beams (i.e., 16) can be formed for each transmit. However, to keep the spatial resolution acceptable in these approaches, spatial compounding of overlapping successive transmits is required. As a result, the effective gain in frame rate is similar to that of a 4MLA system. To date, it remains unclear how conventional 4MLA compares to plane-wave or diverging-wave imaging when operating at similar frame rate. The goal of this study was therefore to directly contrast the performance of these beamforming methods by computer simulation. In this study, the performance of 4 different imaging systems was investigated by quantitatively evaluating the characteristics of their beam profiles. The results showed that the conventional 4MLA and plane wave imaging were very competitive imaging strategies when operating at a similar frame rate. 4MLA performed better in the near field (i.e., 10 to 50 mm), whereas plane-wave imaging had better beam profiles in the far field (i.e., 50 to 90 mm). Although diverging-wave imaging had the poorest performance in the present study, it could potentially be improved by optimizing the settings.

Automatic quantification of regional left ventricular deformation in volumetric ultrasound data remains challenging. Many methods have been proposed to extract myocardial motion, including techniques using block matching, phase-based correlation, differential optical flow methods, and image registration. Our lab previously presented an approach based on elastic registration of subsequent volumes using a B-spline representation of the underlying transformation field. Encouraging results were obtained for the assessment of global left ventricular function, but a thorough validation on a regional level was still lacking. For this purpose, univentricular thick-walled cardiac phantoms were deformed in an experimental setup to locally assess strain accuracy against sonomicrometry as a reference method and to assess whether regions containing stiff inclusions could be detected. Our method showed good correlations against sonomicrometry: r(2) was 0.96, 0.92, and 0.84 for the radial (ε(RR)), longitudinal (ε(LL)), and circumferential (ε(CC)) strain, respectively. Absolute strain errors and strain drift were low for ε(LL) (absolute mean error: 2.42%, drift: -1.05%) and ε(CC) (error: 1.79%, drift: -1.33%) and slightly higher for ε(RR) (error: 3.37%, drift: 3.05%). The discriminative power of our methodology was adequate to resolve full transmural inclusions down to 17 mm in diameter, although the inclusion-to-surrounding tissue stiffness ratio was required to be at least 5:2 (absolute difference of 39.42 kPa). When the inclusion-to-surrounding tissue stiffness ratio was lowered to approximately 2:1 (absolute difference of 22.63 kPa), only larger inclusions down to 27 mm in diameter could still be identified. Radial strain was found not to be reliable in identifying dysfunctional regions.

Left-ventricular remodeling is considered to be an important mechanism of disease progression leading to mechanical dysfunction of the heart. However, the interaction between the physiological changes in the remodeling process and the associated mechanical dysfunction is still poorly understood. Clinically, it has been observed that the left ventricle often undergoes large shape changes, but the importance of left-ventricular shape as a contributing factor to alterations in mechanical function has not been clearly determined. Therefore, the interaction between left-ventricular shape and systolic mechanical function was examined in a computational finite-element study. Hereto, finite-element models were constructed with varying shapes, ranging from an elongated ellipsoid to a sphere. A realistic transmural gradient in fiber orientation was considered. The passive myocardium was described by an incompressible hyperelastic material law with transverse isotropic symmetry. Activation was governed by the eikonal-diffusion equation. Contraction was incorporated using a Hill model. For each shape, simulations were performed in which passive filling was followed by isovolumic contraction and ejection. It was found that the intramyocardial distributions of fiber stress, strain, and stroke work density were shape dependent. Ejection performance was reduced with increasing sphericity, which was regionally related to a reduction in the active fiber stress development, fiber shortening, and stroke work in the midwall and subepicardial region at the midheight level in the left-ventricular wall. Based on these results, we conclude that a significant interaction exists between left-ventricular shape and regional myofiber mechanics, but the importance for left-ventricular remodeling requires further investigation.

Fish swimming is probably one of the most studied and best understood locomotor behaviors in vertebrates. However, many fish also actively exploit sediments. Because of their elongate body shape, anguilliform fishes are not only efficient swimmers but also very maneuverable. Consequently, many species live in complexly structured environments near the bottom and many are known to burrow into the sediment. To better understand burrowing and subsurface locomotion in anguilliform fish we provide descriptive kinematic data on subsurface locomotion in a burrowing eel (Pisodonophis boro) using videofluoroscopy. We also measured the maximal forces that can be exerted by this species during head-first and tail-first burrowing, and explored the implications of head-first burrowing on mechanical stress distribution in the skull. Our data show that P. boro uses lateral undulation to penetrate and move in sandy sediments under water. The kinematics of subsurface locomotion are different from those observed during swimming and are characterized by a very high slip factor. These observations differ considerably from recently published data in terrestrial sand-swimming lizards, and suggest that the sediment behaves like a solid rather than a frictional fluid. Finally, our finite element models show that the cranial shape and structure in the head-first burrowing P. boro is mechanically more suited for head-first burrowing than that of an obligate tail-first burrowing species, Heteroconger hassi.

Anguilliform or eel-like fishes are typically bottom dwellers, some of which are specialized burrowers. Although specializations for burrowing are predicted to affect the kinematics of swimming, it remains unknown to what extent this is actually the case. Here we examine swimming kinematics and efficiency of two burrowing anguilliform species, Pisodonophis boro and Heteroconger hassi, with different degrees of specialization for burrowing. Our data suggest that differences in the swimming kinematics may indeed be related to the differences in burrowing specialization and style between both species. The resemblance between the swimming kinematics of P. boro and previously published data for Anguilla anguilla and Anguilla rostrata may be linked with the relatively limited burrowing specialization of P. boro and suggests an overall stereotypy in anguilliform forward-swimming patterns. The body of H. hassi, in contrast, is more specialized for tail-first burrowing and backward swimming bears a striking resemblance to the backward burrowing motions observed in this species. These motions differ significantly from backward swimming in Anguilla and in P. boro. The kinematics of forward swimming are, however, comparable across species. Thus, our data suggest that specializations for burrowing may affect swimming kinematics in anguilliform fishes, but also that forward swimming and burrowing are not necessarily incompatible. Future studies comparing the kinematics and mechanics of burrowing in these and other anguilliform fishes are needed to better understand how specializations for burrowing constrain backward swimming in H. hassi., (Copyright © 2011 Elsevier GmbH. All rights reserved.)

Myocardial strain quantification in the mouse based on 2-D speckle tracking using real-time ultrasound datasets is feasible but remains challenging. The major difficulty lies in the fact that the frame rate-to-heart rate ratio is relatively low, causing significant decorrelation between subsequent frames. In this setting, regularization is therefore particularly important to discard motion estimates that are improbable. Different regularization methods have been proposed, among which is a class of regularizers based on enforcing preset geometrical characteristics of the motion field. To date, these regularization methods have not been contrasted. The aim of this study was thus to compare the performance of different geometric regularizers in the setting of myocardial motion and strain estimation in murine echocardiography using simulated datasets. In normal models, restricting the spatial curvature of the motion fields resulted in worse radial strain estimates (mean root-mean-square [RMS] error increased from 0.06 to 0.09; p < 0.05), but better circumferential strain estimates (mean RMS error decreased from 0.035 to 0.01; p < 0.05). More accurate circumferential strain estimates were also obtained by convolving a Gaussian function with the lateral motion components (mean RMS error decreased to 0.015; p < 0.05). In infarcted models, no significant differences were found between regularized and nonregularized radial strains. However, for circumferential strain, the curvature method yielded better strain estimates in all regions (mean RMS error decreased from 0.043 to 0.015; p < 0.05), whereas the Gaussian method only improved strain assessment in the remote myocardium (mean RMS error decreased to 0.021; p < 0.05)., (Copyright 2010 World Federation for Ultrasound in Medicine & Biology. Published by Elsevier Inc. All rights reserved.)

Passive filling is a major determinant for the pump performance of the left ventricle and is determined by the filling pressure and the ventricular compliance. In the quantification of the passive mechanical behaviour of the left ventricle and its compliance, focus has been mainly on fiber orientation and constitutive parameters. Although it has been shown that the left-ventricular shape plays an important role in cardiac (patho-)physiology, the dependency on left-ventricular shape has never been studied in detail. Therefore, we have quantified the influence of left-ventricular shape on the overall compliance and the intramyocardial distribution of passive fiber stress and strain during the passive filling period. Hereto, fiber stress and strain were calculated in a finite element analysis of passive inflation of left ventricles with different shapes, ranging from an elongated ellipsoid to a sphere, but keeping the initial cavity volume constant. For each shape, the wall volume was varied to obtain ventricles with different wall thickness. The passive myocardium was described by an incompressible hyperelastic material law with transverse isotropic symmetry along the muscle fiber directions. A realistic transmural distribution in fiber orientation was assumed. We found that compliance was not altered substantially, but the transmural distribution of both passive fiber stress and strain was highly dependent on regional wall curvature and thickness. A low curvature wall was characterized by a maximum in the transmural fiber stress and strain in the mid-wall region, while a steep subendocardial transmural gradient was present in a high curvature wall. The transmural fiber stress and strain gradients in a low and high curvature wall were, respectively, flattened and steepened by an increase in wall thickness., (Copyright (c) 2010 Elsevier Ltd. All rights reserved.)

Left-ventricular shape is an important determinant of regional wall mechanics during passive filling. To examine the influence of left-ventricular shape for the ejection phase, the distribution of active fiber stress at the beginning of ejection was calculated in a finite element study. Hereto, finite element models were constructed with varying left-ventricular shapes, ranging from an elongated ellipsoid to a sphere, but keeping the initial cavity and wall volume constant. A realistic transmural gradient in fiber orientation was assumed. The passive myocardium was described by an incompressible hyperelastic material law with transverse isotropic symmetry along the muscle fiber directions. The activation of the left-ventricular wall was governed by the eikonal-diffusion equation. Active contraction was incorporated using a Hill-like model. For each left-ventricular shape, a simulation was performed in which passive filling was followed by isovolumic contraction. It was found that the transmural gradient of active fiber stress at the beginning of ejection steepens at the mid-height level when the left ventricle becomes more spherical, which was also obtained previously for end-diastolic passive fiber stress and strain.

Passive filling is a major determinant for the pump performance of the left ventricle and is determined by the filling pressure and the ventricular compliance. We quantified the influence of left-ventricular shape on the overall compliance and the distribution of passive fiber stress and strain during the filling period in normal myocardium. Hereto, fiber stress and strain were calculated in a finite element analysis during the inflation of left ventricles of different shape, ranging from an elongated ellipsoid to a sphere, but keeping the initial cavity and wall volume constant. The passive myocardium was described by an incompressible hyperelastic material law with transverse isotropic symmetry along the muscle fiber directions. A realistic transmural gradient in fiber orientation was assumed. While compliance was not altered, the transmural distribution of both passive fiber stress and strain was highly dependent on ventricular shape, where more spherical ventricles exhibited a higher subendocardial gradient in both quantities.

Current ultrasound methods for measuring myocardial strain are often limited to measurements in one or two dimensions. Cardiac motion and deformation however are truly 3-D. With the introduction of matrix transducer technology, 3-D ultrasound imaging of the heart has become feasible but suffers from low temporal and spatial resolution, making 3-D strain estimation challenging. In this paper, it is shown that automatic intensity-based spatio-temporal elastic registration of currently available 3-D volumetric ultrasound data sets can be used to measure the full 3-D strain tensor. The method was validated using simulated 3-D ultrasound data sets of the left ventricle (LV). Three types of data sets were simulated: a normal and symmetric LV with different heart rates, a more realistic asymmetric normal LV and an infarcted LV. The absolute error in the estimated displacement was between 0.47 +/-0.23 and 1.00 +/-0.59 mm, depending on heart rate and amount of background noise. The absolute error on the estimated strain was 9%-21% for the radial strain and 1%-4% for the longitudinal and circumferential strains. No large differences were found between the different types of data sets. The shape of the strain curves was estimated properly and the position of the infarcts could be identified correctly. Preliminary results on clinical data taken in vivo from three healthy volunteers and one patient with an apical aneurism confirmed these findings in a qualitative manner as the strain curves obtained with the proposed method have an amplitude and shape similar to what could be expected.

A local description of the shape of the left ventricle is relevant in assessing the process of adverse ventricular remodeling, associated with most cardiac pathologies, and in monitoring reverse remodeling by therapy. To quantify local shape of the left ventricle, one can calculate the curvature of its epicardial or endocardial surface. The 3D geometry of the heart and especially the ventricles, can typically be described using finite element meshes. From a mathematical point of view these meshes provide a local parametrization of the surface in the 3-dimensional space. We discuss the analytic derivation of the principle curvatures of the left-ventricular surfaces given their smooth finite-element meshes and apply this derivation to assess the regional shape of the normal porcine left ventricle.