高温气体热化学反应的DSMC微观模型分析

热化学耦合的非平衡现象一直是高温气体热化学问题研究的难点,制约了诸如爆轰波胞格结构、低温点火速率等现象的分析.本文以高温氮气离解和氢氧燃烧中的链式置换反应为例,从微观反应概率、振动态指定的反应速率、热力学非平衡态的宏观反应速率、碰撞后的能量再分配等角度,分析了直接蒙特卡罗模拟中的典型化学反应模型(TCE,VFD,QK模型)的微观动力学性质.研究发现,无论是高活化能的高温离解反应还是低活化能的链式置换反应,实际参与反应的分子的振动能概率分布都偏离了平衡态的Boltzmann分布,包含较强振动能额外影响的VFD模型可以很好地模拟高温离解反应,而TCE(VFD的一个特例)和QK模型对活化能较低的链式置换反应的预测效果相对更好.此外,化学反应碰撞后的能量再分配应遵循微观细致平衡原理,细微的偏差都可能造成平动能和振动能难以达到最终的平衡状态.直接蒙特卡罗模拟的应用评估结果表明,化学反应的振动倾向对热化学耦合过程产生了明显的影响,特别是由于高振动能分子更多地参与了化学反应,气体平均振动能的下降将影响后续化学反应的进行.     

高温非平衡流动中的氧分子振动态精细分析

高超声速流动在头激波压缩后常处于高 温条件下的热化学非平衡状态. 本文采用态-态方法和双温度模型计算分析了一维正激波后和高超声速钝体绕流驻点线上的氧气热化学非平衡流动. 态-态方法将氧气的每个振动能级当成独立的组分,通过耦合 Euler 方程或驻点线上的降维 Navier-Stokes 方程,数值求解得 到了高温流动中的精细热化学非平衡状态. 而双温度模型假设氧气的振动能级服从 Boltzmann 分布,通过求解振动能方程得到振动温度. 一维正激波后热化学松弛过程的计算结果表明,态-态计算预测的温度分布和氧原子浓度分布较好地吻合了文 献中的实验结果,而经典的双温度模型的预测结果误差较大,且不同双温度模型的计算结果比较发散. 态-态方法详细地给出了所有振动能级的变化过程. 无论是正激波还是脱体激波后的流场,都是高振动能级首先得到激发;但是数密度大的低振动能级先达到热平衡,而高能级 分子要经过很长距离后才能达到热平衡. 在驻点附近,复合反应生成的氧气分子处于高振动能级,导致高振动能级分子数密度显著高于平衡分布. 计算还发现,经典双温度模型的离解反应速率明显偏离态-态计算结果,无法准确体现振动离解耦合效应对离解反应 速率的影响,但是 Park 双温度模型将离解失去的振动能取为 0.3$\sim 高超声速流动在头激波压缩后常处于高 温条件下的热化学非平衡状态. 本文采用态-态方法和双温度模型计算分析了一维正激波后和高超声速钝体绕流驻点线上的氧气热化学非平衡流动. 态-态方法将氧气的每个振动能级当成独立的组分,通过耦合 Euler 方程或驻点线上的降维 Navier-Stokes 方程,数值求解得 到了高温流动中的精细热化学非平衡状态. 而双温度模型假设氧气的振动能级服从 Boltzmann 分布,通过求解振动能方程得到振动温度. 一维正激波后热化学松弛过程的计算结果表明,态-态计算预测的温度分布和氧原子浓度分布较好地吻合了文 献中的实验结果,而经典的双温度模型的预测结果误差较大,且不同双温度模型的计算结果比较发散. 态-态方法详细地给出了所有振动能级的变化过程. 无论是正激波还是脱体激波后的流场,都是高振动能级首先得到激发;但是数密度大的低振动能级先达到热平衡,而高能级 分子要经过很长距离后才能达到热平衡. 在驻点附近,复合反应生成的氧气分子处于高振动能级,导致高振动能级分子数密度显著高于平衡分布. 计算还发现,经典双温度模型的离解反应速率明显偏离态-态计算结果,无法准确体现振动离解耦合效应对离解反应 速率的影响,但是 Park 双温度模型将离解失去的振动能取为 0.3$\sim $0.5 倍分子离解能是比较合理的.     

高温非平衡流动中的氧分子振动态精细分析

高超声速流动在头激波压缩后常处于高温条件下的热化学非平衡状态.本文采用态-态方法和双温度模型计算分析了一维正激波后和高超声速钝体绕流驻点线上的氧气热化学非平衡流动.态-态方法将氧气的每个振动能级当成独立的组分,通过耦合Euler方程或驻点线上的降维Navier-Stokes方程,数值求解得到了高温流动中的精细热化学非平衡状态.而双温度模型假设氧气的振动能级服从B oltzmann分布,通过求解振动能方程得到振动温度.一维正激波后热化学松弛过程的计算结果表明,态-态计算预测的温度分布和氧原子浓度分布较好地吻合了文献中的实验结果,而经典的双温度模型的预测结果误差较大,且不同双温度模型的计算结果比较发散.态-态方法详细地给出了所有振动能级的变化过程.无论是正激波还是脱体激波后的流场,都是高振动能级首先得到激发;但是数密度大的低振动能级先达到热平衡,而高能级分子要经过很长距离后才能达到热平衡.在驻点附近,复合反应生成的氧气分子处于高振动能级,导致高振动能级分子数密度显著高于平衡分布.计算还发现,经典双温度模型的离解反应速率明显偏离态-态计算结果,无法准确体现振动离解耦合效应对离解反应速率的影响,但是Park双温度模型将离解失去的振动能取为0.3～0.5倍分子离解能是比较合理的.     

高温热化学非平衡气动热试验与仿真技术研究进展

临近空间新型飞行器向全空域、更高马赫数发展,面临的气动热环境会越发恶劣,高温流场气动热预测技术是该类飞行器发展的关键技术之一.高超声速气流通过激波压缩或黏性阻滞减速,分子动能转化为内能,产生了高温.高温引起体分子振动、电子激发,伴随离解、电离反应等一系列复杂气动物理现象,其流场气动热预测面临诸多挑战.文章对高温热化学非平衡气动热预测技术的发展情况进行了分析探讨.首先,阐述了国内外高温气动热地面试验技术的发展历程,重点介绍分析了气动热风洞试验设备的模拟能力及目前试验测试技术的研究水平;然后,调研和讨论了高温气动热数值模拟研究现状,分别从热化学模型、辐射输运和壁面催化/烧蚀等多个角度探讨了热化学非平衡流场气动热数值模拟规律;最后,对气动热预测技术的发展趋势进行了讨论,提出了高温气动热试验与仿真技术后续应重点解决的问题.     

高超声速三维热化学非平衡流场的数值模拟

对三维高超声速热化学非平衡流场进行数值模拟,采用双温度热化学非平衡、11组元空气模型,考虑振动-离解耦合．差分格式采用沈清博士提出的“迎风型NND”格式,用熵修正方法消除了高超声速流数值模拟中的“carbuncle现象”．与LU-SGS方法结合,提高了单步计算效率和收敛性．数值模拟结果与文献结果进行了对比,并在弹道靶中进行了钢质圆球的弓形激波位置实验验证．计算结果与文献、实验的对比说明,三维热化学非平衡流计算程序可以精确地捕捉到强弓形激波,得到合理的空气动力系数．     

凝聚态物质振动-平动能量弛豫过程的分子动力学模拟

本文利用分子动力学计算机模拟方法,研究了稠密态双原子分子振动-平动弛豫速率与分子离解能、密度和温度的关系。发现振动弛豫速率随着分子离解能的增高而下降。这一现象与由光谱数据得到的结果是一致的。它可以用振动频率的下降来解释;分子振动弛豫速率随密度增大而加快,在我们所作的范围内,似乎看不到弛豫速率与温度有关。     

凝聚态物质振动-平动能量弛豫过程的分子动力学模拟

本文利用分子动力学计算机模拟方法,研究了稠密态双原子分子振动-平动弛豫速率与分子离解能、密度和温度的关系。发现振动弛豫速率随着分子离解能的增高而下降。这一现象与由光谱数据得到的结果是一致的。它可以用振动频率的下降来解释;分子振动弛豫速率随密度增大而加快,在我们所作的范围内,似乎看不到弛豫速率与温度有关。     

气固两相流中颗粒碰撞的Monte-Carlo数值模拟

利用颗粒碰撞动力学模型和颗粒几何碰撞率模型,采用Monte—Carlo算法来模拟颗粒之间碰撞,把该算法与求解雷诺应力-概率密度函数模型的有限差分-Monte Carlo算法耦合起来,对轴对称突扩通道内的两相旋流场进行了数值模拟,模拟结果表明,由于颗粒碰撞使颗粒的动能和湍动能在三个坐标方向上进行了再分配,从而导致颗粒的动能和湍动能在三个坐标方向上趋于各向同性;另外,由于颗粒碰撞破坏了颗粒-颗粒、颗粒-流体微团之间的速度关联,从而造成颗粒湍动能及两相速度脉动关联的降低。     

气体化学反应流动的DSMC／EPSM混合算法研究

发展了平衡粒子模拟方法(EPSM)，建立了与高温气体化学反应动力学理论相匹配的：EPSM耦合模型，并通过混合参数进行流区的自动识别，将：EPSM方法与蒙特卡罗直接模拟方法(OSMC)结合，构造了可模拟化学反应流动的DSMC／EPSM混合算法。应用该算法对汲及化学反应的二维高超音速竖板绕流流场进行模拟，将结果与DSMC方法的结果进行比较，验证了新算法对求解化学反应流动的可行性。将混合算法的计算效率与DSMC方法的计算效率进行比较，发现混合算法能够大大提高计算效率。     

水动力作用下流冰撞击闸墩的动力响应研究

高寒地区河冰撞击河道的闸墩结构会产生极端冰载荷和冰激振动,水的动力效应使得碰撞过程更加复杂。采用任意拉格朗日-欧拉流固耦合方法,考虑作用在流冰和闸墩表面的流体力,建立了水-冰-闸墩耦合模型,探究了偶然极端条件下冰-闸墩碰撞的力学特性,设计了冰-砼碰撞实验。结果表明：冰-砼碰撞实验中,撞击力的模拟结果与实验结果吻合良好;对流固耦合的水动力效应分析发现,水-冰-闸墩耦合模型能够体现水的流体特性,在流冰撞击闸墩近场逼近过程中,初始时刻水的动力效应能够增加流冰的动能,撞击楔入闸墩过程中,水介质形成一个瞬态高压力场,产生水垫效应吸收冰体部分动能,从而抑制流冰运动;在不同流冰体积和压缩强度工况下,闸墩结构所承受的冰力随着流冰体积的增大而增大,流冰压缩强度对冰力的影响较小,流冰损伤与闸墩结构响应主要集中在碰撞接触区,流冰撞击闸墩结构引起冰激振动,流冰体积对闸墩振动加速度的影响较大,相同体积的流冰随着压缩强度的增大,振动幅值差异不明显,表明流冰体积是影响冰-闸墩碰撞的关键参数。     

高雷诺数下多柱绕流特性研究

采用改进的延迟分离涡方法数值模拟了高雷诺数下的柱体绕流,包括单圆柱绕流、单方柱绕流、串列双圆柱绕流和串列双方柱绕流,研究了不同雷诺数下圆柱绕流与方柱绕流的水动力特性. 计算结果与实验数据及其他文献的数值计算结果吻合良好,研究表明,单方柱绕流在范围内未出现类似于单圆柱绕流的阻力危机现象,其平均阻力系数\overline{C_d}、升力系数均方根{C\mathrm{\text{'}}}_l及斯特劳哈尔数\mathit{St}维持在一定范围内波动. 串列双圆柱绕流与串列双方柱绕流中,均选取L/D=2.0, 2.5, 3.0, 3.5和4.0这五中间距比进行计算. 串列双圆柱绕流中,当\mathit{Re}=2.2\times 10^4时,在内存在一临界间距比(L_c/D)使得L_c/D前后上下游圆柱的升阻力系数发生跳跃性变化,且当时,下游圆柱的阻力系数为负数. 而当\mathit{Re}=3.0\times 10^6时,则不存在临界间距比,且下游圆柱的阻力系数始终为正数. 串列双方柱绕流在\mathit{Re}=1.6\times 10^4和\mathit{Re}=1.0\times 10^6两种工况下的临界间距比分别处于和区间内,且当时,两个雷诺数下的下游方柱阻力系数均为负数.      

应力诱发界面迁移下晶内孔洞的演化

随着微电子技术的迅猛发展, 集成电路中内连导线的失效问题引起广泛关注. 内连导线内部孔洞萌生、长大、漂移和失稳变形成狭长裂纹, 从而导致电路的开路失效. 这是内连导线失效的常见形式. 而界面迁移是导致微结构形态演化的主要机制之一. 本文基于界面迁移下微结构演化的经典理论和弱解描述, 建立了应力诱发界面迁移下微结构演化的有限单元法, 并验证了算法的可靠性. 对铜内连导线中晶内孔洞的演化进行了数值模拟, 详细分析了应力、线宽及形态比对晶内孔洞演化的影响. 研究结果表明, 椭圆形晶内孔洞存在生长和收缩两种演化分叉趋势. 通过大量数值分析得到了晶内孔洞演化的临界应力{\stackrel{\text{?}}{\sigma }}_c 、临界线宽{\stackrel{\text{?}}{h}}_c 和临界形态比{\beta }_c . 当\stackrel{\text{?}}{h}\mathrm{?}{\stackrel{\text{?}}{h}}_c, \stackrel{\text{?}}{\sigma }或\stackrel{\text{?}}{h} 时, 晶内孔洞会沿长轴长大; 反之, 晶内孔洞会收缩甚至愈合. 此外, 应力\beta 越大、线宽\stackrel{\text{?}}{\sigma } 越小或形态比\stackrel{\text{?}}{h} 越大, 晶内孔洞越易发生长大, 且孔洞面积增大速度越快; \beta 越小、\stackrel{\text{?}}{\sigma } 越大或\stackrel{\text{?}}{h} 越小, 晶内孔洞越易发生收缩, 且孔洞面积减小速度越快.      

一种考虑液相惯性的一维准静力固结模型

经典Terzaghi一维固结理论不考虑孔隙流体惯性影响,且该理论在不同时期模型推导和表述结果差别较大,导致当前仍存在诸多困惑甚至认识混乱的现象. 在笔者前期研究大变形动力固结理论框架内,忽略固相惯性而重点考虑液相惯性影响,经过合理简化建立反映孔隙流体惯性的一维小变形固结波动模型. 该固结波模型具有频散和耗散特性. 采用分离变量法,可得到单面排水和瞬时加载条件下无量纲形式固结波解析解答. 算例分析结果表明:固结波发展规律受无量纲数D_{\mathrm{c}}变化影响而呈现不同性态;D_{\mathrm{c}}数值较大时固结波响应会出现阶跃和正负波动现象;当D_{\mathrm{c}}值较小时,可能出现Mandel-Cryer效应等特殊现象. 通过对早期和后期Terzaghi固结模型的分析和对比,初步探明Terzaghi固结理论模型内部的矛盾性,在普通土体坐标和固相体积坐标两种不同解读条件下,早期Terzaghi (1923,1925)模型可以分别诠释为具有小变形和大变形属性的不同固结模型. 在经典一维固结理论模型的不同诠释背景下,固结波模型也可以据此作出相应拓展和表述. 固结波理论揭示缩尺固结试验中土体物理力学参数与固结波响应两种因素之间存在一种不确定性矛盾,据此建议微观土力学研究重视尺度效应. 固结波模型的意义还在于,可为Terzaghi经典固结模型理论精度分析提供新的依据.      

扩展腔对方波型微混合器混合性能的影响研究

微混合器凭借节约试剂、混合强度高和易于集成等优点,在材料合成、医药制备和生化检测等领域中具有广泛的应用. 为了进一步提高混合性能,保证混合过程的安全性及生化反应结果的准确性,设计了一种带扩展腔的新型方波型微混合器. 在综合考虑混合强度和压降的前提下,通过实验研究和数值模拟分析了窄缝宽度、窄缝长度和扩展腔高度对微混合器混合性能的影响并得到了不同雷诺数\mathit{Re}条件下的最优结构参数. 与方波型微混合器的混合性能进行比较,发现\mathit{Re}=20时,带扩展腔的方波型微混合器的混合强度更高,其中\mathit{Re}时两者混合强度相差最多,可达12%. 在相同\mathit{Re}下,带扩展腔的方波型微混合器的压降要低于方波型微混合器. 对带扩展腔的方波型微混合器进行内部流场分析,发现扩展腔结构能在流体层流状态的基础上引入涡流,使通道中流体的流动状态发生改变、对流增强,进而混合性能提高.      

高超声速激波湍流边界层干扰直接数值模拟研究

高超声速激波与湍流边界层干扰会导致飞行器表面出现局部热流峰值,严重影响飞行器气动性能和飞行安全. 针对高马赫数激波干扰问题,以往数值研究多采用雷诺平均方法,而在直接数值模拟方面的相关工作较为少见. 开展高超声速激波与湍流边界层干扰的直接数值模拟研究,有助于进一步提升对其复杂流动机理认识和理解,同时也将为现有湍流模型和亚格子应力模型的改进提供理论依据. 采用直接数值模拟方法对来流马赫数6.0,34°压缩拐角内激波与湍流边界层的干扰问题进行了研究. 基于雷诺应力各向异性张量,分析了高超声速湍流边界层在压缩拐角内的演化特性. 通过对湍动能输运方程的逐项分析,系统地研究了可压缩效应对湍动能及其输运的影响机制. 采用动态模态分解方法,探讨了干扰流场的非定常运动历程. 研究结果表明,随着湍流边界层往下游发展,近壁湍流的雷诺应力状态由两组元轴对称状态逐渐演化为两组元状态,外层区域则由轴对称膨胀趋近于各向同性. 干扰流场内存在强内在压缩性效应(声效应),其对湍动能输运的影响主要体现在压力--膨胀项,而对膨胀--耗散项影响较小. 高超声速下压缩拐角内的非定常运动仍存在以分离泡膨胀/收缩为特征的低频振荡特性,其物理机制与分离泡剪切层密切相关.      

等离子体激励器诱导射流的湍流特性研究

为了进一步掌握等离子体流动控制机理, 完善等离子体激励器数学模型, 提升等离子体激励器扰动能力, 采用粒子图像测速技术, 在静止空气下开展了介质阻挡放电等离子体激励器诱导射流特性研究. 实验时, 将非对称布局激励器布置在平板模型上, 随后将带有激励器的模型放置在有机玻璃箱内, 从而避免环境气流对测试结果的影响. 基于激励器诱导流场, 分析了激励电压对诱导射流特性的影响, 揭示了较高电压下诱导射流近壁区的拟序结构, 获得了卷起涡、二次涡等拟序结构的演化发展过程, 计算了卷起涡脱落频率, 阐述了卷起涡与启动涡的区别, 初步探索了卷起涡的耗散机制. 结果表明: (1)层流射流不能完全概括等离子体诱导射流特性, 激励电压是影响射流特性的重要参数. 当电压较低时, 诱导射流为层流射流; 当电压较高时, 诱导射流的雷诺数提高, 射流剪切层不稳定, 层流射流逐渐发展为湍流射流. (2)等离子体诱导湍流射流包含着卷起涡、二次涡等拟序结构; 在固定电压下, 这些涡结构存在恒定的卷起频率. (3)当激励电压较高时, 流动不稳定使得卷起涡发生了拉伸、变形, 引起了流场湍动能增大, 从而加速了卷起涡的耗散. 研究结果为全面认识激励器射流特性, 进一步挖掘激励器卷吸掺混能力, 提升激励器控制能力积累基础.      

PMMA膨胀环动态拉伸碎裂实验研究

在强动载作用下, 脆性材料的碎裂问题是一个重要的研究课题, 而脆性材料在冲击拉伸载荷下的力学行为的实验研究相对较匮乏. 提出了一种动态拉伸断(碎)裂的液压膨胀环实验技术, 可用于准脆性/脆性材料的动态拉伸. 利用该技术对有机玻璃(PMMA)圆环试件进行了不同膨胀速度下的动态碎裂实验研究. 从回收碎片的断口形貌和碎片内部残余裂纹观察可知试件的破碎由环向拉伸应力造成, 碎片断口处发出的稀疏波会将周围的拉伸应力卸载, 从而抑制其他裂纹的进一步发展. 利用超高速相机记录了试件的膨胀碎裂过程, 利用DISAR激光速度干涉仪获得了试件外表面粒子的径向膨胀速度历史, 通过试件上的应变片获得了试件的应变历史和断裂应变. 实验结果表明: 在拉伸应变率150\mathrm{~}500s^{-1}范围, 材料的动态断裂应变低于准静态加载下的断裂应变, 体现出“动脆”现象; 随着加载应变率的提高, PMMA 材料的碎片尺寸减小; 无量纲化的PMMA圆环的平均碎片尺寸介于韧性碎裂模型和脆性碎裂模型的预测数值之间, 反映出材料的准脆性特性.      

功能梯度材料动力学问题的POD模型降阶分析

为了快速分析非均质材料结构在复杂载荷作用下的动态响应, 提出一种模型降阶方法, 只需计算结构在简单均质材料情况下的动力学问题, 进而用其计算结果对非均质材料结构进行分析. 首先, 采用结构内部任意一点处的材料参数值作为整个结构的材料参数, 利用有限元分析软件计算该均质材料结构在动态载荷作用下的位移场建立数据库, 该数据库包含计算模型各个节点(自由度为L)在某时间段内L个时刻的位移; 其次, 对数据库中的信息按照时间离散的特定方式组集成瞬像矩阵, 并利用特征正交分解方法对其进行分解, 得到该模型的H个特征正交基底, 选取其中能反应模型主要特征的个(其中\mathrm{~})作为一组最优基底, 通过这组基底建立模型的低阶离散控制方程; 最后, 求解低阶离散微分方程组, 得到功能梯度材料结构在复杂载荷作用下的位移场. 文中分别给出二维和三维算例, 比较了降阶模型和全阶模型计算结果, 验证了该方法的有效性, 并且计算效率能提高1{\mathrm{}}^{1)}2个数量级.      

错列角度对双圆柱涡激振动影响的数值模拟研究

为研究错列角度\alpha对双圆柱涡激振动问题的影响,采用自主研发的基于CIP (constrained interpolation profile)方法的数值模型,对雷诺数\mathit{Re}=100、错列角度\alpha =0^{°}\mathrm{~}90^{°} (间隔15 #x00B0;)的等直径双圆柱涡激振动问题进行数值模拟. 模型在笛卡尔网格系统下建立,采用具有三阶精度的 CIP 方法求解 N-S (Navier--Stokes)方程,采用浸入边界法处理流--固耦合问题,避免了任意拉欧方法下的网格畸变和重叠动网格技术中的大量信息交换问题,保证了模型的计算效率. 重点分析不同错列角度\alpha上下游圆柱的升阻力系数、位移响应、涡脱频率和尾涡模态等. 结果表明:折合速度U_r=2.0\mathrm{~}3.0时,上下游圆柱升阻力随错列角度的增大基本呈单调增大的趋势;U_{\mathrm{r}}=5.0\mathrm{~}8.0时,随错列角度的增大,上下游圆柱阻力变化较小,升力呈“上凸”趋势,在\alpha =15^{°}\mathrm{~}30^{°}取得最大值;U_r=10.0\mathrm{~}13.0时,随错列角度的增大,上下游圆柱阻力变化较小,升力呈“下凹”趋势,在\alpha =30^{°}\mathrm{~}45^{°}取得最小值,且柱体横流向振幅和升力没有明显的对应关系. 最后,结合尾涡模态对以上规律的成因进行分析. 研究结果可为相关海洋工程设计提供参考.      

复杂加载条件下的砂土本构模型

试验表明,饱和砂土的应力应变关系具有显著的密度以及压力依存性,上述两点构成了描述砂土静力加载下变形特性无法忽视的因素. 此外,在循环加载等复杂加载作用下,砂土还会表现出明显的应力诱导各向异性以及相变转换特性. 基于在e--p空间中存在唯一的临界状态线这一基本假定,通过在e--p空间中引入当前状态点与临界状态线的距离R来作为反映密度与压力依存特性的状态参量, 将变相应力比以及峰值应力比表达为状态参量的指数函数,将上述应力比参量引入到统一硬化参量中可准确地反映初始状态下围压、密度 对于单调加载下应力应变关系的影响规律,能描述砂土剪缩、剪胀,应变软化、硬化等特性. 采用非相关联流动法则,p--q空间中采用水滴型屈服面,塑性势面为椭圆面,松砂在单调加载下的静态液化现象也可描述. 为反映循环加载下塑性体积应变的累积特性以及塑形偏应变的滞回特性,在循环加载下将状态参量R表达为应力比参量,并在硬化参数中引入描述应力诱导各向异性特性的旋转硬化部分,所提模型可有效地描述循环加载下剪切模量的衰减特性、刚度衰化性质、强度减小特性,在不排水约束作用下,则会产生往返活动性现象. 通过一系列的模型模拟与试验结果对比,验证了本构模型的有效性及适用性.