Modeling the Viscoelastic Behavior of Quartz and Clay Minerals in Shale by Nanoindentation Creep Tests

The viscoelastic behavior of minerals in shales is important in predicting the macroscale creep behavior of heterogeneous bulk shale. In this study, in situ indentation measurements of two major constitutive minerals (i.e., quartz and clay) in Longmaxi Formation shale from the Sichuan Basin, South China, were conducted using a nanoindentation technique and high-resolution optical microscope. Firstly, quartz and clay minerals were identified under an optical microscope based on their morphology, surface features, reflection characteristics, particle shapes, and indentation responses. Three viscoelastic models (i.e., three-element Voigt, Burger’s, and two-dashpot Kelvin models) were then used to fit the creep data for both minerals. Finally, the effects of peak load on the viscoelastic behavior of quartz and clay minerals were investigated. Our results show that the sizes of the residual imprints on clay minerals were larger than that of quartz for a specific peak load. Moreover, the initial creep rates and depths in clay minerals were higher than those in quartz. However, the creep rates of quartz and clay minerals displayed similar trends, which were independent of peak load. In addition, all three viscoelastic models produced good fits to the experimental data. However, due to the poor fit in the initial holding stage of the three-element Voigt model and instability of the two-dashpot Kelvin model, Burger’s model is best in obtaining the regression parameters. The regression results indicate that the viscoelastic parameters obtained by these models are associated with peak load, and that a relatively small peak load is more reliable for the determination of viscoelastic parameters. Furthermore, the regression values for the viscoelastic parameters of clay minerals were lower than those of quartz and the bulk shale, suggesting the former facilitates the viscoelastic deformation of shale. Our study provides a better understanding of the nanoscale viscoelastic properties of shale, which can be used to predict the time-dependent deformation of shale.

An unified formulation of strong non-local elasticity with fractional order calculus

The research of a formulation to model non-local interactions in the mechanical behavior of matter is currently an open problem. In this context, a strong non-local formulation based on fractional calculus is provided in this paper. This formulation is derived from an analogy with long-memory viscoelastic models. Specifically, the same kind of power-law time-dependent kernel used in Boltzmann integral of viscoelastic stress-strain relation is used as kernel in the Fredholm non-local relation. This non-local formulation leads to stress-strain relation based on the space Riesz integral and derivative of fractional order. For unbounded domain, proposed model can be defined in stress- and in strain-driven formulation and in both cases the stress–strain relation represent a strong non-local model. Also, the proposed strain driven and stress driven formulations defined in terms of Riesz operators are proved to be fully consistent each another. Moreover, the proposed model posses a mechanical meaning and for unbounded non-local rod is described and discussed in detail.

A new method for obtaining model-free viscoelastic material properties from atomic force microscopy experiments using discrete integral transform techniques

Viscoelastic characterization of materials at the micro- and the nanoscale is commonly performed with the aid of force–distance relationships acquired using atomic force microscopy (AFM). The general strategy for existing methods is to fit the observed material behavior to specific viscoelastic models, such as generalized viscoelastic models or power-law rheology models, among others. Here we propose a new method to invert and obtain the viscoelastic properties of a material through the use of the Z-transform, without using a model. We present the rheological viscoelastic relations in their classical derivation and their z-domain correspondence. We illustrate the proposed technique on a model experiment involving a traditional ramp-shaped force–distance AFM curve, demonstrating good agreement between the viscoelastic characteristics extracted from the simulated experiment and the theoretical expectations. We also provide a path for calculating standard viscoelastic responses from the extracted material characteristics. The new technique based on the Z-transform is complementary to previous model-based viscoelastic analyses and can be advantageous with respect to Fourier techniques due to its generality. Additionally, it can handle the unbounded inputs traditionally used to acquire force–distance relationships in AFM, such as ramp functions, in which the cantilever position is displaced linearly with time for a finite period of time.

Characterising the linear viscoelastic behaviour of an injection moulding grade polypropylene polymer

The linear viscoelastic behaviour of an injection moulding grade polypropylene is studied using theoretical and computational methods. Polypropylene has a variety of engineering applications as a component. However, it commonly exhibits viscoelastic deformations. This paper analyses the creep and recovery responses of the BJ368MO polypropylene copolymer using the Burgers and generalised Maxwell models. Within the linear viscoelastic regime, an experimental creep strain at $20\ \text{MPa}$ 20 MPa is used to determine the rheological constants of the models. These constants (springs and dashpots) are determined using a nonlinear least-squares curve fitting of the experimental creep. Then they are used to predict the creep and recovery responses of the polymer at three different stresses, $10\ \text{MPa}$ 10 MPa , $12.5\ \text{MPa}$ 12.5 MPa and $15\ \text{MPa}$ 15 MPa . The experiments are made using tensile specimens designed according to the ASTM D638-14standard. The theoretical evaluations are made using the creep and recovery equations derived from their constitutive. Whereas COMSOL Multiphysics software is used during the finite element (FE) analyses. The results of the theoretical and FE calculations are verified using creep and recovery experiments. Based on the validation analyses, both viscoelastic models showed lower deviations from the experimental results when a computational approach is used. In addition, the viscoelastic models are compared by evaluating the residuals of the creep and recovery strain predictions. The theoretical analyses showed better predictions at $12.5\ \text{MPa}$ 12.5 MPa and $15\ \text{MPa}$ 15 MPa stresses when the generalised Maxwell model is used. However, the improvements are attributed to the recovery predictions. When FE is used, the Burgers model showed lower mean absolute percentage errors (MAPEs) in all creep and recovery predictions. The model has a minimum of 6.37% error at the $10\ \text{MPa}$ 10 MPa stress and a maximum of 8.23% error at the $15\ \text{MPa}$ 15 MPa . By comparison, the generalised Maxwell model showed a minimum of 9.24% error at $12.5\ \text{MPa}$ 12.5 MPa and a maximum of 12.8% error at $15\ \text{MPa}$ 15 MPa stresses. The novelty of this paper is on predicting the creep and recovery behaviour of the polymer using the FE and theoretical approaches in the linear viscoelastic regime. The findings suggest that the FE analyses using the Burgers viscoelastic material model provide better predictions, with all calculated errors falling below 10%.

Analytical Investigation of Viscoelastic Stagnation-Point Flows with Regard to Their Singularity

Singularities in the stress field of the stagnation-point flow of a viscoelastic fluid have been studied for various viscoelastic constitutive models. Analyzing the analytical solutions of these models is the most effective way to study this problem. In this paper, exact analytical solutions of two-dimensional steady wall-free stagnation-point flows for the generic Oldroyd 8-constant model are obtained for the stress field using different material parameter relations. For all solutions, compatibility with the conservation of momentum is considered in our analysis. The resulting solutions usually contain arbitrary functions, whose choice has a crucial effect on the stress distribution. The corresponding singularities are discussed in detail according to the choices of the arbitrary functions. The results can be used to analyze the stress distribution and singularity behavior of a wide spectrum of viscoelastic models derived from the Oldroyd 8-constant model. Many previous results obtained for simple viscoelastic models are reproduced as special cases. Some previous conclusions are amended and new conclusions are drawn. In particular, we find that all models have singularities near the stagnation point and most of them can be avoided by appropriately choosing the model parameters and free functions. In addition, the analytical solution for the stress tensor of a near-wall stagnation-point flow for the Oldroyd-B model is also obtained. Its compatibility with the momentum conservation is discussed and the parameters are identified, which allow for a non-singular solution.

A New Method for Obtaining Model-Free Viscoelastic Material Properties from Atomic Force Microscopy Experiments Using Discrete Integral Transform Techniques

Viscoelastic characterization of materials at the micro- and nanoscales is commonly performed with the aid of force-distance relationships acquired using atomic force microscopy (AFM). The general strategy for existing methods is to fit the observed material behavior to specific viscoelastic models, such as generalized viscoelastic models or power-law rheology models, among others. Here we propose a new method to invert and obtain the viscoelastic properties of a material through the use of the Z-transform, without using a model. We present the rheological viscoelastic relations in their classical derivation and their Z-domain correspondence. We illustrate the proposed technique on a model experiment involving a traditional ramp-shaped force-distance AFM curve, demonstrating good agreement between the viscoelastic characteristics extracted from the simulated experiment and the theoretical expectations. We also provide a path for calculating standard viscoelastic responses from the extracted material characteristics. The new technique based on the Z-transform is complementary to previous model-based viscoelastic analyses and can be advantageous with respect to Fourier techniques due to its generality. Additionally, it can handle the unbounded inputs traditionally used to acquire force-distance relationships in AFM, such as “ramp” functions, in which the cantilever position is displaced linearly with time for a finite period of time.

Fractional-order viscoelastic model of musculoskeletal tissues: correlation with fractals

Self-similar fractals are widely obtained from biomaterials within the human musculoskeletal system, and their viscoelastic behaviours can be described by fractional-order derivatives. However, existing viscoelastic models neglect the internal correlation between the fractal structure of biomaterials and their fractional-order temporal responses. We further expanded the fractal hyper-cell (FHC) viscoelasticity theory to investigate this spatio-temporal correlation. The FHC element was first compared with other material elements and spring–dashpot viscoelastic models, thereby highlighting its discrete and fractal nature. To demonstrate the utility of an FHC, tree-like, ladder-like and triangle-like FHCs were abstracted from human cartilage, tendons and muscle cross-sections, respectively. The duality and symmetry of the FHC element were further discussed, where operating the duality transformation generated new types of FHC elements, and the symmetry breaking of an FHC altered its fractional-order viscoelastic responses. Thus, the correlations between the staggering patterns of FHCs and their rheological power-law orders were established, and the viscoelastic behaviour of the multi-level FHC elements fitted well in stress relaxation experiments at both the macro- and nano-levels of the tendon hierarchy. The FHC element provides a theoretical basis for understanding the connections between structural degeneration of bio-tissues during ageing or disease and their functional changes.