Solution existence conditions for elastoplastic constitutive models of granular materials

In this paper, the conditions of solution existence for stress rates under given strain rates are investigated. The focus of the solution existence investigation is on the non-associated flow rule and elastic stress–strain relationship. Granular materials characterized with strong non-associated plastic flows are used as a particular example for analysis. Various flow rules for granular materials are analyzed, including Rowe’s, Roscoe’s flow rules and their modified versions. In the elastic stress–strain relationships of materials, the effects of Poisson’s ratio on solution existence are investigated. Both isotropic and anisotropic elasticity are considered. Given a granular material and its states, it is found that there exists a critical Poisson’s ratio for a particular non-associated flow rule. When the Poisson’s ratio of a material is above this critical Poisson’s ratio, its constitutive model is susceptible to solution non-existence. It is suggested that special attentions should be paid to the selection of material Poisson’s ratio and non-associated flow rule to ensure the existence of elastoplastic solutions.     

Thermo-Chemo-Electro-Mechanical Formulation of Saturated Charged Porous Solids

A thermo-chemo-electro-mechanical formulation of quasi-static finite deformation of swelling incompressible porous media is derived from a mixture theory including the volume fraction concept. The model consists of an electrically charged porous solid saturated with an ionic solution. Incompressible deformation is assumed. The mixture as a whole is assumed locally electroneutral. Different constituents following different kinematic paths are defined: solid, fluid, anions, cations and neutral solutes. Balance laws are derived for each constituent and for the mixture as a whole. A Lagrangian form of the second law of thermodynamics for incompressible porous media is used to derive the constitutive restrictions of the medium. The material properties are shown to be contained in one strain energy function and a matrix of frictional tensors. A principle of reversibility results from the constitutive restrictions. Existing theories of swelling media should be evaluated with respect to this principle.     

A systematic procedure for constructing critical state models in three dimensions

A general procedure for developing constitutive models for frictional materials possessing a critical state is developed in a three-dimensional context. The procedure starts from the laws of thermo-dynamics, so that the first and second laws of thermo-dynamics are automatically satisfied. There is hence no need to invoke any extraneous stability postulates. The models involve a number of parameters, which can be interpreted in terms of micro-mechanical energy storage and dissipative mechanisms. In most cases non-associated flow rules are predicted and in some cases the yield surfaces are seen to have concave segments. The procedure is more general than that traditionally used for materials with non-associated flow rules, in that plastic potentials are not needed and not presumed to exist. In illustration, examples of families of models are given in which the critical state surface is either the Drucker–Prager or the Matsuoka–Nakai cone.     

A generalized continuum theory for granular materials

A generalized continuum theory for granular media is formulated by allowing for the possibility of rotation of granules. The basic balance laws are presented and based on thermodynamical consideration a set of constitutive equations are derived. The theory naturally gives rise to the generation of antisymmetric stress tensor and existence of couple stresses. The basic equations of motion are derived and it is shown that the theory contains Mohr-Coulomb criterion of limiting equilibrium as a special case. The problem of coupled porosity and microrotational wave propagation is investigated and the rectilinear shear flow of granular materials is discussed.     

Ideal nonsymmetric 4D-medium as a model of invertible dynamic thermoelasticity

The space-time continuum (4D-medium) is considered, and a generalized model of reversible dynamic thermoelasticity is constructed as a theory of elasticity of an ideal (defect-free) nonsymmetric 4D-medium that is transversally-isotropic with respect to the time coordinate. The definitions of stresses and strains for the space-time continuum are introduced. The constitutive equations of the medium model relating the components of nonsymmetric stress and distortion 4D-tensors are stated. Physical interpretations of all tensor components of the thermomechanical properties are given. The Lagrangian of the generalized model of coupled dynamic thermoelasticity is presented, and the Euler equations are analyzed. It is shown that the three Euler equations are generalized equations of motion of the dynamic classical thermoelasticity, and the last, fourth, equation is a generalized heat equation which allows one to predict the wave properties of heat. An energy-consistent version of thermoelasticity is constructed where the Duhamel-Neumann and Maxwell-Cattaneo laws (a nonclassical generalization of the Fourier law for the heat flow) are direct consequences of the constitutive equations.     

Quasi-static and dynamical analyses of a thermoviscoelastic Timoshenko beam using the differential quadrature method

The quasi-static and dynamic responses of a thermoviscoelastic Timoshenko beam subject to thermal loads are analyzed. First, based on the small geometric deformation assumption and Boltzmann constitutive relation, the governing equations for the beam are presented. Second, an extended differential quadrature method(DQM)in the spatial domain and a differential method in the temporal domain are combined to transform the integro-partial-differential governing equations into the ordinary differential equations. Third, the accuracy of the present discrete method is verified by elastic/viscoelastic examples, and the effects of thermal load parameters, material and geometrical parameters on the quasi-static and dynamic responses of the beam are discussed. Numerical results show that the thermal function parameter has a great effect on quasi-static and dynamic responses of the beam. Compared with the thermal relaxation time, the initial vibrational responses of the beam are more sensitive to the mechanical relaxation time of the thermoviscoelastic material.     

The fundamental role of nonlocal and local balance laws of material forces in finite elastoplasticity and damage mechanics

In this paper the fundamental role of independent balance laws of material forces acting on dislocations and microdefects is shown. They enable a thermodynamically consistent formulation of dissipative deformation processes of continua with dislocation motion and defect evolution in the material space on meso- and microlevel.The balance laws of material forces together with the classical balance laws of physical forces and couples, first and second laws of thermodynamics for physical and material space and general constitutive equations are the basis to develop a thermodynamically consistent framework of nonlocal finite elastoplasticity and brittle and ductile damage.It is shown that a weakly-nonlocal formulation of the balance laws of material forces leads to gradient theories, where local theories are obtained, if all gradient contributions are assumed to be small. In this case the local balance laws of material forces together with the constitutive equations represent evolution laws of the material forces. In the classical approach of internal variables they are assumed from the outset with the result that there is a large number of different propositions in the literature.The well-known splitting test of a circular cylinder of concrete is simulated numerically, where the process of deformation in the physical space and defect and plastic evolution in the material space is represented.     

Dynamic strain aging and related instabilities: experimental,theoretical and numerical aspects

Experimental data from uniaxial tensile tests on smooth and notched specimens of aluminium alloy 5083-H116 show that the material exhibits negative strain-rate sensitivity for strain rates within a certain range. The negative strain-rate dependence, which is attributed to dynamic strain aging, leads to serrated stress–strain curves, discontinuous plastic flow and propagating deformation bands during plastic straining (also denoted as the Portevin–Le Chatelier effect). Band analysis and linear perturbation analysis are performed using simple elastic-viscoplastic constitutive equations that include negative strain-rate sensitivity in a simplified manner. The negative strain-rate sensitivity allows for jumps in the plastic strain rate, which in turn permits the existence of localisation bands for the elastic-viscoplastic model. The simple elastic-viscoplastic constitutive model has been implemented in LS-DYNA, and non-linear finite element simulations of smooth and notched tensile test specimens are performed, allowing more detailed investigations into the effects of the negative strain-rate sensitivity on the material's behaviour.     

DYNAMIC ANALYSIS OF A SPATIAL COUPLED TIMOSHENKO ROTATING SHAFT WITH LARGE DISPLACEMENTS

ＩｎｔｒｏｄｕｃｔｉｏｎＲｏｔａｔｉｎｇｓｈａｆｔｓａｒｅｔｈｅｍｏｓｔｖｉｔａｌｃｏｍｐｏｎｅｎｔｓｏｆｍｏｄｅｒｎｉｎｄｕｓｔｒｉａｌａｎｄｐｏｗｅｒｇｅｎｅｒａｔｉｏｎｆａｃｉｌｉｔｉｅｓ.ＤｕｅｔｏｔｈｅｉｍｐｏｒｔａｎｃｅｏｆｔｈｅｓｅｃｏｍｐｏｎｅｎｔｓｔｈｅｒｅｗｅｒｅｗｉｄｅｌｙｓｔｕｄｉｅｓｏｎｔｈｅｖｉｂｒａｔｉｏｎｂｅｈａｖｉｏｒｏｆＥｕｌｅｒ＿Ｂｅｒｎｏｕｌｌｉｒｏｔａｔｉｎｇｓｈａｆｔｓｕｓｉｎｇａｎａｌｙｔｉｃａｌａｎｄｎｕｍｅｒｉｃａｌｍｅｔｈｏｄｓ[1- 4 ].Ｈｏｗｅ…     

Consequences of the second law for a turbulent fluid flow

Second Law statements in thermomechanics applicable to turbulent fluid flow, in which the internal energy in a macroscopic field theory includes contributions both from molecular vibrations and from turbulent fluctuations, are discussed. In the absence of turbulence, these statements naturally reduce to the known and accepted Second Law statements for a nonturbulent medium. The usual version of the Second Law statements — which deny the existence of perpetual motion and place restrictions on the constitutive equations —is extended here in the presence of turbulence; and an additional statement is introduced associated with the tendency of turbulent fluctuations to decay in the absence of external work or the addition of thermal heat. The mathematical representations of various Second Law statements are then used to derive several restrictions on the response variables of the macroscopic turbulence theory. Examples of such variables include the rates of production and dissipation of turbulent fluctuations, the rate of thermal entropy production, internal energy (involving constitutive coefficients which may be taken to be the thermal and turbulent specific heats), turbulent viscosity coefficients and other response functions which control the degree of flow anisotropy in the medium. These Second Law restrictions are then applied to a recent theory of macroscopic turbulent flow by the present authors in which fairly general constitutive equations are presented for the dependent variables of the theory. It is found that not only is the range of values of several constitutive coefficients limited by these Second Law restrictions, but the presence of a number of terms in the constitutive equations is entirely denied.     

A micromechanics-based elastoplastic damage model for granular materials at low confining pressure

This paper is devoted to the formulation of a micromechanics-based constitutive model for granular materials under relatively low confining pressure. The constitutive formulation is performed within the general framework of homogenization for granular materials. However, new rigorous stress localization laws are proposed. Some local constitutive relations are established under the consideration of irreversible thermodynamics. Macroscopic plastic deformation is obtained by considering local plastic sliding in a limit number of families of contact planes. The plastic sliding at each contact plane is described by a non-associated plastic flow rule, taking into account pressure sensitivity and normal dilatancy. Nonlinear elastic deformation related to progressive compaction of contacts is also taken into account. Material softening is described by involving damage process related to degradation of microstructure fabric. The proposed model is applied to some typical granular materials (sands). The numerical predictions are compared with experimental data.     

Mathematical simulation of two-phase fluid flow through a water-flooded porous reservoir with sediment formation

The results of numerical investigations of the problem of flow through a porous reservoir in the process of its water flooding with addition of a sediment-forming component are given. The mathematical model is based on the mass conservation laws for each of the considered phases and components supplemented with the equations of motion and constitutive relations necessary to close the system of equations. In solving the problem, an empirical dependence of the sediment formation intensity on the content of the sediment-forming component in the aqueous solution with allowance for variation in the effective porosity of the medium is used. The main features in solving the sediment-formation problem are distinguished using the empirical dependence and a contrastive analysis of the effect of choosing this dependence on the solution results is carried out. It is shown that the neglect of the experimental results in the mathematical formulation can lead to not only unjustified overestimated results in realizing the method but also give a distorted pattern of the entire process of sediment formation in fluid flow through a water-flooded porous reservoir.     

Simple shear of a strain-hardening elastoplastic hollow circular cylinder

Stress and deformation analysis of the simple shear at finite strain of a strain-hardening elastoplastic hollow circular cylinder is given. Both isotropic and anistotropic hardening models are considered. In the case of isotropic hardening, there is a closed from analytical solution. No normal stresses exist in this case. Purely kinematic hardening with a Mises-type yield condition is utilized as a model of anisotropic hardening. Conventional (average) spin is taken to construct the objective Jaumann derivative needed in the structure of the corresponding constitutive laws. Governing partial differential equations are derived and solved numerically to give stress and deformation distribution following the advance of plastic flow. The extent or range of the appropriateness of the considered constitutive model is also established.     

HOMOTOPY PERTURBATION SOLUTION AND PERIODICITY ANALYSIS OF NONLINEAR VIBRATION OF THIN RECTANGULAR FUNCTIONALLY GRADED PLATES

In this paper nonlinear analysis of a thin rectangular functionally graded piate is formulated in terms of von-Karman＇s dynamic equations. Functionaily Graded Material （FGM） properties vary through the constant thickness of the plate at ambient temperature. By expansion of the solution as a series of mode functions, we reduce the governing equations of motion to a Duffing＇s equation. The homotopy perturbation solution of generated Duffing＇s equation is also obtained and compared with numerical solutions. The sufficient conditions for the existence of periodic oscillatory behavior of the plate are established by using Green＇s function and Schauder＇s fixed point theorem.     

On the simultaneous elongation and inflation of a tubular membrane of BKZ fluid

This work considers a viscoelastic fluid membrane which is initially tubular and bonded at each end to a rigid circular disc. The membrane is subjected to prescribed elongational and internal pressure histories causing it to undergo quasi-static axisymmetric deformation. This example is intended to simulate an experiment which has been recently proposed for the determination of constitutive properties for viscoelastic fluids as well as some polymer sheet forming process.The constitutive equation is presumed to be of integral type. The formulation of the problem leads to a basic system of equations which is intended for numerical solution. It has the structure of a two-point boundary value problem for a system ordinary differential equations at each time. The formulation has the advantage that the equations do not have to be rederived if the constitutive equation is changed. A change in the sub-program for computing stress from stretch history is all that is needed.A numerical method of solution is presented. In a numerical example, the material is taken to be polyisobutylene, modeled as a BKZ fluid.     

The existence of solutions in nonlinear elastodynamics

I.Intr0ductionIntheclassicalelasticity,thelinearconstitutivelawisassumed,wherethehyper-potentialisapositiveorsemi-positivedefinitequadraticform.Forsuchsystems,theexistence,uniquenessandrelativepr0pertiesofsolutionstostaticanddynamicproblemscanbefound[l'2'…     

Second-order bifurcation in elastic-plastic solidS

Second-order rate constitutive equations are formulated for a time-independent elastic-plastic material, obeying the normality flow rule with a smooth yield surface. Under specified regularity restrictions imposed on the involved fields, the regular second-order rate boundary value problem with quasistatic accelerations as unknowns is posed. It is shown that every solution of this generally non-linear rate problem is governed by a variational principle and that the corresponding functional reaches a strict absolute minimum, provided the solution satisfies a sufficient uniqueness condition. With the same incrementally linear comparison solid, Hill's exclusion condition rules out not only a first- but also a second-order bifurcation. The criticality of the exclusion condition is discussed and conditions are indicated under which a second-order bifurcation becomes possible, while the first-order rate problem is still uniquely solvable.     

Three-dimensional limit analysis of rigid blocks assemblages. Part II: Load-path following solution procedure and validation

A novel solution procedure for the non-associated limit analysis of rigid blocks assemblages is proposed. This proposal produces better solutions than previously proposed procedures and it is also able to provide an insight into the structural behaviour prior to failure. The limit analysis model proposed in Part I of this paper and the solution procedure are validated through illustrative examples in three-dimensional masonry piers and walls. The use of limit analysis for three-dimensional problems incorporating non-associated flow rules and a coupled yield surface is novel in the literature.     

Blowup of solutions for the planar motions of rotating nonlinearly elastic rods

This paper treats the motion of flexible, extensible, shearable nonlinearly elastic rods, described by a geometrically exact theory, when they are confined to a plane rotating about a fixed axis at constant angular speed and when they are confined to a fixed plane with one end rotating at a constant angular speed about an axis perpendicular to the fixed plane. The paper gives restrictions on the constitutive equations and initial conditions that ensure that motions become unbounded at rapid rates as time becomes infinite. The analysis of these constitutive restrictions employs the theory of characteristics for single first-order semilinear partial differential equations.     

Three-dimensional limit analysis of rigid blocks assemblages. Part I: Torsion failure on frictional interfaces and limit analysis formulation

The governing equations of plastic torsion of, arbitrarily shaped, frictional interfaces are established and a general solution is proposed. The case of rectangular interfaces is used to study the interactions of torsion strength with bending moments and shear forces. In order to solve limit analysis problems, a piecewise linear approximation of the yield function for rectangular interfaces is proposed. A model for the limit analysis of three-dimensional block assemblages interacting through no-tension frictional interfaces is presented including the proposal for the torsion failure mode. This model takes into account non-associated flow rules and limited compressive stresses at the interfaces.