skip to main content


Title: Finite deformation and fractional order viscoelasticity of an auxetic foam

Auxetic foams exhibit novel mechanical properties due to their unique microstructure for improved energy-absorption and cavity expansion applications that have fascinated the scientific community since their inception. Given the advancements in material processing and performance of polymer open cell auxetic foams, there is a strong desire to fully understand the nonlinear rate-dependent deformation of these materials. The influence of nonlinear compressibility is introduced here along with relaxation effects to improve model predictions for different stretch rates and finite deformation regimes. The viscoelastic behavior of the material is analyzed by comparing fractional order and integer order calculus models. All results are statistically validated using maximum entropy methods to obtain Bayesian posterior densities for the hyperelastic, auxetic, and viscoelastic parameters. It is shown that fractional order viscoelasticity provides [Formula: see text]–[Formula: see text] improvement in prediction over integer order viscoelastic models when the model is calibrated at higher stretch rates where viscoelasticity is more significant.

 
more » « less
NSF-PAR ID:
10368902
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
SAGE Publications
Date Published:
Journal Name:
Journal of Intelligent Material Systems and Structures
Volume:
33
Issue:
14
ISSN:
1045-389X
Page Range / eLocation ID:
p. 1846-1861
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Fractional models and their parameters are sensitive to intrinsic microstructural changes in anomalous materials. We investigate how such physics-informed models propagate the evolving anomalous rheology to the nonlinear dynamics of mechanical systems. In particular, we study the vibration of a fractional, geometrically nonlinear viscoelastic cantilever beam, under base excitation and free vibration, where the viscoelasticity is described by a distributed-order fractional model. We employ Hamilton's principle to obtain the equation of motion with the choice of specific material distribution functions that recover a fractional Kelvin–Voigt viscoelastic model of order α. Through spectral decomposition in space, the resulting time-fractional partial differential equation reduces to a nonlinear time-fractional ordinary differential equation, where the linear counterpart is numerically integrated through a direct L1-difference scheme. We further develop a semi-analytical scheme to solve the nonlinear system through a method of multiple scales, yielding a cubic algebraic equation in terms of the frequency. Our numerical results suggest a set of α-dependent anomalous dynamic qualities, such as far-from-equilibrium power-law decay rates, amplitude super-sensitivity at free vibration, and bifurcation in steady-state amplitude at primary resonance. 
    more » « less
  2. Abstract Entropy dynamics is a Bayesian inference methodology that can be used to quantify time-dependent posterior probability densities that guide the development of complex material models using information theory. Here, we expand its application to non-Gaussian processes to evaluate how fractal structure can influence fractional hyperelasticity and viscoelasticity in elastomers. We investigate how kinematic constraints on fractal polymer network deformation influences the form of hyperelastic constitutive behavior and viscoelasticity in soft materials such as dielectric elastomers, which have applications in the development of adaptive structures. The modeling framework is validated on two dielectric elastomers, VHB 4910 and 4949, over a broad range of stretch rates. It is shown that local fractal time derivatives are equally effective at predicting viscoelasticity in these materials in comparison to nonlocal fractional time derivatives under constant stretch rates. We describe the origin of this accuracy that has implications for simulating large-scale problems such as finite element analysis given the differences in computational efficiency of nonlocal fractional derivatives versus local fractal derivatives. 
    more » « less
  3. Toughness of soft materials such as elastomers and gels depends on their ability to dissipate energy and to reduce stress concentration at the crack tip. The primary energy dissipation mechanism is viscoelasticity. Most analyses and models of fracture are based on linear viscoelastic theory (LVT) where strains are assumed to be small and the relaxation mechanisms are independent of stress or strain history. A well-known paradox is that the size of the dissipative zone predicted by LVT is unrealistically small. Here we use a physically based nonlinear viscoelastic model to illustrate why the linear theory breaks down. Using this nonlinear model and analogs of crack problems, we give a plausible resolution to this paradox. In our model, viscoelasticity arises from the breaking and healing of physical cross-links in the polymer network. When the deformation is small, the kinetics of bond breaking and healing are independent of the strain/stress history and the model reduces to the standard linear theory. For large deformations, localized bond breaking damages the material near the crack tip, reducing stress concentration and dissipating energy at the same time. The damage zone size is a new length scale which depends on the strain required to accelerate bond breaking kinetics. These effects are illustrated by considering two cases with stress concentrations: the evolution of spherical damage in a viscoelastic body subjected to internal pressure, and a zero degree peel test. 
    more » « less
  4. Unlike traditional structural materials, soft solids can often sustain very large deformation before failure, and many exhibit nonlinear viscoelastic behavior. Modeling nonlinear viscoelasticity is a challenging problem for a number of reasons. In particular, a large number of material parameters are needed to capture material response and validation of models can be hindered by limited amounts of experimental data available. We have developed a Gaussian Process (GP) approach to determine the material parameters of a constitutive model describing the mechanical behavior of a soft, viscoelastic PVA hydrogel. A large number of stress histories generated by the constitutive model constitute the training sets. The low-rank representations of stress histories by Singular Value Decomposition (SVD) are taken to be random variables which can be modeled via Gaussian Processes with respect to the material parameters of the constitutive model. We obtain optimal material parameters by minimizing an objective function over the input set. We find that there are many good sets of parameters. Further the process reveals relationships between the model parameters. Results so far show that GP has great potential in fitting constitutive models. 
    more » « less
  5. We have formulated the Fréchet kernel computation using the adjoint-state method based on a fractional viscoacoustic wave equation. We first numerically prove that the 1/2- and the 3/2-order fractional Laplacian operators are self-adjoint. Using this property, we find that the adjoint wave propagator preserves the dispersion and compensates the amplitude, whereas the time-reversed adjoint wave propagator behaves identically to the forward propagator with the same dispersion and dissipation characters. Without introducing rheological mechanisms, this formulation adopts an explicit [Formula: see text] parameterization, which avoids the implicit [Formula: see text] in the conventional viscoacoustic/viscoelastic full-waveform inversion ([Formula: see text]-FWI). In addition, because of the decoupling of operators in the wave equation, the viscoacoustic Fréchet kernel is separated into three distinct contributions with clear physical meanings: lossless propagation, dispersion, and dissipation. We find that the lossless propagation kernel dominates the velocity kernel, whereas the dissipation kernel dominates the attenuation kernel over the dispersion kernel. After validating the Fréchet kernels using the finite-difference method, we conduct a numerical example to demonstrate the capability of the kernels to characterize the velocity and attenuation anomalies. The kernels of different misfit measurements are presented to investigate their different sensitivities. Our results suggest that, rather than the traveltime, the amplitude and the waveform kernels are more suitable to capture attenuation anomalies. These kernels lay the foundation for the multiparameter inversion with the fractional formulation, and the decoupled nature of them promotes our understanding of the significance of different physical processes in [Formula: see text]-FWI. 
    more » « less