On secondary loops in LAOS via self-intersection of Lissajous–Bowditch curves

When the shear stress measured in large amplitude oscillatory shear (LAOS) deformation is represented as a 2-D Lissajous–Bowditch curve, the corresponding trajectory can appear to self-intersect and form secondary loops. This self-intersection is a general consequence of a strongly nonlinear material response to the imposed oscillatory forcing and can be observed for various material systems and constitutive models. We derive the mathematical criteria for the formation of secondary loops, quantify the location of the apparent intersection, and furthermore suggest a qualitative physical understanding for the associated nonlinear material behavior. We show that when secondary loops appear in the viscous projection of the stress response (the 2-D plot of stress vs. strain rate), they are best interpreted by understanding the corresponding elastic response (the 2-D projection of stress vs. strain). The analysis shows clearly that sufficiently strong elastic nonlinearity is required to observe secondary loops on the conjugate viscous projection. Such a strong elastic nonlinearity physically corresponds to a nonlinear viscoelastic shear stress overshoot in which existing stress is unloaded more quickly than new deformation is accumulated. This general understanding of secondary loops in LAOS flows can be applied to various molecular configurations and microstructures such as polymer solutions, polymer melts, soft glassy materials, and other structured fluids.     

A numerical model for the nonlinear interaction of elastic waves with cracks

Microstructures such as cracks and microfractures play a significant role in the nonlinear interactions of elastic waves, but the precise mechanism of why and how this works is less clear. Here we simulate wave propagation to understand these mechanisms, complementing existing theoretical and experimental works.We implement two models, one of homogeneous nonlinear elasticity and one of perturbations to cracks, and then use these models to improve our understanding of the relative importance of cracks and intrinsic nonlinearity. We find, by modeling the perturbations in the speed of a low-amplitude P-wave caused by the propagation of a large-amplitude S-wave that the nonlinear interactions of P- and S-waves with cracks are significant when the particle motion is aligned with the normal to the crack face, resulting in a larger magnitude crack dilation. This improves our understanding of the relationship between microstructure orientations and nonlinear wave interactions to allow for better characterization of fractures for analyzing processes including earthquake response, reservoir properties, and non-destructive testing.     

位移幅值对铜镁合金微动磨损行为的影响

在研制的多功能微动磨损试验机上,开展了不同位移幅值下铜镁合金微动磨损试验,以研究位移幅值对铜镁合金微动磨损行为的影响. 微动过程中记录摩擦系数曲线与F_t-D-N曲线,利用光学显微镜(OM)、扫描电镜(SEM)、能量色散X射线光谱仪(EDS)及三维形貌仪对损伤区域进行了微观分析. 结果显示:随着位移幅值的增加,铜镁合金微动运行状态由部分滑移进入完全滑移,未发现混合滑移状态;部分滑移区中呈现由弹性变形协调逐渐向塑性变形协调转变的趋势. 磨损体积随位移幅值的增加而增加,在完全滑移区中体积损失非常严重. 在弹性变形协调的部分滑移状态下,接触表面损伤轻微,而由塑性变形协调的部分滑移状态下,接触中心出现较大切应力,疲劳裂纹扩展至接触表面导致材料剥落,接触边缘有磨粒磨损和氧化磨损的痕迹. 在完全滑移状态下,接触表面损伤主要为疲劳剥层,磨粒磨损和氧化磨损.      

The singularity in weakly nonlinear fracture mechanics

Material failure by crack propagation essentially involves a concentration of large displacement-gradients near a crack's tip, even at scales where no irreversible deformation and energy dissipation occurs. This physical situation provides the motivation for a systematic gradient expansion of general nonlinear elastic constitutive laws that goes beyond the first order displacement-gradient expansion that is the basis for linear elastic fracture mechanics (LEFM). A weakly nonlinear fracture mechanics theory was recently developed by considering displacement-gradients up to second order. The theory predicts that, at scales within a dynamic lengthscale ℓ from a crack's tip, significant logr displacements and 1/r displacement-gradient contributions arise. Whereas in LEFM the 1/r singularity generates an unbalanced force and must be discarded, we show that this singularity not only exists but is also necessary in the weakly nonlinear theory. The theory generates no spurious forces and is consistent with the notion of the autonomy of the near-tip nonlinear region. The J-integral in the weakly nonlinear theory is also shown to be path-independent, taking the same value as the linear elastic J-integral. Thus, the weakly nonlinear theory retains the key tenets of fracture mechanics, while providing excellent quantitative agreement with measurements near the tip of single propagating cracks. As ℓ is consistent with lengthscales that appear in crack tip instabilities, we suggest that this theory may serve as a promising starting point for resolving open questions in fracture dynamics.     

On using many-particle interatomic potentials to compute elastic properties of graphene and diamond

The elastic properties of diatomic crystals are considered. An approach is proposed that permits calculating the elastic characteristics of crystals by using the interatomic interaction parameters specified as many-particle potentials, i.e., potentials that take into account the effect of the environment on the diatomic interaction. The many-particle interaction is given in the general form obtained in the framework of linear elastic deformation. It is shown that, by expanding in series in small deformation parameters, a group of nonlinear potentials frequently used to model covalent structures can be reduced to this general form. An example of graphene and diamond lattices is used to determine how adequately these potentials describe the elastic characteristics of crystals.     

The solid–fluid transition in a yield stress shear thinning physical gel

We present an experimental investigation of the solid–fluid transition in a yield stress shear thinning physical gel (Carbopol^® 940) under shear. Upon a gradual increase of the external forcing, we observe three distinct deformation regimes: an elastic solid-like regime (characterized by a linear stress–strain dependence), a solid–fluid phase coexistence regime (characterized by a competition between destruction and reformation of the gel), and a purely viscous regime (characterized by a power law stress-rate of strain dependence). The competition between destruction and reformation of the gel is investigated via both systematic measurements of the dynamic elastic moduli (as a function of stress, the amplitude, and temperature) and unsteady flow ramps. The transition from solid behavior to fluid behavior displays a clear hysteresis upon increasing and decreasing values of the external forcing. We find that the deformation power corresponding to the hysteresis region scales linearly with the rate at which the material is being forced (the degree of flow unsteadiness). In the asymptotic limit of small forcing rates, our results agree well with previous steady state investigations of the yielding transition. Based on these experimental findings, we suggest an analogy between the solid–fluid transition and a first-order phase transition, e.g., the magnetization of a ferro-magnet where irreversibility and hysteresis emerge as a consequence of a phase coexistence regime. In order to get further insight into the solid–fluid transition, our experimental findings are complemented by a simple kinetic model that qualitatively describes the structural hysteresis observed in our rheological experiments. The model is fairly well validated against oscillatory flow data by a partial reconstruction of the Pipkin space of the material’s response and its nonlinear spectral behavior.     

Invertibility and Weak Continuity of the Determinant for the Modelling of Cavitation and Fracture in Nonlinear Elasticity

In this paper, we present and analyze a variational model in nonlinear elasticity that allows for cavitation and fracture. The main idea in unifying the theories of cavitation and fracture is to regard both cavities and cracks as phenomena of the creation of a new surface. Accordingly, we define a functional that measures the area of the created surface. This functional has relationships with the theory of Cartesian currents. We show that the boundedness of that functional implies sequential weak continuity of the determinant of the deformation gradient, and that the weak limit of one-to-one almost everywhere deformations is also one-to-one almost everywhere. We then use these results to obtain the existence of minimizers of variational models that incorporate elastic energy and this created surface energy, taking into account orientation-preserving and non-interpenetration conditions.     

Nonlinear laws of dry friction in contact problems of linear theory of elasticity

A problem of elastic body deformation with conditions of dry friction imposed at the boundary is considered. Various friction laws are studied, including linear oneparameter and nonlinear twoparameter laws. A general view of a nonlinear function with two constants is suggested, which determines the friction force as a function of normal pressure. The problem of elastic plate compression by rough infinite plates for a variety of friction conditions on the contact surfaces is solved. The plate equations are employed, which make it possible to specify arbitrary conditions on the front faces without reducing the order of differential equations. The unknown boundary of ideal contact zones and sliding zones is determined. Solutions obtained by using various friction conditions on the contact surfaces are compared.     

An equation of motion for a thick viscoelastic plate

In this paper an equation of motion is presented for a general thick viscoelastic plate, including the effects of shear deformation, extrusion deformation and rotatory inertia. This equation is the generalization of equations of motion for the corresponding thick elastic plate, and it can be degenerated into several types of equations for various special cases.     

考虑粗糙结合面弹塑性接触的黏滑摩擦建模

提出一种同时考虑粗糙面上微凸体弹性变形和塑性接触的切向黏滑摩擦建模方法。采用Hertz弹性理论和Mindlin解描述弹性接触微凸体的切向载荷和相对变形的关系;采用AF(Abbott-Firstone)塑性理论和Fujimoto模型描述塑性接触微凸体切向载荷和相对变形的关系。再利用GW(Greenwood-Williamson)模型统计分析方法建立粗糙表面切向载荷和相对变形之间的关系。将模型与仅考虑微凸体弹性接触情况的模型进行对比,并研究了不同塑性指数对切向载荷和相对变形关系的影响。结果表明：与完全弹性接触模型相比,本文模型引入了塑性接触理论,能够更好地描述粗糙表面切向载荷和相对变形关系,并且考虑不同接触条件下弹性变形微凸体和塑性变形微凸体对切向接触载荷的贡献,在微滑移阶段,主要由弹性接触变形影响,而在进入宏观滑移阶段之后,切向行为主要由塑性变形影响。界面切向载荷由黏着和滑移接触作用共同决定,随着切向变形的增加,滑移接触力逐渐增加,而黏着接触力先增加后减少,反映了界面由微滑移逐渐向宏滑移演化的过程。随着塑性指数的增加,粗糙面上发生塑性接触的微凸体数目逐渐增加,切向黏滑行为主要受到塑性接触特征的控制。     

具有非线性homologous变形约束的桁架结构形态分析

研究了具有非线性homologous变形约束条件的桁架结构形态分析问题。在已有的线性homologous变形约束桁架形态分析的基础上,将结构的节点分成三类:homologous变形约束节点,形状可变节点和边界点。运用Moore-Penrose广义逆矩阵性质,将基础方程组解的存在条件表示为包含形状可变节点未知坐标的非线性方程组,为采用Newton-Raphson方法求解非线性方程组,对AA (A为任意矩阵,A 为A的Moore-Penrose广义逆矩阵)求偏导数,找到了满足保型要求的形态,给出的桁架算例说明了本文方法的有效性。     

On the asymptotics of the stress and strain fields near a crack tip

The problem is solved under the plane strain conditions for a crack of general form, which in general is neither a mode I nor a mode II crack. We assume that the strains are small and the material is nonlinearly elastic. The mathematical statement of the problem is reduced to the eigenvalue problem for a system of ordinary nonlinear differential equations. Its solution is obtained numerically. We show that, for an incompressible material with power-law relations between the stress and strain deviators, the solution (the well-known HRR-asymptotics [1, 2]) exists only for mode I and II cracks. In the general case, we can only speak of approximate solutions. A similar conclusion can be made for different-modulus materials. We analyze the results of the preceding papers [1–7], where specific cases of the problem were considered.     

Finite amplitude,periodic motion of a body supported by arbitrary isotropic,elastic shear mountings

The undamped, finite amplitude, periodic motion of a load supported symmetrically by arbitrary isotropic, elastic shear mountings is investigated. Conditions on the shear response function sufficient to guarantee periodic motions for finite shearing with arbitrary initial data are provided. Some general results applicable for all simple shearing oscillators in the class are derived and illustrated graphically. The mechanical response of the general nonlinear shearing oscillator is compared with the response of a certain linear oscillator of comparable design. As consequence, certain static and dynamic aspects of the motion of an arbitrary nonlinear oscillator supported by shear springs are compared with those of a simple, linear oscillator for which the response is well-known and readily determined for the same initial data. The effect of a finite static shear deformation on the frequency equation for superimposed, small amplitude vibrations of the load is examined. The general analysis is applied to a class of hyperelastic biological tissues; and the frequency relation for finite amplitude oscillations of a load supported by soft tissue is derived. The finite amplitude oscillatory shearing of a general isotropic elastic continuum is described; and three universal relations connecting the stress and the oscillatory shearing deformation for every isotropic elastic material are presented.     

Complex unloading behavior: Nature of the deformation and its consistent constitutive representation

Complex (nonlinear) unloading behavior following plastic straining has been reported as a significant challenge to accurate springback prediction. More fundamentally, the nature of the unloading deformation has not been resolved, being variously attributed to nonlinear/reduced modulus elasticity or to inelastic/“microplastic” effects. Unloading-and-reloading experiments following tensile deformation showed that a special component of strain, deemed here “Quasi-Plastic-Elastic” (“QPE”) strain, has four characteristics. (1) It is recoverable, like elastic deformation. (2) It dissipates work, like plastic deformation. (3) It is rate-independent, in the strain rate range 10⁻⁴-10⁻²/s, contrary to some models of anelasticity to which the unloading modulus effect has been attributed. (4) To first order, the evolution of plastic properties occurs during QPE deformation. These characteristics are as expected for a mechanism of dislocation pile-up and relaxation. A consistent, general, continuum constitutive model was derived incorporating elastic, plastic, and QPE deformation. Using some aspects of two-yield-function approaches with unique modifications to incorporate QPE, the model was implemented in a finite element program with parameters determined for dual-phase steel and applied to draw-bend springback. Significant differences were found compared with standard simulations or ones incorporating modulus reduction. The proposed constitutive approach can be used with a variety of elastic and plastic models to treat the nonlinear unloading and reloading of metals consistently for general three-dimensional problems.     

Finite strain solutions in compressible isotropic elasticity

Three classes of compressible isotropic elastic solids are introduced, for each of which the strain energy, expressed as a function of the three principal invariants of the stretch tensors, is linear in two of its arguments and nonlinear in the third argument. One of these is the class of harmonic materials. Several deformation fields are examined, for which the governing equations, for general compressible isotropic elastic response, reduce to a nonlinear ordinary differential equation. For the three special classes of materials, this differential equation may be solved in closed form, giving a deformation field which is independent of the material response function, or its solution may be reduced to a single quadrature, involving the nonlinear material response function.     

Some Remarks on the Effects of Inertia and Viscous Dissipation in the Onset of Cavitation in Rubber

Through direct comparisons with experiments, Lefèvre et al. (Int. J. Frac. 192:1–23, 2015) have recently confirmed the prevailing belief that the nonlinear elastic properties of rubber play a significant role in the so-called phenomenon of cavitation—that is, the sudden growth of inherent defects in rubber into large enclosed cavities/cracks in response to external stimuli. These comparisons have also made it plain that cavitation in rubber is first and foremost a fracture process that may possibly depend, in addition to the nonlinear elastic properties of the rubber, on inertial effects and/or on the viscous dissipation innate to rubber. This is because the growth of defects into large cavities/cracks is locally in time an extremely fast process.The purpose of this Note is to provide insight into the relevance of inertial and viscous dissipation effects on the onset of cavitation in rubber. To this end, leaving fracture properties aside, we consider the basic problem of the radially symmetric dynamic deformation of a spherical defect embedded at the center of a ball made up of an isotropic incompressible nonlinear viscoelastic solid that is subjected to external hydrostatic loading. Specifically, the defect is taken to be vacuous and the viscoelastic behavior of the solid is characterized by a fairly general class of constitutive relations given in terms of two thermodynamic potentials—namely, a free energy function describing the nonlinear elasticity of the solid and a dissipation potential describing its viscous dissipation—which has been shown to be capable to describe the mechanical response of a broad variety of rubbers over wide ranges of deformations and deformations rates. It is found that the onset of cavitation is not affected by inertial effects so long as the external loads are not applied at a high rate. On the other hand, even when the external loads are applied quasi-statically, viscous dissipation can greatly affect the critical values of the applied loads at which cavitation ensues.     

Effective elastic properties of flexible chiral honeycomb cores including geometrically nonlinear effects

Flexible chiral honeycomb cores generally exhibit nonlinear elastic properties in response to large geometric deformation, which are suited for the design of morphing aerospace structures. However, owing to their complex structure, it is standard to replace the actual core structure with a homogenized core material presenting reasonably equivalent elastic properties in an effort to increase the speed and efficiency of analyzing the mechanical properties of chiral honeycomb sandwich structures. As such, a convenient and efficient method is required to evaluate the effective elastic properties of flexible chiral honeycomb cores under conditions of large deformation. The present work develops an analytical expression for the effective elastic modulus based on a deformable cantilever beam under large deformation. Firstly, Euler–Bernoulli beam theory and micropolar theory are used to analyze the deformation characteristics of chiral honeycombs, and to calculate the effective elastic modulus under small deformation. On that basis, the expression for the effective elastic modulus is improved by including the stretching deformation of the chiral honeycomb structure for a unit cell under conditions of large deformation. The effective elastic moduli calculated by the respective analytical expressions are compared with the results of finite element analysis. The results indicate that the analytical expression obtained under consideration of the geometric nonlinearity is more suitable than the linear expressions for flexible chiral honeycomb cores under conditions of high strain and low elastic modulus.     

A three dimensional field formulation,and isogeometric solutions to point and line defects using Toupin's theory of gradient elasticity at finite strains

We present a field formulation for defects that draws from the classical representation of the cores as force dipoles. We write these dipoles as singular distributions. Exploiting the key insight that the variational setting is the only appropriate one for the theory of distributions, we arrive at universally applicable weak forms for defects in nonlinear elasticity. Remarkably, the standard, Galerkin finite element method yields numerical solutions for the elastic fields of defects that, when parameterized suitably, match very well with classical, linearized elasticity solutions. The true potential of our approach, however, lies in its easy extension to generate solutions to elastic fields of defects in the regime of nonlinear elasticity, and even more notably for Toupin's theory of gradient elasticity at finite strains (Toupin Arch. Ration. Mech. Anal., 11 (1962) 385). In computing these solutions we adopt recent numerical work on an isogeometric analytic framework that enabled the first three-dimensional solutions to general boundary value problems of Toupin's theory (Rudraraju et al. Comput. Methods Appl. Mech. Eng., 278 (2014) 705). We first present exhaustive solutions to point defects, edge and screw dislocations, and a study on the energetics of interacting dislocations. Then, to demonstrate the generality and potential of our treatment, we apply it to other complex dislocation configurations, including loops and low-angle grain boundaries.     

Deformations of an elastic,internally constrained material. Part 1: Homogeneous deformations

The nonlinear elastic response of a class of materials for which the deformation is subject to an internal material constraint described in experiments by James F. Bell on the finite deformation of a variety of metals is investigated. The purely kinematical consequences of the Bell constraint are discussed, and restrictions on the full range of compatible deformations are presented in geometrical terms. Then various forms of the constitutive equation relating the stress and stretch tensors for an isotropic elastic Bell material are presented. Inequalities on the mechanical response functions are introduced. The importance of these in applications is demonstrated in several examples throughout the paper.This paper focuses on homogeneous deformations. In a simple illustration of the theory, a generalized form of Bell's empirical rule for uniaxial loading is derived, and some peculiarities in the response under all-around compressive loading are discussed. General formulae for universal relations possible in an isotropic elastic, Bell constrained material are presented. A simple method for the determination of the left stretch tensor for essentially plane problems is illustrated in the solution of the problem of pure shear of a materially uniform rectangular block. A general formula which includes the empirical rule found in pure shear experiments by Bell is derived as a special case. The whole apparatus is then applied in the solution of the general problem of a homogeneous simple shear superimposed on a uniform triaxial stretch; and the great variety of results possible in an isotropic, elastic Bell material is illustrated. The problem of the finite torsion and extension of a thin-walled cylindrical tube is investigated. The results are shown to be consistent with Bell's data for which the rigid body rotation is found to be quite small compared with the gross deformation of the tube. Several universal formulas relating various kinds of stress components to the deformation independently of the material response functions are derived, including a universal rule relating the axial force to the torque.Constitutive equations for hyperelastic Bell materials are derived. The empirical work function studied by Bell is introduced; and a new constitutive equation is derived, which we name Bell's law. On the basis of this law, we then derive exactly Bell's parabolic laws for uniaxial loading and for pure shear. Also, form Bell's law, a simple constitutive equation relating Bell's deviatoric stress tensor to his finite deviatoric strain tensor is obtained. We thereby derive Bell's invariant parabolic law relating the deviatoric stress intensity to the corresponding strain intensity; and, finally, Bell's fundamental law for the work function expressed in these terms is recovered. This rule is the foundation for all of Bell's own theoretical study of the isotropic materials cataloged in his finite strain experiments on metals, all consistent with the internal material constraint studied here.