A constitutive equation for filled rubber under cyclic loading

This paper describes experiments and the development of constitutive equations to predict the steady-state response of filled rubber under cyclic loading. An MTS servo-hydraulic machine was used to obtain the dynamic hysteresis curves for a filled rubber compound in uniaxial tension-compression. The material tests were performed with mean strains from −0.1 to 0.1, strain amplitudes ranging from 0.02 to 0.1, and strain rates between 0.01 and 10 s⁻¹. Temporary material set, the Payne effect and rate-dependence were observed from the experimental results. A hyper-viscoelastic constitutive model was developed to characterize the dynamic response of the rubber. A cornerstone of this constitutive modeling was to devise a scheme to evaluate material set and a finite strain, non-linear viscoelastic law from the test data. Predictions of the dynamic hysteresis curves using the proposed constitutive equation were found to be in good agreement with the uniaxial test results.     

A new constitutive equation for elastoviscoplastic fluid flows

From thermodynamic theory, a new three-dimensional model for elastoviscoplastic fluid flows is presented. It extends both the Bingham viscoplastic and the Oldroyd viscoelastic models. Fundamental flows are studied: simple shear flow, uniaxial elongation and large amplitude oscillatory shear. The complex moduli (G^′,G^″)$(G^{\prime }\text{,}{G}^{″})$ are found to be in qualitative agreement with experimental data for materials that present microscopic network structures and large scale rearrangements. Various fluids of practical interest, such as liquid foams, droplet emulsions or blood, present such elastoviscoplastic behavior: at low stress, the material behaves as a viscoelastic solid, whereas at stresses above a yield stress, the material behaves as a fluid.     

A plastic constitutive equation incorporating strain,strain-rate,and temperature

An empirical plasticity constitutive form describing the flow stress as a function of strain, strain-rate, and temperature has been developed, fit to data for three dual-phase (DP) steels, and compared with independent experiments outside of the fit domain. Dubbed the “H/V model” (for “Hollomon/Voce”), the function consists of three multiplicative functions describing (a) strain hardening, (b) strain-rate sensitivity, and (c) temperature sensitivity. Neither the multiplicative structure nor the choice of functions (b) or (c) is novel. The strain hardening function, (a), has two novel features: (1) it incorporates a linear combination coefficient, α, that allows representation of Hollomon (power law) behavior (α = 1), Voce (saturation) behavior (α = 0) or any intermediate case (0 < α < 1, and (2) it allows incorporation of the temperature sensitivity of strain hardening rate in a natural way by allowing α to vary with temperature (in the simplest case, linearly). This form therefore allows a natural transition from unbounded strain hardening at low temperatures toward saturation behavior at higher temperatures, consistent with many observations. Hollomon, Voce, H/V models and others selected as representative from the literature were fit for DP590, DP780, and DP980 steels by least-squares using a series of tensile tests up to the uniform strain conducted over a range of temperatures. Jump-rate tests were used to probe strain rate sensitivity. The selected laws were then used with coupled thermo-mechanical finite element (FE) modeling to predict behavior for tests outside the fit range: non-isothermal tensile tests beyond the uniform strain at room temperatures, isothermal tensile tests beyond the uniform strain at several temperatures and hydraulic bulge tests at room temperature. The agreement was best for the H/V model, which captured strain hardening at high strain accurately as well as the variation of strain hardening with temperature. The agreement of FE predictions up to the tensile failure strain illustrates the critical role of deformation-induced heating in high-strength/high ductility alloys, the importance of having a constitutive model that is accurate at large strains, and the implication that damage and void growth are unlikely to be determinant factors in the tensile failure of these alloys. The new constitutive model may have application for a wide range of alloys beyond DP steels, and it may be extended to larger strain rate and temperature ranges using alternate forms of strain rate sensitivity and thermal softening appearing in the literature.     

A cyclic plasticity model for single crystals

The Armstrong–Frederick type kinematic hardening rule was invoked to capture the Bauschinger effect of the cyclic plastic deformation of a single crystal. The yield criterion and flow rule were built on individual slip systems. Material memory was introduced to describe strain range dependent cyclic hardening. The experimental results of copper single crystals were used to evaluate the cyclic plasticity model. It was found that the model was able to accurately describe the cyclic plastic deformation and properly reflect the dislocation substructure evolution. The well-known three distinctive regimes in the cyclic stress–strain curve of the copper single crystals oriented for single slip can be reproduced by using the model. The model can predict the enhanced hardening for crystals oriented for multislip, showing the model's ability to describe anisotropic cyclic plasticity. For a given loading history, the model was able to capture not only the saturated stress–strain response but also the detailed transient stress–strain evolution. The model was used to predict the cyclic plasticity under a high–low loading sequence. Both the stress–strain responses and the microstructural evolution can be appropriately described through the slip system activation.     

The key role of dislocation dissociation in the plastic behaviour of single crystal nickel-based superalloy with low stacking fault energy: Three-dimensional discrete dislocation dynamics modelling

To model the deformation of single crystal nickel based superalloys (SCNBS) with low stacking fault energy (SFE), three-dimensional discrete dislocation dynamics (3D-DDD) is extended by incorporating dislocation dissociation mechanism. The present 3D-DDD simulations show that, consistent with the existing TEM observation, the leading partial can enter the matrix channel efficiently while the trailing partial can hardly glide into it when the dislocation dissociation is taken into account. To determine whether the dislocation dissociation can occur or not, a critical percolation stress (CPS) based criterion is suggested. According to this CPS criterion, for SCNBS there exists a critical matrix channel width. When the channel width is lower than this critical value, the dislocation tends to dissociate into an extended configuration and vice versa. To clarify the influence of dislocation dissociation on CPS, the classical Orowan formula is improved by incorporating the SFE. Moreover, the present 3D-DDD simulations also show that the yielding stress of SCNBSs with low SFE may be overestimated up to 30% if the dislocation dissociation is ignored. With dislocation dissociation being considered, the size effect due to the width of γ matrix channel and the length of γ′ precipitates on the stress–strain responses of SCNBS can be enhanced remarkably. In addition, due to the strong constraint effect by the two-phase microstructure in SCNBS, the configuration of formed junctions is quite different from that in single phase crystals such as Cu. The present results not only provide clear understanding of the two-phase microstructure levelled microplastic mechanisms in SCNBSs with low SFE, but also help to develop new continuum-levelled constitutive laws for SCNBSs.     

An empirical constitutive equation of integral type for viscoelastic liquids

A 5-constant constitutive equation is proposed. The analytical form for the relaxation modulus as a function of flow conditions was chosen based on experimental data for stress-relaxation in solid polymers. The resulting formulae for the material functions in simple and oscillatory shear flow fulfil the empirical Cox-Merz rule as well as other phenomenological relations formulated by Coleman and Markowitz. The theoretical results are compared with experimental data obtained by Han for various polymer melts. Good agreement between theory and experiment is found.     

Influence of single crystal orientation on homogeneous dislocation nucleation under uniaxial loading

Atomistic simulations are used to investigate how the stress required for homogeneous nucleation of partial dislocations in single crystal copper under uniaxial loading changes as a function of crystallographic orientation. Molecular dynamics is employed based on an embedded-atom method potential for Cu at 10 and 300 K. Results indicate that non-Schmid parameters are important for describing the calculated dislocation nucleation behavior for single crystal orientations under tension and compression. A continuum relationship is presented that incorporates Schmid and non-Schmid terms to correlate the nucleation stress over all tensile axis orientations within the stereographic triangle. Simulations investigating the temperature dependence of homogeneous dislocation nucleation yield activation volumes of ≈0.5- and activation energies of . For uniaxial compression, full dislocation loop nucleation is observed, in contrast to uniaxial tension. One of the main differences between uniaxial tension and compression is how the applied stress is resolved normal to the slip plane on which dislocations nucleate—in tension, this normal stress is tensile, and in compression, it is compressive. Last, the tension-compression asymmetry is examined as a function of loading axis orientation. Orientations with a high resolved stress normal to the slip plane on which dislocations nucleate have a larger tension-compression asymmetry with respect to dislocation nucleation than those orientations with a low resolved normal stress. The significance of this research is that the resolved stress normal to the slip plane on which dislocations nucleate plays an important role in partial (and full) dislocation loop nucleation in FCC Cu single crystals.     

Plane constrained shear of single crystal strip with two active slip systems

Within continuum dislocation theory the plane constrained shear of a single crystal strip with two active slip systems is considered. An analytical solution is found for symmetric double slip which exhibits the energetic and dissipative thresholds for dislocation nucleation, the Bauschinger translational work hardening, and the size effects. Comparison with discrete dislocation simulations shows good agreement between the discrete and continuum approaches. Numerical procedures in the general case of non-symmetric double slip are proposed.     

Influence of misfit stresses on dislocation glide in single crystal superalloys: A three-dimensional discrete dislocation dynamics study

In the characteristic $\gamma /\gamma \prime$ microstructure of single crystal superalloys, misfit stresses occur due to a significant lattice mismatch of those two phases. The magnitude of this lattice mismatch depends on the chemical composition of both phases as well as on temperature. Furthermore, the lattice mismatch of γ and $\gamma \prime$ phases can be either positive or negative in sign. The internal stresses caused by such lattice mismatch play a decisive role for the micromechanical processes that lead to the observed macroscopic athermal deformation behavior of these high-temperature alloys. Three-dimensional discrete dislocation dynamics (DDD) simulations are applied to investigate dislocation glide in γ matrix channels and shearing of $\gamma \prime$ precipitates by superdislocations under externally applied uniaxial stresses, by fully taking into account internal misfit stresses. Misfit stress fields are calculated by the fast Fourier transformation (FFT) method and hybridized with DDD simulations. For external loading along the crystallographic [001] direction of the single crystal, it was found that the different internal stress states for negative and positive lattice mismatch result in non-uniform dislocation movement and different dislocation patterns in horizontal and vertical γ matrix channels. Furthermore, positive lattice mismatch produces a lower deformation rate than negative lattice mismatch under the same tensile loading, but for an increasing magnitude of lattice mismatch, the deformation resistance always diminishes. Hence, the best deformation performance is expected to result from alloys with either small positive, or even better, vanishing lattice mismatch between γ and $\gamma \prime$ phase.     

Bending of single crystal microcantilever beams of cube orientation: Finite element model and experiments

The aim of this work is to investigate the microstructure evolution, stress-strain response and strain hardening behavior of microscale beams. For that purpose, two single crystal cantilever beams in the size dependent regime were manufactured by ion beam milling and beams were bent with an indenter device. A crystal plasticity material model for large deformations was implemented in a finite element framework to further investigate the effect of boundary constraints. Simulations were performed using bulk material properties of single crystal copper without any special treatment for the strain gradients. The difference between the slopes of the experimental and the simulated force displacement curves suggested negligible amount of strain gradient hardening compared to the statistical hardening mechanisms.     

A two-dimensional dislocation dynamics model of the plastic deformation of polycrystalline metals

Two-dimensional dislocation dynamics (2D-DD) simulations under fully periodic boundary conditions are employed to study the relation between microstructure and strength of a material. The material is modeled as an elastic continuum that contains a defect microstructure consisting of a preexisting dislocation population, dislocation sources, and grain boundaries. The mechanical response of such a material is tested by uniaxially loading it up to a certain stress and allowing it to relax until the strain rate falls below a threshold. The total plastic strain obtained for a certain stress level yields the quasi-static stress-strain curve of the material. Besides assuming Frank-Read-like dislocation sources, we also investigate the influence of a pre-existing dislocation density on the flow stress of the model material. Our results show that - despite its inherent simplifications - the 2D-DD model yields material behavior that is consistent with the classical theories of Taylor and Hall-Petch. Consequently, if set up in a proper way, these models are suited to study plastic deformation of polycrystalline materials.     

Multiscale dislocation dynamics simulations of shock compression in copper single crystal

Ultra short pulse shock wave propagation, plastic deformation and evolution of dislocations in copper single crystals with (0 0 1), (0 1 1) and (1 1 1) orientations are investigated using multiscale dislocation dynamics plasticity analyses. The effects of peak pressure, pulse duration, crystal anisotropy and the nonlinear elastic properties on the interaction between shock wave and dislocations are investigated. The results of our calculations show that the dislocation density has a power law dependence on pressure with a power of 1.70 and that the dislocation density is proportional to pulse duration and sensitive to crystal orientation. These results are in very good agreement with the analytical predications of Meyers et al. [Meyers, M.A., Gregori, F., Kad, B.K., Schneider, M.S., Kalantar, D.H., Remington, B.A., Ravichandran G., Boehly, T., Wark, J., 2003. Laser-induced shock compression of monocrystalline copper: characterization and analysis. Acta Materialia 51, 1211–1228] and the experimental results of Murr [Murr, L.E., 1981. Residual microstructure-mechanical property relationships in shock loaded metals and alloys. In: Meyers, M.A., Murr, L.E. (Eds.), Shock Waves and High Strain Rate Phenomena in Metals. Plenum, New York, pp. 607–673]. It is shown also that incorporating the effect of crystal anisotropy in the elastic properties results in orientation dependent wave speed and peak pressure. The relaxed configurations of dislocation microstructures show the formation of microbands coincident with the slip planes.     

A multi-scale computational model of crystal plasticity at submicron-to-nanometer scales

Plastic flow in crystal at submicron-to-nanometer scales involves many new interesting problems. In this paper, a unified computational model which directly combines 3D discrete dislocation dynamics (DDD) and continuum mechanics is developed to investigate the plastic behaviors at these scales. In this model, the discrete dislocation plasticity in a finite crystal is solved under a completed continuum mechanics framework: (1) an initial internal stress field is introduced to represent the preexisting stationary dislocations in the crystal; (2) the external boundary condition is handled by finite element method spontaneously; and (3) the constitutive relationship is based on the finite deformation theory of crystal plasticity, but the discrete plastic strains induced by the slip of the newly nucleated or propagating dislocations are calculated by dislocation dynamics methodology instead of phenomenological evolution equations used in conventional crystal plasticity. These discrete plastic strains are then localized to the continuum material points by a Burgers vector density function proposed by us. Various processes, such as loop dislocation evolution, dislocation junction formation etc., are simulated to verify the reliability of this computational model. Specifically, a uniaxial compression test for micro-pillars of Cu is simulated by this model to investigate the ‘dislocation starvation hardening’ observed in the recent experiment.     

On plastic anisotropy of constitutive model for rate-dependent single crystal

An algorithm for single crystals was developed and implemented to simulate plastic anisotropy using a rate-dependent slip model. The proposed procedure was a slightly modified form of single crystal constitutive model of Sarma and Zacharia. Modified Euler method, together with Newton-Raphson method was used to integrate this equation which was stable and efficient. The model together with the developed algorithm was used to study three problems. First, plastic anisotropy was examined by simulating the crystal deformation in tension and plane strain compression, respectively. Secondly, the orientation effect of some material parameters in the model and applied strain rate on plastic anisotropy for single crystal also is investigated. Thirdly, the influence of loading direction on the active slip system was discussed.     

A continuum model for intermittent deformation of single crystal micropillars

Strain bursts are often observed during compression tests of single crystal micropillars. In this work, we formulate a new continuum model that accounts for the strain bursts within the framework of crystal plasticity. The strain bursts are separated from the loading stage (nearly elastic loading) by introducing a dimensionless constant in the continuum model, and are detected by load serrations. The boundary conditions in the context of micropillar compression are studied and they are shown to be changing and unpredictable as plastic deformation proceeds. To evaluate the validity of our model, finite element simulations of the uniaxial compression tests on nickel micropillars are performed. Our simulations produce clearly visible strain bursts during the plastic flow and the produced intermittent flows are comparable with the experimental observations. For the bulk crystal, a series of strain bursts is identified in the course of plastic flow, despite an apparently smooth stress–strain response. We also show that the intermittent flow is intensified in the micrometer-scale due to both increasing numbers of the successive strain bursts and increasing amplitude of the strain burst, when the specimen size decreases. Finally, we show that the occurrences of the strain bursts are always associated with negative values of the second-order work.     

Experimental lower bounds on geometrically necessary dislocation density

A single nickel crystal is indented with a wedge indenter such that a two-dimensional deformation state with three effective plane strain slip systems is induced. The in-plane lattice rotation of the crystal lattice is measured with a three micrometer spatial resolution using Orientation Imaging Microscopy (OIM). All non-zero components of the Nye dislocation density tensor are calculated from the lattice rotation field. A rigorous analytical expression is derived for the lower bound of the total Geometrically Necessary Dislocation (GND) density. Existence and uniqueness of the lower bound are demonstrated, and the apportionment of the total GND density onto the effective individual slip systems is determined. The lower bound solution reduces to the exact solution under circumstances in which only one or two of the effective slip systems are known to have been activated. The results give insight into the active slip systems as well as the dislocation structures formed in the nickel crystal as a result of the wedge indentation.     

Selection of kernel functions used in a new constitutive equation for viscoelastic fluids

A selection of kernel functions is given to be used in a new integral constitutive equation proposed by Piau whereby the deviatoric stress is calculated from the integral of the history of the past intrinsic rate of rotation and rate of deformation tensors through a representation theorem. Piau has demonstrated the objectivity of a frame moving with a given particle whose axis are directed along the eigenvectors of the rate of deformation tensor. The use of such a framework provides a new approach in the attempt to reduce the computational difficulties associated with conventional constitutive equations written in co-deformational or co-rotational reference frames.The shear and primary normal-stress material functions and the extensional (elongational) stress growth function are defined for the proposed integral constitutive equation. These material functions are used to calculate the kernel functions using steady state, stress relaxation and stress growth data of Attané in simple shear flow for monodisperse polystyrene solutions. The shear and extensional stress growth data of Meissner for a polyethylene melt are also used to show the flexibility of the rheological model.The material functions are first written in terms of five monotonically decreasing functions of the time lag between the past and the present time. Then kernel functions are chosen such that when substituted in the new integral constitutive equation they yield the functions used to describe the data. A further condition imposed on the normalized kernel functions is that they be decreasing functions of time lag.     

Continuum framework for dislocation structure,energy and dynamics of dislocation arrays and low angle grain boundaries

We present a continuum framework for dislocation structure, energy and dynamics of dislocation arrays and low angle grain boundaries that are allowed to be nonplanar or nonequilibrium. In our continuum framework, we define a dislocation density potential function on the dislocation array surface or grain boundary to describe the orientation dependent continuous distribution of dislocations in a very simple and accurate way. The continuum formulations incorporate both the long-range dislocation interaction and the local dislocation line energy, and are derived from the discrete dislocation model. The continuum framework recovers the classical Read–Shockley energy formula when the long-range elastic fields of the low angle grain boundaries are canceled out. Applications of our continuum framework in this paper are focused on dislocation structures on static planar and nonplanar low angle grain boundaries and misfitting interfaces. We present two methods under our continuum framework for this purpose, including the method based on the Frank׳s formula and the energy minimization method. We show that for any (planar or nonplanar) low angle grain boundary, the Frank׳s formula holds if and only if the long-range stress field in the continuum model is canceled out, and it does not necessarily hold for a total energy minimum dislocation structure.     

Analysis of heterogeneous deformation and dislocation dynamics in single crystal micropillars under compression

The size dependent deformation of Cu single crystal micropillars with thickness ranging from 0.2 to 2.5 μm subjected to uniaxial compression is investigated using a Multi-scale Dislocation Dynamics Plasticity (MDDP) approach. MDDP is a hybrid elasto-viscoplastic simulation model which couples discrete dislocation dynamics at the micro-scale (software micro3d) with the macroscopic plastic deformation. Our results show that the deformation field in these micropillars is heterogeneous from the onset of plastic flow and is confined to few deformation bands, leading to the formation of ledges and stress concentrations at the surface of the specimen. Furthermore, the simulation yields a serrated stress–strain behavior consisting of discrete strain bursts that correlates well with experimental observations. The intermittent operation and stagnation of discrete dislocation arms is identified as the prominent mechanism that causes heterogeneous deformation and results in the observed macroscopic strain bursts. We show that the critical stress to bow an average maximum dislocation arm, whose length changes during deformation due to pinning events, is responsible for the observed size dependent response of the single crystals. We also reveal that hardening rates, similar to that shown experimentally, occur under relatively constant dislocation densities and are linked to dislocation stagnation due to the formation of entangled dislocation configuration and pinning sites.     

A constitutive equation for a non-Newtonian liquid

The K-BKZ constitutive equation for non-Newtonian liquids usually uses as invariants the coefficients of the characteristic equation of a certain finite strain tensor. We show that the use of the roots of this equation, the squares of the principal stretches, is more efficient. We use a product of the linear relaxation spectrum and a simple sum of powers of the principal stretches as the energy function in the K-BKZ theory. The predictions of this theory agree with Meissner's data to a good approximation over many decades.Dedicated to Professor Arthur S. Lodge on the occasion of his 70th birthday and his retirement from the University of Wisconsin.