An anisotropic, hyperelastic model for skin: experimental measurements, finite element modelling and identification of parameters for human and murine skin

Insights and mechanics-driven modeling of human cutaneous impact injuries

Cutaneous damage mechanisms related to dynamic fragment impacts are dependent on the impact angle, impact energy, and fragment characteristics including shape, volume, contact friction, and orientation. Understanding the cutaneous injury mechanism and its relationship to the fragment parameters is lacking compromising damage classification, treatment, and protection. Here we develop a high-fidelity dynamic mechanics-driven model for partial-thickness skin injuries and demonstrate the influence of fragment parameters on the injury mechanism and damage sequence. The model quantitatively predicts the wound shape, area, and depth into the skin layers for selected impact angles, kinetic energy density, and the fragment projectile type including shape and material. The detailed sequence of impact damage including epidermal tearing that occurs ahead of the fragments initial contact location, subsequent stripping of the epidermal/dermal junction, and crushing of the underlying dermis are revealed. We demonstrate that the fragment contact friction with skin plays a key role in redistributing impact energy affecting the extent of epidermal tearing and dermal crushing. Furthermore, projectile edges markedly affect injury severity dependent on the orientation of the edge during initial impact. The model provides a quantitative framework for understanding the detailed mechanisms of cutaneous damage and a basis for the design of protective equipment.

Biomechanical Response of Head Surrogate With and Without the Helmet

Measurements of brain deformations under injurious loading scenarios are actively sought. In this work, we report experimentally measured head kinematics and corresponding dynamic, two-dimensional brain simulant deformations in head surrogates under a blunt impact, with and without a helmet. Head surrogates used in this work consisted of skin, skull, dura, falx, tentorium, and brain stimulants. The head surrogate geometry was based on the global human body models consortium's head model. A base head surrogate consisting of skin-skull-brain was considered. In addition, the response of two other head surrogates, skin-skull-dura-brain, and skin-skull-dura-brain-falx-tentorium, was investigated. Head surrogate response was studied for sagittal and coronal plane rotations for impactor velocities of 1 and 3 m/s. Response of head surrogates was compared against strain measurements in PMHS. The strain pattern in the brain simulant was heterogenous, and peak strains were established within ∼30 ms. The choice of head surrogate affect the spatiotemporal evolution of strain. For no helmet case, peak MPS of ∼50-60% and peak MSS of ∼35-50% were seen in brain simulant corresponding to peak rotational accelerations of ∼5000-7000 rad/s2. Peak head kinematics and peak MPS have been reduced by up to 75% and 45%, respectively, with the conventional helmet and by up to 90% and 85%, respectively, with the helmet with antirotational pads. Overall, these results provide important, new data on brain simulant strains under a variety of loading scenarios-with and without the helmets.

Impact of Indenter Size and Microrelief Anisotropy on the Tribological Behavior of Human Skin

Everyday, we interact with screens, sensors, and many other devices through contact with the skin. Experimental efforts have increased our knowledge of skin tribology but are challenged by the fact that skin has a complex structure, undergoes finite deformations, has nonlinear material response, and has properties that vary with anatomical location, age, sex, and environmental conditions. Computational models are powerful tools to dissect the individual contribution of these variables to the overall frictional response. Here, we present a three-dimensional high-fidelity multilayer skin computational model including a detailed surface topography or skin microrelief. Four variables are explored: local coefficient of friction (COF), indenter size, mechanical properties of the stratum corneum, and displacement direction. The results indicate that the global COF depends nonlinearly on the local COF, implying a role for skin deformation on the friction response. The global COF is also influenced by the ratio of the indenter size to the microrelief features, with larger indenters smoothing out the role of skin topography. Changes in stiffness of the uppermost layer of skin associated with humidity have a substantial effect on both the contact area and the reaction forces, but the overall changes in the COF are small. Finally, for the microrelief tested, the response can be considered isotropic. We anticipate that this model and results will enable the design of materials and devices for a desired interaction against skin.

Sensory neuron activation from topical treatments modulates the sensorial perception of human skin

Neural signaling of skin sensory perception from topical treatments is often reported in subjective terms such as a sensation of skin "tightness" after using a cleanser or "softness" after applying a moisturizer. However, the mechanism whereby cutaneous mechanoreceptors and corresponding sensory neurons are activated giving rise to these perceptions has not been established. Here, we provide a quantitative approach that couples in vitro biomechanical testing and detailed computational neural stimulation modeling along with a comprehensive in vivo self-assessment survey to demonstrate how cutaneous biomechanical changes in response to treatments are involved in the sensorial perception of the human skin. Strong correlations are identified between reported perception up to 12 hours post treatment and changes in the computed neural stimulation from mechanoreceptors residing deep under the skin surface. The study reveals a quantitative framework for understanding the biomechanical neural activation mechanism and the subjective perception by individuals.

Poroelastic behavior and water permeability of human skin at the nanoscale

Topical skin care products and hydrating compositions (moisturizers or injectable fillers) have been used for years to improve the appearance of, for example facial wrinkles, or to increase "plumpness". Most of the studies have addressed these changes based on the overall mechanical changes associated with an increase in hydration state. However, little is known about the water mobility contribution to these changes as well as the consequences to the specific skin layers. This is important as the biophysical properties and the biochemical composition of normal stratum corneum, epithelium, and dermis vary tremendously from one another. Our current studies and results reported here have focused on a novel approach (dynamic atomic force microscopy-based nanoindentation) to quantify biophysical characteristics of individual layers of ex vivo human skin. We have discovered that our new methods are highly sensitive to the mechanical properties of individual skin layers, as well as their hydration properties. Furthermore, our methods can assess the ability of these individual layers to respond to both compressive and shear deformations. In addition, since human skin is mechanically loaded over a wide range of deformation rates (frequencies), we studied the biophysical properties of skin over a wide frequency range. The poroelasticity model used helps to quantify the hydraulic permeability of the skin layers, providing an innovative method to evaluate and interpret the impact of hydrating compositions on water mobility of these different skin layers.

The assessment of natural human skin tension orientation and its variation according to age for two body areas: Forearm and thigh

Human skin has a complex multilayer structure consisting of non-homogeneous, non-linear, viscoelastic and anisotropic materials subjected to in vivo natural pre-tension. This natural tension stems from networks of collagen and elastin fibers. The 3D organization of the collagen and elastin fibers underpins the multidirectional natural tensions in the skin volume while the state of the networks formed influences the surface topography of the skin. This topography depends on the area of the body and on the age of the person. Experiments reported in the literature have been performed ex vivo or on cadavers. By contrast, this work proposes the characterization of the anisotropic natural tension of the human skin in vivo. Experimental tests were performed on the forearms and thighs of 42 female volunteers representing two age groups [20 - 30] and [45-55] years old. Non-contact impact tests and skin-folding tests were conducted using devices developed at the LTDS (Lyon, France). The impact test generated a Rayleigh wave that spread in the skin. The speed of this wave was measured in 7 directions to study the anisotropy of the skin tension. The image of the skin relief at rest and during the skin folding test was reconstructed by optical confocal microscopy and provided the density of the skin lines printed on the outer surface of the skin. Skin folding test enables the clinician's manual procedure to be instrumented to identify tension lines i.e., Langer lines, for better healing during a surgical procedure. The main directions of natural skin tension deduced from the measured wave speed and the densities of skin lines were [40°-60°] for the forearm and [0°-20°] for the thigh, considering that the longitudinal axis of the body is situated at 90° and the transversal axis at 0°. This method shows the remarkable effect of age and body area on the mechanical behavior of human skin in vivo. The elastic properties and natural tension of the skin decrease with age. This decrease is greater in the directions orthogonal to the skin's tension lines, leading to the accentuation of the anisotropic behavior of the cutaneous tissue. The main direction of skin tension is highly dependent on the area of the body and is directed towards a preferred direction which corresponds to the main direction of skin tension.

Mechano-biological and bio-mechanical pathways in cutaneous wound healing

Injuries to the skin heal through coordinated action of fibroblast-mediated extracellular matrix (ECM) deposition, ECM remodeling, and wound contraction. Defects involving the dermis result in fibrotic scars featuring increased stiffness and altered collagen content and organization. Although computational models are crucial to unravel the underlying biochemical and biophysical mechanisms, simulations of the evolving wound biomechanics are seldom benchmarked against measurements. Here, we leverage recent quantifications of local tissue stiffness in murine wounds to refine a previously-proposed systems-mechanobiological finite-element model. Fibroblasts are considered as the main cell type involved in ECM remodeling and wound contraction. Tissue rebuilding is coordinated by the release and diffusion of a cytokine wave, e.g. TGF-β, itself developed in response to an earlier inflammatory signal triggered by platelet aggregation. We calibrate a model of the evolving wound biomechanics through a custom-developed hierarchical Bayesian inverse analysis procedure. Further calibration is based on published biochemical and morphological murine wound healing data over a 21-day healing period. The calibrated model recapitulates the temporal evolution of: inflammatory signal, fibroblast infiltration, collagen buildup, and wound contraction. Moreover, it enables in silico hypothesis testing, which we explore by: (i) quantifying the alteration of wound contraction profiles corresponding to the measured variability in local wound stiffness; (ii) proposing alternative constitutive links connecting the dynamics of the biochemical fields to the evolving mechanical properties; (iii) discussing the plausibility of a stretch- vs. stiffness-mediated mechanobiological coupling. Ultimately, our model challenges the current understanding of wound biomechanics and mechanobiology, beside offering a versatile tool to explore and eventually control scar fibrosis after injury.

A review of the characterizations of soft tissues used in human body modeling: Scope, limitations, and the path forward

Soft tissue material properties are vital to human body models that evaluate interactions between the human body and its environment. Such models evaluate internal stress/strain responses in soft tissues to investigate issues like pressure injuries. Numerous constitutive models and parameters have been used to represent mechanical behavior of soft tissues in biomechanical models under quasi-static loading. However, researchers reported that generic material properties cannot accurately represent specific target populations due to large inter-individual variability. Two challenges that exist are experimental mechanical characterization and constitutive modeling of biological soft tissues and personalization of constitutive parameters using non-invasive, non-destructive bedside testing methods. It is imperative to understand the scope and appropriate applications for reported material properties. Thus, the goal of this paper was to compile studies from which soft tissue material properties were obtained and categorize them by source of tissue samples, methods used to quantify deformation, and material models used to describe tissues. The collected studies displayed wide ranges of material properties, and factors that affected the properties included whether tissue samples were in vivo or ex vivo, from humans or animals, the body region tested, body position during in vivo studies, deformation measurements, and material models used to describe tissues. Because of the factors that affected reported material properties, it is clear that much progress has been made in understanding soft tissue responses to loading, yet there is a need to broaden the scope of reported soft tissue material properties and better match reported properties to appropriate human body models.

A Micro-Scale Non-Linear Finite Element Model to Optimize the Mechanical Behavior of Bioprinted Constructs

Extrusion-based bioprinting is an enabling biofabrication technique that is used to create heterogeneous tissue constructs according to patient-specific geometries and compositions. The optimization of bioinks as per requirements for specific tissue applications is an essential exercise in ensuring clinical translation of the bioprinting technologies. Most notably, optimum hydrogel polymer concentrations are required to ensure adequate mechanical properties of bioprinted constructs without causing significant shear stresses on cells. However, experimental iterations are often tedious for optimizing the bioink properties. In this work, a nonlinear finite element modeling approach has been undertaken to determine the effect of different bioink parameters such as composition, concentration on the range of stresses being experienced by the cells in the bioprinted construct. The stress distribution of the cells at different parts of the constructs has also been modeled. It is found that both bioink chemical compositions and concentrations can substantially alter the stress effects experienced by the cells. Concentrated regions of softer cells near pore regions were found to increase stress concentrations by almost three times compared with stress generated in cells away from the pores. The study provides a method for rapid optimization of bioinks, design of bioprinted constructs, as well as toolpath plans for fabricating constructs with homogenous properties.

Bayesian Inference With Gaussian Process Surrogates to Characterize Anisotropic Mechanical Properties of Skin From Suction Tests

One of the intrinsic features of skin and other biological tissues is the high variation in the mechanical properties across individuals and different demographics. Mechanical characterization of skin is still a challenge because the need for subject-specific in vivo parameters prevents us from utilizing traditional methods, e.g., uniaxial tensile test. Suction devices have been suggested as the best candidate to acquire mechanical properties of skin noninvasively, but capturing anisotropic properties using a circular probe opening-which is the conventional suction device-is not possible. On the other hand, noncircular probe openings can drive different deformations with respect to fiber orientation and therefore could be used to characterize the anisotropic mechanics of skin noninvasively. We propose the use of elliptical probe openings and a methodology to solve the inverse problem of finding mechanical properties from suction measurements. The proposed probe is tested virtually by solving the forward problem of skin deformation by a finite element (FE) model. The forward problem is a function of the material parameters. In order to solve the inverse problem of determining skin properties from suction data, we use a Bayesian framework. The FE model is an expensive forward function, and is thus substituted with a Gaussian process metamodel to enable the Bayesian inference problem.

The Finite Element Analysis Research on Microneedle Design Strategy and Transdermal Drug Delivery System

Microneedles (MNs) as a novel transdermal drug delivery system have shown great potential for therapeutic and disease diagnosis applications by continually providing minimally invasive, portable, cost-effective, high bioavailability, and easy-to-use tools compared to traditional parenteral administrations. However, microneedle transdermal drug delivery is still in its infancy. Many research studies need further in-depth exploration, such as safety, structural characteristics, and drug loading performance evaluation. Finite element analysis (FEA) uses mathematical approximations to simulate real physical systems (geometry and load conditions). It can simplify complex engineering problems to guide the precise preparation and potential industrialization of microneedles, which has attracted extensive attention. This article introduces FEA research for microneedle transdermal drug delivery systems, focusing on microneedle design strategy, skin mechanics models, skin permeability, and the FEA research on drug delivery by MNs.

Ultrasonic microrheology for ex vivo skin explants monitoring: A proof of concept

As an answer to alternative non-animal testing, biosensors dedicated to the ex vivo skin explants monitoring are a challenge to study physiological-like behavior and optimize new topical products. Because of the skin viscoelastic behavior, mechanical tests are commonly based on macroscopic measurement and give global descriptors of its state. Other techniques, including photoacoustic ones, are more focused on the molecular scale. There is a gap to fill in the mesoscopic range to get information about the microstructure of the skin. This article presents the proof-of-concept of a biosensor coupling a thickness shear-mode transducer with human skin explants kept in life-like state for a week. Thanks to a multifrequency analysis of the transducer impedance, this biosensor is able to monitor the viscoelastic properties of the skin. To extract the complex shear modulus and the microstructural evolutions, a mechanical model based on fractional calculus is used. As a preliminary results, the sensitivity of the sensor to probe the skin viscoelasticity in lifelike state and the impact of its culture medium are presented. A suitable microstructural coefficient is also extracted in order to identify mechanical breaches in the skin barrier after the application of peeling products.

A Review on Damage and Rupture Modelling for Soft Tissues

Computational modelling of damage and rupture of non-connective and connective soft tissues due to pathological and supra-physiological mechanisms is vital in the fundamental understanding of failures. Recent advancements in soft tissue damage models play an essential role in developing artificial tissues, medical devices/implants, and surgical intervention practices. The current article reviews the recently developed damage models and rupture models that considered the microstructure of the tissues. Earlier review works presented damage and rupture separately, wherein this work reviews both damage and rupture in soft tissues. Wherein the present article provides a detailed review of various models on the damage evolution and tear in soft tissues focusing on key conceptual ideas, advantages, limitations, and challenges. Some key challenges of damage and rupture models are outlined in the article, which helps extend the present damage and rupture models to various soft tissues.

A hyperelastic model to capture the mechanical behaviour and histological aspects of the soft tissues

It is well established that the soft connective tissues show a nonlinear elastic response that comes from their microstructural arrangement. Tissues' microstructure alters with various physiological conditions and may affect their mechanical responses. Therefore, the accurate prediction of tissue's mechanical response is crucial for clinical diagnosis and treatments. Thus, a physically motivated and mathematically simplified model is required for the accurate prediction of tissues' mechanical and structural responses. This study explored the 'Exp-Ln' hyperelastic model (Khajehsaeid et al., 2013) to capture soft tissues' mechanical and histological behaviour. In this work, uniaxial tensile test data for the belly and back pig skin were extracted from the experiments performed in our laboratory, whereas uniaxial test data for other soft tissues (human skin, tendon, ligament, and aorta) were extracted from the literature. The 'Exp-Ln; and other hyperelastic models (e.g. Money Rivlin, Ogden, Yeoh, and Gent models) were fitted with these experimental data, and obtained results were compared between the models. These results show that the 'Exp-Ln' model could capture the mechanical behaviour of soft tissues more accurately than other hyperelastic models. This model was also found numerically stable for all modes and ranges of deformation. This study also investigated the link between 'Exp-Ln' material parameters and tissue's histological parameters. The histological parameters such as collagen content, fibre free length, crosslink density, and collagen arrangement were measured using staining and ATR-FTIR techniques. The material parameters were found statistically correlated with the histological parameters. Further, 'Exp-Ln' model was implemented in ABAQUS through the VUMAT subroutine, where the mechanical behaviour of various soft tissues was simulated for different modes of deformation. The finite element analysis results obtained using the 'Exp-Ln' model agreed with the experiments and were more accurate than other hyperelastic models. Overall, these results demonstrate the capability of 'Exp-Ln' model to predict the mechanical and structural responses of the soft tissues.

Changes in Tissue Composition and Load Response After Transtibial Amputation Indicate Biomechanical Adaptation

Despite the potential for biomechanical conditioning with prosthetic use, the soft tissues of residual limbs following lower-limb amputation are vulnerable to damage. Imaging studies revealing morphological changes in these soft tissues have not distinguished between superficial and intramuscular adipose distribution, despite the recognition that intramuscular fat levels indicate reduced tolerance to mechanical loading. Furthermore, it is unclear how these changes may alter tissue tone and stiffness, which are key features in prosthetic socket design. This study was designed to compare the morphology and biomechanical response of limb tissues to mechanical loading in individuals with and without transtibial amputation, using magnetic resonance imaging in combination with tissue structural stiffness. The results revealed higher adipose infiltrating muscle in residual limbs than in intact limbs (residual: median 2.5% (range 0.2-8.9%); contralateral: 1.7% (0.1-5.1%); control: 0.9% (0.4-1.3%)), indicating muscle atrophy and adaptation post-amputation. The intramuscular adipose content correlated negatively with daily socket use, although there was no association with time post-amputation. Residual limbs were significantly stiffer than intact limbs at the patellar tendon site, which plays a key role in load transfer across the limb-prosthesis interface. The tissue changes following amputation have relevance in the clinical understanding of prosthetic socket design variables and soft tissue damage risk in this vulnerable group.

Mechanical modeling and characterization of human skin: A review

This paper reviews the advances made in recent years on modeling approaches and experimental techniques to characterize the mechanical properties of human skin. The skin is the largest organ of the human body that has a complex multi-layered structure with different mechanical behaviors. The mechanical properties of human skin play an important role in distinguishing between healthy and unhealthy skin. Furthermore, knowing these mechanical properties enables computer simulation, skin research, clinical studies, as well as diagnosis and treatment monitoring of skin diseases. This paper reviews the recent efforts on modeling skin using linear, nonlinear, viscoelastic, and anisotropic materials. The work also focuses on aging effects, microstructure analysis, and non-invasive methods for skin testing. A detailed explanation of the skin structure and numerical models, such as finite element models, are discussed in this work. This work also compares different experimental methods that measure the mechanical properties of human skin. The work reviews the experimental results in the literature and shows how the mechanical properties of human skin vary with the skin sites, the layers, and the structure of human skin. The paper also discusses how state-of-the-art technology can advance skin research.

Mechanical properties of whole-body soft human tissues: a review

The mechanical properties of soft tissues play a key role in studying human injuries and their mitigation strategies. While such properties are indispensable for computational modelling of biological systems, they serve as important references in loading and failure experiments, and also for the development of tissue simulants. To date, experimental studies have measured the mechanical properties of peripheral tissues (e.g. skin)in-vivoand limited internal tissuesex-vivoin cadavers (e.g. brain and the heart). The lack of knowledge on a majority of human tissues inhibit their study for applications ranging from surgical planning, ballistic testing, implantable medical device development, and the assessment of traumatic injuries. The purpose of this work is to overcome such challenges through an extensive review of the literature reporting the mechanical properties of whole-body soft tissues from head to toe. Specifically, the available linear mechanical properties of all human tissues were compiled. Non-linear biomechanical models were also introduced, and the soft human tissues characterized using such models were summarized. The literature gaps identified from this work will help future biomechanical studies on soft human tissue characterization and the development of accurate medical models for the study and mitigation of injuries.

Effect of collagen degradation on the mechanical behavior and wrinkling of skin

Chronological skin aging is a complex process that is controlled by numerous intrinsic and extrinsic factors. One major factor is the gradual degradation of the dermal collagen fiber network. As a step toward understanding the mechanistic importance of dermal tissue in the process of aging, this study employs analytical and multiscale computational models to elucidate the effect of collagen fiber bundle disintegration on the mechanical properties and topography of skin. Here, human skin is modeled as a soft composite with an anisotropic dermal layer. The anisotropy of the tissue is governed by collagen fiber bundles with varying densities, average fiber alignments, and normalized alignment distributions. In all finite element models examined, collagen fiber bundle degradation results in progressive decreases in dermal and full-thickness composite stiffness. This reduction is more profound when collagen bundles align with the compression axis. Aged skin models with low collagen fiber bundle densities under compression exhibit notably smaller critical wrinkling strains and larger critical wavelengths than younger skin models, in agreement with in vivo wrinkling behavior with age. The propensity for skin wrinkling can be directly attributable to the degradation of collagen fiber bundles, a relationship that has previously been assumed but unsubstantiated. While linear-elastic analytical models fail to capture the postbuckling behavior in skin, nonlinear finite element models can predict the complex bifurcations of the compressed skin with different densities of collagen bundles.

Changes in the three-dimensional microscale topography of human skin with aging impact its mechanical and tribological behavior

Human skin enables interaction with diverse materials every day and at all times. The ability to grasp objects, feel textures, and perceive the environment depends on the mechanical behavior, complex structure, and microscale topography of human skin. At the same time, abrasive interactions, such as sometimes occur with prostheses or textiles, can damage the skin and impair its function. Previous theoretical and computational efforts have shown that skin's surface topography or microrelief is crucial for its tribological behavior. However, current understanding is limited to adult surface profiles and simplified two-dimensional simulations. Yet, the skin has a rich set of features in three dimensions, and the geometry of skin is known to change with aging. Here we create a numerical model of a dynamic indentation test to elucidate the effect of changes in microscale topography with aging on the skin's response under indentation and sliding contact with a spherical indenter. We create three different microrelief geometries representative of different ages based on experimental reports from the literature. We perform the indentation and sliding steps, and calculate the normal and tangential forces on the indenter as it moves in three distinct directions based on the characteristic skin lines. The model also evaluates the effect of varying the material parameters. Our results show that the microscale topography of the skin in three dimensions, together with the mechanical behavior of the skin layers, lead to distinctive trends on the stress and strain distribution. The major finding is the increasing role of anisotropy which emerges from the geometric changes seen with aging.

A method to analyze the influence of mechanical strain on dermal collagen morphologies

Collagen fibers and their orientation play a major role in the mechanical behavior of soft biological tissue such as skin. Here, we present a proof-of-principle study correlating mechanical properties with collagen fiber network morphologies. A dedicated multiphoton stretching device allows for mechanical deformations in combination with a simultaneous analysis of its collagen fiber network by second harmonic generation imaging (SHG). The recently introduced Fiber Image Network Evaluation (FINE) algorithm is used to obtain detailed information about the morphology with regard to fiber families in collagen network images. To demonstrate the potential of our method, we investigate an isotropic and an anisotropic ex-vivo dorsal pig skin sample under quasi-static cyclic stretching and relaxation sequences. Families of collagen fibers are found to form a partially aligned collagen network under strain. We find that the relative force uptake is accomplished in two steps. Firstly, fibers align within their fiber families and, secondly, fiber families orient in the direction of force. The maximum alignment of the collagen fiber network is found to be determined by the largest strain. Isotropic and anisotropic samples reveal a different micro structural behavior under repeated deformation leading to a similar force uptake after two stretching cycles. Our method correlates mechanical properties with morphologies in collagen fiber networks.

Numerical analysis of the strain distribution in skin domes formed upon the application of hypobaric pressure

BACKGROUND

Suction cups are widely used in applications such as in measurement of mechanical properties of skin in vivo, in drug delivery devices or in acupuncture treatment. Understanding mechanical response of skin under hypobaric pressure is of great importance for users of suction cups. The aim of this work is to predict the hypobaric pressure induced 3D stretching of the skin.

METHODS

Experimental skin tensile tests were carried out for mechanical property characterization. Both linear elasticity and hyperelasticity parameters were determined and implemented in Finite Element modelling. Skin suction tests were performed in both experiments and FEM simulations for model validation. 3D skin stretching is then visualized in detail in FEM simulations.

RESULTS

The simulations showed that the skin was compressed consistently along the thickness direction, leading to reduced thickness. At the center of the dome, the radial and angular strain decreases from the top surface to the bottom surface, although always in tension. Hyperelasticity modelling showed superiority over linear elasticity modelling while predicting the strain distribution because the stretch ratio reaches values exceeding the initial linear elastic stage of the stress-strain curve for skin.

CONCLUSION

Hyperelasticity modelling is an effective approach to predict the 3D strain distribution, which paves a way to accurately design safe commercial products that interface with the skin.

On the gular sac tissue of the brown pelican: Structural characterization and mechanical properties

The brown pelican (Pelecanus occidentalis) wields one of the largest bills of any bird and is distinguished by the deployable throat pouch of extensible tissue used to capture prey. Here we report on mechanical properties and microstructure of the pouch skin. It exhibits significant anisotropy, with the transverse direction having maximum nominal tensile strains of 200% to 300%, triple the value in the longitudinal direction. This is a higher extensibility than most conventional skin and is the result of the requirement of the sac to net fish; it should expand laterally, with controlled longitudinal stretch. Transmission electron microscopy provides microstructural evidence of the directionality of the collagen fibers and reveals the individual collagen fibrils with a bimodal diameter distribution having peaks at 100 and 170 nm. These dimensions are similar to collagen in mammal skin. In the lateral direction, the fibers form a curvy pattern with a radius of approximately 2 µm wherein the fibrils reorient, straighten, slide, and stretch elastically under tensile load. A second mechanism operates in the transverse direction; the membrane forms a corrugated pattern that, upon straightening of collagen fibrils, confers additional extensibility. This elicits the anisotropic response observed in tensile testing. This work focuses on the mechanical characterization based on the effect of relative bird age, sample location on the pouch, and strain rate. Anterior-posterior location and strain rate are not major influencers on exhibited strengths and extensibilities. However, bird age and dorsal-ventral location are found to affect the mechanical response of the pouch significantly. A physically-based constitutive model is developed for the middle layer of the gular sac, based on observations, which predicts maximum stresses, strains, and the shape of the stress-strain curve consistent with the experimental results.

In vivo panoramic human skin shape and deformation measurement using mirror-assisted multi-view digital image correlation

Panoramic shape and deformation measurements of human skin in vivo may provide important information for biomechanical analysis, exercise guidance and medical diagnosis. This work proposes the application of an advanced mirror-assisted multi-view digital image correlation (DIC) method for dynamic measurements of 360-deg shape and deformation of human body parts in vivo. The main advantage of this method consists in its capabilities to perform full-panoramic non-contact measurements with a single pair of synchronized cameras and two planar mirrors thus representing a lean yet effective alternative to conventional multi-camera DIC systems in 'surrounding' configuration. We demonstrate the capabilities of this method by measuring the full-panoramic shape of a plastic human head, the deformation of a woman face and the principal strain distribution over the full-360-deg surface of a forearm during fist clenching. The applications of this method can be the most disparate but, given the possibility to determine the full-field strains and derived information (e.g. skin tension lines), we envisage a great potential for the study of skin biomechanics in vivo.

An open source pipeline for design of experiments for hyperelastic models of the skin with applications to keloids

The aim of this work is to characterize the mechanical parameters governing the in-plane behavior of human skin and, in particular, of a keloid-scar. We consider 2D hyperelastic bi-material model of a keloid and the surrounding healthy skin. The problem of finding the optimal model parameters that minimize the misfit between the model observations and the in vivo experimental measurements is solved using our in-house developed inverse solver that is based on the FEniCS finite element computational platform. The paper focuses on the model parameter sensitivity quantification with respect to the experimental measurements, such as the displacement field and reaction force measurements. The developed tools quantify the significance of different measurements on different model parameters and, in turn, give insight into a given model's ability to capture experimental measurements. Finally, an a priori estimate for the model parameter sensitivity is proposed that is independent of the actual measurements and that is defined in the whole computational domain. This estimate is primarily useful for the design of experiments, specifically, in localizing the optimal displacement field measurement sites for the maximum impact on model parameter inference.

Frequency dependent inelastic response of collagen architecture of pig dermis under cyclic tensile loading: An experimental study

The evaluation of collagen architecture of the dermis in response to mechanical stimulation is important as it affects the macroscopic mechanical properties of the dermis. A detailed understanding of the processes involved in the alteration of the collagen structure is required to correlate the mechanical stimulation with tissue remodeling. This study investigated the effect of cyclic frequencies i.e. low (0.1 Hz), medium (2.0 Hz), and high (5.0 Hz) (physiological range) in the alteration of pig dermis collagen structure and its correlation with the macroscopic mechanical response of the dermis. The assessment of the collagen structure of virgin and mechanical tested specimens at tropocollagen, collagen fibril, and fiber level was performed using Fourier-transform infrared-attenuated total reflection (FTIR-ATR), atomic force microscopy (AFM), and scanning electron microscopy (SEM) respectively. After 10³ cycles, a significantly higher alteration in collagen structure with discrete plastic-type damage was found for low frequency. This frequency dependent alteration of the collagen structure was found in correlation with the dermis macroscopic response. The value of inelastic strain, stress softening, damage parameter (reduction in elastic modulus), and reduction in energy dissipation were observed significantly large for slow frequency. A power-law based empirical relations, as a function of frequency and number of cycles, were proposed to predict the value of inelastic strain and damage parameter. This study also suggests that hierarchical structural response against the mechanical stimulation is time-dependent rather than cycle-dependent, may affect the tissue remodeling.

Is there any objective and independent characterization and modeling of soft biological tissues?

The characterization of soft tissue raises several difficulties. Indeed, soft biological tissues usually shrink when dissected from their in vivo location. This shrinkage is characteristic of the release of residual stresses, since soft tissues are indeed often pre-stressed in their physiological configuration. During experimental loading, large extension at very low level of force are expected and assumed to be related to the progressive recruitment and stretching of fibers. However, the first phase of the mechanical test is also aiming at recovering the pre-stressed in vivo behavior. As a consequence, the initial phase, corresponding to the recovering of prestress and/or recruitment of fiberes, is questionable and frequently removed. One of the preferred methods to erase it consists in applying a preforce or prestress to the sample: this allows to easily get rid of the sample retensioning range. However this operation can impact the interpretation of the identified mechanical parameters. This study presents an evaluation of the impact of the data processing on the mechanical properties of a numerically defined material. For this purpose, a finite element simulation was performed to replicate a uniaxial tensile test on a biological soft tissue sample. The influence of different pre-stretches on the mechanical parameters of a second order Yeoh model was investigated. The Yeoh mechanical parameters, or any other strain energy density, depend strongly on any pre- and post-processing choices: they adapt to compensate the error made when choosing an arbitrary level of prestretch or prestress. This observation spreads to any modeling approach used in soft tissues. Mechanical parameters are indeed naturally bound to the choice of the pre-stretch (or pre-stress) through the elongation and the constitutive law. Regardless of the model, it would therefore be pointless to compare mechanical parameters if the conditions for the processing of experimental raw data are not fully documented.

Hyperelastic Ex Vivo Cervical Tissue Mechanical Characterization

This paper presents the results of the comparison between a proposed Fourth Order Elastic Constants (FOECs) nonlinear model defined in the sense of Landau's theory, and the two most contrasted hyperelastic models in the literature, Mooney-Rivlin, and Ogden models. A mechanical testing protocol is developed to investigate the large-strain response of ex vivo cervical tissue samples in uniaxial tension in its two principal anatomical locations, the epithelial and connective layers. The final aim of this work is to compare the reconstructed shear modulus of the epithelial and connective layers of cervical tissue. According to the obtained results, the nonlinear parameter A from the proposed FOEC model could be an important biomarker in cervical tissue diagnosis. In addition, the calculated shear modulus depended on the anatomical location of the cervical tissue (μepithelial = 1.29 ± 0.15 MPa, and μconnective = 3.60 ± 0.63 MPa).

Numerical Investigation of Auxetic Textured Soft Strain Gauge for Monitoring Animal Skin

Recent advances in hyperelastic materials and self-sensing sensor designs have enabled the creation of dense compliant sensor networks for the cost-effective monitoring of structures. The authors have proposed a sensing skin based on soft polymer composites by developing soft elastomeric capacitor (SEC) technology that transduces geometric variations into a measurable change in capacitance. A limitation of the technology is in its low gauge factor and lack of sensing directionality. In this paper, we propose a corrugated SEC through surface texture, which provides improvements in its performance by significantly decreasing its transverse Poisson's ratio, and thus improving its sensing directionality and gauge factor. We investigate patterns inspired by auxetic structures for enhanced unidirectional strain monitoring. Numerical models are constructed and validated to evaluate the performance of textured SECs, and to study their performance at monitoring strain on animal skin. Results show that the auxetic patterns can yield a significant increase in the overall gauge factor and decrease the stress experienced by the animal skin, with the re-entrant hexagonal honeycomb pattern outperforming all of the other patterns.

Constitutive description of skin dermis: Through analytical continuum and coarse-grained approaches for multi-scale understanding

Although there are many successful descriptions of the mechanical response of dermis at different levels of complexity and incorporating varying degrees of the physical phenomena involved in deformation, observations indicate that the unraveling of fibers involves a complex three-dimensional process in which they interact in ways that resemble a braided pattern. Here we develop two complementary treatments to gain a better understanding of the mechanical response of dermis: a) an analytical treatment incorporating fibril stiffness, interfibrillar frictional sliding, and the effect of lateral fibers on the extension of a primary fiber; b) a coarse-grained molecular dynamics model comprised of an array of parallel curved fibrils simulating a fiber. Interfibrillar frictional sliding and stiffness are also captured. Both analytical and molecular dynamics models operate at a scale compatible with the wavelength of collagen fibers (~10 µm). The constitutive description presented here incorporates important physical processes taking place during deformation of dermis and thus represents an advance in our understanding of these phenomena. STATEMENT OF SIGNIFICANCE: Microstructural observations of the dermis of skin during tensile deformation indicate that the unraveling of fibers involves a complex three-dimensional process which replicates the effects of braiding. Two complementary constitutive modeling treatments were developed to gain a better understanding of the mechanical response of dermis: an analytical treatment incorporating fibril stiffness, interfibrillar sliding, and the effect of transverse fibers; and a coarse-grained molecular dynamics model describing the fibril bundling effect. An important novel aspect of the current contribution is the recognition that tridimensional collagen fiber arrangements play an important role in the mechanical response. The constitutive description presented here incorporates physical processes taking place during deformation of the dermis and thus represents an advance in our understanding of these phenomena.

Characterization of human female breast and abdominal skin elasticity using a bulge test

Characterization of material properties of human skin is required to develop a physics-based biomechanical model that can predict deformation of female breast after cosmetic and reconstructive surgery. In this paper, we have adopted an experimental approach to characterize the biaxial response of human skin using bulge tests. Skin specimens were harvested from breast and abdominal skin of female subjects who underwent mastectomy and/or reconstruction at The University of Texas MD Anderson Cancer Center and who provided informed consent. The specimens were tested within 2 h of harvest, and after freezing for different time periods but not exceeding 6 months. Our experimental results show that storage in a freezer at -20 °C for up to about 40 days does not lead to changes in the mechanical response of the skin beyond statistical variation. Moreover, displacement at the apex of the bulged specimen versus applied pressure varies significantly between different specimens from the same subject and from different subjects. The bulge test results were used in an inverse optimization procedure in order to calibrate two different constitutive material models - the angular integration model proposed by Lanir (1983) and the generalized structure tensor formulation of Gasser et al. (2006). The material parameters were estimated through a cost function that penalized deviations of the displacement and principal curvatures at the apex. Generally, acceptable fits were obtained with both models, although the angular integration model was able to fit the curvatures slightly better than the Gasser et al. model. The range of the model parameters has been extracted for use in physics-based biomechanical models of the breast.

Propagation of uncertainty in the mechanical and biological response of growing tissues using multi-fidelity Gaussian process regression

A key feature of living tissues is their capacity to remodel and grow in response to environmental cues. Within continuum mechanics, this process can be captured with the multiplicative split of the deformation gradient into growth and elastic contributions. The mechanical and biological response during tissue adaptation is characterized by inherent variability. Accounting for this uncertainty is critical to better understand tissue mechanobiology, and, moreover, it is of practical importance if we aim to develop predictive models for clinical use. However, the current gold standard in computational models of growth and remodeling remains the use of deterministic finite element (FE) simulations. Here we focus on tissue expansion, a popular technique in which skin is stretched by a balloon-like device inducing its growth. We construct FE models of tissue expansion with various levels of detail, and show that a sufficiently broad set of FE simulations from these models can be used to train an accurate and efficient multi-fidelity Gaussian process (GP) surrogate. The approach is not limited to simulation data, rather, it can fuse different kinds of data, including from experiments. The main appeal of the framework relies on the common experience that highly detailed models (or experiments) are more accurate but also more costly, while simpler models (or experiments) can be easily evaluated but are bound to have some error. In these situations, doing uncertainty analysis tasks with the high fidelity models alone is not feasible and, conversely, relying solely on low fidelity approximations is also undesirable. We show that a multi-fidelity GP outperforms the high fidelity GP and low fidelity GP when tested against the most detailed FE model. In turn, having trained the multi-fidelity GP model, we showcase the propagation of uncertainty from the mechanical and biological response parameters to the spatio-temporal growth outcomes. We expect that the methods and applications in this paper will enable future research in parameter calibration under uncertainty and uncertainty propagation in real clinical scenarios involving tissue growth and remodeling.

The regional-dependent biaxial behavior of young and aged mouse skin: A detailed histomechanical characterization, residual strain analysis, and constitutive model

Skin fulfills several vital functions, many of which are dependent on its mechanical properties. Therefore, as mice have become an invaluable model for skin research, determining murine skin's mechanical properties is important. Specifically, skin's mechanical properties are important for functional tests as well as for prognostic and diagnostic purposes. Additionally, computational simulations of skin behavior are becoming commonplace, rendering accurate models of murine skin's constitutive behavior necessary. To date, our knowledge of mouse skin mechanics shows significant gaps. For example, there are no comprehensive reports correlating skin's mechanical properties with region, age, and direction. Moreover, mouse skin's residual strain behavior has not been reported on. In our current work, we set out to fill these gaps. Based on histology, 2-photon microscopy, and planar biaxial testing, while accurately tracking various reference configurations, we report on differences in gross structure, microstructural organization, and constitutive response of skin, and cast those properties into a versatile Fung-type hyperelastic constitutive law for three reference configurations. Our data is the most comprehensive report contrasting the mechanical properties of young (12 weeks) and aged (52 weeks) mouse skin and will, thus, be valuable to basic science as control data, and provide accurate constitutive laws for mouse skin modeling. STATEMENT OF SIGNIFICANCE: Our findings are significant as they fill several gaps in our understanding of mouse skin mechanics. This is particularly important as mouse skin is becoming a frequent and critical model of human skin for cosmetic and medical science. Specifically, we quantified how mechanical properties of mice skin vary with age, with location, and with direction. Additionally, we cast our findings into constitutive models that can be used by others for predictive computer simulations of skin behavior.

Structural and Biomechanical Characterization of a Scaffold-Free Skin Equivalent Model via Biophysical Methods

AIMS

Among in vitro skin models, the scaffold-free skin equivalent (SFSE), without exogenous material, is interesting for pharmacotoxicological studies. Our aim was to adapt in vivo biophysical methods to study the structure, thickness, and extracellular matrix of our in vitro model without any chemical fixation needed as for histology.

METHODS

We evaluated 3 batches of SFSE and characterized them by histology, transmission electron microscopy (TEM), and immunofluorescence. In parallel, we investigated 3 biophysical methods classically used for in vivo evaluation, optical coherence tomography (OCT), and laser scanning microscopy (LSM) imaging devices as well as the cutometer suction to study the biomechanical properties.

RESULTS

OCT allowed the evaluation of SFSE total thickness and its different compartments. LSM has a greater resolution enabling an evaluation at the cell scale and the orientation of collagen fibers. The viscoelasticity measurement by cutometry was possible on our thin skin model and might be linked with mature collagen bundles visible in TEM and LSM and with elastic fibers seen in immunofluorescence.

CONCLUSION

Our data demonstrated the simplicity and sensitivity of these different in vivo biophysical devices on our thin skin model. These noninvasive tools allow to study the morphology and the biomechanics of in vitro models.

Maximum tensile stress and strain of skin of the domestic pig-differences concerning pigs from organic and non-organic farming

The purpose of this work has been to determine differences in biomechanical properties of porcine skin from organic and non-organic farming as porcine skin is widely used as a model for human skin. A test apparatus was used, using gravity to stretch and finally tear a dumbbell-shaped specimen of prepared abdominal skin with a testing surface area of 25 × 4 mm. A total of 32 specimens were taken from seven individual pigs, three from organic and four from non-organic farming, in different orientations with respect to the Langer's lines. The tests were performed at a dynamic speed of around 1.66 m/s (corresponding to a nominal strain rate of 67 s^-1). Engineering strain at rupture was higher in pig skin from non-organic farming with values up to 321% as opposed to 90% in organic pig skin. The maximum tensile stress found in non-organic pig skin was lower than in pig skin from organic farming with maximum values of 34 MPa as opposed to 58 MPa. The reason for the difference in biomechanical properties is unclear; the effect of sunlight is discussed as well as other factors like age and exercise. It seems that the biomechanical properties of porcine skin from organic farming are more similar to those of human skin.

Linking microvascular collapse to tissue hypoxia in a multiscale model of pressure ulcer initiation

Pressure ulcers are devastating injuries that disproportionately affect the older adult population. The initiating factor of pressure ulcers is local ischemia, or lack of perfusion at the microvascular level, following tissue compression against bony prominences. In turn, lack of blood flow leads to a drop in oxygen concentration, i.e, hypoxia, that ultimately leads to cell death, tissue necrosis, and disruption of tissue continuity. Despite our qualitative understanding of the initiating mechanisms of pressure ulcers, we are lacking quantitative knowledge of the relationship between applied pressure, skin mechanical properties as well as structure, and tissue hypoxia. This gap in our understanding is, at least in part, due to the limitations of current imaging technologies that cannot simultaneously image the microvascular architecture, while quantifying tissue deformation. We overcome this limitation in our work by combining realistic microvascular geometries with appropriate mechanical constitutive models into a microscale finite element model of the skin. By solving boundary value problems on a representative volume element via the finite element method, we can predict blood volume fractions in response to physiological skin loading conditions (i.e., shear and compression). We then use blood volume fraction as a homogenized variable to couple tissue-level skin mechanics to an oxygen diffusion model. With our model, we find that moderate levels of pressure applied to the outer skin surface lead to oxygen concentration contours indicative of tissue hypoxia. For instance, we show that applying a pressure of 60 kPa at the skin surface leads to a decrease in oxygen partial pressure from a physiological value of 65 mmHg to a hypoxic level of 31 mmHg. Additionally, we explore the sensitivity of local oxygen concentration to skin thickness and tissue stiffness, two age-related skin parameters. We find that, for a given pressure, oxygen concentration decreases with decreasing skin thickness and skin stiffness. Future work will include rigorous calibration and validation of this model, which may render our work an important tool toward developing better prevention and treatment tools for pressure ulcers specifically targeted toward the older adult patient population.

Non-invasive in vivo quantification of human skin tension lines

Human skin is a composite tissue that exhibits anisotropic mechanical properties. This anisotropy arises primarily from the alignment of collagen and elastin fibers in the dermis, which causes the skin to exhibit greater tension in one direction, making it appear stiffer. A diverse number of skin tension guidelines have been developed to assist surgeons in making incisions that produce the least conspicuous scars. However, skin anisotropy is believed to vary from subject to subject, and no single guideline is universally recognized as the best to implement for surgical applications. To date, no system exists that can rapidly and non-invasively measure lines of skin tension in vivo. In this article, we evaluate the ability of a new aspiration system to measure the anisotropy of human skin. The device painlessly applies a radial stress of 17 kPa to a region of skin, and captures radially asymmetric skin deformations via a dermal camera. These deformations are used to quantify orientations of strain extrema and the direction of greatest skin stiffness. The ratio of these asymmetric strains varies between 1 and -0.75. A simple 2D transverse isotropic model captures this behavior for multiple anatomical sites. Clinical trials reveal that skin tension line orientations are comparable with existing skin tension maps and generally agree across subjects, however orientations statistically differ between individuals. As such, existing guidelines appear to provide only approximate estimates of skin tension orientation. STATEMENT OF SIGNIFICANCE: Skin tension lines (STL) in human skin arise primarily from collagen fiber alignment in the dermis. These lines are used by surgeons to guide incisions that produce the least conspicuous scars. While numerous anatomical STL maps exist, no single guideline is universally recognized as the most reliable. Moreover, manual methods of quantifying STL are imprecise. For the first time, we have developed a device capable of rapidly and non-invasively measuring STL orientations in vivo, using a single test. Our results are used to establish a simple constitutive model of mechanical skin anisotropy. Clinical trials further reveal STL orientations are comparable with existing maps, but statistically differ between individuals. Existing guidelines therefore appear to provide only approximate estimates of STL orientation.

Tensile behavior and structural characterization of pig dermis

Skin, the outermost layer of the body, fulfills a broad range of functions, protecting internal organs from damage and infection, while regulating the body's temperature and water content via the exchange of heat and fluids. It must be able to withstand and recover from extensive deformation and damage that can occur during growth, movement, and potential injuries. A detailed investigation of the evolution of the collagen architecture of the dermis as a function of deformation is conducted, which reveals new aspects that help us to understand the mechanical response of skin. Juvenile pig is used as a model material because of its similarity to human skin. The dermis is found to have a tridimensional woven structure of collagen fibers, which evolves with deformation. After failure, we observe that the fibers have straightened and aligned in the direction of tension. The effects of strain-rate change, cyclic loading, stress relaxation, and orientation are quantitatively established. Digital image correlation techniques are implemented to quantify skin's anisotropy; measurements of the Poisson ratio are reported. This is coupled with transmission electron microscopy which enables obtaining quantitative strain parameters evaluated through the orientation and curvature of the collagen fibers and their changes, for the first time in all three dimensions of the tissue. A model experiment using braided human hair in tension exhibits a similar J-curve response to skin, and we propose that this fiber configuration is at least partially responsible for the monotonic increase of the tangent modulus of skin with strain. The obtained results are intended to serve as a basis for structurally-based models of skin. STATEMENT OF SIGNIFICANCE: Our study reveals a new aspect of the dermis: it is comprised of a tridimensional woven structure of collagen fibers, which evolves with deformation. This is enabled by primarily two techniques, transmission electron microscopy on three perpendicular planes and confocal images with second harmonic generation fluorescence of collagen, captured at different intervals of depth. After failure, the fibers have straightened and aligned in the direction of tension. Digital image correlation techniques are implemented to quantify skin's anisotropy; measurements of the Poisson ratio are reported. A model experiment using braided human hair in tension exhibits a similar J-curve response to skin, and we propose that this fiber configuration is at least partially responsible for the monotonic increase of the tangent modulus of skin with strain.

On the compressibility and poroelasticity of human and murine skin

A total of 37 human and 33 murine skin samples were subjected to uniaxial monotonic, cyclic, and relaxation experiments. Detailed analysis of the three-dimensional kinematic response showed that skin volume is significantly reduced as a consequence of a tensile elongation. This behavior is most pronounced in monotonic but persists in cyclic tests. The dehydration associated with volume loss depends on the osmolarity of the environment, so that tension relaxation changes as a consequence of modifying the ionic strength of the environmental bath. Similar to ex vivo observations, complementary in vivo stretching experiments on human volar forearms showed strong in-plane lateral contraction. A biphasic homogenized model is proposed which allows representing all relevant features of the observed mechanical response.

Biomechanical Modeling of Wounded Skin

Skin injury is the most common type of injury, which manifests itself in the form of wounds and cuts. A minor wound repairs itself within a short span of time. However, deep wounds require adequate care and sometime clinical interventions such as surgical suturing for their timely closure and healing. In literature, mechanical properties of skin and other tissues are well known. However, the anisotropic behavior of wounded skin has not been studied yet, specifically with respect to localized overstraining and possibilities of rupture. In the current work, the biomechanics of common skin wound geometries were studied with a biofidelic skin phantom, using uniaxial mechanical testing and Digital Image Correlation (DIC). Global and local mechanical properties were investigated, and possibilities of rupture due to localized overstraining were studied across different wound geometries and locations. Based on the experiments, a finite element (FE) model was developed for a common elliptical skin wound geometry. The fidelity of this FE model was evaluated with simulation of uniaxial tension tests. The induced strain distributions and stress-stretch responses of the FE model correlated very well with the experiments (R2 > 0.95). This model would be useful for prediction of the mechanical response of common wound geometries, especially with respect to their chances of rupture due to localized overstraining. This knowledge would be indispensable for pre-surgical planning, and also in robotic surgeries, for selection of appropriate wound closure techniques, which do not overstrain the skin tissue or initiate tearing.

Modelling skin wound healing angiogenesis: A review

The occurrence of wounds is a main health concern in Western society due to their high frequency and treatment cost. During wound healing, the formation of a functional blood vessel network through angiogenesis is an essential process. Angiogenesis allows the reestablishment of the normal blood flow, the sufficient exchange of oxygen and nutrients and the removal of metabolic waste, necessary for cell proliferation and viability. Mathematical and computational models provide new tools to improve the healing process. In fact, over the last thirty years, in silico models have been continuously formulated to describe the effect of several biological and mechanical factors in angiogenesis during wound healing. Additionally, with different levels of complexity, these models allow coupling the human skin structure, to distinct cell types and growth factors, to study extracellular matrix composition and to understand its deformation. This paper discusses how in silico models, which are more economical and less time-consuming comparatively to laboratory methodologies, can help test new strategies to promote/optimize angiogenesis. The continuum, cell-based and hybrid mathematical models of wound healing angiogenesis are reviewed in the present paper, in order to identify possible improvements. Accordingly, the development of higher dimension models incorporating multiscale analysis at molecular, cellular and tissue level remains a challenge that future models should consider.

Biomechanical Modeling of Human Skin Tissue Surrogates

Surrogates, which precisely simulate nonlinear mechanical properties of the human skin at different body sites, would be indispensable for biomechanical testing applications, such as estimating the accurate load response of skin implants and prosthetics to study the biomechanics of static and dynamic loading conditions on the skin, dermatological and sports injuries, and estimating the dynamic load response of lethal and nonlethal ballistics. To date, human skin surrogates have been developed mainly with materials, such as gelatin and polydimethylsiloxane (PDMS), based on assumption of simplified mechanical properties, such as an average elastic modulus (estimated through indentation tests), and Poisson’s ratio. In addition, pigskin and cowhides, which have widely varying mechanical properties, have been used to simulate human skin. In the current work, a novel elastomer-based material system is developed, which precisely mimics the nonlinear stress–stretch behavior, elastic modulus at high and low strains, and fracture strengths of the natural human skin at different body sites. The manufacturing and fabrication process of these skin surrogates are discussed, and mechanical testing results are presented.

A model of human skin under large amplitude oscillatory shear

Skin mechanics is of importance in various fields of research when accurate predictions of the mechanical response of skin is essential. This study aims to develop a new constitutive model for human skin that is capable of describing the heterogeneous, nonlinear viscoelastic mechanical response of human skin under shear deformation. This complex mechanical response was determined by performing large amplitude oscillatory shear (LAOS) experiments on ex vivo human skin samples. It was combined with digital image correlation (DIC) on the cross-sectional area to assess heterogeneity. The skin is modeled as a one-dimensional layered structure, with every sublayer behaving as a nonlinear viscoelastic material. Heterogeneity is implemented by varying the stiffness with skin depth. Using an iterative parameter estimation method all model parameters were optimized simultaneously. The model accurately captures strain stiffening, shear thinning, softening effect and nonlinear viscous dissipation, as experimentally observed in the mechanical response to LAOS. The heterogeneous properties described by the model were in good agreement with the experimental DIC results. The presented mathematical description forms the basis for a future constitutive model definition that, by implementation in a finite element method, has the capability of describing the full 3D mechanical behavior of human skin.

Ultrasound elastography reliably identifies altered mechanical properties of burned soft tissues

Although burn injury to the skin and subcutaneous tissues is common in both civilian and military scenarios, a significant knowledge gap exists in quantifying changes in tissue properties as a result of burns. In this study, we present a noninvasive technique based on ultrasound elastography which can reliably assess altered nonlinear mechanical properties of a burned tissue. In particular, ex vivo porcine skin tissues have been exposed to four different burn conditions: (i) 200°F for 10s, (ii) 200°F for 30s, (iii) 450°F for 10s, and (iv) 450°F for 30s. A custom-developed instrument including a robotically controlled ultrasound probe and force sensors has been used to compress the tissue samples to compute two parameters (C₁₀ and C₂₀) of a reduced second-order polynomial hyperelastic material model. The results indicate that while the linear model parameter (C₁₀) does not show a statistically significant difference between the test conditions, the nonlinear model parameter (C₂₀) reliably identifies three (ii-iv) of the four cases (p<0.05) when comparing burned with unburned tissues with a classification accuracy of 60-87%. Additionally, softening of the tissue is observed because of the change in structure of the collagen fibers. The ultrasound elastography-based technique has potential for application under in vivo conditions, which is left for future work.

Elastic compliance of fibrillar assemblies in type I collagen

Fibrillary assemblies of Type I collagen find important applications in tissue engineering and as matrices for biophysical studies. The mechanical and structural properties of these structures are governed by factors such as protein concentration, temperature, pH and ionic strength. This study reports on an impedance based analysis of the elastic compliance of fibrillary assemblies of Type I collagen using quartz crystal microbalance with dissipation (QCM-D) at a fundamental frequency of 5 MHz and overtones (n = 3,5,7,9,11). Here, In situ partial fibrillation of the adsorbing collagen followed by its fibrillary assemblies on hydrophilic gold coated quartz surface have been crosslinked using Gallic acid (GA), Chromium (III) gallate (Cr-GA), Catechin (Cat), Tetrakis(hydroxymethyl)phosphonium sulfate (THPS) and Oxazolidine (Ox). This approach allows direct comparison of how viscoelastic properties track the structural evolution of the fiber and network length scales. The collagen crosslinking shows significant positive impact on the protein's mechanical behaviour and on the type of crosslinking agents used. The elastic modulus increases as collagen <GA < THPS < Cr-GA < Cat < Ox. Atomic force microscopic studies on the adsorbed collagen after cross linking confirmed the presence of fibrous assemblies. The results indicate stabilization and reinforcement through strong physical entanglement between the molecules of collagen as well as chemical interaction between collagen matrix and fibrils during cross linking. The elastic compliance evaluated from ΔDissipation/Δfreq. from QCM-D showed that cross linking with GA, Cr-GA and Ox resulted in flexible fibrillary network while agents like THPS and Cat showed elastic moduli similar to that of pure collagen. Results suggest that optimal collagen-crosslinking agent ratio and degree of crosslinking of collagen can help tailor the mechanical properties for specific applications in design of bio-materials of these composites.

Implementation and validation of constitutive relations for human dermis mechanical response

Finite element models in conjunction with adequate constitutive relations are pivotal in several physiological and medical applications related to both native and engineered tissues, allowing to predict the tissue response under various loading states. In order to get reliable results, however, the validation of the constitutive models is crucial. Therefore, the main purpose of this work is to provide an experimental-computational approach to the biomechanical investigation of soft tissues such as the dermis. This is accomplished by implementing and validating three widely adopted hyperelastic constitutive models (the Ogden, the Holzapfel, and the Gasser-Ogden-Holzapfel laws) supposed to be adequate to reproduce human reticular dermis mechanical behavior. Biaxial experimental data have represented the basis for the determination of the respective material parameters identified thanks to the definition of a cost function accounting for the discrepancy between experimental and predicted data. Afterwards, the experimental tests have been reproduced through finite element simulations. Hence, the constitutive laws have been validated comparing experimental and numerical outcomes in terms of displacements of four reference points and stress-strain relations. Hence, an experimental-numerical framework is proposed for the investigation of collagenous tissues, which could become more accurate with larger and independent experimental datasets. Graphical abstract ᅟ.

Technologies to monitor the health of loaded skin tissues

There are many situations where the skin and underlying soft tissues are compromised by mechanical loading in the form or pressure, or pressure in combination with shear. If sustained, this can lead to damage in the tissues particularly adjacent to bony prominences, resulting in chronic wounds. An array of bioengineering technologies have been adopted to assess the integrity of loaded soft tissues. This paper aims to review these approaches for the quantification, simulation and early detection of mechanically-induced skin damage. The review considers different measurements at the interface between the skin and support surface/medical device, involving pressure, shear, friction and the local microclimate. The potential of the techniques to monitor the physiological response of the skin to these external stimuli including biophysical measurement devices and sampling of biofluids are critically analysed. In addition, it includes an analysis of medical imaging technologies and computational modelling to provide a means by which tissue deformation can be quantified and thresholds for tissue damage defined. Bioengineering measurement and imaging technologies have provided an insight into the temporal status of loaded skin. Despite the advances in technology, to date, the translation to clinical tools which are robust and cost effective has been limited. There is a need to adapt existing technologies and simulation platforms to enable patients, carers and clinicians to employ appropriate intervention strategies to minimise soft tissue damage.