Simulation and interpretation of diffuse mode bifurcation of elastoplastic solids

The diffuse mode bifurcation of elastoplastic solids at finite strain is investigated. The multiplicative decomposition of deformation gradient and the hyperelasto-plastic constitutive relationship are adapted to the numerical bifurcation analysis of the elastoplastic solids. First, bifurcation analyses of rectangular plane strain specimens subjected to uniaxial compression are conducted. The onset of the diffuse mode bifurcations from a homogeneous state is detected; moreover, the post-bifurcation states for these modes are traced to arrive at localization to narrow band zones, which look like shear bands. The occurrence of diffuse mode bifurcation, followed by localization, is advanced as a possible mechanism to create complex deformation and localization patterns, such as shear bands. These computational diffuse modes and localization zones are shown to be in good agreement with the associated experimental ones observed for sand specimens to ensure the validity of this mechanism. Next, the degradation of horizontal sway stiffness of a rectangular specimen due to plane strain uniaxial compression is pointed out as a cause of the bifurcation of the first antisymmetric diffuse mode, which triggers the tilting of the specimen. Last, circular and punching failures of a footing on a foundation are simulated.     

Axisymmetric bifurcation analysis in soils by the tangential-subloading surface model

This article discusses localized bifurcation modes corresponding to shear band formation and diffuse bifurcation modes corresponding to bulge formation for cylindrical soil specimen subjected to an axisymmetric load under undrained conditions. We employ the tangential-subloading surface model, which exhibits the characteristic regimes of the governing equations: elliptic, hyperbolic and parabolic. Also, conditions for shear band formation, shear band inclination, diffuse bulging formation, and the long and short wavelength limits of diffuse bulging modes are discussed in relation to material properties and their state of stress, i.e. the stress ratio and the normal-yield ratio. Tangential-plastic strain rate term is required for the analyses of shear band and diffuse bulging. The shear band and the diffuse bulging are generated in not only normal-yield but also subyield states and they are severely affected by the normal-yield ratio describing the degree of approach to the normal-yield state.     

Non-normality and bifurcation in plane strain tension and compression

The bifurcations of a rectangular block subject to plane strain tension or compression are investigated. The block material is taken to be incompressible and is characterized by an incrementally linear constitutive law for which “normality” does not necessarily hold. The consequences of non-normality regarding bifurcation are given primary emphasis here. The characteristic regimes of the governing equations (elliptic, parabolic and hyperbolic) are detennined. In each of these regimes both symmetric and antisymmetric diffuse bifurcation modes are available. Additionally, in the hyperbolic and parabolic regimes, bifurcation into a localized shear band mode is also possible. Particular attention is given to the limiting cases of long wavelength and soon wavelength diffuse bifurcation modes. The range of parameter values is identified for which bifurcation into some localized mode may precede bifurcation into a long wavelength diffuse mode. Some difficulties associated with employing a linear incremental solid in a bifurcation analysis, when primary interest is in the bifurcation of an underlying elastic-plastic solid, are also discussed.     

Necking of biaxially stretched elastic-plastic circular plates

The necking of an elastic-plastic circular plate under uniform radial tensile loading is investigated both within the framework of the three-dimensional theory and within the context of the plane-stress approximation. Attention is restricted to axisymmetric deformations of the plate. The material behavior is described by two different constitutive laws. One is a finite-strain version of the simplest flow-theory of plasticity and the other is a finite-strain generalization of the simplest deformationtheory, which is employed as a simple model of a solid with a vertex on its yield surface. For an initially uniform plate made of an incompressible material, bifurcation from the uniformly stretched state is studied analytically. The regimes of stress and moduli where the governing axisymmetric three-dimensional equations are elliptic, parabolic or hyperbolic are identified. The plane-stress local-necking mode emerges as the appropriate limiting mode from the bifurcation modes available in the elliptic regime. In the elliptic regime, the main qualitative features of the bifurcation behavior are revealed by the plane-stress analysis, although three-dimensional effects delay the onset of necking somewhat. For the deformation theory employed here, the first bifurcation modes are encountered in the parabolic regime if the hardening-rate is sufficiently high. These bifurcations are not revealed by a plane-stress analysis. For a plate with an initial inhomogeneity, the growth of an imperfection is studied by a perturbation method, by a plane-stress analysis of localized necking, and by numerical computations within the framework of the three-dimensional theory. When bifurcation of the corresponding perfect plate takes place in the elliptic regime, the finite element results show that the plane-stress analysis gives reasonably good agreement with the numerical results. When bifurcation of the corresponding perfect plate first occurs in the parabolic regime, then a bifurcation of the imperfect plate is encountered, that is, the finite element stiffness matrix ceases to be positive definite.     

On evolution of thermo-plastic shear band

The present paper briefly reviews analytical studies of the evolution of thermoplastic shear band, i.e. emergence from uniform deformation, post-instability growth and late stage behaviour. The case studied is the simple shear of temperature and rate-dependent materials with heat transfer. Uniform mode exists before a critical state, if no heat flows out of testpiece. Upon reaching the critical state, bifurcation appears as a result of disturbances, which leads to instability and the formation of narrow shear band. Initially, the band, due to temperature disturbance, can shrink with increasing temperature and strain rate owing to unsteady flow. Then heat conduction dominates and causes the shear band to expand. The postmortem appearance of thermo-plastic shear band manifests itself as balance of plastic work rate and heat diffusion. Melting may also take place within the band.     

Microscopic symmetric bifurcation condition of cellular solids based on a homogenization theory of finite deformation

In this paper, we establish a homogenization framework to analyze the microscopic symmetric bifurcation buckling of cellular solids subjected to macroscopically uniform compression. To this end, describing the principle of virtual work for infinite periodic materials in the updated Lagrangian form, we build a homogenization theory of finite deformation, which satisfies the principle of material objectivity. Then, we state a postulate that at the onset of microscopic symmetric bifurcation, microscopic velocity becomes spontaneous, yet changing the sign of such spontaneous velocity has no influence on the variation in macroscopic states. By applying this postulate to the homogenization theory, we derive the conditions to be satisfied at the onset of microscopic symmetric bifurcation. The resulting conditions are verified by analyzing numerically the in-plane biaxial buckling of an elastic hexagonal honeycomb. It is thus shown that three kinds of experimentally observed buckling modes of honeycombs i.e., uniaxial, biaxial and flower-like modes, are attained and classified as microscopic symmetric bifurcation. It is also shown that the multiplicity of bifurcation gives rise to the complex cell-patterns in the biaxial and flower-like modes.     

Investigation on mechanism of plastic deformation by digital speckle pattern interferometry

An experimental scheme of digital speckle pattern interferometry using two electromagnetic shutters and a three-frame color image board is presented to measure two in-plane components of incremental displacement field of diffuse objects. The deformation of the entire process from initial elastic deformation to fracture can be observed by using a Magneto Optical (MO) disk that has enough memory to record the image data of the entire process. By this system, some tensile experiments are achieved to investigate the mechanism of plastic deformation. The experimental results show that plastic deformation of tensile deformation is nonuniform. The material deforms in the distribution of domains. In the domains, the main deformation is uniform shear and rotation. At the boundary of the domains, the deformation is very large normal strain. With load, both the range and mode of domains are changing. At the end of plastic deformation, only two domains become dominant, and fracture happens at the boundary of the two domains.     

Influence of void nucleation on ductile shear fracture at a free surface

An approximate continuum model of a ductile, porous material is used to study the influence of the nucleation and growth of micro-voids on the formation of shear bands and the occurrence of surface shear fracture in a solid subject to plane strain tension. Bifurcation into diffuse modes is analysed for a plane strain tensile specimen described by these constitutive relations, which account for a considerable plastic dilatancy due to void growth and for the possibility of non-normality of the plastic flow law. In particular, bifurcation into surface wave modes and the possible influence of such modes triggering shear bands is investigated. For solids with initial imperfactions such as a surface undulation, a local material inhomogeneity on an inclusion colony, the inception and growth of plastic flow localization is analysed numerically. Both the formation of void-sheets and the final growth of cracks in the shear bands is described numerically. Some special features of shear band development in the solid obeying non-normality are studied by a simple model problem.     

Bifurcation phenomena in the plane compression test

A basic, compression, bifurcation problem is studied by methods similar to those used by R. Hill and J. W. Hutchinson (1975) for the corresponding tension problem. Bifurcations from a state of homogeneous in-plane compression loading are investigated for a rectangular block of incompressible material constrained to undergo plane deformations. The sides of the block are tractionfree, and it is loaded compressively by a uniform, shear-free, relative displacement of its ends. For a wide class of incrementally-linear time-independent materials only two instantaneous moduli enter into the analysis. Diffuse modes of both symmetric and antisymmetric bifurcation are examined in the elliptic regime of the governing equations, and the possibility of localized modes is considered both inside and outside this regime. Lowest bifurcation stresses are computed for essentially the entire range of possible combinations of material properties and geometry, and these are compared with results obtained by Hill and Hutchinson for the tension problem. The limiting case of large thickness (the semi-infinite block) is considered, confirming the results of M. A. Biot (1965).     

Formation of shear bands in plane sheet

The formation of shear bands in plane sheet is studied, both analytically and experimentally, to enhance the fundamental understanding of this phenomenon and to develop a capability for predicting material failure. The evolution of voids is measured and its interaction with the process of shear banding is examined. The evolving dilatancy in plasticity is shown to have a vital role in analysing the shear-band type of bifurcation, and tremendously reduces the theoretical value of critical stresses. The analyses, referring to both localized and diffuse modes of bifurcation, fairly explain the corresponding observations obtained through testing a dual-phase steel sheet and provide a justification of the constitutive model used.     

Two-mode combination resonances of an in-extensional rotating shaft with large amplitude

In this paper, two-mode combination resonances of a simply supported rotating shaft are investigated. The shaft is modeled as an in-extensional spinning beam with large amplitude. Rotary inertia and gyroscopic effects are included, but shear deformation is neglected. The equations of motion are derived with the aid of the Hamilton principle and then transformed to the complex form. The method of harmonic balance is applied to obtain analytical solutions. Frequency-response curves are plotted for the combination resonances of the first and the second modes. The effects of eccentricity and external damping are investigated on the steady state response of the rotating shaft. The loci of saddle node bifurcation points are plotted as functions of external damping and eccentricity. The results are validated with numerical simulations.     

Buckling of a flush-mounted plate in simple shear flow

The buckling of an elastic plate with arbitrary shape flush-mounted on a rigid wall and deforming under the action of a uniform tangential load due to an overpassing simple shear flow is considered. Working under the auspices of the theory of elastic instability of plates governed by the linear von Kármán equation, an eigenvalue problem is formulated for the buckled state resulting in a fourth-order partial differential equation with position-dependent coefficients parameterized by the Poisson ratio. The governing equation also describes the deformation of a plate clamped around the edges on a vertical wall and buckling under the action of its own weight. Solutions are computed analytically for a circular plate by applying a Fourier series expansion to derive an infinite system of coupled ordinary differential equations and then implementing orthogonal collocation, and numerically for elliptical and rectangular plates by using a finite-element method. The eigenvalues of the resulting generalized algebraic eigenvalue problem are bifurcation points in the solution space, physically representing critical thresholds of the uniform tangential load above which the plate buckles and wrinkles due to the partially compressive developing stresses. The associated eigenfunctions representing possible modes of deformation are illustrated, and the effect of the Poisson ratio and plate shape is discussed.     

Bifurcations in granular media: macro- and micro-mechanics approaches

Various failure modes related to different kinds of bifurcations occur in nonassociated elastoplastic materials such as geomaterials. After presenting experimental evidence, we study this question by means of phenomenological constitutive relations and direct numerical simulations based on the discrete element method. The second-order work criterion related to diffuse failure modes is particularly considered within the framework of continuum and discrete mechanics. The equations of the bifurcation domain boundary and unstable stress direction cones are established. Diffuse failure is simulated numerically by perturbing bifurcation states. To cite this article: F. Darve et al., C. R. Mecanique 335 (2007).     

Analysis of deformation localization based on proposed theory of ultrasonic wave velocity propagating in plastically deformed solids

The localization of plastic deformation is discussed as “stationary discontinuity” characterized by a vanishing velocity of an acceleration wave derived using the author’s proposed theory of ultrasonic wave velocities propagating in plastically deformed solids. To formulate the proposed theory, the elasto-plastic coupling effect was introduced to consider the elastic stiffness degradation due to the plastic deformation. The driving force of the deformation localization is caused by the yield vertex effect, which introduces a pronounced softening of the shear modulus, and geometrical softening due to double slip caused by lattice rotations. In the present paper, it is examined theoretically and experimentally that the diagonal terms of the introduced elasto-plastic coupling tensor represent a slight hardening followed by a pronounced softening of the elastic modulus induced by the point defect development caused by cross slides among dislocations at multiple slip stages similar to the yield vertex effects. The off-diagonal terms represent geometrical softening induced by lattice rotations such as texture evolution. Then, based on the coincidence of the onset strains between localization and acceleration waves of vanishing velocity, the diagrams of diffuse necking, localized necking and forming limit are analyzed by applying the proposed acoustic tensor, which is based on the generalized Christoffel tensor derived by the author, and solving cut off conditions of the quasi-longitudinal wave to determine the onset strains of deformation localization and localization modes. As a result, diagrams of diffuse necking, localized necking and forming limit were obtained. Moreover, the localization modes were determined and distinguished as the SH-mode, SV-mode, tearing mode and splitting mode.     

Analytical solutions of tsunamis generated by underwater earthquakes

Tsunamis induced by underwater earthquakes are theoretically analyzed by applying the linear potential theory. Special attention is placed on the initial state of tsunami. For instantaneous seabed deformations, analytical wave solutions induced by three fundamental seabed deformations at initial stage are derived rather than integral expressions in past studies. These analytical solutions constitute a fundamental base for analyzing waves generated by arbitrary seabed displacement with the help of Fourier analysis. Tsunamis induced by non-instantaneous seabed deformation are analyzed as well. For the sake of examining the contributions of all wave components involved in the tsunami waveform, the amplitude density is proposed to examine the effects of deformation width, water depth, harmonic mode and rising time on waveforms. Results show that a larger ratio of water depth to deformation width results in a greater difference between initial waveform and seabed deformation, and the effect of the rising time is significant in deeper-water configuration. For cosine and sine seabed liftings, the effects of higher harmonic modes might be ignored.     

Dynamics of inelastic deformation of porous rocks and formation of localized compaction zones studied by numerical modeling

The paper presents a numerical analysis of the inelastic deformation process in porous rocks during different stages of its development and under non-equiaxial loading. Although numerous experimental studies have already investigated many aspects of plasticity in porous rocks, numerical modeling gives valuable insight into the dynamics of the process, since experimental methods cannot extract detailed information about the specimen structure during the test and have strong limitations on the number of tests. The numerical simulations have reproduced all different modes of deformation observed in experimental studies: dilatant and compactive shear, compaction without shear, uniform deformation, and deformation with localization. However, the main emphasis is on analysis of the compaction mode of plastic deformation and compaction localization, which is characteristic for many porous rocks and can be observed in other porous materials as well. The study is largely inspired by applications in petroleum industry, i.e. surface subsidence and reservoir compaction caused by extraction of hydrocarbons and decrease of reservoir pressure. Special attention is given to the conditions, evolution, and characteristic patterns of compaction localization, which is often manifested in the form of compaction bands. Results of the study include stress-strain curves, spatial configurations and characteristics of localized zones, analysis of bifurcation of stress paths inside and outside localized zones and analysis of the influence of porous rocks properties on compaction behavior. Among other results are examples of the interplay between compaction and shear modes of deformation.To model the evolution of plastic deformation in porous rocks, a new constitutive model is formulated and implemented, with the emphasis on selection of adequate functions defining evolution of yield surface with deformation. The set of control parameters of the model is kept as short as possible; the parameters are carefully selected to have simple and intuitive physical interpretation whenever possible. Results demonstrate that evolution of the yield surface with deformation has major influence on the resulting characteristics of deformation patterns, which is not sufficiently acknowledged in the literature.     

Buckling of a beam subjected to electromagnetic force and its stabilization by controlling the perturbation of the bifurcation

For a beam subjected to electromagnetic force, magnetoelastic buckling due to the increase of such force is theoretically investigated by taking account of the nonlinearity of the electromagnetic force and the elastic force of the beam. Using Liapunov-Schmidt method and center manifold theory, the equilibrium space, the bifurcation set and the bifurcation diagram are theoretically derived. Also, the effect of the higher modes other than the buckling mode on the mode shape of the postbuckling state is discussed. Furthermore, a control method to stabilize the magnetoelastic buckling is proposed, and the unstable equilibrium state of the beam in the postbuckling state, i.e., the straight position of the beam, is stabilized by controlling the perturbation of the bifurcation.     

A Finite Strain Analysis of Cavity Formation at a Rigid Inhomogeneity

The formation of a cavity by inclusion-matrix interfacial separation is examined by analyzing the response of a plane rigid inclusion embedded in an unbounded incompressible matrix subject to remote equibiaxial dead load traction. A vanishingly thin interfacial cohesive zone, characterized by normal and tangential interface force-separation constitutive relations, is assumed to govern separation behavior. Rotationally symmetric cavity shapes (circles) are shown to be solutions of an interfacial integral equation depending on the strain energy density of the matrix, the interface force constitutive relation and the remote loading. Nonsymmetrical cavity formation, under rotationally symmetric conditions of geometry and loading, is treated within the theory of infinitesimal strain superimposed on a given finite strain state. Rotationally symmetric and nonsymmetric bifurcations are analyzed and detailed results, for the Mooney–Rivlin strain energy density and for an exponential interface force-separation law, are presented. For the nonsymmetric rigid body displacement mode, a simple formula for the critical load is presented. The effect on bifurcation behavior of interfacial shear stiffness and other interface parameters is treated as well. In particular we demonstrate that (i) for the smooth interface nonsymmetric bifurcation always precedes rotationally symmetric bifurcation, (ii) unlike rotationally symmetric bifurcation, there is no threshold value of interface parameter for which nonsymmetric bifurcation will not occur and (iii) interfacial shear may significantly delay the onset of nonsymmetric bifurcation. Also discussed is the range of validity of a nonlinear infinitesimal strain theory previously presented by the author (Levy [1]). This revised version was published online in July 2006 with corrections to the Cover Date.     

Bifurcation phenomena in the rigid inclusion power law matrix composite sphere

Bifurcation of interface separation related to cavity nucleation is analyzed for a radially loaded composite sphere consisting of a rigid inclusion separated from a power law matrix by a uniform, non-linear cohesive zone. Equations for the spherically symmetric and non-symmetric problems are obtained from a hyperelastic finite strain theory by a limiting process that preserves non-linear matrix and interface response at infinitesimal strain. A complete solution to the symmetric problem is presented including bifurcation load, stresses, and evolution of elasto-plastic boundary and interface separation. An analysis of non-symmetric bifurcation, under symmetric conditions of geometry and loading, yields the bifurcation load and first non-symmetric mode shape associated with rigid inclusion displacement. An energy analysis is carried out for both symmetric and non-symmetric problems in order to assess stability of spherically symmetric states to spherically symmetric and non-symmetric “rigid body mode” perturbations.Results are provided for an interface force law that captures interface failure in normal mode and linear response in shear mode. For the symmetric problem, (i) there are threshold parameter values above which bifurcation will generally not occur, (ii) threshold values below which there do not exist equilibria in the post bifurcation regime, (iii) bifurcation occurs after attainment of the maximum interface strength. For the non-symmetric problem, (i) bifurcation always occurs, although it can be delayed by interfacial shear, (ii) for the smooth interface, non-symmetric bifurcation occurs after attainment of the maximum interface strength and always precedes symmetric bifurcation.     

Stability and bifurcation of doubly curved shallow panels under quasi-static uniform load

The focus of this paper is on the investigation of the mathematical nature of buckling from the point of view of bifurcation theory. For the doubly curved orthotropic panels subjected to quasi-static uniform load and with hinged boundary conditions, the solution to the non-linear partial differential equation is partitioned into two parts and projected onto the complete space spanned by the eigenfunctions of the linear operator of the governing equation. Furthermore, the fundamental branch, from which a new solution will emanate, is approximated by the first single mode pair which is close to the real membrane state. Whereas the ensuing bifurcated branch is approximated by the other single mode pair, under the assumption that the coupling between modes can be neglected. The present analysis could give a deep insight into the mechanism of the instability of panel structures, and show that there exists a mode transition at the critical point and the snap-through, then results from saddle-node bifurcation on the bifurcated branch. As a conclusion, the buckling of the system studied can be stated as: a bifurcated branch emanates from the fundamental branch at a critical point, and a saddle-node bifurcation, behaving as jumping, then occurs on the ensuing bifurcated branch.