Hyperelastic constitutive relationship for the strain-rate dependent behavior of shoulder and other joint cartilages

Non-linear finite deformations of articular cartilages under physiological loading conditions can be attributed to hyperelastic behavior. This paper contains experimental results of indentation tests in finite deformation and proposes an empirical based new generalized hyperelastic constitutive model to account for strain-rate dependency for humeral head cartilage tissues. The generalized model is based on existing hyperelastic constitutive relationships that are extensively used to represent biological tissues in biomechanical literature. The experimental results were obtained for three loading velocities, corresponding to low (1x10-3 s-1), moderate and high strain-rates (1x10-1 s-1), which represent physiological loading rates that are experienced in daily activities such as lifting, holding objects and sporting activities. Hyperelastic material parameters were identified by non linear curve fitting procedure. Analysis demonstrated that the material behavior of cartilage can be effectively decoupled into strain-rate independent(elastic) and dependent parts. Further, experiments conducted using different indenters indicated that the parameters obtained are significantly affected by the indenter size, potentially due to structural inhomogeneity of the tissue. The hyperelastic constitutive model developed in this paper opens a new avenue for the exploration of material properties of cartilage tissues.

Local stress-strain distribution and load transfer across cartilage matrix at micro-scale using combined microscopy-based finite element method

Articular cartilage is the load-bearing tissue that consists of proteoglycan macromolecules entrapped between collagen fibrils in a three-dimensional architecture. To date, the drudgery of searching for mathematical models to represent the biomechanics of such a system continues without providing a fitting description of its functional response to load at micro-scale level. We believe that the major complication arose when cartilage was first envisaged as a multiphasic model with distinguishable components and that quantifying those and searching for the laws that govern their interaction is inadequate. To the thesis of this paper, cartilage as a bulk is as much continuum as is the response of its components to the external stimuli. For this reason, we framed the fundamental question as to what would be the mechano-structural functionality of such a system in the total absence of one of its key constituents-proteoglycans. To answer this, hydrated normal and proteoglycan depleted samples were tested under confined compression while finite element models were reproduced, for the first time, based on the structural microarchitecture of the cross-sectional profile of the matrices. These micro-porous in silico models served as virtual transducers to produce an internal noninvasive probing mechanism beyond experimental capabilities to render the matrices micromechanics and several others properties like permeability, orientation etc. The results demonstrated that load transfer was closely related to the microarchitecture of the hyperelastic models that represent solid skeleton stress and fluid response based on the state of the collagen network with and without the swollen proteoglycans. In other words, the stress gradient during deformation was a function of the structural pattern of the network and acted in concert with the position-dependent compositional state of the matrix. This reveals that the interaction between indistinguishable components in real cartilage is superimposed by its microarchitectural state which directly influences macromechanical behavior.

Nonlinear pre-fire and post-fire analysis of steel frames

A numerical procedure based on the plastic hinge concept for study of the structural behaviour of steel framed structures exposed to fire is described. Most previous research on fire analysis considered the structural performance due to rising temperature. When strain reversal occurs during the cooling phase, the stress–strain curve is different. The plastic deformation is incorporated into the stress–strain curve to model the strain reversal effect in which unloading under elastic behaviour is allowed. This unloading response is traced by the incremental–iterative Newton–Raphson method. The mechanical properties of the steel member in the present fire analysis follows both Eurocode 3 Part 1.2 and BS5950 Part 8, which implicitly allow for thermal creep deformation. This paper presents an efficient fire analysis procedure for predicting thermal and cooling effects on an isolated element and a multi-storey frame. Several numerical and experimental examples related to structural behaviour in cooling phase are studied and compared with results obtained by other researchers. The proposed method is effective in the fire safety design and analysis of a building in a real fire scenario. The scope of investigation is of great significance since a large number of rescuers would normally enter a fire site as soon as the fire is extinguished and during the cooling phase, so a structural collapse can be catastrophic.

Performance-based structural fire design of steel frames using conventional computer software

Fire incident in buildings is common, so the fire safety design of the framed structure is imperative, especially for the unprotected or partly protected bare steel frames. However, software for structural fire analysis is not widely available. As a result, the performance-based structural fire design is urged on the basis of using user-friendly and conventional nonlinear computer analysis programs so that engineers do not need to acquire new structural analysis software for structural fire analysis and design. The tool is desired to have the capacity of simulating the different fire scenarios and associated detrimental effects efficiently, which includes second-order P-D and P-d effects and material yielding. Also the nonlinear behaviour of large-scale structure becomes complicated when under fire, and thus its simulation relies on an efficient and effective numerical analysis to cope with intricate nonlinear effects due to fire. To this end, the present fire study utilizes a second order elastic/plastic analysis software NIDA to predict structural behaviour of bare steel framed structures at elevated temperatures. This fire study considers thermal expansion and material degradation due to heating. Degradation of material strength with increasing temperature is included by a set of temperature-stress-strain curves according to BS5950 Part 8 mainly, which implicitly allows for creep deformation. This finite element stiffness formulation of beam-column elements is derived from the fifth-order PEP element which facilitates the computer modeling by one member per element. The Newton-Raphson method is used in the nonlinear solution procedure in order to trace the nonlinear equilibrium path at specified elevated temperatures. Several numerical and experimental verifications of framed structures are presented and compared against solutions in literature. The proposed method permits engineers to adopt the performance-based structural fire analysis and design using typical second-order nonlinear structural analysis software.

Strain engineering of selective chemical adsorption on monolayer MoS2

Nanomaterials are prone to influence by chemical adsorption because of their large surface to volume ratios. This enables sensitive detection of adsorbed chemical species which, in turn, can tune the property of the host material. Recent studies discovered that single and multi-layer molybdenum disulfide (MoS2) films are ultra-sensitive to several important environmental molecules. Here we report new findings from ab inito calculations that reveal substantially enhanced adsorption of NO and NH3 on strained monolayer MoS2 with significant impact on the properties of the adsorbates and the MoS2 layer. The magnetic moment of adsorbed NO can be tuned between 0 and 1 μB; strain also induces an electronic phase transition between half-metal and metal. Adsorption of NH3 weakens the MoS2 layer considerably, which explains the large discrepancy between the experimentally measured strength and breaking strain of MoS2 films and previous theoretical predictions. On the other hand, adsorption of NO2, CO, and CO2 is insensitive to the strain condition in the MoS2 layer. This contrasting behavior allows sensitive strain engineering of selective chemical adsorption on MoS2 with effective tuning of mechanical, electronic, and magnetic properties. These results suggest new design strategies for constructing MoS2-based ultrahigh-sensitivity nanoscale sensors and electromechanical devices.

Exploration of mechanisms underlying the strain-rate-dependent mechanical property of single chondrocytes

Based on the characterization by Atomic Force Microscopy (AFM), we report that the mechanical property of single chondrocytes has dependency on the strain-rates. By comparing the mechanical deformation responses and the Young’s moduli of living and fixed chondrocytes at four different strain-rates, we explore the deformation mechanisms underlying this dependency property. We found that the strain-rate-dependent mechanical property of living cells is governed by both of the cellular cytoskeleton (CSK) and the intracellular fluid when the fixed chondrocytes is mainly governed by their intracellular fluid which is called the consolidation-dependent deformation behavior. Finally, we report that the porohyperelastic (PHE) constitutive material model which can capture the consolidation-dependent behavior of both living and fixed chondrocytes is a potential candidature to study living cell biomechanics.

Effect of constraint on void growth near a blunt crack tip

The growth of a single cylindrical hole ahead of a blunt crack tip was studied using large deformation finite element analysis in three-point bend specimens with different precrack depth. The effect of small second phase particles was taken into account by incorporating Gurson’s constitutive equation. The effects of strain hardening and the initial distance from the hole to the crack tip were also investigated. The results show that the variation of crack tip opening displacement with load is not sensitive to constraint level. The effects of constraint on the growth of hole and ductile initiation toughness are diminished with decreasing initial distance from the hole to the blunt crack tip.

Numerical investigation on stable crack growth in plane stress

Large deformation finite element analysis has been carried out to investigate the stress-strain fields ahead of a growing crack for compact tension .a=W D 0:5/ and three-point bend .a=W D 0:1 and 0:5/ specimens under plane stress condition. The crack growth is controlled by the experimental J -integral resistance curves measured by Sun et al. The results indicate that the distributions of opening stress, equivalent stress and equivalent strain ahead of a growing crack are not sensitive to specimen geometry. For both stationary and growing cracks, similar distributions of opening stress and triaxiality can be found along the ligament. During stable crack growth, the crack-tip opening displacement (CTOD) resistance curve and the cohesive fracture energy in the fracture process zone are independent of specimen geometry and may be suitable criteria for characterizing stable crack growth in plane stress.

Swelling induced finite strain flexure in a rectangular block of an isotropic elastic material

The deformation of a rectangular block into an annular wedge is studied with respect to the state of swelling interior to the block. Nonuniform swelling fields are shown to generate these flexure deformations in the absence of resultant forces and bending moments. Analytical expressions for the deformation fields demonstrate these effects for both incompressible and compressible generalizations of conventional hyperelastic materials. Existing results in the absence of a swelling agent are recovered as special cases.

Determination of strain-rate-dependent mechanical behavior of living and fixed osteocytes and chondrocytes using atomic force microscopy and inverse finite element analysis

The aim of this paper is to determine the strain-rate-dependent mechanical behavior of living and fixed osteocytes and chondrocytes, in vitro. Firstly, Atomic Force Microscopy (AFM) was used to obtain the force-indentation curves of these single cells at four different strain-rates. These results were then employed in inverse finite element analysis (FEA) using Modified Standard neo-Hookean Solid (MSnHS) idealization of these cells to determine their mechanical properties. In addition, a FEA model with a newly developed spring element was employed to accurately simulate AFM evaluation in this study. We report that both cytoskeleton (CSK) and intracellular fluid govern the strain-rate-dependent mechanical property of living cells whereas intracellular fluid plays a predominant role on fixed cells’ behavior. In addition, through the comparisons, it can be concluded that osteocytes are stiffer than chondrocytes at all strain-rates tested indicating that the cells could be the biomarker of their tissue origin. Finally, we report that MSnHS is able to capture the strain-rate-dependent mechanical behavior of osteocyte and chondrocyte for both living and fixed cells. Therefore, we concluded that the MSnHS is a good model for exploration of mechanical deformation responses of single osteocytes and chondrocytes. This study could open a new avenue for analysis of mechanical behavior of osteocytes and chondrocytes as well as other similar types of cells.

Temperature and strain-rate dependent fracture strength of graphynes

Graphyne is an allotrope of graphene. The mechanical properties of graphynes (α-, β-, γ- and 6,6,12-graphynes) under uniaxial tension deformation at different temperatures and strain rates are studied using molecular dynamics simulations. It is found that graphynes are more sensitive to temperature changes than graphene in terms of fracture strength and Young's modulus. The temperature sensitivity of the different graphynes is proportionally related to the percentage of acetylenic linkages in their structures, with the α-graphyne (having 100% of acetylenic linkages) being most sensitive to temperature. For the same graphyne, temperature exerts a more pronounced effect on the Young's modulus than fracture strength, which is different from that of graphene. The mechanical properties of graphynes are also sensitive to strain rate, in particular at higher temperatures.

Thermal conductivity of configurable two-dimensional carbon nanotube architecture and strain modulation

We reported the thermal conductivity of the two-dimensional carbon nanotube (CNT)-based architecture, which can be constructed through welding of single-wall CNTs by electron beam. Using large-scale nonequilibrium molecular dynamics simulations, the thermal conductivity is found to vary with different junction types due to their different phonon scatterings at the junction. The strong length and strain dependence of the thermal conductivity suggests an effective avenue to tune the thermal transport properties of the CNT-based architecture, benefiting the design of nanoscale thermal rectifiers or phonon engineering.

A nonconservative Lagrangian framework for statistical fluid registration - SAFIRA

In this paper, we used a nonconservative Lagrangian mechanics approach to formulate a new statistical algorithm for fluid registration of 3-D brain images. This algorithm is named SAFIRA, acronym for statistically-assisted fluid image registration algorithm. A nonstatistical version of this algorithm was implemented, where the deformation was regularized by penalizing deviations from a zero rate of strain. In, the terms regularizing the deformation included the covariance of the deformation matrices Σ and the vector fields (q). Here, we used a Lagrangian framework to reformulate this algorithm, showing that the regularizing terms essentially allow nonconservative work to occur during the flow. Given 3-D brain images from a group of subjects, vector fields and their corresponding deformation matrices are computed in a first round of registrations using the nonstatistical implementation. Covariance matrices for both the deformation matrices and the vector fields are then obtained and incorporated (separately or jointly) in the nonconservative terms, creating four versions of SAFIRA. We evaluated and compared our algorithms' performance on 92 3-D brain scans from healthy monozygotic and dizygotic twins; 2-D validations are also shown for corpus callosum shapes delineated at midline in the same subjects. After preliminary tests to demonstrate each method, we compared their detection power using tensor-based morphometry (TBM), a technique to analyze local volumetric differences in brain structure. We compared the accuracy of each algorithm variant using various statistical metrics derived from the images and deformation fields. All these tests were also run with a traditional fluid method, which has been quite widely used in TBM studies. The versions incorporating vector-based empirical statistics on brain variation were consistently more accurate than their counterparts, when used for automated volumetric quantification in new brain images. This suggests the advantages of this approach for large-scale neuroimaging studies.

Best individual template selection from deformation tensor minimization

We study the influence of the choice of template in tensor-based morphometry. Using 3D brain MR images from 10 monozygotic twin pairs, we defined a tensor-based distance in the log-Euclidean framework [1] between each image pair in the study. Relative to this metric, twin pairs were found to be closer to each other on average than random pairings, consistent with evidence that brain structure is under strong genetic control. We also computed the intraclass correlation and associated permutation p-value at each voxel for the determinant of the Jacobian matrix of the transformation. The cumulative distribution function (cdf) of the p-values was found at each voxel for each of the templates and compared to the null distribution. Surprisingly, there was very little difference between CDFs of statistics computed from analyses using different templates. As the brain with least log-Euclidean deformation cost, the mean template defined here avoids the blurring caused by creating a synthetic image from a population, and when selected from a large population, avoids bias by being geometrically centered, in a metric that is sensitive enough to anatomical similarity that it can even detect genetic affinity among anatomies.

Formation of highly conductive composite coatings and their applications to broadband antennas and mechanical transducers

Tight networks of interwoven carbon nanotube bundles are formed in our highly conductive composite. The composite possesses propertiessuggesting a two-dimensional percolative network rather than other reported dispersions displaying three-dimensional networks. Binding nanotubes into large but tight bundles dramatically alters the morphology and electronic transport dynamics of the composite. This enables itto carry higher levels of charge in the macroscale leading to conductivities as high as 1600 S/cm. We now discuss in further detail, the electronic and physical properties of the nanotube composites through Raman spectroscopy and transmission electron microscopy analysis. When controlled and usedappropriately, the interesting properties of these composites reveal their potential for practical device applications. For instance, we used this composite to fabricate coatings, whic improve the properties of an electromagnetic antenna/amplifier transducer. The resulting transducer possesses a broadband range up to GHz frequencies. A strain gauge transducer was also fabricated using changes in conductivity to monitor structural deformations in the composite coatings.