Constitutive modeling of discontinuities by means of discrete and continuum approximations and damage models

In this paper the mathematical modeling of discontinuities using the discrete approximation and the continuum approximation with weak discontinuities is presented. First, the kinematics of discontinuities is discussed, then two constitutive models based on the continuum damage mechanics theory are developed. The first model is an isotropic damage model and is used in the discrete approximation. The second model is an anisotropic damage model and is used in the continuum approximation. These models are characterized for weighing the mode of failure in the failure criterion. An energy analysis is proposed to establish the equations that relate the parameters of both constitutive models; the fulfillment of the involved equations guarantee that both models are energetically equivalent. It is concluded that the proposed models are suitable to reproduce the constitutive behavior of discontinuities.     

Dynamic steady crack growth in elastic-plastic solids—propagation of strong discontinuities

The near-tip field of a mode I crack growing steadily under plane strain conditions is studied. A key issue is whether strong discontinuities can propagate under dynamic conditions. Theories which impose rather restrictive assumptions on the structure of an admissible deformation path through a dynamically propagating discontinuity have been proposed recently. Asymptotic solutions for dynamic crack growth, based on such theories, do not contain any discontinuities. In the present work a broader family of deformation paths is considered and we show that a discontinuity can propagate dynamically without violating any of the mechanical constitutive relations of the material. The proposed theory for the propagation of strong discontinuities is corroborated by very detailed finite element calculations. The latter shows a plane of strong discontinuity emanating from the crack tip (with its normal pointing in the direction of crack advance) and moving with the tip. Elastic unloading ahead of and/or behind the plane of discontinuity and behind the crack tip have also been observed.The numerical investigation is performed within the framework of a boundary layer formulation whereby the remote loading is fully specified by the first two terms in the asymptotic solution of the elasto-dynamic crack tip field, characterized by K₁, and T. It is shown that the family of near-tip fields, associated with a given crack speed, can be arranged into a one-parameter field based on a characteristic length, L_g, which scales with the smallest dimension of the plastic zone. This extends a previous result for quasi-static crack growth.     

Propagation of discontinuity waves in complex rheological media with viscous properties

The propagation of discontinuity waves of various order in rheological media is examined. It is assumed that the region of discontinuity of values can be represented by an intermediate layer of infinitesimal thickness. By means of this representation, results can be obtained for a rather wide class of continuous media with viscous properties, which generalize Duhem's results. The first integrals of the laws of momentum and energy conservation are obtained, which hold inside the intermediate layer at a shock wave.It is shown that when viscosity elements are introduced in a special way into the rheological model of a continuous medium, discontinuity waves of any order are propagated in the medium, and that at the surface of a strong discontinuity in a heat-conducting medium, the temperature is continuous. Additional conditions for strain discontinuities at the viscosity elements are obtained. For certain inclusions of the viscosity elements into the rheological model discontinuity waves do not propagate; instead there is merely a weak discontinuity surface which acts as an interface between the flow region of the continuous medium and the region in the state of rest. Contact discontinuities can occur in any continuous medium.The possible existence of a geometrical discontinuity surface in a viscous gas was examined first by Duhem [1]. He established that singluar strong-discontinuity surfaces cannot take place in a viscous gas. However, if one assumes that the velocity and temperature are continuous in the passage through a singular surface, only contact discontinuities are possible [2].     

On the equivalence between traction- and stress-based approaches for the modeling of localized failure in solids

This work investigates systematically traction- and stress-based approaches for the modeling of strong and regularized discontinuities induced by localized failure in solids. Two complementary methodologies, i.e., discontinuities localized in an elastic solid and strain localization of an inelastic softening solid, are addressed. In the former it is assumed a priori that the discontinuity forms with a continuous stress field and along the known orientation. A traction-based failure criterion is introduced to characterize the discontinuity and the orientation is determined from Mohr's maximization postulate. If the displacement jumps are retained as independent variables, the strong/regularized discontinuity approaches follow, requiring constitutive models for both the bulk and discontinuity. Elimination of the displacement jumps at the material point level results in the embedded/smeared discontinuity approaches in which an overall inelastic constitutive model fulfilling the static constraint suffices. The second methodology is then adopted to check whether the assumed strain localization can occur and identify its consequences on the resulting approaches. The kinematic constraint guaranteeing stress boundedness and continuity upon strain localization is established for general inelastic softening solids. Application to a unified stress-based elastoplastic damage model naturally yields all the ingredients of a localized model for the discontinuity (band), justifying the first methodology. Two dual but not necessarily equivalent approaches, i.e., the traction-based elastoplastic damage model and the stress-based projected discontinuity model, are identified. The former is equivalent to the embedded and smeared discontinuity approaches, whereas in the later the discontinuity orientation and associated failure criterion are determined consistently from the kinematic constraint rather than given a priori. The bi-directional connections and equivalence conditions between the traction- and stress-based approaches are classified. Closed-form results under plane stress condition are also given. A generic failure criterion of either elliptic, parabolic or hyperbolic type is analyzed in a unified manner, with the classical von Mises (J₂), Drucker–Prager, Mohr–Coulomb and many other frequently employed criteria recovered as its particular cases.     

Critical distances for the formation of strong discontinuities in fluid-saturated solids

An increase in the stiffness of a solid in compression is known to lead to the steepening of the profiles of compression waves and, as a consequence, to the formation of strong discontinuities from continuous waves propagating in the solid. In this paper, the critical distance required for a continuous wave to turn into a shock wave is calculated from the evolution equation for a weak discontinuity (acceleration wave) propagating into a quiescent region. Infinite growth of the amplitude of an acceleration wave in a finite time signifies the transition to a strong discontinuity. Relations between the critical distances for plane, cylindrical and spherical waves are established. Numerical examples are presented for a particular case of the pressure-dependent stiffness typical of granular solids such as sand or soil, with emphasis placed on the influence of a small amount of free gas in the pore fluid.     

Micromechanics based second gradient continuum theory for shear band modeling in cohesive granular materials following damage elasticity

Gradient theories, as a regularized continuum mechanics approach, have found wide applications for modeling strain localization failure process. This paper presents a second gradient stress–strain damage elasticity theory based upon the method of virtual power. The theory considers the strain gradient and its conjugated double stresses. Instead of introducing an intrinsic material length scale into the constitutive law in an ad hoc fashion, a microstructural granular mechanics approach is applied to derive the higher-order constitutive coefficients such that the internal length scale parameter reflects the natural granularity of the underlying material microstructure. The derivations of the required damage constitutive relationships, the strong form governing equations as well as its weak form for the second gradient model are described. The recently popularized Element-Free Galerkin (EFG) method is then employed to discretize the weak form equilibrium equation for accommodating the resultant higher-order continuity requirements and further handling the mesh sensitivity problem. Numerical examples for shear band simulations show that the proposed second gradient continuum model can produce stable, accurate as well as mesh-size independent solutions without a priori assumption of the shear band path.     

Numerical investigation of the slope discontinuities in large deformation finite element formulations

In many multibody system applications, the system components are made of structural elements that can have different orientations, leading to slope discontinuities. In this paper, a numerical investigation of a new procedure that can be used to model structures with slope discontinuities in the finite element absolute nodal coordinate formulation (ANCF) is presented. This procedure can be applied to model slope discontinuities in the case of commutative rotations of gradient deficient elements that are used for modeling thin beam and plate structures. An important special case to which the proposed procedure can be applied is the case of all planar gradient deficient ANCF finite elements. The use of the proposed method leads to a constant orthogonal element transformation that describes an arbitrary initial configuration. As a consequence, one obtains, in the case of large commutative rotations and large deformations, a constant mass matrix for structures which have complex geometry. The procedure used in this investigation to model slope discontinuities requires the use of the concept of the intermediate finite element coordinate system. For each finite element, a new set of gradient coordinates that define, at the discontinuity node, the element deformation with respect to the intermediate element coordinate system is introduced. These new gradient coordinates are assumed to be equal for the two finite elements at the point of intersection. That is, the change of the gradients of two elements at the intersection point from their respective intermediate initial reference configuration is assumed to be the same. This procedure leads to a set of linear algebraic equations that define the orthogonal transformation matrix for the finite element. Numerical examples are presented in order to demonstrate the use of the proposed procedure for modeling slope discontinuities.     

Closed form solutions of Euler–Bernoulli beams with singularities

The problem of the integration of the static governing equations of the uniform Euler–Bernoulli beam with discontinuities is studied. In particular, two types of discontinuities have been considered: flexural stiffness and slope discontinuities. Both the above mentioned discontinuities have been modeled as singularities of the flexural stiffness by means of superimposition of suitable distributions (generalized functions) to a uniform one dimensional field. Closed form solutions of governing differential equation, requiring the knowledge of the boundary conditions only, are proposed, and no continuity conditions are enforced at intermediate cross-sections where discontinuities are shown. The continuity conditions are in fact embedded in the flexural stiffness model and are automatically accounted for by the proposed integration procedure. Finally, the proposed closed form solution for the cases of slope discontinuity is compared with the solution of a beam having an internal hinge with rotational spring reproducing the slope discontinuity.     

Weak discontinuities in a chemically reacting atmosphere

Summary The growth and decay of a weak discontinuity headed by a singular surface of arbitrary shape in three dimensions is investigated in a chemically reacting atmosphere, in the absence of dissipative mechanisms such as viscosity, diffusion and heat conduction. The combined effects of the disequilibrium due to the chemical reaction and a wave front curvature on the propagation of discontinuities have been examined and discussed. It has been observed that the chemical disequilibrium, with its Arrhenius rate dependence, causes the compression wave to steepen more swiftly that it does in an inert atmosphere. The critical values of the initial discontinuity, and time for shock formation, in cases of diverging and converging waves, have been determined.     

The numerical simulations of explosion and implosion in air: use of a modified Harten's TVD scheme

Numerical simulations of explosion and implosion in air are carried out with a modified Harten's TVD scheme. The new scheme has a high resolution for contact discontinuities in addition to maintaining the good features of Harten's TVD scheme. In the numerical experiment of spherical explosion in air, the second shock wave (which does not exist in the one‐dimensional shock tube problem) and its subsequent implosion on the origin have been successfully captured. The positions of the main shock wave, the contact discontinuity and the second shock wave have shown satisfactory agreement with those predicted from previous analysis. The numerical results are also compared with those obtained experimentally. Finally, simulations of a cylindrical explosion and implosion in air are carried out. Results of the cylindrical implosion in air are compared with those of previous work, including the interaction of the reflected main shock wave with the contact discontinuity and the formation of a second shock wave. All these attest to the successful use of the modified Harten's TVD scheme for the simulations of shock waves arising from explosion and implosion. Copyright © 1999 John Wiley & Sons, Ltd.     

Analysis of an interfacial crack in a piezoelectric bi-material via the extended Green’s functions and displacement discontinuity method

Based on the extended Stroh formalism, we first derive the extended Green’s functions for an extended dislocation and displacement discontinuity located at the interface of a piezoelectric bi-material. These include Green’s functions of the extended dislocation, displacement discontinuities within a finite interval and the concentrated displacement discontinuities, all on the interface. The Green’s functions are then applied to obtain the integro-differential equation governing the interfacial crack. To eliminate the oscillating singularities associated with the delta function in the Green’s functions, we represent the delta function in terms of the Gaussian distribution function. In so doing, the integro-differential equation is reduced to a standard integral equation for the interfacial crack problem in piezoelectric bi-material with the extended displacement discontinuities being the unknowns. A simple numerical approach is also proposed to solve the integral equation for the displacement discontinuities, along with the asymptotic expressions of the extended intensity factors and J-integral in terms of the discontinuities near the crack tip. In numerical examples, the effect of the Gaussian parameter on the numerical results is discussed, and the influence of different extended loadings on the interfacial crack behaviors is further investigated.     

Damping numerical oscillations in hybrid solvers through detection of Gibbs phenomenon

A Gibbs phenomenon detector that is useful in damping numerical oscillations in hybrid solvers for compressible turbulence is proposed and tested. It is designed to function in regions away from discontinuities where commonly used discontinuity sensors are ineffective. Using this Gibbs phenomenon detector in addition to a discontinuity sensor for combining central and shock capturing schemes provides an integrated way of dealing with numerical oscillations generated by shock waves and contact lines that are normal to the flow. When complete suppression of numerical oscillations is not possible, they are sufficiently localized. Canonical tests and large eddy simulations show that inclusion of the proposed detector does not cause additional damping of ‘well‐resolved’ physical oscillations. Copyright © 2017 John Wiley & Sons, Ltd.     

Superposition of non-ordinary state-based peridynamics and finite element method for material failure simulations

Superposition of non-ordinary state-based peridynamics and finite element method for material failure simulations, including crack propagation and strain localization is developed. By this approach, a peridynamic model capable of effectively treating strong and weak discontinuities is superimposed in the critical regions over an underlying finite element mesh placed over the entire problem domain. A rigorous variational framework of coupling local finite element and nonlocal peridynamics approximations that is free of blending parameters is developed. Several numerical examples involving mixed-model fracture, three-dimensional adaptive crack propagation and strain localization induced ductile failure demonstrate the rational and efficiency of the proposed superposition-based coupling approach.

     

Function Spaces for Liquid Crystals

We consider the relationship between three continuum liquid crystal theories: Oseen–Frank, Ericksen and Landau–de Gennes. It is known that the function space is an important part of the mathematical model and by considering various function space choices for the order parameters s, n, and Q, we establish connections between the variational formulations of these theories. We use these results to justify a version of the Oseen–Frank theory using special functions of bounded variation. This proposed model can describe both orientable and non-orientable defects. Finally we study a number of frustrated nematic and cholesteric liquid crystal systems and show that the model predicts the existence of point and surface discontinuities in the director.     

Interference of stationary and non-stationary shock waves

The mathematical model of a gasdynamic discontinuity is used in the area of study concerning gas flows with large gradients of gasdynamic functions. Gasdynamic functions before and after the discontinuity meet non-linear algebraic equations called the dynamic compatibility conditions on the discontinuities. Different modes of shock wave structures forming as a result of regular or irregular interference of the incoming discontinuities of different types are described. Ranges of the initial flow parameters definition such that either shock wave structures of different modes take place or interference equations have no solutions are determined. Most attention is given to arbitrary triple shock-wave configurations. Their classification is proposed. Differential characteristics of the steady flow are studied. The notion “differential characteristics” includes first derivatives of the fundamental gasdynamic parameters with respect to natural coordinates and curvatures of the discontinuities surfaces. Effect of unsteadiness on the triple-shock configuration is examined. Some problems arising at creation of complete local theory of steady and propagating gasdynamic discontinuities interference are formulated.     

Some observations of volumetric instabilities in soils

The direct shear apparatus was developed for soil testing because it reproduced the shear failure surfaces that formed a part of the failure mechanism of many geotechnical systems. Typical shear tests on dense sands show a softening load:displacement response which is associated with significant volumetric expansion. While shear localisations can be detected from discontinuities in marker layers, the change in density can be detected using radiography. The progressive formation of shear localisations is readily observed in situations which impose a discontinuity of boundary displacement and these can naturally be interpreted as precursors to a failure mechanism. However, more subtle patterns of volumetric strain or density localisation can be observed in situations where no such obvious boundary displacement discontinuity exists but the sand body is subjected to a more general shearing. Such patterns have a structure which is clearly related to the size of the sand particles. Several examples of such patterns are presented and implications for soil testing and for model tests on soils are discussed.     

Restrictions on dynamically propagating surfaces of strong discontinuity in elastic-plastic solids

For dynamic three-dimensional deformations of elastic-plastic materials, we elicit conditions necessary for the existence of propagating surfaces of strong discontinuity (across which components of stress, strain or material velocity jump). This is accomplished within a small-displacement-gradient formulation of standard weak continuum-mechanical assumptions of momentum conservation and geometrical compatibility, and skeletal constitutive assumptions which permit very general elastic and plastic anisotropy including yield surface vertices and anisotropic hardening. In addition to deriving very explicit restrictions on propagating strong discontinuities in general deformations, we prove that for anti-plane strain and incompressible plane strain deformations, such strong discontinuities can exist only at elastic wave speeds in generally anisotropic elastic-ideally plastic materials unless a material's yield locus in stress space contains a linear segment. The results derived seem essential for correct and complete construction of solutions to dynamic elastic-plastic boundary-value problems.     

Eulerian formulation of transport equations for three-dimensional shock waves in simple elastic solids

The propagation of a three-dimensional shock wave in an elastic solid is studied. The material is assumed to be a simple elastic solid in which the Cauchy stress depends on the deformation gradient only. It is shown that the growth or decay of a discontinuity ψ depends on (i) an unknown quantity φ^? behind the shock wave, (ii) the two principal curvatures of the shock surface, (iii) the gradient on the shock surface of the shock wave speeds and (iv) the inhomogeneous term which depends on the motion ahead of the shock surface and vanishes when the motion ahead of the shock surface is uniform. If a proper choice is made of the propagation vectorb along which the growth or decay of the discontinuity is measured, the dependence on item (iii) can be avoided. However,b assumes different directions depending on the choice of discontinuity ψ with which one is concerned and the unknown quantity φ^? behind the shock wave on which one chooses to depend. As in the case of one-dimensional shock waves, the growth (or decay) of one discontinuity may not be accompanied by the growth (or decay) of other discontinuities. A universal equation relating the growth or decay of discontinuities in the normal stress, normal velocity and specific volume is also presented.     

Broadband elastic metamaterial with single negativity by mimicking lattice systems

The narrow bandwidth is a significant limitation of elastic metamaterials for practical engineering applications. In this paper, a broadband elastic metamaterial with single negativity (negative mass density or Young's modulus) is proposed by mimicking lattice systems. It has two stop bands and the bandwidth of the second one is infinite theoretically. The effect of the relevant parameters on band gaps is discussed. A continuum model is proposed and the selection of materials is discussed in detail. It shows that continuum metamaterials can be described accurately by using the lattice model, and the second stopband can be ultra-broad but not infinite. This discrepancy is investigated and a method is provided to calculate the upper limit of the second stopband for a continuum metamaterial. As a verification, the proposed metamaterial is used for wave mitigation over broadband frequency ranges. Moreover, the present method is extended to design 2D anisotropic elastic metamaterials, and a device to control the direction of elastic wave transmission is proposed as an example.     

Tangential discontinuities of parameters of a polar fluid under shear strain

Possible formation of tangential discontinuities of parameters of a deformable polar fluid is examined by the example of glycerin. It is experimentally established that glycerin under weak shear loads possesses the properties of a non-Newtonian elastoviscoplastic fluid, and formation of tangential discontinuities in viscosity is possible. In the discontinuity region, glycerin has the properties of a low-viscosity fluid, and the structure of the medium is reconstructed after unloading. A rheological equation of the examined fluid is derived, which allows one to analyze the behavior of the medium in different modes of its deformation, including the formation of a local region with reduced viscosity and a tensile stress field. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 3, pp. 41–49, May–June, 2005.