Tissue viscoelasticity is related to tissue composition but may not fully predict the apparent-level viscoelasticity in human trabecular bone – An experimental and finite element study

A Review: Biomechanical Aspects of the Fallopian Tube Relevant to its Function in Fertility

The fallopian tube (FT) plays a crucial role in the reproductive process by providing an ideal biomechanical and biochemical environment for fertilization and early embryo development. Despite its importance, the biomechanical functions of the FT that originate from its morphological aspects, and ultrastructural aspects, as well as the mechanical properties of FT, have not been studied nor used sufficiently, which limits the understanding of fertilization, mechanotrasduction, and mechanobiology during embryo development, as well as the replication of the FT in laboratory settings for infertility treatments. This paper reviews and revives valuable information on human FT reported in medical literature in the past five decades relevant to the biomechanical aspects of FT. In this review, we summarized the current state of knowledge concerning the morphological, ultrastructural aspects, and mechanical properties of the human FT. We also investigate the potential arising from a thorough consideration of the biomechanical functions and exploring often neglected mechanical aspects. Our investigation encompasses both macroscopic measurements (such as length, diameter, and thickness) and microscopic measurements (including the height of epithelial cells, the percentage of ciliated cells, cilia structure, and ciliary beat frequency). Our primary focus has been on healthy women of reproductive age. We have examined various measurement techniques, encompassing conventional metrology, 2D histological data as well as new spatial measurement techniques such as micro-CT.

Bone Organoids: Recent Advances and Future Challenges

Bone defects stemming from tumorous growths, traumatic events, and diverse conditions present a profound conundrum in clinical practice and research. While bone has the inherent ability to regenerate, substantial bone anomalies require bone regeneration techniques. Bone organoids represent a new concept in this field, involving the 3D self-assembly of bone-associated stem cells guided in vitro with or without extracellular matrix material, resulting in a tissue that mimics the structural, functional, and genetic properties of native bone tissue. Within the scientific panorama, bone organoids ascend to an esteemed status, securing significant experimental endorsement. Through a synthesis of current literature and pioneering studies, this review offers a comprehensive survey of the bone organoid paradigm, delves into the quintessential architecture and ontogeny of bone, and highlights the latest progress in bone organoid fabrication. Further, existing challenges and prospective directions for future research are identified, advocating for interdisciplinary collaboration to fully harness the potential of this burgeoning domain. Conclusively, as bone organoid technology continues to mature, its implications for both clinical and research landscapes are poised to be profound.

Efficient nonlinear homogenization of bones using a cluster-based model order reduction technique

We present a reduced order model for efficient nonlinear homogenization of bones, accounting for strength difference effects and containing some well-known plasticity models (like von Mises or Drucker-Prager) as special cases. The reduced order homogenization is done by using a cluster-based model order reduction technique, called cluster-based nonuniform transformation field analysis. For an offline phase, a space-time decomposition is performed on the mesoscopic plastic strain fields, while a clustering analysis is employed for a spatial decomposition of the mesoscale RVE model. A volumetric-deviatoric split is additionally introduced to capture the enriched characteristics of the mesoscopic plastic strain fields. For an online analysis, the reduced order model is formulated in a unified minimization problem, which is compatible with a large variety of material models. Both cortical and trabecular bones are considered for numerical experiments. Compared to conventional FE-based RVE computations, the developed reduced order model renders a considerable acceleration rate beyond10 3 , while maintaining a sufficient accuracy level.

Effective Governing Equations for Viscoelastic Composites

We derive the governing equations for the overall behaviour of linear viscoelastic composites comprising two families of elastic inclusions, subphases and/or fibres, and an incompressible Newtonian fluid interacting with the solid phases at the microscale. We assume that the distance between each of the subphases is very small in comparison to the length of the whole material (the macroscale). We can exploit this sharp scale separation and apply the asymptotic (periodic) homogenization method (AHM) which decouples spatial scales and leads to the derivation of the new homogenised model. It does this via upscaling the fluid-structure interaction problem that arises between the multiple elastic phases and the fluid. As we do not assume that the fluid flow is characterised by a parabolic profile, the new macroscale model, which consists of partial differential equations, is of Kelvin-Voigt viscoelastic type (rather than poroelastic). The novel model has coefficients that encode the properties of the microstructure and are to be computed by solving a single local differential fluid-structure interaction (FSI) problem where the solid and the fluid phases are all present and described by the one problem. The model reduces to the case described by Burridge and Keller (1981) when there is only one elastic phase in contact with the fluid. This model is applicable when the distance between adjacent phases is smaller than the average radius of the fluid flowing in the pores, which can be the case for various highly heterogeneous systems encountered in real-world (e.g., biological, or geological) scenarios of interest.

The Relationship of Whole Human Vertebral Body Creep to Bone Density and Texture via Clinically Available Imaging Modalities

Creep deformation of human vertebrae accumulates under physiological levels of load and is understood to contribute to the progression toward clinically observable vertebral fracture. However, little information is available in terms of clinically measurable predictors of creep behavior in human vertebrae. In this study, creep tests were performed on 22 human cadaveric T12 vertebrae (13 male, 9 female; age 41-90). Areal and volumetric bone density parameters were measured from the same specimens using dual x-ray absorptiometry and high resolution computed tomography. Image textural analyses (which probe the organization of image intensities within the cancellous bone in low resolution clinical imaging) were performed using digital tomosynthesis (DTS) images. Multiple regression models were constructed to examine the relationship between creep properties and bone density and DTS image textural parameters. For the standard clinical imaging configuration, models including DTS derived image textural parameters alone were generally more explanatory (adjusted R²: 0.14-0.68) than those with bone density parameters forced in the models (adjusted R²: 0.17-0.61). Metrics of textural heterogeneity and anisotropy presented as the most explanatory imaging markers for creep deformation and recovery from creep. These metrics of image texture may help provide, independent from bone mass, important clinically measurable indicators of the time dependent deformation of human vertebrae.

Embracing Mechanobiology in Next Generation Organ-On-A-Chip Models of Bone Metastasis

Bone metastasis in breast cancer is associated with high mortality. Biomechanical cues presented by the extracellular matrix play a vital role in driving cancer metastasis. The lack of in vitro models that recapitulate the mechanical aspects of the in vivo microenvironment hinders the development of novel targeted therapies. Organ-on-a-chip (OOAC) platforms have recently emerged as a new generation of in vitro models that can mimic cell-cell interactions, enable control over fluid flow and allow the introduction of mechanical cues. Biomaterials used within OOAC platforms can determine the physical microenvironment that cells reside in and affect their behavior, adhesion, and localization. Refining the design of OOAC platforms to recreate microenvironmental regulation of metastasis and probe cell-matrix interactions will advance our understanding of breast cancer metastasis and support the development of next-generation metastasis-on-a-chip platforms. In this mini-review, we discuss the role of mechanobiology on the behavior of breast cancer and bone-residing cells, summarize the current capabilities of OOAC platforms for modeling breast cancer metastasis to bone, and highlight design opportunities offered by the incorporation of mechanobiological cues in these platforms.

Characterization and comparison of hyper-viscoelastic properties of normal and osteoporotic bone using stress-relaxation experiment

Bone tissue behavior under various loads is nonlinear elastic due to irreversible energy absorption. Also, viscoelasticity is one of the most important properties of bone which is very important in dynamic analyses and helps a lot in making artificial bone. In this study, rat tibia bone specimens were subjected to compression stress-relaxation test for normal (n = 5) and osteoporotic (n = 5) groups in order to characterize their mechanical properties using finite element modeling coupled with an optimization algorithm. Using this method, the structural equation parameters for the Neo-Hookean model and the Prony series coefficients were used to describe the hyper-elastic and the viscoelastic behavior of specimens, respectively; moreover, the properties of materials including the bulk, shear and Young's moduli for both groups were obtained and compared. The shear modulus was also gained as a function of time. In addition, the percentage of stress reduction and its relation to the initial stress were investigated for specimens. Finally, the effect of changes in each of the parameters of the hyper-viscoelastic structural equation on the force response was determined. Results of this study can be used in predicting the transient response and dynamic analysis of the bone.

The effects of despeckling filters on pore size measurements in collagen scaffold micro-CT data

Micro-CT is often used to assess the characteristics of porous structures such as tissue engineering scaffolds and trabecular bone. Prior to analysis, micro-CT images can be thresholded and filtered to remove noise. Scaffold pore size affects mechanical properties and biological cell behaviour and is a frequently assessed parameter. This paper identifies and characterizes an artefact affecting a commonly used filter which erroneously increases mean pore size. The 3D sweep despeckling filter removes all but the largest object within a volume of interest, and therefore deletes any disconnected objects located at the periphery, increasing measured mean pore size. This artefact is characterized, and effective methods to mitigate its effects are devised, involving despeckling a sufficiently large volume of interest, then reducing the volume of interest in size to remove the error prior to analysis. Techniques to effectively apply this method to other data sets are described. This method eliminates the artefact but is time-consuming and computationally expensive. Alternative, more economical filters which remove objects below a specified size are also assessed but are shown to be affected by the same artefact. These results will help to guide the implementation of future studies investigating the effects of pore size.

Influence of non-enzymatic glycation on the mechanical properties of cortical bone

Poor bone quality induced by non-enzymatic glycation (NEG) of bone tissue in patients with type 2 diabetes mellitus (T2DM) is regarded as the major factor of bone fragility and affecting bone mechanical properties. A comprehensive and systemic mechanical investigation for evaluating the effect of NEG on bone was still lacking. In order to provide additional information for the bone quality of T2DM, the effects of NEG on mechanical properties of cortical bone were investigated in terms of elastoplasticity, fracture toughness and viscoelasticity. All samples of cortical bone, including the samples of strength test (n = 20), fracture toughness test (n = 40, quasi-static and fall-like conditions with displacement rates of 10^-3 mm/s and 10 mm/s, respectively) and stress relaxation test (n = 20), were harvested from bovine tibiae. The samples of each test were equally divided into incubated-control group and ribose-incubated group. All mechanical tests were performed after incubating all samples for 15 days. Post-yield strain (p = 0.014), post-yield energy (p < 0.0001) and damage fraction (p = 0.040) of ribose-incubated group were significantly lower than those of incubated-control group, but secant modulus (p = 0.029) of ribose-incubated group was significantly higher than that of incubated-control group. In quasi-static condition, the plastic contribution J_pl of fracture toughness (p = 0.043) of ribose-incubated group was significantly lower than that of incubated-control group. In fall-like condition, there were no differences in J_pl, elastic contribution J_el and J-integral in both two groups. The quasi-static J_el (p < 0.0001, p < 0.0001) of incubated-control and ribose-incubated groups and J-integral (p = 0.007) of incubated-control group were all significantly higher than those of fall-like condition. In stress relaxation test, initial modulus E₀ (p = 0.040) and equilibrium modulus (p = 0.029) of ribose-incubated group were significantly higher than those of incubated-control group. Reductions of relaxation modulus, which were the differences between two adjacent time points within 700 s-3000 s for ribose-incubated group, were significantly lower than those of incubated-control group. NEG could decrease the post-yield properties and quasi-static facture toughness of cortical bone, especially the plastic contribution of quasi-static fracture toughness. It could also decrease the viscoelasticity of cortical bone. The present study confirmed the negative effects of NEG on the mechanical properties of cortical bone in terms of elastoplasticity, fracture toughness and viscoelasticity, but NEG had no significant effect on the fracture toughness of cortical bone at fall-like loading. These results provided more evidence for increased fragility of cortical bone in patients with T2DM.

Identifying Regional Viscoelastic Properties of the Superficial Muscular Aponeurotic System

BACKGROUND

Suspension of the superficial muscular aponeurotic system (SMAS) is generally believed to be necessary in facelift surgery. Although many techniques have been suggested, all rely on the viscoelastic properties of the SMAS.

OBJECTIVES

The aim of this study was to determine the viscoelastic properties of bursting strength, stress relaxation, and creep in the lateral, mid-cheek, and medial regions of the SMAS.

METHODS

The viscoelastic properties of the SMAS were determined in 12 cadaveric hemifaces. Lateral SMAS was classified as the SMAS overlying the parotid gland; mid-cheek SMAS as anterior to the parotid and overlying the masseter muscle; and medial SMAS as including tissue extending medial from the lateral canthus and ending at the nasolabial fold.

RESULTS

The 3 SMAS regions showed significantly different bursting strengths: 38.9 N for the lateral SMAS, 26.7 N for the mid-cheek SMAS, and 11.9 N for the medial SMAS (P < 0.0001). Stress relaxation was similar in all vertical regions with measurements of 54% in the lateral, 48% in the mid-cheek, and 59% in the medial SMAS. Creep was found to be similar in the lateral and mid-cheek SMAS with values of 18% and 19%, respectively. The medial SMAS was noted to have a higher creep at 22%.

CONCLUSIONS

The lateral SMAS has a stronger bursting strength than the mid-cheek and medial SMAS. Creep appears to be lower in the lateral and mid-cheek SMAS. Stress relaxation appears to be similar in all 3 vertical regions.

A Review on Recent Advances in the Constitutive Modeling of Bone Tissue

PURPOSE OF REVIEW

Image-based finite element analysis (FEA) to predict and understand the biomechanical response has become an essential methodology in musculoskeletal research. An important part of such simulation models is the constitutive material model of which recent advances are summarized in this review.

RECENT FINDINGS

The review shows that existing models from other fields were introduced, such as cohesion zone (cortical bone) or phase-field models (trabecular bone). Some progress has been made in describing cortical bone involving physical mechanisms such as microcracks. Problems with validations at different length scales remain a problem. The improvement of recent constitutive models is partially obscured by uncertainties that affect overall predictions, such as image quality and calibration or boundary conditions. Nevertheless, in vivo CT-based FEA simulations based on a sophisticated constitutive behavior are a very valuable tool for clinical-related osteoporosis research.

A clinically friendly viscoelastic finite element analysis model of the mandible with Herbst appliance

INTRODUCTION

As a powerful numerical approximation tool, finite element analysis (FEA) has been widely used to predict stress and strain distributions in facial bones generated by orthodontic appliances. Previous FEA models were constructed on the basis of a linear elastic phase of the bone response (eg, elastic bone strains to loading). However, what is more useful for clinical understanding would be predicting long-term strains and displacements of bone-segments responding to loading, yet tissue responses are (1) not promptly observable and (2) hard to predict in nature.

METHODS

Viscoelastic property of the mandibular bone was incorporated into FEA models to visualize long-term, time-dependent stress and strain patterns in the mandible after being exposed to orthopedic stress. A mandible under loading by a Herbst appliance was modeled, and outcomes of the constructed elastic and viscoelastic models were compared.

RESULTS

Patterns and magnitudes of the displacement throughout the mandible predicted by the viscoelastic model were exhibited in accordance with previous clinical outcomes of Herbst appliance therapy. The elastic models exhibited similar displacement patterns; however, the magnitude of the displacements in the models was invariably small (approximately 1 per 100) compared with those outputs of corresponding viscoelastic models. The corresponding maximum stress level in our viscoelastic mandible subjected to the Herbst appliance with the same loading was considerably low and relaxed in various regions when compared with the elastic model.

CONCLUSIONS

We suggest that a viscoelastic model of the mandible mimics our general prediction of orthopedic treatment outcomes better than those by elastic models.

Patient-Specific Bone Multiscale Modelling, Fracture Simulation and Risk Analysis-A Survey

This paper provides a starting point for researchers and practitioners from biology, medicine, physics and engineering who can benefit from an up-to-date literature survey on patient-specific bone fracture modelling, simulation and risk analysis. This survey hints at a framework for devising realistic patient-specific bone fracture simulations. This paper has 18 sections: Section 1 presents the main interested parties; Section 2 explains the organzation of the text; Section 3 motivates further work on patient-specific bone fracture simulation; Section 4 motivates this survey; Section 5 concerns the collection of bibliographical references; Section 6 motivates the physico-mathematical approach to bone fracture; Section 7 presents the modelling of bone as a continuum; Section 8 categorizes the surveyed literature into a continuum mechanics framework; Section 9 concerns the computational modelling of bone geometry; Section 10 concerns the estimation of bone mechanical properties; Section 11 concerns the selection of boundary conditions representative of bone trauma; Section 12 concerns bone fracture simulation; Section 13 presents the multiscale structure of bone; Section 14 concerns the multiscale mathematical modelling of bone; Section 15 concerns the experimental validation of bone fracture simulations; Section 16 concerns bone fracture risk assessment. Lastly, glossaries for symbols, acronyms, and physico-mathematical terms are provided.

Microstructure-based numerical computational method for the insertion torque of dental implant

The bone quality has a significant effect on the insertion torque of dental implant. In most clinical studies, bone density is used as a gold standard in predicting insertion torque. By contrast, trabecular microstructure is ignored. In this study, a microstructure-based numerical computational method with high accuracy and efficiency for the insertion torque of dental implant was proposed by introducing two microscopic variables, namely, volume fraction and fabric tensor. First, two kinds of 3D microstructural solid models with same volume fraction and fabric tensor were established on the basis of the microstructural topology of six reference specimens. Second, a new numerical simulation method based on homogenous theory was used to explore the material models of these 3D microstructural solid models at the microscopic scale. Then, the anisotropic material models of specimens were developed on the basis of the mixture rule. Thereafter, a numerical simulation based on the anisotropic finite element (FE) model was carried out to acquire the insertion torque. To demonstrate the efficiency and accuracy of the simulation based on the anisotropic FE model, numerical simulations based on isotropic FE model and micro-computer tomography (micro-CT) FE models were also implemented as comparisons. Comparison of the simulated peak insertion torques of the anisotropic, isotropic, and micro-CT FE models with insertion experiments demonstrated the feasibility and potential of the proposed method. The anisotropic FE model reduced the time consumption by 91.85% and enhanced the accuracy by 11.82% compared with the micro-CT and isotropic FE models, respectively.

Analysis of Trabecular Bone Mechanics Using Machine Learning

"Bone remodeling" is a dynamic process, and mutliphase analysis incorporated with the forecasting algorithm can help the biologists and orthopedics to interpret the laboratory generated results and to apply them in improving applications in the fields of "drug design, treatment, and therapy" of diseased bones. The metastasized bone microenvironment has always remained a challenging puzzle for the researchers. A multiphase computational model is interfaced with the artificial intelligence algorithm in a hybrid manner during this research. Trabecular surface remodeling is presented in this article, with the aid of video graphic footage, and the associated parametric thresholds are derived from artificial intelligence and clinical data.

Standardizing compression testing for measuring the stiffness of human bone

Objectives

The ability to determine human bone stiffness is of clinical relevance in many fields, including bone quality assessment and orthopaedic prosthesis design. Stiffness can be measured using compression testing, an experimental technique commonly used to test bone specimens in vitro. This systematic review aims to determine how best to perform compression testing of human bone.

Methods

A keyword search of all English language articles up until December 2017 of compression testing of bone was undertaken in Medline, Embase, PubMed, and Scopus databases. Studies using bulk tissue, animal tissue, whole bone, or testing techniques other than compression testing were excluded.

Results

A total of 4712 abstracts were retrieved, with 177 papers included in the analysis; 20 studies directly analyzed the compression testing technique to improve the accuracy of testing. Several influencing factors should be considered when testing bone samples in compression. These include the method of data analysis, specimen storage, specimen preparation, testing configuration, and loading protocol.

Conclusion

Compression testing is a widely used technique for measuring the stiffness of bone but there is a great deal of inter-study variation in experimental techniques across the literature. Based on best evidence from the literature, suggestions for bone compression testing are made in this review, although further studies are needed to establish standardized bone testing techniques in order to increase the comparability and reliability of bone stiffness studies.

Cite this article: S. Zhao, M. Arnold, S. Ma, R. L. Abel, J. P. Cobb, U. Hansen, O. Boughton. Standardizing compression testing for measuring the stiffness of human bone. Bone Joint Res 2018;7:524–538. DOI: 10.1302/2046-3758.78.BJR-2018-0025.R1.