基于离散元的南海软黏土剪切变形模拟

软黏土剪切变形是决定海洋构筑物、管线安全存在的重要因素,离散元能有效模拟黏土剪切变形。通过对南海神狐海域原状软黏土开展固结不排水三轴试验,得到离散元模拟所需宏观剪切强度参数指标;采用试错法调整细观参数使得双轴压缩模拟结果与室内试验相吻合,标定出一组适用于南海软黏土的细观参数组合。离散元模拟结果得到,在应变达到5.3%时,土体发生破坏,形成一条与X轴反方向为49°的剪切带;应变增加到10%时,倾角增大到57°。测量圆对剪切带内外的应变进行监测得到,剪切带内土体初始应变较大,在达到4%应变后不再继续发生变形;剪切带边缘区域,土体应变持续增大;剪切带外应变较小。一般的室内三轴试验无法得到剪切带内外土体应变变化定量数据,离散元模拟是对室内三轴试验的重要补充。     

基于离散元弹性边界模拟的砂土剪切行为分析

采用二维离散元方法开展数值双轴试验研究砂土剪切力学行为。将室内试验中橡皮膜对砂土试样的柔性包裹效应抽象为作用在试样边界颗粒上的指定法向压力,对传统离散元数值双轴试验进行了改进。改进后的离散元双轴试验能有效模拟室内试验中砂土试样剪切破坏后的侧向不均匀变形现象,并以此开展参数研究,探明了围压、颗粒摩擦因数等参数对砂土剪切行为的影响规律,最后从细观角度得到了剪切带形成过程与相关机理。研究结果表明:砂土应变软化与剪缩现象均随围压的增加而增强,剪胀效应则随围压增加而减弱;砂土变形模量随围压增大而增大,但增幅不明显;砂土残余强度基本不受围压影响,但峰值强度受围压影响明显,砂土峰值强度和对应的轴应变均随颗粒摩擦因数增大而增大,围压和颗粒摩擦因数对试样峰值强度的贡献是有限的;剪切带两侧的砂土颗粒均沿相反方向平移,接触力链主要沿加载方向分布,表现出卸载特性;剪切带内部的砂土颗粒出现相互错动,同时绕自身形心转动角位移较大,接触力链相应发生水平向偏转,出现明显的应变局部化现象,应变局部化首先发生在试样中心,此后迅速贯穿整个试样,同时砂土试样到达峰值强度。剪切带内部砂土颗粒结构疏松,颗粒间接触减弱,颗粒配位数减小,土体孔隙率增大。     

非饱和粗、细颗粒混合土在动荷载作用下的细观力学响应

为了研究非饱和粗、细颗粒混合土在动荷载作用下的细观力学性质,进行了非饱和粗、细颗粒混合土动三轴的离散元数值分析。根据室内动三轴试验结果标定PFC2D接触模型参数,进行非饱和粗、细颗粒混合土动三轴离散元模拟,分析其细观力学机理。结果表明:离散元模拟结果与室内试验结果曲线拟合度基本一致,能较好地反映非饱和粗、细颗粒混合土在动荷载作用下力学特征。相同含水率条件下,粗、细颗粒混合土的累积动应变随动应力幅值的增加而增大;随围压的增加先小幅度增大后再大幅度减小。孔隙度随动应力幅值的增大而减小;随围压的增大也减小。试样颗粒在剪切过程中,上下端运动速度最大,整体呈现对称递减状,中部运动速度最小。     

格宾网加筋土强度和变形特性的三维离散元细观模拟

为揭示格宾网加筋土强度和变形特性的细观力学本质,利用PFC3D三维离散元计算平台的FISH语言编写数值伺服程序实现大三轴试验环境和土颗粒模型生成算法,模拟大三轴素土试样应力-应变宏观响应校正土颗粒细观参数;然后,引入平行黏结构造格宾网颗粒流模型,通过模拟其室内拉伸试验校正格宾网模型细观参数;在此基础上,建立格宾网加筋土大三轴试样的三维颗粒数值模型模拟格宾网加筋土应力-应变关系,揭示格宾网加筋土颗粒之间接触力和内部剪切位移场的渐进发展规律。研究结果表明:格宾网加筋土颗粒模型数值试验能对其力学机理做出定量和定性预测;在试样两端和两层格宾网之间,分别形成"八"字形和""形倾斜剪切带;随着试样围压增加,试样内部剪切带范围缩小;加筋层数较少时,剪切滑移带范围较大;土颗粒接触力在整个格宾网截面内分布较均匀,在两层格宾网之间集中到试样截面的中部。     

苏北粉土非线性特性离散元模拟与参数分析

离散元法对颗粒系统进行动态分析,在离散颗粒单元接触分析与宏观响应间建立联系,使其成为试验手段分析细观特性的有益补充.采用三维颗粒离散元程序YADE进行了颗粒细观参数的参数分析,分析细观参数对三轴试验宏观响应的影响.研究表明：弹性模量和泊松比不仅与接触法向刚度和刚度比有关,同时随粒间摩擦角变化;而固结和剪切两个阶段粒间摩擦角控制了试样的剪胀水平和孔隙率向临界孔隙率的发展趋势.根据参数分析结果,提出了改进的三轴试验细观参数的标定步骤,并对室内试验进行模拟.模拟结果与室内试验结果较为吻合,表明该数值方法应用于三轴非线性变形特性研究具有可行性.     

基于颗粒流离散元的粘性土三轴剪切试验数值模拟

为了研究粘性土体的工程力学性质,本文通过对粘性试样进行室内三轴固结排水试验,得到了试样土体的应力—应变关系曲线以及强度指标,基于颗粒流离散元法的原理,利用PFC3D软件建立三轴试样模型,并设定颗粒间的接触为接触黏结模型,通过颗粒流数值试验对细观参数进行标定,从细观角度对粘性土的工程力学特性做了一定的研究。对比了不同围压下数值试验与室内试验的应力应变关系曲线,得到了较好的模拟效果。分析细观参数的变化对粘性土的应力—应变关系曲线的影响,得出边界条件、孔隙率、粘结强度、摩擦系数、接触刚度等细观参数的取值对粘性土的宏观力学性质有着显著的影响并对影响规律做了探讨,为今后的三轴试验颗粒流数值模拟提供了有价值的参考。     

基于颗粒流离散元的黏性土三轴剪切试验数值模拟

通过对黏性试样进行室内三轴固结排水试验,研究黏性土体的工程力学性质,得到了试样土体的应力-应变关系曲线以及强度指标,基于颗粒流离散元法的原理,利用PFC3D软件建立三轴试样模型,并设定颗粒间的接触为接触黏结模型,通过颗粒流数值试验对细观参数进行标定,从细观角度对黏性土的工程力学特性做了一定的研究。对比了不同围压下数值试验与室内试验的应力-应变关系曲线,得到了较好的模拟效果。分析细观参数的变化对黏性土的应力-应变关系曲线的影响,得出边界条件、孔隙率、黏结强度、摩擦系数、接触刚度等细观参数的取值对黏性土的宏观力学性质有着显著的影响并对影响规律做了探讨,为今后的三轴试验颗粒流数值模拟提供了有价值的参考。     

非饱和土力学行为的三维离散元分析

基于离散单元法颗粒流理论,土体颗粒单元间采用Hill接触模型来考虑非饱和土中吸力的作用,建立了非饱和土的颗粒流模型。根据室内试验结果对模型细观参数进行标定,开展了固结排水三轴剪切试验离散元数值模拟,揭示了在不同围压和吸力下,模型试样的应力-应变关系与强度特性等宏观力学性质以及细观信息平均配位数的变化规律。根据不同吸力下的莫尔-库伦破坏包线,建立了吸力与抗剪强度参数之间的关系。此外,通过标定室内试验结果中不同含水率下的应力-应变曲线,建立了Hill模型中吸力与含水率的关系。研究结果表明：Hill接触模型可以较好的模拟非饱和土的力学行为,可为非饱和土的离散元模拟提供参考。     

垃圾炉渣力学性质的颗粒流模拟

基于离散单元法理论和PFC3D软件,采用平行胶结模型来模拟炉渣中的胶结作用,建立了城市垃圾炉渣的细观颗粒模型。将模型细观参数和材料宏观力学性质联系起来,标定了模型的细观力学参数,将数值模拟结果和室内试验数据进行了对比。结果表明:离散元方法可以较好的模拟炉渣的应力-应变变形关系;细部结构调整的变化速率是影响模型受力变形机制的主要因素;改变边界刚度,分析了边界刚度对剪切强度参数和压缩模量的影响。从细观角度揭示影响炉渣强度的因素,对离散单元法的应用有一定的价值。     

基于柔性边界的非饱和土三轴试验及离散元分析

非饱和土三轴试验的离散元模拟,对于从细观角度揭示其力学特性至关重要,但目前研究多采用刚性边界条件.使用黏结球颗粒代替橡胶膜,实现围压柔性加载.在此基础上,提出非饱和土三轴数值试验模拟流程,并进行不同围压下的三轴数值试验.从宏观破坏形态和细观剪切带、配位数、接触力等方面研究不同围压下非饱和土的力学特征.结果表明：非饱和土力学特性受围压影响较大,峰值应力和弹性模量随围压增大而增大,围压越大,剪胀性越小;试样主要表现为鼓胀破坏,破坏后剪切带在形态上呈对称的“X”形,剪切带的形成开始于颗粒的转动,到达峰值应力时,内部转动颗粒发生贯通,剪切带基本成型,围压越大,剪切带上转动颗粒数目和转动角度越小;配位数先增加后减少,与体应变变化趋势一致,且随围压增大而增大;剪切使法向接触力概率密度函数峰值左移,并且使法向接触力不均匀性增大,围压越大,法向接触力分布越均匀.该研究成果可为从细观认识非饱和土的力学行为提供途径.     

温度对季冻区粉质黏土强度及变形特性的影响

为研究温度变化对天然状态及饱和状态下粉质黏土强度以及变形的影响,以典型季冻区广泛分布的粉质黏土为研究对象,通过GDS非饱和土三轴测试系统对天然状态下非饱和土以及GDS温控式静/动三轴测试系统对饱和粉质黏土,在不同温度条件下进行三轴试验。对试验结果进行研究分析可得：非饱和以及饱和粉质黏土的应力-应变曲线均呈现出应变硬化的特性。非饱和粉质黏土的黏聚力随温度降低表现出不断增加的趋势,而内摩擦角随温度降低逐渐减小,但整体变化趋势较小。饱和土的黏聚力与非饱和土变化趋势相同,随温度的降低其黏聚力逐渐增加,但相同温度时其黏聚力小于非饱和土,其内摩擦角则呈先减小后增加的转变趋势。无竖向压力作用时,相同围压条件下,非饱和与饱和粉质黏土轴向变形量随温度降低而增加,相同温度条件下,两者的轴向变形量随围压的增加会有所减小。试验成果可为季冻区粉质黏土地层工程的设计施工提供理论依据。     

土压平衡盾构改性砂土离散元模型参数反演方法研究

首次将遗传神经网络与三轴试验离散元数值模拟有机结合,用于改性砂土颗粒离散元接触模型参数反演.反演目标是使三轴试验离散元模拟曲线与真实实验曲线误差最小,采用的求解策略是基于遗传神经网络的参数识别.三轴试验的离散元数值模拟为网络提供训练样本,遗传算法映射网络输入与输出样本之间的复杂非线性关系,改性土三轴试验的真实测量曲线为参数反演提供依据.以反演结果为接触模型参数的三轴试验离散元模拟曲线与真实实验曲线相吻合,为改性土离散元接触模型参数的确定提供了有效和准确的方法,为进一步的盾构密封舱压力分析奠定了基础.     

新型倒T式防波堤结构长期变形特性预测方法研究

新型倒T式防波堤结构是可供深水环境软土地基选用的一种新型港工结构型式。深水环境中结构所受波浪力、水流力等长期荷载作用显著,但其长期变形特性尚不清楚。为探讨该类结构长期变形特性,本文基于波浪等循环荷载长期作用下软粘土变形与循环次数之间相互关系同静荷载作用下土体蠕变与时间之间相互关系相似的特性,借鉴蠕变理论,采用幂指数函数拟合饱和软粘土循环三轴试验结果,结果显示二者吻合较好,说明基于蠕变理论采用幂指数函数预测饱和软粘土长期变形的方法是可行的。进一步,利用ABAQUS二次开发平台Creep子程序实现了这一方法的有限元计算手段,建立了预测饱和软粘土长期变形的有限元实用简化计算方法,最后,采用该方法探讨了新型倒T式防波堤结构的长期变形特性。     

长期循环荷载作用下饱和软粘土应变速率的研究

为了加深对长期循环荷载作用下饱和黏土变形性状的认识,对珠江三角洲的典型淤泥质饱和软黏土进行了室内循环三轴试验,着重研究了饱和黏土的残余应变速率。通过试验结果分析可得：（1）在长期循环荷载作用下,黏土的应变速率随时间的增大而减小,应变速率对数与时间对数间的关系可用直线描述;（2）初始静偏应力与动偏应力的耦合作用以及不同的排水条件对黏土残余应变速率有较大程度的影响。这些结论可为长期循环荷载作用下饱和黏土变形机理的进一步研究提供理论基础。     

土工格栅软基路堤工程变形特征数值分析

土工格栅作为一种新型土工合成材料,具有高强度、低延伸率等特点,被用作路堤加筋材料。运用有限元通用程序AD INA对土工格栅软基路堤工程进行二维有限元分析,土体本构采用Mohr-Coulomb模型,土工格栅采用杆单元模拟,接触单元采用平面应变2-D单元。通过计算得到软基路堤工程各部位与各阶段的沉降、侧向位移变化特征,有限元计算结果与实测数据基本吻合。     

复合土钉支护的软土基坑开挖有限元模拟分析

基于平面应变假设,采用自行编制的有限元程序建立数值计算模型,考察某复合土钉支护工程下的软土基坑开挖变形性状及土钉的受力情况.结果表明,深搅桩-土钉/锚杆支护软土基坑中,锚杆张拉力的作用对基坑水平位移及受力影响较大;复合土钉支护结构并无特别针对坑底土体隆起的加固措施,须采用考虑时空效应的分层分块的开挖模式,以确保坑底隆起...     

地铁振动荷载作用下饱和软粘土性状微观研究

位于或穿越饱和软粘土层中的地铁隧道在运营之后，沿轴线产生较大的沉降量，很大程度上与饱和软粘土在地铁振动荷载作用下土体微观结构的变形和破坏、发生残余变形有关．通过扫描电镜试验对饱和软粘土颗粒的微观性状进行了研究，从土的微结构基本单元的形状和大小、基本单元体的接触状态、基本单元体问的联结形式、孔隙的形状和大小以及土体变形3个阶段的力学特性等方面进行较深入的研究．另外，将原状土和循环三轴试验破坏后两种状态的土样作了对比，从微观的角度解释了土体的变形机理．研究成果对有效控制、避免或减轻隧道经过地区地面沉降灾害的发生，指导地铁工程的设计、施工以及地铁运行后沉降的预测评价和防治等提供了有价值的参考．