J Biomech Eng. 2000;122(4):297-303. doi:10.1115/1.1286559.

A computational model that accounts for blood–tissue interaction under physiological flow conditions was developed and applied to a thin-walled model of the left heart. This model consisted of the left ventricle, left atrium, and pulmonary vein flow. The input functions for the model included the pulmonary vein driving pressure and time-dependent relationship for changes in chamber tissue properties during the simulation. The Immersed Boundary Method was used for the interaction of the tissue and blood in response to fluid forces and changes in tissue pathophysiology, and the fluid mass and momentum conservation equations were solved using Patankar’s Semi-Implicit Method for Pressure Linked Equations (SIMPLE). This model was used to examine the flow fields in the left heart under abnormal diastolic conditions of delayed ventricular relaxation, delayed ventricular relaxation with increased ventricular stiffness, and delayed ventricular relaxation with an increased atrial contraction. The results obtained from the left heart model were compared to clinically observed diastolic flow conditions, and to the results from simulations of normal diastolic function in this model 1. Cases involving impairment of diastolic function were modeled with changes to the input functions for fiber relaxation/contraction of the chambers. The three cases of diastolic dysfunction investigated agreed with the changes in diastolic flow fields seen clinically. The effect of delayed relaxation was to decrease the early filling magnitude, and this decrease was larger when the stiffness of the ventricle was increased. Also, increasing the contraction of the atrium during atrial systole resulted in a higher late filling velocity and atrial pressure. The results show that dysfunction can be modeled by changing the relationships for fiber resting-length and/or stiffness. This provides confidence in future modeling of disease, especially changes to chamber properties to examine the effect of local dysfunction on global flow fields. [S0148-0731(00)00104-7]

Commentary by Dr. Valentin Fuster
J Biomech Eng. 2000;122(4):304-309. doi:10.1115/1.1287171.

Clinical studies using transcranial Doppler ultrasonography in patients with mechanical heart valves (MHV) have detected gaseous emboli. The relationship of gaseous emboli release and cavitation on MHV has been a subject of debate in the literature. To study the influence of cavitation and gas content on the formation and growth of stable gas bubbles, a mock circulatory loop, which employed a Medtronic-Hall pyrolytic carbon disk valve in the mitral position, was used. A high-speed video camera allowed observation of cavitation and gas bubble release on the inflow valve surfaces as a function of cavitation intensity and carbon dioxide (CO2) concentration, while an ultrasonic monitoring system scanned the aortic outflow tract to quantify gas bubble production by calculating the gray scale levels of the images. In the absence of cavitation, no stable gas bubbles were formed. When gas bubbles were formed, they were first seen a few milliseconds after and in the vicinity of cavitation collapse. The volume of the gas bubbles detected in the aortic track increased with both increased CO2 and increased cavitation intensity. No correlation was observed between O2 concentration and bubble volume. We conclude that cavitation is an essential precursor to stable gas bubble formation, and CO2, the most soluble blood gas, is the major component of stable gas bubbles. [S0148-0731(00)00204-1]

Commentary by Dr. Valentin Fuster
J Biomech Eng. 2000;122(4):310-320. doi:10.1115/1.1287157.

To evaluate the local hemodynamic implications of coronary artery balloon angioplasty, computational fluid dynamics (CFD) was applied in a group of patients previously reported by [Wilson et al. (1988), 77 , pp. 873–885] with representative stenosis geometry post-angioplasty and with measured values of coronary flow reserve returning to a normal range (3.6±0.3). During undisturbed flow in the absence of diagnostic catheter sensors within the lesions, the computed mean pressure drop Δp̃ was only about 1 mmHg at basal flow, and increased moderately to about 8 mmHg for hyperemic flow. Corresponding elevated levels of mean wall shear stress in the midthroat region of the residual stenoses, which are common after angioplasty procedures, increased from about 60 to 290 dynes/cm2 during hyperemia. The computations (R̃ee≃100–400;αe=2.25) indicated that the pulsatile flow field was principally quasi-steady during the cardiac cycle, but there was phase lag in the pressure drop−mean velocity (Δp−ū) relation. Time-averaged pressure drop values, Δp̃, were about 20 percent higher than calculated pressure drop values, Δps, for steady flow, similar to previous in vitro measurements by Cho et al. (1983). In the throat region, viscous effects were confined to the near-wall region, and entrance effects were evident during the cardiac cycle. Proximal to the lesion, velocity profiles deviated from parabolic shape at lower velocities during the cardiac cycle. The flow field was very complex in the oscillatory separated flow reattachment region in the distal vessel where pressure recovery occurred. These results may also serve as a useful reference against catheter-measured pressure drops and velocity ratios (hemodynamic endpoints) and arteriographic (anatomic) endpoints post-angioplasty. Some comparisons to previous studies of flow through stenoses models are also shown for perspective purposes. [S0148-0731(00)00304-6]

Commentary by Dr. Valentin Fuster
J Biomech Eng. 2000;122(4):321-326. doi:10.1115/1.1286560.

The purpose of this study was to determine the hemolytic potentials of discrete bubble cavitation and attached cavitation. To generate controlled cavitation events, a venturi-geometry hydrodynamic device, called a Cavitation Susceptibility Meter (CSM), was constructed. A comparison between the hemolytic potential of discrete bubble cavitation and attached cavitation was investigated with a single-pass flow apparatus and a recirculating flow apparatus, both utilizing the CSM. An analytical model, based on spherical bubble dynamics, was developed for predicting the hemolysis caused by discrete bubble cavitation. Experimentally, discrete bubble cavitation did not correlate with a measurable increase in plasma-free hemoglobin (PFHb), as predicted by the analytical model. However, attached cavitation did result in significant PFHb generation. The rate of PFHb generation scaled inversely with the Cavitation number at a constant flow rate, suggesting that the size of the attached cavity was the dominant hemolytic factor. [S0148-0731(00)00404-0]

Commentary by Dr. Valentin Fuster
J Biomech Eng. 2000;122(4):327-335. doi:10.1115/1.1287158.

We have formulated the first constitutive model to describe the complete measured planar biaxial stress–strain relationship of the native and glutaraldehyde-treated aortic valve cusp using a structurally guided approach. When applied to native, zero-pressure fixed, and low-pressure fixed cusps, only three parameters were needed to simulate fully the highly anisotropic, and nonlinear in-plane biaxial mechanical behavior. Differences in the behavior of the native and zero- and low-pressure fixed cusps were found to be primarily due to changes in the effective fiber stress–strain behavior. Further, the model was able to account for the effects of small (<10 deg) misalignments in the cuspal specimens with respect to the biaxial test axes that increased the accuracy of the model material parameters. Although based upon a simplified cuspal structure, the model underscored the role of the angular orientation of the fibers that completely accounted for extreme mechanical anisotropy and pronounced axial coupling. Knowledge of the mechanics of the aortic cusp derived from this model may aid in the understanding of fatigue damage in bioprosthetic heart valves and, potentially, lay the groundwork for the design of tissue-engineered scaffolds for replacement heart valves. [S0148-0731(00)00504-5]

Commentary by Dr. Valentin Fuster
J Biomech Eng. 2000;122(4):336-346. doi:10.1115/1.1286316.

The main objective of this study is to determine the nature of electric fields inside articular cartilage while accounting for the effects of both streaming potential and diffusion potential. Specifically, we solve two tissue mechano-electrochemical problems using the triphasic theories developed by Lai et al. (1991, ASME J. Biomech Eng., 113 , pp. 245–258) and Gu et al. (1998, ASME J. Biomech. Eng., 120 , pp. 169–180) (1) the steady one-dimensional permeation problem; and (2) the transient one-dimensional ramped-displacement, confined-compression, stress-relaxation problem (both in an open circuit condition) so as to be able to calculate the compressive strain, the electric potential, and the fixed charged density (FCD) inside cartilage. Our calculations show that in these two technically important problems, the diffusion potential effects compete against the flow-induced kinetic effects (streaming potential) for dominance of the electric potential inside the tissue. For softer tissues of similar FCD (i.e., lower aggregate modulus), the diffusion potential effects are enhanced when the tissue is being compressed (i.e., increasing its FCD in a nonuniform manner) either by direct compression or by drag-induced compaction; indeed, the diffusion potential effect may dominate over the streaming potential effect. The polarity of the electric potential field is in the same direction of interstitial fluid flow when streaming potential dominates, and in the opposite direction of fluid flow when diffusion potential dominates. For physiologically realistic articular cartilage material parameters, the polarity of electric potential across the tissue on the outside (surface to surface) may be opposite to the polarity across the tissue on the inside (surface to surface). Since the electromechanical signals that chodrocytes perceive in situ are the stresses, strains, pressures and the electric field generated inside the extracellular matrix when the tissue is deformed, the results from this study offer new challenges for the understanding of possible mechanisms that control chondrocyte biosyntheses. [S0148-0731(00)00604-X]

Commentary by Dr. Valentin Fuster
J Biomech Eng. 2000;122(4):347-353. doi:10.1115/1.1286561.

A microstructural model of cartilage was developed to investigate the relative contribution of tissue matrix components to its elastostatic properties. Cartilage was depicted as a tensed collagen lattice pressurized by the Donnan osmotic swelling pressure of proteoglycans. As a first step in modeling the collagen lattice, two-dimensional networks of tensed, elastic, interconnected cables were studied as conceptual models. The models were subjected to the boundary conditions of confined compression and stress–strain curves and elastic moduli were obtained as a function of a two-dimensional equivalent of swelling pressure. Model predictions were compared to equilibrium confined compression moduli of calf cartilage obtained at different bath concentrations ranging from 0.01 to 0.50 M NaCl. It was found that a triangular cable network provided the most consistent correspondence to the experimental data. The model showed that the cartilage collagen network remained tensed under large confined compression strains and could therefore support shear stress. The model also predicted that the elastic moduli increased with increasing swelling pressure in a manner qualitatively similar to experimental observations. Although the model did not preclude potential contributions of other tissue components and mechanisms, the consistency of model predictions with experimental observations suggests that the cartilage collagen network, prestressed by proteoglycan swelling pressure, plays an important role in supporting compression. [S0148-0731(00)00704-4]

Commentary by Dr. Valentin Fuster
J Biomech Eng. 2000;122(4):354-363. doi:10.1115/1.1288207.

Current brain deformation models have predominantly reflected solid constitutive relationships generated from empirical ex vivo data and have largely overlooked interstitial hydrodynamic effects. In the context of a technique to update images intraoperatively for image-guided neuronavigation, we have developed and quantified the deformation characteristics of a three-dimensional porous media finite element model of brain deformation in vivo. Results have demonstrated at least 75–85 percent predictive capability, but have also indicated that interstitial hydrodynamics are important. In this paper we investigate interstitial pressure transient behavior in brain tissue when subjected to an acute surgical load consistent with neurosurgical events. Data are presented from three in vivo porcine experiments where subsurface tissue deformation and interhemispheric pressure gradients were measured under conditions of an applied mechanical deformation and then compared to calculations with our three-dimensional brain model. Results demonstrate that porous-media consolidation captures the hydraulic behavior of brain tissue subjected to comparable surgical loads and that the experimental protocol causes minimal trauma to porcine brain tissue. Working values for hydraulic conductivity of white and gray matter are also reported and an assessment of transient pressure gradient effects with respect to deformation is provided. [S0148-0731(00)00804-9]

Commentary by Dr. Valentin Fuster
J Biomech Eng. 2000;122(4):364-371. doi:10.1115/1.1287160.

The mechanical properties of the adult human skull are well documented, but little information is available for the infant skull. To determine the age-dependent changes in skull properties, we tested human and porcine infant cranial bone in three-point bending. The measurement of elastic modulus in the human and porcine infant cranial bone agrees with and extends previous published data [McPherson, G. K., and Kriewall, T. J. (1980), J. Biomech., 13 , pp. 9–16] for human infant cranial bone. After confirming that the porcine and human cranial bone properties were comparable, additional tensile and three-point bending studies were conducted on porcine cranial bone and suture. Comparisons of the porcine infant data with previously published adult human data demonstrate that the elastic modulus, ultimate stress, and energy absorbed to failure increase, and the ultimate strain decreases with age for cranial bone. Likewise, we conclude that the elastic modulus, ultimate stress, and energy absorbed to failure increase with age for sutures. We constructed two finite element models of an idealized one-month old infant head, one with pediatric and the other adult skull properties, and subjected them to impact loading to investigate the contribution of the cranial bone properties on the intracranial tissue deformation pattern. The computational simulations demonstrate that the comparatively compliant skull and membranous suture properties of the infant brain case are associated with large cranial shape changes, and a more diffuse pattern of brain distortion than when the skull takes on adult properties. These studies are a fundamental initial step in predicting the unique mechanical response of the pediatric skull to traumatic loads associated with head injury and, thus, for defining head injury thresholds for children. [S0148-0731(00)00904-3]

Commentary by Dr. Valentin Fuster
J Biomech Eng. 2000;122(4):372-379. doi:10.1115/1.1288208.

Thermal therapy of benign prostatic hyperplasia requires accurate prediction of the temperature distribution induced by the heating within the prostatic tissue. In this study, the Pennes bioheat transfer equation was used to model the transient heat transfer inside the canine prostate during transurethral microwave thermal therapy. Incorporating the specific absorption rate of microwave energy in tissue, a closed-form analytical solution was obtained. Good agreement was found between the theoretical predictions and in-vivo experimental results. Effects of blood perfusion and the cooling at the urethral wall on the temperature rise were investigated within the prostate during heating. The peak intraprostatic temperatures attained by application of 5, 10, or 15 W microwave power were predicted to be 38°C,41°C, and 44°C. Results from this study will help optimize the thermal dose that can be applied to target tissue during the therapy. [S0148-0731(00)01004-9]

Commentary by Dr. Valentin Fuster
J Biomech Eng. 2000;122(4):380-386. doi:10.1115/1.1286562.

In order to assess the significance of the dynamics of neural control signals for the rise time of muscle moment, simulations of isometric and dynamic plantar flexion contractions were performed using electromyographic signals (EMG signals) of m. triceps surae as input. When excitation dynamics of the muscle model was optimized for an M-wave of the medial head of m. gastrocnemius (GM), the model was able to make reasonable predictions of the rise time of muscle moment during voluntary isometric plantar flexion contractions on the basis of voluntary GM EMG signals. The rise time of muscle moment in the model was for the greater part determined by the amplitude of the first EMG burst. For dynamic jumplike movements of the ankle joint, however, no relationship between rise time of muscle moment in the experiment and muscle moment predicted by the model on the basis of GM EMG signals was found. Since rise time of muscle moment varied over a small range for this movement, it cannot be completely excluded that stimulation dynamics plays a role in control of these simple single-joint movements. [S0148-0731(00)01104-3]

Commentary by Dr. Valentin Fuster
J Biomech Eng. 2000;122(4):387-393. doi:10.1115/1.1287161.

Although it is well accepted that bone tissue metabolism is regulated by external mechanical loads, it remains unclear to what load-induced physical signals bone cells respond. In this study, a novel computer-controlled stretch device and parallel plate flow chamber were employed to investigate cytosolic calcium (Ca2+i) mobilization in response to a range of dynamic substrate strain levels (0.1–10 percent, 1 Hz) and oscillating fluid flow (2 N/m2, 1 Hz). In addition, we quantified the effect of dynamic substrate strain and oscillating fluid flow on the expression of mRNA for the bone matrix protein osteopontin (OPN). Our data demonstrate that continuum strain levels observed for routine physical activities (<0.5 percent) do not induce Ca2+i responses in osteoblastic cells in vitro. However, there was a significant increase in the number of responding cells at larger strain levels. Moreover, we found no change in osteopontin mRNA level in response to 0.5 percent strain at 1 Hz. In contrast, oscillating fluid flow predicted to occur in the lacunar–canalicular system due to routine physical activities (2 N/m2, 1 Hz) caused significant increases in both Ca2+i and OPN mRNA. These data suggest that, relative to fluid flow, substrate deformation may play less of a role in bone cell mechanotransduction associated with bone adaptation to routine loads. [S0148-0731(00)01204-8]

Commentary by Dr. Valentin Fuster
J Biomech Eng. 2000;122(4):394-401. doi:10.1115/1.1286563.

A new cell-level finite element formulation is presented and used to investigate how epithelia and other planar collections of viscous cells might deform during events such as embryo morphogenesis and wound healing. Forces arising from cytoskeletal components, cytoplasm viscosity, and cell-cell adhesions are included. Individual cells are modeled using multiple finite elements, and cell rearrangements can occur. Simulations of cell-sheet stretching indicate that the initial stages of sheet stretching are characterized by changes in cell shape, while subsequent stages are governed by cell rearrangement. Inferences can be made from the simulations about the forces that act in real cell sheets when suitable experimental data are available. [S0148-0731(00)01404-7]

Commentary by Dr. Valentin Fuster
J Biomech Eng. 2000;122(4):402-407. doi:10.1115/1.1288205.

Finite element–based computer simulations are used to investigate a number of phenomena, including tissue engulfment, cell sorting, and checkerboard-pattern formation, exhibited by heterotypic cell aggregates. The simulations show that these phenomena can be driven by a single equivalent force, namely a surface (or interfacial) tension, that results from cytoskeletal components and cell–cell adhesions. They also reveal that tissue engulfment, cell sorting, and checkerboard-pattern formation involve several discernible mechanical features or stages. With the aid of analytical arguments, we identify the conditions necessary for each of these phenomena. These findings are consistent with previous experimental investigations and computer simulations, but pose significant challenges to current theories of cell sorting and tissue engulfment. [S0148-0731(00)01304-2]

Commentary by Dr. Valentin Fuster
J Biomech Eng. 2000;122(4):408-415. doi:10.1115/1.1286564.

Polymer-based composites are widely used in restorative dentistry as alternatives to metals and ceramics to fill cavities in teeth. They adhere to the walls of the cavity in the tooth, thus forming a composite body consisting of dentine, enamel, and composite resin. Geometric discontinuities along the interfaces between these materials can induce singularities in the stress field, which in turn lead to premature failure of the restoration. In the present investigation, a complex stress function technique is employed to derive the order of the stress singularity. It is shown that the order of the singularity depends on both the material properties of the restorative material and the local geometry of the cavity. It is also shown that the singularity in the stress field can be avoided through careful design of the cavity shape. The results presented correlate well with experimental results reported in the literature. [S0148-0731(00)01504-1]

Commentary by Dr. Valentin Fuster
J Biomech Eng. 2000;122(4):416-422. doi:10.1115/1.1286317.

The paper presents a novel method for recording amplitude and phase of 6D-vibrations of a spatial pendulum over a wide frequency range (10 Hz up to 20 kHz). The six degrees of freedom of the pendulum mass were monitored by three electrodynamic stereo pickups. At rest, the tips of the needles and the pendulum’s center of mass defined the reference system with respect to which the oscillations of the mass were recorded in terms of their amplitudes and phases. Its small dimensions, constant transfer characteristics, linearity, high dynamics, and virtual lack of reaction onto the moving system over the entire frequency range provided the advantages of the measuring system. This method was used to analyze the spatial 6D-vibrations of the head of a cemented femoral hip endoprosthesis when the femur was stimulated to bending vibrations. The head of the prosthesis carried out axial rotational vibrations at every frequency used to stimulate the femur. The amplitudes of the axial rotations of the cortical bone were small in comparison to the ones of the prosthesis head, indicating that axial rotational vibrations following femur bending vibrations mainly stressed the spongiosa and the cement layer. This was observed over the entire frequency range, including at the low frequencies relevant for gait. Over the low-frequency range, as well as at some of the higher resonance frequencies, stationary instantaneous helical axes characterized the vibrations. The measurements suggest the mechanism that the interface “implant-bone” may already be stressed by axial torsional loads when the femur is loaded by bending impacts that are known to occur during walking. [S0148-0731(00)01604-6]

Commentary by Dr. Valentin Fuster
J Biomech Eng. 2000;122(4):423-429. doi:10.1115/1.1287162.

Rosette strain gages indicate shear and principal strains at specific points, whereas photoelastic coatings provide shear strain information over a broad area. Information regarding bone loading and load transfer from a prosthetic implant to adjacent bone can be obtained using either strain-measuring technique on loaded femora. This study compared proximal femoral strains derived from photoelastic coatings to those obtained from rosette strain gages applied directly to the bone in order to determine the relationships between photoelastic shear strains and rosette shear and principal strains. Photoelastic shear strains underestimated rosette shear strains and exceeded the larger of the rosette principal strains. Principal strains derived from photoelastic coatings augmented with strain separator gages underestimated their rosette counterparts in most instances. Correlation was strong and nearly linear for all measures, indicating that photoelastic coatings can accurately express proportional strain changes despite imperfect agreement in absolute strain magnitudes. The best agreement between absolute strain magnitudes occurred in the proximal medial, or calcar, region. Understanding the relationships between the various measures obtained using the two strain measurement methods will allow more accurate estimates of actual strains to be made from photoelastic coatings. [S0148-0731(00)01704-0]

Commentary by Dr. Valentin Fuster
J Biomech Eng. 2000;122(4):430-436. doi:10.1115/1.1286318.

Glenoid component loosening is the dominant cause of failure in total shoulder arthroplasty. It is presumed that loosening in the glenoid is caused by high stresses in the cement layer. Several anchorage systems have been designed with the aim of reducing the loosening rate, the two major categories being “keeled” fixation and “pegged” fixation. However, no three-dimensional finite element analysis has been performed to quantify the stresses in the cement or to compare the different glenoid prosthesis anchorage systems. The objective of this study was to determine the stresses in the cement layer and surrounding bone for glenoid replacement components. A three-dimensional model of the scapula was generated using CT data for geometry and material property definition. Keeled and pegged designs were inserted into the glenoid, surrounded by a 1-mm layer of bone cement. A 90 deg arm abduction load with a full muscle and joint load was applied, following van der Helm (1994). Deformations of the prosthesis, stresses in the cement, and stresses in the bone were calculated. Stresses were also calculated for a simulated case of rheumatoid arthritis (RA) in which bone properties were modified to reflect that condition. A maximum principal stress-based failure model was used to predict what quantity of the cement is at risk of failure at the levels of stress computed. The prediction is that 94 percent (pegged prosthesis) and 68 percent (keeled prosthesis) of the cement has a greater than 95 percent probability of survival in normal bone. In RA bone, however, the situation is reversed where 86 percent (pegged prosthesis) and 99 percent (keeled prosthesis) of the cement has a greater than 95 percent probability of survival. Bone stresses are shown to be not much affected by the prosthesis design, except at the tip of the central peg or keel. It is concluded that a “pegged” anchorage system is superior for normal bone, whereas a “keeled” anchorage system is superior for RA bone. [S0148-0731(00)01804-5]

Commentary by Dr. Valentin Fuster
J Biomech Eng. 2000;122(4):437-445. doi:10.1115/1.1286677.

Joint moment estimation using the traditional inverse dynamics analysis presents two challenging problems, which limit its reliability. First, the quality of the computed moments depends directly on unreliable estimates of the segment accelerations obtained numerically by differentiating noisy marker measurements. Second, the representation of joint moments from combined video and force plate measurements belongs to a class of ill-posed problems, which does not possess a unique solution. This paper presents a well-posed representation derived from an embedded constraint equation. The proposed method, referred to as the embedded constraint representation (ECR), provides unique moment estimates, which satisfy all measurement constraints and boundary conditions and require fewer acceleration components than the traditional inverse dynamics method. Specifically, for an n-segment open chain planar system, the ECR requires n−3 acceleration components as compared to 3(n−1) components required by the traditional (from ground up) inverse dynamics analysis. Based on a simulated experiment using a simple three-segment model, the precision of the ECR is evaluated at different noise levels and compared to the traditional inverse dynamics technique. At the lowest noise levels, the inverse dynamics method is up to 50 percent more accurate while at the highest noise levels the ECR method is up to 100 percent more accurate. The ECR results over the entire range of noise levels reveals an average improvement on the order 20 percent in estimating the moments distal to the force plate and no significant improvement in estimating moments proximal to the force plate. The new method is particularly advantageous in a combined video, force plate, and accelerometery sensing strategy. [S0148-0731(00)01904-X]

Commentary by Dr. Valentin Fuster
J Biomech Eng. 2000;122(4):446-452. doi:10.1115/1.1286678.

Bicycle pedaling has been studied from both a motor control and an equipment setup and design perspective. In both cases, although the dynamics of the bicycle drive system may have an influence on the results, a thorough understanding of the dynamics has not been developed. This study pursued three objectives related to developing such an understanding. The first was to identify the limitations of the inertial/frictional drive system model commonly used in the literature. The second was to investigate the advantages of an inertial/frictional/compliant model. The final objective was to use these models to develop a methodology for configuring a laboratory ergometer to emulate the drive system dynamics of road riding. Experimental data collected from the resulting road-riding emulator and from a standard ergometer confirmed that the inertial/frictional model is adequate for most studies of road-riding mechanics or pedaling coordination. However, the compliant model was needed to reproduce the phase shift in crank angle variations observed experimentally when emulating the high inertia of road riding. This finding may be significant for equipment setup and design studies where crank kinematic variations are important or for motor control studies where fine control issues are of interest. [S0148-0731(00)02004-5]

Commentary by Dr. Valentin Fuster


J Biomech Eng. 2000;122(4):453-454. doi:10.1115/1.1286565.

Several types of mechanical cardiac prostheses have been constructed with Delrin occluders, a material that is subject to osmotic swelling. The leaflets are designed to expand to specific tolerances when immersed in blood. The synthetic blood analogs commonly used in vitro contain hydrophilic compounds that can alter the osmotic expansion of the Delrin occluders. A static leak test chamber was employed to illustrate the effects of various test fluids on the sustained regurgitation phase of Delrin valves. [S0148-0731(00)02104-X]

Commentary by Dr. Valentin Fuster
J Biomech Eng. 2000;122(4):454-456. doi:10.1115/1.1288210.

Quantifying the stress distribution through the arterial wall is essential to studies of arterial growth and disease. Previous studies have shown that both residual stress, as measured by opening angle, and differing material properties for the media-intima and the adventitial layers affect the transmural circumferential stress θ) distribution. Because a lack of comprehensive data on a single species and artery has led to combinations from multiple sources, this study determined the sensitivity of σθ to published variations in both opening angle and layer thickness data. We fit material properties to previously published experimental data for pressure–diameter relations and opening angles of rabbit carotid artery, and predicted σθ through the arterial wall at physiologic conditions. Using a one-layer model, the ratio of σθ at the internal wall to the mean σθ decreased from 2.34 to 0.98 as the opening angle increased from 60 to 130 deg. In a two-layer model using a 95 deg opening angle, mean σθ in the adventitia increased (112 percent for 25 percent adventitia) and mean σθ in the media decreased (47 percent for 25 percent adventitia). These results suggest that both residual stress and wall layers have important effects on transmural stress distribution. Thus, experimental measurements of loading curves, opening angles, and wall composition from the same species and artery are needed to accurately predict the transmural stress distribution in the arterial wall. [S0148-0731(00)02204-4]

J Biomech Eng. 2000;122(4):457-460. doi:10.1115/1.1287163.

This article describes the design and development of a system that is capable of quantifying the thermal comfort of bicycle helmets. The motivation for the development of the system stems from the desire both to increase helmet use and to provide the designer with a quantitative method of evaluating the thermal comfort of a helmet. The system consists of a heated mannequin head form, a heated reference sphere, a small wind tunnel, and a data acquisition system. Both the head form and the reference sphere were instrumented with thermocouples. The system is capable of simulating riding speeds ranging from 4.5–15.5 m/s. A cooling effectiveness, C1, that is independent of both ambient conditions and wind velocity is defined as a measure of how well the helmet ventilates as compared to the reference sphere. The system was validated by testing six commercially available bicycle helmets manufactured between approximately 1992 and 1998. [S0148-0731(00)02304-9]


Commentary by Dr. Valentin Fuster

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In