岩体中饱和渗流应力耦合模型研究进展

在大坝渗流与控制、各种水工隧洞和交通隧道以及地下各种洞室开挖、采矿工程与油藏开采中的水力劈裂、 岩坡和坝基的稳定性研究等领域中经常会遇到岩体的饱和渗流应力耦合问题.本文对这些领域中的渗流应力耦合问题分为6种力学模型: 等效连续介质模型、 裂隙网络模型、 双重介质模型、 断裂力学模型、 连续损伤力学模型和统计模型.前3种主要从经典的黏弹塑性本构着手考虑, 着重对渗流场的处理, 后3种则是从考虑岩体在耦合作用下发生损伤破裂行为方面切入, 着重于岩体结构内部发生质的改变带来的更加复杂的耦合效应.对这些模型进行了详细的介绍, 指出了这些模型的优缺点与适用范围, 给出了一些有代表性的研究成果, 并对该研究领域中未来的研究方向进行了展望.     

The Three-Stage Model Based on Strain Strength Distribution for the Tensile Failure Process of Rock and Concrete Materials

A three-stage model is introduced to describe the tensile failure process of rock and concrete materials.Failure of the material is defined to contain three stages in the model,which include elastic deformation stage,body damage stage and localization damage stage.The failure mode change from uniform body damage to localization damage is expressed.The heterogeneity of material is described with strain strength distribution.The fracture factor and intact factor,defined as the distribution function of strain strength,are used to express the fracture state in the failure process.And the distributive parameters can be determined through the experimental stress-strain curve.     

A stochastic model for the permeability characteristics of saturated cemented porous media undergoing freezing

As the temperature of a saturated porous medium drops, the water in the pores starts to freeze. Since the temperature at which the phase change takes place is dependent on the pore size, the permeability of the medium changes continuously. Simultaneously, due to the expansion of water on freezing, it is forced to migrate through the pore body thus inducing stresses in material matrix. The stresses developed and the consequent frost damage are therefore dependent on the change in the permeability characteristics of the medium on freezing. This paper deals with the numerical prediction of permeability characteristics of porous cemented media saturated with water undergoing progressive freezing.A bond percolation model is used to generate the pore structure according to an assumed poresize distribution. Permeability of the medium at various temperatures is computed by solving the network problem. The computed results are compared with other analytical and experimental results. The proposed model predicts a threshold temperature below which permeability drops to zero. This phenomenon is crucial in developing a deeper understanding of the mechanism of frost damage to cemented porous materials such as bricks, stone, concrete, etc.     

Scattered fracture of porous materials with brittle skeleton

A model of damage accumulation in a porous medium with a brittle skeleton saturated with a compressible fluid is formulated in the isothermal approximation. The model takes account of the skeleton elastic energy transformation into the surface energy of microcracks. In the case of arbitrary deformations of an anisotropic material, constitutive equations are obtained in a general form that is necessary and sufficient for the objectivity and thermodynamic consistency principles to be satisfied. We also formulate the kinetics equation ensuring that the scattered fracture dissipation is nonnegative for any loading history. For small deviations from the initial state, we propose an elastic potential which permits describing the principal characteristics of the behavior of a saturated porous medium with a brittle skeleton. We study the acoustic properties of the material under study and find their relationship with the strength criterion depending on the accumulated damage and the material current deformation. We consider the problem of scattered fracture of a saturated porous material in a neighborhood of a spherical cavity. We show that the cavity failure occurs if the Hadamard condition is violated.     

单轴拉伸条件下脆性岩石微裂纹损伤模型研究

利用断裂力学、损伤力学和均匀化原理,对脆性岩石单轴拉伸条件下的力学特性进行分析,建立了脆性岩石的微裂纹损伤本构模型.首先对岩石内部微裂纹的统计分布规律进行分析,给出了理论分析过程中微裂纹分布的假设条件,在此基础上,参考已有研究成果,得到含细长微裂纹脆性岩石有效弹性参数的计算公式.然后,对岩石内部单一微裂纹进行断裂力学和损伤力学分析,得到了扩展裂纹尖端的应力强度因子计算公式,在一定微裂纹断裂扩展准则和断裂扩展速率的假设基础上,利用积分原理,得到了岩石整体的损伤变量和损伤演化方程,由此建立单轴拉伸条件下脆性岩石的微裂纹损伤本构模型.最后,通过一花岗岩的单轴拉伸试验结果对微裂纹损伤本构模型进行了验证.     

强地震荷载作用下含断层岩质边坡稳定性分析

基于岩石材料的弹塑性损伤理论,引入两个独立的损伤变量分别控制岩石的拉压行为,通过对地震荷载作用下边坡场地响应的有限无数值模拟,总结含裂隙的岩质边坡发生局部损伤软化、破坏的区域的扩展、演化过程和规律.结合对某大型水电坝坝址左岸的岩质边坡在强地震动载荷作用下(百年一遇、超越概率2%)发生破坏的数值模拟,揭示这类边坡滑坡的机理,并根据模拟的结果,对整治边坡破坏提出参考依据.     

STUDY ON NONLINEAR CONSTITUTIVE RELATIONSHIP AND FRACTURE BEHAVIOR OF STITCHED C/SiC COMPOSITES$^{\bf 1)}$

C/SiC复合材料具有高比强度、高比模量和优良的热稳定性能等一系列优点, 广泛应用于航空航天领域中. 裂纹扩展进而引起的脆性断裂是其主要失效形式之一, 因而材料的断裂性能分析对材料的结构设计和应用有重要的指导意义. 本文开展了缝合式C/SiC复合材料简单力学试验和断裂试验, 研究了材料在不同载荷下的力学响应及断裂特征. 基于缝合式C/SiC复合材料简单力学试验, 建立了材料宏观非线性损伤本构方程, 并模拟了缝合式C/SiC复合材料单边切口梁和双悬臂梁的断裂行为. 本构方程采用简单函数描述了材料在复杂应力状态下的非线性应力-应变曲线, 并考虑了反向加载过程中造成的裂纹闭合. 基于商业有限元软件ABAQUS, 通过编写UMAT子程序实现非线性损伤本构方程, 采用单个单元验证了建立的本构方程的有效性. 在此基础上, 采用线弹性损伤本构和非线性损伤本构分别模拟了缝合式C/SiC复合材料单边切口梁和双悬臂梁的断裂行为. 采用非线性损伤本构方程模拟的力-位移曲线结果与试验结果更为吻合, 非线性损伤本构预测的失效载荷与试验失效载荷更为接近, 验证了所建立的非线性损伤本构方程的准确性, 为C/SiC复合材料断裂行为的研究提供了借鉴, 为缝合式C/SiC复合材料结构的设计和应用提供了理论基础.     

缝合式C/SiC复合材料非线性本构关系及断裂行为研究

C/SiC复合材料具有高比强度、高比模量和优良的热稳定性能等一系列优点, 广泛应用于航空航天领域中. 裂纹扩展进而引起的脆性断裂是其主要失效形式之一, 因而材料的断裂性能分析对材料的结构设计和应用有重要的指导意义. 本文开展了缝合式C/SiC复合材料简单力学试验和断裂试验, 研究了材料在不同载荷下的力学响应及断裂特征. 基于缝合式C/SiC复合材料简单力学试验, 建立了材料宏观非线性损伤本构方程, 并模拟了缝合式C/SiC复合材料单边切口梁和双悬臂梁的断裂行为. 本构方程采用简单函数描述了材料在复杂应力状态下的非线性应力-应变曲线, 并考虑了反向加载过程中造成的裂纹闭合. 基于商业有限元软件ABAQUS, 通过编写UMAT子程序实现非线性损伤本构方程, 采用单个单元验证了建立的本构方程的有效性. 在此基础上, 采用线弹性损伤本构和非线性损伤本构分别模拟了缝合式C/SiC复合材料单边切口梁和双悬臂梁的断裂行为. 采用非线性损伤本构方程模拟的力-位移曲线结果与试验结果更为吻合, 非线性损伤本构预测的失效载荷与试验失效载荷更为接近, 验证了所建立的非线性损伤本构方程的准确性, 为C/SiC复合材料断裂行为的研究提供了借鉴, 为缝合式C/SiC复合材料结构的设计和应用提供了理论基础.      

Constitutive relationship of brittle rock subjected to dynamic uniaxial tensile loads with microcrack interaction effects

A micromechanical model is proposed to describe both stable and unstable damage evolution in microcrack-weakened brittle rock material subjected to dynamic uniaxial tensile loads. The basic idea of the present model is to classify the constitution relationship of rock material subjected to dynamic uniaxial tensile loads into four stages including some of the stages of linear elasticity, pre-peak nonlinear hardening, rapid stress drop, and strain softening, and to investigate their corresponding micromechanical damage mechanisms individually. Special attention is paid to the transition from structure rearrangements on microscale to the macroscopic inelastic strain, to the transition from distribution damage to localization of damage and the transition from homogeneous deformation to localization of deformation. The influence of all microcracks with different sizes and orientations are introduced into the constitutive relation by using the statistical average method. Effects of microcrack interaction on the complete stress-strain relation as well as the localization of damage for microcrack-weakened brittle rock material are analyzed by using effective medium method. Each microcrack is assumed to be embedded in an approximate effective medium that is weakened by uniformly distributed microcracks of the statistically-averaged length depending on the actual damage state. The elastic moduli of the approximate effective medium can be determined by using the dilute distribution method. Micromechanical kinetic equations for stable and unstable growth characterizing the ‘process domains’ of active microcracks are taken into account. These ‘process domains’ together with ‘open microcrack domains’ completely determine the integration domains of ensemble averaged constitutive equations relating macro-strain and macro-stress. Theoretical predictions have shown to consistent with the experimental results.     

非贯通裂隙岩体损伤演化率相关性及变形特征

含非贯通裂隙岩体是自然界中岩体的主要赋存形式，其裂隙几何特征对岩体的强度及变形均产生显著影响。应变率对岩体的损伤演化及黏滞效应也具有显著的率相关性。首先，运用模型元件的方法，将非贯通裂隙岩体动态破坏过程视为具复合损伤、静态弹性特性、动态黏滞特性的非均质点组成，对黏弹性响应的Maxwell体进行改进，将细观损伤体与裂隙损伤演化的宏观损伤体根据等效应变假设并联组成宏细观复合损伤体，构建综合考虑岩体宏细观缺陷的动态损伤模型；其次，基于断裂力学及应变能理论，对岩体宏观裂隙动态扩展的能量机制进行分析，综合考虑初始裂隙应变能、裂隙动态损伤演化过程应变能、裂隙闭合应变能，得到裂隙岩体宏观动态损伤变量计算公式；最后，将模型计算结果与实验结果进行比较，模型计算结果与实验结果吻合较好，证明了模型的合理性，同时利用模型讨论了裂隙倾角、应变率、岩石性质对岩体变形特征的影响规律。     

Localization of deformation and stress–strain relation for mesoscopic heterogeneous brittle rock materials under unloading

Stress redistribution induced by excavation of underground engineering and slope engineering results in the unloading zone in parts of surrounding rock masses. The mechanical behaviors of crack-weakened rock masses under unloading are different from those of crack-weakened rock masses under loading. A micromechanics-based model has been proposed for brittle rock material undergoing irreversible changes of their microscopic structures due to microcrack growth when axial stress is held constant while lateral confinement is reduced. The basic idea of the present model is to classify the constitution relation of rock material into four stages including some of the stages of linear elasticity, pre-peak nonlinear hardening, rapid stress drop, and strain softening, and to investigate their corresponding micromechanical damage mechanisms individually. Special attention is paid to the transition from structure rearrangements on microscale to the macroscopic inelastic strain, to the transition from distribution damage to localization of damage and the transition from homogeneous deformation to localization of deformation. The closed-form explicit expression for the complete stress–strain relation of rock materials containing cracks under unloading is obtained. The results show that the complete stress–strain relation and the strength of rock materials under unloading depend on the crack spacing, the fracture toughness of rock materials, orientation of the cracks, the crack half-length and the crack density parameter.     

修正金属本构模型在超高速撞击模拟中的应用

为更加准确地计算93钨合金弹超高速撞击Q345钢板问题,构建了修正的金属本构模型。引入GRAY三相物态方程描述材料相态变化,采用Johnson-Cook强度模型描述撞击后期材料的力学行为。结合封加波损伤演化模型以及Johnson-Cook失效模型描述不同应力三轴度下材料的拉伸、剪切失效行为;引入曹祥提出的断裂演化模型,描述材料失效后应力归零的过程。通过对比超高速撞击数值模拟结果与实验结果,验证了本构模型的适用性,并进一步分析了典型弹靶撞击条件下破片群的空间分布特征。研究结果表明:基于修正金属本构模型获得的超高速撞击靶板穿孔直径、弹体侵蚀长度、破片群扩展速度结果与实验结果一致;GRAY三相物态方程能够相对准确地给出弹体撞击首层靶板以及剩余弹体、破片群撞击第2层靶板时弹靶材料的熔化情况;封加波损伤演化模型能够准确判断超高速撞击过程中靶板是否产生层裂破坏;综合封加波损伤演化模型、Johnson-Cook失效模型以及曹祥提出的断裂演化模型后,数值模拟获得的破片群撞击后效靶板的穿孔面积与累积数量的统计曲线结果与实验结果一致;获得了典型条件下的柱形93钨弹体超高速撞击Q345靶板破片群空间分布结果,破片群的前端具有较高的质量、轴向动量以及横向动量（绝对值）。     

Dynamic void nucleation and growth in solids: A self-consistent statistical theory

We present a framework for a self-consistent theory of spall fracture in ductile materials, based on the dynamics of void nucleation and growth. The constitutive model for the material is divided into elastic and “plastic” parts, where the elastic part represents the volumetric response of a porous elastic material, and the “plastic” part is generated by a collection of representative volume elements (RVEs) of incompressible material. Each RVE is a thick-walled spherical shell, whose average porosity is the same as that of the surrounding porous continuum, thus simulating void interaction through the resulting lowered resistance to further void growth. All voids nucleate and grow according to the appropriate dynamics for a thick-walled sphere made of incompressible material. The macroscopic spherical stress in the material drives the response in all volume elements, which have a distribution of critical stresses for void nucleation, and the statistically weighted sum of the void volumes of all RVEs generates the global porosity. Thus, macroscopic pressure, porosity, and a distribution of growing microscopic voids are fully coupled dynamically. An example is given for a rate-independent, perfectly plastic material. The dynamics of void growth gives rise to a rate effect in the macroscopic material even though the parent material is rate independent.     

Fracture behaviour in tension of viscoplastic porous metallic materials saturated with liquid

The present work extends the investigation which has been initiated in Parts I and II of this study (Martin, C.L., Favier, D., Suéry, M., 1997a. Viscoplastic behaviour of porous metallic materials saturated with liquid, part I: constitutive equations. Int. J. Plasticity 13, 215–235; Martin, C.L., Favier, D. Suéry, M., 1997b. Viscoplastic behaviour of porous metallic materials saturated with liquid, Part II: experimental identification on a Sn–Pb model alloy. Int. J. Plasticity 13, 237–259) to the tensile behaviour of viscoplastic porous metallic materials saturated with liquid. Simple tensile experiments together with ring extension tests are carried out to study the fracture behaviour of this class of material. Ring tests consist in applying an internal pressure on a specimen with a ring shape. A Sn–Pb model alloy with a dendritic microstructure is used to characterise the behaviour of the material up to fracture. The liquid presence is accounted for to derive the intrinsic behaviour of the solid skeleton. The collected data are then incorporated in the model framework presented in Part I. A simple modification of the model allows the treatment of the strong asymmetry between tension and compression which is exhibited by these materials.     

Microcrack interaction brittle rock subjected to uniaxial tensile loads

A micromechanics-based model is proposed to describe unstable damage evolution in microcrack-weakened brittle rock material. The influence of all microcracks with different sizes and orientations are introduced into the constitutive relation by using the statistical average method. Effects of microcrack interaction on the complete stress–strain relation as well as the localization of damage for microcrack-weakened brittle rock material are analyzed by using effective medium method. Each microcrack is assumed to be embedded in an approximate effective medium that is weakened by uniformly distributed microcracks of the statistically-averaged length depending on the actual damage state. The elastic moduli of the approximate effective medium can be determined by using the dilute distribution method. Micromechanical kinetic equations for stable and unstable growth characterizing the ‘process domains’ of active microcracks are taken into account. These ‘process domains’ together with ‘open microcrack domains’ completely determine the integration domains of ensemble averaged constitutive equations relating macro-strain and macro-stress. Theoretical predictions have shown to be consistent with the experimental results.     

准脆性材料的弹塑性各向异性损伤耦合本构模型研究

针对准脆性材料的非线性特征：强度软化和刚度退化、单边效应、侧限强化和拉压软化、不可恢复变形、剪胀及非弹性体胀,在热动力学框架内,建立了准脆性材料的弹塑性与各向异性损伤耦合的本构关系。对准脆性材料的变形机理和损伤诱发的各向异性进行了诠释,并给出了损伤构形和有效构形中各物理量之间的关系。在有效应力空间内,建立了塑性屈服准则、拉压不同的塑性随动强化法则和各向同性强化法则。在损伤构形中,采用应变能释放率,建立了拉压损伤准则、拉压不同的损伤随动强化法则和各向同性强化法则。基于塑性屈服准则和损伤准则,构建了塑性势泛函和损伤势泛函,并由正交性法则,给出了塑性和损伤强化效应内变量的演化规律,同时,联立塑性屈服面和损伤加载面,给出了塑性流动和损伤演化内变量的演化法则。将损伤力学和塑性力学结合起来,建立了应变驱动的应力-应变增量本构关系,给出了本构数值积分的要点。以单轴加载-卸载往复试验识别和校准了本构材料常数,并对单轴单调试验、单轴加载-卸载往复试验、二轴受压、二轴拉压试验和三轴受压试验进行了预测,并与试验结果作了比较,结果表明,所建本构模型对准脆性材料的非线性材料性能有良好的预测能力。     

基于微结构损伤的珠光体材料本构行为研究

珠光体材料由许多随机取向的珠光体团组成，每一珠光体团又由铁素体片和渗碳体片交替 叠合而成. 实验表明具有较小片层间距的珠光体材料具有较高的强度和寿命. 通过分析构成 珠光体团的各相及界面的损伤特征，利用损伤耗散功的概念，得到了铁素体、渗碳体与界面 损伤的具有统一形式的演化规律，将其嵌入基于细片层微结构和非经典塑性理论的珠光体团 的弹塑性本构方程，利用Hill自洽方法，得到了珠光体材料的损伤本构描述. 所得到的损 伤本构描述显含表征微结构特征的片层间距，容易证明具有较小片层间距的珠光体材料具有 较好的综合力学性能. 对珠光体材料的拉伸和循环拉压特性进行了分析，得到了与实验一 致的结果.     

岩土介质的分形孔隙和分形粒子

岩土介质是晶粒状材料,存在大量的孔隙。这些孔隙的存在严重地影响岩土介质的力学、物理和化学性能。大量的研究表明,岩土介质的孔隙几何从原子尺度到晶粒尺寸范围内均表现出分形特征。目前分形几何已被广泛应用来研究岩土的孔隙率,输运特性和渗透性等等。本文从5个方面来讨论岩土介质的分形孔隙和分形粒子:①分形几何简介;②孔隙介质的分形模型;③岩土介质的分形孔隙特性和它们的分形量测方法;④岩土的分形粒子;⑤水土保持估计中的分形毛细管模型。      

多相孔隙介质的本构描述

在混合物理论的框架内,提出了一个处理和和非饱和的多相孔隙介质（工程材料）变形和强度本构建模的新方法,并给出了相应的具体描述,该方法将多相工程材料看成由两个层次的３种混合物模型组成,即第一层次的“固相”和“液相”混合物模型和第二层次的由“固相”和“液相”组成的混合物模型,用该模型来模拟要研究的多相工程材料,这里定义的“固相”和“液相”实际上也是由各自组分组成的混合物,用以反映“固相”介质和“液相”介     

Phase-field modeling of hydraulic fracture

In this work a theoretical framework implementing the phase-field approach to fracture is used to couple the physics of flow through porous media and cracks with the mechanics of fracture. The main modeling challenge addressed in this work, which is a challenge for all diffuse crack representations, is on how to allow for the flow of fluid and the action of fluid pressure on the aggregate within the diffuse damage zone of the cracks. The theory is constructed by presenting the general physical balance laws and conducting a consistent thermodynamic analysis to constrain the constitutive relationships. Constitutive equations that reproduce the desired responses at the various limits of the phase-field parameter are proposed in order to capture Darcy-type flow in the intact porous medium and Stokes-type flow within open cracks. A finite element formulation for the solution of the governing model equations is presented and discussed. Finally, the theoretical and numerical model is shown to compare favorably to several important analytical solutions. More complex and interesting calculations are also presented to illustrate some of the advantageous features of the approach.