Strain localisation modelling and pore pressure in saturated sand samples

Dynamic strain localisation theory together with a multiphase material model is used to simulate shear band dominated processes in fully saturated sand samples. The fluid-saturated medium is viewed as multi-phase continuum consisting of a solid skeleton and pores filled by fluids. The governing equations are based on the general framework of averaging theories. A generalised plasticity constitutive model for fully saturated soils is adopted in the computational process. Both samples of medium-loose and dense sands are studied. Negative water pressures, which are important in localisation phenomena of fully saturated dilatant geomaterials, are obtained for dense sands, while positive water pressures result for medium-loose sands.     

饱和砂土基于相变状态的不排水本构模型

饱和砂土基于相变状态的改进本构模型能较好地描述砂土不排水情况下的应力-应变关系。通过改进模型中相变状态参量的计算方法,将e-lgp°表达式加以改进,通过详细的阐述,提出e_pt-lgp的关系式,并通过三轴压缩试验的数据对提出的关系式加以验证。应用新的关系式计算不排水条件下试样的相变孔隙比e_pt,进而计算出模型中的相变状态参量。改进后的不排水本构模型的计算结果与不排水三轴压缩的试验结果能较好地吻合,即验证了本构模型的有效性,也验证了e_pt-lgp关系式的合理性,使原有的边界面本构模型更好地描述不排水的情况下、不同密度和围压状态下饱和砂土的力学行为。     

基于SMP破坏准则的土体弹塑性动力本构模型

利用土体的塑性流动理论,提出了用于描述饱和砂土动力反应性质的弹塑性本构模型。土体总的变形由三部分组成:即弹性应变、与体积屈服机制相关的塑性应变和与剪切屈服机制相关的塑性应变。土体在初始加载与卸载和重新加载阶段性质的差别通过采用不同的模型参数加以反映。该模型能够较为准确地描述饱和砂土在单调加载和循环加载条件下的反应性质。     

Numerical differentiation for use in integrating unified constitutive equations

The NONSS method for numerically integrating time-dependent constitutive equations requires evaluation of the partial derivatives of the constitutive equations. A numerical method for evaluating these partial derivatives has been developed, thus avoiding the cumbersome task of analytically re-deriving the partial derivatives when improvements to the constitutive model are evaluated. This method is based upon the adaptive finite-difference algorithm published by Stepleman and Winarsky. Modifications of their algorithm are described that allow its efficient use within an extended version of NONSS, called NONSS-ND. Using the MATMOD constitutive equations, the capabilities of NONSS-ND are demonstrated by comparing the results of a wide variety of deformation simulations with those of NONSS, for which the derivatives are computed from analytical solutions. The numerical differentiation within NONSS-ND is shown to be sufficiently accurate for use in integrating time-dependent constitutive equations, but results in a 30–100 per cent increase in computation time. Multi-element structural calculations are therefore found to be uneconomical, but NONSS-ND is well suited for testing alternative constitutive models owing to its accuracy and its general applicability to highly non-linear systems of differential equations.     

A mixture model for the compaction of saturated sand

Free draining water saturated sands and dry sands compact progressively under cyclic shear loading, and the rate of compaction increases as the shear strain amplitude increases, independent of the confining stress magnitude. As the sand grains are relatively incompressible, this compaction is essentially an irreversible porosity decrease induced by rearrangement of the granular structure. We present a constitutive model of differential type for porosity variation which has the minimal ingredients necessary to reflect the observed cyclic loading phenomena, and determine the associated material functions of a particular form by correlation with cyclic loading data. A hypoelastic shear response is also correlated to cyclic loading data to complete an isotropic constitutive model. Pore pressure generation under cyclic loading in undrained conditions is evaluated to illustrate the predictions of the model.     

A new perspective on variational methods for stability analysis of columns

Summary The classical linear stability equations for columns take the form of equilibrium equations in terms of displacements. From a variational point of view these equations emerge as conditions of extremum from an energy or virtual work functionals. In principle it should be possible to arrive at the governing equations for stability of columns from a complementary energy or complementary virtual work functional involving force quantities alone. This is rarely done. Both the energy and the complementary energy approaches depend upon the relevant constitutive equations. In this study a formulation is employed wherein focus is maintained on the constitutive equations. It is shown that under certain admissibility conditions a least squares functional, forcing the satisfaction of the constitutive equations, yields the energy and the so-called pure complementary energy functionals as subsets.By way of illustration, a number of examples of conservative and non conservative column buckling problems are analysed by the procedures outlined.With 3 Figures     

A generalized midpoint algorithm for the integration of a generalized plasticity model for sands

This paper describes a method of performing the integration of generalized plasticity models, in which, unlike classical elastoplasticity, the yield surface is not explicitly defined. The algorithm is based on a generalized midpoint scheme and is applied to a specific generalized plasticity model for sands, in which a hyperelastic formulation is introduced to describe the reversible component of the soil response instead of the hypoelastic approach originally proposed. In this way, an efficient integration scheme is developed in the elastic strain space. The consistent, algorithmic tangent operator is derived. Isoerror maps are generated to study the local accuracy of the numerical integration algorithm. Results from a series of numerical examples based on the simulation of drained triaxial tests are given to illustrate the accuracy and convergence properties of the algorithm, both at the local and at the global level. Finally an example is given of the simulation of a cyclic triaxial test to illustrate the improvement on accuracy caused by the use of a hyperelastic law into the constitutive equations, as opposed to the hypoelastic formulation initially adopted in the model. Copyright © 2008 John Wiley & Sons, Ltd.     

Higher-order fractional constitutive equations of viscoelastic materials involving three different parameters and their relaxation and creep functions

Two higher-order fractional viscoelastic material models consisting of the fractional Voigt model (FVM) and the fractional Maxwell model (FMM) are considered. Their higher-order fractional constitutive equations are derived due to the models’ constructions. We call them the higher-order fractional constitutive equations because they contain three different fractional parameters and the maximum order of equations is more than one. The relaxation and creep functions of the higher-order fractional constitutive equations are obtained by Laplace transform method. As particular cases, the analytical solutions of standard (integer-order) quadratic constitutive equations are contained. The generalized Mittag–Leffler function and H-Fox function play an important role in the solutions of the higher-order fractional constitutive equations. Finally, experimental data of human cranial bone are used to fit with the models given by this paper. The fitting plots show that the models given in the paper are efficient in describing the property of viscoelastic materials.     

各向同性弹性介质非线性本构方程

从张量函数出发,围绕共轭应力、应变变量,研究了各向同性非线性弹性介质各种形式的本构方程以及各种形式方程之间的关系。推导出用张量不变量,标量不变量表示的两种形式非线性Green弹性介质本构方程。证明了方程是完备的,不可约的。作为应用举例,研究了橡胶材料的工程应用问题。     

Assessment of plate theories for multilayered angle-ply plates

The accuracy of laminated plate theories for the analysis of angle-ply multilayered plates is investigated. Classical Lamination Theory, First Order and Reddy's High Order Shear Deformation Theories are considered, together with a new two-dimensional theory recently proposed by the second author (2D-theory). Deflections and through-the-thickness stress distributions are compared with the exact 3D elasticity solution. Energy errors of stress fields obtained from the various theories through both constitutive and equilibrium equations are also evaluated. It is shown that equilibrium equations are required to obtain sufficiently accurate stress distributions from CLT, FSDT and HSDT, whereas stresses can be obtained directly from constitutive equations for the 2D-theory.     

Elasto-viscoplastic constitutive equations for rock-type materials (finite deformation)

In the present paper the elasto-viscoplastic frame indifferent constitutive equations for rock and rock-type materials are proposed within the frame-work of large deformations. The constitutive equations are written with respect to an arbitrary configuration at time t as reference configurations. Our model describe the dilatancy or compressibility and creep for rock-type materials as well as the existence of the irreversible part of strain even when the applied stresses are relatively small. Constitutive hypotheses are based on the experimental evidences which reveal the complex behaviour or rock-type materials. All the constitutive functions and moduli involved in the model can be determined from a complete set of experimental data.     

大块非晶合金复杂应力状态塑性本构方程

为得到复杂应力状态下大块非晶合金的塑性本构方程,对Zr41.2Ti13.8Ni10Cu12.5Be22.5大块非晶合金进行简单拉伸、压缩和扭转实验,得到相应的一维本构关系,并应用弹塑性力学的相关理论,将一维本构方程延拓到多维应力空间.研究结果表明:该材料符合Mises屈服条件并具有强化特性;利用简单压缩和扭转实验结果可推导出材料在复杂应力状态下的塑性本构方程,且这两个塑性本构方程是一致的.     

Applications of tensor functions to the formulation of constitutive equations involving damage and initial anisotropy

In this paper some results of the tensor function theory are applied to the formulation of constitutive equations of isotropic and anisotropic materials in the secondary and tertiary creep stage. The creep process, in its tertiary phase, is characterized by a damage tensor. Because of its microscopic nature, damage has, in general, an anisotropic character even in cases where the material was originally isotropic, i.e. isotropic in its virgin state. Fissure orientation and length cause anisotropic macroscopic behaviour. In the first part of the paper some possible ways of representing constitutive equations involving (initial) anisotropy of the material (e.g. from rolling) and involving anisotropic creep-damage are dealt with. The formulations of such equations are based upon theorems concerning tensor-valued functions. Furthermore, some simplified constitutive equations for more practical use are discussed. The main problem of this part is: to find an irreducible set of tensor generators. Besides the problem of finding such tensor generators it is very important to determine the scalar coefficients in constitutive equations as functions of the invariants and experimental data. The second part of the paper is concerned with the determination of the scalar functions. This can be done by using tensorial interpolation methods as pointed out in detail.     

Thermal aspects of shear localization in microporous viscoplastic solids

A set of phenomenological state variable constitutive equations for large deformations of isotropic, thermo-elasto-viscoplastic, porous materials is considered, and a new semi-implicit, incrementally objective time-integration procedure for these constitutive equations is developed. The constitutive equations and the time-integration procedure are implemented in a finite element program, and used to study the shear band formation and ductile fracture initiation triggered by flow softening due to the combined effects of void growth and deformation heating, in a plane strain tension test.     

Generalised constitutive equations for densification of metal matrix coated fibre composites

Abstract

The present paper completes a study of constitutive equations for the consolidation processing of continuous fibre reinforced metal matrix composite materials. It builds on an earlier paper in which physically based constitutive equations were derived for the case of symmetrical, isostatic loading. In the present paper, constitutive equations are developed for in plane, general stress states. The total deformation of the consolidating composite is expressed as the sum of a conventional deviatoric creep term, together with a dilatational term, which was derived using a variational method previously published. The equations contain only two material parameters, which are the conventional creep coefficient and exponent for the fibre coating material (in this case, Ti-6Al-4V). The resulting equations have been implemented into finite element software enabling the simulation of practical consolidation processes. The model has been verified by comparing predicted results with those obtained from independent micromechanical models. A number of experimental tests have been carried out, and the model is used to predict the rates of densification for a range of experimental pressure and temperature histories. Good comparisons have been achieved.     

含损伤复合材料的湿热弹性响应

依据不可逆热力学理论, 未引入任何附加假设, 建立了湿热弹性各向异性损伤复合材料的一般理论。应用建立损伤本构方程的本构泛函展开法, 推导出湿热弹性损伤材料全部本构方程的一般形式, 其中包括比自由能密度表达式、 应力-应变关系、 熵密度方程、 损伤应变能释放率表达式、 吸湿对偶力表达式、 湿-热-固-损伤耦合的热传导方程和损伤演化方程。研究表明, 在本构方程中含有若干损伤效应函数, 表征损伤对材料宏观力学性能与湿、 热性能的影响, 其具体形式可由细观力学解确定, 从而使连续损伤力学与细观损伤力学有机结合在一起。最后, 从细观力学与实验观测两个角度, 举例说明损伤效应函数与系数张量的确定方法, 为分析变温变湿环境下复合材料的损伤问题提供重要的理论依据。      

Granular element method (GEM): linking inter-particle forces with macroscopic loading

We present a new method capable of inferring, for the first time, inter-particle contact forces in irregularly-shaped natural granular materials (e.g., sands), using basic Newtonian mechanics and balance of linear momentum at the particle level. The method furnishes a relationship between inter-particle forces and corresponding average particle stresses, which can be inferred, for instance, from measurements of average particle strains emanating from advanced experimental techniques (e.g., 3D X-ray diffraction). Inter-particle forces are the missing link in understanding how forces are transmitted in complex granular structures and the key to developing physics-based constitutive models. We present two numerical examples to verify the method and showcase its promise.     

Hyper-elastic constitutive equations of conjugate stresses and strain tensors for the Seth-Hill strain measures

The use of hypo-elastic constitutive equations for large strains in nonlinear finite element applications usually requires special considerations. For example, the strain does not tend to zero upon unloading in some elastic loading-unloading closed cycles. Furthermore, these equations are based on objective material time rate tensors, which require incrementally objective algorithms for numerical applications and integration. Hyper-elastic constitutive equations on the other hand do not require such considerations. However, their behaviour for large elastic strains is important and may differ in tension and compression. In the present work, Hyper-elastic constitutive equations for the Seth-Hill strains and their conjugate stresses are explored as a natural generalisation of Hook’s law for finite elastic deformations. Based on the uniaxial and simple shear tests, the response of the material for different constitutive equations is examined. Together with an objective rate model, the effect of different constitutive laws on Cauchy stress components is compared. It is shown that the constitutive equation based on logarithmic strain and its conjugate stress gives results closer to that of the rate model. In addition, the use of Biot stress-strain pairs for a bar element results in an elastic spring which obeys the Hook’s law even for large deformations and has the same behaviour in both tension and compression. The effect of the constitutive equation on the volume change of the material has also been considered here.     

An approach to finite static axially symmetric deformation of a hyperelastic membrane

Summary A simple approach to the problem of finite quasi-static axially symmetric deformation of an isotropic incompressible hyperelastic membrane is presented. Lagrangian type equilibrium equations, expressed in terms of the Biot stresses, are used along with constitutive relations expressing the principal components of Biot stress in terms of the principal stretches. Numerical results, obtained from the application of a finite element method to the governing equilibrium equations and constitutive relations, are presented for two problems, and for two different strain energy functions. It is shown how the proposed equilibrium equations can be obtained from a variational principle, and the variational principle is also used to obtain approximate solutions.     

Acceleration waves in simple and linear viscoelastic micropolar materials

In this work the conditions of propagation of acceleration waves are studied in different types of micropolar media: (a) in simple media whose constitutive equations are given in functional form, (b) in linear viscoelastic media of which the constitutive equations are obtained in the limit of small deformations from a simple micropolar medium obeying the axiom of fading memory. In contrast to the first case, two kinds of acceleration waves can propagate independently in the second case.