Perturbation growth and localization in fluid-saturated inelastic porous media under quasi-static loadings

Instabilities in inelastic saturated porous media are investigated here for general three-dimensional states under quasi-static loadings using a perturbation approach and focussing in particular on the two limiting cases of the onset of growth and of blowing up of perturbations.For associative flow rules for the skeleton, both onset of growth and blowing up of perturbations depend only on the underlying drained properties. Unbounded growth is obtained when the condition of localization for the underlying drained deformation (singularity of the drained acoustic tensor) is approached or just passed. Onset of growth has always a divergence growth character and critical conditions are always associated to the shortwavelength regime leading to the fact that the failure mode is expected to be a localized one.For non-associative behaviour of the skeleton we show in contrast that the onset of growth and unbounded growth may be defined either by the drained or undrained properties. One or the other depends on the details of the constitutive behaviour but also on the type of loadings. In particular, unbounded growth occurs when either the condition of localization under drained or undrained conditions is first passed. Transition from decaying to growing behaviour may have a divergence character or flutter-type character. Here the critical conditions are associated either to the shortwavelength or to the longwavelength regimes and therefore the failure mode may be localized or diffuse.The hierarchy between criticality of drained and undrained properties is analysed for a general class of constitutive equations and the results are fully and explicitly illustrated for saturated porous media with skeleton obeying Drucker-Prager like constitutive model.     

3-D DYNAMIC RESPONSE OF TRANSVERSELY ISOTROPIC SATURATED SOILS

Introduction Thedynamicanalysisofsaturatedporoelasticmediaisofbroadmeaninginmanyaspectssuchasseismology,earthquakeengineering,soilmechanicsandgeophysics,etc.Thetransverselyisotropicsaturatedsoilsissolid fluidcoupledtwophasemediainwhichthe soilskeletondisp…     

Basic formulation of static and dynamic behaviours of soil and other porous media

This paper treats the soil and rock or concrete as a two-phase medium composed of a solid skeleton and an interstitial fluid. Under the unsaturated condition of interstitial fluid, the third phase is taken into consideration. In this paper, we consider the constitutive relation, dynamic and kinematic relations, and also various approximations with their limits of validity, including those of most practical engineering situations such as the consolidation problem and the undrained behaviour. The finite element discretization and the time dependent solution of various classes of soil problem are also discussed.     

Two-dimensional stress wave analysis in incompressible elastic solids

Two-dimensional stress waves in a general incompressible elastic solid are investigated. First, basic equations for simple waves and shock waves are presented for a general strain energy function. Then the characteristic wave speeds and the associated characteristic vectors are deduced. It is shown that there usually exist two simple waves and two shock waves. Finally, two examples are given for the case of plane strain deformation and antiplane strain deformation, respectively. It is proved that, in the case of plane strain deformation the oblique reflection problem of a plane shock is not solvable in general.     

Discontinuous bifurcation analysis of thermodynamically consistent gradient poroplastic materials

In this work, analytical and numerical solutions of the condition for discontinuous bifurcation of thermodynamically consistent gradient-based poroplastic materials are obtained and evaluated. The main aim is the analysis of the potentials for localized failure modes in the form of discontinuous bifurcation in partially saturated gradient-based poroplastic materials as well as the dependence of these potentials on the current hydraulic and stress conditions. Also the main differences with the localization conditions of the related local theory for poroplastic materials are evaluated to perfectly understand the regularization capabilities of the non-local gradient-based one. Firstly, the condition for discontinuous bifurcation is formulated from wave propagation analyses in poroplastic media. The material formulation employed in this work for the spectral properties evaluation of the discontinuous bifurcation condition is the thermodynamically consistent, gradient-based modified Cam Clay model for partially saturated porous media previously proposed by the authors. The main and novel feature of this constitutive theory is the inclusion of a gradient internal length of the porous phase which, together with the characteristic length of the solid skeleton, comprehensively defined the non-local characteristics of the represented porous material. After presenting the fundamental equations of the thermodynamically consistent gradient based poroplastic constitutive model, the analytical expressions of the critical hardening/softening modulus for discontinuous bifurcation under both drained and undrained conditions are obtained. As a particular case, the related local constitutive model is also evaluated from the discontinuous bifurcation condition stand point. Then, the localization analysis of the thermodynamically consistent non-local and local poroplastic Cam Clay theories is performed. The results demonstrate, on the one hand and related to the local poroplastic material, the decisive role of the pore pressure and of the volumetric non-associativity degree on the location of the transition point between ductile and brittle failure regimes in the stress space. On the other hand, the results demonstrate as well the regularization capabilities of the non-local gradient-based poroplastic theory, with exception of a particular stress condition which is also evaluated in this work. Finally, it is also shown that, due to dependence of the characteristic lengths for the pore and skeleton phases on the hydraulic and stress conditions, the non-local theory is able to reproduce the strong reduction of failure diffusion that takes place under both, low confinement and low pore pressure of partially saturated porous materials, without loosing, however, the ellipticity of the related differential equations.     

An efficient method of analysis for the drained and undrained behaviour of an elastic soil

In this paper a method of analysis for the drained and undrained behaviour of a saturated soil with an elastic anisotropic skeleton is proposed. The method depends upon embedding the problem in a series of problems depending upon a parameter, then expanding the solution in terms of that parameter. The method has the advantage that both the drained and undrained behaviour can be found simultaneously and the problems of a singular stiffness matrix usually associated with undrained behaviour of an elastic soil is avoided.     

饱和土中球空腔的瞬态动力响应

采用工程上通用的饱和土力学模型,考虑流体与固全之间的耦合作用,利用Ｌａｐｌａｃｅ变换求解了饱和土中球空腔的瞬态动力响应问题,得到了变换域内的解析解,借助数值Ｌａｐｌａｃｅ换求解了饱和土中球空腔瞬态动力响应的位移、应力及孔压的变化规律,为分析地下结构瞬态动力响应提供了一种有效的方法,模型符合工程实际,有一定的工程应用价值。     

川藏公路地质环境与整治改建方案的思考

川藏公路由于地质环境复杂、建设标准低、后遗病害多,抗灾能力差,泥石流、滑坡、山崩、雪害、水毁等自然灾害频繁发生,公路阻车断道严重。国家投入巨资进行整治改建,并取得了明显的效果,但由于自然环境特殊、影响因素复杂,许多特大型、大型工程地质病害问题还没有可行、可靠的解决方案。本文通过分析川藏公路沿线的地质环境和灾害特点,总结历年整治改建和经验的教训,提出川藏公路建设的途径、可能达到的目标和应采用的原则。     

Modelling of two-phase transient flow and combustion of granular propellants

Gas-permeable solid propellants possess great potential for producing high thrusts during extremely short time intervals. The paper presents a theoretical model describing the important physical phenomena taking place in both gaseous and solid phases in such two-phase, unsteady, reactive-flow systems. The governing equations were derived in the form of coupled, non-linear partial differential equations and the resulting formulation was compared with those used by previous authors. A stable, fast-converging, fully-implicit numerical method, first proposed by Spalding and implemented in a general-purpose computer program, was used to solve the equations, so that effects of importance can be studied. Solutions are presented that predict the pressure wave build-up and accelerating flame front for beds of granulated solid propellants fixed in a rigid enclosure, and in gun barrels with an accelerating projectile. The results show that the flame front accelerates and the rate of pressurization increases substantially in the downstream direction. Some discussion is given of the sensitivity of these predictions to the assumed constitutive relations for interphase heat and momentum transfer and the solid burning-rate law. It is argued that the controversy surrounding the hyperbolicity of the equations, which has delayed progress towards theoretical prediction of two-phase flows, is ill-founded, since the equations are shown to possess unique solutions, as indeed do the physical systems they represent.     

成层饱和介质平面波斜入射问题的一维化时域方法

地震波斜入射下自由场的输入是大型结构抗震分析中亟待解决的问题之一,尤其是成层饱和多孔介质自由场问题,由于问题的复杂性,目前研究甚少. 本文基于Biot提出的饱和多孔介质动力方程,建立了一种新的求解平面波斜入射下基岩上覆饱和多孔介质成层场地自由场分析的一维化时域计算方法. 该方法首先根据Snell定律将饱和多孔介质二维空间问题转化为一维时域问题,通过对深度方向的有限元离散,得到饱和多孔介质波动问题的一维化有限元方程,然后采用单相弹性介质精确人工边界条件模拟基岩半空间的波动辐射和输入特征,通过考虑基岩与饱和多孔介质间透水或不透水边界条件以及不同饱和多孔介质交界面边界条件,形成基岩上覆成层饱和介质系统的整体有限元方程,最后采用中心差分法与Newmark平均加速度近似格式相结合的方法对时间进行离散,得到节点的动力时程的显式表达. 典型场地的地震反应分析表明,本文方法的计算结果与传递矩阵法结合傅里叶变换的计算结果完全吻合,证明了其有效性.      

Acceleration wave speeds in a hypoplastic constitutive model

The paper presents a theoretical investigation of acceleration waves in a plastic material described by an incrementally non-linear hypoplastic constitutive equation. Speeds of plane acceleration waves and their dependence on the stress state are calculated. The spectrum of possible wave speeds is found to be continuous, which is in contrast to discrete wave speed spectra in incrementally linear models. Two types of ill-posedness are revealed, known as flutter ill-posedness and stationary discontinuity. The wave speed analysis is also performed by the characteristic method, leading to different equations compared to the acceleration wave approach. It is proved that for the considered hypoplastic constitutive equation both approaches give identical wave speed spectra.     

On the Propagation of Long‐Wave Perturbations in a Two‐Layer Free‐Boundary Rotational Fluid

A mathematical model for the propagation of longwave perturbations in a freeboundary shear flow of an ideal stratified twolayer fluid is considered. The characteristic equation defining the velocity of perturbation propagation in the fluid is obtained and studied. The necessary hyperbolicity conditions for the equations of motion are formulated for flows with a monotonic velocity profile over depth, and the characteristic form of the system is calculated. It is shown that the problem of deriving the sufficient hyperbolicity conditions is equivalent to solving a system of singular integral equations. The limiting cases of weak and strong stratification are studied. For these models, the necessary and sufficient hyperbolicity conditions are formulated, and the equations of motion are reduced to the Riemann integral invariants conserved along the characteristics.     

Material instabilities in inelastic saturated porous media under dynamic loadings

Instabilities in inelastic saturated porous media are investigated here for general three-dimensional states under dynamic loadings using a perturbation approach.Under quasi-static conditions, unbounded growth of perturbations is related to the emergence of stationary discontinuities under drained or undrained conditions, while under dynamic conditions, unbounded growth is related either to the emergence of stationary discontinuities (and these are set by drained conditions) or to the appearance of the flutter phenomenon (acceleration waves).For associative behaviour unbounded growth always corresponds to localization under drained conditions and the onset of growth of perturbations occurs here only through divergence growth. It is only for non-associative flow that unbounded growth may correspond to undrained localization in quasi-static conditions and to flutter under dynamic conditions.     

Shallow waves in a two-layer vortex fluid under a lid

A mathematical model of the vortex motion of an ideal two-layer fluid in a narrow straight channel is considered. The fluid motion in the Eulerian-Lagrangian coordinate system is described by quasilinear integrodifferential equations. Transformations of a set of the equations of motion which make it possible to apply the general method of studying integrodifferential equations of shallow-water theory, which is based on the generalization of the concepts of characteristics and the hyperbolicity for systems with operator functionals, are found. A characteristic equation is derived and analyzed. The necessary hyperbolicity conditions for a set of equations of motion of flows with a monotone-in-depth velocity profile are formulated. It is shown that the problem of sufficient hyperbolicity conditions is equivalent to the solution of a certain singular integral equation. In addition, the case of a strong jump in density (a heavy fluid in the lower layer and a quite lightweight fluid in the upper layer) is considered. A modeling that results in simplification of the system of equations of motion with its physical meaning preserved is carried out. For this system, the necessary and sufficient hyperbolicity conditions are given. Novosibirsk State University, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 3, pp. 68–80, May–June, 1999.     

Quasi-static versus dynamic failure instabilities in fluid-saturated porous mediaModes de rupture sous conditions quasi-statiques et dynamiques dans les milieux poreux saturés

Using a linear perturbation approach, we show that under quasi-static conditions, unbounded growth of perturbations coincides with localization under drained or undrained conditions. Under dynamic loadings, unbounded growth is related either to the emergence of stationary discontinuities (and these are set by drained conditions) or to the appearance of the flutter phenomenon (acceleration waves). For associative behaviour the inception of unbounded growth is always set (under both static and dynamic conditions) by the singularity of the drained acoustic tensor. It is only for non-associative flow that unbounded growth may correspond to undrained localization in quasi-static conditions and to flutter under dynamic conditions. To cite this article: A. Benallal, C. Comi, C. R. Mecanique 330 (2002) 339–345.     

Model of Nonlinear Evolution of Long-Wave Perturbations on an Ideally Conducting Jet with Current in a Longitudinal Magnetic Field. Collision of Magnetized Jets

Within the framework of the magnetohydrodynamic approach, a system of equations is derived for nonlinear evolution of long-wave axisymmetric perturbations on a conducting fluid jet with surface electric current, located along the axis of a conducting solid cylinder in a longitudinal magnetic field. The fluid is assumed to be inviscid, incompressible, and ideally conducting, like the cylinder walls. It is shown that, if the longitudinal field is uniform and the axial flow is shear-free, this system can be either hyperbolic or elliptic-hyperbolic, depending on problem parameters. The boundaries of hyperbolicity and ellipticity regions in the space of solutions are determined. In the hyperbolicity region, equations of characteristics and conditions on them are obtained. The problem of the decay of velocity discontinuity on the jet is considered. Conditions are found for the existence of a continuous self-similar solution in the hyperbolicity region, corresponding to collision of jets.     

Second sound and internal energy in solids

A non-linear thermodynamic model of heat-conducting anisotropic solid is elaborated which turns out to be in a conservative form. Then, through the associated main field variables, the symmetry and the hyperbolicity properties are investigated. As outstanding applications, the analysis of the grow of the wave discontinuities and the evaluation of the critical time are performed. Finally, the Rankine-Hugoniot conditions for the system of equations are given in detail.     

不可压饱和粘弹性多孔介质的变分原理及其无网格模拟

基于饱和多孔介质理论,在固相和液相微观不可压,固相骨架小变形且满足线性粘弹性积分型本构关系的假定下,建立了流体饱和粘弹性多孔介质动力响应的若干Gurtin型变分原理,包括Hu-Washizu变分原理.利用所建立的变分原理,导出了流体饱和粘弹性多孔介质动力响应无网格数值模拟的离散控制方程,此方程是一个关于时间的对称微分方程组,便于分析计算.作为数值例子,研究了流体饱和粘弹性多孔柱体的一维动力响应,数值结果揭示了流体饱和粘弹性多孔柱体中波的传播特性以及固相粘性的影响.     

Axisymmetric Thermoelastoplastic Stress State of Isotropic Solids of Revolution Under Impulsive Loading

Consideration is given to the solution of a dynamic problem for a solid of revolution with an arbitrary meridional section under impulsive thermomechanical loading inducing elastoplastic strains. The theory of small elastoplastic deformations is used. The constitutive equations are linearized by the variable-parameter method. The unloading process is described by a linear law. The solution technique involves the finite-element approximation in spatial coordinates and the finite-difference representation of time derivatives. Based on the principle of linear summation, recurrent relations are derived for successive evaluation of nodal displacements by an explicit scheme in time. Solution for cylinders and disks are presented to illustrate the influence of elastoplastic deformations on wave processes