一种有限元转化为颗粒离散元的方法及其应用研究     

冯春  李世海  刘晓宇 《岩土力学》2015,36(4):1027-1034

为了充分发挥有限元与颗粒离散元各自的优势,提出了一种由有限元转化为颗粒流的方法。数值模型首先用较粗的有限元网格进行离散,并在单元上引入连续介质本构模型。力学计算开始后,实时跟踪各单元的应力状态。一旦某单元的应力满足Mohr-Coulomb准则或最大拉应力准则,删除该单元,同时创建具有一定数目、随机分布且微嵌入的颗粒簇。其后,该单元所在区域的非连续变形及失稳断裂由颗粒簇演化获得。各颗粒的质量、材料参数、速度、位移、接触力等信息根据插值从有限元单元中继承。为了实现有限元与颗粒流接触面的耦合计算,引入了点-棱（二维）及点-面（三维）接触模型,通过法向及切向弹簧实现接触力的计算。颗粒球与有限元板的碰撞分析、单轴压缩、岩石切割等案例展示了上述方法的精确性及合理性。  相似文献   

基于岩土材料软化特性的滑坡多级滑动面分析方法研究     

薛海斌  党发宁  尹小涛  雷曼  杨超 《岩土力学》2015,36(11):3235-3242

在实际工程中,滑坡经常呈现出多级破坏的特征,而一般计算和设计中仅关注最危险的滑动面及对应的最小安全系数,这往往遗留安全隐患。在考虑岩土材料软化特性的基础上,借助FLAC3D和Matlab软件平台构建了无需人工干预便可有效模拟与评价滑坡多级破坏的理论框架。以旬阳县党家坝廉租房小区滑坡为例,通过对塑性剪应变、塑性拉应变、剪切应变增量等特征变量的渐进发展规律分析发现：滑坡多级滑动面的形成,其时间和空间顺序并不一定相对应。各级滑动面在时间上按产生的先后顺序依次为第1级主滑面、第2级主滑面、次级滑面;而在空间上按照从前到后的顺序依次为第1级主滑面、次级滑面、第2级主滑面。从计算所得的滑坡最终破坏形态中发现,坡体内部出现的滑动面条数与现场采集的拉裂缝数基本保持一致;第1级主滑面的入口位置与滑坡前缘拉裂缝位置基本吻合;而第2级主滑面及次级滑面的位置与现场勘查到的拉裂缝位置出现偏差。从材料参数的渐进发展规律中发现,滑面上强度参数的分布及大小均随着滑面的产生、发展而逐渐变化,其变化区间为峰值强度到残余强度,这是此方法可以有效模拟多级滑动面形成过程与各级滑动面之间相互影响的核心所在。借助矢量和法成功地实现了基于滑面上强度参数渐进发展规律的安全系数演化过程的确定。通过对安全系数演化过程的分析发现,在滑坡多级滑动面的形成过程中安全系数大小顺序呈现出3种不同的状态,这很好地揭示了各级滑动面在形成过程中的主、被动关系。  相似文献   

基于数字图像测量系统的饱和细砂渐进破坏特性研究     

高军程  郭莹  贾金青  涂兵雄 《岩土力学》2016,37(5):1343-1350

基于亚像素角点检测的数字图像测量系统能够记录三轴试样表面方块角点的位移,从而获得试样表面每一时刻的局部应变及应变场。分析了试样不同位置的局部径向应变在剪切过程中的变化,获得了以不同初始成样含水率制备的松砂和密砂试样在不同特征时刻的应力与应变的特征值。通过轴向应变场的变化分析了从应变局部化出现到剪切带发育、形成的这一完整的渐进破坏过程,总结了剪切带形成时的局部最大轴向应变特征值,并定性地分析了剪切带内、外土体在渐进破坏过程中不同的轴向应变增长率。试验结果表明：在应变场中试样应变局部化明显,并可以依此确定应变局部化的出现、剪切带的形成;对于密砂及初始含水率为0%制备的松砂试样,应力在应变局部化出现之初即达到峰值,剪切带形成时应力已经开始下降,进入应变软化阶段;以初始含水率为6%和12%制备的松砂试样在达到应力峰值时剪切带已经形成;剪切带内土体的局部轴向应变增长幅度比剪切带外的土体大得多。试样整体轴向应变的增大主要是由剪切带内土体剪切破坏产生的较大轴向应变所致。  相似文献   

干密度对粗粒料力学特性的影响     

姜景山  程展林  左永振  丁红顺 《岩土力学》2018,39(2):507-514

粗粒料的力学特性不仅取决于应力状态,还与粗粒料本身的松密程度密切相关,即与粗粒料的材料状态相关。通过不同初始干密度的粗粒料常规大型三轴各向等压固结排水剪切试验研究了干密度对粗粒料力学特性的影响。试验结果表明：粗粒料应力-应变曲线的形态取决于密度和围压的共同作用,破坏状态之前密度和围压共同决定了粗粒料的强度和变形,而当应变足够大处于渐进状态或临界状态时,粗粒料的应力和体积应变受初始干密度的影响逐渐减小直至消失。对于软化型曲线,相变状态的应力小于渐进状态的应力,而渐进状态的应力小于破坏状态的应力。无论围压大小如何,随着初始干密度的增大,应力-应变曲线的硬化性逐渐减弱,而软化性则逐渐增强。  相似文献   

水压影响岩石渐进破裂过程的试验研究   总被引：2，自引：0，他引：2  

彭俊  荣冠  周创兵  王小江  侯迪 《岩土力学》2013,34(4):941-946

岩石的压缩过程伴随着裂纹扩展,在低、中围压条件下扩展的裂纹主要受到张拉作用而产生,张拉破坏是首要的作用机制。研究发现,岩石的渐进破裂过程受岩石的结构和构造影响,比如矿物成分、颗粒大小以及胶结情况等,另外,围压以及开挖扰动等外界因素对岩石的渐进破裂过程也有重要影响。基于试验方法探讨水压对岩石渐进破坏过程的影响,首先,阐述岩石渐进破裂过程各个阶段的特征;然后,对岩石渐进破裂指标--岩石的启裂强度 和损伤强度 的确定方法进行总结;最后,选取细粒的石英砂岩为研究对象,通过试验研究水压对岩石渐进破坏过程的影响。试验结果表明,相同围压条件下,随着岩样两端水压的增大,岩石的 有逐渐增大的趋势,而岩石的 和峰值强度 逐渐变小;随着围压逐渐增大,岩石的启裂强度 、损伤强度 及峰值强度 均逐渐增大。  相似文献   

地层渐进成拱对浅埋隧道上覆土压力影响研究     

汪大海  贺少辉  刘夏冰  张嘉文  姚文博 《岩土力学》2019,40(6):2311-2322

砂土及完整性较差、黏聚强度较小的破碎岩体中,浅埋隧道地层拱作用机制随地层变形发展而变化,受此影响隧道松动土压力也相应变化。常规方法忽略了地层拱不同阶段力学机制的不同,同时未考虑剪切面转动与大主应力旋转之间的相互关系,因此,不能解决浅埋隧道地层能否成拱、地层拱何时贯通至地表以及地层拱发展对隧道松动土压力影响等问题。针对这一情况,提出渐进地层拱力学模型以反映不同阶段地层拱的力学机制;其次,同时考虑主应力旋转、剪切面转动及二者相互关系,确定拱内土体应力分布;随后,优化了传统主应力偏转与地层差异沉降间的数学模型。在此基础上确定了渐进地层拱对隧道松动土压力的影响。改进方法结果与传统方法结果及试验结果的对比验证了改进方法的有效性与适用性。通过参数分析研究了隧道初始松动压力、随地层变形发展的松动压力以及地层拱贯通至地表时的极限变形等关键参数。最后,对下北山超大跨浅埋隧道的研究说明了改进方法的实用性。结论显示：（1）初始松动压力为初始松动区内土体重力,初始松动区范围不受覆跨比影响,而受地层强度影响,随内摩擦角增加而减小;（2）最大拱效应阶段以后,松动土压力随地层变形发展而增加,深埋、小跨度隧道（ ）增长速率较慢,反之较快;（3）极限变形随覆跨比、内摩擦角增加而增加,深埋、小跨度隧道地层拱效应更明显;（4）对于下北山隧道,初始地层拱存在,初始松动土压力为0.37 ,极限松动土压力为0.41 ,最终松动土压力为0.54 ,隧道变形应控制在5.7%以下避免地层拱贯通至地表。  相似文献   

考虑张拉-剪切渐进破坏的边坡矢量和法研究     

王伟  陈国庆  郑水全  张广泽  王栋 《岩土力学》2019,40(Z1):468-476

针对边坡渐进破坏过程中岩土体强度参数不断劣化现象,采用应变软化模型代替传统的理想塑性模型,提出一种考虑张拉?剪切渐进破坏的边坡矢量和分析方法。从安全系数表达式与最危险滑动面搜索两方面进行研究,首先基于滑动和力的矢量特性,改进矢量和法滑动趋势与安全系数表达式,进而基于改进自适应遗传算法搜索边坡最危险滑动面。算例表明,采用应变软化模型时的矢量和法安全系数与最危险滑动面形态具有动态变化与调整的特性。该基于张拉-剪切应变软化模型的计算方法进一步扩展了边坡矢量和分析法。  相似文献   

一种考虑雪粒径变化的积雪面积反演算法     

王杰  黄春林  郝晓华 《地球信息科学学报》2017,19(1):101-109

光谱混合分析能够提取亚像元信息,被广泛地应用于遥感影像目标探测之中。本文针对MODIS积雪遥感影像,基于光谱混合分析框架,利用渐进辐射传输模型建立不同粒径大小的雪反射率光谱库,提出了一种考虑端元变化及二次辐射的雪盖面积反演算法。此算法首先利用渐进辐射传输模型建立不同粒径大小积雪的反射率光谱库,然后使用序贯最大角凸锥方法获取植被、土壤与岩石、阴影的光谱库。在建立各种地物反射率光谱库之后,利用均方根误差最小的方法获取最优端元组合。在此基础上,考虑端元独立辐射以及积雪与其它地物的二次辐射过程,利用稀疏光谱混合模型获取积雪面积与雪粒径大小。实验结果表明：此方法能够同时反演雪粒径与积雪面积,反演的雪粒径相比单波段的渐进辐射传输模型小,反演的积雪面积相比MOD10A1产品精度略微提高。  相似文献   

水泥土桩加固海洋土的渐进变形破坏试验与数值分析     

周成  殷建华  房震 《岩土力学》2005,26(Z1):205-208

应用一个高度非线性的三维各向异性粘弹塑性(EVP)模型及有限元法,对一水泥土桩(CDM)加固海洋土的复合地基进行数值模拟。海洋软粘土的渐进变形采用粘弹塑性模型来模拟,水泥土桩破坏前后的应力-应变关系采用一个双曲线模型来模拟,提出了一个简单且非常实用的数值分析方法来解决水泥土桩碎裂软化的数值模拟问题。将数值模拟结果与实测数据进行了比较,二者吻合较好。  相似文献   

RT模式下刚性挡墙土压力计算方法研究   总被引：5，自引：0，他引：5  

龚慈  魏纲  徐日庆 《岩土力学》2006,27(9):1588-1592

针对绕墙顶向外转动的刚性挡土墙,提出一种土压力计算方法。根据土体渐进破坏机理,考虑土拱效应,建立了填土内摩擦角及墙土接触面上外摩擦角的发挥与土体位移的非线性关系,并根据初始应力条件确定初始内摩擦角。采用改进的水平层分析法计算各转角下的土压力分布,并得到土压力合力大小及其作用点的计算公式。通过比较,不同转角下土压力强度、合力大小以及作用点计算值与模型试验实测结果接近。  相似文献