基于Weibull分布的盐岩分数阶蠕变本构模型研究

采用流变程控仪对湖北江汉油田盐岩试样进行了三轴蠕变实验.基于Weibull随机分布获得岩石损伤变量随加载时间的变化形式,构建了变黏性系数Abel粘壶本构关系.基于分数阶微分理论,采用将变黏性系数Abel粘壶替代西原正夫模型中Newton粘壶的方法,建立了盐岩蠕变损伤本构模型,并采用最小二乘法拟合获得了该本构模型中的参数.通过拟合曲线对比发现,论文提出的盐岩分数阶蠕变本构模型在描述盐岩的蠕变全过程,尤其是加速蠕变过程中具有较明显的优势,并且验证了西元正夫模型为论文所提出的本构模型的一个特例.     

一种改进的R-L分数阶导数定义在固体推进剂粘弹性本构模型中的应用

针对Riemann-Liouville分数阶导数定义的不足之处进行了改进,利用改进过的分数阶导数定义,建立了类Kelvin体粘弹性本构模型,并应用于某固体推进剂上,对本构方程中的三个参数进行了求解,与经典的prony级数模型进行比较,采用分数阶导数的类Kelvin体粘弹性本构模型与实验结果能很好地吻合.     

软岩的动态力学本构模型

利用分离式霍布金森压杆实验系统(SHPB)对2种典型的软岩砂质、泥岩进行动态力学性能测试,测试的应力应变曲线表现出显著的应变硬化和塑性流动等复杂的动态力学特性。基于实验结果,在朱-王-唐模型的基础上,并且考虑软岩本身结构缺陷的影响,建立了一种适应软岩材料的黏弹性统计损伤模型。该模型由2个Maxwell体和1个损伤体并联组成,2个Maxwell体用来描述对软岩高分子材料的高低应变率响应以及软岩材料在动载作用下的塑性变形特性,考虑软岩材料本身缺陷的影响,用一个损伤体代替朱-王-唐模型中的非弹性弹簧。利用黏弹性统计损伤模型对不同应变率下软岩材料的动态应力应变曲线进行拟合,拟合曲线与测试曲线一致性良好。     

结构整数阶模型的分数阶控制方法研究

为了把分数阶控制方法引入到控制策略中,假设原整数阶模型为分数阶模型,将其表示为状态方程形式;通过系统输出相同的方式建立分数阶和整数阶系统之间的联系;为保证分数阶系统控制器能够用于整数阶系统,采用了基于系统输出的控制策略。算例结果表明:在相同控制能耗条件下,与整数阶控制策略相比,采用本文方法能使位移降低30%。     

分数阶导数粘弹性模型的有限元法

在分析分数阶导数三元件模型理论的基础上,把分数阶导数三元件模型引入有限元模型中,推导出具有分数阶导数三元件本构关系的粘弹性结构动力学有限元格式。同时,应用分数阶导数型粘弹性结构动力学方程的数值算法求解了该有限元格式的数值解。并以二维沥青路面结构为例进行了路面动态粘弹性响应分析。算例分析表明,该方法能够正确有效地进行路面动态粘弹性分析。     

老化粘弹性蠕变本构方程中的热力学条件

针对具有"老化"特性的粘弹性材料,从不可逆热力学角度给出了其本构关系应该满足的条件;以混凝土材料为例,对某些蠕变关系式进行探讨,指出了老化材料蠕变方程应该满足的关系.     

基于扰动状态概念的软岩蠕变本构模型与试验验证

软岩材料的蠕变过程从初始加载阶段开始即容易产生塑性变形,这与硬岩存在较大差异,而传统的元件模型仅在应力超过某一阈值后才开始描述其粘塑性行为。通过分析岩石加卸载时效变形过程中弹塑性状态的演变规律,以扰动状态理论为基础,确定了以塑性变形为变量的扰动因子函数,并通过塑性变形随时间的变化特征进一步建立了以时间为自变量的蠕变扰动因子演化方程。在此基础上,针对泥质页岩的蠕变变形过程,选取Burgers和Bingham模型分别描述扰动状态理论中的相对完整状态和完全调整状态,并选取蠕变扰动因子为权重函数建立了基于扰动状态理论的蠕变本构模型。通过泥质页岩的室内蠕变特性试验对模型进行有效性检验表明,理论曲线与实测曲线的逼近程度较高,蠕变扰动函数能随时间的发展持续调整相对完整状态向完全调整状态转换的过程,较高应力状态下软岩进入完全调整状态的速率较快,相对传统模型该模型具有更好的非线性及阶段协调性,可较好地描述软岩的初始、稳定和加速蠕变阶段。     

岩石流变分数阶模型研究

以Kelvin流变模型为研究对象,提出了一种分数阶Kelvin流变模型。首先,把Kelvin模型中的整数阶导数改为分数阶导数,考虑到岩石材料的频率通常不超过1000 Hz,在分数阶拟合时,拟合频段选取为[0 1000],进而利用Oustalop滤波算法把分数阶表示为整数阶模式;其次,利用试验数据对分数阶模型进行参数识别,考虑到分数阶Kelvin模型具有强非线性的特点,引入了Levenberg-Marquardt优化算法来确定未知参数;最后,对于频域表示的流变方程,利用Laplace逆变换获得流变精确表达式。仿真实例表明本文方法可以很好地反映岩石流变特性。     

分数阶微分流变模型圆形硐室黏弹塑性分析

基于分数阶微分黏弹塑性流变模型,同时考虑受硐室掌子面影响的围岩应力释放效应,根据Laplace变换、分数阶微积分理论、Mittag-Leffler函数,推导了圆形硐室黏塑性区半径、应力、位移的理论解。将理论解与西原本构模型解进行对比,证明了分数阶模型解的合理性。分析结果表明:黏塑性区半径、洞壁环向应力的分数阶微分模型解和西原模型解随时间最终都趋于稳定,两种解法在稳定前有一定差异,稳定后大小相同;一定范围内,掌子面对围岩变形具有一定抑制作用,越靠近掌子面,黏塑性区半径和洞壁径向位移越小;流变特征对硐室围岩黏塑性区环向应力的影响较大,离硐室圆心越近,影响效果越明显,而黏塑性区径向应力和黏弹性区环向应力、径向应力受围岩流变的影响较小。     

人工冻结深部软岩单轴力学性能试验及蠕变模型

为了了解深部软岩在冻结条件下的单轴力学性能,以东北地区的原状泥砂岩为试验对象,利用自行研制的WDT-100型人工冻土试验仪器,对其进行不同温度下的人工冻土单轴抗压强度试验和单轴蠕变试验,得到泥砂岩单轴压缩应力-应变关系曲线,各温度下试样的单轴抗压强度以及蠕变曲线.单轴压缩试验结果表明：试样在给定温度和加载速率条件下,单轴压缩应力-应变关系曲线都有较为明显的屈服点,并且都在屈服点后,强度有所提高,出现硬化现象.单轴蠕变试验结果表明：单轴压缩蠕变曲线有非线性特征,单轴压缩蠕变的等时应力-应变曲线向应变轴靠拢;单轴压缩时蠕变模量随时间的增长而降低.最后采用遗传算法优化模型参数,得出泥砂岩蠕变经验方程.与试验结果对比,发现拟合情况较好.     

基于饱依丁-汤姆逊模型的软岩损伤蠕变模型

为更加准确地描述深部软岩的蠕变全过程,以饱依丁-汤姆逊模型为基础,将损伤变量引入到蠕变方程中,并构造了可描述加速蠕变阶段的模型元件,通过串联得到了新的蠕变模型。结果表明:该模型能够较好地描述软岩的衰减蠕变、稳态蠕变、加速蠕变阶段。用该模型对泥质页岩的蠕变试验结果进行数据拟合对比,结果显示该模型的拟合曲线与实验曲线基本吻合,拟合度为0.99733,能够很好地反映出软岩蠕变的特性,适合用于描述深部软岩的蠕变行为。     

有限应变下蠕变—弹塑性—损伤耦合本构模型的数值模拟

本文提出了一个模拟有限元应变下各向同性蠕变－弹塑性－损伤耦合本构行为的数值模型。模型基于不可逆过程的热动力学和内状态变量理论、以及应力就变本构关系的超弹性形式描述和变形梯度的弹性及非弹性部分的乘式分解。     

三轴加载下盐岩蠕变损伤特性的研究

基于不可逆热力学,采用内变量刻划岩石材料的不可逆变形历史,引入四阶损伤张量,建立了盐岩的蠕变损伤本构模型,对盐岩在复杂应力条件下的蠕变行为进行了描述,与实验结果比较表明该模型能较好的描述盐岩不同阶段的蠕变行为及其主要因素对盐岩蠕变特性的影响.     

三种分形和分数阶导数阻尼振动模型的比较研究

标准的整数阶导数方程不能准确描述粘弹性材料的记忆性参考文献[1]和阻尼的分数次幂频率依赖[2],因此分形导数、分数阶导数及正定分数阶导数被用于描述粘弹性介质中的阻尼振动.该文通过分析模型和数值模拟,比较了三种模型描述的振动过程.结果显示,当p小于约O.75或大于约1.9时(p为非整数阶导数的阶数),分形导数模型衰减最快;当P大于约0.75且小于约1.9时,正定分数阶导数模型衰减最快,衰减最慢的分别为分数阶导数模型(p1).且正定分数阶导数模型衰减快于分数阶导数模型,当p接近2时,两种模型较为相近.     

黏弹性材料等效分数阶微观结构标准线性固体模型

从黏弹性材料微观链结构出发,以橡胶基黏弹性材料超弹性理论分子网链高斯(Gauss)统计模型和黏滞流动理论为基础,研究黏弹性材料的微观结构、填料等对黏弹性性能的影响.用温频等效原理描述温度对黏弹性材料力学性能的影响,建立了可以有效描述黏弹性材料耗能特性的等效分数阶微观结构标准线性固体模型.采用动态热机械分析仪(DMA)对高聚物黏弹性材料力学性能、耗能能力进行测试.试验表明:在低温区域,储能模量较大,随着温度的升高,储能模量下降显著;能量损耗因子在高温和低温区域数值较小,在玻璃化转变温度附近数值较高.根据测试数据对所提等效分数阶微观结构标准线性固体模型进行验证,该力学模型能够较好地描述黏弹性材料储能模量和能量损耗因子随温度的变化趋势.用9050A和ZN22黏弹性材料对模型的有效性进一步验证,结果表明:9050A和ZN22黏弹性材料具有较好的耗能能力,所提出的等效分数阶微观结构标准线性固体模型能够准确地描述微观结构和填料对黏弹性材料宏观性能的影响,能够准确地描述黏弹性材料在不同温度和频率下的动态力学性能.     

丙烯酸弹性体的率相关分数阶黏弹性模型研究

丙烯酸弹性体VHB 4910作为一种重要的介电弹性体,在软体机器人、致动器、俘能器和智能隔振器等领域有很好的应用前景.但材料的非线性黏弹性对其力学行为有显著影响.近来分数阶模型在复杂材料的建模中取得了成功.本文基于分数阶有限变形Kelvin-Voigt流变学模型建立弹性体的三维张量本构,并进一步推导单向拉伸情况下的本构...     

网纹红土分数阶应力松弛模型

网纹红土松弛特性具有非线性。基于分数阶微积分理论,探讨了能更准确描述网纹红土松弛非线性特性和全过程的FVMS(Fractional Voigt and Maxwell model in series)松弛模型和FVMP(Fractional Voigt and Maxwell model in parallel)松弛模型及其理论解,进而应用提出的模型对三轴松弛试验实测数据进行反演,讨论了分数阶阶数的敏感性,并与西原模型和Burgers模型进行对比分析。研究结果表明,建立的四元件分数阶松弛本构模型应用于网纹红土应力松弛特性分析是有效可行的,模型灵活且精度更高,参数确定简便,发现分数阶阶数对应力松弛量的影响较大,但其对FVMS模型和FVMP模型松弛速率的影响不同,为实际工程长期稳定性分析提供了参考。     

人工冻结深部软岩单轴力学性能试验及蠕变模型1）

为了了解深部软岩在冻结条件下的单轴力学性能,以东北地区的原状泥砂岩为试验对象,利用自行研制的WDT-100型人工冻土试验仪器,对其进行不同温度下的人工冻土单轴抗压强度试验和单轴蠕变试验,得到泥砂岩单轴压缩应力-应变关系曲线,各温度下试样的单轴抗压强度以及蠕变曲线.单轴压缩试验结果表明：试样在给定温度和加载速率条件下,单轴压缩应力-应变关系曲线都有较为明显的屈服点,并且都在屈服点后,强度有所提高,出现硬化现象.单轴蠕变试验结果表明：单轴压缩蠕变曲线有非线性特征,单轴压缩蠕变的等时应力-应变曲线向应变轴靠拢;单轴压缩时蠕变模量随时间的增长而降低.最后采用遗传算法优化模型参数,得出泥砂岩蠕变经验方程.与试验结果对比,发现拟合情况较好.     

主应力轴旋转下黏土累积变形的分数阶元件模型

在长期交通载荷作用下土体塑性累积变形本构模型对路基沉降计算至为关键.元件组合模型可以计算岩土体循环累积应变,但现有的各类元件模型未能反映饱和软黏土的主应力轴循环旋转现象.在对饱和软黏土进行等向固结条件下的主应力轴循环旋转加载试验及非等向固结下的循环扭剪试验基础上,将Abel黏壶代替Burgers模型中的Newton黏壶,得到分数阶Burgers模型;利用遗传算法优化循环塑性累积应变的Burgers模型和分数阶Burgers模型的参数,通过对比两组模型的计算值与试验值,发现分数阶模型更适合模拟计算循环载荷下饱和软黏土的累积变形.     

基于分数阶模型的牡蛎壳动力学特性研究

贝壳、牡蛎等天然材料因其轻质高强的力学特性在材料设计等领域受到了广泛的关注,但由于材料本身结构的复杂性,对其力学行为的研究十分困难。近年来,分数阶模型在研究材料的力学特性上取得了成功,相比传统模型,分数阶模型可以更好地表征复杂介质的应力或应变与时间的关系。因此,本文从波传播理论出发,以分数阶模型作为材料本构,得到了复杂介质的波传播控制方程。通过Laplace变换得到了控制方程的解析解,并通过Laplace数值逆变换分析了波的衰减对分数阶模型中参量的敏感性,讨论了不同于材料弹性、黏性的材料“惯性”特性。接着,基于解析解和多种实验测试信号,给出了得到分数阶模型参数的拟合式子。以牡蛎材料作为研究对象,利用CO₂脉冲激光器进行小试样的冲击加载、应用两点激光干涉测速系统(laser interferometer velocimetry system, VISAR)对表面粒子的速度进行测量,得到了4种密度下不同厚度的牡蛎壳试样的粒子速度时程曲线,再结合上述理论方法分析得到了牡蛎壳试样的Abel模型和分数阶Maxwell模型的参数,模型参数反映了牡蛎壳试样的细微观结构特征。结果...