Showing papers on "Viscoplasticity published in 2017"

Numerical simulations of complex yield-stress fluid flows

Abstract: Viscoplasticity is characterized by a yield stress, below which the materials will not deform and above which they will deform and flow according to different constitutive relations. Viscoplastic models include the Bingham plastic, the Herschel-Bulkley model and the Casson model. All of these ideal models are discontinuous. Analytical solutions exist for such models in simple flows. For general flow fields, it is necessary to develop numerical techniques to track down yielded/unyielded regions. This can be avoided by introducing into the models a regularization parameter, which facilitates the solution process and produces virtually the same results as the ideal models by the right choice of its value. This work reviews several benchmark problems of viscoplastic flows, such as entry and exit flows from dies, flows around a sphere and a bubble and squeeze flows. Examples are also given for typical processing flows of viscoplastic materials, where the extent and shape of the yielded/unyielded regions are clearly shown. The above-mentioned viscoplastic models leave undetermined the stress and elastic deformation in the solid region. Moreover, deviations have been reported between predictions with these models and experiments for flows around particles using Carbopol, one of the very often used and heretofore widely accepted as a simple “viscoplastic” fluid. These have been partially remedied in very recent studies using the elastoviscoplastic models proposed by Saramito.

A modern look on yield stress fluids

Abstract: A concept of viscoplasticity advanced exactly one century ago by Bingham appears very fruitful because there are many natural and artificial materials that demonstrate viscoplastic behavior, i.e., they are able to pass from a solid to a liquid state under the influence of applied stress. However, although this transition was originally considered as a jump-like phenomenon occurring at a certain stress—the yield stress—numerous subsequent studies have shown that the real situation is more complicated. A long-term discussion about the possibility of flow at low stresses less than the yield stress came to today’s conclusion denying this possibility as being opposite to the existence of the maximal Newtonian viscosity in viscoelastic polymeric fluids. So, there is a contradiction between the central dogma of rheology which says that “everything flows” and the alleged impossibility for flow at a solid-like state of viscoplastic fluids. Then, the concept of the fragile destruction of an inner structure responsible for a solid-like state at the definite (yield) stress was replaced by an understanding of the yielding as a transition extending over some stress range and occurring in time. So, instead of the yield stress, yielding is characterized by the dependence of durability (or time-to-break) on the applied stress. In this review, experimental facts and the new understanding of yielding as a kinetic process are discussed. Besides, some other alternative methods for measuring the yield stress are considered.

The yielding and the linear-to-nonlinear viscoelastic transition of an elastoviscoplastic material

Abstract: Elastoviscoplastic materials present a transition from a gel-like to a liquidlike state induced by shearing: While the first is primarily elastic, the second is predominantly viscous. The point that characterizes this transition is usually known as the yield point, which is associated to critical quantities such as yield stress and/or yield strain. Another characteristic of elastoviscoplastic materials is the transition from linear to nonlinear viscoelasticity. In the current work, a commercial hair gel, which is an elastoviscoplastic material, was tested in two rotational rheometers in order to evaluate these two transition points. Stress oscillatory amplitude sweeps at different frequencies were performed and a Fourier-Transform analysis was applied to the results in order to determine the linear viscoelastic limit. The linear viscoelastic limit stresses and strains at different frequencies were then compared to quantities that are usually associated to the yield point: The extrapolated zero-shear-rate ...

Time-Dependent Mechanical Response of APbX3 (A = Cs, CH3NH3; X = I, Br) Single Crystals

Abstract: The ease of processing hybrid organic-inorganic perovskite (HOIPs) films, belonging to a material class with composition ABX3 , from solution and at mild temperatures promises their use in deformable technologies, including flexible photovoltaic devices, sensors, and displays. To successfully apply these materials in deformable devices, knowledge of their mechanical response to dynamic strain is necessary. The authors elucidate the time- and rate-dependent mechanical properties of HOIPs and an inorganic perovskite (IP) single crystal by measuring nanoindentation creep and stress relaxation. The observation of pop-in events and slip bands on the surface of the indented crystals demonstrate dislocation-mediated plastic deformation. The magnitudes of creep and relaxation of both HOIPs and IPs are similar, negating prior hypothesis that the presence of organic A-site cations alters the mechanical response of these materials. Moreover, these samples exhibit a pronounced increase in creep, and stress relaxation as a function of indentation rate whose magnitudes reflect differences in the rates of nucleation and propagation of dislocations within the crystal structures of HOIPs and IP. This contribution provides understanding that is critical for designing perovskite devices capable of withstanding mechanical deformations.

A finite difference method for off-fault plasticity throughout the earthquake cycle

Abstract: We have developed an efficient computational framework for simulating multiple earthquake cycles with off-fault plasticity. The method is developed for the classical antiplane problem of a vertical strike-slip fault governed by rate-and-state friction, with inertial effects captured through the radiation-damping approximation. Both rate-independent plasticity and viscoplasticity are considered, where stresses are constrained by a Drucker–Prager yield condition. The off-fault volume is discretized using finite differences and tectonic loading is imposed by displacing the remote side boundaries at a constant rate. Time-stepping combines an adaptive Runge–Kutta method with an incremental solution process which makes use of an elastoplastic tangent stiffness tensor and the return-mapping algorithm. Solutions are verified by convergence tests and comparison to a finite element solution. We quantify how viscosity, isotropic hardening, and cohesion affect the magnitude and off-fault extent of plastic strain that develops over many ruptures. If hardening is included, plastic strain saturates after the first event and the response during subsequent ruptures is effectively elastic. For viscoplasticity without hardening, however, successive ruptures continue to generate additional plastic strain. In all cases, coseismic slip in the shallow sub-surface is diminished compared to slip accumulated at depth during interseismic loading. The evolution of this slip deficit with each subsequent event, however, is dictated by the plasticity model. Integration of the off-fault plastic strain from the viscoplastic model reveals that a significant amount of tectonic offset is accommodated by inelastic deformation ( ∼ 0.1 m per rupture, or ∼ 10% of the tectonic deformation budget).

Arbitrary order 2D virtual elements for polygonal meshes: Part II, inelastic problem

Abstract: The present paper is the second part of a twofold work, whose first part is reported in [3], concerning a newly developed Virtual Element Method (VEM) for 2D continuum problems. The first part of the work proposed a study for linear elastic problem. The aim of this part is to explore the features of the VEM formulation when material nonlinearity is considered, showing that the accuracy and easiness of implementation discovered in the analysis inherent to the first part of the work are still retained. Three different nonlinear constitutive laws are considered in the VEM formulation. In particular, the generalized viscoplastic model, the classical Mises plasticity with isotropic/kinematic hardening and a shape memory alloy (SMA) constitutive law are implemented. The versatility with respect to all the considered nonlinear material constitutive laws is demonstrated through several numerical examples, also remarking that the proposed 2D VEM formulation can be straightforwardly implemented as in a standard nonlinear structural finite element method (FEM) framework.

Yield limit analysis of particle motion in a yield-stress fluid

Abstract: A theoretical and numerical study of yield-stress fluid creeping flow about a particle is presented. Yield-stress fluids can hold rigid particles statically buoyant if the yield stress is large enough. In addressing sedimentation of rigid particles in viscoplastic fluids, we should know this critical ‘yield number’ beyond which there is no motion. As we get close to this limit, the role of viscosity becomes negligible in comparison to the plastic contribution in the leading order, since we are approaching the zero-shear-rate limit. Admissible stress fields in this limit can be found by using the characteristics of the governing equations of perfect plasticity (i.e. the sliplines). This approach yields a lower bound of the critical plastic drag force or equivalently the critical yield number. Admissible velocity fields also can be postulated to calculate the upper bound. This analysis methodology is examined for three families of particle shapes (ellipse, rectangle and diamond) over a wide range of aspect ratios. Numerical experiments of either resistance or mobility problems in a viscoplastic fluid validate the predictions of slipline theory and reveal interesting aspects of the flow in the yield limit (e.g. viscoplastic boundary layers). We also investigate in detail the cases of high and low aspect ratio of the particles.

Multiscale testing-analysis of asphaltic materials considering viscoelastic and viscoplastic deformation

Abstract: Fine aggregate matrix (FAM) is a phase consisting of asphalt binder, air voids, fine aggregates and fillers It acts as a primary phase in evaluating the damage and deformation of entire asphalt concrete mixtures The simplicity, repeatability and efficiency of the FAM testing make it a very attractive specification-type approach for evaluating the performance characteristics of the entire asphalt concrete mixtures This study explores a linkage in the deformation characteristics between the two length scales: asphalt concrete mixture scale and its corresponding FAM scale To that end, a simple creep-recovery test was conducted for both mixtures (ie asphalt concrete mixture and its corresponding FAM phase) at various stress levels Test results were compared and analysed using Schapery’s single-integral viscoelastic theory and Perzyna-type viscoplasticity with a generalised Drucker–Prager yield surface In particular, stress-dependent nonlinear viscoelastic and viscoplastic behaviours were chara

A dislocation-based crystal plasticity framework for dynamic ductile failure of single crystals

Abstract: A framework for dislocation-based viscoplasticity and dynamic ductile failure has been developed to model high strain rate deformation and damage in single crystals. The rate-dependence of the crystal plasticity formulation is based on the physics of relativistic dislocation kinetics suited for extremely high strain rates. The damage evolution is based on the dynamics of void growth, which are governed by both micro-inertia as well as dislocation kinetics and dislocation substructure evolution. An averaging scheme is proposed in order to approximate the evolution of the dislocation substructure in both the macroscale as well as its spatial distribution at the microscale. Additionally, a concept of a single equivalent dislocation density that effectively captures the collective influence of dislocation density on all active slip systems is proposed here. Together, these concepts and approximations enable the use of semi-analytic solutions for void growth dynamics developed in (Wilkerson and Ramesh, 2014), which greatly reduce the computational overhead that would otherwise be required. The resulting homogenized framework has been implemented into a commercially available finite element package, and a validation study against a suite of direct numerical simulations was carried out.

Comparison of two creep degradation modeling approaches for soft structured soils

Abstract: The creep degradation is a common phenomenon for soft structured clays. In this paper, the creep degradation behavior for soft structured clays is first studied by combining intrinsic creep behavior and the structure indicator. A creep-implicit model and a creep-explicit model corresponding to a stress-based and a creep-based structure indicators are developed, respectively, under one-dimensional condition. Parameters determination for both models is straightforward from oedometer tests. Coupled with consolidation theory, both models are used to simulate oedometer tests with different structural levels and load durations on three clays. The predictive ability of the two models on creep behavior, creep degradation behavior and evolution of structure indicator is analyzed. The relationship between the two structure indicators is discussed based on experimental results. The comparison between experimental and numerical results demonstrates that both models can accurately describe the creep degradation behavior of soft structured clay under one-dimensional loading.

Twin-induced hardening in extruded Mg alloy AM30

Abstract: The mechanical anisotropy of extruded Mg alloy AM30 was recently investigated, both experimentally and numerically (Oppedal et al., 2013 [1]). The authors highlighted the need to include cases of intermediate levels of twinning as a critical validation, in addition to the limiting cases where the contribution of twinning is dominant or negligible. However, further scrutiny of the results revealed that the experimental data employed in that study were still inadequate to constrain all necessary parameters in the crystal plasticity models. In particular, the parameters describing the prismatic slip mechanism were under-constrained. In the present study, additional experimental data were obtained to provide the necessary constraint. Based upon these experimental results, clearer conclusions can be drawn about the requirements of a crystal plasticity model, which must accurately account for strong tension-compression asymmetry as well as anisotropy in both strength and strain hardening behavior. In addition to the previously employed Viscoplastic Self-Consistent (VPSC) model, with the Predominant Twin Reorientation (PTR) scheme, the Elastic Visco-Plastic Self-Consistent (EVPSC) model, with the recently developed Twinning and De-Twinning (TDT) description, is applied to simulate the uniaxial response along arbitrary directions. It is demonstrated that accounting for the initial texture and calibrating the EVPSC-TDT model using uniaxial tension and compression along the extrusion direction permits prediction of the strength anisotropy and strain hardening behavior along arbitrary straining directions.