考虑损伤的内变量黏弹-黏塑性本构方程

基于Rice 不可逆内变量热力学框架,在约束构型空间中讨论材料的蠕变损伤问题. 通过给定具体的余能密度函数和内变量演化方程推导出考虑损伤的内变量黏弹-黏塑性本构方程. 通过模型相似材料单轴蠕变加卸载试验对一维情况下的本构方程进行参数辨识和模型验证,本构方程能很好地描述黏弹性变形和各蠕变阶段.不同的蠕变阶段具有不同的能量耗散特点. 受应力扰动后,不考虑损伤的材料系统能自发趋于热力学平衡态或稳定态. 在考虑损伤的整个蠕变过程中,材料系统先趋于平衡态再背离平衡态发展. 能量耗散率可作为材料系统热力学状态偏离平衡态的测度;能量耗散率的时间导数可用于表征系统的演化趋势;两者的域内积分值可作为结构长期稳定性的评价指标.      

考虑损伤的内变量黏弹--黏塑性本构方程

基于Rice不可逆内变量热力学框架,在约束构型空间中讨论材料的蠕变损伤问题.通过给定具体的余能密度函数和内变量演化方程推导出考虑损伤的内变量黏弹--黏塑性本构方程.通过模型相似材料单轴蠕变加卸载试验对一维情况下的本构方程进行参数辨识和模型验证,本构方程能很好地描述黏弹性变形和各蠕变阶段.不同的蠕变阶段具有不同的能量耗散特点.受应力扰动后,不考虑损伤的材料系统能自发趋于热力学平衡态或稳定态.在考虑损伤的整个蠕变过程中,材料系统先趋于平衡态再背离平衡态发展.能量耗散率可作为材料系统热力学状态偏离平衡态的测度;能量耗散率的时间导数可用于表征系统的演化趋势;两者的域内积分值可作为结构长期稳定性的评价指标.     

黏弹流体流动的数值模拟研究进展

综述了黏弹流体流动数值模拟的研究进展,突出介绍近十年来有限元法在黏弹流体流动数值模拟研究中取得的成果,通过动量方程的适当变形和本构方程离散权函数的合理选择,可以显著增强数值计算的稳定性。得到较高Weissenberg数下的解,同时文中对黏弹流体流动数值模拟中本构方程的应用、非等温情况和三维空间下的研究进行了介绍。     

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

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

一个反映土体复杂加卸荷路径的弹黏塑性模型

基于Perzyna 过应力理论,构造了一个土体各向异性弹黏塑性本构模型. 模型引入Wheeler 旋转硬化法则,能够较好的描述土体的初始各向异性和应力诱发各向异性. 借助ABAQUS 软件中UMAT 子程序接口,模型采用隐式积分算法——图形返回算法实现. 通过对一组复杂加卸荷路径的三轴不排水剪切试验(HKMD) 及一组分级加载三轴不排水蠕变试验(Sachville clay) 的模拟,表明本模型能够合理反映土体的率效应、应力松弛、蠕变及土体各向异性现象,验证了模型的有效性.     

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

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

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

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

用于设计灵敏度分析的弹粘塑性边界元法

基于一致切线算子概念的弹粘塑性隐式边界元方法,进行了弹粘塑性设计灵敏度分析.采用了Perzyna经典粘塑性本构模型,针对包含各向同性硬化和运动硬化的混合硬化模型,导出了弹粘塑性灵敏度分析的径向返回算法和一致切线算子.利用直接微分的方法,建立了设计灵敏度分析的弹粘塑性边界元增量方程,导出了弹粘塑性径向返回的灵敏度公式.给出了在不同粘塑性流动参数下的三个典型算例的结果,与ANSYS结果相吻合,证明了方法是正确的.     

基于黏弹性的微注射成型流动模拟

应用有限元方法研究了微注射成型中瞬态、可压缩、非牛顿熔体流动的黏弹性对流动前沿及流动平衡的影响。基于Phan-Thien-Tanner模型建立了熔体流动的本构方程,利用Hele-Shaw假设和简化建立了瞬态、可压缩、非牛顿熔体流动的连续性方程、动量方程、能量方程;为了有效地描述微注射成型的尺寸效应,采用了边界滑移和表面张力边界条件。通过分部积分和待定系数法导出了带有边界信息的变分方程和求解应力分量的半解析公式,构造了有限元离散求解及超松驰迭代算法。模拟结果表明:熔体的黏弹性对浇口附近的压力和后续的熔体流动前沿有重要影响;与黏性模型相比,黏弹性模型可以控制模拟压力的快速增长,减少不同型腔之间的充填差异,与短射实验结果也更吻合。     

聊聊动态塑性和黏塑性

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

SnPb钎料合金的粘塑性Anand本构方程

采用统一型粘塑性本构 Anand方程描述了电子封装焊点 Sn Pb钎料合金的非弹性变形行为 ,基于 Sn Pb 合金的弹塑性蠕变本构方程和实验数据 ,确定了6 2 Sn36 Pb2 Ag、6 0 Sn40 Pb、96 .5 Sn3.5 Ag和 97.5 Pb2 .5 Sn四种钎料合金 Anand方程的材料参数 ,验证了粘塑性 Anand本构方程对 Sn Pb合金在恒应变速率和稳态塑性流动条件下应力应变行为的预测能力。结果表明 ,Anand方程能有效描述 Sn Pb钎料的粘塑性本构行为 ,并可应用于电子封装 Sn Pb焊点的可靠性模拟和失效分析     

Some basic solutions for nonlinear creep

The aim of the paper is to derive the exact analytical expressions for torsion and bending creep of rods that obey the Norton–Bailey, Prandtl–Garofalo and Naumenko–Altenbach–Gorash constitutive models. The common secondary creep constitutive model is the Norton–Bailey law which gives a power law relationship between creep rate and stress. The closed form solutions for fractional Norton–Bailey creep law are derived. The analytical formulas express the torque and bending moment as functions of the time for the period of relaxation. Other formulas express the twist rate and curvature as functions of the time for the duration of engineering creep experiment. The derived formulas are suitable for the practically important problems of machinery. Namely, the formulas are relevant for calculation of hereditary effects for helical, leaf and disk springs and twisted shafts.     

Finite deformation of a polymer: experiments and modeling

The response of a polymer (polytetrafluoroethylene) to quasi-static and dynamic loading is determined and modeled. The polytetrafluoroethylene is extremely ductile and highly nonlinear in elastic as well as plastic behaviors including elastic unloading. Constitutive model developed earlier by Khan, Huang and Liang (KHL) is extended to include the responses of polymeric materials. The strain rate hardening, creep, and relaxation behaviors of polytetrafluoroethylene were determined through extensive experimental study. Based on the observation that both viscoelastic and viscoplastic deformation of polytetrafluoroethylene are time dependent and nonlinear, a phenomenalogical viscoelasto–plastic constitutive model is presented by a series connection of a viscoelastic deformation module (represented by three elements standard solid spring dashpot model), and a viscoplastic deformation module represented by KHL model. The KHL module is affected only when the stress exceeds the initial yield stress. The comparison between the predictions from the extended model and experimental data for uniaxial static and dynamic compression, creep and relaxation demonstrate that the proposed constitutive model is able to represent the observed time dependent mechanical behavior of polytetrafluoroethylene polytetrafluoroethylene qualitatively and quantitatively.     

Experimental stress analysis for process of nonsteady inelastic deformation by photoviscoplasticity

The proposed photoviscoplasticity is a method for determining the stress distribution in nonsteady inelastic deformation during creep. The governing fundamental relations are derived by considering the effect of time as involved in the viscoplastic strain rate and the stress rate, and by considering the distinctive deformation properties of celluloid as model material. The validity of the fundamental relations are evaluated according to calibration test by using celluloid. The proposed method is applied for two practical applications of the compression of blocks by elastic punches and the compression of strip having a circular hole or semicircular notches. The time-dependent variation of stress states can be analyzed with the same accuracy as in the photoelasticity.     

Transient and steady-state nanoindentation creep of polymeric materials

The transient and steady-state nanoindentation creep of polymeric materials was investigated. The creep model is used to explain the experimental data of transient and steady-state creep dominated by viscoelastic deformation and power-law creep deformation, respectively. The Burgers viscoelastic model was used to interpret the transient creep in polymers under nano-indentation. Explicit expression for the displacement of transient creep was derived using the correspondence principle of linear viscoelasticity theory. The power law of strain rate-stress relation was used to explain the creep displacement during the steady state. Three polymers of poly(methyl methacrylate), hydroxyethyl methacrylate copolymer, and the fast-cure acrylic resin were used to measure the nanoindentation creep. The transient creep data are in good agreement with the predictions from the Burgers viscoelastic model. The creep displacement is mainly attributed to the viscous flow of the Kelvin element, and the computed values of viscosities (η_1,c, η_2,c) increase with decreasing preloading rate. By comparing the steady-state creep data with the power law of strain rate-stress relation, the stress exponents for the above polymeric materials were quantitatively determined.     

Explicit and implicit viscoplastic models for polymeric composite

The viscoplastic behavior of a carbon-fiber/polymer matrix composite was investigated through two different modeling efforts. The first model is phenomenological in nature and utilizes the tensile and stress relaxation experiments to predict the creep strain. The phenomenological model was constructed based on the overstress viscoplastic model. In the second model, the composite viscoplastic behavior is captured via neural networks formulation. The neural networks model was constructed directly from the experimental results obtained via creep tests performed at various stress–temperature conditions. The neural network was trained to predict the creep strain based on the stress–temperature–time values. The performance of the neural model is evaluated through the mean squared error between the neural network prediction and the experimental creep strain results. To minimize this error, several optimization techniques were examined. The minimization of the error when carried out by the Truncated Newton method outperforms the standard back-propagation and the conjugate gradient method in terms of convergence rate and accuracy. Using neural network with truncated Newton training algorithm, the prediction of the creep strain was very satisfactory compared to the phenomenological model.     

The viscoelastic and viscoplastic behavior of low-density polyethylene

Three series of tensile tests with constant cross-head speeds (ranging from 5 to 200 mm/min), tensile relaxation tests (at strains from 0.03 to 0.09) and tensile creep tests (at stresses from 2.0 to 6.0 MPa) are performed on low-density polyethylene at room temperature. Constitutive equations are derived for the time-dependent response of semicrystalline polymers at isothermal deformation with small strains. A polymer is treated as an equivalent heterogeneous network of chains bridged by temporary junctions (entanglements, physical cross-links and lamellar blocks). The network is thought of as an ensemble of meso-regions linked with each other. The viscoelastic behavior of a polymer is modelled as thermally-induced rearrangement of strands (separation of active strands from temporary junctions and merging of dangling strands with the network). The viscoplastic response reflects mutual displacement of meso-domains driven by macro-strains. Stress–strain relations for uniaxial deformation are developed by using the laws of thermodynamics. The governing equations involve five material constants that are found by fitting the observations. Fair agreement is demonstrated between the experimental data and the results of numerical simulation. It is shown that observations in conventional creep tests reflect not only the viscoelastic, but also the viscoplastic behavior of an ensemble of meso-regions.     

高应力高水压下砂岩三轴蠕变特性试验研究

对砂岩进行高围压高水压条件下的三轴压缩蠕变试验。试验表明,在整个加载过程中,孔隙水压力主要起到增强轴向变形和横向变形的作用;但在加载的初始阶段,孔隙水压力在一定程度上抑制了轴向变形。当应力水平高于屈服应力时,横向蠕变速率明显大于轴向蠕变速率,且横向蠕变率先进入加速蠕变阶段。本文提出一个新的非线性黏性元件,并引入一个能判断是否进入加速蠕变阶段的计时器元件,组建一个非线性黏塑性加速蠕变启动模型,将该黏塑性模型与Burgers模型串联,构建一个新的非线性黏弹塑性蠕变模型,推导了该模型在常规三轴应力状态下的本构方程。基于试验结果,通过对非线性优化算法(BFGS)的Matlab编程,实现对本文提出蠕变模型的参数识别,识别效果比较理想。对比试验曲线与拟合曲线,二者相当吻合,验证了新提出的非线性黏弹塑性蠕变模型的正确性。