非饱和多孔介质中混合元法的化学-热-渗流-力学耦合的本构模拟

在文献[1]中建立的多孔介质中化学-热-渗流-力学(CTHM)本构模型基础上,针对文献[2]建立的非饱和多孔介质中热-渗流-力学耦合分析的混合有限元方法,发展了非饱和多孔介质中混合元的化学-热-渗流-力学(CTHM)耦合本构模拟算法。采用非关联流动多重屈服准则模拟非饱和多孔介质的材料非线性行为。推导了u-pw-pa-T形式的包含了耦合率本构方程积分的向后欧拉映射算法和一致性弹塑性切线模量矩阵(单元刚度矩阵)的混合元一致性算法。本文给出了临界状态线(CSL)和状态边界面(SBS)两个屈服准则的一致性算法。数值结果显示了本文所发展的混合元耦合本构模拟算法在模拟由热、化学、力学荷载共同引起的多孔介质中化学-热-渗流-力学(CTHM)耦合行为的能力和有效性。     

脆性固体损伤和局部软化的有限元分析

本文研究混凝土、岩石一类工程中常用的应变软化材料的有限元分析方法。在作者以往有关粘塑性损伤本构模型的工作基础上，给出了一组便于有限元计算的本构方程表达式。包括损伤弹性矩阵和局部损伤软化矩阵，分别运用于计算硬化和软化阶段的有限元刚度矩阵；对所提出的本构方程的实验验证计算和一些算例的有限元数值分析，表明文中给出的本构方程是可行的，相应的有限元算法能较好地对损伤固体的局部软化效应进行数值分析，并可成功地追踪应力应变响应的软化曲线     

黏弹-Perzyna黏塑性有限元法应力更新隐式算法

黏弹-黏塑性耦合模型的黏弹性部分由弹簧、黏壶和Kelvin链串联而成,黏塑性部分为双曲线型DruckerPrager屈服函数、各向同性硬化和Perzyna黏塑性流动模型。基于黏弹性蠕变柔度,通过定义与弹性问题相对应的与时间增量相关的黏弹性剪切模量和体积模量,导出增量递推形式的本构方程。为保证算法的收敛和稳定性,把Perzyna黏塑性流动方程转化为与弹塑性相似的一致性条件,建立黏塑性增量因子单侧逼近其收敛值的N-R迭代算法。最后,给出应力更新完全隐式算法和最终计算公式。分别采用黏弹性、黏弹-塑性和黏弹-黏塑性本构关系对一地基蠕变模型进行三维有限元分析和比较,结果表明,本文算法具有较高的计算效率和稳定性。     

一种黏弹塑性统一本构模型

经过对大量有关统一本构模型的文献资料分析，指出了现有统一本构模型存在的问题，并通过对材料实验数据的分析，指出了黏弹性现象在实验中的表现，并据此将黏弹性引入到弹性黏塑性统一本构模型之中，建立了黏弹塑性统一本构模型，通过模型的数值模拟证明：模型计算结果无论在变形趋势上，还是在数值精度上都与实验数据符合得很好，克服了此前统一本构模型存在的问题。黏弹塑性统一本构模型的产生将统一本构模型的产生将统一本构理论的内涵扩大到黏弹性范围，进而构造了一个黏弹塑性理论的新框架。     

复杂载荷作用下的结构弹塑性有限元分析

本文首先根据描写复杂载荷作用下材料特征行为的需要, 提出结合Dafalias双面模型和广义Hodge混合硬化模型两者特点的新的本构模型,并给出此模型本构方程的似近解析积分表达式。针对复杂载荷作用下的结构弹塑性分析的特点,本文还提出若干算法措施,以保证计算机执行过程的可靠性和有效性。实际算例的结果表明本文所提出的本构模型及算法措施是可行的。     

高温条件下混凝土墙体的非线性分析研究

本文建立了一个弹塑性-损伤耦合本构模型用于数值模拟高温下混凝土的真实破坏过程。导出了一个利用Newton-Raphson迭代的一般的直接应力返回映射算法。同时求解应力向量和塑性、损伤的内状态变量。并推导了用于化学-热-湿-力学耦合分析的全局守恒方程Newton-Raphson迭代过程的一致性切线模量矩阵。建议了一个用于弹塑性-损伤耦合分析的两级求解过程。给出的数值例题结果显示了所提出的数法和公式的正确性,表明了所发展的弹塑性-损伤耦合本构模型在模拟高温下混凝土墙体中复杂破坏过程的能力。     

结构性软土弹塑性模型的隐式算法实现

对于考虑软土结构性的高度非线性弹塑性本构模型,在采用Newton-CPPM隐式算法对模型进行数值实现的过程中容易出现Jacobian矩阵奇异和不收敛问题。为此,本文提出了两种改进隐式算法。考虑到Newton-CPPM隐式算法是局部收敛性算法,因此引入大范围收敛的同伦延拓算法对Newton-CPPM算法的迭代初值进行改进,形成了同伦-Newton-CPPM算法。考虑到Newton-CPPM隐式算法单个迭代步的计算量过大,因此借鉴显式算法的思想提出一种两阶段迭代算法,第一阶段先求出一致性参数,第二阶段采用类似于显示算法的方法进行回代得出状态变量的值。然后,以考虑软土结构性的SANICLAY模型为例,从弹塑性本构模型的组成和算法的特点两个角度分析了引起Jacobian矩阵奇异和不收敛问题的原因,并且在单单元计算的基础上,对全显式算法、传统隐式算法和两种改进隐式算法在计算收敛性、计算精度和计算效率方面进行了对比。最后,将同伦-Newton-CPPM算法和传统隐式算法用于地基承载力多单元计算中,结果表明该算法能够有效地解决Jacobian矩阵奇异和不收敛问题。      

一种混凝土损伤模型和数值方法

Ｏｔｏｓｅｎ准则在准静态时和实验结果有很好的一致性，但是混凝土的动态力学性质和准静态相比有明显的变化。此时，需要用应变率相关的本构模型来描述混凝土的力学行为。从Ｏｔ－ｔｏｓｅｎ的四参数混凝土破坏准则出发，考虑损伤、静水压和应变率对本构关系的影响，建立了混凝土的粘塑性本构模型。同时给出了基于该本构模型的混凝土的有限元计算方法:在积分内变量时采用改进的龙格－库塔格式，在时间方向上的积分使用带有步长控制的ｇｅｎｅｒａｌｉｚｅｄ－方法，有效地保证了积分的稳定和精度，给混凝土的进一步研究提供了方便。     

非比例循环载荷下塑性模量的讨论

对1Cr18Ni9Ti不锈钢进行菱形应变路径的非比例循环本构实验,根据实验结果及分析,对双面本构模型给出了一个极限面定义,以保证屈服面不与极限面相交,在这个极限面定义及实验结果的基础上,对不同本构模型所定义的距离与塑性模量的相关性进行了讨论,表明Mroz和Chen的距离与塑性模量具有较好的相关性。     

????????????????????????????

基于群体与适应度的概念，应用改进的PSO算法，从随机解出发，提出了基于PSO 算法的本构模型参数识别方法. 该方法解决了橡胶类材料大应变时硬化现象的本构模型参数 的确定这一难题. 首先通过单轴拉伸本构模型实验，在针对硬化实验曲线存在拐点的情 况下，应用PSO算法进行拟合，最后利用简单剪切实验进行验证. 结果表明该方法科学可行， 且具有速度快、精度高、易于收敛等优点. 并且有效地解决了本构模型参数识别的困难，可 广泛应用于各种复杂材料.     

饱和多孔介质中的混合有限元法和有限应变下应变局部化分析

对基于Biot理论的饱和多孔介质中动力-渗流耦合分析提出了一个耦合场混合元．固相位移、应变和有效应力以及流相压力、压力梯度和Darcy速度在单元内均处理为独立变量分别插值．基于胡海昌-Washizu三变量广义变分原理给出的饱和多孔介质动力-渗流耦合问题控制方程的单元弱形式，导出了单元公式．进一步导出了考虑压力相关非关联塑性的非线性单元公式和发展了相应的一致性算法．对几何非线性分析，采用了共旋公式途径．数值结果例题显示所发展耦合场混合元模拟大应变下由应变软化引起以应变局部化为特征的渐进破坏现象的性能．     

基于显式有限元方法的两相介质弹塑性动力反应计算分析

针对增量形式的流体饱和两相多孔介质弹塑性波动方程组,运用基于显式逐步积分格式的时域显式有限元方法对该波动方程组进行求解,并应用基于SMP破坏准则的弹塑性动力本构模型描述两相介质的动力反应性质,对两相介质在输入地震波作用下的弹塑性动力反应进行计算和分析,将计算结果与相应的弹性动力反应的计算结果进行对比;对本文应用的弹塑性...     

CONSTITUTIVE RELATION OF UNSATURATED SOIL BY USE OF THE MIXTURE THEORY (Ⅱ)—LINEAR CONSTITUTIVE EQUATIONS AND FIELD EQUATIONS

IntroductionInpart(Ⅰ )ofthework[1],byuseofmixturetheory ,thenonlinearconstitutiveequationsandthefieldequationsofunsaturatedsoilwereconstructed ,andthecompleteequationsforthethermodynamicsystemofunsaturatedsoilwasformed .Inthispart,thelinearconstitutiveequationsandfieldequationsofunsaturatedsoilareobtainedthroughlinearizingnonlinearequations,andthelinearequationsarewrittenintheformssimilartoBiot’sequationsforsaturatedporousmedia .ItisprovedthatDarcy’slawissuitabletodescribethemotionofliquid…     

A thermodynamical gradient theory for deformation and strain localization of porous media

In this work, a thermodynamically consistent gradient formulation for partially saturated cohesive-frictional porous media is proposed. The constitutive model includes a classical or local hardening law and a softening formulation with state parameters of non-local character based on gradient theory. Internal characteristic length in softening regime accounts for the strong shear band width sensitivity of partially saturated porous media regarding both governing stress state and hydraulic conditions. In this way the variation of the transition point (TP) of brittle-ductile failure mode can be realistically described depending on current confinement condition and saturation level. After describing the thermodynamically consistent gradient theory the paper focuses on its extension to the case of partially saturated porous media and, moreover, on the formulation of the gradient-based characteristic length in terms of stress and hydraulic conditions. Then the localization indicator for discontinuous bifurcation is formulated for both drained and undrained conditions.     

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.     

A micropolar formulation of the Desai hierarchical model for elastoplastic porous media

In this contribution, the Desai hierarchical model is extended to the case of 3D elasto-plastic two-phase micropolar continuum. An unconditionally stable implicit Euler backward algorithm for integration of the constitutive relations is developed and presented in detail. The regularizing effect of the introduction of a length parameter to the classical Desai model is demonstrated by means of several case studies of single phase and two-phase porous media.     

STUDY ON ANALYSIS METHOD OF “THREE BOX” IN HEAT AND MASS TRANSFER PROCESS OF RESERVOIR MEDIA

油藏多孔介质孔隙组成及结构变化多样,一些特性参数很难全部获得,精确描述和分析困难;另外,多孔介质内渗流过程水力条件和作用机理复杂,存在热流固耦合作用,目前的一些分析方法和研究模型具有一定的局限性.提出了油藏多孔介质的表征单元体（representative elementary volume,REV）描述表征方法;基于表征单元体建立了多孔介质的黑箱模型、灰箱模型和白箱模型,据此提出了多孔介质的“黑箱→灰箱→白箱”分析过程.基于黑箱模型和灰箱模型推导了REV导热系数计算公式、给出了REV热质传递过程的热平衡方程.结合中国油藏热采情况,对多孔介质导热系数变化规律和蒸汽驱热质传递特性进行了分析,得到了一些有意义的结果.该工作为多孔介质热质传递过程分析提供了新思路和新方法.     

ONSET CONDITION OF STRAIN LOCALIZATION IN MATRIX OF SATURATED POROUS MEDIA

Introduction Strainlocalizationofgeomaterialsisoneofmostpopularfailuretypesinnature,which canbeshowedaslandslidesandmudflowsinmountainousareasunderincessantorheavy raining,especiallythevegetationisseverelydamagedbywoodsharvest;pipingeffect,a typeoflocalfa…     

Stationary discontinuity and flutter instability of wave propagation in elasto-plastic saturated porous media

In the present work a model based on the Biot theory for simulating coupled hydrodynamic behavior in saturated porous media is utilized with integration of the inertial coupling effect between the solid-fluid phases of the media into the model. The non-associated Drucker-Prager criterion to describe nonlinear constitutive behavior of pressure dependent elasto-plasticity for the media is particularly considered. With no consideration of compressibility of solid grains and the pore fluid, the discontinuity and instability of the wave propagation in saturated porous media are analyzed for the plane strain problems in detail. The critical conditions of stationary discontinuity and flutter instability in the wave propagation are given. The results and conclusions obtained by the present work will provide some bases or clues for overcoming the difficulties in numerical modeling of wave propagation in the media subjected to dynamic loading. The project supported by the National Natural Science Foundation of China (19832010)