分级加载条件下紫色泥岩三轴蠕变特性研究

流变性质是岩石材料的重要力学特性,以紫色泥岩为例,在分级加载条件下研究了其三轴蠕变特性。采用自行研制的重力杠杆式岩石蠕变实验机,并配备三轴压力室,对紫色泥岩实施分级加载三轴蠕变实验。实验结果表明,紫色泥岩在相同轴压条件下,围压越大,瞬时弹性变形和蠕变变形越小;在相同围压条件下,轴压越大,二者变形越大。紫色泥岩的蠕变速率具有相同的变化趋势,且在相同轴压条件下,围压越大,稳定蠕变阶段的蠕变速率越不明显。另外,应力与瞬时应变线性拟合结果表明,瞬时轴向应变与轴向荷载近似成比例增长。通过分析紫色泥岩蠕变曲线变化趋势,认为采用H/M模型进行模拟较为合理。并且采用MATLAB软件非线性回归分析进行模型参数拟合,最后得到实验曲线与理论曲线的最大误差在2e-5左右。证明采用的H/M模型能较好地描述紫色泥岩的蠕变特性。     

岩石蠕变特性试验研究

针对大安山煤矿+550水平西二石门围岩蠕变变形严重的问题,采用室内蠕变试验与理论分析相结合的方法研究岩石蠕变特性。选取泥质粉砂岩作为软岩试样,通过自制三轴蠕变仪对所取岩样进行恒围压分级增轴压的三轴压缩蠕变试验,结果表明:随应力水平的提高,起始蠕变速率增大,进入稳态蠕变阶段用时延长,稳态蠕变速率从0h-1增大到3.12e-4·h-1,瞬时应变以线性关系增加,围压9MPa条件下的起始蠕变应力阈值为6MPa,当应力水平超过岩石长期强度时出现加速蠕变。基于蠕变试验规律,考虑起始蠕变应力阈值,将一种可用阶跃函数表示的开关元件与广义K体并联,并引入损伤元件对CVISC蠕变模型进行改进。采用1st-Opt中的Levenberg-Marquardt+通用全局优化法计算蠕变参数并进行反演。反演结果表明理论曲线与试验曲线的吻合效果较好,证明了所建模型的合理性。上述研究成果可为岩石工程的长期稳定性分析提供试验和理论依据。     

界面特性对短纤维金属基复合材料蠕变行为的影响

基于短纤维增强金属基复合材料(MMC)的单纤维三维模型(三相)，利用粘弹性有限元分析方法对影响金属基复合材料的蠕变行为的因素进行了较为系统的分析。研究中主要讨论了界面特性和纤维取向角对金属基复合材料的蠕变性能的影响。研究结果发现，界面特性诸如厚度、模量和应力指数都对纤维最大轴应力和稳定蠕变率产生影响：稳态蠕变率随界面模量的增大而逐渐减小，当高于基体模量时基本保持不变；纤维轴应力的变化与蠕变率正好相反。稳态蠕变率随界面厚度、应力指数的增加而增大；而轴应力则随之减小。同时不同的纤维取向也影响金属基复合材料蠕变时的轴应力分布和稳态蠕变率。     

砂岩短时单轴直接拉伸蠕变特性试验研究

采用自行设计加工的挂重型岩石材料直接拉伸装置,对重庆市某地红砂岩进行了短时分级加载单轴直接拉伸蠕变试验,并对该砂岩的蠕变特性进行了分析。结果表明,砂岩在该试验条件下具有明显的蠕变特性,其蠕变表现为衰减蠕变和稳态蠕变两个阶段,砂岩蠕变量与稳态蠕变速率随着荷载的增加而增加。应用三参量广义Kelvin模型与Burgers蠕变模型描述该砂岩的蠕变规律,结果表明,两个模型均可以准确描述该砂岩的蠕变特性,其中,Burgers蠕变模型效果更佳。本文所得结论可为从事相关研究的工作人员提供一定的借鉴。     

结构面泥岩非线性蠕变特性试验研究

利用自制蠕变试验装置对含有平直闭合结构面的泥岩试样进行单轴蠕变试验研究。结果表明:轴向应力为8~12MPa时,结构面泥岩蠕应变与总应变比值呈线性增加;随着轴向应力增加,蠕应变与总应变比值突变呈上凹型增加,继而呈上凸型变化;结构面泥岩蠕变速率随轴向应力增加衰减显著,与时间为负指数关系。为描述不同轴向应力下结构面泥岩蠕变规律,提出一种非线性元件与时间呈负指数关系的非线性黏滞体,构建结构面泥岩法向非线性蠕变模型,并进行参数辨识及拟合分析。分析发现改进的Burgers模型能很好地描述结构面泥岩的法向非线性蠕变特性。     

含单裂隙岩体蠕变破裂裂纹扩展演化规律研究

岩体内部存在大量不同尺寸裂隙，形成不同结构形式的裂隙岩体，裂隙的存在会降低岩体的力学性能，对于岩体工程的整体稳定性具有很大的影响。本研究基于滑动翼裂纹模型对恒定载荷下含单裂隙岩体蠕变变形进行了理论推导和分析，基于MAPLE软件直观给出恒定载荷下裂隙角度、长度对蠕变变形的影响，采用连续-非连续数值分析软件GDEM,对预制单裂隙在不同倾角、不同长度情况下岩体的蠕变破裂裂纹扩展进行数值模拟。通过理论分析和数值模拟可以看出，岩体蠕变速率随着裂隙长度的增大而增大；随着裂隙角度的增加出现先增大后减小的趋势，其中裂隙倾角为30°时蠕变速率最大。含单裂隙岩体蠕变破裂裂纹扩展规律的研究可为岩体工程蠕变破裂规律研究及蠕变治理提供一定的理论依据和参考。     

稳定阶段非线性蠕变的不可逆热力学解释

文中首先介绍了作者提出的反耗散概念和由之引出的耗散条件，在此基础上把材料参数当作遵守一定分布的随机变量，建立了稳定阶段蠕变的概率本构模型。它包含了线性模型和幂织经模型等，并能解释变速率随应力增加而增加的现象。     

蠕变压痕实验的计算机模拟

采用有限元的方法对双相材料的蠕变压痕实验进行了数值模拟，在有限元数值解的基础上，定义了相应于传统单轴蠕变实验的“等效应力”和“等效应变”，建立了蠕变压痕实验同传统单轴蠕变实验之间的关系，给出了确定薄膜蠕变应力指数和蠕变常数的方法；同时数值解的结果表明，实验中通过控制压痕深度不超过薄膜厚度的 ５％～１０％，忽略基体的硬化指数对确定薄膜性能的影响存在一定的误差，但基体的弹性模量对确定薄膜的蠕变性能影响不大。     

砂岩峰后蠕变失稳特性试验研究

为揭示岩石峰后蠕变失稳特性,利用MTS815电液伺服岩石力学试验系统,对红砂岩进行了峰后蠕变试验,采用线性回归的方法估算了蠕变失稳时间,探讨了损伤、初始围压以及围压增量对峰后蠕变特性的影响规律。结果表明:随着损伤的加大,蠕变变形速率增大,蠕变失稳时间缩短。同一围压下,各级应力水平下蠕变量随应力水平的增大而增大,且服从指数函数关系。增大初始围压,各级应力水平下的蠕变变形速率均有一定程度的降低,且对等速蠕变阶段蠕变速率的影响较加速蠕变阶段蠕变速率大,蠕变失稳时间延长。施加围压增量,大大减小了蠕变变形速率,促使了加速蠕变阶段向等速蠕变阶段转化,明显延长了蠕变失稳时间,减小了蠕变形变量。研究结果为工程围岩的长期稳定性分析提供了参考依据。     

蠕变应变的予测

蠕变分析仅确定稳态解,由于多数低合金钢进入加速阶段较早,结果与实际情况差别较大,需研究能更精确予测蠕变应变的方法.文中介绍了制订新方法的途径和算法,推荐了一种带应力函数的蠕变速率综合方程,可满意地予测蠕变各阶段的应变,并举例说明它在蠕变分析中的应用.     

Non-linear creep damage under one-dimensional variable tensile stress

A non-linear damage relation, containing the axial strain history and a time integral over the stress history, is proposed for the case of one-dimensional time dependent tensile stress. Non-linear steady and transient creep terms are included in the axial strain relation, and elastic and creep Poisson's ratios are introduced into the lateral strain relation. Using these relations, complete damage solutions are obtained for the constant stress rate, step stress, relaxation and constant load tests. Observations are made concerning the associated rupture times.     

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.     

ANALYSIS ON PSEUDO-STEADY INDENTATION CREEP

Theoretical analysis and finite element （FE） simulation have been carried out for a constant specific load rate （CSLR） indentation creep test. Analytical results indicate that both the representative stress and the indentation strain rate become constant after a transient period. Moreover, the FE simulation reveals that both the contours of equivalent stress and equivalent plastic strain rate underneath the indenter evolve with geometrical self-similarity. This suggests that pseudo-steady indentation creep occurs in the region beneath the indenter. The representative points in the region are defined as the ones with the equivalent stress equal to the representative stress. In addition, it is revealed that the proportionality between indentation strain rate and equivalent plastic strain rate holds at the representative points during the pseudo-steady indentation creep of a power law material. A control volume （CV） beneath the indenter, which governs the indenter velocity, is identified. The size of the CV at the indented surface is approximately 2.5 times the size of the impression. The stress exponent for creep can be obtained from the pseudosteady indentation creep data. These results demonstrate that the CSLR testing technique can be used to evaluate creep parameters with the same accuracy as conventional uniaxial creep tests.     

An analytical and computational study of crack initiation under transient creep conditions

An analytical method has been developed to predict creep crack initiation (CCI), based on the accumulation of a critical level of damage at a critical distance. The method accounts for the re-distribution of stress from the elastic or elastic–plastic field, experienced on initial loading, to a steady state creep stress distribution, via a transient creep region. The method has been applied to predict CCI times in a fracture specimen of type 316H stainless steel at 550 °C. The failure model has been also been implemented into a finite element (FE) framework. Reasonable and conservative predictions of CCI time can be obtained from the analytical solution relative to FE solutions. Conservative predictions of experimental CCI times are obtained when stress redistribution is taking into account. However, CCI times predicted from a steady state creep model are found to be non-conservative.     

应用于高速轮轨滚动接触的蠕滑理论算法对比研究

运用多体动力学仿真软件UM建立了车辆-轨道多体动力学模型，并基于轮轨滚动蠕滑理论对比分析了FASTSIM算法与CONTACT算法两者在稳态和瞬态滚动状态中轮轨动态作用力之间的差异.计算结果表明:列车以稳态滚动状态行驶在直线轨道和曲线轨道上时通过FASTSIM算法计算得出的横、纵向轮轨蠕滑力与CONTACT算法的计算值存在的差异值较小，达到了10%~15%，而在瞬态滚动状态下尤其是在大蠕滑以及考虑风阻等比较恶劣的工况下行驶时计算得出的横、纵向蠕滑力与CONTACT算法的计算值存在较大的误差，差值达到了25%，并且随着车速的增加，最大差异值达到了30%，这表明CONTACT算法在瞬态滚动工况下进行蠕滑力计算上相比于FASTSIM算法更加精准，更适用于作出对列车运行时的安全性方面的评价.      

炭质泥岩泥化夹层的流变特性及长期强度

通过分析炭质泥岩泥化夹层的直剪流变实验成果,表明炭质泥岩泥化夹层的法向应力和剪应力为线性关系,存在统一的蠕变加速阶段来临时的极限应变量。基于老化理论,建立了适合弱层流变特性的流变力学模型,该模型可以描述初始蠕变、等速蠕变和加速蠕变三个阶段过程。依据蠕变加速阶段来临时的极限应变量得到了长期强度方程,从而为边坡变形分析和滑坡中长期预报提供了可靠的依据。     

Effect of constraint induced by crack depth on creep crack-tip stress field in CT specimens

In this paper, the effect of constraint induced by the crack depth on creep crack-tip stress field in compact tension (CT) specimens is examined by finite element analysis, and the effect of creep deformation and damage on the Hutchinson–Rice–Rosengren (HRR) singularity stress field are discussed. The results show the constraint induced by crack depth causes the difference in crack-tip opening stress distributions between the specimens with different crack depth at the same C*. The maximum opening stress appears at a distance from crack tips, and the stress singularity near the crack tips does not exist due to the crack-tip blunting caused by the large creep deformation in the vicinity of the crack tips. The actual stress calculated by the finite element method (FEM) in front of crack tip is significantly lower than that predicted by the HRR field. Based on the reference stress field in the deep crack CT specimen with high constraint, a new constraint parameter R is defined and the constraint effect in the shallow crack specimen is examined at different distances ahead of the crack tip from transient to steady-state creep conditions. During the early stages of creep constraint increases with time, and then approaches a steady state value as time increases. With increasing the distance from crack tips and applied load, the negative R increases and the constraint decreases.     

STUDY OF CREEP PARAMETERS OF GRANITE IN STEADY STATE UNDER CONDITIONS OF HIGH TEMPERATURE AND HYDROSTATIC PRESSURE

高温高压花岗岩钻孔实验表明,温度低于500℃,静水应力低于150MPa 状态下,岩体的钻孔变形均属于稳态蠕变变形阶段. 该文选取了广义开尔文模型来反映其特征,通过拉普拉斯变换及逆变换,详细推演出了钻孔径向位移解析解,并且考虑温度-应力的耦合效应,给出了模型参数随温度及应力变化的关系式. 利用该关系进行拟合计算,说明广义开尔文模型来表达高温高压环境中的花岗岩稳态蠕变变形特性,寻求蠕变参数是合理可靠的. 该文对于高温岩体地热资源开发中的钻孔施工与维护、钻孔变形量预测等方面,具有实际的指导意义.     

Recoverable creep deformation and transient local stress concentration due to heterogeneous grain-boundary diffusion and sliding in polycrystalline solids

Numerical simulations are used to investigate the influence of heterogeneity in grain-boundary diffusivity and sliding resistance on the creep response of a polycrystal. We model a polycrystal as a two-dimensional assembly of elastic grains, separated by sharp grain boundaries. The crystal deforms plastically by stress driven mass transport along the grain boundaries, together with grain-boundary sliding. Heterogeneity is idealized by assigning each grain boundary one of two possible values of diffusivity and sliding viscosity. We compute steady state and transient creep rates as functions of the diffusivity mismatch and relative fractions of grain boundaries with fast and slow diffusion. In addition, our results show that under transient conditions, flux divergences develop at the intersection between grain boundaries with fast and slow diffusivity, which generate high local stress concentrations. The stress concentrations develop at a rate determined by the fast diffusion coefficient, and subsequently relax at a rate determined by the slow diffusion coefficient. The influence of the mismatch in diffusion coefficient, loading conditions, and material properties on the magnitude of this stress concentration is investigated in detail using a simple model problem with a planar grain boundary. The strain energy associated with these stress concentrations also makes a small fraction of the plastic strain due to diffusion and sliding recoverable on unloading. We discuss the implications of these results for conventional polycrystalline solids at high temperatures and for nanostructured materials where grain-boundary diffusion becomes one of the primary inelastic deformation mechanisms even at room temperature.