The twist-to-bend compliance of the Rheum rhabarbarum petiole: integrated computations and experiments

Plant petioles can be considered as hierarchical cellular structures, displaying geometric features defined at multiple length scales. Their macroscopic mechanical properties are the cumulative outcome of structural properties attained at each level of the structural hierarchy. This work appraises the compliance of a rhubarb stalk by determining the stalk’s bending and torsional stiffness both computationally and experimentally. In our model, the irregular cross-sectional shape of the petiole and the layers of the constituent tissues are considered to evaluate the stiffness properties at the structural level. The arbitrary shape contour of the petiole is generated with reasonable accuracy by the Gielis superformula. The stiffness and architecture of the constituent layered tissues are modeled by using the concept of shape transformers so as to obtain the computational twist-to-bend ratio for the petiole. The rhubarb stalk exhibits a ratio of flexural to torsional stiffness 4.04 (computational) and 3.83 (experimental) in comparison with 1.5 for isotropic, incompressible, circular cylinders, values that demonstrate the relative structural compliance to flexure and torsion.     

掌叶大黄根茎大黄多糖的贮藏和分布特征

采用组织化学方法和苯酚硫酸比色法研究了大黄多糖在掌叶大黄(Rheum palmatum)根茎中的贮藏分布特征及其含量变化规律。结果表明: 大黄多糖在根茎中呈多位点贮藏, 在根茎的周皮、皮层、次生维管组织的薄壁细胞以及髓内不同程度地贮藏和积累了一定数量的大黄多糖, 次生木质部的木薄壁细胞以及次生韧皮部的韧皮薄壁细胞是大黄多糖的主要贮藏和积累部位; 不同发育时期根茎中大黄多糖含量的变化规律为: 随着植物的生长, 根茎及其各组织中大黄多糖的总含量表现为逐渐增高的趋势, 但在发育的后期略有下降; 与木薄壁细胞相比, 韧皮薄壁细胞贮藏大黄多糖量相对较多, 大黄多糖的贮藏和积累方式为逐渐累积。次生维管组织为大黄多糖贮藏和积累的主要组织。     

掌叶大黄根茎大黄多糖的贮藏和分布特征

采用组织化学方法和苯酚硫酸比色法研究了大黄多糖在掌叶大黄（Rheum palmatum）根茎中的贮藏分布特征及其含量变化规律。结果表明：大黄多糖在根茎中呈多位点贮藏,在根茎的周皮、皮层、次生维管组织的薄壁细胞以及髓内不同程度地贮藏和积累了一定数量的大黄多糖,次生木质部的木薄壁细胞以及次生韧皮部的韧皮薄壁细胞是大黄多糖的主要贮藏和积累部位;不同发育时期根茎中大黄多糖含量的变化规律为：随着植物的生长,根茎及其各组织中大黄多糖的总含量表现为逐渐增高的趋势,但在发育的后期略有下降;与木薄壁细胞相比,韧皮薄壁细胞贮藏大黄多糖量相对较多,大黄多糖的贮藏和积累方式为逐渐累积。次生维管组织为大黄多糖贮藏和积累的主要组织。     

Tooth Eruption Results from Bone Remodelling Driven by Bite Forces Sensed by Soft Tissue Dental Follicles: A Finite Element Analysis

Intermittent tongue, lip and cheek forces influence precise tooth position, so we here examine the possibility that tissue remodelling driven by functional bite-force-induced jaw-strain accounts for tooth eruption. Notably, although a separate true ‘eruptive force’ is widely assumed, there is little direct evidence for such a force. We constructed a three dimensional finite element model from axial computerized tomography of an 8 year old child mandible containing 12 erupted and 8 unerupted teeth. Tissues modelled included: cortical bone, cancellous bone, soft tissue dental follicle, periodontal ligament, enamel, dentine, pulp and articular cartilage. Strain and hydrostatic stress during incisive and unilateral molar bite force were modelled, with force applied via medial and lateral pterygoid, temporalis, masseter and digastric muscles. Strain was maximal in the soft tissue follicle as opposed to surrounding bone, consistent with follicle as an effective mechanosensor. Initial numerical analysis of dental follicle soft tissue overlying crowns and beneath the roots of unerupted teeth was of volume and hydrostatic stress. To numerically evaluate biological significance of differing hydrostatic stress levels normalized for variable finite element volume, ‘biological response units’ in Nmm were defined and calculated by multiplication of hydrostatic stress and volume for each finite element. Graphical representations revealed similar overall responses for individual teeth regardless if incisive or right molar bite force was studied. There was general compression in the soft tissues over crowns of most unerupted teeth, and general tension in the soft tissues beneath roots. Not conforming to this pattern were the unerupted second molars, which do not erupt at this developmental stage. Data support a new hypothesis for tooth eruption, in which the follicular soft tissues detect bite-force-induced bone-strain, and direct bone remodelling at the inner surface of the surrounding bony crypt, with the effect of enabling tooth eruption into the mouth.     

Static indentation of anisotropic biomaterials using axially asymmetric indenters--a computational study

Indentation has historically been used by biomechanicians to extract the small strain elastic or viscoelastic properties of biological tissues. Because of the axisymmetry of indenters used in these studies however, analysis of the results requires the assumption of material isotropy and often yields an "effective" elastic modulus. Since most biological tissues such as bone and myocardium are known to be anisotropic, the use of conventional indentation techniques for estimating material properties is therefore limited. The feasibility of using an axially asymmetric indenter to determine material directions and in-plane material properties for anisotropic tissue is explored here using finite element analysis. The load versus displacement curves as would be measured by an indenter depend on the orientation of the indenter cross section relative to the in-plane material axes, thus suggesting a method for determining the underlying material directions. Additionally, the stiffness of the tissue response to indentation is sensitive to the values of the in-plane anisotropic material properties and prestretches, and thus test results can be used to back out relevant constitutive parameters.     

Nonlinear incompressible finite element for simulating loading of cardiac tissue--Part I: Two dimensional formulation for thin myocardial strips

A two-dimensional incompressible plane-stress finite element is formulated for the simulation of the passive-state mechanics of thin myocardial strips. The formulation employs a total Lagrangian and materially nonlinear approach, being based on a recently proposed structural material law, which is derived from the histological composition of the tissue. The ensuing finite element allows to demonstrate the mechanical properties of a single myocardial layer containing uniformly directed fibers by simulating various loading cases such as tension, compression and shear. The results of these cases show that the fiber direction is considerably stiffer than the cross-fiber direction, that there is significant coupling between these two directions, and that the shear stiffness of the tissue is lower than its tensile and compressive stiffness.     

Anisotropy of fibrous tissues in relation to the distribution of tensed and buckled fibers

Fibrous tissues are characterized by a much higher stiffness in tension than compression. This study uses microstructural modeling to analyze the material symmetry of fibrous tissues undergoing tension and compression, to better understand how material symmetry relates to the distribution of tensed and buckled fibers. The analysis is also used to determine whether the behavior predicted from a microstructural model can be identically described by phenomenological continuum models. The analysis confirms that in the case when all the fibers are in tension in the current configuration, the material symmetry of a fibrous tissue in the corresponding reference configuration is dictated by the symmetry of its fiber angular distribution in that configuration. However, if the strain field exhibits a mix of tensile and compressive principal normal strains, the fibrous tissue is represented by a material body which consists only of those fibers which are in tension; the material symmetry of this body may be deduced from the superposition of the planes of symmetry of the strain and the planes of symmetry of the angular fiber distribution. Thus the material symmetry is dictated by the symmetry of the angular distribution of only those fibers which are in tension. Examples are provided for various fiber angular distribution symmetries. In particular, it is found that a fibrous tissue with isotropic fiber angular distribution exhibits orthotropic symmetry when subjected to a mix of tensile and compressive principal normal strains, with the planes of symmetry normal to the principal directions of the strain. This anisotropy occurs even under infinitesimal strains and is distinct from the anisotropy induced from the finite rotation of fibers. It is also noted that fibrous materials are not stable under all strain states due to the inability of fibers to sustain compression along their axis; this instability can be overcome by the incorporation of a ground matrix. It is concluded that the material response predicted using a microstructural model of the fibers cannot be described exactly by phenomenological continuum models. These results are also applicable to nonbiological fiber-composite materials.     

Comparison of the linear finite element prediction of deformation and strain of human cancellous bone to 3D digital volume correlation measurements

The mechanical properties of cancellous bone and the biological response of the tissue to mechanical loading are related to deformation and strain in the trabeculae during function. Due to the small size of trabeculae, their motion is difficult to measure. To avoid the need to measure trabecular motions during loading the finite element method has been used to estimate trabecular level mechanical deformation. This analytical approach has been empirically successful in that the analytical models are solvable and their results correlate with the macroscopically measured stiffness and strength of bones. The present work is a direct comparison of finite element predictions to measurements of the deformation and strain at near trabecular level. Using the method of digital volume correlation, we measured the deformation and calculated the strain at a resolution approaching the trabecular level for cancellous bone specimens loaded in uniaxial compression. Smoothed results from linearly elastic finite element models of the same mechanical tests were correlated to the empirical three-dimensional (3D) deformation in the direction of loading with a coefficient of determination as high as 97% and a slope of the prediction near one. However, real deformations in the directions perpendicular to the loading direction were not as well predicted by the analytical models. Our results show, that the finite element modeling of the internal deformation and strain in cancellous bone can be accurate in one direction but that this does not ensure accuracy for all deformations and strains.     

掌叶大黄根多糖的积累分布特征

采用组织化学方法和苯酚硫酸比色法研究了掌叶大黄（Rheum palmatum）根中大黄多糖的贮藏分布特征和含量变化规律。结果表明：大黄多糖在根内的贮藏是多位点的,在根周皮的栓内层、次生维管组织的薄壁细胞内不同程度地贮藏和积累了一定数量的大黄多糖,次生木质部的木薄壁细胞和次生韧皮部的韧皮薄壁细胞是主要贮藏和积累的部位;不同发育时期根中大黄多糖含量的变化规律为,随着植物的不断成熟,根及其各组织中大黄多糖的总含量表现为逐渐增高的趋势,但在发育的后期略有下降;韧皮薄壁细胞与木薄壁细胞相比,贮藏大黄多糖的含量相对较多,大黄多糖的贮藏积累方式为逐渐累积的方式。     

A large-strain finite element formulation for biological tissues with application to mitral valve leaflet tissue mechanics

This paper presents a finite element formulation suitable for large-strain modeling of biological tissues and uses this formulation to implement an accurate finite element model for mitral valve leaflet tissue. First, an experimentally derived strain energy function is obtained from literature. This function is implemented in finite elements using the mixed pressure-displacement formulation. A modification is made to aid in maintaining positive definiteness of the stiffness matrix at low strains. The numerical implementation is shown to be accurate in representing the analytical model of material behavior. The mixed formulation is useful for modeling of soft biological tissues in general, and the model presented here is applicable to finite element simulation of mitral valve mechanics.     

Microplane constitutive model and computational framework for blood vessel tissue

This paper presents a nonlinearly elastic anisotropic microplane formulation in 3D for computational constitutive modeling of arterial soft tissue in the passive regime. The constitutive modeling of arterial (and other biological) soft tissue is crucial for accurate finite element calculations, which in turn are essential for design of implants, surgical procedures, bioartificial tissue, as well as determination of effect of progressive diseases on tissues and implants. The model presented is defined at a lower scale (mesoscale) than the conventional macroscale and it incorporates the effect of all the (collagen) fibers which are anisotropic structural components distributed in all directions within the tissue material in addition to that of isotropic bulk tissue. It is shown that the proposed model not only reproduces Holzapfel's recent model but also improves on it by accounting for the actual three-dimensional distribution of fiber orientation in the arterial wall, which endows the model with advanced capabilities in simulation of remodeling of soft tissue. The formulation is flexible so that its parameters could be adjusted to represent the arterial wall either as a single material or a material composed of several layers in finite element analyses of arteries. Explicit algorithms for both the material subroutine and the explicit integration with dynamic relaxation of equations of motion using finite element method are given. To circumvent the slow convergence of the standard dynamic relaxation and small time steps dictated by the stability of the explicit integrator, an adaptive dynamic relaxation technique that ensures stability and fastest possible convergence rates is developed. Incompressibility is enforced using penalty method with an updated penalty parameter. The model is used to simulate experimental data from the literature demonstrating that the model response is in excellent agreement with the data. An experimental procedure to determine the distribution of fiber directions in 3D for biological soft tissue is suggested in accordance with the microplane concept. It is also argued that this microplane formulation could be modified or extended to model many other phenomena of interest in biomechanics.     

Trabecular bone's mechanical properties are affected by its non-uniform mineral distribution

The bone remodeling process takes place at the surface of trabeculae and results in a non-uniform mineral distribution. This will affect the mechanical properties of cancellous bone, because the properties of bone tissue depend on its mineral content. We investigated how large this effect is by simulating several non-uniform mineral distributions in 3D finite element models of human trabecular bone and calculating the apparent stiffness of these models. In the ‘linear model’ we assumed a linear relation between mineral content and Young's modulus of the tissue. In the ‘exponential model’ we included an empirical exponential relation in the model. When the linear model was used the mineral distribution slightly changed the apparent stiffness, the difference varied between an 8% decrease and a 4% increase compared to the uniform model with the same BMD. The exponential model resulted in up to 20% increased apparent stiffness in the main load-bearing direction. A thin less mineralized surface layer (28 μm) and highly mineralized interstitial bone (mimicking mineralization resulting from anti-resorptive treatment) resulted in the highest stiffness. This could explain large reductions in fracture risk resulting from small increases in BMD. The non-uniform mineral distribution could also explain why bone tissue stiffness determined using nano-indentation is usually higher than finite element (FE)-determined stiffness. We conclude that the non-uniform mineral distribution in trabeculae does affect the mechanical properties of cancellous bone and that the tissue stiffness determined using FE-modeling could be improved by including detailed information about mineral distribution in trabeculae in the models.     

Multiscale mechanobiology: Coupling models of adhesion kinetics and nonlinear tissue mechanics

The mechanical behavior of tissues at the macroscale is tightly coupled to cellular activity at the microscale. Dermal wound healing is a prominent example of a complex system in which multiscale mechanics regulate restoration of tissue form and function. In cutaneous wound healing, a fibrin matrix is populated by fibroblasts migrating in from a surrounding tissue made mostly out of collagen. Fibroblasts both respond to mechanical cues, such as fiber alignment and stiffness, as well as exert active stresses needed for wound closure.Here, we develop a multiscale model with a two-way coupling between a microscale cell adhesion model and a macroscale tissue mechanics model. Starting from the well-known model of adhesion kinetics proposed by Bell, we extend the formulation to account for nonlinear mechanics of fibrin and collagen and show how this nonlinear response naturally captures stretch-driven mechanosensing. We then embed the new nonlinear adhesion model into a custom finite element implementation of tissue mechanical equilibrium. Strains and stresses at the tissue level are coupled with the solution of the microscale adhesion model at each integration point of the finite element mesh. In addition, solution of the adhesion model is coupled with the active contractile stress of the cell population. The multiscale model successfully captures the mechanical response of biopolymer fibers and gels, contractile stresses generated by fibroblasts, and stress-strain contours observed during wound healing. We anticipate that this framework will not only increase our understanding of how mechanical cues guide cellular behavior in cutaneous wound healing, but will also be helpful in the study of mechanobiology, growth, and remodeling in other tissues.     

Three-dimensional finite element analysis of the foot during standing--a material sensitivity study

Information on the internal stresses/strains in the human foot and the pressure distribution at the plantar support interface under loading is useful in enhancing knowledge on the biomechanics of the ankle-foot complex. While techniques for plantar pressure measurements are well established, direct measurement of the internal stresses/strains is difficult. A three-dimensional (3D) finite element model of the human foot and ankle was developed using the actual geometry of the foot skeleton and soft tissues, which were obtained from 3D reconstruction of MR images. Except the phalanges that were fused, the interaction among the metatarsals, cuneiforms, cuboid, navicular, talus, calcaneus, tibia and fibula were defined as contact surfaces, which allow relative articulating movement. The plantar fascia and 72 major ligaments were simulated using tension-only truss elements by connecting the corresponding attachment points on the bone surfaces. The bony and ligamentous structures were embedded in a volume of soft tissues. The encapsulated soft tissue was defined as hyperelastic, while the bony and ligamentous structures were assumed to be linearly elastic. The effects of soft tissue stiffening on the stress distribution of the plantar surface and bony structures during balanced standing were investigated. Increases of soft tissue stiffness from 2 and up to 5 times the normal values were used to approximate the pathologically stiffened tissue behaviour with increasing stages of diabetic neuropathy. The results showed that a five-fold increase in soft tissue stiffness led to about 35% and 33% increase in the peak plantar pressure at the forefoot and rearfoot regions, respectively. This corresponded to about 47% decrease in the total contact area between the plantar foot and the horizontal support surface. Peak bone stress was found at the third metatarsal in all calculated cases with a minimal increase of about 7% with soft tissue stiffening.     

Cell mechanics studied by a reconstituted model tissue

Tissue models reconstituted from cells and extracellular matrix (ECM) simulate natural tissues. Cytoskeletal and matrix proteins govern the force exerted by a tissue and its stiffness. Cells regulate cytoskeletal structure and remodel ECM to produce mechanical changes during tissue development and wound healing. Characterization and control of mechanical properties of reconstituted tissues are essential for tissue engineering applications. We have quantitatively characterized mechanical properties of connective tissue models, fibroblast-populated matrices (FPMs), via uniaxial stretch measurements. FPMs resemble natural tissues in their exponential dependence of stress on strain and linear dependence of stiffness on force at a given strain. Activating cellular contractile forces by calf serum and disrupting F-actin by cytochalasin D yield "active" and "passive" components, which respectively emphasize cellular and matrix mechanical contributions. The strain-dependent stress and elastic modulus of the active component were independent of cell density above a threshold density. The same quantities for the passive component increased with cell number due to compression and reorganization of the matrix by the cells.     

A fibre reorientation model for orthotropic multiplicative growth

The main goal of this contribution consists in the development of a remodelling framework for orthotropic continua whereby the underlying symmetry group is incorporated via two fibre families. Special emphasis is placed on the modelling of biological tissues at finite deformations. Besides the incorporation of a referential mass source, anisotropic growth is addressed by means of a multiplicative decomposition of the overall deformation gradient into an elastic and a growth distortion. Projected quantities of a configurational growth stress tensor are advocated as driving forces for time-dependent saturation–type evolution of the principal values of the growth distortion. Moreover, the reorientation of both fibre families, which directly affects the strain energy as well as the growth distortion itself, is guided by analyzing critical energy points. In particular, a time-dependent formulation is developed which aligns the fibre directions according to the principal stretch directions. Finally, the proposed framework is embedded into a finite element context so that representative numerical examples, examining growth and resorption in volume and density together with fibre reorientation, close this study.     

A transversely isotropic viscoelastic constitutive equation for brainstem undergoing finite deformation

The objective of this study was to define the constitutive response of brainstem undergoing finite shear deformation. Brainstem was characterized as a transversely isotropic viscoelastic material and the material model was formulated for numerical implementation. Model parameters were fit to shear data obtained in porcine brainstem specimens undergoing finite shear deformation in three directions: parallel, perpendicular, and cross sectional to axonal fiber orientation and determined using a combined approach of finite element analysis (FEA) and a genetic algorithm (GA) optimizing method. The average initial shear modulus of brainstem matrix of 4-week old pigs was 12.7 Pa, and therefore the brainstem offers little resistance to large shear deformations in the parallel or perpendicular directions, due to the dominant contribution of the matrix in these directions. The fiber reinforcement stiffness was 121.2 Pa, indicating that brainstem is anisotropic and that axonal fibers have an important role in the cross-sectional direction. The first two leading relative shear relaxation moduli were 0.8973 and 0.0741, respectively, with corresponding characteristic times of 0.0047 and 1.4538 s, respectively, implying rapid relaxation of shear stresses. The developed material model and parameter estimation technique are likely to find broad applications in neural and orthopaedic tissues.     

Muscle-driven finite element simulation of human foot movements

This paper describes a finite element scheme for realistic muscle-driven simulation of human foot movements. The scheme is used to simulate human ankle plantar flexion. A three-dimensional anatomically detailed finite element model of human foot and lower leg is developed and the idea of generating natural foot movement based entirely on the contraction of the plantar flexor muscles is used. The bones, ligaments, articular cartilage, muscles, tendons, as well as the rest soft tissues of human foot and lower leg are included in the model. A realistic three-dimensional continuum constitutive model that describes the biomechanical behaviour of muscles and tendons is used. Both the active and passive properties of muscle tissue are accounted for. The materials for bones and ligaments are considered as homogeneous, isotropic and linearly elastic, whereas the articular cartilage and the rest soft tissues (mainly fat) are defined as hyperelastic materials. The model is used to estimate muscle tissue deformations as well as stresses and strains that develop in the lower leg muscles during plantar flexion of the ankle. Stresses and strains that develop in Achilles tendon during such a movement are also investigated.