非等轴颗粒付林图解

非等轴颗粒付林图解由直线a=1、b=1、a=b和曲线a2b=1、ab2=1分割成10个微区,并对应于不同的e_x、e_y、e_z组合和应变类型。这些直线和曲线分别对应于特殊的e_x、e_y、e_z组合和岩石应变类型。岩石的应变类型有18种。函数a.b的表征方程有4种类型,并分别表征直线a=1/k₁和b=1/k₂所分割的4个区间。a,b表征方程的确定取决于岩石总的变形行为。此图解可用于非等轴颗粒标志体应变测量和应变类型的确定和图示,文章给出了应用实例.      

煤样全应力应变过程中的渗透系数—应变方程

借助现代化的电液伺服岩石力学试验系统,以数控瞬态渗透法进行了全应力应变过程的软煤样渗透特性试验。实验中首次设置了环向应变传感器,得出煤样渗透性与主应力差、轴应变、体积应变关系曲线,并拟合出相应方程。从中可知,煤样全应力应变过程对应的渗透系数是体积应变的双值函数,体积缩小时为2次多项式,体积膨胀时为5次多项式。该方程用在应力场—渗流场耦合的数值分析中,可使计算结果更符合工程实际。      

真空-堆载预压作用下软土蠕变特性试验研究

采用改进的、可施加负压的三轴仪开展了真空预压、堆载预压以及真空-堆载联合预压作用下软土的固结蠕变试验,描述了加载率、应力比和时间等对软土蠕变特性的影响,分析了轴向应变（率）、体积应变（率）及偏应变（率）与应力比和时间之间的关系。结果表明：轴向应变率和体积应变率与时间的对数关系并非线性,但经过若干天以后可近似认为是线性关系;在不同应力比n下,体应变与时间的关系可用双曲线方程来表示;偏应变、偏应变率与时间的关系符合双曲线方程,通过温州软黏土样的蠕变试验结果验证了其有效性。     

循环荷载下饱和软黏土的累积塑性应变试验研究

通过以交通荷载为背景的饱和重塑软黏土室内不排水动三轴试验,研究了循环荷载作用下饱和重塑软黏土的累积塑性应变发展形态,可分为3种类型：稳定型、破坏型和临界型。根据稳定型累积塑性应变发展曲线特点,提出了饱和软黏土的稳定累积塑性应变方程,并通过试验结果分析了该方程中各拟合参数与动应力幅值、固结围压和静偏应力的变化规律,同时提出了求解无静偏应力条件下软黏土临界动应力的解析方法。通过动三轴试验结果,提出了含动应力幅值、固结围压、静偏应力和循环周次等影响因素的累积塑性应变拟合模型。     

隔离非线性平面应变单元模型及其在Drucker-Prager模型中的应用

提出了一种基于隔离非线性理论的平面应变单元分析模型。对于线性的四边形等参单元,采用2×2个高斯积分点作为单元的非线性应变插值点。考虑到平面应变单元面外无应变分量存在,故仅在面内的3个方向建立非线性应变场以及相应的控制方程。采用Woodbury公式与组合近似法联合求解控制方程,使得整个非线性分析过程仅为多次的初始弹性刚度矩阵回代计算以及稀疏矩阵与向量的乘积,极大地提高了控制方程的求解效率。基于时间复杂度的计算效率分析表明：本文算法的计算效率相对传统变刚度法显著提高,且相对于精确的Woodbury公式则大大提高了非线性自由度的临界比例。将Drucker- Prager准则的转移应力解析解法用于典型平面应变模型的非线性分析,数值算例验证了平面应变单元模型的正确性以及算法的高效性,为平面应变类型的非线性分析提供了一种新思路。     

隔离非线性平面应变单元模型及其在Drucker-Prager模型中的应用

提出了一种基于隔离非线性理论的平面应变单元分析模型。对于线性的四边形等参单元,采用2×2个高斯积分点作为单元的非线性应变插值点。考虑到平面应变单元面外无应变分量存在,故仅在面内的3个方向建立非线性应变场以及相应的控制方程。采用Woodbury公式与组合近似法联合求解控制方程,使得整个非线性分析过程仅为多次的初始弹性刚度矩阵回代计算以及稀疏矩阵与向量的乘积,极大地提高了控制方程的求解效率。基于时间复杂度的计算效率分析表明:所提算法的计算效率相对传统变刚度法显著提高,且相对于精确的Woodbury公式,则大大提高了非线性自由度的临界比例。将Drucker-Prager准则的转移应力解析解法用于典型平面应变模型的非线性分析,数值算例验证了平面应变单元模型的正确性以及算法的高效性,为平面应变类型的非线性分析提供了一种新思路。     

非饱和土渗透函数方程的间接确定

利用应力应变控制式非饱和土三轴仪,测定了非饱和土的基质吸力和体积含水量,由此得出粉质粘土一系列的土-水特征曲线,并由非线性曲线拟合得出土-水特征曲线方程,根据其与渗透函数的关系,推导出Mualem渗透函数方程的具体形式,由其计算得出的非饱和土的渗透系数和直接测量法得到的渗透系数较为一致,从而进一步证明了Mualem方程的具体形式的合理性。     

原状Q2黄土三轴剪切细观结构演化定量研究

基于非饱和原状Q2黄土CT-三轴剪切试验得到的结构演化清晰的CT图像及CT数据,提出利用CT数随偏应力的关系确定屈服应力的新方法。基于CT数定义了结构性参数和结构演化变量;分析了结构演化变量与净围压和吸力的关系,相对于应变（包括偏应变和体应变）,当吸力相同时,净围压越大,结构演化变量值越大;当净围压相同时,吸力越大,结构演化变量值越小。建立了非饱和原状Q2黄土三轴剪切过程的结构演化方程,该方程为结构演化变量与偏应变和体应变的关系,能同时反映净围压和吸力等对结构演化的影响。利用提出的结构演化方程对试验数据进行拟合,拟合结果较理想。     

基于动态CT识别的冻土单轴压缩损伤本构模型

采用连续介质力学与热力学方法,建立了冻土单轴压缩损伤本构模型。以兰州冻结黄土为例,基于冻土单轴压缩动态CT试验,对冻土附加损伤的2个阶段,即塑性损伤和微裂纹扩展损伤分别采用硬化曲线法和动态CT识别方法进行了损伤计算,在此基础上给出了冻土附加损伤与弹性应变的对应关系,建立了基于试验的冻土损伤演化方程,并进行了有效应力的计算。     

极射赤平投影中的数学方法

本文提供了严格依据赤平投影原理并利用空间解析几何知识推导出的公式计算方法，可以写出已知产状的线或面的赤平投影方程，也能计算两线所共的面，两面的交线，线或面之间的夹角以及它们的平分线或平分面，可以进行真倾斜和视倾斜的换算，倾斜和侧伏的换算；模拟线面绕轴的旋转以及恢复节理的原始产状；最后阐述用计算机进行节理统计和绘制极点等密图的原理。所用的向量表示方法直观表达了相关构造要素之间的空间关系，可用于计算机辅助教学，也能用于生产实践，精度可达0.01°。     

层状地基群桩沉降计算的剪切位移解析算法

考虑群桩的“束缚作用”,基于剪切位移法的理论,提出了竖向荷载作用下用于层状地基大规模群桩沉降分析的简捷实用的解析算法。以单桩位移积分方程为基础,导出了桩顶位移与轴力和桩底位移与轴力之间的关系,考虑桩-桩相互作用,得出了计算群桩沉降的柔度矩阵方程。推导过程中,桩被分成任意n段,因此该方法可以用于地基土任意分层的群桩沉降计算。算例分析表明,该方法与边混合法和界元法有较好的一致性。     

脉动孔隙水压下三轴加压对煤岩损伤变形的影响

为研究脉动孔隙水压作用下低透性松软煤岩在不同三轴加压条件下的损伤变形规律,采用RLW-2000M煤岩流变仪进行了不同轴压和围压系列下的试验研究。结果表明：煤岩破坏的应变曲线具有先疏后密再疏的变化规律,且累积应变和残余应变随着围压的增加而减小,而随着轴压的增加而增大;平均单次耗散能随着围压的增加而减小,随着轴压的增加而增加,其关系符合指数函数变化特征;采用累积耗散能计算了煤岩破坏过程中的损伤变化特征,并采用Logistic方程逆函数对其进行了拟合分析,得到了煤岩损伤演化模型,为预测煤岩破坏周期提供了依据。研究结果为具有不同地应力条件下的低透气性松软煤层施工脉动水力压裂提供了一定的理论依据。     

循环荷载下考虑主应力旋转的软土力学响应研究

目前用于轨道路堤分析计算的参数多是通过三轴试验获得,而列车荷载经过时土单元应力路径与循环三轴试验加载路径的显著差异可能导致预测失真而引发工程问题。针对某城市轨道进行了三维计算,并分析了移动列车荷载下土单元应力幅值、循环周数、主应力轴旋转的分布规律。在该基础上,采用空心圆柱循环扭剪试验对该复杂应力路径进行了模拟,研究了饱和软黏土的孔压及变形累积特性。结果表明,列车荷载经过时地基土单元大主应力将在(－90°, 90°)内发生旋转;该主应力轴旋转将显著促进软黏土孔压和应变的累积,竖向应力幅值为15 kPa时,循环扭剪试验产生的孔压值比循环三轴试验高77%,累积应变增大了近50%;随着循环应力水平提高,二者累积孔压及累积应变的差值还会进一步增大,甚至出现循环三轴下仅产生较小应变而循环扭剪下已破坏的本质性差异。     

考虑冻融循环的季冻区高等级公路路基填土本构模型研究

在季冻区高等级公路路基场地取土,制成重塑土试件进行冻融循环试验。而后对经历不同冻融循环次数的土体进行三轴固结排水（CD）剪切试验,测得应力、应变数据。在大量试验数据的基础上,遵循岩土本构变化规律,以岩土弹塑性力学为理论基础,在p-q空间,以椭圆方程拟合体变屈服面,以抛物线方程拟合剪切屈服面,采用塑性关联流动,得到考虑冻融循环条件下融土的弹塑性本构方程。通过编程计算后,对比计算与试验结果,证明该双屈服面本构模型能够较为准确地预估冻融循环土体的应力-应变关系,该本构模型对季冻区路基土体长期稳定性分析和工程预测具有重要参考意义。     

循环荷载下饱和重塑黏质粉土的动力特性研究

选取宁波饱和重塑黏质粉土,开展了黏质粉土的动三轴试验,研究了围压、动应力、排水条件、温度及频率等因素对土动力特性的影响,提出了采用相对动应力来分析不同围压、不同动应力下土的累积塑性应变变化规律,并建立了孔压—累积塑性应变关系的经验公式,结果表明:累积塑性应变随着相对动应力的增大而增大;温度越高,累积塑性应变值越小,且相同温度增量（△T=15℃）条件下的累积塑性应变增量也随之减小;建立的孔压-累积塑性应变经验公式可用于长期振动荷载作用后黏质粉土的孔压计算。     

基于Lade-Duncan准则的平面应变强度公式

在弹性-理想塑性和相关联流动法则的假设下,基于Lade-Duncan准则,得到了平面应变破坏的中主应力条件,提出一种新的平面应变条件下的强度公式。该公式将平面应变条件和常规三轴压缩条件下的抗剪强度指标联系起来,只须通过常规三轴试验,就可确定平面应变条件下的内摩擦角和黏聚力,可以有效反映中主应力对平面应变内摩擦角的提高效应。证明了等效内摩擦角达到极大值的条件即为平面应变条件下无黏性土的破坏条件,以及绘制π平面上3种准则的几何关系,阐明平面应变强度公式的数学意义。最后,通过地基承载力试验结果验证了该公式的有效性。在使用规范方法或有限元软件进行平面应变计算时,该公式可考虑中主应力的影响,选取更为有效的抗剪强度指标。     

基于虚内键模型的岩石单轴压缩全过程曲线模拟

利用虚内键模型理论（virtual internal bond model,简称VIB）,由岩石单轴压缩损伤本构方程和虚内键联接法则方程的相似性,提出了岩石单轴压缩破坏的虚内键密度演化函数 ,通过含虚内键密度演化方程的弹性张量 ,实现了岩石单轴压缩应力-应变全过程曲线的数值模拟。数值模拟结果表明,通过参数的合理选取,虚内键模型理论可以模拟不同岩样的全过程曲线。     

三峡泄滩滑坡滑动带土的改进Singh-Mitchell蠕变方程

进行了三峡泄滩滑坡体滑动带土的三轴排水蠕变试验研究,在此基础上给出了该滑动带土的 Singh-Mitchell 蠕变方程,即剪应力-应变关系采用指数函数,应变-时间关系采用幂函数来描述土体的流变性质。并根据分段拟合的思想,提出了改进的 Singh-Mitchell 蠕变方程,即在应变-时间幂函数中使用不同的幂值λ ,来分别模拟滑动带土的衰减、匀速和恒定蠕变。通过对比可知,改进 Singh-Mitchell 蠕变方程较 Singh-Mitchell 蠕变方程能更好地模拟流变性较小的滑动带粘土。     

非饱和土弹塑性模型参数的试验确定及有限元法

简要介绍广泛应用的非饱和土Alonso模型的原理与发展,并通过3组共7个试样的非饱和土三轴试验,得到模型的11个参数。在饱和土Cam本构关系的基础上,推导出一个基于非饱和土Alonso模型的应力-应变增量方程的计算公式,该方程与饱和土的本构方程形式相同。编制了能够考虑净应力和吸力二者对土体硬化规律影响的有限元程序,从而为非饱和土弹塑性计算提供一条途径,便于进一步应用于实际工程。     

三轴试验中砂土刚度的统计和回归分析及应用

砂土的密度和应力状态对其刚度有很大的影响。计算岩土工程中许多硬化土体模型都是基于邓肯−张模型得出的,没有考虑到密度对砂土刚度的影响。而在极致密或松散的砂土的三轴压缩过程中,剪切应变的上升会引起密度的显著变化。为了评估粒径分布、密度及应力状态对砂土刚度的影响,使用统计和回归方法对来自莫斯科和明斯克的15个建筑工地的962个土壤样本的各向同性三轴试验数据进行分析。基于密度和应力状态参数的影响,提出了评估不同粒径砂土刚度的经验方程。对来自欧洲、印度和美国的冲积土和陆地土试验的比较分析表明,其砂土的刚度与莫斯科和明斯克的砂土在同一范围。所提出的方程可用于初步估计有限元法计算里的刚度参数,也可应用于岩土工程模型（允许考虑刚度的变化、水平和垂直分布）。此外,还提出了基于邓肯−张模型的半经验关系。当密度的变化影响土的刚度时,该半经验关系可为受到大变形和（或）复杂加载路径影响的松散和非常致密的砂土提供更为真实的结果。一般来说,岩土工程师可将获得的结果应用于更为复杂的土体模型设计中。