一种新的岩石非线性黏弹塑性蠕变模型研究

为了克服传统元件组合模型不能描述岩石蠕变过程中非线性特征的缺陷，首先根据加速蠕变阶段的应变和应变率随蠕变时间急剧增大的特点，建立黏塑性应变与蠕变时间的指数函数关系并提出非线性黏塑性体．将该非线性黏塑性体与广义Burgers蠕变模型串联，建立可以描述岩石全蠕变过程的非线性黏弹塑性蠕变模型，根据叠加原理得到一维应力状态下的轴向蠕变方程．然后基于塑性力学理论指出岩石三维蠕变本构方程建立过程中的不足之处，并给出非线性黏弹塑性蠕变模型合理的三维蠕变方程．最后采用不同应力水平下砂岩轴向蠕变试验对模型合理性进行验证，结果表明：拟合曲线与试验曲线吻合度较高，所建蠕变模型能够很好地描述砂岩在不同应力水平下的蠕变变形规律，尤其对加速蠕变阶段的非线性特征描述效果很好，验证了模型的合理性．     

I型定常扩展裂纹尖端的弹黏塑性场

考虑材料在扩展裂纹尖端的黏性效应，假设黏性系数与塑性应变率的幂次成反比，对幂硬化材料中平面应变扩展裂纹尖端场进行了弹黏塑性渐近分析，得到了不含间断的连续解，并讨论了I型裂纹数值解的性质随各参数的变化规律. 分析表明应力和应变均具有幂奇异性，并且只有在线性硬化时，尖端场的弹、黏、塑性才可以合理匹配. 对于I型裂纹，裂尖场不含弹性卸载区. 当裂纹扩展速度趋于零时，动态解趋于准静态解，表明准静态解是动态解的特殊形式；如果进一步考虑硬化系数为零的极限情况，便可退化为Hui和Riedel的非线性黏弹性解.     

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

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

Ⅱ型动态扩展裂纹尖端的弹黏塑性渐近场

考虑材料的黏性效应建立了Ⅱ型动态扩展裂纹尖端的力学模型,假设黏性系数与塑性等效应变率的幂次成反比,通过分析使尖端场的弹、黏、塑性得到合理匹配,并给出边界条件作为扩展裂纹定解的补充条件,对理想塑性材料中平面应变扩展裂纹尖端场进行了弹黏塑性渐近分析,得到了不含间断的连续解,并讨论了Ⅱ型裂纹数值解的性质随各参数的变化规律.分析表明应力和应变均具有幂奇异性,对于Ⅱ型裂纹,裂尖场不含弹性卸载区.引入Airy应力函数,求得了Ⅱ型准静态裂纹尖端场的控制方程,并进行了数值分析,给出了裂纹尖端的应力应变场.当裂纹扩展速度(M→0)趋于零时,动态解趋于准静态解,表明准静态解是动态解的特殊形式.     

混凝土黏塑性动力损伤本构关系

从静力弹塑性损伤本构关系的基本框架出发, 综合考虑塑性应变与损伤 演化的率敏感性, 建立了能够较为全面地描述混凝土在动力加载条件下非线性性能的混凝土 黏塑性动力损伤本构模型. 为了考虑塑性应变的率敏感性, 基于Perzyna理论推导了有效应 力空间黏塑性力学基本公式, 采用改进的Perzyna型动力演化方程, 将损伤静力演化方程推 广到动力加载情形. 基于并联弹簧模型, 从概率论的角度推导给出了一维损伤静力演化方程, 并基于能量等效应变的基本概念将其推广到多维损伤演化. 利用数值模拟, 计算得到了 混凝土在不同应变率下的应力应变全曲线, 同时得到了一维动力提高因子和二维动力强度包 络图, 数值结果与试验结果的对比表明了该模型的有效性.     

聊聊动态塑性和黏塑性

固体力学研究者致力于具有应力-应变本构关系（以下简称为形变型本构关系）的变形体的力学响应研究，而流体力学研究者致力于具有应力-应变率本构关系（以下简称为流动型本构关系）的流动体的力学响应研究。当涉及结构和材料的动态塑性时，到底应该用“塑性变形”还是“塑性流动”来表示？本文从宏观塑性本构理论和微观位错动力学机理两个角度，分别讨论并指出塑性本构关系属于流动型黏塑性率相关本构关系，且同时适用于加载和卸载；因而不应该用应力-应变图来描述塑性加-卸载过程。弹塑性本构关系则是一种形变型和流动型本构关系的耦合。     

NEPE推进剂低高应变率下改进的黏-超弹本构模型

为研究低高应变率条件下NEPE推进剂的力学特性,通过电子万能试验机和分离式霍普金森杆装置,对NEPE推进剂进行了准静态和冲击实验,得到了不同应变率下（1.667×10?4～4 500 s?1）的应力-应变曲线。实验结果表明NEPE推进剂具有明显的非线性弹性和应变率敏感性,随着应变率的增加,材料的强度、屈服应力和弹性模量显著增加,与低应变率相比,高应变率条件下材料的应变率敏感性更高。在高速冲击下材料内部瞬间产生大量热量无法及时散发出去,使得材料内部温度升高,导致材料出现软化效应,力学性能降低。本文建立了一个非线性黏超弹本构模型,其中采用Rivlin应变能函数来描述稳态超弹响应部分,采用积分型本构模型来描述材料的动态黏弹性响应部分,考虑到松弛时间具有应变率相关性,本文采用了一个率相关松弛函数来替代传统的Prony级数形式。使用极慢速压缩实验数据对本构模型中的超弹部分进行拟合获得超弹参数,然后用准静态和动态实验数据对本构模型进行拟合得出其他参数。不同应变率下的预测曲线与实验曲线具有较好的重合度,证明了该模型可以很好地描述低高应变率下NEPE推进剂的力学特性。     

双材料界面Ⅰ型扩展裂纹尖端的弹黏塑性场

双材料界面中存在材料黏性效应, 对界面裂纹尖端场的分布和界面本身性能 的变化起着重要的影响. 考虑裂纹尖端的奇异性, 建立了双材料界面扩展裂纹尖端的弹黏塑 性控制方程. 引入界面裂纹尖端的位移势函数和边界条件, 对刚性-弹黏塑性界面I型界面 裂纹进行了数值分析, 求得了界面裂纹尖端应力应变场, 并讨论了界面裂纹尖端场随各影响 参数的变化规律. 计算结果表明, 黏性效应是研究界面扩展裂纹尖端场时的一个主要因素, 界面裂纹尖端为弹黏塑性场, 其场受材料的黏性系数、马赫数和奇异性指数控制.     

黏弹性板的大挠度蠕变屈曲

研究了具有初始小挠度受轴向压载黏弹性板的蠕变屈曲问题，在建立控制方程时，利用了von Karman非线性应变-位移关系，并考虑了初始挠度，用标准线性固体模型描述材料的黏弹性特性，在求解非线性积分方程时，利用梯形公式计算记忆积分式，将非线性积分方程化为非线性代数方程进行数值求解，得到了结构的蠕变变形过程，又将问题退化到小挠度情况进行研究，得到了挠度随时间扩展的解析解，分析了瞬时失稳临界载荷、持久临界载荷的物理意义，讨论了考虑几何非线性对黏弹性板蠕变屈曲的影响。     

黏弹性模型微分型本构方程的矩阵形式

针对多个弹性元件和黏性元件以任意连接方式组成的线性黏弹性模型,本文探究了其本构方程的通用矩阵形式表述。首先将研究问题扩展为由Maxwell基本单元构成的标准模型,然后转化为有向图,根据独立路径和闭包围的形式表征出基本应力方程和应变方程,进一步推导得到了任意线性黏弹性模型的微分型本构方程的一般矩阵形式。论文最后总结并建立了一套适合计算机编程的固定范式,利用Python编程实现了该算法、获得了一些数值计算结果。     

A coupled temperature and strain rate dependent yield function for dynamic deformations of bcc metals

A coupled temperature and strain rate microstructure physically based yield function is proposed in this work. It is incorporated along with the Clausius–Duhem inequality and an appropriate free energy definition in a general thermodynamic framework for deriving a three-dimensional kinematical model for thermo-viscoplastic deformations of body centered cubic (bcc) metals. The evolution equations are expressed in terms of the material time derivatives of the elastic strain, accumulated plastic strain (isotropic hardening), and the back stress conjugate tensor (kinematic hardening). The viscoplastic multipliers are obtained using both the Consistency and Perzyna viscoplasticity models. The athermal yield function is employed instead of the static yield function in the case of the Perzyna viscoplasticity model. It is found that the static strain rate value, at which the material shows rate-independent behavior, varies with the material deformation temperature. Computational aspects of the proposed model are addressed through the finite element implementation with an implicit stress integration algorithm. Finite element simulations are performed by implementing the proposed viscoplasticity constitutive models in the commercial finite element program ABAQUS/Explicit [ABAQUS, 2003. User Manual, Version 6.3. Habbitt, Karlsson and Sorensen Inc., Providence, RI] via the user material subroutine coded as VUMAT. Numerical implementation for a simple compression problem meshed with one element is used to validate the proposed model implementation with applications to tantalum, niobium, and vanadium at low and high strain rates and temperatures. The analysis of a tensile shear banding is also investigated to show the effectiveness and the performance of the proposed framework in describing the strain localizations at high velocity impact. Results show mesh independency as a result of the viscoplastic regularization used in the proposed formulation.     

Modeling deformation behavior of polymers with viscoplasticity theory based on overstress

The nonlinear strain rate sensitivity, multiple creep and recovery behavior of polyphenylene oxide (PPO), which were explored through strain rate-controlled experiments at ambient temperature by Khan [The deformation behavior of solid polymers and modeling with the viscoplasticity theory based overstress, Ph.D. Thesis, Rensselaer Polytechnic Institute, New York], are modeled using the modified viscoplasticity theory based on overstress (VBO). In addition, VBO used by Krempl and Ho [An overstress model for solid polymer deformation behavior applied to Nylon 66, ASTM STP 1357, 2000, p. 118] and the classical VBO are used to demonstrate the improved modeling capabilities of VBO for solid polymer deformation. The unified model (VBO) has two tensor valued state variables, the equilibrium and kinematic stresses and two scalar valued states variables, drag and isotropic stresses. The simulations include monotonic loading and unloading at various strain rates, multiple creep and recovery at zero stress. Since creep behavior has been found to be profoundly influenced by the level of the stress, the tests are performed at different stresses above and below the yield point. Numerical results are compared to experimental data. It is shown that nonlinear rate sensitivity, nonlinear unloading, creep and recovery at zero stress can be reproduced using the modified viscoplasticity theory based on overstress.     

The parametric variational principle for Perzyna model in viscoplasticity

This paper presents the parametric variational principle for Perzyna model which isone of the main constitutive relations of viscoplasticity.The principle,by which the potentialenergy function is minimized under a constrained condition transformed by the constitutiverelations of viscoplasticity,is free from the bound of Drucker’s postulate of plastic flow andconsequently suitable for solving the nonassociated plastic flow problems.Furthermore,the paper has proven the presented principle and discussed the creep problem.     

The initiation and growth of adiabatic shear bands

A simple version of thermo/viscoplasticity theory is used to model the formation of adiabatic shear bands in high rate deformation of solids. The one dimensional shearing deformation of a finite slab is considered. For the constitutive assumptions made in this paper, homogeneous shearing produces a stress/strain response curve that always has a maximum when strain and rate hardening, plastic heating, and thermal softening are taken into account. Shear bands form if a perturbation is added to the homogeneous fields just before peak stress is obtained with these new fields being used as initial conditions. The resulting initial/boundary value problem is solved by the finite element method for one set of material parameters. The shear band grows slowly at first, then accelerates sharply, until finally the plastic strain rate in the center reaches a maximum, followed by a slow decline. Stress drops rapidly throughout the slab, and the central temperature increases rapidly as the peak in strain rate develops.     

Classes of flow rules for finite viscoplasticity

Classes of flow rules for finite viscoplasticity are defined by assuming that certain measures for plastic strain rate and plastic spin depend on the state variables but not on the plastic deformation. It is shown that three of these classes are mutually exclusive for finite elastic strains. For small elastic shear strains, two of the three classes are approximately equivalent. A number of exact and approximate kinematic relations between the various measures for plastic strain rate and plastic spin are derived. Some inconsistent flow rules encountered in the literature are also discussed. Throughout the paper, arbitrarily anisotropic materials are considered, and some of the simplifications resulting from the assumption of isotropy are noted.     

A thermo-viscoelastic–viscoplastic–viscodamage constitutive model for asphaltic materials

A temperature-dependent viscodamage model is proposed and coupled to the temperature-dependent Schapery’s nonlinear viscoelasticity and the temperature-dependent Perzyna’s viscoplasticity constitutive model presented in Abu Al-Rub et al., 2009, Huang et al., in press in order to model the nonlinear constitutive behavior of asphalt mixes. The thermo-viscodamage model is formulated to be a function of temperature, total effective strain, and the damage driving force which is expressed in terms of the stress invariants of the effective stress in the undamaged configuration. This expression for the damage force allows for the distinction between the influence of compression and extension loading conditions on damage nucleation and growth. A systematic procedure for obtaining the thermo-viscodamage model parameters using creep test data at different stress levels and different temperatures is presented. The recursive-iterative and radial return algorithms are used for the numerical implementation of the nonlinear viscoelasticity and viscoplasticity models, respectively, whereas the viscodamage model is implemented using the effective (undamaged) configuration concept. Numerical algorithms are implemented in the well-known finite element code Abaqus via the user material subroutine UMAT. The model is then calibrated and verified by comparing the model predictions with experimental data that include creep-recovery, creep, and uniaxial constant strain rate tests over a range of temperatures, stress levels, and strain rates. It is shown that the presented constitutive model is capable of predicting the nonlinear behavior of asphaltic mixes under different loading conditions.     

An experimental study of uniaxial creep,cyclic creep and relaxation of aisi type 304 stainless steel at room temperature

Following previous work (Krempl, 1979), a servocontrolled testing machine and strain measurement at the gage length were used to study the uniaxial rate(time)-dependent behavior of AISI Type 304 stainless steel at room temperature. The test results show that the creep strain accumulated in a given period of time depends strongly on the stress-rate preceding the creep test. In constant stress-rate zero-to-tension loading the creep strain accumulated in a fixed time-period at a given stress level is always higher during loading than during unloading. Continued cycling causes an exhaustion of creep ratchetting which depends on the stress-rate. Periods of creep and relaxation introduced during completely reversed plastic cycling show that the curved portions of the hysteresis loop exhibit most of the inelasticity. In the straight portions, creep and relaxation are small and there exists a region commencing after unloading where the behavior is similar to that at the origin for virgin materials. This region does not extend to zero stress.The results are at variance with creep theory and with viscoplasticity theories which assume that the yield surface expands with the stress. They support the theory of viscoplasticity based on total strain and overstress.     

THE OPTIMAL CONTROL VARIATIONAL PRINCIPLE AND FINITEELEMENTS ANALYSIS FOR VISCOPLASTIC DYNAMICS

I.Viscoplastic-DynamicsConstitutiveRelationThestructurematerialwouldarisesimultaneouslytheelasticity,viscosityandplasticityindynamichighstraincondition.Atpresent,theovcrstressmodel,viscoplasticitymodel,quasi-linearconstitutiverelationtheoryandintrinsicvar…