Bone Volume Fraction and Fabric Anisotropy Are Better Determinants of Trabecular Bone Stiffness than Other Morphological Variables

As our population ages, more individuals suffer from osteoporosis. This disease leads to impaired trabecular architecture and increased fracture risk. It is essential to understand how morphological and mechanical properties of the cancellous bone are related. Morphologyelasticity relationships based on bone volume fraction (BV/TV) and fabric anisotropy explain up to 98% of the variation in elastic properties. Yet, other morphological variables such as individual trabeculae segmentation (ITS) and trabecular bone score (TBS) could improve the stiffness predictions. A total of 743 micro-computed tomography reconstructions of cubic trabecular bone samples extracted from femur, radius, vertebrae and iliac crest were analysed. Their morphology was assessed via 25 variables and their stiffness tensor (inline image) was computed from six independent load cases using micro finite element analyses. Variance inflation factors were calculated to evaluate collinearity between morphological variables and decide upon their inclusion in morphology-elasticity relationships. The statistically admissible morphological variables were included in a multi-linear regression modelling the dependent variable inline image. The contribution of each independent variable was evaluated (ANOVA). Our results show that BV/TV is the best determinant of inline image (inline image=0.889), especially in combination with fabric (inline image=0.968). Including the other independent predictors hardly affected the amount of variance explained by the model (inline image=0.975). Across all anatomical sites, BV/TV explained 87% of the variance of the bone elastic properties. Fabric further described 10% of the bone stiffness, but the improvement in variance explanation by adding other independent factors was marginal (<1%). These findings confirm that BV/TV and fabric are the best determinants of trabecular bone stiffness and show, against common belief, that other morphological variables do not bring any further contribution. These overall conclusions remain to be confirmed for specific bone diseases and post-elastic properties.

Image-based biomechanical assessment of vertebral body and intervertebral disc in the human lumbar spine

Life expectancy continuously increases but our society faces age-related conditions. Among musculoskeletal diseases, osteoporosis associated with risk of vertebral fracture and degenerative intervertebral disc (IVD) are painful pathologies responsible for tremendous healthcare costs. Hence, reliable diagnostic tools are necessary to plan a treatment or follow up its efficacy. Yet, radiographic and MRI techniques, respectively clinical standards for evaluation of bone strength and IVD degeneration, are unspecific and not objective. Increasingly used in biomedical engineering, CT-based finite element (FE) models constitute the state-of-art for vertebral strength prediction. However, as non-invasive biomechanical evaluation and personalised FE models of the IVD are not available, rigid boundary conditions (BCs) are applied on the FE models to avoid uncertainties of disc degeneration that might bias the predictions. Moreover, considering the impact of low back pain, the biomechanical status of the IVD is needed as a criterion for early disc degeneration. Thus, the first FE study focuses on two rigid BCs applied on the vertebral bodies during compression test of cadaver vertebral bodies, vertebral sections and PMMA embedding. The second FE study highlights the large influence of the intervertebral disc’s compliance on the vertebral strength, damage distribution and its initiation. The third study introduces a new protocol for normalisation of the IVD stiffness in compression, torsion and bending using MRI-based data to account for its morphology. In the last study, a new criterion (Otsu threshold) for disc degeneration based on quantitative MRI data (axial T2 map) is proposed. The results show that vertebral strength and damage distribution computed with rigid BCs are identical. Yet, large discrepancies in strength and damage localisation were observed when the vertebral bodies were loaded via IVDs. The normalisation protocol attenuated the effect of geometry on the IVD stiffnesses without complete suppression. Finally, the Otsu threshold computed in the posterior part of annulus fibrosus was related to the disc biomechanics and meet objectivity and simplicity required for a clinical application. In conclusion, the stiffness normalisation protocol necessary for consistent IVD comparisons and the relation found between degeneration, mechanical response of the IVD and Otsu threshold lead the way for non-invasive evaluation biomechanical status of the IVD. As the FE prediction of vertebral strength is largely influenced by the IVD conditions, this data could also improve the future FE models of osteoporotic vertebra.

A numerical model for inert gas transport in the lung based on fractal airway morphology

Patients suffering from cystic fibrosis (CF) show thick secretions, mucus plugging and bronchiectasis in bronchial and alveolar ducts. This results in substantial structural changes of the airway morphology and heterogeneous ventilation. Disease progression and treatment effects are monitored by so-called gas washout tests, where the change in concentration of an inert gas is measured over a single or multiple breaths. The result of the tests based on the profile of the measured concentration is a marker for the severity of the ventilation inhomogeneity strongly affected by the airway morphology. However, it is hard to localize underlying obstructions to specific parts of the airways, especially if occurring in the lung periphery. In order to support the analysis of lung function tests (e.g. multi-breath washout), we developed a numerical model of the entire airway tree, coupling a lumped parameter model for the lung ventilation with a 4th-order accurate finite difference model of a 1D advection-diffusion equation for the transport of an inert gas. The boundary conditions for the flow problem comprise the pressure and flow profile at the mouth, which is typically known from clinical washout tests. The natural asymmetry of the lung morphology is approximated by a generic, fractal, asymmetric branching scheme which we applied for the conducting airways. A conducting airway ends when its dimension falls below a predefined limit. A model acinus is then connected to each terminal airway. The morphology of an acinus unit comprises a network of expandable cells. A regional, linear constitutive law describes the pressure-volume relation between the pleural gap and the acinus. The cyclic expansion (breathing) of each acinus unit depends on the resistance of the feeding airway and on the flow resistance and stiffness of the cells themselves. Special care was taken in the development of a conservative numerical scheme for the gas transport across bifurcations, handling spatially and temporally varying advective and diffusive fluxes over a wide range of scales. Implicit time integration was applied to account for the numerical stiffness resulting from the discretized transport equation. Local or regional modification of the airway dimension, resistance or tissue stiffness are introduced to mimic pathological airway restrictions typical for CF. This leads to a more heterogeneous ventilation of the model lung. As a result the concentration in some distal parts of the lung model remains increased for a longer duration. The inert gas concentration at the mouth towards the end of the expirations is composed of gas from regions with very different washout efficiency. This results in a steeper slope of the corresponding part of the washout profile.

Experimental validation of a nonlinear μ FE model based on cohesive-frictional plasticity for trabecular bone

Trabecular bone is a porous mineralized tissue playing a major load bearing role in the human body. Prediction of age-related and disease-related fractures and the behavior of bone implant systems needs a thorough understanding of its structure-mechanical property relationships, which can be obtained using microcomputed tomography-based finite element modeling. In this study, a nonlinear model for trabecular bone as a cohesive-frictional material was implemented in a large-scale computational framework and validated by comparison of μFE simulations with experimental tests in uniaxial tension and compression. A good correspondence of stiffness and yield points between simulations and experiments was found for a wide range of bone volume fraction and degree of anisotropy in both tension and compression using a non-calibrated, average set of material parameters. These results demonstrate the ability of the model to capture the effects leading to failure of bone for three anatomical sites and several donors, which may be used to determine the apparent behavior of trabecular bone and its evolution with age, disease, and treatment in the future.

Evaluation of HA-fibrin-based hydrogel for restoration of degenerated intervertebral disc

Introduction: Intervertebral disc degeneration is associated with loss of nucleus pulposus (NP) tissue and reduced disc height[1]. A number of therapies, including synthetic and natural biomaterials, have been developed to restore full disc function and to minimize the pain and disability caused by this disease. Fibrin-based biomaterials are used as a replacement for NP or as a cell carrier for tissue engineering approaches[2]. While the behavior of such gels is well-characterized from a material point of view, little is known about their contribution to intervertebral disc (IVD) restoration under dynamic loads. The aim of the present study is the evaluation of a hyaluronic acid fibrin-based hydrogel (ProCore) used to repair an in vitro model of disc degeneration under dynamic loading. Methods: In vitro model of disc degeneration was induced in intact coccygeal bovine IVD by papain digestion of the NP as previously described[3]. In order to characterize fibrin hydrogels, four experimental groups were considered: 1) intact IVD (control), 2) IVD injected with PBS, 3) injection of hydrogels in degenerative IVD and 4) injection of hydrogels in combination with human bone marrow-derived mesenchymal stem cells (MSC) in degenerative IVD. All of the groups were subjected to dynamic loading protocols consisting of 0.2MPa static compression superimposed with ±2° torsion at 0.2Hz for 8h per day and maintained for 7 days. Additionally, one group consisted of degenerative IVD injected with hydrogel and subjected to static compression. Disc heights were monitored after the duration of the loading and compared to the initial disc height. The macrostructure of the formed tissue and the cellular distribution was evaluated by histological means. Results: After one week of loading, the degenerative IVD filled with hydrogel in combination with MSC (dynamic load), hydrogels (dynamic load) and hydrogels (static load) showed a reduction in height by 30%, 15% and 20%, respectively, as compared to their initial disc height. Histological sections showed that the HA-fibrin gel fully occupied the nucleotomized region of the disc and that fibrin was effective in filling the discontinuities of the cavity region. Furthermore, the cells were homogenously distributed along the fibrin hydrogels after 7 days of loading. Discussion: In this study, we showed that fibrin hydrogels showed a good integration within the papain-induced model of disc degeneration and can withstand the applied loads. Fibrin hydrogels can contribute to disc restoration by possibly maintaining adequate stiffness of the tissue and thus preventing disorganization of the surrounding IVD. References: 1. Jarman, J.P., Arpinar, V.E., Baruah, D., Klein, A.P., Maiman, D.J., and Tugan Muftuler, L. (2014). Intervertebral disc height loss demonstrates the threshold of major pathological changes during degeneration. Eur Spine J . 2. Colombini, A., Ceriani, C., Banfi, G., Brayda-Bruno, M., and Moretti, M. (2014). Fibrin in intervertebral disc tissue engineering. Tissue Eng Part B Rev . 3. Chan, S.C., Bürki, A., Bonél, H.M., Benneker, L.M., and Gantenbein-Ritter, B. (2013). Papain-induced in vitro disc degeneration model for the study of injectable nucleus pulposus therapy. Spine J 13, 273-283. Acknowledgement We thank the Swiss National Science Foundation SNF #310030_153411 for funding.

1D.03: Inactive Matrix Gla Protein is Assocated with Renal Resistive Index in a Population-Based Study

OBJECTIVE Renal resistive index (RRI) varies directly with renal vascular stiffness and pulse pressure. RRI correlates positively with arteriolosclerosis in damaged kidneys and predicts progressive renal dysfunction. Matrix Gla-protein (MGP) is a vascular calcification inhibitor that needs vitamin K to be activated. Inactive MGP, known as desphospho-uncarboxylated MGP (dp-ucMGP), can be measured in plasma and has been associated with various cardiovascular (CV) markers, CV outcomes and mortality. In this study we hypothesize that increased RRI is associated with high levels of dp-ucMGP. DESIGN AND METHOD We recruited participants via a multi-center family-based cross-sectional study in Switzerland exploring the role of genes and kidney hemodynamics in blood pressure regulation. Dp-ucMGP was quantified in plasma samples by sandwich ELISA. Renal doppler sonography was performed using a standardized protocol to measure RRIs on 3 segmental arteries in each kidney. The mean of the 6 measures was reported. Multiple regression analysis was performed to estimate associations between RRI and dp-ucMGP adjusting for sex, age, pulse pressure, mean pressure, renal function and other CV risk factors. RESULTS We included 1035 participants in our analyses. Mean values were 0.64 ± 0.06 for RRI and 0.44 ± 0.21 (nmol/L) for dp-ucMGP. RRI was positively associated with dp-ucMGP both before and after adjustment for sex, age, body mass index, pulse pressure, mean pressure, heart rate, renal function, low and high density lipoprotein, smoking status, diabetes, blood pressure and cholesterol lowering drugs, and history of CV disease (P < 0.001). CONCLUSIONS RRI is independently and positively associated with high levels of dp-ucMGP after adjustment for pulse pressure and common CV risk factors. Further studies are needed to determine if vitamin K supplementation can have a positive effect on renal vascular stiffness and kidney function.

Rotation-Induced Breakdown of Torsional Quantum Control

Control of the torsional angles of nonrigid molecules is key for the development of emerging areas like molecular electronics and nanotechnology. Based on a rigorous calculation of the rotation-torsion-Stark energy levels of nonrigid biphenyl-like molecules, we show that, unlike previously believed, instantaneous rotation-torsion-Stark eigenstates of such molecules, interacting with a strong laser field, present a large degree of delocalization in the torsional coordinate even for the lowest energy states. This is due to a strong coupling between overall rotation and torsion leading to a breakdown of the torsional alignment. Thus, adiabatic control of changes on the planarity of this kind of molecule is essentially impossible unless the temperature is on the order of a few Kelvin.

Magnetic properties of nanostructured systems based on TbFe2

The aim of this work is to study the magnetic properties of annealed [Fe3Ga/TbFe2]n heterostructures grown by sputtering at room temperature. The interest of investigating multilayers comprised of TbFe2 and Fe3Ga is their complementary properties in terms of coercivity and magnetostriction. We have studied the thickness combination which optimizes the magnetic and magnetostrictive properties of the annealed multilayers. The crystallization of the Laves phase upon the thermal treatment in heterostructures with thick TbFe2 layers promotes the increase of the coercivity. This crystallization seems to be prevented by the low mechanical stiffness of the Fe3Ga. [Fe3Ga/TbFe2]n heterostructures show promising characteristics, λ of 340 ppm and a HC of 220 Oe, for the development of new magnetostrictive devices.

Acoustic displacement triangle based on the individual element test

A three node, displacement based, acoustic element is developed. In order to avoid spurious rotational modes, a higher order stiffness is introduced. The higher order stiffness is developed from an incompatible strain field which computes element volume changes under nodal rotational displacements fields. The higher order strain satisfies the IET requirements, non affecting convergence. The higher order stiffness is modulated, element by element, with a factor. Thus, the displacement based formulation is capable of placing the spurious rotational modes over the range of physical compressional modes that can be accurately captured by the mesh.

An interactive web-based learning assistant for structural analysis of continuous beams

There are significant levels of concern about the relevance and the difficulty of learning some issues on Strength of Materials and Structural Analysis. Most students of Continuum Mechanics and Structural Analysis in Civil Engineering usually point out some key learning aspects as especially difficult for acquiring specific skills. These key concepts entail comprehension difficulties but ease access and applicability to structural analysis in more advanced subjects. Likewise, some elusive but basic structural concepts, such as flexibility, stiffness or influence lines, are paramount for developing further skills required for advanced structural design: tall buildings, arch-type structures as well as bridges. As new curricular itineraries are currently being implemented, it appears appropriate to devise a repository of interactive web-based applications for training in those basic concepts. That will hopefully train the student to understand the complexity of such concepts, to develop intuitive knowledge on actual structural response and to improve their preparation for exams. In this work, a web-based learning assistant system for influence lines on continuous beams is presented. It consists of a collection of interactive user-friendly applications accessible via Web. It is performed in both Spanish and English languages. Rather than a “black box” system, the procedure involves open interaction with the student, who can simulate and virtually envisage the structural response. Thus, the student is enabled to set the geometric, topologic and mechanic layout of a continuous beam and to change or shift the loading and the support conditions. Simultaneously, the changes in the beam response prompt on the screen, so that the effects of the several issues involved in structural analysis become apparent. The system is performed through a set of web pages which encompasses interactive exercises and problems, written in JavaScript under JQuery and DyGraphs frameworks, given that their efficiency and graphic capabilities are renowned. Students can freely boost their self-study on this subject in order to face their exams more confidently. Besides, this collection is expected to be added to the "Virtual Lab of Continuum Mechanics" of the UPM, launched in 2013 (http://serviciosgate.upm.es/laboratoriosvirtuales/laboratorios/medios-continuos-en-construcci%C3%B3n)

Hierarchical and hybrid polymer nanocomposites based on carbon nanotubes and inorganic fullerene type nanoparticles

The influence of singlewalled carbon nanotubes (SWCNT) and inorganic fullerenelike tungsten disulfide nanoparticles (IFWS2) on the morphology and thermal, mechanical and electrical performance of multifunctional fibrereinforced polymer composites has been investigated. Significant improvements were observed in stiffness, strength and toughness in poly (ether ether ketone) (PEEK) / (SWCNT) / glass fibre (GF) laminates when a compatibilizer was used for wrapping the CNTs. Hybrid poly(phenylene sulphide) (PPS)/IFWS2/ carbon fibre (CF) reinforced polymer composites showed improved mechanical and tribological properties attributed to a synergetic effect between the IF nanoparticles and CF.

Nonlinear seismic behaviour of concrete arch bridges

Arch bridge structural solution has been known for centuries, in fact the simple nature of arch that require low tension and shear strength was an advantage as the simple materials like stone and brick were the only option back in ancient centuries. By the pass of time especially after industrial revolution, the new materials were adopted in construction of arch bridges to reach longer spans. Nowadays one long span arch bridge is made of steel, concrete or combination of these two as "CFST", as the result of using these high strength materials, very long spans can be achieved. The current record for longest arch belongs to Chaotianmen bridge over Yangtze river in China with 552 meters span made of steel and the longest reinforced concrete type is Wanxian bridge which also cross the Yangtze river through a 420 meters span. Today the designer is no longer limited by span length as long as arch bridge is the most applicable solution among other approaches, i.e. cable stayed and suspended bridges are more reasonable if very long span is desired. Like any super structure, the economical and architectural aspects in construction of a bridge is extremely important, in other words, as a narrower bridge has better appearance, it also require smaller volume of material which make the design more economical. Design of such bridge, beside the high strength materials, requires precise structural analysis approaches capable of integrating the combination of material behaviour and complex geometry of structure and various types of loads which may be applied to bridge during its service life. Depend on the design strategy, analysis may only evaluates the linear elastic behaviour of structure or consider the nonlinear properties as well. Although most of structures in the past were designed to act in their elastic range, the rapid increase in computational capacity allow us to consider different sources of nonlinearities in order to achieve a more realistic evaluations where the dynamic behaviour of bridge is important especially in seismic zones where large movements may occur or structure experience P - _ effect during the earthquake. The above mentioned type of analysis is computationally expensive and very time consuming. In recent years, several methods were proposed in order to resolve this problem. Discussion of recent developments on these methods and their application on long span concrete arch bridges is the main goal of this research. Accordingly available long span concrete arch bridges have been studied to gather the critical information about their geometrical aspects and properties of their materials. Based on concluded information, several concrete arch bridges were designed for further studies. The main span of these bridges range from 100 to 400 meters. The Structural analysis methods implemented in in this study are as following: Elastic Analysis: Direct Response History Analysis (DRHA): This method solves the direct equation of motion over time history of applied acceleration or imposed load in linear elastic range. Modal Response History Analysis (MRHA): Similar to DRHA, this method is also based on time history, but the equation of motion is simplified to single degree of freedom system and calculates the response of each mode independently. Performing this analysis require less time than DRHA. Modal Response Spectrum Analysis (MRSA): As it is obvious from its name, this method calculates the peak response of structure for each mode and combine them using modal combination rules based on the introduced spectra of ground motion. This method is expected to be fastest among Elastic analysis. Inelastic Analysis: Nonlinear Response History Analysis (NL-RHA): The most accurate strategy to address significant nonlinearities in structural dynamics is undoubtedly the nonlinear response history analysis which is similar to DRHA but extended to inelastic range by updating the stiffness matrix for every iteration. This onerous task, clearly increase the computational cost especially for unsymmetrical buildings that requires to be analyzed in a full 3D model for taking the torsional effects in to consideration. Modal Pushover Analysis (MPA): The Modal Pushover Analysis is basically the MRHA but extended to inelastic stage. After all, the MRHA cannot solve the system of dynamics because the resisting force fs(u; u_ ) is unknown for inelastic stage. The solution of MPA for this obstacle is using the previously recorded fs to evaluate system of dynamics. Extended Modal Pushover Analysis (EMPA): Expanded Modal pushover is a one of very recent proposed methods which evaluates response of structure under multi-directional excitation using the modal pushover analysis strategy. In one specific mode,the original pushover neglect the contribution of the directions different than characteristic one, this is reasonable in regular symmetric building but a structure with complex shape like long span arch bridges may go through strong modal coupling. This method intend to consider modal coupling while it take same time of computation as MPA. Coupled Nonlinear Static Pushover Analysis (CNSP): The EMPA includes the contribution of non-characteristic direction to the formal MPA procedure. However the static pushovers in EMPA are performed individually for every mode, accordingly the resulted values from different modes can be combined but this is only valid in elastic phase; as soon as any element in structure starts yielding the neutral axis of that section is no longer fixed for both response during the earthquake, meaning the longitudinal deflection unavoidably affect the transverse one or vice versa. To overcome this drawback, the CNSP suggests executing pushover analysis for governing modes of each direction at the same time. This strategy is estimated to be more accurate than MPA and EMPA, moreover the calculation time is reduced because only one pushover analysis is required. Regardless of the strategy, the accuracy of structural analysis is highly dependent on modelling and numerical integration approaches used in evaluation of each method. Therefore the widely used Finite Element Method is implemented in process of all analysis performed in this research. In order to address the study, chapter 2, starts with gathered information about constructed long span arch bridges, this chapter continuous with geometrical and material definition of new models. Chapter 3 provides the detailed information about structural analysis strategies; furthermore the step by step description of procedure of all methods is available in Appendix A. The document ends with the description of results and conclusion of chapter 4.

Torsional rigidity of single actin filaments and actin–actin bond breaking force under torsion measured directly by in vitro micromanipulation

Knowledge of the elastic properties of actin filaments is crucial for considering its role in muscle contraction, cellular motile events, and formation of cell shape. The stiffness of actin filaments in the directions of stretching and bending has been determined. In this study, we have directly determined the torsional rigidity and breaking force of single actin filaments by measuring the rotational Brownian motion and tensile strength using optical tweezers and microneedles, respectively. Rotational angular fluctuations of filaments supplied the torsional rigidity as (8.0 ± 1.2) × 10−26 Nm2. This value is similar to that deduced from the longitudinal rigidity, assuming the actin filament to be a homogeneous rod. The breaking force of the actin–actin bond was measured while twisting a filament through various angles using microneedles. The breaking force decreased greatly under twist, e.g., from 600–320 pN when filaments were turned through 90°, independent of the rotational direction. Our results indicate that an actin filament exhibits comparable flexibility in the rotational and longitudinal directions, but breaks more easily under torsional load.

Numerical and experimental investigation of structural stiffness influence on ratcheting convection cell in granular soils under cyclic loading

Lateral cyclic loaded structures in granular soils can lead to an accumulation of irreversible strains by changing their mechanical response (densification) and forming a closed convective cell in the upper layer of the bedding. In the present thesis the convective cell dimension, formation and grain migration inside this closed volume have been studied and presented in relation to structural stiffness and different loads. This relation was experimentally investigated by applying a cyclic lateral force to a scaled flexible vertical element embedded in dry granular soil. The model was monitored with a camera in order to derive the displacement field by means of the PIV technique. Modelling large soil deformation turns out to be difficult, using mesh-based methods. Consequently, a mesh-free approach (DEM) was chosen in order to investigate the granular flow with the aim of extracting interesting micromechanical information. In both the numerical and experimental analyses the effect of different loading magnitudes and different dimensions of the vertical element were considered. The main results regarded the different development, shape and dimensions of the convection cell and the surface settlements. Moreover, the Discrete Element Method has proven to give satisfactory results in the modelling of large deformation phenomena such as the ratcheting convective cell.

Contributions to Physics-Based Aeroservoelastic Uncertainty Analysis

Thesis (Ph.D.)--University of Washington, 2016-06