Search

We develop a fractional return-mapping framework for power-law visco-elasto-plasticity. In our approach, the fractional viscoelasticity is accounted through canonical combinations of Scott-Blair elements to construct a series of well-known fractional linear viscoelastic models, such as Kelvin-Voigt, Maxwell, Kelvin-Zener and Poynting-Thomson. We also consider a fractional quasi-linear version of Fung's model to account for stress/strain nonlinearity. The fractional viscoelastic models are combined with a fractional visco-plastic device, coupled with fractional viscoelastic models involving serial combinations of Scott-Blair elements. We then develop a general return-mapping procedure, which is fully implicit for linear viscoelastic models, and semi-implicit for the quasi-linear case. We find that, in the correction phase, the discrete stress projection and plastic slip have the same form for all the considered models, although with different property and time-step dependent projection terms. A series of numerical experiments is carried out with analytical and reference solutions to demonstrate the convergence and computational cost of the proposed framework, which is shown to be at least first-order accurate for general loading conditions. Our numerical results demonstrate that the developed framework is more flexible, preserves the numerical accuracy of existing approaches while being more computationally tractable in the visco-plastic range due to a reduction of $50\%$ in CPU time. Our formulation is especially suited for emerging applications of fractional calculus in bio-tissues that present the hallmark of multiple viscoelastic power-laws coupled with visco-plasticity.

We introduce a data-driven fractional modeling framework aimed at complex materials, and particularly bio-tissues. From multi-step relaxation experiments of distinct anatomical locations of porcine urinary bladder, we identify an anomalous relaxation character, with two power-law-like behaviors for short/long long times, and nonlinearity for strains greater than 25%. The first component of our framework is an existence study, to determine admissible fractional viscoelastic models that qualitatively describe the linear relaxation behavior. After the linear viscoelastic model is selected, the second stage adds large-strain effects to the framework through a fractional quasi-linear viscoelastic approach, given by a multiplicative kernel decomposition of the selected relaxation function and a nonlinear elastic response for the bio-tissue of interest. From single-relaxation data of the urinary bladder, a fractional Maxwell model captures both short/long-term behaviors with two fractional orders, being the most suitable model for small strains at the first stage. For the second stage, multi-step relaxation data under large strains were employed to calibrate a four-parameter fractional quasi-linear viscoelastic model, that combines a Scott-Blair relaxation function and an exponential instantaneous stress response, to describe the elastin/collagen phases of bladder rheology. Our obtained results demonstrate that the employed fractional quasi-linear model, with a single fractional order in the range $\alpha = 0.25 - 0.30$, is suitable for the porcine urinary bladder, producing errors below 2% without need for recalibration over subsequent applied strains. We conclude that fractional models are attractive tools to capture the bladder tissue behavior under small-to-large strains and multiple time scales, therefore being potential alternatives to describe multiple stages of bladder functionality.

From their birth in the manufacturing process, materials inherently contain defects that affect the mechanical behavior across multiple length and time-scales, including vacancies, dislocations, voids and cracks. Understanding, modeling, and real-time simulation of the underlying stochastic micro-structure defect evolution is therefore vital towards multi-scale coupling and propagating numerous sources of uncertainty from atomistic to eventually aging continuum mechanics. We develop a graph-based surrogate model of dislocation glide for computation of dislocation mobility. We model an edge dislocation as a random walker, jumping between neighboring nodes of a graph following a Poisson stochastic process. The network representation functions as a coarse-graining of a molecular dynamics simulation that provides dislocation trajectories for an empirical computation of jump rates. With this construction, we recover the original atomistic mobility estimates, with remarkable computational speed-up and accuracy. Furthermore, the underlying stochastic process provides the statistics of dislocation mobility associated to the original molecular dynamics simulation, allowing an efficient propagation of material parameters and uncertainties across the scales.

Fractional models and their parameters are sensitive to changes in the intrinsic micro-structures of anomalous materials. We investigate how such physics-informed models propagate the evolving anomalous rheology to the nonlinear dynamics of mechanical systems. In particular, we analyze the vibration of a fractional, geometrically nonlinear viscoelastic cantilever beam, under base excitation and free vibration, where the viscoelastic response is general through a distributed-order fractional model. We employ Hamilton's principle to obtain the corresponding equation of motion with the choice of specific material distribution functions that recover a fractional Kelvin-Voigt viscoelastic model of order $\alpha$. Through spectral decomposition in space, the resulting time-fractional partial differential equation reduces to a nonlinear time-fractional ordinary differential equation, in which the linear counterpart is numerically integrated by employing a direct L1-difference scheme. We further develop a semi-analytical scheme to solve the nonlinear system through a method of multiple scales, which yields a cubic algebraic equation in terms of the frequency. Our numerical results suggest a set of $\alpha$-dependent anomalous dynamic qualities, such as far-from-equilibrium power-law amplitude decay rates, super-sensitivity of amplitude response at free vibration, and bifurcation in steady-state amplitude at primary resonance., Comment: arXiv admin note: substantial text overlap with arXiv:1909.02142

We develop a thermodynamically consistent, fractional visco-elasto-plastic model coupled with damage for anomalous materials. The model utilizes Scott-Blair rheological elements for both visco-elastic/plastic parts. The constitutive equations are obtained through Helmholtz free-energy potentials for Scott-Blair elements, together with a memory-dependent fractional yield function and dissipation inequalities. A memory-dependent Lemaitre-type damage is introduced through fractional damage energy release rates. For time-fractional integration of the resulting nonlinear system of equations, we develop a first-order semi-implicit fractional return-mapping algorithm. We also develop a finite-difference discretization for the fractional damage energy release rate, which results into Hankel-type matrix-vector operations for each time-step, allowing us to reduce the computational complexity from $\mathcal{O}(N^3)$ to $\mathcal{O}(N^2)$ through the use of Fast Fourier Transforms. Our numerical results demonstrate that the fractional orders for visco-elasto-plasticity play a crucial role in damage evolution, due to the competition between the anomalous plastic slip and bulk damage energy release rates.

Efficient long-time integration of nonlinear fractional differential equations is significantly challenging due to the integro-differential nature of the fractional operators. In addition, the inherent non-smoothness introduced by the inverse power-law kernels deteriorates the accuracy and efficiency of many existing numerical methods. We develop two efficient first- and second-order implicit-explicit (IMEX) methods for accurate time-integration of stiff/nonlinear fractional differential equations with fractional order $\alpha \in (0,1]$ and prove their convergence and linear stability properties. The developed methods are based on a linear multi-step fractional Adams-Moulton method (FAMM), followed by the extrapolation of the nonlinear force terms. In order to handle the singularities nearby the initial time, we employ Lubich-like corrections to the resulting fractional operators. The obtained linear stability regions of the developed IMEX methods are larger than existing IMEX methods in the literature. Furthermore, the size of the stability regions increase with the decrease of fractional order values, which is suitable for stiff problems. We also rewrite the resulting IMEX methods in the language of nonlinear Toeplitz systems, where we employ a fast inversion scheme to achieve a computational complexity of $\mathcal{O}(N \log N)$, where $N$ denotes the number of time-steps. Our computational results demonstrate that the developed schemes can achieve global first- and second-order accuracy for highly-oscillatory stiff/nonlinear problems with singularities.

We investigate the nonlinear vibration of a fractional viscoelastic cantilever beam, subject to base excitation, where the viscoelasticity takes the general form of a distributed-order fractional model, and the beam curvature introduces geometric nonlinearity into the governing equation. We utilize the extended Hamilton principle to derive the governing equation of motion for specific material distribution functions that lead to fractional Kelvin-Voigt viscoelastic model. By spectral decomposition in space, the resulting governing fractional PDE reduces to nonlinear time-fractional ODEs. We use direct numerical integration in the decoupled system, in which we observe the anomalous power-law decay rate of amplitude in the linearized model. We further develop a semi-analytical scheme to solve the nonlinear equations, using method of multiple scales as a perturbation technique. We replace the expensive numerical time integration with a cubic algebraic equation to solve for frequency response of the system. We observe the super sensitivity of response amplitude to the fractional model parameters at free vibration, and bifurcation in steady-state amplitude at primary resonance.

Solutions to fractional models inherently exhibit non-smooth behavior, which significantly deteriorates the accuracy and therefore efficiency of existing numerical methods. We develop a two-stage data-infused computational framework for accurate time-integration of single- and multi-term fractional differential equations. In the first stage, we formulate a self-singularity-capturing scheme, given available/observable data for diminutive time. In this approach, the fractional differential equation provides the necessary knowledge/insight on how the hidden singularity can bridge between the initial and the subsequent short-time solution data. We develop a new self-singularity-capturing finite-difference algorithm for automatic determination of the underlying power-law singularities nearby the initial data, employing gradient descent optimization. In the second stage, we can utilize the multi-singular behavior of solution in a variety of numerical methods, without resorting to making any ad-hoc/uneducated guesses for the solution singularities. Particularly, we employed an implicit finite-difference method, where the captured singularities, in the first stage, are taken into account through some Lubich-like correction terms, leading to an accuracy of order $\mathcal{O}(\Delta t^{3-\alpha})$. Our computational results demonstrate that the developed framework can either fully capture or successfully control the solution error in the time-integration of fractional differential equations, especially in the presence of strong multi-singularities.