Suction based mechanical characterization of superficial facial soft tissues

On Mechanical Behavior and Characterization of Soft Tissues

The growth and advancements done in solid mechanics and metallurgy have come up with various characterization techniques that help in prediction of elastic properties of different types of materials-isotropic, anisotropic, transverse isotropic, etc. Soft tissues which refer to fibrous tissues, fat, blood vessels, muscles and other tissues that support the body were found to have some control over its mechanical properties. This mechanical behavior of soft tissues has recently shifted the attention of many researchers to develop methods to characterize and describe the mechanical response of soft tissues. The paper discusses the biomechanical nature of soft tissues and the work done to characterize their elastic properties. The paper gives a review of the behavior and characteristics of soft tissues extracted from various experimental tests employed in their characterization. Soft tissues exhibit complex behavior and various complexities are involved in their experimental testing due to their small size and fragile nature. The paper focuses on the conventionally used tensile and compression tests and the difficulties encountered in soft tissue characterization. It also describes the utility of ultrasound technique which is a non-destructive method to characterize soft tissues. Tensile and compression test used to characterize materials are destructive in nature. Ultrasound technique can provide a better way to characterize material in a non-destructive manner.

Second harmonic generation microscopy, biaxial mechanical tests and fiber dispersion models in human skin biomechanics

Advanced numerical simulations of the mechanical behavior of human skin require thorough calibration of the material's constitutive models based on experimental ex vivo mechanical tests along with images of tissue microstructure for a variety of biomedical applications. In this work, a total of 14 human healthy skin samples and 4 additional scarred skin samples were experimentally analyzed to gain deep insights into the biomechanics of human skin. In particular, second harmonic generation (SHG) microscopy was used to extract detailed images of the distribution of collagen fibers, which were subsequently processed using a three-dimensional Fourier transform-based method recently proposed by the authors to quantify the distribution of fiber orientations. Mechanical tests under both biaxial and uniaxial loading were performed to calibrate the relevant mechanical parameters of two widely used constitutive models of soft fiber-reinforced biological tissues that account for non-symmetrical fiber dispersion. The calibration of the models allowed us to identify correlations between the mechanical parameters of the constitutive models considered. STATEMENT OF SIGNIFICANCE: Constitutive models for soft collagenous tissues can accurately reproduce the complex nonlinear and anisotropic mechanical behavior of skin. However, a comprehensive analysis of both microstructural and mechanical parameters is still missing for human skin. In this study, these parameters are determined by combining biaxial mechanical tests and SHG stacks of collagen fibers on ex vivo healthy human skin samples. The constitutive parameters are provided for two widely used hyperelastic models and enable accurate characterization of skin mechanical behavior for advanced numerical simulations.

Skin viscoelasticity effects on the periodic mechanical stimuli propagation between skin layers

Our daily lives are constantly surrounded by dynamic stimuli, and our skin is deformed in a time-dependent manner. Although skin plays an important role in transmitting stimuli received at the surface to mechanoreceptors, few studies have investigated how differences in skin viscoelasticity affect the mechanical stimuli propagation in the skin. Therefore, using a finite element model, we evaluated the effects and trends of changes in the stiffness and viscoelasticity of the skin on the propagation of mechanical quantities between skin layers where mechanoreceptors are present when subjected to periodic stimuli. First, we constructed a new, sophisticated mathematical model of skin viscoelasticity based on the history-dependent deformation behavior of human skin obtained experimentally. We were able to construct a skin model that thoroughly reproduced the actual human skin deformation behavior at oscillations as fast as 10 Hz by setting viscoelastic parameters with a short time constant (0.001-0.006 s). Then, we calculated how skin material parameters affect the propagation of the mechanical quantities in the skin during the history-dependent skin deformation response to periodic stimuli. The finite element analysis showed that not only stiffness but also viscoelasticity markedly affected the mechanical stimuli propagation in the skin, and the effect differed depending on the layer. In particular, greater immediate responsiveness of the dermis contributed to greater propagation of the mechanical stimulus. Our results indicate that more attention needs to be given to the differences in the time-dependent intradermal mechanical stimuli propagation caused by individual's skin viscoelasticity.

Stratum corneum compliance enhances tactile sensitivity through increasing skin deformation: A study protocol for a randomized controlled trial

BACKGROUND

Tactile sensation plays a crucial role in object manipulation, communication, and even emotional well-being. It has been reported that the deformability of skin (also described as skin compliance) that shows a large mechanical response to stimuli is associated with high tactile sensitivity. However, although the compliance of the stratum corneum, the outermost layer of skin, can change daily due to skin care and environmental factors, few studies have quantified the effect of the stratum corneum on tactile sensation.

AIMS

We investigated the changes in tactile sensitivity resulting from skin hydration and identified corresponding alterations in the compliance of the stratum corneum.

METHODS

A randomized controlled trial was conducted. Participants were randomly assigned to an intervention group (n = 20) that had a moisturizing cream applied to their cheeks or a control group (n = 19) that had Milli-Q water applied to their cheeks. Tactile discrimination performance was assessed using psychophysical techniques before and after application. The water content, mechanical response characteristics, and penetration of PEG/PPG-17/4 dimethyl ether from the cream in the stratum corneum were evaluated to identify hydration effects. Skin deformations occurring during tactile sensation were measured concurrently using a suction device employed for tactile stimulation.

RESULTS

Tactile sensitivity was increased in participants who had cream applied to the skin surface, while no significant change was observed in participants who received Milli-Q water. The improved discrimination of tactile stimulus intensity was directly related to the magnitude of skin displacement. The higher water content of the stratum corneum due to cream application decreased the dynamic modulus of elasticity of the stratum corneum and increased the skin's extensibility in response to tactile stimuli.

CONCLUSIONS

Hydrating the stratum corneum significantly enhances tactile sensitivity and is accompanied by an increase in skin extensibility, a factor in tactile intensity perception. The compliance of the thin stratum corneum layer plays a crucial role in tactile experiences that involve skin stretching.

Characterization of the skin-to-bone mechanical interaction on porcine scalp: A combined experimental and computational approach

Fasciae are soft tissues permitting a large but finite sliding between organs, but also between skin and its underlying elements. The contribution of fasciae has been seldomly reported in the literature, and is usually neglected or overly simplified within simulations. In the present contribution, we propose to use peeling tests in order to quantify the skin-to-bone interaction associated with a simple computational approach based on a geometrical modeling of the skin-to-bone interface. To this aim, a new experimental set up combined with a computational model to characterize the skin-to-bone interaction were proposed. The current work is devoted to the porcine scalp complex since it constitutes a common mechanical surrogate for the human scalp complex. The ad hoc computational approach and peeling set up were firstly evaluated on a validation material, before being used to characterize the skin-to-bone interaction within 6 porcine specimens harvested from the scalp. Our experimental setup allowed to measure the peeling response of porcine scalp, showing a three-regimes response including a plateau force. The computational approach satisfyingly reproduced the peeling response based uniquely on experimental-based parameters and on a discrete modeling of skin-to-bone interface. The presented methodology is a first attempt to propose a computationally efficient geometrically based model able to take into account the skin-to-bone interaction up to failure and corroborated by experimental data, and may be largely extended to the modeling of soft interactions between biological human tissues in the future.

A comparison of commercially available synthetic skin substitutes for surgical simulation training

Objective

Simulation training provides an important opportunity to accelerate surgical skills acquisition whilst safeguarding patients. This study compares the suitability of different synthetic skin substitutes for use in surgical simulation training.

Design

Data was collected for eight commercially available synthetic skin substitutes and included cost, delivery time, subjective assessment of fidelity by surgeons and trainees, and objective comparison with the biomechanics of human skin was made through cutometry and durometry measurements. Cutometry and durometry data was collected from three healthy adults from the forearm, forehead and back, with measurements being repeated in triplicate. Subjective assessment of skin pad quality was collected using an 8-criteria questionnaire, graded using a 5-point Likert scale for fidelity to normal skin.

Results

The questionnaire assessment was completed by 30 trainees and practitioners. Overall, felt pads received the poorest outcomes in all criteria; cutometry and durometry results demonstrate poor similarity to skin, and felt received the lowest scores in the questionnaire, although the cheapest. Foam dressings were similar in both cutometric and durometric properties to skin of the face, back and arm. Clinical outcomes of foam dressings were similar to the most expensive commercial skin pad.

Conclusions

Bilaminar foam-based dressings provide a low cost, high fidelity non-biological simulation of skin for surgical training, which is non-inferior to more expensive specifically designed products. Many products designed to act as skin substitutes for surgical simulation fail to adequately replicate the anatomical and mechanical properties of skin.

Identifiability of soft tissue constitutive parameters from in-vivo macro-indentation

Reliable identification of soft tissue material parameters is frequently required in a variety of applications, particularly for biomechanical simulations using finite element analysis (FEA). However, determining representative constitutive laws and material parameters is challenging and often comprises a bottleneck that hinders the successful implementation of FEA. Soft tissues exhibit a nonlinear response and are commonly modeled using hyperelastic constitutive laws. In-vivo material parameter identification, for which standard mechanical tests (e.g., uniaxial tension and compression) are inapplicable, is commonly achieved using finite macro-indentation test. Due to the lack of analytical solutions, the parameters are commonly identified using inverse FEA (iFEA), in which simulated results and experimental data are iteratively compared. However, determining what data must be collected to accurately identify a unique parameter set remains unclear. This work investigates the sensitivities of two types of measurements: indentation force-depth data (e.g., measured using an instrumented indenter) and full-field surface displacements (e.g., using digital image correlation). To eliminate model fidelity and measurement-related errors, we employed an axisymmetric indentation FE model to produce synthetic data for four 2-parameter hyperelastic constitutive laws: compressible Neo-Hookean, and nearly incompressible Mooney-Rivlin, Ogden, and Ogden-Moerman models. For each constitutive law, we computed the objective functions representing the discrepancies in the reaction force, the surface displacement, and their combination, and visualized them for hundreds of parameter sets, spanning a representative range as found in the literature for the bulk soft tissue complex in human lower limbs. Moreover, we quantified three identifiability metrics, which provided insights into the uniqueness (or lack thereof) and the sensitivities. This approach provides a clear and systematic evaluation of the parameter identifiability, which is independent of the selection of the optimization algorithm and initial guesses required in iFEA. Our analysis indicated that the indenter's force-depth data, despite being commonly used for parameter identification, was insufficient for reliably and accurately identifying both parameters for all the investigated material models and that the surface displacement data improved the parameter identifiability in all cases, although the Mooney-Rivlin parameters remained poorly identifiable. Informed by the results, we then discuss several identification strategies for each constitutive model. Finally, we openly provide the codes used in this study, to allow others to further investigate the indentation problem according to their specifications (e.g., by modifying the geometries, dimensions, mesh, material models, boundary conditions, contact parameters, or objective functions).

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.

Can machine learning accelerate soft material parameter identification from complex mechanical test data?

Identifying the constitutive parameters of soft materials often requires heterogeneous mechanical test modes, such as simple shear. In turn, interpreting the resulting complex deformations necessitates the use of inverse strategies that iteratively call forward finite element solutions. In the past, we have found that the cost of repeatedly solving non-trivial boundary value problems can be prohibitively expensive. In this current work, we leverage our prior experimentally derived mechanical test data to explore an alternative approach. Specifically, we investigate whether a machine learning-based approach can accelerate the process of identifying material parameters based on our mechanical test data. Toward this end, we pursue two different strategies. In the first strategy, we replace the forward finite element simulations within an iterative optimization framework with a machine learning-based metamodel. Here, we explore both Gaussian process regression and neural network metamodels. In the second strategy, we forgo the iterative optimization framework and use a stand alone neural network to predict the entire material parameter set directly from experimental results. We first evaluate both approaches with simple shear experiments on blood clot, an isotropic, homogeneous material. Next, we evaluate both approaches against simple shear and uniaxial loading experiments on right ventricular myocardium, an anisotropic, heterogeneous material. We find that replacing the forward finite element simulations with metamodels significantly accelerates the parameter identification process with excellent results in the case of blood clot, and with satisfying results in the case of right ventricular myocardium. On the other hand, we find that replacing the entire optimization framework with a neural network yielded unsatisfying results, especially for right ventricular myocardium. Overall, the importance of our work stems from providing a baseline example showing how machine learning can accelerate the process of material parameter identification for soft materials from complex mechanical data, and from providing an open access experimental and simulation dataset that may serve as a benchmark dataset for others interested in applying machine learning techniques to soft tissue biomechanics.

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.

Anisotropic damage model for collagenous tissues and its application to model fracture and needle insertion mechanics

The analysis of tissue mechanics in biomedical applications demands nonlinear constitutive models able to capture the energy dissipation mechanisms, such as damage, that occur during tissue deformation. Furthermore, implementation of sophisticated material models in finite element models is essential to improve medical devices and diagnostic tools. Building on previous work toward microstructure-driven models of collagenous tissue, here we show a constitutive model based on fiber orientation and waviness distributions for skin that captures not only the anisotropic strain-stiffening response of this and other collagen-based tissues, but, additionally, accounts for tissue damage directly as a function of changes in the microstructure, in particular changes in the fiber waviness distribution. The implementation of this nonlinear constitutive model as a user subroutine in the popular finite element package Abaqus enables large-scale finite element simulations for biomedical applications. We showcase the performance of the model in fracture simulations during pure shear tests, as well as simulations of needle insertion into skin relevant to auto-injector design.

A computer based facial flaps simulator using projective dynamics

BACKGROUND AND OBJECTIVES

Interactive surgical simulation using the finite element method to model human skin mechanics has been an elusive goal. Mass-spring networks, while fast, do not provide the required accuracy.

METHODS

This paper presents an interactive, cognitive, facial flaps simulator based on a projective dynamics computational framework. Projective dynamics is able to generate rapid, stable results following changes to the facial soft tissues created by the surgeon, even in the face of sudden increases in skin resistance as its stretch limit is reached or collision between tissues occurs. Our prior work with the finite element method had been hampered by these considerations. Surgical tools are provided for; skin incision, undermining, deep tissue cutting, and excision. A spring-like "skin hook" is used for retraction. Spring-based sutures can be placed individually or automatically placed as a row between cardinal sutures.

RESULTS

Examples of an Abbe/Estlander lip reconstruction, a paramedian forehead flap to the nose, a retroauricular flap reconstruction of the external ear, and a cervico-facial flap reconstruction of a cheek defect are presented.

CONCLUSIONS

Projective dynamics has significant advantages over mass-spring and finite element methods as the physics backbone for interactive soft tissue surgical simulation.

Comparison of material properties of heel pad between adults with and without type 2 diabetes history: An in-vivo investigation during gait

OBJECTIVE

This study was aimed to compare the material properties of heel pad between diabetes patients and healthy adults, and investigate the impact of compressive loading history and length of diabetes course on the material properties of heel pad.

METHODS

The dual fluoroscopic imaging system (DFIS) and dynamic foot-ground contact pressure-test plate were used for measuring the material properties, including primary thickness, peak strain, peak stress, stiffness, viscous modulus and energy dissipation ratio (EDR), both at time zero and following continuous loading. Material properties between healthy adults and DM patients were compared both at time zero and following continuous weight bearing. After then, comparison between time-zero material properties and properties following continuous loading was performed to identify the loading history-dependent biomechanical behaviour of heel pad. Subgroup-based sensitivity analysis was then conducted to investigate the diabetes course (<10 years vs. ≥10 years) on the material properties of heel pad.

RESULTS

Ten type II DM subjects (20 legs), aged from 59 to 73 (average: 67.8 ± 4.9), and 10 age-matched healthy adults (20 legs), aged from 59 to 72 (average: 64.4 ± 3.4), were enrolled. Diabetes history was demonstrated to be associated with significantly lower primary thickness (t=3.18, p=0.003**), higher peak strain (t=2.41, p=0.021*), lower stiffness (w=283, p=0.024*) and lower viscous modulus (w=331, p<0.001***) at time zero, and significantly lower primary thickness (t=3.30, p=0.002**), higher peak strain (w=120, p=0.031*) and lower viscous modulus (t=3.42, p=0.002**) following continuous loading. The continuous loading was found to be associated with significantly lower primary thickness (paired-w=204, p<0.001***) and viscous modulus (paired-t=5.45, p<0.001***) in healthy adults, and significantly lower primary thickness (paired-w=206, p<0.001***) and viscous modulus (paired-t=7.47, p<0.001***) in diabetes group. No any significant difference was found when conducting the subgroup analysis based on length of diabetes course (<10 years vs. ≥10 years), but the regression analysis showed that the length of diabetes history was positively associated with the peak strain, at time zero (r=0.506, p<0.050) and following continuous loading (r=0.584, p<0.010).

CONCLUSIONS

Diabetes patients were found to be associated with decreased primary thickness and viscous modulus, and increased peak strain, which may contribute to the vulnerability of heel pad to injury and ulceration. Pre-compression history-dependent behaviour is observable in soft tissue of heel pad, with lowered primary thickness and viscous modulus.

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.

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.

Non-Invasive in Vivo Quantification of Directional Dependent Variation in Mechanical Properties for Human Skin

Skin is the body’s largest organ, and it shows non-linear and anisotropic behavior under the deformation. This behavior of the skin is due to the waviness and preferred orientation (in a particular direction) of collagen fibers. This preferred orientation of collagen fibers results in natural pre-tension and anisotropy of the skin. The knowledge of natural skin pre-tension and anisotropy is essential during incisions and surgery. The available suction-based devices quantify the anisotropy through the displacement field and cannot measure the stress-strain relation in particular directions. Therefore, in the current study, an in vivo full-field measurement suction apparatus was developed to measure the stress and strain of skin in all planar directions through a single experiment. First, this apparatus was tested on silicone substrates of known properties, and then it was used to test the skin of 12 human forearms. Further, to check the effect of hand stability on the measurements, the obtained results of the skin were compared with the results of a standard test performed in the same skin using a steady setup. The consistency between these two results confirms that the stability of the hand does not influence the measurements of skin properties. Furthermore, using the developed apparatus, the skin’s anisotropy and its relation with the Kraissl’s lines orientation was quantified by measuring the toe and linear moduli at an interval of one degree. The minimum and maximum values of the toe and linear moduli were 0.52 ± 0.09 and 0.59 ± 0.11 MPa, and 3.09 ± 0.47 and 5.52 ± 1.13 MPa, respectively. Also, the direction of maximum moduli was found almost similar to Kraissl’s lines’ orientation. These results confirm the contribution of skin pre-tension on the anisotropy of the skin. The present apparatus mimics the tissue expansion procedure, where observation of the test may be helpful in the selection of size and shape of the expander.

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.

Improving reconstructive surgery design using Gaussian process surrogates to capture material behavior uncertainty

To produce functional, aesthetically natural results, reconstructive surgeries must be planned to minimize stress as excessive loads near wounds have been shown to produce pathological scarring and other complications (Gurtner et al., 2011). Presently, stress cannot easily be measured in the operating room. Consequently, surgeons rely on intuition and experience (Paul et al., 2016; Buchanan et al., 2016). Predictive computational tools are ideal candidates for surgery planning. Finite element (FE) simulations have shown promise in predicting stress fields on large skin patches and in complex cases, helping to identify potential regions of complication. Unfortunately, these simulations are computationally expensive and deterministic (Lee et al., 2018a). However, running a few, well selected FE simulations allows us to create Gaussian process (GP) surrogate models of local cutaneous flaps that are computationally efficient and able to predict stress and strain for arbitrary material parameters. Here, we create GP surrogates for the advancement, rotation, and transposition flaps. We then use the predictive capability of these surrogates to perform a global sensitivity analysis, ultimately showing that fiber direction has the most significant impact on strain field variations. We then perform an optimization to determine the optimal fiber direction for each flap for three different objectives driven by clinical guidelines (Leedy et al., 2005; Rohrer and Bhatia, 2005). While material properties are not controlled by the surgeon and are actually a source of uncertainty, the surgeon can in fact control the orientation of the flap with respect to the skin's relaxed tension lines, which are associated with the underlying fiber orientation (Borges, 1984). Therefore, fiber direction is the only material parameter that can be optimized clinically. The optimization task relies on the efficiency of the GP surrogates to calculate the expected cost of different strategies when the uncertainty of other material parameters is included. We propose optimal flap orientations for the three cost functions and that can help in reducing stress resulting from the surgery and ultimately reduce complications associated with excessive mechanical loading near wounds.

A biphasic multilayer computational model of human skin

The present study investigates the layer-specific mechanical behavior of human skin. Motivated by skin’s histology, a biphasic model is proposed which differentiates between epidermis, papillary and reticular dermis, and hypodermis. Inverse analysis of ex vivo tensile and in vivo suction experiments yields mechanical parameters for each layer and predicts a stiff reticular dermis and successively softer papillary dermis, epidermis and hypodermis. Layer-specific analysis of simulations underlines the dominating role of the reticular dermis in tensile loading. Furthermore, it shows that the observed out-of-plane deflection in ex vivo tensile tests is a direct consequence of the layered structure of skin. In in vivo suction experiments, the softer upper layers strongly influence the mechanical response, whose dissipative part is determined by interstitial fluid redistribution within the tissue. Magnetic resonance imaging-based visualization of skin deformation in suction experiments confirms the deformation pattern predicted by the multilayer model, showing a consistent decrease in dermal thickness for large probe opening diameters.

On the Reliability of Suction Measurements for Skin Characterization

in vivo skin characterization methods were shown to be useful in the detection of microstructural alterations of the dermis due to skin diseases. Specifically, the diagnostic potential of skin suction has been widely explored, yet measurement uncertainties prevented so far its application in clinical assessment. In this work, we analyze specific factors influencing the reliability of suction measurements. We recently proposed a novel suction device, called Nimble, addressing the limitations of existing instruments, and applied it in clinical trials quantifying mechanical differences between healthy skin and scars. Measurements were performed with the commercial device Cutometer and with the new device. A set of new suction measurements was carried out on scar tissue and healthy skin, and FE-based inverse analysis was applied to determine corresponding parameters of a hyperelastic-viscoelastic material model. FE simulations were used to rationalize differences between suction protocols and to analyze specific factors influencing the measurement procedure. Tissue stiffness obtained from Cutometer measurements was significantly higher compared to the one from Nimble measurements, which was shown to be associated with the higher deformation levels in the Cutometer and the nonlinear mechanical response of skin. The effect of the contact force exerted on skin during suction measurements was quantified, along with an analysis of the effectiveness of a corresponding correction procedure. Parametric studies demonstrated the inherently higher sensitivity of displacement- over load-controlled suction measurements, thus rationalizing the superior ability of the Nimble to distinguish between tissues.

How cosmetic tightening products modulate the biomechanics and morphology of human skin

The active and passive mechanical behavior of a cosmetic tightening product for skin anti-aging is investigated based on a wide range of in vivo and in vitro measurements. The experimental data are used to inform a numerical model of the attained cosmetic effect, which is then implemented in a commercial finite-element framework and used to analyze the mechanisms that regulate the biomechanical interaction between the native tissue and the tightening film. Such a film reduces wrinkles and enhances skin consistency by increasing its stiffness by 48-107% and reducing inelastic, non-recoverable deformations (-47%). The substrate deformability influences both the extent of tightening and the reduction of wrinkle amplitude. The present findings allow, for the first time, to rationalize the mechanisms of action of cosmetic products with a tightening action and provide quantitative evidence for further optimization of this fascinating class of biomaterials.

Personalized Computational Models of Tissue-Rearrangement in the Scalp Predict the Mechanical Stress Signature of Rotation Flaps

OBJECTIVE

To elucidate the mechanics of scalp rotation flaps through 3D imaging and computational modeling. Excessive tension near a wound or sutured region can delay wound healing or trigger complications. Measuring tension in the operating room is challenging, instead, noninvasive methods to improve surgical planning are needed.

DESIGN

Multi-view stereo allows creation of 3D patient-specific geometries based on a set of photographs. The patient-specific 3D geometry is imported into a finite element (FE) platform to perform a virtual procedure. The simulation is compared with the clinical outcome. Additional simulations quantify the effect of individual flap parameters on the resulting tension distribution.

PARTICIPANTS

Rotation flaps for reconstruction of scalp defects following melanoma resection in 2 cases are presented. Rotation flaps were designed without preoperative FE preparation.

MAIN OUTCOME MEASURE

Tension distribution over the operated region.

RESULTS

The tension from FE shows peaks at the base and distal ends of the scalp rotation flap. The predicted geometry from the simulation aligns with postoperative photographs. Simulations exploring the flap design parameters show variation in the tension. Lower tensions were achieved when rotation was oriented with respect to skin tension lines (horizontal tissue fibers) and smaller rotation angles.

CONCLUSIONS

Tension distribution following rotation of scalp flaps can be predicted through personalized FE simulations. Flaps can be designed to reduce tension using FE, which may greatly improve the reliability of scalp reconstruction in craniofacial surgery, critical in complex cases when scalp reconstruction is essential for coverage of hardware, implants, and/or bone graft.

Activin-mediated alterations of the fibroblast transcriptome and matrisome control the biomechanical properties of skin wounds

Matrix deposition is essential for wound repair, but when excessive, leads to hypertrophic scars and fibrosis. The factors that control matrix deposition in skin wounds have only partially been identified and the consequences of matrix alterations for the mechanical properties of wounds are largely unknown. Here, we report how a single diffusible factor, activin A, affects the healing process across scales. Bioinformatics analysis of wound fibroblast transcriptome data combined with biochemical and histopathological analyses of wounds and functional in vitro studies identify that activin promotes pro-fibrotic gene expression signatures and processes, including glycoprotein and proteoglycan biosynthesis, collagen deposition, and altered collagen cross-linking. As a consequence, activin strongly reduces the wound and scar deformability, as identified by a non-invasive in vivo method for biomechanical analysis. These results provide mechanistic insight into the roles of activin in wound repair and fibrosis and identify the functional consequences of alterations in the wound matrisome at the biomechanical level.

The relationship between histopathology, gene expression, and biochemical and mechanical properties of wounds is largely unknown. Here, the authors show that activin A alters wound healing at multiple levels by promoting pro-fibrotic gene expression and matrix deposition, thereby affecting biomechanical properties of skin wounds.

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.

In vivo stress relaxation of human scalp

OBJECTIVE

Conduct a first in vivo study on the large deformation stress relaxation behavior of the human scalp.

METHODS

This study was conducted during Mohs micrographic surgery of the scalp of 14 patients aged 59-90 with wounds initially ranging from 9 to 41 mm wide. The initial wound diameter was measured under zero applied force. Then, the force required to close each wound using a single size 1 nylon suture and a SUTUREGARD suture retention device was measured, after which the suture was then locked in the retention device at fixed displacement. At time points of 300 s, 600 s, and 1800 s, the suture retention device was released, and the wound opening was again recorded at zero force, and the force required to close the wound was recorded.

RESULTS

The average wound closure force relaxed by 44% and 65% after 300 s and 1800 s, respectively. Average wound width decreased 30% and 42%, after 300 s and 1800 s, respectively, due to creep deformation. Furthermore, all wounds relaxed to be below 15 N of closure force after 600 s, which is considered the maximum clinically acceptable force. A relaxation time of ∼270 s and a threshold force for creep of ∼5 N was found.

SIGNIFICANCE

Results of this study provide the first quantitative clinical guidance for efficient scalp closure of large wounds by creep deformation and stress relaxation. Furthermore, the methodology developed here can be used as a basis for future in vivo studies of the stress relaxation and creep deformation of human scalp, which in turn can provide data for the development and validation of constitutive models for scalp deformation.

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.

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.

Surface deformation tracking and modelling of soft materials

Many computer vision algorithms have been presented to track surface deformations, but few have provided a direct comparison of measurements with other stereoscopic approaches and physics-based models. We have previously developed a phase-based cross-correlation algorithm to track dense distributions of displacements over three-dimensional surfaces. In the present work, we compare this algorithm with one that uses an independent tracking system, derived from an array of fluorescent microspheres. A smooth bicubic Hermite mesh was fitted to deformations obtained from the phase-based cross-correlation data. This mesh was then used to estimate the microsphere locations, which were compared to stereo reconstructions of the microsphere positions. The method was applied to a 35 mm × 35 mm × 35 mm soft silicone gel cube under indentation, with three square bands of microspheres placed around the indenter tip. At an indentation depth of 4.5 mm, the root-mean-square (RMS) differences between the reconstructed positions of the microspheres and their identified positions for the inner, middle, and outer bands were 60 µm, 20 µm, and 19 µm, respectively. The usefulness of the strain-tracking data for physics-based finite element modelling of large deformation mechanics was then demonstrated by estimating a neo-Hookean stiffness parameter for the gel. At the optimal constitutive parameter estimate, the RMS difference between the measured microsphere positions and their finite element model-predicted locations was 143 µm.

A novel ultra-light suction device for mechanical characterization of skin

Suction experiments have been extensively applied for skin characterization. In these tests the deformation behavior of superficial tissue layers determines the elevation of the skin surface observed when a predefined negative (suction) pressure history is applied. The ability of such measurements to differentiate between skin conditions is limited by the variability of the elevation response observed in repeated experiments. The scatter was shown to be associated with the force exerted by the observer when holding the instrument against the skin. We have developed a novel suction device and a measurement procedure aiming at a tighter control of mechanical boundary conditions during the experiments. The new device weighs only 3.5 g and thus allows to minimize the force applied on the skin during the test. In this way, it is possible to reliably characterize the mechanical response of skin, also in case of low values of suction pressure and deformation. The influence of the contact force is analyzed through experiments on skin and synthetic materials, and rationalized based on corresponding finite element calculations. A comparative study, involving measurements on four body locations in two subjects by three observers, showed the good performance of the new procedure, specific advantages, and limitations with respect to the Cutometer^®, i.e. the suction device most widely applied for skin characterization. As a byproduct of the present investigation, a correction procedure is proposed for the Cutometer measurements, which allows to partially compensate for the influence of the contact force. The characteristics of the new suction method are discussed in view of future applications for diagnostic purposes.

DISPOSABLE SYSTEM FOR IN-VIVO MECHANICAL CHARACTERIZATION OF SOFT TISSUES BASED ON VOLUME MEASUREMENT

In-vivo characterization of soft tissues is a key step toward biomechanical simulation and planning of intra-operative assisted surgery. To achieve this, aspiration method is a standard technique: tissue is aspirated through a hole while measuring the pressure and associated apex height. An inverse problem is then solved to identify the material mechanical properties. In the literature, the apex height is usually measured using a camera, which induces design difficulties, in particular in regards on the required sterilization process for in-vivo measurements. In this paper, the idea is to replace the apex height optical measurement by the measurement of the aspirated tissue volume. The proposed method enables to reduce the system head to a simple tube: sterilizations becomes easy and the system is disposable after use. The proposed system is thus the simplest, lightest and cheapest one could achieve. It is also to the authors knowledge the first time ever in aspiration method that the aspired volume is the extracted data. As the data signal-to-noise ratio is the main factor impacting any applied inverse method when extracting the mechanical properties, the aim of this work is to assess and compare the experimental signal-to-noise ratio in the raw volume measurements obtained either optically or with the method proposed. Explicit results of inverse methods using volumes as input data are not presented in this paper for concision purpose. The effects on accuracy of various experimental parameters has been investigated and quantified: the volume measurement has proved to present a same order or even better signal-to-noise ratio compared to optical measurements.

Improving tissue expansion protocols through computational modeling

Tissue expansion is a common technique in reconstructive surgery used to grow skin in vivo for correction of large defects. Despite its popularity, there is a lack of quantitative understanding of how stretch leads to growth of new skin. This has resulted in several arbitrary expansion protocols that rely on the surgeon's personal training and experience rather than on accurate predictive models. For example, choosing between slow or rapid expansion, or small or large inflation volumes remains controversial. Here we explore four tissue expansion protocols by systematically varying the inflation volume and the protocol duration in a porcine model. The quantitative analysis combines three-dimensional photography, isogeometric kinematics, and finite growth theory. Strikingly, all four protocols generate similar peak stretches, but different growth patterns: Smaller filling volumes of 30 ml per inflation did not result in notable expander-induced growth neither for the short nor for the long protocol; larger filling volumes of 60 ml per inflation trigger skin adaptation, with larger expander-induced growth in regions of larger stretch, and more expander-induced growth for the 14-day compared to the 10-day expansion protocol. Our results suggest that expander-induced growth is not triggered by the local stretch alone. While stretch is clearly a driver for growth, the local stretch at a given point is not enough to predict the expander-induced growth at that location. From a clinical perspective, our study suggests that longer expansion protocols are needed to ensure sufficient growth of sizable skin patches.

Multi-view stereo in the operating room allows prediction of healing complications in a patient-specific model of reconstructive surgery

Excessive mechanical stress leads to wound healing complications following reconstructive surgery. However, this knowledge is not easily applicable in clinical scenarios due to the difficulty in measuring stress contours during complex tissue rearrangement procedures. Computational tools have been proposed as an alternative to address this need, but obtaining patient specific geometries with an affordable and flexible setup has remained a challenge. Here we present a methodology to determine the stress contours from a reconstructive procedure on a patient-specific finite element model based on multi-view stereo (MVS). MVS is a noninvasive technology that allows reconstruction of 3D geometries using a standard digital camera, making it ideal for the operating room. Finite element analysis can then be used on the patient-specific geometry to perform a virtual surgery and predict regions at risk of complications. We applied our approach to the case of a 7-year-old patient who was treated to correct a cranial contour deformity and resect two large areas of scalp scarring. The simulation showed several zones of high stress, particularly near the suture lines at the distal ends of the flaps. The patient did show delayed healing and partial flap tip necrosis at one of such predicted regions at the 30-day follow up visit. Our results further establish the application of computational tools in individualized medical scenarios to advance preoperative planing and anticipate regions of concern immediately after surgery.

Multidirectional In Vivo Characterization of Skin Using Wiener Nonlinear Stochastic System Identification Techniques

A triaxial force-sensitive microrobot was developed to dynamically perturb skin in multiple deformation modes, in vivo. Wiener static nonlinear identification was used to extract the linear dynamics and static nonlinearity of the force-displacement behavior of skin. Stochastic input forces were applied to the volar forearm and thenar eminence of the hand, producing probe tip perturbations in indentation and tangential extension. Wiener static nonlinear approaches reproduced the resulting displacements with variances accounted for (VAF) ranging 94-97%, indicating a good fit to the data. These approaches provided VAF improvements of 0.1-3.4% over linear models. Thenar eminence stiffness measures were approximately twice those measured on the forearm. Damping was shown to be significantly higher on the palm, whereas the perturbed mass typically was lower. Coefficients of variation (CVs) for nonlinear parameters were assessed within and across individuals. Individual CVs ranged from 2% to 11% for indentation and from 2% to 19% for extension. Stochastic perturbations with incrementally increasing mean amplitudes were applied to the same test areas. Differences between full-scale and incremental reduced-scale perturbations were investigated. Different incremental preloading schemes were investigated. However, no significant difference in parameters was found between different incremental preloading schemes. Incremental schemes provided depth-dependent estimates of stiffness and damping, ranging from 300 N/m and 2 Ns/m, respectively, at the surface to 5 kN/m and 50 Ns/m at greater depths. The device and techniques used in this research have potential applications in areas, such as evaluating skincare products, assessing skin hydration, or analyzing wound healing.

Location-specific mechanical response and morphology of facial soft tissues

The facial tissue of 9 healthy volunteers (m/f; age: 23-60y) is characterized at three different locations using a procedure combining suction measurements and 18MHz ultrasound imaging. The time-dependent and multilayered nature of skin is accounted for by adopting multiple loading protocols which differ with respect to suction probe opening size and rate of tissue deformation. Over 700 suction measurements were conducted and analyzed according to location-specific mechanical and morphological characteristics. All corresponding data are reported and made available for facial tissue analysis and biomechanical modeling. Higher skin stiffness is measured at the forehead in comparison to jaw and parotid; these two regions are further characterized by lower creep deformation. Thicker tissue regions display a tendency towards a more compliant and less dissipative response. Comparison of superficial layer thickness and corresponding mechanical measurements suggests that connective tissue density determines the resistance to deformation in suction experiments.

Mathematical and computational modelling of skin biophysics: a review

The objective of this paper is to provide a review on some aspects of the mathematical and computational modelling of skin biophysics, with special focus on constitutive theories based on nonlinear continuum mechanics from elasticity, through anelasticity, including growth, to thermoelasticity. Microstructural and phenomenological approaches combining imaging techniques are also discussed. Finally, recent research applications on skin wrinkles will be presented to highlight the potential of physics-based modelling of skin in tackling global challenges such as ageing of the population and the associated skin degradation, diseases and traumas.

Quantification of Strain in a Porcine Model of Skin Expansion Using Multi-View Stereo and Isogeometric Kinematics

Tissue expansion is a popular technique in plastic and reconstructive surgery that grows skin in vivo for correction of large defects such as burns and giant congenital nevi. Despite its widespread use, planning and executing an expansion protocol is challenging due to the difficulty in measuring the deformation imposed at each inflation step and over the length of the procedure. Quantifying the deformation fields is crucial, as the distribution of stretch over time determines the rate and amount of skin grown at the end of the treatment. In this manuscript, we present a method to study tissue expansion in order to gain quantitative knowledge of the deformations induced during an expansion process. This experimental protocol incorporates multi-view stereo and isogeometric kinematic analysis in a porcine model of tissue expansion. Multi-view stereo allows three-dimensional geometric reconstruction from uncalibrated sequences of images. The isogeometric kinematic analysis uses splines to describe the regional deformations between smooth surfaces with few mesh points. Our protocol has the potential to bridge the gap between basic scientific inquiry regarding the mechanics of skin expansion and the clinical setting. Eventually, we expect that the knowledge gained with our methodology will enable treatment planning using computational simulations of skin deformation in a personalized manner.