颗粒离散元模拟堆积碎石土变形的参数灵敏度分析

采用颗粒离散单元方法,实现任意形状块石的模拟,基于室内试验数据标定滑坡坡体物质的堆(残、坡)积碎石土的细观参数,并考虑试样尺度效应,对影响堆积碎石土宏观变形特征的细观参数进行虚拟试验敏感度分析。研究表明,1以三轴剪切试验获得测试数据并标定颗粒离散元细观参数具有相对误差在5%以内的可靠性;2虚拟模型1(101mm×200mm)与模型2(300mm×600mm)对标定的细观参数具有尺度效应,但其相对误差控制在9%以内;3离散元颗粒间的摩擦系数与碎石土的内摩擦角、抗剪强度及残余强度为非线性正相关,摩擦系数每增加0.1,峰值主应力差平均增加118.85kPa,残余强度平均增加90.44kPa;4围压越大,材料的剪胀性越弱,围压在100kPa~500kPa时,剪胀特征值K在3~6之间变化,随着围压增加,模型破坏时的粘结力以近线性增加;5杨氏模量越大,碎石土的抗剪强度越大,但两者之间并不成线性关系,且不同杨氏模量对材料的残余强度没有显著影响。     

土体含水率与滑坡风险相关性实验研究

降雨入渗容易造成边坡滑坡自然灾害,含水率是影响土体抗剪强度的关键因素之一。本文以滑坡区红黏土为研究对象,制作了5组不同含水率的土样,并在不同垂直压力下进行快速剪切实验,得到每组对应的剪切强度参数,研究了含水率与剪切强度参数之间的关系。利用实验所得到的土工参数进行高地下水位边坡稳定性数值模拟,拟合得到边坡坡顶含水率与边坡安全系数的关系。结果表明:含水率从10%增加到25%,红黏土强度指标粘聚力、内摩擦角约分别下降50%和30%;含水率与粘聚力和内摩擦角间均呈半对数关系;边坡坡顶含水率增加导致边坡安全系数近似成二次曲线下降趋势。     

黏性土宏细观参数相关性研究

运用颗粒流程序中的bond原理结合Mohr-Coulomb破坏准则对黏性土双轴试验进行了数值模拟,探讨了黏性土宏观强度指标与其对应的细观参数之间的关系.分析结果表明,黏性土内摩擦角与颗粒摩擦角呈线性关系,黏性土黏聚力-颗粒接触强度两者拟合曲线的斜率与颗粒间摩擦系数呈幂函数关系,且两者的相关性都很好.通过该宏细观关系,得出上海第②层褐黄色粉质黏土细观参数,对其进行了双轴试验的颗粒流模拟,计算得到的应力应变曲线与实测结果一致.研究成果对于黏性土细观参数的取值及颗粒流模拟具有指导意义.     

含水率与粗料含量对砾石土力学性能影响的试验研究

诸多隧道围岩由宽级配的砾石土组成,砾石土含水率和粗料含量变化时,其力学性能存在较大差异。本文选取集呼高速金盆湾隧道浅埋段具有代表性的砾石土样本进行了室内大型直剪试验,研究了砾石土的力学性能随含水率和粗料含量的变化规律。试验结果表明:当含水率在5.37%~19.19%之间变化时,粘聚力随粗料含量的变化呈相同的变化规律,随着粗料含量增加,粘聚力先迅速下降,后显著增高,在粗料含量为45%左右时取得最小值;当含水率为10%~15%之间的某值时,细料颗粒处于最优含水率状态,当含水率超过最优含水率时,粘聚力显著下降;内摩擦角随粗料含量增加逐渐增加,随含水率增加逐渐降低;结合隧道围岩砾石土的颗粒级配特征,当含水率大于10%时,粘聚力和内摩擦角显著降低,隧道围岩稳定性也随之降低。     

不同含水率下尾矿与土工格栅界面特性的试验研究

含水率是影响尾矿加筋结构稳定性的重要因素。通过室内拉拔试验,对不同含水率(1.5%、9.7%、19.3%)下尾矿与土工格栅的界面特性进行了研究。结果表明,土工格栅的最大剪应力与法向应力呈线性相关,能够反映筋-尾矿界面参数似粘聚力和似摩擦角;似粘聚力随含水率增大先增大后减小,似摩擦角随含水率增大先减小后增大;在含水率增大过程中存在界限含水率,大小约为10%;在界限含水率下,似粘聚力达到最大值,似摩擦角达到最小值;似摩擦系数随法向应力增大呈非线性递减,可以二次抛物线拟合反映二者关系;在不同含水率条件下,似摩擦系数随含水率增大呈现先增大后减小规律,与对应界面参数似粘聚力的变化一致。上述试验分析成果对应用土工格栅进行尾矿坝加筋设计有一定的指导价值。     

不同含水率和压实度铅锌尾砂抗剪强度试验研究

尾砂的抗剪强度是分析尾矿坝失稳问题的重要参数。为了探究含水率和压实度对铅锌尾砂抗剪强度的影响,开展三种压实度(80%、90%、100%)情况下五组不同含水率(3%、6%、9%、12%、15%)铅锌尾砂试样的抗剪强度试验。结果表明:铅锌尾砂的剪应力随压实度的增大而增大,随含水率的增大先增大后减小;随着含水率的增大,铅锌尾砂的内摩擦角先减小后增大,粘聚力先增大后减小;在含水率为6%时,粘聚力最大;随着压实度的增大,铅锌尾砂内摩擦角和粘聚力均呈线性增长,含水率为6%时压实度对粘聚力的影响最大,含水率超过12%,压实度对粘聚力的影响逐渐减小;保持尾砂含水率在6%左右和提高尾矿坝的压实度,能够有效提高尾砂抗剪强度,降低尾矿坝溃坝危险性。研究结果对尾矿坝的安全管理和稳定性评价具有一定意义。     

金沙江中游凝灰岩强度参数综合识别研究

凝灰岩是一种典型的分布不规律、易崩解、难取样火山碎屑岩,其强度参数难以准确评估。本文依托华丽高速公路金沙江特大桥工程,判定丽江岸为反倾坡,为陡倾结构面或卸荷裂隙带剪入、缓倾结构面剪出组合滑移模式,破坏规模有限;华坪岸为顺层坡,为陡倾结构面或卸荷裂隙带剪入、凝灰岩层主滑、缓倾结构面剪出组合滑移模式,破坏规模较大。T3和T2凝灰岩软弱层为关键层,其强度参数对岸坡稳定性影响很大。针对此种情况提出了经验、试验、数值反分析相结合综合识别凝灰岩强度参数的方法。利用这种方法,在工程岩体分级标准GBT50218-2014、水利水电岩石力学参数手册、室内外试验、数值反分析的基础上,综合确定了凝灰岩带泥化参数为内摩擦角24°、粘聚力50k Pa、泥化率60%;凝灰岩层内摩擦角29°、粘聚力90k Pa,在工程经验取值的范围之内;玄武岩内摩擦角36°、粘聚力600k Pa,在水利水电工程岩石力学手册推荐值范围之内;IV陡倾结构面强度参数内摩擦角35°、粘聚力150k Pa,符合一般结构面试验经验值。经计算复核,埋深浅、局部分布的T3凝灰岩主滑模式下安全系数为1.65~2.15;埋深大、分布连续的T2凝灰岩主滑的安全系数为1.32~1.67,T2主控下的安全系数要比局部分布T3情况小约0.33~0.48,平均降低约21%。虽然T2主控下的最不利工况安全系数也大于1.30,但是在实际施工和运行过程中仍需关注T2凝灰岩层可能形成的深大滑移问题。类似和临近工程可参考该凝灰岩强度参数或强度参数综合识别方法。     

PFC2D平节理模型细观参数标定方法

细观参数的快速准确确定对PFC正确模拟岩石力学性能非常重要.目前,平节理模型细观参数的标定研究均未考虑对泊松比和压拉比有很大影响的黏结比例φb和影响岩石强度及裂纹产生难易程度的黏结强度变异系数Rsd这两个细观参数,且未考虑作为判断围岩损伤范围重要指标的宏观参数起裂强度σci在参数中新增φb,Rsd和σci,采用正交数值试验方法研究了平节理模型8个细观参数与岩石6个宏观参数之间的相关性,确定了各宏观参数与主要细观参数之间的拟合关系,并分析了宏细观参数间的趋势关系.在此基础上,提出了平节理模型细观参数的标定流程,然后根据片麻花岗岩物理试验的宏观力学参数确定了模型的细观参数并进行数值试验,模拟结果与物理试验结果基本吻合,说明标定流程可行.     

参变量变分原理求解土的变形模量与压缩模量间的关系

运用参变量最小势能原理推导了弹性及塑性状态下的土的变形模量与压缩模量间的关系,认为两者比值与土的内摩擦角Ф、泊松比μ及粘聚力c等土性指标有关,提出了塑性影响因子a及临界应力p^*的概念,指出当试验时土体应力大于土体的临界应力值时,变形模量与压缩模量的理论关系式应引入塑性影响因子a,从而克服了以往基于弹性理论得出的变形模量理论值远小于试验值的不足。     

临近边坡的基础的极限承载力

本文推导了了靠近边坡的基础的极限承载力计算公式,它由地基内摩擦角,粘聚力,土壤比生,以及几何参数：基础宽,坡顶距和坡角等完全决定。     

爆炸产物发生化学反应对药室压力影响的研究

研究了爆炸产物发生化学反应对约束爆炸准静态压力的影响,推导了考虑爆炸产物发生化学反应时,约束空间内准静态压力的计算公式.以TNT炸药发生约束爆炸为例,分析了不同药量体积比下爆炸产物的化学反应动力学过程.编写了计算TNT炸药约束爆炸后准静态压力的程序,将程序计算结果与实验测量结果和ConWep程序计算结果进行对比,验证了计算方法的准确性,证明了爆炸产物发生化学反应对约束爆炸准静态压力有重要的影响.考虑爆炸产物发生化学反应对压力的影响,利用LS-DYNA对TNT炸药在一长方体药室内的爆炸进行了数值模拟,模拟结果表明,考虑爆炸产物发生化学反应时,药室会出现局部塑性应变与应力集中现象,最大塑性应变达到0.8‰.     

掺气分流墩墩头的动水压强特性研究

对掺气分流墩墩头的动水压强特性进行了试验研究,以求了解墩头的抗空蚀性能及动水压强特性。试验结果得出：时均压强P／Pman沿墩高和Cp沿横断面的分布规律;墩头轴线与来流有夹角时Cp沿横断面的分布规律;墩头壁压脉动压强特性,在墩后充分供气时,墩头的空化特性。     

孔隙度对PBX炸药冲击起爆影响的实验研究

为了研究孔隙度(装药密度)对PBX炸药冲击起爆爆轰成长的影响,采用炸药冲击起爆锰铜压阻一维拉格朗日实验测试系统,测量了不同孔隙度的PBXC03炸药(HMX的质量分数为87%,TATB的质量分数为7%,黏结剂的质量分数为6%)冲击起爆过程不同拉格朗日位置的压力-时间历史。结果显示:在本文装药范围和加载条件下,孔隙度对PBX炸药冲击起爆爆轰过程的影响不单调,中等密度的炸药冲击起爆和爆轰成长最快,这是热点点火过程与燃烧反应过程共同作用的结果。     

缺陷对石墨烯摩擦性能影响的分子动力学研究

采用分子动力学方法模拟硅探针在空位缺陷和Stone-Wales(SW)缺陷石墨烯上的滑移过程,研究空位缺陷和SW缺陷对石墨烯摩擦力的影响.研究结果表明:两种缺陷石墨烯摩擦力大于完美石墨烯,空位缺陷使石墨烯界面势垒增大导致能量耗散增加,摩擦力增大;SW缺陷使石墨烯表面形成凸起,阻碍探针滑移,摩擦力增大.空位缺陷石墨烯平均摩擦力随缺陷浓度的增加而增加,Y向空位缺陷石墨烯平均摩擦力大于X向,这都是由空位陷处能量势垒和缺陷与探针切向作用距离共同决定的.SW2型缺陷石墨烯摩擦力大于SW1型,X向SW2型缺陷石墨烯摩擦力大于Y向SW2型,因为存在相邻五边形碳原子环结构的石墨烯表面更容易产生凸起,摩擦力较大.以上研究结果完善了缺陷石墨烯的摩擦机制,对设计和开发石墨烯微纳器件提够了理论依据和指导.     

基于NPF-CKF的捷联惯导系统动基座初始对准技术

当海况不佳时,水下航行器大幅晃动,捷联惯导系统无法快速完成自主初始对准,因此提出了利用多普勒计程仪提供的速度信息进行运动中辅助对准。针对在非线性对准中扩展卡尔曼滤波存在精度低,且需要计算雅可比矩阵等不足,提出了一种基于非线性预测滤波的求容积卡尔曼滤波算法。该滤波算法将惯性器件测量误差作为模型误差使用非线性预测滤波器进行实时预测,然后再利用求容积卡尔曼滤波对模型误差补偿后的系统进行状态估计。仿真结果表明,与扩展卡尔曼滤波和求容积卡尔曼滤波算法相比,该滤波算法能够不仅提高失准角特别是方位失准角的估计精度,其精度约为45″,而且加快了收敛速度。同时由于该滤波算法降低了系统状态的维数,因此也大大减少了计算量。     

二烷基二硫代磷酸锌在铁表面等温吸附的研究

本文对二烷基二硫化磷酸锌在不同温度下自液体石蜡溶液中于铁表面的等温吸附进行了试验研究。结果表明，在２８５－３５３Ｋ的温度范围内，恒温搅拌３小时就基本达到了吸附平衡。在２８５Ｋ和３０３Ｋ这两种温度下的吸附等温线均属Ｓ型曲线，这表明二烷基二硫代磷酸锌主要是以分子形态和表面反胶的形式被吸附；而其在较高温度（３３３Ｋ和３５３Ｋ）下的吸附等温线却发生了显著的变化，这是二烷基二硫代磷酸锌在铁表面催化分解和氧化     

OSTR试验下感度分布参数的推断

研究了OSTR试验数据的统计处理方法。给出了推断感度分布参数的两种方法 :条件推断与非条件推断方法。通过研究与模拟 ,得出条件推断方法适用于推断机械撞击作用火工品的一些R响应点的结论。     

基于三维数字图像相关的三维形貌拼接研究

三维数字图像相关(DIC)方法已在众多工程技术领域中得到了日益广泛的应用.当该方法用来测试表面曲率较大的物体的变形时,受双相机之间夹角的限制,不能实现全场的三维形貌测量.采用在两个(及以上)不同位置拍摄,分别得到被测物的局部形貌,然后用三维DIC算法计算局部坐标系之间的坐标变换关系,将局部形貌转换至同一坐标系下并拼接得到被测物全场三维形貌.实验结果验证了该方法的可行性和可靠性.     

基于LMI的框架结构分数阶控制方法研究

以框架结构为研究对象,采用分数阶状态方程进行控制器设计,提出了基于LMI的分数阶控制器设计方法。首先,把原系统振动方程改为阶次为1～2的分数阶状态方程,同时,构造同阶次分数阶控制器状态方程,进而组装为一个整体状态方程;其次,根据分数阶系统稳定性条件确定满足系统渐进稳定性的不等式。为了获得可行解,在不等式中附加了两个参数,同时,令控制器参数矩阵为对称;再次,从仿真结果来看,控制效果过好,但控制力偏大,为此,控制力计算时附加了调节因子,以期在满足控制效果的基础上,降低所需控制力;最后给出了一个算例说明本文方法的可行性。     

Spacing effects on hydrodynamics of heave plates on offshore structures

The nonlinear viscous flow problem associated with a heaving vertical cylinder with two heave plates in the form of two circular disks attached is solved using a finite difference method. Numerical experiments are carried out to investigate the spanwise length effects on the hydrodynamic properties, such as added mass and damping coefficients. Over a Keulegan–Carpenter (KC) number range from 0.1 to 0.6 at frequency parameter (β=7.869×10⁷), calculations indicate that a KC-dependent critical value of spanwise length L/D_d (L—the distance between the disks, and D_d—the diameter of the disks) exists. A significant influence of L/D_d on the vortex shedding patterns around the disks and the hydrodynamic coefficients is revealed when L/D_d is smaller than the critical value due to strong interaction between vortices of different disks. Beyond that limit, however, both added mass and damping coefficients are shown to be rather stable, indicating the independence of the vortex shedding processes of different disks.