改进宾汉姆流变模型及其参数辨识

传统宾汉姆流变模型不适于描述岩石蠕变的非线性衰减和加速阶段.本文提出一种非线性的改进模型,引入作用函数表征衰减蠕变阶段,引入非线性黏性元件与传统模型串联表征岩石的加速蠕变阶段.通过对绿片岩和砂岩蠕变实验数据进行非线性拟合得到模型参数,与实验曲线的对比分析表明了改进的效果及改进模型的合理性.     

基于西原模型的非线性黏弹塑性流变模型

在西原模型的基础上,引入能反映加速蠕变起始点的非线性黏塑性元件和瞬时塑性变形的塑性元件,建立了一个能完整反映蠕变全过程的非线性黏弹塑性流变模型. 在其一维流变模型结构的基础上建立三维的流变模型,推导三维流变方程,接着探讨了在三轴压缩试验中确定三维蠕变参数的方法,并对参数进行了敏感性分析. 用三轴压缩试验结果验证模型的有效性,结果表明模型能很好地模拟蠕变的三阶段,并且在多级偏应力水平下可以采用同一组蠕变参数.     

黄土损伤与流变耦合模型及参数研究

以试验为基础,对不同含水量Q3黄土的蠕变特征进行了分析研究,结合黄土结构损伤特性及其对黄土力学特征影响,通过考虑结构损伤对变形模量的劣化来反映黄土结构性对蠕变规律的影响,提出了黄土损伤与流变耦合作用的本构模型。该模型不仅具有参数少的优点,而且能够准确描述黄土的突发性破坏特征。利用结构性软土流变的研究成果,给出了黄土流变模型参数的确定方法,并对模型参数随非饱和黄土含水量变化的规律进行了研究。数值模拟结果表明,该模型能够很好地描述黄土的加速蠕变阶段的变形特征。     

岩石三轴蠕变试验黏弹性解析及参数识别

基于岩石三轴蠕变试验,推导了相应的黏弹性解析解,并对三参量H-K流变模型给出了黏弹性解析解的显式以及流变参数的识别方法. 在此基础上,针对花岗岩及片岩两种岩石的室内蠕变分级加载试验曲线进行了流变参数的识别与反演. 分析表明,同一岩样在不同加载级别下的弹性及黏弹性参数在量值上均有所不同,其中黏弹性模量和黏滞系数的量值随载荷的增加出现损伤衰减,基于试验反演所获参数的解析蠕变曲线与试验蠕变曲线吻合较好.     

改进的西原模型及其在冻结壁上的应用

通过对现有流变模型的分析及结合室内冻黏土大量的单轴蠕变试验,试验结果表明可以忽略冻土的瞬时变形并采用简单的屈服函数对线性牛顿体进行简化,从而对西原模型进行了改进,并导出了改进后模型在三向应力条件下冻土的非线性蠕变方程．采用改进的西原模型推导了冻结壁应力场和位移场的理论计算公式,探讨了防止冻结管断裂的有效措施．计算和工程实例基本相符,可以应用于工程计算．     

含分数阶Bingham模型的阻尼减振系统时滞半主动控制

针对基于磁流变液阻尼器的半主动控制系统中存在的时滞问题, 采用了一种将可控的时滞变量引入半主动控制切换条件的控制策略, 研究了考虑时滞的天棚阻尼控制切换条件对半主动阻尼减振系统的影响, 分析了含有分数阶Bingham模型的线性刚度系统在基础激励下的振动特性. 利用平均法得到了系统在含时滞半主动控制策略下主共振响应的近似解析解, 根据Lyapunov理论分析了系统的稳定性. 通过数值解验证了近似解析解的准确性, 二者具有较好的一致性. 利用近似解析解分析了固定激励频率下时滞对系统幅频响应特性的影响, 以及主共振峰值响应和共振频率随时滞变化的特性规律. 结果表明, 含时滞的半主动控制系统存在一个小时滞区间, 使得系统的振幅在主共振峰对应的频率附近低于不考虑时滞时系统的振幅, 且存在最优时滞使得系统的振幅大幅度降低; 而大时滞的引入会加剧系统的振动, 导致系统的颤振. 确定了基于分数阶Bingham模型的线性刚度系统在天棚阻尼半主动控制下的时滞选取原则, 为振动系统半主动阻尼控制中的时滞选取提供了参考.      

一种带梯度惩罚项的离散力学最优控制方法

离散力学最优控制方法(discrete mechanics and optimal control, DMOC) 是最新提出的一种有效的优化控制模型. 由于原始方法在求解优化问题时容易产生控制力的震荡,不利于控制力的稳定输出. 该文首先对产生这种数值现象的原因进行了简单的分析,并提出了可行的解决路径——加入惩罚梯度项. 最后采用数值算例对该方法进行验证,数值算例表明:改进的方法可以在保证力量最小化的前提下,使得控制力连续平滑.     

基于改进西原模型巷道围岩流变参数反演分析

当应力达到一定水平时,岩石经过衰减和稳定蠕变过程之后会发生加速蠕变破坏,但是传统西原体模型无法更有效地反映该加速蠕变阶段。为了解决这一问题,引入一个非线性流变元件,并将其串联在西原体模型上,组成一个新的六元件模型,该模型可充分反映岩石蠕变的三个阶段。基于改进的西原体模型本构方程,应用Laplace变换,得到了盾构圆形巷道围岩的径向位移变化规律;根据巷道围岩收敛变形现场实测数据,反演得到了岩石黏弹塑性蠕变参数;由巷道围岩径向位移计算结果和实测数据对比分析表明:两者吻合较好;当t ≤25 d时,巷道围岩处于加速变形阶段,其径向位移变化速率随时间的增长而不断增加,当t>25 d时,巷道围岩进入稳定变形阶段,随时间的推移,其径向围岩位移变化速率逐渐趋于稳定。     

冻结黏土力学特性分析及改进S-M蠕变模型1)

为合理描述人工冻结法施工的矿井土层力学特性,对某矿井土样的重塑土进行单轴抗压强度实验和单轴压缩蠕变实验,得到重塑冻结黏土在不同温度及载荷加载等级下的应力-应变曲线和蠕变曲线。实验结果表明：冻结温度越低,土样的单轴强度越大;相同冻结温度下,土样的蠕变变形随着应力水平的升高而增大。单轴压缩蠕变的等时应力-应变曲线随时间发展向应变轴靠拢;土样经历初始蠕变和等速蠕变两个阶段,在较高应力水平下有进入加速蠕变的趋势;对S-M模型中各参数的意义进行修正并考虑温度的影响,得到人工冻结黏土改进S-M蠕变显式模型,然后采用粒子群算法对人工冻结黏土S-M蠕变模型参数进行优化。改进S-M蠕变显式模型理论计算值与实验值吻合良好,表明改进S-M蠕变显式模型能较好模拟人工冻结黏土的蠕变特性。改进S-M蠕变显式模型为人工冻土蠕变计算提供一种新方法。     

A RHEOLOGICAL CONSOLIDATION MODEL FOR SOFT CLAY WITH CONSIDERATION OF TEMPERATURE AND THE ANALYTICAL

基于四元件流变模型,考虑温度影响,建立了宁波软黏土流变固结模型,并利用拉普拉斯变换得到瞬时加载条件下的考虑温度 影响的饱和软黏土流变固结解析解;利用不同温度下软黏土的流变固结试验结果,拟合获得了四元件流变模型参数,进行编程得到理论值, 并与试验值进行对比. 结果表明：该模型能较好反映宁波软黏土流变固结特性,计算结果与试验结果较为吻合.温度的升高导致渗透系数增大,在相同的时间内孔压消散 越快, 固结 越快.     

STUDY ON THE PARAMETER VALUE DOMAIN OF MACRO-MICRO CONSTITUTIVE MODEL

宏微观耦合本构模型的参数识别往往通过反分析方法进行,为了使参数识别结果具有高的置信度,需要确定合适的参数取值范围. 基于动态再结晶过程的微观机理以及相应本构方程的数学特征,提出一个确定参数取值范围的方法. 首先详细给出考虑动态再结晶的黏塑性本构模型,并根据模型构造物理机理,提出通过6 步确定该模型参数取值范围的方法;其次,对300M 低合金钢进行不同温度、应变速率下的热变形试验,测试宏观的流动应力-应变数据及微观的组织数据;然后应用提出的方法,依据试验数据,确定参数取值范围;最后,基于确定参数取值范围中获得的知识,对模型进行局部修改,使模型模拟结果更接近实验结果.     

七方程可压缩多相流模型的HLLC格式及应用

针对Saurel和Abgrall提出的两速度两压力的七方程可压缩多相流模型,改进了其数值解法并应用于模拟可压缩多介质流动问题.在Saurel等的算子分裂法基础上,根据Abgrall的多相流系统应满足速度和压力的均匀性不随时间改变的思想,推导了与HLLC格式一致的非守恒项离散格式以及体积分数发展方程的迎风格式.进一步,通过改变分裂步顺序,构造了稳健的结合算子分裂的三阶TVD龙格-库塔方法.最后通过几个一维和二维高密度比高压力比气液两相流算例,显示了该方法在计算精度和稳健性上的改进效果.      

A COMPLEX PARAMETER MODEL UPDATING METHOD FOR DAMPING CHARACTERISTICS OF FINITE ELEMENT MODELS

阻尼对于结构动力学响应具有重要的影响,但有限元模型一般很难对阻尼特性进行精确建模.基于实测频响函数,研究了一种有限元模型阻尼特性的复参数修正方法.以待修正区域各单元质量、刚度矩阵的比例修正系数为复修正参数,建立了单元矩阵比例修正的灵敏度方程直接算法,并对比分析了复修正参数与不同阻尼特性之间的数学关系.以六自由度集中参数模型和25杆平面桁架模型为例,验证了复参数修正方法在阻尼特性修正中的有效性.     

DYNAMIC ANALYSIS OF BOUNDED DOMAINS BY SBFE AND THE IMPROVED CONTINUED-FRACTION EXPANSION

基于比例边界有限元理论框架,通过采用连分式展开和引入辅助变量,将有限域的动力刚度矩阵和质量矩阵采用高阶的矩阵表示. 采用改进的连分式法求解比例边界有限元方程中的动力刚度矩阵. 通过增加连分式展开的阶数,该求解方法能包含动力分析的主要频率范围. 针对结构自由度较多的系统当连分式阶数逐渐增大时,原连分式算法可能会造成矩阵运算病态的问题,提出采用改进的连分式算法能有效地提高数值计算稳定性.通过对一正八边形的自由振动分析及矩形平面的时域分析,算例结果表明改进算法的鲁棒性更强,适合大规模系统的动力分析.     

深水套管切割力学模型

针对深水套管切割力学模型单一、切割泵冲选择困难等问题, 利用长杆铣削模型, 建立了深水套管切割力学模型, 计算出钻柱切割扭矩和切割允许的最大泵冲, 得出了切割扭矩占井口驱动扭矩的比例与作业水深之间的关系, 对深水弃井作业中套管割断判断指标的形成具有指导意义. 结合南海某两口井套管切割的现场数据, 验证了模型的正确性, 其计算结果可以为深水弃井作业中合理的泵冲参数选择提供理论参考, 为深水弃井作业提供技术保证, 进一步提高切割效率, 具有很好的应用价值.     

基于压缩蠕变试验的流变参数理论解析反演

大坝坝基的长期稳定性是大坝运营安全的重要保障, 而坝基岩体中的软弱夹层是影响其变形和稳定的重要因素. 为研究大岗山水电站坝基中辉绿岩脉软弱夹层在长期载荷作用下的变形机制, 在坝基边坡的试验平硐内垂直于软弱夹层进行了现场大型圆形刚性承压板压缩蠕变试验. 辨识得到可较准确表示其蠕变特性的5参量广义Kelvin模型, 克服了3 参量广义Kelvin模型收敛过快的缺陷. 基于弹性力学中的布辛涅斯克问题, 通过黏弹性理论中的拉普拉斯变换及逆变换, 推导出了刚性承压板下部岩体的5参量广义Kelvin模型表示的黏弹性变形公式, 并以此为基础反演得到流变参数.     

煤矿泥岩流变模型辨识及支护方法的确定

结合高家梁煤矿斜井泥岩的流变变形特征, 对泥岩进行分析建模, 设计有效的泥岩斜井支护措施. 通过对泥岩进行蠕变加卸载试验, 采用统一流变力学模型理论建立了能表现泥岩流变特性的弹性--黏性--黏弹性--黏弹塑性模型, 并确定了流变模型参数. 结合泥岩试验测得的泥岩力学特性和现场观测资料, 设计采用锚喷网联合支护方法, 并采用大型有限元软件ADINA进行数值模拟分析, 该支护方法能有效控制泥岩变形, 采用的锚喷网支护方法合理可靠.     

改进的非定常Bingham岩石蠕变模型及参数辨识

传统Bingham模型适用于描述蠕变的前两个阶段,为了更加准确地描述岩体蠕变的三个阶段,本文提出了一种改进的非定常Bingham蠕变模型,并推导出其蠕变本构方程．该模型在传统Bingham模型的基础上串联一个非定常黏壶元件,并考虑岩体弹性模量及原有黏壶元件蠕变参数的非定常化．基于该模型对泥岩的蠕变试验曲线进行拟合及参数辨识,结果表明该模型具有较高的拟合精度,可以准确地描述泥岩蠕变的三个阶段,并实现了在数值软件中的应用,为实际工程提供了一种简单有效的蠕变本构模型．     

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

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

INTERNAL SOLITARY WAVE LOADS EXPERIMENTS AND ITS THEORETICAL MODEL FOR A CYLINDRICAL STRUCTURE

在大型重力式密度分层水槽中, 对内孤立波与圆柱型结构的相互作用特性开展了系列实验. 基于两层流体中 内孤立波的KdV,eKdV和MCC理论, 建立了圆柱型结构内孤立波载荷的理论预报模型, 给出了该载荷理论预报模型中3类内孤立波理论的适用性条件.研究表明, 圆柱型结构内孤立波水平载荷包括水平Froude-Krylov力、附加质量力和拖曳力3个部分, 可以由Morison公式计算, 而内孤立波垂向载荷主要为垂向Froude-Krylov力, 可以由内孤立波诱导动压力计算.系列实验结果表明, 附加质量系数可以取为常数1.0, 拖曳力系数与内孤立波诱导速度场的雷诺数之间为指数函数关系, 而且基于理论预报模型的数值结果与系列实验结果吻合.