Fracture Energy Research Articles

Multi-Modal Sensing Ionogels with Tunable Mechanical Properties and Environmental Stability for Aquatic and Atmospheric Environments.

Ionogels have garnered significant interest due to their great potential in flexible iontronic devices. However, their limited mechanical tunability and environmental intolerance have posed significant challenges for their integration into next-generation flexible electronics in different scenarios. Herein, the synergistic effect of cation-oxygen coordination interaction and hydrogen bonding is leveraged to construct a 3D supramolecular network, resulting in ionogels with tunable modulus, stretchability, and strength, achieving an unprecedented elongation at break of 10800%. Moreover, the supramolecular network endows the ionogels with extremely high fracture energy, crack insensitivity, and high elasticity. Meanwhile, the high environmental stability and hydrophobic network of the ionogels further shield them from the unfavorable effects of temperature variations and water molecules, enabling them to operate within a broad temperature range and exhibit robust underwater adhesion. Then, the ionogel is assembled into a wearable sensor, demonstrating its great potential in flexible sensing (temperature, pressure, and strain) and underwater signal transmission. This work can inspire the applications of ionogels in multifunctional sensing and wearable fields.

Fracture behavior of 3D-printed thermoplastics—A dilemma of fracture toughness versus fracture energy

The Material Extrusion (ME) technique in 3D printing facilitates the production of thermoplastic materials with diverse cellular or lattice-like infill geometries, enabling the creation of lightweight, high-performance materials. This study investigates the influence of infill geometry and relative density on the fracture behavior of 3D-printed thermoplastics produced through ME. Polylactic acid (PLA) is used as the model material, with unit cell geometries-square, triangle, and Kagome-and relative densities ranging from 0.4 to 0.8 explored. Tensile tests are first conducted to determine in-plane elastic moduli in two orthogonal directions, accounting for unit cell anisotropy. These results are employed in effective fracture energy calculations. Compact tension fracture tests follow, quantitatively characterize crack growth characteristics in both directions. In both tests, digital image correlation method is used to measure full field strain distributions. Three significant findings emerge. First, there is a power scaling relationship between Elastic Modulus and relative density across all unit cells, with Kagome exhibiting the highest correlation constant. Triangle and Kagome unit cells show negligible orthotropic responses in perpendicular directions for varying relative densities. Second, the relationship between nondimensional fracture toughness and relative density follows a similar scaling law, with minor variations depending on unit cell type. The correlation factor slightly dependens on unit cell geometry but remains around 3 across the the examined relative density range. Third, considering the dissipation of inelastic fracture energy, both triangle and Kagome unit cell geometries demonstrate an order of magnitude higher effective fracture energy (i.e., crack growth resistance) compared to the square unit cell. This is associated with toughening mechanisms such as crack deflection and bifurcation. These findings highlight the disparity between fracture toughness and effective fracture energy in assessing the fracture resistance of 3D-printed thermoplastics, clearly demonstrating that fracture toughness alone is insufficient for a comprehensive assessment. The proposed scaling laws provide predictive insights into the fracture toughness of polymeric materials produced via the ME method.

Mechanical Properties of Permeable Concrete Reinforced with Polypropylene Fibers for Different Water–Cement Ratios

Permeable concrete is a material that allows water filtration, reduces surface runoff, and maintains the natural water cycle. Previous studies have shown that its mechanical properties, particularly its compressive and flexural tensile strengths, are generally lower than those of conventional concrete, with significant variability observed among similar tests. This study investigates the compressive strength, flexural strength, and permeability of polypropylene fiber-reinforced permeable concrete specimens at two water–cement ratios (0.30 and 0.35). The mix design was conducted using the ACI 522R-10 method. Forty-eight cylinders measuring 200 mm × 100 mm were fabricated for permeability and compression tests. Additionally, 12 beams measuring 100 mm × 100 mm × 350 mm were produced and subjected to simple flexural testing in accordance with ASTM C78 guidelines. Compressive strength versus permeability and load versus deflection graphs were plotted, and the fracture energy was calculated for various deflections. The results indicate that the addition of fibers increased permeability and tensile strength but decreased compressive strength. Furthermore, an increase in the water–cement ratio led to higher compressive and flexural tensile strengths.

Research on the Gap Effect of Circular Concrete-Filled Steel Tubes Using the Improved Cohesive Zone Model

Understanding the influence of gap distribution characteristics on the mechanical properties of circular concrete-filled steel tubes (CCFSTs) under bending load is important for stability and support design in engineering projects. In this study, the improved cohesive zone model considering friction was used to describe the mechanical behavior of mortar interfaces. Meanwhile, the concrete damage plastic model and isotropic elastoplastic model were applied for core concrete and steel tubes. The improved cohesive zone model has a unified potential function that governs the Mode I and Mode II failure processes of mortar interfaces to realize the mechanical interaction between concrete and steel. A smooth frictional function was utilized in the elastic stage to calculate the accurate frictional effect. Furthermore, the capability of the model in addressing unloading and reloading was verified, and the fracture energy varied accordingly during the cyclic loading. Then, the mechanical response of CCFSTs was investigated under bending loads by setting different gap sizes and angles between the gap and loading direction. The results show that under three-point bending, the equivalent plastic strains at the middle part of CCFSTs are much larger and the peak bearing forces are much lower than the other degrees when the angles between the coronal gap axis and loading direction equal 0° and 180°. In addition, the order of the peak bearing forces, from highest to lowest, is when the height of the coronal-cap gap increases from 0.0 mm to 2.5 mm, 5.0 mm, and 7.5 mm. The significant effect makes it inappropriate to ignore the weakening of the structural performance caused by coronal gaps in structural design.

Tensile performance prediction of CFRPs with voids using multiscale analysis and neural networks

The overall work provides a novel multiscale modeling strategy applicable to the mechanical response prediction of arbitrarily laid porous laminates under uniaxial tension by establishing representative defective volumes (RDVs) for the inter-ply and intra-ply regions. A representative defective volume modeling method based on void characterization results using X-ray computed tomography for the interlaminar region was proposed. Elastic properties, tensile strength, and fracture energies of porous and poreless microscale models were obtained and used in the macroscale model. Defective properties were assigned to partial elements in the macroscale model according to the preset distribution mode. The Extended Finite Element Method was taken to simulate the failure process. Tensile properties of porous laminates with stacking sequences of [90]8, [90/0/90/0]2S, and [45/90/-45/0]2S were experimentally tested and numerically simulated. The predicted tensile properties and crack propagation behavior coincide well with the test results. Besides, the effects of micro-voids and macroscopic void distribution were discussed in detail. Further, a back-propagation neural network was used to achieve a fast prediction of the tensile properties and save computational costs. The computation time (<1 s) greatly improved compared with multiscale modeling (>38 h). This work provides an efficient tool for composite performance prediction, design, and optimization.

Dynamic damage behavior of auxetic textile reinforced concrete under impact loading

To investigate the dynamic compression behavior of Textile-reinforced concrete (TRC) under impact load, low speed and high speed impact tests were performed using a ball drop test and a split Hopkinson pressure bar (SHPB) with ϕ50mm, respectively. In view of the advantages of auxetic materials in impact resistance and energy absorption, the auxetic fiber preforms with negative Poisson's ratio effect were woven into meshes and added to the cement-based materials to obtain the auxetic textile-reinforced concrete. Thus, the auxetic textile reinforced cement-based material is obtained. In this paper, TRC with different textile types, laying layers and surface treatments were prepared, and impact tests with different impact speeds were carried out to evaluate the impact resistance of TRC. The impact damage behavior of TRC was analyzed through the method of digital scattering analysis, and then the impact damage mechanism of TRC was investigated. The results show that the auxetic effect of auxetic textiles can still be stabilized under impact loads of different velocities, which provides better control of deformation and energy consumption in the matrix material. Under different impact heights of falling balls, the impact fracture energy of cementitious materials reinforced with auxetic textiles was increased by up to 22.5 and 20 times, respectively, compared with that of blank samples.

Shrinkage cracking law and anti-crack inverse design in early-age concrete: A novel perspective on the development of crack resistance properties

Early-age cracking in massive concrete has long attracted research focus, essentially governed by the game between crack driving and resistance. However, cracking resistance aspect received inadequate attention. This study examines the impact of crack resistance properties development processes on shrinkage cracking law in early-age concrete, employing chemo-thermo-mechanical coupling phase-field model. Two quantitative evaluation indicators, namely the index of damage development and average growth rate of maximum damage, are introduced to characterize the process of damage and cracking. The findings demonstrate that the development processes of elastic modulus and tensile strength exert a considerable influence on shrinkage cracking process in early-age concrete, whereas the impact of fracture energy development process is insignificant. Multifactor analysis revealed substantial coupling effects among parameters. Additionally, an anti-crack reverse design method is proposed based on the global optimization analysis, which can guide the optimization design of shrinkage crack prevention and control in early-age concrete from both perspectives of crack driving and crack resistance.

Phase-field finite element modelling of creep crack growth in martensitic steels

Creep fracture presents a major concern for structural materials and critical components operating at elevated temperatures, thus requiring effective computational models. This study presents a phase-field framework for modelling creep crack growth and fracture behaviour of modern high-temperature materials such as Creep Strength Enhanced Ferritic (CSEF) martensitic steels. The model was formulated using the thermodynamic principles of the variational phase field theory of fracture and considering some physical aspects of creep fracture. Within the modelling framework, a dissipation potential dependent on creep damage is introduced to capture the effect of creep cavitation on the fracture energy and the creep crack growth resistance of the solid in a phenomenological manner. An elastoplastic power-law creep model is coupled to the phase-field formulations to account for the non-linear deformation processes ahead of the crack tip due to inelastic deformations, as well as their contribution to fracture at high temperatures. The capability of the proposed model is assessed against experimental data from compact tension (CT) creep tests conducted on P91 and P92 steels. Good agreement was obtained between the FE-predicted creep crack growth behaviour and the experimental measurements, showcasing the model’s feasibility. Further, numerical experiments were conducted using the proposed model to elucidate some key aspects influencing the fracture behaviour of martensitic steels. The proposed computational framework not only demonstrated good capability but was also able to offer improved mechanistic insights into the influence of material tendency to develop creep cavities on crack growth behaviour. This work contributes valuable insights into understanding the fracture process of CSEF steels at elevated temperatures and further demystifies the role of creep ductility.

Criterion for hydraulic fracture propagation behaviour at coal measure composite reservoir interface based on energy release rate theory

In this study, a discriminative model of fracture expansion through a composite reservoir interface was established. During fracture transformation, the model, based on the theory of energy release rate at the fracture tip, helps clarify the critical conditions and assess the law of hydraulic fracture expansion through coal composite reservoir interfaces. This model considers shear and tensile processes during fracturing at the coal–rock interface and the effects of interface material property differences on fracture extension. In addition, fracture extension pattern differentiation is achieved by comparing the hydraulic fracture tip penetration and bending energy release rates (Gp/Gd) with the ratio of reservoir fracture and interface fracture energies (ΓR/ΓF). The accuracy of this discriminant criterion was verified by comparing the predicted results of the discriminant model with fracturing test results under identical conditions. The effects of various factors on Gp/Gd and hydraulic fracture penetration extension were investigated by using the sand-mudstone composite reservoir of the Shibox Formation in the Linxing Block, China, as the geological background. The results indicate that for certain values of ΓR/ΓF, Gp/Gd increases logarithmically with an increase in water pressure in the wellbore. Moreover, the difference in ground stress and elastic modulus between the layers exponentially increases with increasing Poisson’s ratio difference and the angle between fracture and interface and decreases logarithmically with increasing fracture height. These results indicate that hydraulic fractures are likely to penetrate from reservoirs with a high elastic modulus and Poisson’s ratio into reservoirs with a low elastic modulus and Poisson’s ratio. Furthermore, hydraulic fractures are likely to penetrate formations with large differences in pinch angle and ground stress. Under specific geological conditions, high water pressure in the wellbore and small seam height are favourable for fracture penetration.

Squid-Inspired Anti-Salt Skin-Like Elastomers With Superhigh Damage Resistance for Aquatic Soft Robots.

Cephalopod skins evolve multiple functions in response to environmental adaptation, encompassing nonlinear mechanoreponse, damage tolerance property, and resistance to seawater. Despite tremendous progress in skin-mimicking materials, the integration of these desirable properties into a single material system remains an ongoing challenge. Here, drawing inspiration from the structure of reflectin proteins in cephalopod skins, a long-term anti-salt elastomer with skin-like nonlinear mechanical properties and extraordinary damage resistance properties is presented. Cation-π interaction is incorporated to induce the geometrically confined nanophases of hydrogen bond domains, resulting in elastomers with exceptional true tensile strength (456.5 ± 68.9MPa) and unprecedently high fracture energy (103.7 ± 45.7kJ m-2). Furthermore, the cation-π interaction effectively protects the hydrogen bond domains from corrosion by high-concentration saline solution. The utilization of the resultant skin-like elastomer has been demonstrated by aquatic soft robotics capable of grasping sharp objects. The combined advantages render the present elastomer highly promising for salt enviroment applications, particularly in addressing the challenges posed by sweat, in vivo, and harsh oceanic environments.

Mechanisms for Recycling E-Glass Waste to Enhance The Fracture Energy of 3d Printed Structural Dune Sand Mortar

Mechanisms for Recycling E-Glass Waste to Enhance The Fracture Energy of 3d Printed Structural Dune Sand Mortar

Mechanism of Rock Mass Detachment Using Undercutting Anchors: A Numerical Finite Element Method (FEM) Analysis.

Undercutting anchors are structural elements used in construction and geotechnics to stabilize both structures and soils. Their main applications include stabilizing slopes and embankments, reinforcing foundations, and providing support during tunnel construction and other underground works. The authors propose the use of these anchors in rock mass detachment technology. This article presents a comprehensive analysis of the mechanism behind rock mass detachment using an undercutting anchor. Particular attention is given to the influence of parameters such as the fracture energy of the medium and the coefficient of friction between the medium and the anchor head on the detachment process of rock elements during anchor expansion in the drilled hole. Numerical FEM analysis was employed to model the effect of changes in the shape and size of failure cones under varying simulation conditions. The discussed problem is crucial for evaluating the effectiveness of this anchor design under non-standard conditions, particularly in the unconventional destruction of rock media.

Effect of key design parameters on the physical / mechanical properties of coarse aggregated UHPC and its structural response subjected to bending

This study experimentally investigates the effect of key design parameters on the physical, mechanical and structural performance of coarse aggregated ultra-high performance concrete (CA-UHPC), with an aim to provide design guidance for on-site engineering applications. The variables of maximum aggregate size, aggregate content, fiber dosage and curing condition are selected for investigation. The results showed that increasing aggregate size does not significantly affect the compressive strength, but somewhat increases the elastic modulus and fracture energy; the increase of aggregate content results in a reasonable increase of compressive strength and elastic modulus. In particular, the 480 kg/m3 coarse aggregated UHPC shows a compressive strength of ∼150 MPa, elastic modulus of ∼50 GPa, autogenous shrinkage of <200 µε, and outstanding durability. The structural specimens fabricated with coarse aggregated UHPC exhibit excellent load-carrying capacity and ductility when subjected to four-point bending. More importantly, it was found that the effect of curing condition on the structural performance is insignificant. The results of this study can provide insights into the development of advanced UHPC materials and structures and promote their applications in practice.

Adhesive and cohesive fracture of blood clots: Experiments and modeling

Blood clots represent living materials composed of a polymer network and an abundance of cells. They might fracture within the bulk material of the clot (cohesive fracture), at the interface between the clot and the surrounding tissue (adhesive fracture), or through a combination of both modes (hybrid fracture). The clot fracture within vascular systems and injury sites could lead to life-threatening conditions. Despite the significance, understanding and modeling the fracture behaviors of blood clots, including their dependence on mechanical loading and cellular components, remain in a nascent stage. In this study, we employ an integrated experimental-computational approach to comprehensively investigate the fracture behaviors of bovine blood clots. We explore various mechanical factors, substrates, and cellular components such as red blood cells (RBCs) and platelets. Our findings reveal that among various tissue substrates, blood clots exhibit the highest interfacial adhesion energy with muscle, and the lowest to the inner arterial lining, consistent with their biological function. Both interfacial adhesion energy and bulk fracture energy are rate-dependent, although they exhibit different dependencies. Also, RBCs and platelets have different effects on clot fracture. An increase in RBC content tends to toughen both adhesion and fracture of blood clots. However, an increase in platelet content enhances interfacial adhesion energy but lowers the bulk fracture energy. The platelet content also governs the shift from adhesive fracture to hybrid fracture. To model clot fracture, we developed two finite element models incorporating a coupled cohesive-zone and Mullins-effect approach to simulate pure shear fracture and peeling of blood clots. These models, validated through experimental data, elucidate the interplay between intrinsic fracture toughness, interfacial strength, and bulk energy dissipation during clot fracture. This study significantly advances our understanding of clot mechanics, providing valuable insights into the mechanics of similar living materials and the management of clot-related disorders such as hemorrhage and thrombosis.

Revamp of the eco-friendly grouting materials by mixing basalt fiber for earthen sites conservation: Compressive performance and constitutive models

Reinforcement techniques, including anchorage and sealing with grout, are frequently employed to protect earthen sites and minimize the impact of fissures, the effectiveness of these approaches relies heavily on the mechanical properties of grout. To enhance the strength and ductility of the conventional grout, an eco-friendly composite grouting material mixture basalt fibers of varying lengths and contents is prepared. The DIC-based UCS tests were conducted to evaluate its compressive performance. Constitutive models of the basalt fiber-reinforced grout, including elastic perfectly-plastic (EPP) and median nonlinear (MNL) constitutive models, were established respectively based on the equal energy rule and cross-sectional data analysis for simplification application and precision evaluation purposes. The results show that the mixture of basalt fibers resulted in a distinct semi-brittle failure mode for the grout, with a significant increase in its compressive strength, ductility, compressive fracture energy, and residual bearing capacity. However, the elastic modulus will be lower when the fibers are relatively longer. Considering the differences in the relations between performance parameters and fiber content, and length, the highest compressive strength was reached at the fiber length of 6 mm and the content of 0.2 wt%, while maximum ductility was achieved at 18 mm and 0.4 wt%. The fracture energy exhibits an approximate positive correlation with fiber content. In addition, the EPP and MNL constitutive models can effectively characterize the compressive mechanical behavior of semi-brittle materials to satisfy different requirements on the precision of engineering analysis. This study contributes valuable insights into the development of high-ductility grouting materials of earthen sites, and the establishment of constitutive models for semi-brittle materials.

90-degree peeling of elastic thin films from elastic soft substrates: Theoretical solutions and experimental verification

Peeling of thin films has been widely used in adhesion measurement, film transfer and bio-inspired design. Most previous studies focused on the peeling of thin films from rigid substrates, but soft substrates are common in practical applications. Herein, we propose a two-dimensional model based on the bilinear cohesive law to characterize the 90-degree peeling of elastic thin films from elastic soft substrates, and obtain theoretical solutions expressed in terms of the Chebyshev series. The theoretical solutions match well with the finite element method results, including the load-displacement curves and the bulging deformation of soft substrates. We find that with decreasing substrate modulus, the maximum peeling force (Pmax) decreases but the steady-state peeling force remains unchanged. With the present solutions, the interfacial strength and fracture energy can be extracted simultaneously from the 90-degree peeling experiments of thin film/soft substrate systems, and then the experimentally measured Pmax for different film thicknesses can be well predicted. Furthermore, we obtain a new power scaling law of Pmax, where the scaling exponent depends on substrate elasticity. These results can help us measure the interfacial properties of thin film/soft substrate systems via peel tests, and regulate their peeling behaviors by interface design.

Improving mode I fracture response of glass fiber-epoxy composite adhesive joints using sanding treatment and MWCNTs inclusion

This study investigates the enhancement of fracture behavior in laminated composites using a combination of surface treatments and nanomaterial incorporation. The adherent material used is glass fiber-reinforced polymer (GFRP), and the surface treatment was applied to the GFRP surface. The adhesive used in the joints is an epoxy-based adhesive, which was reinforced with Multi-Walled Carbon Nanotubes (MWCNTs) to improve its mechanical properties. Initially, mechanical sanding techniques were explored, revealing significant improvements of 55% in fracture energy and 38% in load-bearing capacity with optimal 240-grit sanding. Subsequently, MWCNTs were introduced at varying weight percentages (0.1%, 0.3%, and 0.5%) into the adhesive. This combined approach aimed to synergistically examine the impact of surface pretreatments and nanoparticle integration on the fracture behavior of the GFRP joints. Results highlight the potential for fracture resistance and enhanced load-carrying capabilities in reinforced specimens, with maximum load and fracture energy improvements of up to 92% and 50%, respectively, compared to unreinforced specimens and those subjected solely to sanding treatments. Moreover, the shift in failure mode from adhesive failure to cohesive substrate failure was observed and is attributed to the enhanced bonding area created by the 240-grit sanding and the increased fracture resistance from the incorporation of 0.3 wt% MWCNTs. The inclusion of MWCNTs provides additional toughening mechanisms such as crack bridging and pull-out, which significantly increase the fracture energy. This comprehensive investigation underscores the intricate interplay between mechanical preparation and nanomaterial incorporation in significantly improving the performance of laminated composites.

Analysis of the impact of exfoliated graphene oxide on the mechanical performance and in-plane fracture resistance of epoxy-based nanocomposite

Graphene oxide (GO) is a versatile material, derivative of graphene, and it has gained significant attention in the field of polymer nanocomposites due to its exceptional properties and unique structural features. These characteristics make it a successful secondary reinforcing agent for enhancing mechanical performance, fatigue resistance, and other various physical and thermal properties of nanocomposites. However, the low fracture toughness of nanocomposite has restricted the overall structural applications. In this research work, the epoxy nanocomposites were prepared with different weights of exfoliated graphene oxide (E-GO) nanoparticles by using a dual mixing probe ultra-sonication for homogeneous distribution and proper dispersion. A substantial enhancement in mechanical performance and fracture toughness of nanocomposite was observed. The good dispersion and interfacial adhesion among E-GO and epoxy matrix were investigated through the critical analysis of the fracture surfaces of the nanocomposite. The enhancement of mechanical performance of E-GO nanofiller-reinforced epoxy composite was observed superior at 1 wt. % of GO60 particle concentration. The maximum increment of mechanical properties such as tensile strength, flexural strength, and work of fracture were 40.763, 39.23, and 12.106 % respectively, whereas maximum tensile modulus and flexural modulus were 19.3724, 27.63 % at 1.5 wt. % of GO60 as compared to the neat composite. However, in the case of fracture toughness and energy, maximum improvement of 1.256 Mpa.m1/2 and 0.472 KJ/m2 respectively was observed at 1 wt. % of GO60. Such enhancement was primarily due to the bolstering of in-plane crack propagation resistance in the nanocomposites. Additionally, the highlight point in thermomechanical properties of the polymer nanocomposite is optimal storage modulus and damping factor such as 4066.75 MPa, and 0.46 respectively, for the GO60 at 1 wt.% nanocomposite. Moreover, GO60 at 1 wt.% demonstrates greater thermal stability, withstanding up to 50% material degradation at 413.65°C.

Sustainably sourced and DOPO-derived triols as hardeners for epoxy thermosets: A promising solution to synchronously improve flame retardancy, mechanical strength and toughness

The flammability and brittleness of the diglycidyl ether of bisphenol A (DGEBA)-type epoxy resin greatly restrict its application in electronic appliances and structural materials. However, the synchronous improvement in the flame retardancy, mechanical strength, and toughness of DGEBA-type epoxy thermosets remains a significant challenge. To overcome this challenge, in this study, a sustainably sourced and DOPO-derived triol (DOPO-GL-DCDA) was synthesized from a bio-based material, glycerol 1, 3-diglycerolate diacrylate (GL-DCDA). A series of epoxy thermosets with different phosphorus contents were prepared by introducing DOPO-GL-DCDA into the DGEBA matrix. EP/DOPO-GL-DCDA demonstrated excellent flame retardancy and mechanical properties without compromising transparency. Specifically, EP/DOPO-GL-DCDA-1.0P exhibited a limit oxygen index (LOI) of 34.0 % and a UL-94 V-0 rating, and the total heat release (THR) and peak heat release rate (PHRR) of EP/DOPO-GL-DCDA-1.5P decreased by 36.6 % and 46.8 %, respectively, compared with those of pure EP. In addition, compared with those of pure EP, the tensile strength, impact strength, fracture toughness (KIC), and fracture energy (GIC) of EP/DOPO-GL-DCDA-1.5P increased by 62.8 %, 69.9 %, 62.9 %, and 99.0 %, respectively. Moreover, EP/DOPO-GL-DCDA had better dielectric properties than pure EP. Under the premise of maintaining transparency, this work provided a promising solution to synchronously improve the flame retardancy, mechanical strength, and toughness of epoxy thermosets.

Universal phase-field mixture representation of thermodynamics and shock-wave mechanics in porous soft biologic continua.

A continuum mixture theory is formulated for large deformations, thermal effects, phase interactions, and degradation of soft biologic tissues suitable at high pressures and low to very high strain rates. Tissues consist of one or more solid and fluid phases and can demonstrate nonlinear anisotropic elastic, viscoelastic, thermoelastic, and poroelastic physics. Under extreme deformations or shock loading, tissues may fracture, tear, or rupture. Existing models do not account for all physics simultaneously, and most poromechanics and soft-tissue models assume incompressibility of some or all constituents, generally inappropriate for modeling shock waves or extreme compressions. Motivated by these prior limitations, a thermodynamically consistent formulation that combines a continuum theory of mixtures, compressible nonlinear anisotropic thermoelasticity, viscoelasticity, and phase-field mechanics of fracture is constructed to resolve the pertinent physics. A metric tensor of generalized Finsler space supplies geometric insight on effects of rearrangements of microstructure, for example degradation, growth, and remodeling. Shocks are modeled as singular surfaces. Hugoniot states and shock decay are analyzed: Solutions account for concurrent viscoelasticity, fracture, and interphase momentum and energy exchange not all contained in previous analyses. Suitability of the framework for representing blood, skeletal muscle, and liver is demonstrated by agreement with experimental data and observations across a range of loading rates and pressures. Insight into previously unresolved physics is obtained, for example importance of rate sensitivity of damage and quantification of effects of dissipation from viscoelasticity and phase interactions on shock decay.