Showing papers in "Journal of Biomechanical Engineering-transactions of The Asme in 1999"

The effects of the surface topography of micromachined titanium substrata on cell behavior in vitro and in vivo.

Abstract: Surface properties, including topography and chemistry, are of prime importance in establishing the response of tissues to biomaterials. Microfabrication techniques have enabled the production of precisely controlled surface topographies that have been used as substrata for cells in culture and on devices implanted in vivo. This article reviews aspects of cell behavior involved in tissue response to implants with an emphasis on the effects of topography. Microfabricated grooved surfaces produce orientation and directed locomotion of epithelial cells in vitro and can inhibit epithelial downgrowth on implants. The effects depend on the groove dimensions and they are modified by epithelial cell-cell interactions. Fibroblasts similarly exhibit contact guidance on grooved surfaces, but fibroblast shape in vitro differs markedly from that found in vivo. Surface topography is important in establishing tissue organization adjacent to implants, with smooth surfaces generally being associated with fibrous tissue encapsulation. Grooved topographies appear to have promise in reducing encapsulation in the short term, but additional studies employing three-dimensional reconstruction and diverse topographies are needed to understand better the process of connective-tissue organization adjacent to implants. Microfabricated surfaces can increase the frequency of mineralized bone-like tissue nodules adjacent to subcutaneously implanted surfaces in rats. Orientation of these nodules with grooves occurs both in culture and on implants. Detailed comparisons of cell behavior on micromachined substrata in vitro and in vivo are difficult because of the number and complexity of factors, such as population density and micromotion, that can differ between these conditions.

Effect of Impact Load on Articular Cartilage: Cell Metabolism and Viability, and Matrix Water Content

Abstract: Significant evidence exists that trauma to a joint produced by a single impact load below that which causes subchondral bone fracture can result in permanent damage to the cartilage matrix, including surface fissures, loss of proteoglycan, and cell death. Limited information exists, however, on the effect of a varying impact stress on chondrocyte biophysiology and matrix integrity. Based on our previous work, we hypothesized that a stress-dependent response exists for both the chondrocyte's metabolic activity and viability and the matrix's hydration. This hypothesis was tested by impacting bovine cartilage explants with nominal stresses ranging from 0.5 to 65 MPa and measuring proteoglycan biosynthesis, cell viability, and water content immediately after impaction and 24 hours later. We found that proteoglycan biosynthesis decreased and water content increased with increasing impact stress. However, there appeared to be a critical threshold stress (15-20 MPa) that caused cell death and apparent rupture of the collagen fiber matrix at the time of impaction. We concluded that the cell death and collagen rupture are responsible for the observed alterations in the tissue's metabolism and water content, respectively, although the exact mechanism causing this damage could not be determined.

A validated three-dimensional computational model of a human knee joint.

Abstract: This paper presents a three-dimensional finite element tibio-femoral joint model of a human knee that was validated using experimental data. The geometry of the joint model was obtained from magnetic resonance (MR) images of a cadaveric knee specimen. The same specimen was biomechanically tested using a robotic/universal force-moment sensor (UFS) system and knee kinematic data under anterior-posterior tibial loads (up to 100 N) were obtained. In the finite element model (FEM), cartilage was modeled as an elastic material, ligaments were represented as nonlinear elastic springs, and menisci were simulated by equivalent-resistance springs. Reference lengths (zero-load lengths) of the ligaments and stiffness of the meniscus springs were estimated using an optimization procedure that involved the minimization of the differences between the kinematics predicted by the model and those obtained experimentally. The joint kinematics and in-situ forces in the ligaments in response to axial tibial moments of up to 10 Nm were calculated using the model and were compared with published experimental data on knee specimens. It was also demonstrated that the equivalent-resistance springs representing the menisci are important for accurate calculation of knee kinematics. Thus, the methodology developed in this study can be a valuable tool for further analysis of knee joint function and could serve as a step toward the development of more advanced computational knee models.

Analysis of indentation: implications for measuring mechanical properties with atomic force microscopy.

Abstract: Indentation using the atomic force microscope (AFM) has potential to measure detailed micromechanical properties of soft biological samples. However, interpretation of the results is complicated by the tapered shape of the AFM probe tip, and its small size relative to the depth of indentation. Finite element models (FEMs) were used to examine effects of indentation depth, tip geometry, and material nonlinearity and heterogeneity on the finite indentation response. Widely applied infinitesimal strain models agreed with FEM results for linear elastic materials, but yielded substantial errors in the estimated properties for nonlinear elastic materials. By accounting for the indenter geometry to compute an apparent elastic modulus as a function of indentation depth, nonlinearity and heterogeneity of material properties may be identified. Furthermore, combined finite indentation and biaxial stretch may reveal the specific functional form of the constitutive law--a requirement for quantitative estimates of material constants to be extracted from AFM indentation data.

A fibril-network-reinforced biphasic model of cartilage in unconfined compression.

Abstract: Cartilage mechanical function relies on a composite structure of a collagen fibrillar network entrapping a proteoglycan matrix. Previous biphasic or poroelastic models of this tissue, which have approximated its composite structure using a homogeneous solid phase, have experienced difficulties in describing measured material responses. Progress to date in resolving these difficulties has demonstrated that a constitutive low that is successful for one test geometry (confined compression) is not necessarily successful for another (unconfined compression). In this study, we hypothesize that an alternative fibril-reinforced composite biphasic representation of cartilage can predict measured material responses and explore this hypothesis by developing and solving analytically a fibril-reinforced biphasic model for the case of uniaxial unconfined compression with frictionless compressing platens. The fibrils were considered to provide stiffness in tension only. The lateral stiffening provided by the fibril network dramatically increased the frequency dependence of disk rigidity in dynamic sinusoidal compression and the magnitude of the stress relaxation transient, in qualitative agreement with previously published data. Fitting newly obtained experimental stress relaxation data to the composite model allowed extraction of mechanical parameters from these tests, such as the rigidity of the fibril network, in addition to the elastic constants and the hydraulic permeability of the remaining matrix. Model calculations further highlight a potentially important difference between homogeneous and fibril-reinforced composite models. In the latter type of model, the stresses carried by different constituents can be dissimilar, even in sign (compression versus tension) even though strains can be identical. Such behavior, resulting only from a structurally physiological description, could have consequences in the efforts to understand the mechanical signals that determine cellular and extracellular biological responses to mechanical loads in cartilage.

Molding of Deep Polydimethylsiloxane Microstructures for Microfluidics and Biological Applications

Abstract: Here we demonstrate the microfabrication of deep (> 25 microns) polymeric microstructures created by replica-molding polydimethylsiloxane (PDMS) from microfabricated Si substrates. The use of PDMS structures in microfluidics and biological applications is discussed. We investigated the feasibility of two methods for the microfabrication of the Si molds: deep plasma etch of silicon-on-insulator (SOI) wafers and photolithographic patterning of a spin-coated photoplastic layer. Although the SOI wafers can be patterned at higher resolution, we found that the inexpensive photoplastic yields similar replication fidelity. The latter is mostly limited by the mechanical stability of the replicated PDMS structures. As an example, we demonstrate the selective delivery of different cell suspensions to specific locations of a tissue culture substrate resulting in micropatterns of attached cells.

Rapid, automated nucleic acid probe assays using silicon microstructures for nucleic acid concentration.

Abstract: A system for rapid point-of-use nucleic acid (NA) analysis based on PCR techniques is described. The extraction and concentration of DNA from test samples has been accomplished utilizing silicon fluidic microchips with high surface-area-to-volume ratios. Short (500 bp) and medium size (48,000 bp) DNA have been captured, washed, and eluted using the silicon dioxide surfaces of these chips. Chaotropic (GuHCl) salt solutions were used as binding agents. Wash and elution agents consisted of ethanol-based solutions and water, respectively. DNA quantities approaching 40 ng/cm2 of binding area were captured from input solutions in the 100-1000 ng/mL concentration range. For dilute samples of interest for pathogen detection, PCR and gel electrophoresis were used to demonstrate extraction efficiencies of about 50 percent, and concentration factors of about 10x using bacteriophage lambda DNA as the target. Rapid, multichannel PCR thermal cycling modules with integrated solid-state detection components have also been demonstrated. These results confirm the viability of utilizing these components as elements of a compact, disposable cartridge system for the detection of NA in applications such as clinical diagnostics, biowarfare agent detection, food quality control, and environmental monitoring.

Microcontact printing for precise control of nerve cell growth in culture.

Abstract: Microcontact printing, facilitated by silane linker chemistry and high-relief stamps, creates precise patterns of proteins, which in turn control growth of hippocampal neurons in culture. This additive, multi-mask technique permits several different molecules to be patterned on the same substrate. The covalent linker technology permits relatively long-term (two-week) compliance of neurons to the stamped pattern against a polyethylene glycol background. When polylysine was stamped adjacent to a laminin/polylysine mixture, neural somata and dendrites preferred the polylysine while axons prefer the mixture or the border between the two.

A fluid--structure interaction finite element analysis of pulsatile blood flow through a compliant stenotic artery.

Abstract: A new model is used to analyze the fully coupled problem of pulsatile blood flow through a compliant, axisymmetric stenotic artery using the finite element method. The model uses large displacement and large strain theory for the solid, and the full Navier-Stokes equations for the fluid. The effect of increasing area reduction on fluid dynamic and structural stresses is presented. Results show that pressure drop, peak wall shear stress, and maximum principal stress in the lesion all increase dramatically as the area reduction in the stenosis is increased from 51 to 89 percent. Further reductions in stenosis cross-sectional area, however, produce relatively little additional change in these parameters due to a concomitant reduction in flow rate caused by the losses in the constriction. Inner wall hoop stretch amplitude just distal to the stenosis also increases with increasing stenosis severity, as downstream pressures are reduced to a physiological minimum. The contraction of the artery distal to the stenosis generates a significant compressive stress on the downstream shoulder of the lesion. Dynamic narrowing of the stenosis is also seen, further augmenting area constriction at times of peak flow. Pressure drop results are found to compare well to an experimentally based theoretical curve, despite the assumption of laminar flow.

Convergence behavior of high-resolution finite element models of trabecular bone

Abstract: The convergence behavior of finite element models depends on the size of elements used, the element polynomial order, and on the complexity of the applied loads. For high-resolution models of trabecular bone, changes in architecture and density may also be important. The goal of this study was to investigate the influence of these factors on the convergence behavior of high-resolution models of trabecular bone. Two human vertebral and two bovine tibial trabecular bone specimens were modeled at four resolutions ranging from 20 to 80 microns and subjected to both compressive and shear loading. Results indicated that convergence behavior depended on both loading mode (axial versus shear) and volume fraction of the specimen. Compared to the 20 microns resolution, the differences in apparent Young's modulus at 40 microns resolution were less than 5 percent for all specimens, and for apparent shear modulus were less than 7 percent. By contrast, differences at 80 microns resolution in apparent modulus were up to 41 percent, depending on the specimen tested and loading mode. Overall, differences in apparent properties were always less than 10 percent when the ratio of mean trabecular thickness to element size was greater than four. Use of higher order elements did not improve the results. Tissue level parameters such as maximum principal strain did not converge. Tissue level strains converged when considered relative to a threshold value, but only if the strains were evaluated at Gauss points rather than element centroids. These findings indicate that good convergence can be obtained with this modeling technique, although element size should be chosen based on factors such as loading mode, mean trabecular thickness, and the particular output parameter of interest.

Dynamic lumped element response of the human fingerpad.

Abstract: The dynamic response of the fingerpad plays an important role in the tactile sensory response and precision manipulation, as well as in ergonomic design. This paper investigates the dynamic lumped element response of the human fingerpad in vivo to a compressive load. A flat probe indented the fingerpad at a constant velocity, then held a constant position. The resulting force (0-2 N) increased rapidly with indentation, then relaxed during the hold phase. A quasilinear viscoelastic model successfully explained the experimental data. The instantaneous elastic response increased exponentially with position, and the reduced relaxation function included three decaying exponentials ( with time constants of approximately 4 ms, 70 ms, and 1.4 s) plus a constant. The model was confirmed with data from sinusoidal displacement trajectories.

A method for planar biaxial mechanical testing that includes in-plane shear.

Abstract: A limitation in virtually all planar biaxial studies of soft tissues has been the inability to include the effects of in-plane shear. This is due to the inability of current mechanical testing devices to induce a state of in-plane shear, due to the added cost and complexity. In the current study, a straightforward method is presented for planar biaxial testing that induces a combined state of in-plane shear and normal strains. The method relies on rotation of the test specimen's material axes with respect to the device axes and on rotating carriages to allow the specimen to undergo in-plane shear freely. To demonstrate the method, five glutaraldehyde treated bovine pericardium specimens were prepared with their preferred fiber directions (defining the material axes) oriented at 45 deg to the device axes to induce a maximum shear state. The test protocol included a wide range of biaxial strain states, and the resulting biaxial data re-expressed in material axes coordinate system. The resulting biaxial data was then fit to the following strain energy function W: [equation: see text] where E'ij is the Green's strain tensor in the material axes coordinate system and c and Ai are constants. While W was able to fit the data very well, the constants A5 and A6 were found not to contribute significantly to the fit and were considered unnecessary to model the shear strain response. In conclusion, while not able to control the amount of shear strain independently or induce a state of pure shear, the method presented readily produces a state of simultaneous in-plane shear and normal strains. Further, the method is very general and can be applied to any anisotropic planar tissue that has identifiable material axes.

Dynamic contact of the human fingerpad against a flat surface.

Abstract: This paper investigates the dynamic, distributed pressure response of the human fingerpad in vivo when it first makes contact with an object. A flat probe was indented against the fingerpad at a 20 to 40 degree angle. Ramp-and-hold and sinusoidal displacement trajectories were applied to the fingerpad within a force range of 0-2 N. The dynamic spatial distribution of the pressure response was measured using a tactile array sensor. Both the local pressure variation and the total force exhibited nonlinear stiffness (exponential with displacement) and significant temporal relaxation. The shape of the contact pressure distribution could plausibly be described by an inverted paraboloid. A model based on the contact of a rigid plane (the object) and a linear viscoelastic sphere (the fingerpad), modified to include a nonlinear modulus of elasticity, can account for the principal features of the distributed pressure response.

Mechanical properties of collagen fascicles from the rabbit patellar tendon.

Abstract: Tensile and viscoelastic properties of collagen fascicles of approximately 300 μm in diameter, which were obtained from rabbit patellar tendons, were studied using a newly designed micro-tensile tester. Their cross-sectional areas were determined with a video dimension analyzer combined with a CCD camera and a low magnification microscope. There were no statistically significant differences in tensile properties among the fascicles obtained from six medial-to-lateral locations of the patellar tendon. Tangent modulus, tensile strength, and strain at failure of the fascicles determined at about 1.5 percent/s strain rate were 216 ± 68 MPa, 17.2 ± 4.1 MPa, and 10.9 ± 1.6 percent (mean ± S.D.), respectively. These properties were much different from those of bulk patellar tendons; for example, the tensile strength and strain at failure of these fascicles were 42 percent and 179 percent of those of bulk tendons, respectively. Tangent modulus and tensile strength of collagen fascicles determined at 1 percent/s strain rate were 35 percent larger than those at 0.01 percent/s. The strain at failure was independent of strain rate. Relaxation tests showed that the reduction of stress was approximately 25 percent at 300 seconds. These stress relaxation behavior and strain rate effects of collagen fascicles differed greatly from those of bulk tendons. The differences in tensile and viscoelastic properties between fascicles and bulk tendons may be attributable to ground substances, mechanical interaction between fascicles, and the difference of crimp structure of collagen fibrils.

Remodeling of a Collagenous Tissue at Fixed Lengths

Abstract: Mature tissues can often adapt to changes in their chemical, mechanical, or thermal environment. For example, in response to sustained increases or decreases in mechanical loads, some tissues grow and remodel so as to restore the stress or strain to its homeostatic state. Whereas most previous work addresses gross descriptors of tissue growth, this paper introduces a possible cell-mediated mechanism by which remodeling may occur in a soft connective tissue--that the kinetics of collagen deposition and degradation is similar regardless of the configuration of the body at which it occurs. The proposed theoretical framework applies to three-dimensional settings, but it is illustrated by focusing on the remodeling of a uniaxial collagenous tissue that is maintained at a fixed length for an extended period. It is shown that qualitative features expected of such remodeling (e.g., an increased compliance and increased stress-free length when remodeling occurs at an extended length) are easily realized. Growth and remodeling are complex phenomena, however, and are likely accomplished via multiple complementary mechanisms. There is a need, therefore, to identify other candidate mechanisms and, of course, to collect experimental data suitable for testing and refining the possible theories.

In vivo tracking of the human patella using cine phase contrast magnetic resonance imaging.

Abstract: Improper patellar tracking is often considered to be the cause of patellar-femoral pain. Unfortunately, our knowledge of patellar-femoral-tibial (knee) joint kinematics is severely limited due to a lack of three-dimensional, noninvasive, in vivo measurement techniques. This study presents the first large-scale, dynamic, three-dimensional, noninvasive, in vivo study of nonimpaired knee joint kinematics during volitional leg extensions. Cine-phase contrast magnetic resonance imaging was used to measure the velocity profiles of the patella, femur, and tibia in 18 unimpaired knees during leg extensions, resisted by a 34 N weight. Bone displacements were calculated through integration and then converted into three-dimensional orientation angles. We found that the patella displaced laterally, superiorly, and anteriorly as the knee extended. Further, patellar flexion lagged knee flexion, patellar tilt was variable, and patellar rotation was fairly constant throughout extension.

A Quantitative Investigation of Structure-Function Relationships in a Tendon Fascicle Model

Abstract: These studies sought to investigate quantitative relationships between the complex composite structure and mechanical properties of tendon. The isolated mouse tail tendon fascicle was chosen as an appropriate model for these so-called "structure-function" investigations. Specifically, collagen fibril diameters and mechanical properties were measured in fascicles from immature (3 week) control, adult (8 week) control, and adult (8 week) MovI3 transgenic mice. Results demonstrated a moderate correlation between mean fibril diameter and fascicle stiffness (r = 0.73, p = 0.001) and maximum load (r = 0.75, p < 0.001), whereas a weak correlation with fascicle modulus (r = 0.39, p = 0.11) and maximum stress (r = 0.48, p = 0.04). An analysis of pooled within-group correlations revealed no strong structure-function trends evidenced at the local or group level, indicating that correlations observed in the general structure-function analyses were due primarily to having three different experimental groups, rather than significant correlations of parameters within the groups.

Influence of cell deformation on leukocyte rolling adhesion in shear flow.

Abstract: Blood cell interaction with vascular endothelium is important in microcirculation, where rolling adhesion of circulating leukocytes along the surface of endothelial cells is a prerequisite for leukocyte emigration underflow conditions. HL-60 cell rolling adhesion to surface-immob ilized P-selectin in shear flow was investigated using a side-view flow chamber, which permitted measurements of cell deformation and cell-substrate contact length as well as cell rolling velocity. A two-dimensiona l model was developed based on the assumption that fluid energy input to a rolling cell was essentially distributed into two parts: cytoplasmic viscous dissipation, and energy needed to break adhesion bonds between the rolling cell and its substrate. The flow flelds of extracellular fluid and intracellular cytoplasm were solved using finite element methods with a deformable cell membrane represented by an elastic ring. The adhesion energy loss was calculated based on receptor-ligand kinetics equations. It was found that, as a result of shear-flow-induced cell deformation, cell-substrate contact area under high wall shear stresses (20 dyn/crrf) could be as much as twice of that under low stresses (0.5 dyn/cm'). An increase in contact area may cause more energy dissipation to both adhesion bonds and viscous cytoplasm, whereas the fluid energy input may decrease due to the flattened cell shape. Our model predicts that leukocyte rolling velocity will reach a plateau as shear stress increases, which agrees with both in vivo and in vitro experimental observations.

A numerical study of blood flow patterns in anatomically realistic and simplified end-to-side anastomoses.

Abstract: Purpose: Recently, some numerical and experimental studies of blood flow in large arteries have attempted to accurately replicate in vivo arterial geometries, while others have utilized simplified models. The objective of this study was to determine how much an anatomically realistic geometry can be simplified without the loss of significant hemodynamic information. Method: A human femoral-popliteal bypass graft was used to reconstruct an anatomically faithful finite element model of an end-to-side anastomosis. Nonideal geometric features of the model were removed in sequential steps to produce a series of successively simplified models. Blood flow patterns were numerically computed for each geometry, and the flow and wall shear stress fields were analyzed to determine the significance of each level of geometric simplification. Results: The removal of small local surface features and out-of-plane curvature did not significantly change the flow and wall shear stress distributions in the end-to-side anastomosis. Local changes in arterial caliber played a more significant role, depending upon the location and extent of the change. The graft-to-host artery diameter ratio was found to be a strong determinant of wall shear stress patterns in regions that are typically associated with disease processes. Conclusions: For the specific case of an end-to-side anastomosis, simplified models provide sufficient information for comparing hemodynamics with qualitative or averaged disease locations, provided the primary' geometric features are well replicated. The ratio of the graft-to-host artery diameter was shown to be the most important geometric feature. Secondary geometric features such as local arterial caliber changes, out-of-plane curvature, and small-scale surface topology are less important determinants of the wall shear stress patterns. However, if patient-specific disease information is available for the same arterial geometry, accurate replication of both primary and secondary geometric features is likely required.

Surfaces designed to control the projected area and shape of individual cells.

Abstract: Materials with spatially resolved surface chemistry were designed to isolate individual mammalian cells to determine the influence of projected area on specific cell functions (e.g., proliferation, cytoskeletal organization). Surfaces were fabricated using a photolithographic process resulting in islands of cell binding N-(2-aminoethyl)-3-aminopropyl-trimethoxystian (EDS) separated by a nonadhesive interpenetrating polymer network [poly (acrylamide-co-ethylene glycol); P(AAm-co-EG)]. The surfaces contained over 3800 adhesive islands/cm 2 , allowing for isolation of single cells with projected areas ranging from 100 μm 2 to 10,000 μm 2 . These surfaces provide a useful tool for researching how cell morphology and mechanical forces affect cell function.

Prediction of antagonistic muscle forces using inverse dynamic optimization during flexion/extension of the knee.

Abstract: This paper examined the feasibility of using different optimization criteria in inverse dynamic optimization to predict antagonistic muscle forces and joint reaction forces during isokinetic flexion/extension and isometric extension exercises of the knee. Both quadriceps and hamstrings muscle groups were included in this study. The knee joint motion included flexion/extension, varus/valgus, and internal/external rotations. Four linear, nonlinear, and physiological optimization criteria were utilized in the optimization procedure. All optimization criteria adopted in this paper were shown to be able to predict antagonistic muscle contraction during flexion and extension of the knee. The predicted muscle forces were compared in temporal patterns with EMG activities (averaged data measured from five subjects). Joint reaction forces were predicted to be similar using all optimization criteria. In comparison with previous studies, these results suggested that the kinematic information involved in the inverse dynamic optimization plays an important role in prediction of the recruitment of antagonistic muscles rather than the selection of a particular optimization criterion. Therefore, it might be concluded that a properly formulated inverse dynamic optimization procedure should describe the knee joint rotation in three orthogonal planes.

Interrelation of Creep and Relaxation: A Modeling Approach for Ligaments

Abstract: Experimental data (Thornton et al., 1997) show that relaxation proceeds more rapidly (a greater slope on a log-log scale) than creep in ligament, a fact not explained by linear viscoelasticity. An interrelation between creep and relaxation is therefore developed for ligaments based on a single-integral nonlinear superposition model. This interrelation differs from the convolution relation obtained by Laplace transforms for linear materials. We demonstrate via continuum concepts of nonlinear viscoelasticity that such a difference in rate between creep and relaxation phenomenologically occurs when the nonlinearity is of a strain-stiffening type, i.e., the stress-strain curve is concave up as observed in ligament. We also show that it is inconsistent to assume a Fung-type constitutive law (Fung, 1972) for both creep and relaxation. Using the published data of Thornton et al. (1997), the nonlinear interrelation developed herein predicts creep behavior from relaxation data well (R > or = 0.998). Although data are limited and the causal mechanisms associated with viscoelastic tissue behavior are complex, continuum concepts demonstrated here appear capable of interrelating creep and relaxation with fidelity.

Study on effect of graded facetectomy on change in lumbar motion segment torsional flexibility using three-dimensional continuum contact representation for facet joints

Abstract: Facet joints provide rigidity to the lumbar motion segment and thus protect the disk, particularly against torsional injury. A surgical procedure that fully or partially removes the facet joints (facetectomy) will decrease the mechanical stiffness of the motion segment, and potentially place the disk at risk of injury. Analytical models can be used to understand the effect of facet joints on motion segment stability. Using a facet joint model that represents the contact area as contact between two surfaces rather than as point contact, it was concluded that a substantial sudden change in rotational motion, due to applied torsion moment, was observed after 75 percent of any one of the facet joints was removed. Applied torsional moment loading produced coupled extension motion in the intact motion segment. This coupled motion also experienced a large change following complete unilateral facetectomy. Clinically, the present study showed that surgical intervention in the form of unilateral or bilateral total facetectomy might require fusion to reduce the primary torsion motion.

Optimization algorithm performance in determining optimal controls in human movement analyses.

Abstract: The objective of this study was to evaluate the performance of different multivariate optimization algorithms by solving a "tracking" problem using a forward dynamic model of pedaling. The tracking problem was defined as solving for the muscle controls (muscle stimulation onset, offset, and magnitude) that minimized the error between experimentally collected kinetic and kinematic data and the simulation results of pedaling at 90 rpm and 250 W. Three different algorithms were evaluated: a downhill simplex method, a gradient-based sequential quadratic programming algorithm, and a simulated annealing global optimization routine. The results showed that the simulated annealing algorithm performed for superior to the conventional routines by converging more rapidly and avoiding local minima.

The Effects of External Compression on Venous Blood Flow and Tissue Deformation in the Lower Leg

Abstract: External pneumatic compression of the lower legs is effective as prophylaxis against deep vein thrombosis. In a typical application, inflatable cuffs are wrapped around the patient's legs and periodically inflated to prevent stasis, accelerate venous blood flow, and enhance fibrinolysis. The purpose of this study was to examine the stress distribution within the tissues, and the corresponding venous blood flow and intravascular shear stress with different external compression modalities. A two-dimensional finite element analysis (FEA) was used to determine venous collapse as a function of internal (venous) pressure and the magnitude and spatial distribution of external (surface) pressure. Using the one-dimensional equations governing flow in a collapsible tube and the relations for venous collapse from the FEA, blood flow resulting from external compression was simulated. Tests were conducted to compare circumferentially symmetric (C) and asymmetric (A) compression and to examine distributions of pressure along the limb. Results show that A compression produces greater vessel collapse and generates larger blood flow velocities and shear stresses than C compression. The differences between axially uniform and graded-sequential compression are less marked than previously found, with uniform compression providing slightly greater peak flow velocities and shear stresses. The major advantage of graded-sequential compression is found at midcalf. Strains at the lumenal border are approximately 20 percent at an external pressure of 50 mmHg (6650 Pa) with all compression modalities.

A cellular solid criterion for predicting the axial-shear failure properties of bovine trabecular bone.

Abstract: In a long-term effort to develop a complete multi-axial failure criterion for human trabecular bone, the overall goal of this study was to compare the ability of a simple cellular solid mechanistic criterion versus the Tsai-Wu, Principal Strain, and von Mises phenomenological criteria--all normalized to minimize effects of interspecimen heterogeneity of strength--to predict the on-axis axial-shear failure properties of bovine trabecular bone. The Cellular Solid criterion that was developed here assumed that vertical trabeculae failed due to a linear superposition of axial compression/tension and bending stresses, induced by the apparent level axial and shear loading, respectively. Twenty-seven bovine tibial trabecular bone specimens were destructively tested on-axis without end artifacts, loaded either in combined tension-torsion (n = 10), compression-torsion (n = 11), or uniaxially (n = 6). For compression-shear, the mean (+/- S.D.) percentage errors between measured values and criterion predictions were 7.7 +/- 12.6 percent, 19.7 +/- 23.2 percent, 22.8 +/- 18.9 percent, and 82.4 +/- 64.5 percent for the Cellular Solid, Tsai-Wu, Principal Strain, and von Mises criteria, respectively; corresponding mean errors for tension-shear were -5.2 +/- 11.8 percent, 14.3 +/- 12.5 percent, 6.9 +/- 7.6 percent, and 57.7 +/- 46.3 percent. Statistical analysis indicated that the Cellular Solid criterion was the best performer for compression-shear, and performed as well as the Principal Strain criterion for tension-shear. These data should substantially improve the ability to predict axial-shear failure of dense trabecular bone. More importantly, the results firmly establish the importance of cellular solid analysis for understanding and predicting the multiaxial failure behavior of trabecular bone.

Extraction of Quasi-Linear Viscoelastic Parameters for Lower Limb Soft Tissues From Manual Indentation Experiment

Abstract: A manual indentation protocol was established to assess the quasi-linear viscoelastic (QLV) properties of lower limb soft tissues. The QLV parameters were extracted using a curve-fitting procedure on the experimental indentation data. The load-indentation responses were obtained using an ultrasound indentation apparatus with a hand-held pen-sized probe. Limb soft tissues at four sites of eight normal young subjects were tested in three body postures. Four QLV model parameters were extracted from the experimental data. The initial modulus E0 ranged from 0.22 kPa to 58.4 kPa. The nonlinear factor E1 ranged from 21.7 kPa to 547 kPa. The time constant tau ranged from 0.05 s to 8.93 s. The time-dependent materials parameter alpha ranged from 0.029 to 0.277. Large variations of the parameters were noted among subjects, sites, and postures.

Application of the Tsai–Wu Quadratic Multiaxial Failure Criterion to Bovine Trabecular Bone

Abstract: As a first step toward development of a multiaxial failure criterion for human trabecular bone, the Tsai-Wu quadratic failure criterion was modified as a function of apparent density and applied to bovine tibial trabecular bone. Previous data from uniaxial compressive, tensile, and torsion tests (n = 139 total) were combined with those from new triaxial tests (n = 17) to calibrate and then verify the criterion. Combinations of axial compression and radial pressure were used to produce the triaxial compressive stress states. All tests were performed with minimal end artifacts in the principal material coordinate system of the trabecular network. Results indicated that the stress interaction term F12 exhibited a strong nonlinear dependence on apparent density (r2 > 0.99), ranging from -0.126 MPa-2 at low densities (0.29 g/cm3) to 0.005 MPa-2 at high densities (0.63 g/cm3). After calibration and when used to predict behavior of new-specimens without any curve-fitting, the Tsai-Wu criterion had a mean (+/- SD) error of -32.6 +/- 10.6 percent. Except for the highest density triaxial specimens, most (15/17 specimens) failed at axial stresses close to their predicted uniaxial values, and some reinforcement for transverse loading was observed. We conclude that the Tsai-Wu quadratic criterion, as formulated here, is at best only a reasonable predictor of the multiaxial failure behavior of trabecular bone, and further work is required before it can be confidently applied to human bone.

Strain distribution in the layered wall of the esophagus.

Abstract: The function of the esophagus is to move food by peristaltic motion, which is the result of the interaction of the tissue forces in the esophageal wall and the hydrodynamic forces in the food bolus. To understand the tissue forces in the esophagus, it is necessary to know the zero-stress state of the esophagus, and the stress-strain relationships of the tissues. This article is addressed to the first topic: the representation of zero-stress state of the esophagus by the states of zero stress-resultant and zero bending moment of the mucosa-submucosa and the muscle layers. It is shown that at the states of zero stress-resultant and zero bending moment, these two layers are not tubes of smaller radii but are open sectors whose shapes are approximately cylindrical and more or less circular. When the sectors are approximated by circular sectors, we measured their radii, opening angles, and average thickness around the circumference. Data on the radii, thickness-to-radius ratios, and the opening angles of these sectors are presented. Knowing the zero-stress state of these two layers, we can compute the strain distribution in the wall at any in vivo state, as well as the residual strain in the esophageal wall at the no-load state. The results of the in vivo states are compared to those obtained by a conventional approach, which treats the esophageal wall as a homogeneous material, and to another popular simplification, which ignores the residual strains completely. It is shown that the errors caused by the homogeneous wall assumption are relatively minor, but those caused by ignoring the residual strains completely are severe.

An Evaluation of Pseudoelastic Descriptors Used in Arterial Mechanics

Abstract: Understanding how transmural distributions of stress relate to the mechanisms of vascular growth, remodeling, and disease necessitates computations that are based on a constitutive relation for the arterial wall. Although a number of candidate relations are in the literature, they have not been compared in detail. In this note, three commonly used descriptors of the passive behavior of common carotid arteries are compared using simple "thought experiments." It is shown that two of the three relations are inherently limited in the degree of anisotropy they allow,that each predicts a different anisotropy, and that one yields physically unrealistic predictions given many of the published values of the material parameters. Based on this comparison, it appears that the relation proposed by Chuong and Fung is the best available, though there may be a need to search for an alternate form, particularly for muscular arteries. The methods presented herein are offered as a guide to help the experimentalist identify alternative forms of pseudostrain-energy functions for arteries.