Strain gradient effects on cyclic plasticity

Size effects on the cyclic shear response are studied numerically using a recent higher order strain gradient visco-plasticity theory accounting for both dissipative and energetic gradient hardening. Numerical investigations of the response under cyclic pure shear and shear of a finite slab between rigid platens have been carried out, using the finite element method. It is shown for elastic-perfectly plastic solids how dissipative gradient effects lead to increased yield strength, whereas energetic gradient contributions lead to increased hardening as well as a Bauschinger effect. For linearly hardening materials it is quantified how dissipative and energetic gradient effects promote hardening above that of conventional predictions. Usually, increased hardening is attributed to energetic gradient effects, but here it is found that also dissipative gradient effects lead to additional hardening in the presence of conventional material hardening. Furthermore, it is shown that dissipative gradient effects can lead to both an increase and a decrease in the dissipation per load cycle depending on the magnitude of the dissipative length parameter, whereas energetic gradient effects lead to decreasing dissipation for increasing energetic length parameter. For dissipative gradient effects it is found that dissipation has a maximum value for some none zero value of the material length parameter, which depends on the magnitude of the deformation cycles.     

A multi-phase theory for the detonation of granular explosives containing an arbitrary number of solid components

Multiphase continuum models are commonly used to predict the shock, combustion and detonation behavior of granular energetic mixtures containing solid reactants and gaseous products. These models often include phase interaction terms that formally satisfy the strong form of the Second Law of Thermodynamics and provide flexibility in distributing dissipation between phases arising from non-equilibrium phenomena. This work presents a thermodynamically compatible constitutive theory for reactive systems containing an arbitrary number of solid components. The theory represents a rigorous extension of the two-phase theory formulated by Bdzil et al., based on the well-studied Baer–Nunziato model. Forms of the gas–solid and solid–solid interphase sources suggested by general reactions of type A → B are considered, where the combustion processes discussed in Bdzil et al. are treated as a special case. The model energetics are augmented by supplemental evolutionary equations that track local changes in phase temperatures due to dissipative and transport processes allowing for the identification of dominant energetic processes. This capability provides a mean to identify system parameters (e.g., metal particle size and mass fraction in metalized energetic mixtures) which optimize performance metrics. Detonation predictions are given for mixtures of granular HMX and aluminum to demonstrate model features and to highlight the effect of aluminum particle self-heating by oxidation on detonation. Predicted spatial profiles for mechanical fields, and the heating contributions from individual dissipative processes, illustrate how aluminum particle size can affect the coupling of oxidative heating to the explosive reaction zone.     

Modeling of polycrystals using a gradient crystal plasticity theory that includes dissipative micro-stresses

This study investigates thermodynamically consistent dissipative hardening in gradient crystal plasticity in a large-deformation context. A viscoplastic model which accounts for constitutive dependence on the slip, the slip gradient as well as the slip rate gradient is presented. The model is an extension of that due to Gurtin (Gurtin, M. E., J. Mech. Phys. Solids, 52 (2004) 2545–2568 and Gurtin, M. E., J. Mech. Phys. Solids, 56 (2008) 640–662)), and is guided by the viscoplastic model and algorithm of Ekh et al. (Ekh, M., Grymer, M., Runesson, K. and Svedberg, T., Int. J. Numer. Meths Engng, 72 (2007) 197–220) whose governing equations are equivalent to those of Gurtin for the purely energetic case. In contrast to the Gurtin formulation and in line with that due to Ekh et al., viscoplasticity in the present model is accounted for through a Perzyna-type regularization. The resulting theory includes three different types of hardening: standard isotropic hardening is incorporated as well as energetic hardening driven by the slip gradient. In addition, as a third type, dissipative hardening associated with plastic strain rate gradients is included. Numerical computations are carried out and discussed for the large strain, viscoplastic model with non-zero dissipative backstress.     

Homogenized Interface Model Describing Inhomogeneities Located on a Surface

We study the influence of heterogeneities located near a planar surface on the elastic response of a three-dimensional elastic medium. These heterogeneities can be either reinforcements, like steel reinforcements in concrete, or defects, like micro-cracks periodically distributed. We prove that their influence is of the second order from an energetic viewpoint. Then, we propose an ??up to second order effective model?? in which the influence of the heterogeneities is given by a surface energy contribution involving both the jump of displacement across the surface and the tangential strain components on the surface. The effective coefficients entering in the definition of the surface energy are obtained by solving ??elementary?? elastic problems formulated on an infinite representative cell containing the defects. We analyze this model, in particular the properties of the effective surface coefficients, and establish its coherence with limit models previously described in the literature for stiff or soft interfaces. This approach is finally applied to several kinds of heterogeneities.     

Microcanonical Entropy and Mesoscale Dislocation Mechanics and Plasticity

A methodology is devised to utilize the statistical mechanical entropy of an isolated, constrained atomistic system to define constitutive response functions for the dissipative driving-force and energetic fields in continuum thermomechanics. A thermodynamic model of dislocation mechanics is discussed as an example. Primary outcomes are constitutive relations for the back-stress tensor and the Cauchy stress tensor in terms of the elastic distortion, mass density, polar dislocation density, and the scalar statistical density.     

A theory for grain boundaries with strain-gradient plasticity

In this work, the effect of the material microstructural interface between two materials (i.e., grain boundary in polycrystalls) is adopted into a thermodynamic-based higher order strain gradient plasticity framework. The developed grain boundary flow rule accounts for the energy storage at the grain boundary due to the dislocation pile up as well as energy dissipation caused by the dislocation transfer through the grain boundary. The theory is developed based on the decomposition of the thermodynamic conjugate forces into energetic and dissipative counterparts which provides the constitutive equations to have both energetic and dissipative gradient length scales for the grain and grain boundary. The numerical solution for the proposed framework is also presented here within the finite element context. The material parameters of the gradient framework are also calibrated using an extensive set of micro-scale experimental measurements of thin metal films over a wide range of size and temperature of the samples.     

Computational modeling of size-dependent superelasticity of shape memory alloys

We propose a nonlocal continuum model to describe the size-dependent superelastic responses observed in recent experiments of shape memory alloys. The modeling approach extends a superelasticity formulation based on the martensitic volume fraction, and combines it with gradient plasticity theories. Size effects are incorporated through two internal length scales, an energetic length scale and a dissipative length scale, which correspond to the gradient terms in the free energy and the dissipation, respectively. We also propose a computational framework based on a variational formulation to solve the coupled governing equations resulting from the nonlocal superelastic model. Within this framework, a robust and scalable algorithm is implemented for large scale three-dimensional problems. A numerical study of the grain boundary constraint effect shows that the model is able to capture the size-dependent stress hysteresis and strain hardening during the loading and unloading cycles in polycrystalline SMAs.     

球面聚心气相爆轰波传播过程的数值研究

应用频散可控耗散格式对球面聚心气相爆轰波的传播过程进行了数值模拟 研究. 通过跟踪波阵面上压力和温度变化，分析这些参数在爆轰波传播过程中的演变规律， 及其与几何尺度和初始条件之间的依赖关系. 研究结果表明，仅在远场波面压力的变化近似 只依赖于r/R，对称中心附近则需要考虑初始半径$R$的影响；波面压力与初始压力的变化呈 线性关系；汇聚过程中温度升高比压力慢得多.     

Nonlinear seismic waves: A model for site effects

The aim of this paper is to propose a possible mathematical model of site effects that occur when seismic waves propagate through a sediment filled basin. The model is based on the mechanical properties of the medium (that we consider as a granular material) through which the seismic waves propagate. By looking for asymptotic solutions having the features of a progressive wave, we derive an evolution equation which is a modified Korteweg–deVries–Burgers equation containing also a nonlinear dissipative term. This equation is integrated numerically and the modelled site amplification is evaluated by using the smoothed spectral ratio between the propagated profile of the wave and the initial one.     

碳氢燃料点火燃烧的简化化学反应动力学模型

基于``准稳态'方法建立了一套复杂化学反应动力学模型简化方法和相应的软件SPARCK. 并以3种典型的碳氢燃料------甲烷、乙烯和庚烷为研究对象，从甲烷点火燃烧的GRI2.11详 细基元反应动力学模型出发简化得出了包含14个组分10步总包反应形式的简化化学反应动 力学模型，从乙烯燃烧的51组分365详细基元反应模型出发简化得出了包含20个组分16 步总包反应形式的简化化学反应动力学模型，从庚烷点火燃烧的160组分1540详细基元反 应模型出发简化得出了包含26个组分22步总包反应形式的简化化学反应动力学模型. 通过 对典型激波管试验的结果对比可以看出：得到的简化反应动力学模型能较为有效地再现 详细基元反应模型的反应机理，具有较高的计算精度. 在工程计算中有较好的应用前景.     

Dislocation pile-ups in bicrystals within continuum dislocation theory

Within continuum dislocation theory the plastic deformation of bicrystals under a mixed deformation of plane constrained uniaxial extension and shear is investigated with regard to the nucleation of dislocations and the dislocation pile-up near the phase boundaries of a model bicrystal with one active slip system within each single crystal. For plane uniaxial extension, we present a closed-form analytical solution for the evolution of the plastic distortion and of the dislocation network in the case of symmetric slip planes (i.e. for twins), which exhibits an energetic as well as a dissipative threshold for the dislocation nucleation. The general solution for non-symmetric slip systems is obtained numerically. For a combined deformation of extension and shear, we analyze the possibility of linearly superposing results obtained for both loading cases independently. All solutions presented in this paper also display the Bauschinger effect of translational work hardening and a size effect typical to problems of crystal plasticity.     

Global Classical Solutions for Partially Dissipative Hyperbolic System of Balance Laws

The basic existence theory of Kato and Majda enables us to obtain local-in-time classical solutions to generally quasilinear hyperbolic systems in the framework of Sobolev spaces (in x) with higher regularity. However, it remains a challenging open problem whether classical solutions still preserve well-posedness in the case of critical regularity. This paper is concerned with partially dissipative hyperbolic system of balance laws. Under the entropy dissipative assumption, we establish the local well-posedness and blow-up criterion of classical solutions in the framework of Besov spaces with critical regularity with the aid of the standard iteration argument and Friedrichs’ regularization method. Then we explore the theory of function spaces and develop an elementary fact that indicates the relation between homogeneous and inhomogeneous Chemin–Lerner spaces (mixed space-time Besov spaces). This fact allows us to capture the dissipation rates generated from the partial dissipative source term and further obtain the global well-posedness and stability by assuming at all times the Shizuta–Kawashima algebraic condition. As a direct application, the corresponding well-posedness and stability of classical solutions to the compressible Euler equations with damping are also obtained.     

On torsion of a single crystal rod

The paper aims at calculating the dislocation distribution inside a single crystal rod loaded in torsion within the framework of continuum dislocation theory. We construct an explicit analytical solution of this problem in terms of the modified Bessel and hypergeometric functions. The interesting features of this solution are the energetic and dissipative thresholds for dislocation nucleation, the translational work hardening, and the size effect. The comparison with experimental results shows quite good agreement of the torque-twist curves for small up to moderate twists.     

Discontinuity-induced bifurcations in systems with impacts and friction: Discontinuities in the impact law

This paper concerns the non-smooth dynamics of planar mechanical systems with isolated contact in the presence of Coulomb friction. Following Stronge [Impact Mechanics, Cambridge University Press, Cambridge, 2000], a set of closed-form analytic formulae is derived for a rigid-body impact law based on an energetic coefficient of restitution and a resolution of the impact phase into distinct segments of relative slip and stick. Thus, the impact behavior is consistent both with the assumption of Coulomb friction and with the dissipative nature of impacts. The analysis highlights the presence of boundaries between open regions of initial conditions and parameter values corresponding to distinct forms of the impact law and investigates the smoothness properties of the impact law across these boundaries. It is shown how discontinuities in the impact law are associated with discontinuity-induced bifurcations of periodic trajectories, including non-smooth folds and persistence scenarios. Numerical analysis of an example mechanical model is used to illustrate and validate the conclusions.     

Time-dependent behavior of passive skeletal muscle

An isotropic three-dimensional nonlinear viscoelastic model is developed to simulate the time-dependent behavior of passive skeletal muscle. The development of the model is stimulated by experimental data that characterize the response during simple uniaxial stress cyclic loading and unloading. Of particular interest is the rate-dependent response, the recovery of muscle properties from the preconditioned to the unconditioned state and stress relaxation at constant stretch during loading and unloading. The model considers the material to be a composite of a nonlinear hyperelastic component in parallel with a nonlinear dissipative component. The strain energy and the corresponding stress measures are separated additively into hyperelastic and dissipative parts. In contrast to standard nonlinear inelastic models, here the dissipative component is modeled using an evolution equation that combines rate-independent and rate-dependent responses smoothly with no finite elastic range. Large deformation evolution equations for the distortional deformations in the elastic and in the dissipative component are presented. A robust, strongly objective numerical integration algorithm is used to model rate-dependent and rate-independent inelastic responses. The constitutive formulation is specialized to simulate the experimental data. The nonlinear viscoelastic model accurately represents the time-dependent passive response of skeletal muscle.     

Plane constrained uniaxial extension of a single crystal strip

Within continuum dislocation theory the plane constrained uniaxial extension of a single crystal strip deforming in single or double slip is analyzed. For the single and symmetric double slip, the closed-form analytical solutions are found which exhibits the energetic and dissipative thresholds for dislocation nucleation, the Bauschinger translational work hardening, and the size effect. Numerical solutions for the non-symmetric double slip are obtained by finite element procedures.     

Interfacial energy effects within the framework of strain gradient plasticity

In the framework of strain gradient plasticity, a solid body with boundary surface playing the role of a dissipative boundary layer endowed with surface tension and surface energy, is addressed. Using the so-called residual-based gradient plasticity theory, the state equations and the higher order boundary conditions are derived quite naturally for both the bulk material and the boundary layer. A phenomenological constitutive model is envisioned, in which the bulk material and the boundary layer obey (rate independent associative) coupled plasticity evolution laws, with kinematic hardening laws of differential nature for the bulk material, but of nondifferential nature for the layer. A combined global maximum dissipation principle is shown to hold. The higher order boundary conditions are discussed in details and categorized in relation to some peculiar features of the boundary surface, and their basic role in the coupling of the bulk/layer plasticity evolution laws is pointed out. The case of an internal interface is also studied. An illustrative example relating to a shear model exhibiting energetic size effects is presented. The theory provides a unified view on gradient plasticity with interfacial energy effects.     

A mechanism of transformational plasticity

In order to understand the phenomenon of reversible plasticity exhibited by shape memory alloys and other smart materials, we study an elementary prototypical model. Building on an original idea of Müller and Villaggio [17], we consider an inhomogeneous ensemble of bi-stable elements connected in series and loaded in a soft device. To interpret the fine structure of the hysteresis loops observed experimentally, we assume that the dynamics is maximally dissipative and investigate different evolutiona ry strategies for a “driven” system with external force changing quasi-statically. Our main result is that the inhomogeneity of the elastic properties leads to a distinctive hardening with serrations of a Portevin-Le Chatelier type and produces a realistic memory structure characterized by the “congruency” and “return point memory” properties. Received December 28, 2001 / Published online June 4, 2002 Dedicated to Ingo Müller on the occasion of his 65th birthday Communicated by Kolumban Hutter, Darmstadt     

Grain boundary interface mechanics in strain gradient crystal plasticity

Interactions between dislocations and grain boundaries play an important role in the plastic deformation of polycrystalline metals. Capturing accurately the behaviour of these internal interfaces is particularly important for applications where the relative grain boundary fraction is significant, such as ultra fine-grained metals, thin films and micro-devices. Incorporating these micro-scale interactions (which are sensitive to a number of dislocation, interface and crystallographic parameters) within a macro-scale crystal plasticity model poses a challenge. The innovative features in the present paper include (i) the formulation of a thermodynamically consistent grain boundary interface model within a microstructurally motivated strain gradient crystal plasticity framework, (ii) the presence of intra-grain slip system coupling through a microstructurally derived internal stress, (iii) the incorporation of inter-grain slip system coupling via an interface energy accounting for both the magnitude and direction of contributions to the residual defect from all slip systems in the two neighbouring grains, and (iv) the numerical implementation of the grain boundary model to directly investigate the influence of the interface constitutive parameters on plastic deformation. The model problem of a bicrystal deforming in plane strain is analysed. The influence of dissipative and energetic interface hardening, grain misorientation, asymmetry in the grain orientations and the grain size are systematically investigated. In each case, the crystal response is compared with reference calculations with grain boundaries that are either ‘microhard’ (impenetrable to dislocations) or ‘microfree’ (an infinite dislocation sink).