When aiming at assessing motor control development, natural walking (NW), and tandem walking (TW) are two locomotor tasks that allow analyzing different characteristics of motor control performance. NW is the reference locomotor task, expected to become more and more automatic with age. TW is a nonparadigmatic task used in clinics to highlight eventual impairments and to evaluate how a child deals with a new challenging motor experience. This work aims at investigating motor development in school-aged children, by assessing quantitatively their performance during TW and NW. Eighty children (6–10 years) participated in the study. Trunk acceleration data and nonlinear measures (recurrence quantification analysis (RQA), and multiscale entropy (MSE)) were used to characterize trunk postural control and motor complexity. The results were analyzed with respect to age and standard clinical assessment of TW (number of correct consecutive steps), by means of Spearman correlation coefficients. RQA and MSE allowed highlighting age-related changes in both postural control of the trunk and motor complexity, while classic standard assessment of TW resulted uniformly distributed in the different age groups. The present results suggest this quantitative approach as relevant when assessing the motor development in schoolchildren and complementary to standard clinical tests.

Delivery of biological fluids, such as surfactant solutions, into lungs is a major strategy to treat respiratory disorders including respiratory distress syndrome that is caused by insufficient or dysfunctional natural lung surfactant. The instilled solution forms liquid plugs in lung airways. The plugs propagate downstream in airways by inspired air or ventilation, continuously split at airway bifurcations to smaller daughter plugs, simultaneously lose mass from their trailing menisci, and eventually rupture. A uniform distribution of the instilled biofluid in lung airways is expected to increase the treatments success. The uniformity of distribution of instilled liquid in the lungs greatly depends on the splitting of liquid plugs between daughter airways, especially in the first few generations from which airways of different lobes of lungs emerge. To mechanistically understand this process, we developed a bioengineering approach to computationally design three-dimensional bifurcating airway models using morphometric data of human lungs, fabricate physical models, and examine dynamics of liquid plug splitting. We found that orientation of bifurcating airways has a major effect on the splitting of liquid plugs between daughter airways. Changing the relative gravitational orientation of daughter tubes with respect to the horizontal plane caused a more asymmetric splitting of liquid plugs. Increasing the propagation speed of plugs partially counteracted this effect. Using airway models of smaller dimensions reduced the asymmetry of plug splitting. This work provides a step toward developing delivery strategies for uniform distribution of therapeutic fluids in the lungs.

Kinetics assessment of the human head-arms-trunk (HAT) complex via a multisegment model is a useful tool for objective clinical evaluation of several pathological conditions. Inaccuracies in body segment parameters (BSPs) are a major source of uncertainty in the estimation of the joint moments associated with the multisegment HAT. Given the large intersubject variability, there is currently no comprehensive database for the estimation of BSPs for the HAT. We propose a nonlinear, multistep, optimization-based, noninvasive method for estimating individual-specific BSPs and calculating joint moments in a multisegment HAT model. Eleven nondisabled individuals participated in a trunk-bending experiment and their body motion was recorded using cameras and a force plate. A seven-segment model of the HAT was reconstructed for each participant. An initial guess of the BSPs was obtained by individual-specific scaling of the BSPs calculated from the male visible human (MVH) images. The intersegmental moments were calculated using both bottom-up and top-down inverse dynamics approaches. Our proposed method adjusted the scaled BSPs and center of pressure (COP) offsets to estimate optimal individual-specific BSPs that minimize the difference between the moments obtained by top-down and bottom-up inverse dynamics approaches. Our results indicate that the proposed method reduced the error in the net joint moment estimation (defined as the difference between the net joint moment calculated via bottom-up and top-down approaches) by 79.3% (median among participants). Our proposed method enables an optimized estimation of individual-specific BSPs and, consequently, a less erroneous assessment of the three-dimensional (3D) kinetics of a multisegment HAT model.

Osteoarthritis sufferers commonly have first metatarsophalangeal joint (MTPJ) problems in which articular surfaces are changed permanently due to fatigue. Therefore, medical devices for early diagnosis would increase the opportunity for prevention of disease progression. In previous studies on stiffness of the first MTPJ many details, although functionally of great importance, have not been fully considered including: design and size of the device, tribology consideration, and errors from device. Therefore, the motivation of our research was to enhance the device design by reducing the size of the device, and device design was enhanced by minimizing measurement errors through development of a new ergonomic left and right foot instrument located medial to the first MTPJ (instead of beneath the foot). The first MTPJ stiffness (N mm/kg radian) measurement was taken on 28 subjects with two replicates per subject by the same tester. The first MTPJ stiffness ranged from 3.49 to 14.42 N mm/kg radian with the mean (SD) value of 8.28 (3.15) N mm/kg radian for the left feet and 3.91 to 11.90 N mm/kg radian with the mean (SD) value of 7.65 (2.07) N mm/kg radian for the right feet. Reliability evaluation was measured using intraclass correlation coefficient and described an excellent reliability between two tests.

Rapid advancements in the multidisciplinary field of biomedical engineering (BME) require competitive engineers with skill sets in a broad range of subjects including biology, physiology, mechanics, circuits, and programming. Accordingly, such a need should be reflected in the training of BME students. Among those skills, computer programming is an essential tool that is used in a wide variety of applications. In this paper, we have provided our experience in incorporating project-based learning, a promising approach in active learning, for teaching computer programming to BME students. We describe a low-cost method for using physical, hands-on computing that directly relates to BME. Additionally, we detail our efforts to teach multiple programming languages in one semester and provide a detailed analysis of the outcomes. We also provide basic materials for other instructors to adapt to fit their own needs.

Clinical treatments of skeletal muscle weakness are hindered by a lack of an approach to evaluate individual muscle force. Intramuscular pressure (IMP) has shown a correlation to muscle force in vivo, but patient to patient and muscle to muscle variability results in difficulty of utilizing IMP to estimate muscle force. The goal of this work was to develop a finite element model of whole skeletal muscle that can predict IMP under passive and active conditions to further investigate the mechanisms of IMP variability. A previously validated hypervisco-poroelastic constitutive approach was modified to incorporate muscle activation through an inhomogeneous geometry. Model parameters were optimized to fit model stress to experimental data, and the resulting model fluid pressurization data were utilized for validation. Model fitting was excellent (root-mean-square error or RMSE <1.5 kPa for passive and active conditions), and IMP predictive capability was strong for both passive (RMSE 3.5 mmHg) and active (RMSE 10 mmHg at in vivo lengths) conditions. Additionally, model fluid pressure was affected by length under isometric conditions, as increases in stretch yielded decreases in fluid pressurization following a contraction, resulting from counteracting Poisson effects. Model pressure also varied spatially, with the highest gradients located near aponeuroses. These findings may explain variability of in vivo IMP measurements in the clinic, and thus help reduce this variability in future studies. Further development of this model to include isotonic contractions and muscle weakness would greatly benefit this work.

Wearable sensors embedded with inertial measurement units have become commonplace for the measurement of head impact biomechanics, but individual systems often suffer from a lack of measurement fidelity. While some researchers have focused on developing highly accurate, single sensor systems, we have taken a parallel approach in investigating optimal estimation techniques with multiple noisy sensors. In this work, we present a sensor network methodology that utilizes multiple skin patch sensors arranged on the head and combines their data to obtain a more accurate estimate than any individual sensor in the network. Our methodology visually localizes subject-specific sensor transformations, and based on rigid body assumptions, applies estimation algorithms to obtain a minimum mean squared error estimate. During mild soccer headers, individual skin patch sensors had over 100% error in peak angular velocity magnitude, angular acceleration magnitude, and linear acceleration magnitude. However, when properly networked using our visual localization and estimation methodology, we obtained kinematic estimates with median errors below 20%. While we demonstrate this methodology with skin patch sensors in mild soccer head impacts, the formulation can be generally applied to any dynamic scenario, such as measurement of cadaver head impact dynamics using arbitrarily placed sensors.

In-stent restenosis ails many patients who have undergone stenting. When the stented artery is a bifurcation, the intervention is particularly critical because of the complex stent geometry involved in these structures. Computational fluid dynamics (CFD) has been shown to be an effective approach when modeling blood flow behavior and understanding the mechanisms that underlie in-stent restenosis. However, these CFD models require validation through experimental data in order to be reliable. It is with this purpose in mind that we performed particle image velocimetry (PIV) measurements of velocity fields within flows through a simplified coronary bifurcation. Although the flow in this simplified bifurcation differs from the actual blood flow, it emulates the main fluid dynamic mechanisms found in hemodynamic flow. Experimental measurements were performed for several stenting techniques in both steady and unsteady flow conditions. The test conditions were strictly controlled, and uncertainty was accurately predicted. The results obtained in this research represent readily accessible, easy to emulate, detailed velocity fields and geometry, and they have been successfully used to validate our numerical model. These data can be used as a benchmark for further development of numerical CFD modeling in terms of comparison of the main flow pattern characteristics.

Quantitative computed tomography (qCT) relies on calibrated bone mineral density data. If a calibration phantom is absent from the CT scan, post hoc calibration becomes necessary. Scanning a calibration phantom after-the-fact and applying that calibration to uncalibrated scans has been used previously. Alternatively, the estimated density is known to vary with CT settings, suggesting that it may be possible to predict the calibration terms using CT settings. This study compares a novel CT setting regression method for post hoc calibration to standard and post hoc phantom-only calibrations. Five cadaveric upper limbs were scanned at 11 combinations of peak tube voltage and current (80–140 kV and 100–300 mA) with two calibration phantoms. Density calibrations were performed for the cadaver scans, and scans of the phantoms alone. Stepwise linear regression determined if the calibration equation terms were predictable using peak tube voltage and current. Peak tube voltage, but not current, was significantly correlated with regression calibration terms. Calibration equation slope was significantly related to the type of phantom (p < 0.001), calibration method (p = 0.026), and peak tube voltage (p < 0.001), but not current (p = 1.000). The calibration equation vertical intercept was significantly related to the type of phantom (p < 0.001), and peak tube voltage (p = 0.006), but not calibration method (p = 0.682), or current (p = 0.822). Accordingly, regression can correlate peak tube voltage with density calibration terms. Suggesting that, while standard qCT calibration is preferable, regression calibration may be an acceptable post hoc method when necessary.

The determination of the elastic modulus of bone is important in studying the response of bone to loading and is determined using a destructive three-point bending method. Reference point indentation (RPI), with one cycle of indentation, offers a nondestructive alternative to determine the elastic modulus. While the elastic modulus could be determined using a nondestructive procedure for ex vivo experiments, for in vivo testing, the three-point bending technique may not be practical and hence RPI is viewed as a potential alternative and explored in this study. Using the RPI measurements, total indentation distance (TID), creep indentation distance, indentation force, and the unloading slope, we have developed a numerical analysis procedure using the Oliver–Pharr (O/P) method to estimate the indentation elastic modulus. Two methods were used to determine the area function: (1) Oliver–Pharr (O/P—based on a numerical procedure) and (2) geometric (based on the calculation of the projected area of indentation). The indentation moduli of polymethyl methacrylate (PMMA) calculated by the O/P (3.49–3.68 GPa) and geometric (3.33–3.49 GPa) methods were similar to values in literature (3.5–4 GPa). In a study using femurs from C57Bl/6 mice of different ages and genders, the three-point bending modulus was lower than the indentation modulus. In femurs from 4 to 5 months old TOPGAL mice, we found that the indentation modulus from the geometric (5.61 ± 1.25 GPa) and O/P (5.53 ± 1.27 GPa) methods was higher than the three-point bending modulus (5.28 ± 0.34 GPa). In females, the indentation modulus from the geometric (7.45 ± 0.86 GPa) and O/P (7.46 ± 0.92 GPa) methods was also higher than the three-point bending modulus (7.33 ± 1.13 GPa). We can conclude from this study that the RPI determined values are relatively close to three-point bending values.

Cytoplasmic viscosity-dependent margination of red blood cells (RBC) for flow inside microchannels was studied using numerical simulations, and the results were verified with microfluidic experiments. Wide range of suspension volume fractions or hematocrits was considered in this study. Lattice Boltzmann method for fluid-phase coupled with spectrin-link method for RBC membrane deformation was used for accurate analysis of cell margination. RBC margination behavior shows strong dependence on the internal viscosity of the RBCs. At equilibrium, RBCs with higher internal viscosity marginate closer to the channel wall and the RBCs with normal internal viscosity migrate to the central core of the channel. Same margination pattern has been verified through experiments conducted with straight channel microfluidic devices. Segregation between RBCs of different internal viscosity is enhanced as the shear rate and the hematocrit increases. Stronger separation between normal RBCs and RBCs with high internal viscosity is obtained as the width of a high aspect ratio channel is reduced. Overall, the margination behavior of RBCs with different internal viscosities resembles with the margination behavior of RBCs with different levels of deformability. Observations from this work will be useful in designing microfluidic devices for separating the subpopulations of RBCs with different levels of deformability that appear in many hematologic diseases such as sickle cell disease (SCD), malaria, or cancer.

The inertial properties of a helmet play an important role in both athletic performance and head protection. In this study, we measured the inertial properties of 37 football helmets, a National Operating Committee on Standards for Athletic Equipment (NOCSAE) size 7¼ headform, and a 50th percentile male Hybrid III dummy head. The helmet measurements were taken with the helmets placed on the Hybrid III dummy head. The center of gravity and moment of inertia were measured about six axes (x, y, z, xy, yz, and xz), allowing for a complete description of the inertial properties of the head and helmets. Total helmet mass averaged 1834±231 g, split between the shell (1377±200 g) and the facemask (457±101 g). On average, the football helmets weighed 41±5% as much as the Hybrid III dummy head. The center of gravity of the helmeted head was 1.1±3.0 mm anterior and 10.3±1.9 mm superior to the center of gravity of the bare head. The moment of inertia of the helmeted head was approximately 2.2±0.2 times greater than the bare head about all axes.

Although uniaxial tensile testing is commonly used to evaluate failure properties of vascular tissue, there is no established protocol for specimen shape or gripping method. Large percentages of specimens are reported to fail near the clamp and can potentially confound the studies, or, if discarded will result in sample waste. The objective of this study is to identify sample geometry and clamping conditions that can achieve consistent failure in the midregion of small arterial specimens, even for vessels from older individuals. Failure location was assessed in 17 dogbone specimens from human cerebral and sheep carotid arteries using soft inserts. For comparison with commonly used protocols, an additional 22 rectangular samples were tested using either sandpaper or foam tape inserts. Midsample failure was achieved in 94% of the dogbone specimens, while only 14% of the rectangular samples failed in the midregion, the other 86% failing close to the clamps. Additionally, we found midregion failure was more likely to be abrupt, caused by cracking or necking. In contrast, clamp failure was more likely to be gradual and included a delamination mode not seen in midregion failure. Hence, this work provides an approach that can be used to obtain consistent midspecimen failure, avoiding confounding clamp-related artifacts. Furthermore, with consistent midregion failure, studies can be designed to image the failure process in small vascular samples providing valuable quantitative information about changes to collagen and elastin structure during the failure process.

Much of our understanding of the role of elbow ligaments to overall joint biomechanics has been developed through in vitro cadaver studies using joint motion simulators. The principle of superposition can be used to indirectly compute the force contributions of ligaments during prescribed motions. Previous studies have analyzed the contribution of different soft tissue structures to the stability of human elbow joints, but have limitations in evaluating the loads sustained by those tissues. This paper introduces a unique, hybrid experimental-computational technique for measuring and simulating the biomechanical contributions of ligaments to elbow joint kinematics and stability. in vitro testing of cadaveric joints is enhanced by the incorporation of fully parametric virtual ligaments, which are used in place of the native joint stabilizers to characterize the contribution of elbow ligaments during simple flexion–extension (FE) motions using the principle of superposition. Our results support previously reported findings that the anterior medial collateral ligament (AMCL) and the radial collateral ligament (RCL) are the primary soft tissue stabilizers for the elbow joint. Tuned virtual ligaments employed in this study were able to restore the kinematics and laxity of elbows to within 2 deg of native joint behavior. The hybrid framework presented in this study demonstrates promising capabilities in measuring the biomechanical contribution of ligamentous structures to joint stability.

Percutaneous pedicle screw fixation (PPSF) is a well-known minimally invasive surgery (MIS) employed in the treatment of thoracolumbar burst fractures (TBF). However, hardware failure and loss of angular correction are common limitations caused by the poor support of the anterior column of the spine. Balloon kyphoplasty (KP) is another MIS that was successfully used in the treatment of compression fractures by augmenting the injured vertebral body with cement. To overcome the limitations of stand-alone PPSF, it was suggested to augment PPSF with KP as a surgical treatment of TBF. Yet, little is known about the biomechanical alteration occurred to the spine after performing such procedure. The objective of this study was to evaluate and compare the immediate post-operative biomechanical performance of stand-alone PPSF, stand-alone-KP, and KP-augmented PPSF procedures. Novel three-dimensional (3D) finite element (FE) models of the thoracolumbar junction that describes the fractured spine and the three investigated procedures were developed and tested under mechanical loading conditions. The spinal stiffness, stresses at the implanted hardware, and the intradiscal pressure at the upper and lower segments were measured and compared. The results showed no major differences in the measured parameters between stand-alone PPSF and KP-augmented PPSF procedures, and demonstrated that the stand-alone KP may restore the stiffness of the intact spine. Accordingly, there was no immediate post-operative biomechanical advantage in augmenting PPSF with KP when compared to stand-alone PPSF, and fatigue testing may be required to evaluate the long-term biomechanical performance of such procedures.

Inspiratory flow in a multigeneration pig lung airways was numerically studied at a steady inlet flow rate of 3.2 × 10 −4 m 3 /s corresponding to a Reynolds number of 1150 in the trachea. The model was validated by comparing velocity distributions with previous measurements and simulations in simplified airway geometries. Simulation results provided detailed maps of the axial and secondary flow patterns at different cross sections of the airway tree. The vortex core regions in the airways were visualized using absolute helicity values and suggested the presence of secondary flow vortices where two counter-rotating vortices were observed at the main bifurcation and in many other bifurcations. Both laminar and turbulent flows were considered. Results showed that axial and secondary flows were comparable in the laminar and turbulent cases. Turbulent kinetic energy (TKE) vanished in the more distal airways, which indicates that the flow in these airways approaches laminar flow conditions. The simulation results suggested viscous pressure drop values comparable to earlier studies. The monopodial asymmetric nature of airway branching in pigs resulted in airflow patterns that are different from the less asymmetric human airways. The major daughters of the pig airways tended to have high airflow ratios, which may lead to different particle distribution and sound generation patterns. These differences need to be taken into consideration when interpreting the results of animal studies involving pigs before generalizing these results to humans.

There is lack of investigation capturing the complex mechanical interaction of tissue-engineered intervertebral disk (IVD) constructs in physiologically relevant environmental conditions. In this study, mechanical characterization of anisotropic electrospinning (ES) substrates made of polycaprolactone (PCL) was carried out in wet and dry conditions and viability of human bone marrow derived mesenchymal stem cells (hMSCs) seeded within double layers of ES PCL were also studied. Cyclic compression of IVD-like constructs composed of an agarose core confined by ES PCL double layers was implemented using a bioreactor and the cellular response to the mechanical stimulation was evaluated. Tensile tests showed decrease of elastic modulus of ES PCL as the angle of stretching increased, and at 60 deg stretching angle in wet, the maximum ultimate tensile strength (UTS) was observed. Based on the configuration of IVD-like constructs, the calculated circumferential stress experienced by the ES PCL double layers was 40 times of the vertical compressive stress. Confined compression of IVD-like constructs at 5% and 10% displacement dramatically reduced cell viability, particularly at 10%, although cell presence in small and isolated area can still be observed after mechanical conditioning. Hence, material mechanical properties of tissue-engineered scaffolds, composed of fibril structure of polymer with low melting point, are affected by the testing condition. Circumferential stress induced by axial compressive stimulation, conveyed to the ES PCL double layer wrapped around an agarose core, can affect the viability of cells seeded at the interface, depending on the mechanical configuration and magnitude of the load.

Atomic force microscopy (AFM) has been used to measure cellular stiffness at different osmolarities to investigate the effect of osmotic pressure on cells. However, substantial direct evidence is essential to clarify the phenomena derived from the experimental results. This study used both the single-point and force mapping methods to measure the effective Young's modulus of the cell by using temporal and spatial information. The single-point force measurements confirmed the positive correlation between cellular stiffness and osmolarity. The force mapping measurements provided local stiffness on the cellular surface and identified the cytoskeleton distribution underneath the plasma membrane. At hyper-osmolarity, the cytoskeleton was observed to cover most of the area underneath the plasma membrane, and the effective Young's modulus on the area with cytoskeleton support was determined to be higher than that at iso-osmolarity. The overall increase in cellular Young's modulus confirmed the occurrence of cytoskeleton compression at hyper-osmolarity. On the other hand, although the average Young's modulus at hypo-osmolarity was lower than that at iso-osmolarity, we observed that the local Young's modulus measured on the areas with cytoskeleton support remained similar from iso-osmolarity to hypo-osmolarity. The reduction of the average Young's modulus at hypo-osmolarity was attributed to reduced cytoskeleton coverage underneath the plasma membrane.

Head impact exposure in popular contact sports is not well understood, especially in the youth population, despite recent advances in impact-sensing technology which has allowed widespread collection of real-time head impact data. Previous studies indicate that a custom-instrumented mouthpiece is a superior method for collecting accurate head acceleration data. The objective of this study was to evaluate the efficacy of mounting a sensor device inside an acrylic retainer form factor to measure six-degrees-of-freedom (6DOF) head kinematic response. This study compares 6DOF mouthpiece kinematics at the head center of gravity (CG) to kinematics measured by an anthropomorphic test device (ATD). This study found that when instrumentation is mounted in the rigid retainer form factor, there is good coupling with the upper dentition and highly accurate kinematic results compared to the ATD. Peak head kinematics were correlated with r 2 > 0.98 for both rotational velocity and linear acceleration and r 2 = 0.93 for rotational acceleration. These results indicate that a rigid retainer-based form factor is an accurate and promising method of collecting head impact data. This device can be used to study head impacts in helmeted contact sports such as football, hockey, and lacrosse as well as nonhelmeted sports such as soccer and basketball. Understanding the magnitude and frequency of impacts sustained in various sports using an accurate head impact sensor, such as the one presented in this study, will improve our understanding of head impact exposure and sports-related concussion.

Simulations of soft tissues require accurate and robust constitutive models, whose form is derived from carefully designed experimental studies. For such investigations of membranes or thin specimens, planar biaxial systems have been used extensively. Yet, all such systems remain limited in their ability to: (1) fully prescribe in-plane deformation gradient tensor F 2 D , (2) ensure homogeneity of the applied deformation, and (3) be able to accommodate sufficiently small specimens to ensure a reasonable degree of material homogeneity. To address these issues, we have developed a novel planar biaxial testing device that overcomes these difficulties and is capable of full control of the in-plane deformation gradient tensor F 2 D and of testing specimens as small as ∼4 mm × ∼4 mm. Individual actuation of the specimen attachment points, combined with a robust real-time feedback control, enabled the device to enforce any arbitrary F 2 D with a high degree of accuracy and homogeneity. Results from extensive device validation trials and example tissues illustrated the ability of the device to perform as designed and gather data needed for developing and validating constitutive models. Examples included the murine aortic tissues, allowing for investigators to take advantage of the genetic manipulation of murine disease models. These capabilities highlight the potential of the device to serve as a platform for informing and verifying the results of inverse models and for conducting robust, controlled investigation into the biomechanics of very local behaviors of soft tissues and membrane biomaterials.