基于改进Bingham模型的软岩参数非定常三维非线性黏弹塑性蠕变本构研究

传统线性元件组合模型难以描述岩石非线性加速蠕变阶段特性,给工程应用带来了巨大困难。提出一种新的应变触发非线性黏壶元件,并将其与Bingham模型串联,建立了可以描述岩石等速和加速蠕变特性的一维本构模型,进而推导出相应的三维蠕变本构模型,再将三轴压缩过程中岩石弹性模量的衰减方程引入该三维蠕变本构方程,得到一个能反映岩石蠕变全过程特性的三维非线性黏弹塑蠕变本构模型,根据泥岩蠕变试验数据对该模型参数进行了辨识。模型理论计算值与泥岩蠕变试验数据的对比表明改进的Bingham模型不仅可以充分反映软岩的等速蠕变特性,还可以很好地描述软岩的衰减和加速蠕变规律,且模型元件少,组合形式简单,为软岩非线性蠕变本构模型研究提供了新的借鉴。     

岩石空心圆柱扭剪仪研制的重点问题及研究进展

随着深部岩体工程的发展,开挖卸荷导致的应力主轴旋转问题引起了广泛关注。首先,分析总结了岩体工程中存在的应力主轴旋转问题及其对岩体工程稳定性的影响;其次,介绍了土体空心圆柱扭剪仪的发展现状及其特点,从荷载的施加方式、试样尺寸两个方面指出了研制岩石空心圆柱扭剪仪的关键问题,并与土体空心圆柱扭剪仪进行了对比,提出了新的扭矩施加技术和合理的岩石空心圆柱试样尺寸;最后,总结概括了考虑应力主轴旋转的土体本构模型的研究成果,并对考虑应力主轴旋转的岩石本构关系的研究进行了展望。该工作将对岩石空心圆柱扭剪仪的研制和相关理论的研究提供一定的参考和指导。     

基于旁压试验土体弹塑性本构模型初探

在大量旁压试验数据分析的基础上,通过曲线拟合,建立了旁压试验弹塑性阶段曲线的椭圆方程,利用土体SMP屈服准则和Rowe流动法则,推导出土体塑性阶段应力增量与应变增量间关系矩阵。在弹性阶段,假设土体应力-应变服从广义虎克定律,建立了基于旁压试验的土体弹塑性本构模型。编制了相应的计算程序,将文中模型计算结果与实际旁压试验曲线进行对比,初步验证了模型的正确性。将本构模型编译为ABAQUS自定义材料子程序UMAT,通过有限元对比分析,文中模型计算变形较弹性模型小,较摩尔-库仑模型大,与模型建立的假设一致。基于旁压试验的土体弹塑性本构模型参数少,且易于获取,便于在实际工程中应用。     

西安地区地裂缝带Q3原状黄土流变特性试验研究

汾渭盆地是我国乃至世界上地裂缝最为发育、灾害最为严重的地区。考虑在漫长的地质历史过程中黄土的流变特性对地裂缝破裂扩展模式的影响,采用分级循环加卸载方式,对西安地区长安地裂缝（fc1）带Q3原状黄土进行了不同围压下的固结不排水三轴流变特性试验,得到了地裂缝带Q3原状黄土的蠕变加、卸载曲线和应变速率与时间关系曲线。试验结果表明,不同围压下地裂缝带Q3原状黄土表现为低应力下为减速蠕变,高应力下为减速蠕变和等速蠕变的特征。并基于统一流变力学模型理论,建立了适用于西安地区地裂缝带Q3原状黄土的三维蠕变本构模型,获得了相应的本构模型参数。该成果为研究隐伏地裂缝破裂扩展模式及其力学机制奠定了基础。     

断层地震孕育和发生的不稳定性模型

在以前的工作中,考虑直立走滑型断层地震,假设断层面微元破裂强度遵循Weibull概率分布,由细观力学方法推导出断层面的宏观本构关系是一个非线性函数,表现为弹性-软化塑性特征,在此基础上用稳定性理论研究了地震稳定性问题.而实际断层大多是倾斜的,为此,本文首先建立了由围岩和倾斜断层构成的平面地震力学模型,采用宏观的断层载荷-变形的全过程曲线,详细讨论了倾斜断层地震的不稳定性问题.结果表明,远场一旦施加位移,断层也同时错动,这可能与实际情况不符合.为了更好的模拟断层的初始能量累计过程,进一步对断层本构模型进行改进.考虑断层面破裂强度,采用Coulomb破裂准则,则断层表现为刚塑性本构关系,只有当断层面剪应力达到一个临界值时,断层才开始错动.研究表明,对于倾斜断层地震,与直立走滑型断层地震一样,系统刚度比β（围岩切线刚度与断层刚度最大值之比）是决定地震失稳的重要参数,只有当β<1时才会出现地震失稳,且伴随应力突跳和围岩应变能释放.当β≥1时,仅仅是断层无震滑动,不会发生地震.在远场应施以位移形式边界条件,以致地震失稳发生在平衡路径的位移转向点并伴有应力突跳.     

Bouc-Wen土体动力模型的阈值应变研究

探讨Bouc-Wen土体动力模型的适用条件与参数对滞回圈的影响等。基于热力学基本定律,研究了Bouc-Wen土体动力模型的增量耗散函数表达式,讨论了不同动应变水平下某筑坝土石料的能量耗散机制。发现筑坝土石料的动力特性存在2个阈值应变,即第l和第2阈值应变。其对应的模量比分别为0.94左右及0.50～0.80之间。第2阈值应变与传统意义上以孔压升高或体积变化为标准的门槛应变相当。对进一步认识土体动应力-应变特性有一定意义。     

考虑残余含气量的非饱和土的水力耦合本构模型

大量的非饱和土干湿循环试验表明,当土体处于吸湿过程直至吸力降低为0 kPa时,土体并不能达到完全饱和状态,还存在一定的残余气体。在高饱和度时,由于残余气体以封闭气泡的形式分布在土体中,土体呈现较大的压缩性,使其与饱和土的性质不同。在这种状态下,现有的非饱和土本构模型预测到的土饱和度为100%,与试验结果存在一定的偏差。为了使本构模型在高饱和度状态时具有较高的精度,对非饱和土的毛细滞回和塑性变形耦合本构模型进行了修正,使其能够考虑残余含气量的影响。通过预测与实测结果比较,证明了新模型能够有效地模拟残余含气量对非饱和土力学特性的影响。     

多块体形状堆石体碾压颗粒破碎数值模拟

通过6种典型堆石块体的形状近似,分别采用两种接触本构模型建立了多块体形状堆石体离散元数值模型,研究堆石体在碾压荷载作用下的颗粒破碎过程,建立颗粒破碎的量化计算方法,分析碾压前、后堆石级配曲线的变化,讨论接触本构模型和颗粒形状对块体破碎的影响。模拟结果显示,碾压荷载下堆石颗粒以张拉破碎为主,随着碾压遍数的增加,局部开始出现剪切破碎;提出的颗粒破碎量化计算方法,在大粒径范围对粒径变化幅度预测偏大,但级配曲线整体趋势与实测结果比较吻合;相比较于接触连接模型,平行连接模型与现场碾压试验结果更接近。6种块体形状的数值结果显示,随着形状系数的增加,在其他条件不变的情况下颗粒破碎率逐渐降低,其中类长方形块体在碾压荷载作用下颗粒破碎最明显。相比较于纯圆颗粒或者单一非圆颗粒,采用的6种颗粒形状建立的堆石体振动碾压离散元模型,更接近现场实际情况。     

地铁隧道地震反应非线性分析

本文对动力有限元法中的一系列问题进行了讨论,对土体动力本构模型和土-结构接触面模型的设定提出了建议．在此基础上以天津地铁1、2号线工程为实例,分别对水平地震作用下圆形盾构隧道的动位移、动应力与动应变以及地表位移等动力反应进行了非线性分析,得出了一些有益的规律和结论,从而为地铁隧道的抗震设计与研究提供参考和依据。     

深部人工冻土黏弹塑损伤耦合本构模型研究及其验证

采用改进的Mohr-Coulomb屈服准则(Zienkiewicz-Pande准则), 建立了能描述人工冻土3个阶段蠕变特性的黏弹塑损伤耦合本构模型. 基于该本构模型建立了模型的柔度矩阵, 将本构模型嵌入到ADINA有限元程序中. 计算表明: 数值模拟得到的人工冻土蠕变特性曲线与实测值吻合良好, 最大误差为4.8%, 说明基于改进的Mohr-Coulomb准则建立的黏弹塑损伤耦合本构模型能够准确的描述人工冻土剪切特性和3个阶段的蠕变特性.获得的本构模型具有参数少、容易确定等优点. 深部人工冻土黏弹塑损伤耦合本构模型的建立, 对于人工冻土结构物长期稳定分析和预测具有重大的理论指导意义.