The subsidence by mining car traveling on deep-sea soft bottom based on Burger's creep model
-
摘要:
在陆地矿产资源日渐枯竭的今天,深海矿产资源已成为全球各个国家争相开采与利用的焦点,深海采矿车是实现深海矿产资源开采的重要装备。海底稀软底质是一种承载力与抗剪强度极低的特殊底质,在采矿作业中,深海稀软底质的物理力学特性直接影响采矿车行走的稳定性。文章选取Burger’s接触模型作为深海稀软底质的本构模型,对某海域海底稀软原状土开展室内三轴试验,通过PFC3D颗粒流数值模拟实验对比实际三轴试验,对稀软底质的Burger’s蠕变模型进行参数标定,同时依据标定结果改变相应参数,针对5种不同底质条件的工况,建立海底采矿车的数字仿真模型,模拟各工况下采矿车在不同行驶速度时的下陷深度。结果显示,下陷深度会随行驶速度呈非线性变化,在一定范围内随着行驶速度的增大而减少并逐渐趋于稳定。同时结果还表明,该区域海底稀软底质具有更高的黏粒含量(38.1%~48.4%)、含水率( 88.13%~137.79%)和压缩性(压缩系数:1.86~3.73 MPa−1,压缩模量:1.26 ~2.13 MPa),具有更低的密度(1.3 ~1.5 g/cm3)和强度特性(贯入阻力:0.19 ~1.32 N,黏聚力:3.7~6.9 kPa,内摩擦角:2.4°~3.9°),即承载力较低,蠕变性能较强。本研究在宏观上做了一般的探讨,为类似参数的稀软底质下海底采矿车的运行安全控制提供了较好借鉴与依据。
-
关键词:
- 海底采矿车 /
- 海底稀软底质 /
- Burger’s蠕变模型 /
- PFC3D颗粒流
Abstract:With the depletion of terrestrial mineral resources, the deep-sea mineral resources have become the focus of exploitation and utilization in all countries of the world. A deep-sea mining car is an important equipment for deep-sea mineral resources mining, and a soft marine sediment is a special substrate with very low bearing capacity and shear strength. In mining operations, the physical and mechanical properties of soft marine sediment directly affect the stability of mining vehicles. The Burger's contact model was selected as the constitutive model of deep-sea soft sediment, and a laboratory triaxial test was carried out on seafloor soft undisturbed soil in a certain area. By comparing the actual triaxial test with PFC3D particle flow numerical simulation experiment, parameters of the Burger's creep model of soft sediment were calibrated. Meanwhile, the corresponding parameters were modified according to the calibration results, the digital simulation model of the seabed mining vehicle was established under five working conditions with different soil beds, and the subsidence depth of the mining car under different driving speeds under each working condition was simulated. The results show that the subsidence depth changes nonlinearly with the driving speed, and the subsidence depth decreases with the increasing of the speed in a certain range and gradually tends to be stable. In addition, the soft sediment in the creep area features higher clay content (38.1%~48.4%), water content (88.13%~137.79%), and compressibility (compression coefficient: 1.86~3.73 MPa−1, compression modulus: 1.26~2.13 MPa), and lower density (1.3~1.5 g/cm3) and strength (penetration resistance: 0.19~1.32 N, cohesion: 3.7~6.9 kPa, internal friction angle: 2.4°~3.9°), indicating that the bearing capacity is low and the creep performance is strong. This study provided a reference on a macro level and a theoretical basis for safe operation of seabed mining vehicle on soft bottom with similar parameters.
-
随着对深海的不断探索,丰富的海洋矿产资源成为科学家关注的重点。遍布于各大洋底的多金属结核是锰、铜、镍、钴等多种金属的重要来源,研究表明其具有巨大的储量规模和经济价值[1]。据统计,全球大洋约15%的面积被多金属结核所覆盖[2],而太平洋C-C区、西北太平洋海山盆地、中印度洋海盆等多金属结核富集区大多赋存于低承载力与抗剪强度的稀软底质中[3-5]。伴随着南海海洋勘探深入,我国在南海多个海域发现了多金属结核成矿区,大量陆源物质进入海洋也为南海北部边缘海结核提供了丰富的成矿物质来源,便于南海边缘海结核的快速生长成矿,作为极具潜力的海洋矿产资源,南海北部多金属结核成矿区成为试验目标靶区,是目前国内深水资源勘探、开发的重点区域[6]。但铁锰结核主要分布在海底有松软沉积物的表面或者向下几厘米的表层沉积物中,开采起来较为困难[7]。深海采矿车是实现深海矿产资源开采的重要环节,由于深海沉积物不同于陆地底质,是一种具有高含水率、高液限、高塑性、大孔隙比、低强度等特性的特殊稀软底质[8],在采矿作业中,深海稀软底质的物理力学特性极易影响采矿车行走的稳定性,较低的承载力会使得采矿车具有下陷的风险,导致采矿车不能行走或者丧失采矿功能,严重时可能会使整个采矿系统瘫痪。因此,研究海底稀软底质的物理力学特性对深海采矿车行驶下陷的影响具有重要意义,可为未来实现深海矿产的开采和利用提供必要的经验积累和技术储备[9],同时为类似稀软底质区提供借鉴与参考。
目前,已有学者针对深海采矿车进行了广泛的研究。如Leng等[10]研究了深海采矿车的动态特征,包括动态建模方法、轨迹特征和采矿采集系统。讨论了考虑深海沉积物对深海采矿车性能的相互作用。Schulte等[11]通过一定比例的膨润土与水模拟海底地质的物理力学特性,研究该模拟底质在不同剪切装置下的剪切应力-位移关系,分析了采矿车深海行走性能,计算了细软底质可提供的牵引力以及采矿车的静态与动态沉陷。Hong等[12]建立了履带式采矿车在海底稀软底质上行走的3D瞬态模型,通过压力-沉陷、剪切应力-位移和剪切应力-动态沉陷等关系式描述了稀软底质对履带车的作用力,对履带式采矿车直线行驶与转向行驶性能进行了动力学分析。戴瑜等[13]通过构建海底履带车行驶牵引性能的数学计算模型,探明了履带车海底行走牵引力与滑转率之间的定量关系,研究了履带车接地段参数对其牵引力的影响特性,提出了履带式采矿车海底行走最佳滑转率、滑转率控制范围及履带接地参数的最佳设计值。陈峰等[14]考虑了海底履带机器车的特殊工作环境,以及车体的特殊设计,建立了海底履带机器车的数学模型,使用了MADYMO多刚体动力学仿真软件,建立了深海底履带机器车结构模型,进行了仿真分析,并与实验结果进行了比较,验证了模型的有效性。Qi等[15]建立了描述海底土壤流变特性的分数阶导数Burger’s模型,并通过太平洋C-C矿区模拟土壤蠕变试验确定了其关键参数。预测了DSMV的履带接触地面长度和行走速度,为提高采矿效率提供了安全设计。
深海稀软底质具有多孔性、高含水率和低剪切强度等力学特性且呈流动态,深海采矿车在其床面行走,容易产生随时间变化的沉陷现象,因此,研究中应考虑深海稀软底质的蠕变特性。此外,土体所表现出的物理力学性质,其根本都是内在结构发生变化的外在反映,结构性对土的性质的变化起着重要的决定作用,缺乏对土体结构特征及其演化规律的描述会导致实际的工程计算结果难以模拟土体的实际状态,进而出现较大的计算误差,故而原状土的结构性对研究结果具有重要意义[16]。综上,前人的研究要么没有考虑海底稀软底质的蠕变特性,要么仅采用模拟土样来模拟海底底质的物理性质,缺乏原状土的结构性,因此均具有一定的局限性。
考虑到深海稀软底质的蠕变特性,Burger’s模型作为常用的蠕变模型之一,可反映稀软底质的黏弹性特点。因此选用PFC3D颗粒流数值模拟软件为研究手段,以其中Burger’s接触模型作为海底稀软底质的本构模型,进行模拟研究。软件中模型接触力按时间步长不断进行迭代,模拟应力变形随时间的发展过程,从而迭代计算稀软底质的非恒定荷载作用。以某海域稀软底质原状土为依托,开展土工试验与室内三轴试验,获取其物理力学性质,并通过室内三轴试验与模拟三轴试验对PFC3D中Burger’s接触模型进行参数标定。依据标定结果,改变相应参数从而获得不同底质条件的5种工况,针对5种工况建立了采矿车在海底稀软底质地层上行走的数字仿真简化模型,通过计算模型研究了各工况下海底采矿车在不同行驶速度时的下陷深度,从而分析海底采矿车在稀软底质的行驶稳定性,在宏观上做一般的探讨,为类似参数的海底稀软底质的采矿车行驶提供借鉴与参考。
1. 研究方法
1.1 Burger’s蠕变本构模型
对海底稀软底质选用Burger’s蠕变模型,通过开尔文模型(Kelvin)和麦克斯韦模型(Maxwell)在法向和剪切方向串联来模拟蠕变机制。开尔文模型是线性弹簧和阻尼元件并联作用的组合,麦克斯韦模型是线性弹簧和阻尼元件串联作用的组合。伯格的模型作用于一个很小的区域,因此只传递一种力。
在法向和剪切方向上,Burger’s模型的蠕变组分如图1所示,该模型将开尔文模型和麦克斯韦模型串联起来。
在法线方向上,麦克斯韦模型提供了刚度为$ {E}_{{\mathrm{m}}_{n}} $的线性弹簧和黏性系数为$ {\eta }_{{\mathrm{m}}_{n}} $的阻尼器,开尔文模型提供了刚度为$ {E}_{{\mathrm{k}}_{n}} $的线性弹簧和黏性系数为$ {\eta }_{{\mathrm{k}}_{n}} $的阻尼器。在剪切方向上,麦克斯韦模型提供了一个刚度为$ {E}_{{\mathrm{m}}_{s}} $的线性弹簧和一个黏性系数为$ {\eta }_{{\mathrm{m}}_{s}} $的阻尼器,开尔文模型提供了一个刚度为$ {E}_{{\mathrm{k}}_{s}} $的线性弹簧和一个黏性系数为${\eta }_{{\mathrm{k}}_{n}}$的阻尼器。基于该本构模型,石崇等[17]在《颗粒流(PFC5.0)数值模拟技术及应用》中给出颗粒间接触力。
在每个方向上,Burger’s模型的总位移u为模型的开尔文部分位移$ {u}_{\mathrm{k}} $和麦克斯韦部分位移$ {u}_{{\mathrm{m}}_{E}} $、$ {u}_{{\mathrm{m}}_{\eta }} $三者之和,如公式1所示:
$$ u={u}_{\mathrm{k}}+{u}_{{\mathrm{m}}_{E}}+{u}_{{\mathrm{m}}_{\eta }} $$ (1) 式1关于时间t的一、二阶导数分别如公式2、公式3所示:
$$ \dot{u}={\dot{u}_{\mathrm{k}}}+{\dot{u}_{{\mathrm{m}}_{E}}}+{\dot{u}_{{\mathrm{m}}_{\eta }}} $$ (2) $$ \ddot{u}={\ddot{u}_{\mathrm{k}}}+{\ddot{u}_{{\mathrm{m}}_{E}}}+{\ddot{u}_{{\mathrm{m}}_{\eta }}} $$ (3) 式中,$ \dot{u} $为位移关于时间t的一阶导数;$ \ddot{u} $为位移关于时间t的二阶导数。
接触力F,利用模型中开尔文体的刚度$ {E}_{\mathrm{k}} $、黏性系数$ {\eta }_{\mathrm{k}} $和一阶导数,由公式4给出,同时,其一阶导数$ \dot{F} $如公式5所示:
$$ F=\pm {E}_{\mathrm{k}}{u}_{\mathrm{k}}\pm {\eta }_{\mathrm{k}}{\dot{u}_{\mathrm{k}}} $$ (4) $$ \dot{F}=\pm {E}_{\mathrm{k}}{\dot{u}_{\mathrm{k}}}\pm {\eta }_{\mathrm{k}}{\ddot{u}_{\mathrm{k}}} $$ (5) 式中,符号$ \pm $分别表示法线方向和切线方向,$ + $表示法线方向,$ - $表示切线方向。同样,利用麦克斯韦体的刚度$ {E}_{\mathrm{m}} $和黏性系数$ {\eta }_{\mathrm{m}} $,给出接触力F及其一、二阶导数$ \dot{F} $、$ \ddot{F} $:
$$ F=\pm {E}_{\mathrm{m}}{u}_{{\mathrm{m}}_{E}} $$ $$ \dot{F}=\pm {E}_{\mathrm{m}}{\dot{u}_{{\mathrm{m}}_{E}}} $$ $$ \ddot{F}=\pm {E}_{\mathrm{m}}{\ddot{u}_{{\mathrm{m}}_{E}}} $$ (6) $$ F=\pm {\eta }_{\mathrm{m}}{\dot{u}_{{\mathrm{m}}_{\eta }}} $$ $$ \dot{F}=\pm {\eta }_{\mathrm{m}}{\ddot{u}_{{\mathrm{m}}_{\eta }}} $$ 综上,可求得接触力F的二阶微分方程如公式7:
$$ F+\left(\frac{{\eta }_{\mathrm{k}}}{{E}_{\mathrm{k}}}+\frac{{\eta }_{\mathrm{m}}}{{E}_{\mathrm{k}}}+\frac{{\eta }_{\mathrm{m}}}{{E}_{\mathrm{m}}}\right)\dot{F}+\frac{{\eta }_{\mathrm{k}}{\eta }_{\mathrm{m}}}{{E}_{\mathrm{k}}{E}_{\mathrm{m}}}\ddot{F}=\pm {\eta }_{\mathrm{m}}\dot{u}\pm \frac{{\eta }_{\mathrm{k}}{\eta }_{\mathrm{m}}}{{E}_{\mathrm{k}}}\ddot{u} $$ (7) 由公式7可知,Burger’s蠕变模型的本构模型如公式8所示:
$$ \sigma +\left(\frac{{\eta }_{\mathrm{k}}}{{E}_{\mathrm{k}}}+\frac{{\eta }_{\mathrm{m}}}{{E}_{\mathrm{k}}}+\frac{{\eta }_{\mathrm{m}}}{{E}_{\mathrm{m}}}\right)\dot{\sigma }+\frac{{\eta }_{\mathrm{k}}{\eta }_{\mathrm{m}}}{{E}_{\mathrm{k}}{E}_{\mathrm{m}}}\ddot{\sigma }=+{\eta }_{\mathrm{m}}\dot{\varepsilon }+\frac{{\eta }_{\mathrm{k}}{\eta }_{\mathrm{m}}}{{E}_{\mathrm{k}}}\ddot{\varepsilon } $$ (8) 式中,$ \sigma $为应力;$ \epsilon $为应变。
1.2 基于离散单元法的Burger’s模型
连续介质理论中的Burger’s模型将试样视为连续体,模型作用的对象是整个连续的试样[18]。而在颗粒离散单元法中,当前时步颗粒间的接触力是根据上一时步的接触力值以及当前和上一时步颗粒间发生的位移来更新的,因此,接触力可采用有限差分法进行计算。石崇等分别对模型中开尔文体以及麦克斯韦体的接触力和位移进行了公式推导,同时根据黏弹性线性叠加原理给出总位移[17]。
1.2.1 模型中的开尔文体
对于模型中的开尔文体,其为线性弹簧和阻尼器并联构成,因此,位移与接触力之间的关系可由公式9表示[19]:
$$ {\dot{u}_{\mathrm{k}}}=\frac{-{E}_{\mathrm{k}}{u}_{\mathrm{k}}\pm F}{{\eta }_{\mathrm{k}}} $$ (9) 采用中心差分方式近似求取$ {u}_{\mathrm{k}} $和F的平均值,可得:
$$ \frac{{u}_{\mathrm{k}}^{t+1}-{u}_{\mathrm{k}}^{t}}{\Delta t}=\frac{1}{{\eta }_{\mathrm{k}}}\left[-\frac{{E}_{\mathrm{k}}\left({u}_{\mathrm{k}}^{t+1}+{u}_{\mathrm{k}}^{t}\right)}{2}\pm \frac{{F}^{t+1}+{F}^{t}}{2}\right] $$ (10) 式中,t表示当前时步的计算结果;t+1表示下一时步的计算结果;$ \Delta t $为时间步长。
对公式10进行整理,可得到公式11所示开尔文体位移的计算公式:
$$ {u}_{\mathrm{k}}^{t+1}=\frac{1}{A}\left[B{u}_{\mathrm{k}}^{t}\pm \frac{\Delta t}{2{\eta }_{\mathrm{k}}}\left({F}^{t+1}+{F}^{t}\right)\right] $$ (11) 式中,A、B为关于时间步长的函数,可分别由公式12、公式13求得:
$$ A=1+\frac{{E}_{\mathrm{k}}\Delta t}{2{\eta }_{\mathrm{k}}} $$ (12) $$ B=1-\frac{{E}_{\mathrm{k}}\Delta t}{2{\eta }_{\mathrm{k}}} $$ (13) 1.2.2 模型中的麦克斯韦体
对于模型中的麦克斯韦体,其为线性弹簧和阻尼器串联构成,因此,其位移及位移的一阶导数分别为公式14、15所示:
$$ {u}_{\mathrm{m}}={u}_{{\mathrm{m}}_{E}}+{u}_{{\mathrm{m}}_{\eta }} $$ (14) $$ {\dot{u}_{\mathrm{m}}}={\dot{u}_{{\mathrm{m}}_{E}}}+{\dot{u}_{{\mathrm{m}}_{\eta }}} $$ (15) 将公式6代入公式15可得:
$$ {\dot{u}_{\mathrm{m}}}=\pm \frac{\dot{F}}{{E}_{\mathrm{m}}}\pm \frac{F}{{\eta }_{\mathrm{m}}} $$ (16) 同样地,采用中心差分方式近似求取$ {u}_{\mathrm{m}} $和F的平均值,可得:
$$ \frac{{u}_{\mathrm{m}}^{t+1}-{u}_{\mathrm{m}}^{t}}{\Delta t}=\pm \frac{{F}^{t+1}-{F}^{t}}{{E}_{\mathrm{m}}\Delta t}\pm \frac{{F}^{t+1}+{F}^{t}}{2{\eta }_{\mathrm{m}}} $$ (17) 对公式17进行整理,可得到麦克斯韦体位移的计算公式18:
$$ {u}_{\mathrm{m}}^{t+1}=\pm \frac{{F}^{t+1}-{F}^{t}}{{E}_{\mathrm{m}}}\pm \frac{\Delta t\left({F}^{t+1}+{F}^{t}\right)}{2{\eta }_{\mathrm{m}}}+{u}_{\mathrm{m}}^{t} $$ (18) 1.2.3 整个Burger’s模型
对于整个Burger’s模型,根据黏弹性线性叠加原理,总位移u可由公式19给出:
$$ u={u}_{\mathrm{k}}+{u}_{\mathrm{m}} $$ (19) 对公式19进行中心差分近似,可得:
$$ {u}^{t+1}-{u}^{t}={u}_{\mathrm{k}}^{t+1}-{u}_{\mathrm{k}}^{t}+{u}_{\mathrm{m}}^{t+1}-{u}_{\mathrm{m}}^{t} $$ (20) 综上可得,最终接触力计算公式如公式21所示:
$$ {F}^{t+1}=\pm \frac{1}{C}\left[{u}^{t+1}-{u}^{t}+\left(1-\frac{B}{A}\right){u}_{\mathrm{k}}^{t}\mp D{F}^{t}\right] $$ (21) 式中,C、D为关于时间步长的函数,可分别由公式22、23求得:
$$ C=\frac{\Delta t}{2{\eta }_{\mathrm{k}}A}+\frac{1}{{E}_{\mathrm{m}}}+\frac{\Delta t}{2{\eta }_{\mathrm{m}}} $$ (22) $$ D=\frac{\Delta t}{2{\eta }_{\mathrm{k}}A}-\frac{1}{{E}_{\mathrm{m}}}+\frac{\Delta t}{2{\eta }_{\mathrm{m}}} $$ (23) 2. 微观参数标定
本研究采用PFC3D颗粒流数值模拟软件中的Burger’s接触模型作为海底稀软底质的本构模型进行模拟研究,为得较符合实际的数值模拟结果,数值模型的微观参数应与材料的宏观力学特性相互对应。石崇等[17]给了一种基于三轴试验的Burger’s接触模型参数的标定方法,是一种较为理想并应用广泛的进行颗粒微观参数标定的校准方法。本节通过在PFC3D中模拟常规三轴压缩实验,同时开展室内常规三轴试验,基于微观参数与宏观参数的关系,对微观参数进行调试,从而使数值模拟试验与室内三轴试验所表现出的宏观力学特性保持一致。
2.1 模型参数影响
在颗粒离散单元法中,Burger’s蠕变模型共有9个基本参数,如表1所示。
表 1 Burger’s蠕变模型参数Table 1. The parameters of Burger’s creep model作用方向 Maxwell体 Kelvin体 滑动摩擦系数
f弹性系数Em/MPa 黏性系数ηm/(MP·s) 弹性系数Ek/MPa 黏性系数ηk/(MP·s) 法向 bur_knm bur_cnm bur_knk bur_cnk bur_fric 切向 bur_ksm bur_csm bur_ksk bur_csk 为了明确每个参数对瞬时强度特性以及蠕变特性的影响,杨振伟等[18]在PFC3D中建立了常规三轴压缩实验的数值模型,在模拟试验中采用控制变量法对每个参数进行了分类研究。为方便计算,对各个参数的法向和切向设置了相同的参数值,从而将Burger’s接触模型参数简化为5个,即Em、ηm、Ek、ηk和f。
2.1.1 模型参数对瞬时强度特性的影响
取弹性模量、泊松比、单轴抗压强度等三项物理力学参数作为反应瞬时强度特性的指标,通过5组对比试验得到了不同模型参数取值下相应的瞬时强度结果。模型中各参数对瞬时强度特性的影响关系如表2所示。
表 2 Burger’s蠕变模型参数对瞬时强度特性的影响Table 2. The effect of Burger’s creep model parameters on instantaneous strength characteristics瞬时强度特性 Maxwell体 Kelvin体 滑动摩擦系数f 弹性系数Em 黏性系数ηm 弹性系数Ek 黏性系数ηk 弹性模量 正相关 不相关 不相关 不相关 正相关 泊松比 正相关 不相关 不相关 不相关 负相关 单轴抗压强度 正相关 不相关 不相关 不相关 正相关 由试验结果可知,在Burger’s蠕变接触模型中,仅有Maxwell体弹性系数和摩擦系数对瞬时强度特性存在显著影响,其余参数均对瞬时强度特性无显著影响。因此,在进行Burger’s蠕变模型参数标定调试时,只需调试Maxwell体弹性系数和滑动摩擦系数即可,对其他参数可直接匹配无需调整。
2.1.2 模型参数对蠕变特性的影响
同样地,通过对比试验记录在不同模型参数取值下轴向应变随时间的改变量,得出各模型参数不同取值情况的蠕变加载、卸载曲线。取瞬时应变量、起始应变量、起始应变率等八项参数作为反映蠕变特性的指标,模型中各参数对蠕变特性的影响关系如表3所示。
表 3 Burger’s蠕变模型参数对蠕变特性的影响Table 3. The effect of Burger’s creep model parameters on creep characteristics瞬时强度特性 Maxwell体 Kelvin体 滑动摩擦
系数f弹性系数
Em黏性系数
ηm弹性系数
Ek黏性系数
ηk瞬时应变量 负相关 不相关 不相关 不相关 负相关 起始应变量 不相关 不相关 不相关 不相关 不相关 起始应变率 不相关 负相关 不相关 负相关 不相关 稳定应变量 负相关 负相关 负相关 不相关 不相关 稳定应变率 不相关 不相关 不相关 不相关 不相关 瞬时恢复量 负相关 不相关 不相关 不相关 不相关 弹性后效回复率 不相关 不相关 负相关 负相关 不相关 残余应变量 不相关 负相关 不相关 不相关 不相关 试验表明,材料在经历短暂初始顺势变形后,开始进入等速蠕变阶段,此后应变增加速率随时间增长逐渐变小并趋于0;卸除偏应力,变形将急剧减小,但随着弹性应变的不断恢复,应变减小速率逐渐变小并最终趋于稳定。
不同模型参数下的蠕变特性表明,Maxwell体的弹性模量与瞬时应变量、瞬时恢复量,黏性系数与稳定蠕变率、残余应变量呈负相关关系;Kelvin体的弹性模量与起始蠕变量,黏性系数与起始蠕变率、弹性后效恢复率呈负相关关系;滑动摩擦系数则与瞬时应变量呈负相关关系。
2.2 物理标定试验
2.2.1 试验仪器及材料
试验仪器采用德国APS-WILLE Geotechnik伺服电机控制动三轴试验系统,该试验系统是由电机加载系统、三轴室及升降装置、自动体积压力控制器(VPC)、动态控制器以及计算机等部分组成(图2)。
动三轴试验仪加载活塞能补偿在动态和静态试验期间由于三轴室内部的加载杆的体积变化而导致的围压变化,可以准确控制轴向和围压之间的相位。电机作动器与操作软件相结合,能够生成并施加不同类型的加载波形到试样,最大荷载可达60 kN,荷载精度为0.1% F·S。三轴室采用可更换的上帽和底座,适用于直径为38 ~150 mm的试样。
试验材料为某区域海底底质原状土样,土样采用特定密封容器保存,同时在封口处用保鲜膜困扎密封。运输过程中按原状方向放置,加固容器以防止土样在运输过程中受到振动、冲击或其他影响,同时避免在运输过程中暴露土样于高温、高湿度等不适宜环境中。
采用Mastersizer 3000激光粒度分析仪对海底沉积物土样进行粒度分析,结果表明,沉积物中砂砾占2.3%~5.2%,黏粒占 38.1%~48.4%。采用多功能X射线衍射仪对沉积物进行矿物成分分析,沉积物土样主要由石英、白云母、方解石、钠长石和绿泥石等矿物组成。其中,石英占40.9%,白云母占29.3%,方解石占22.7%,钠长石占6.5%,绿泥石占0.6%。
采用GYS-2型液塑限测定仪、ZJ四联应变控制式直剪仪、WG型单杠杆固结仪等试验仪器,对海底沉积物土样进行天然含水率及天然密度试验、液塑限试验、四联剪切试验、固结压缩试验等室内试验。
根据室内土工试验结果,同时综合现场调查数据,得到所示该土样所处海域海底底质物理力学性质范围(表4)。
表 4 研究海域海底底质物理力学性质范围Table 4. Physical and mechanical properties of seabed sediments in the study area物理性质 结果范围 力学性质 结果范围 天然含水率/% 88.13~137.79 粘聚力
/kPa3.7~6.9 天然密度/(g/cm3) 1.3~1.5 内摩擦角
/(°)2.4~3.9 孔隙比 2.46~3.85 压缩系数/MPa−1 1.86~3.73 液性指数 0.96~1.97 压缩模量/MPa 1.26~2.13 塑性指数 37.2~57.8 贯入阻力
/N0.19~1.32 由试验结果可知,该海域海底底质土样具有较高的含水率、液限、可塑性和压缩性,具有更低的密度和强度特性。
2.2.2 试验内容及结果
常规三轴试验是较为理想的进行颗粒细观参数标定的基准试验,其作用只是为了对数值模拟中的微观参数进行标定,理论上保持数值模拟三轴试验与室内常规三轴试验二者的试验条件相对应,即可进行对微观参数的标定。因此本次研究选取固结排水试验,对3组南海原状土进行了不同围压的三轴压缩试验。试验土样为圆柱形试样,试样尺寸为直径×高度=38 mm×76 mm,剪切类型为固结排水剪切,3组围压分别为100、150、200 kPa。采用等应变速率剪切,考虑到剪切类型为固结排水剪切,根据《土工试验方法标准》(GB/T 50123—2019),剪切速率宜为0.003%/min~0.012%/min。试样轴向尺寸为76 mm,选取加载速率为0.008 mm/min [20]。具体试验方案如表5所示。
表 5 三轴试验方案Table 5. Triaxial test scheme试验次数 试样尺寸
(直径/mm×高度/mm)剪切类型 围压
/kPa加载速率
/
(mm/min)1 38×76 固结排水 100 0.008 2 38×76 固结排水 150 0.008 3 38×76 固结排水 200 0.008 根据《土工试验方法标准》(GB/T 50123—2019),具体试验步骤如下:
(a)该原状土为较软土样,采用切土盘削切成直径为38 mm的圆柱形土样,取试样高度为76 mm,按要求削平上下两端。放入饱和器在水下静止饱和10 h。
(b)在压力室底座依此放置透水石、滤纸、试样、滤纸、透水石以及试样帽,用橡胶膜套在试样外并扎紧在压力室底座及试样帽上。安装压力室罩并注水,对试样施加相应围压进行固结。
(c)固结完成后,按照试验方案选取的加载速率进行轴向加载,直至试样达到设定的轴向应变。
(d)在不同围压条件下,完成三组三轴压缩试验,提取试验数据并分析,得到不同围压条件下的应力应变曲线(图3)。
2.3 微观参数标定
2.3.1 数值模型建立
数值模拟采用PFC3D颗粒流软件,为减少颗粒数量,提高计算速度,三轴试验数值模拟中颗粒的尺寸相比于实际颗粒粒径有所放大。Huang等[21]指出,在三轴试验的离散元模拟中,当尺寸比SR(即试样最小边长与颗粒中值粒径的比值)相对较小时,会引起边界效应,而根据河海大学石崇等专家的研究结果可知,当尺寸比SR>40时,试样尺寸对数值计算结果的影响可忽略不计[22]。
因此,取尺寸比SR为40,根据室内试验测得60~150 μm的粒径级配,数值模型取同比例的粒径级配,生成粒径范围为0.027~0.068 m的颗粒。根据室内物理试验测定的研究海域海底底质物理力学性质范围表(表4),取孔隙率为0.7,生成尺寸为直径×高度=3.8 m×7.6 m的试样模型,等比放大模型尺寸可大大提高数值模型计算速度。在模型试样的侧面以及上下两面生成墙体,生成的常规三轴试验数值模型如图4所示。
侧面作为伺服控制平面,通过控制其收缩和扩张对模型试样施加整体围压并在伺服机制下保持围压的恒定。上下两墙体作为加载板,可以沿模型试样轴向方向上下移动,可通过定义上下加载板的移动速度来实现等应变速率的固结排水剪切试验。
2.3.2 数值模型参数标定
数值模型生成后,根据表4研究海域海底底质物理力学性质范围,定义模型的颗粒密度等已知参数。对模型中颗粒间的接触定义为Burger’s接触模型,为方便计算,同样对各个参数的法向和切向设置了相同的参数值,从而将Burger’s接触模型参数简化为5个。标定之前,同样根据表4中粘聚力、内摩擦角、压缩模量等参数,对5个基本参数初步赋予相应的经验值。
参数赋予后,通过伺服控制平面对模型试样施加相应恒定围压,实现试样的各向同性固结,固结过程中保持排水条件。固结完成后,对加载板施加等同于室内三轴试验的加载速度进行压缩试验,试验结束后输出应力-应变曲线。
根据首次试验得到的应力-应变曲线,并按照本文第2.1节中的参数调试方法对5个基本参数进行调试。经调试后,所采用的模型微观力学参数如表6所示。且不同围压下室内三轴试验与数值模型三轴试验所得应力-应变曲线对比结果如图5所示。
表 6 模型微观力学参数Table 6. Micromechanical parameters of the model参数类型 Maxwell体 Kelvin体 滑动摩擦
系数
f弹性系数
Em/MPa黏性系数
ηm/(MP·s)弹性系数
Ek/MPa黏性系数
ηk/MP·s参数取值 0.2 10 1.5 0.03 0.1 根据对比结果可以看出,数值模拟结果与室内试验结果较为吻合,因此,可以用该参数作为Burger’s接触模型的微观参数取值,同时根据该结果亦可论证选用Burger’s蠕变模型的合理性。
2.4 不同工况下模型参数取值
根据表2、3中蠕变模型参数分别对瞬时强度特性以及蠕变特性的影响机制,依据上述标定结果,改变不同模型参数获得不同底质条件的5种工况,从而研究不同底质条件下采矿车行驶的下陷情况。对应工况下的模型参数取值如表7所示。
表 7 不同工况下的模型微观力学参数Table 7. The micromechanical parameters of the model under different working conditions工况 Maxwell体 Kelvin体 滑动摩擦
系数f弹性系数
Em/MPa黏性系数
ηm/(MP·s)弹性系数
Ek/MPa黏性系数
ηk/(MP·s)1 0.2 10 1.5 0.03 0.1 2 0.1 10 1.5 0.03 0.1 3 0.3 10 1.5 0.03 0.1 4 0.2 5 1.5 0.03 0.1 5 0.2 15 1.5 0.03 0.1 3. 海底采矿车行驶数值模拟
3.1 计算模型
本文采用鲲龙500深海采矿车为模型原型,采矿车整体结构如图6所示,长5.4 m,宽3.5 m,高1.9 m,在陆地上总重9.5 t。
利用PFC3D软件中的Clump模块生成采矿车的简化模型,履带触地尺寸为0.8 m×4.6 m。通过赋予Clump模块体积、密度等参数,从而达到计算模型模拟采矿车整车质量以及通过履带对土体产生相应压力的效果,三维简化模型如图7所示。
按上述三轴试验数值模型的参数生成海底土层模型,由于计算能力的限制,在海底采矿车行驶模拟中,海底底质颗粒粒径同样有所放大。取尺寸比SR为40,按同比例的粒径级配修改数值模型中的粒径级配为0.07~0.18 m,最终生成尺寸为10 m×20 m×10 m的海底地层模型。定义颗粒间的接触为Burger’s接触模型,赋予表7中不同工况下的接触模型参数。
施加流体单元,打开渗流模块和浮力计算选项,读取流体单元的节点文件和单元构成文件,设置流体单元的密度为1040 kg/m3,流体动力粘滞系数为1×10−3 Pa·s,施加流体单元压力,定义静水压力。施加重力场,在重力和浮力作用下进行地层平衡,得到海底地层计算模型(图8)。
调用采矿车三维简化模型,在重力和浮力作用下下沉到地层表面,与地层接触,得到采矿车海底行驶计算模型(图9)。
履带式采矿车在海底移动的实际速度通常为0.5~1.0 m/s[14],因此本研究采用0.5、1.0、1.5以及2.0 m/s的速度进行仿真计算,模拟不同工况下采矿车在不同行驶速度中的下陷情况。同时,取速度为0,模拟静止情况下的下陷情况。
3.2 颗粒与流体耦合作用分析
PFC中颗粒与流体之间的运算是通过差分方法实现的[23-26],其指导思想是将颗粒与流体之间的相间作用力作为固相与液相之间的桥梁,通过该作用力在两相之间建立联系,在整个计算过程中,通过渗透介质变化反映流场的变化,进而不断更新相间作用力。
液相与固相之间的相互耦合过程如图10 所示,整个计算域内,颗粒采用离散介质方法从微观尺度进行研究,流体基于连续介质方法从中尺度层面考虑其平均值[27]。
图10中,首先按孔隙率生成颗粒,引入流场,并通过计算使得颗粒处于稳定状态。在当前时步,按照颗粒的位置决定颗粒之间重叠量,从而计算颗粒间的接触力,并将该接触力施加到颗粒的运动方程;另一方面,根据颗粒的分布确定流体计算单元内的孔隙率,根据初始条件和边界条件以及获得的单元孔隙率将流体方程进行离散求解,计算该时刻的流场,经迭代稳定后,获取各流体单元的压力梯度,用于计算固相与液相之间的相间作用力,将此作用力作为颗粒受到的外力施加到运动方程。颗粒在内外力作用下会产生新的加速度、速度、位移,颗粒位置重新调整,再由力-位移法则得出下一时步颗粒间的接触力,根据颗粒分布确定下一时步流体单元的孔隙率,将其代入流体方程,如此反复,从而实现颗粒与流体之间的耦合计算。
基本模拟过程如下:① 生成墙体和颗粒;② 设置墙体和颗粒属性;③ 读入流体网格;④ 设置流体属性;⑤ 循环求解。
3.3 计算结果与分析
图11所示不同工况下海底采矿车在不同车速下的下陷情况。由下陷曲线得出不同速度工况下采矿车稳定时间与稳定下陷量(表8)。
表 8 不同工况下采矿车各行驶速度的下陷数据Table 8. Subsidence data of a mining car at different driving speeds under different working conditions参数 工况 行驶速度 0 0.5/(m/s) 1.0/(m/s) 1.5/(m/s) 2.0/(m/s) 稳定
时间
/s1 9.5 7.4 6.8 5.8 3.5 2 9.6 7.8 6.5 6.1 3.8 3 9.5 8.0 6.3 5.6 3.2 4 9.7 8.9 7.9 7.0 6.3 5 8.2 6.3 5.6 4.9 2.6 稳定
下陷量/m1 0.331 0.183 0.162 0.132 0.114 2 0.372 0.226 0.184 0.165 0.130 3 0.281 0.163 0.137 0.116 0.091 4 0.335 0.186 0.163 0.130 0.115 5 0.331 0.179 0.159 0.129 0.109 对于土样所处海域,由工况1曲线可知下陷情况分为两个阶段:0~0.5 s时在初始装载期间发生的瞬时下陷,以及0.5 s后发生的蠕变下陷。且由其他工况曲线可知,不同底质条件下的下陷情况均可大致分为以上两个阶段。
进一步分析可知,对于同一种底质条件,海底采矿车在静止时下陷量最大,随着采矿车行驶速度的增加,其下陷量会随之减少并趋于稳定;蠕变随时间的变化而变化,当海底采矿车行驶速度更快时,海底稀软底质的蠕变时间减少,从而导致蠕变下陷量减小。对于不同底质条件,对比图11a、b、c可知,较大弹性系数的底质其下沉量较小,即具有较强的承载力;对比图11a、d、e可知,较大黏性系数的底质其下陷稳定时间较短,蠕变时间较少,即具有较低的蠕变性能。
本次研究中海底采矿车行驶的模拟下陷趋势与Bekker[28]在1969年的研究结论一致,因此,采用基于Burger’s蠕变模型的PFC颗粒流软件是可行的,本次研究具有一定的可靠性。
Xu等[29]于2016年通过制备人工材料来模拟太平洋多金属结核矿区的深海稀软沉积物的性质,进行了不同压应力下的压缩蠕变试验和不同压剪应力下的剪切蠕变试验,从而获得压剪蠕变参数,建立了压沉模型和压剪复合蠕变本构模型,利用RecurDyn软件计算了履带式矿车在不同速度时的压剪复合蠕变下陷。对比其结果,本次研究获得了一致的蠕变下陷趋势,进一步验证了本研究的可靠性。
4. 结论
(1)本文以某海域海底稀软底质原状土样为依托,开展了相关土工试验,综合现场调查资料得到了海底稀软底质的物理力学性质范围:物理性质中,天然含水率为88.13%~137.79%,天然密度为1.3~1.5 g/cm3,孔隙比为2.46~3.85,液性指数和塑性指数分别为0.96~1.97和37.2~57.8;力学性质中,粘聚力为3.7~6.9 kPa,内摩擦角为2.4°~3.9°,压缩系数为1.86~3.73 MPa−1,压缩模量为1.26~2.13 MPa,贯入阻力为0.19~1.32 N。试验结果表明,该海域海底稀软底质具有更高的黏粒含量、含水率和压缩性,具有更低的密度和强度特性,即承载力较低,蠕变性能较强。同时为进一步进行不同底质条件下海底采矿车在不同行驶速度时的下陷深度研究奠定了数据基础。
(2)模拟结果显示,对于同一种底质条件,下陷量会随时间呈非线性的变化,除装载期间发生的瞬时下陷外,约0.5 s后会因海底底质的蠕变性能而发生蠕变下陷,并逐渐趋于稳定;且稳定时间与行驶速度呈负相关,行驶速度越快,蠕变下陷量达到稳定时所用时间越少;就不同行驶速度而言,蠕变下陷量与行驶速度呈负相关,当行驶速度更快时,海底底质的蠕变时间减少,从而会导致蠕变下陷量的减小。
(3)对于不同底质条件,对较大弹性系数的底质其下沉量较小,即具有较强的承载力;较大黏性系数的底质其下陷稳定时间较短,蠕变时间较少,即具有较低的蠕变性能。
-
表 1 Burger’s蠕变模型参数
Table 1 The parameters of Burger’s creep model
作用方向 Maxwell体 Kelvin体 滑动摩擦系数
f弹性系数Em/MPa 黏性系数ηm/(MP·s) 弹性系数Ek/MPa 黏性系数ηk/(MP·s) 法向 bur_knm bur_cnm bur_knk bur_cnk bur_fric 切向 bur_ksm bur_csm bur_ksk bur_csk 表 2 Burger’s蠕变模型参数对瞬时强度特性的影响
Table 2 The effect of Burger’s creep model parameters on instantaneous strength characteristics
瞬时强度特性 Maxwell体 Kelvin体 滑动摩擦系数f 弹性系数Em 黏性系数ηm 弹性系数Ek 黏性系数ηk 弹性模量 正相关 不相关 不相关 不相关 正相关 泊松比 正相关 不相关 不相关 不相关 负相关 单轴抗压强度 正相关 不相关 不相关 不相关 正相关 表 3 Burger’s蠕变模型参数对蠕变特性的影响
Table 3 The effect of Burger’s creep model parameters on creep characteristics
瞬时强度特性 Maxwell体 Kelvin体 滑动摩擦
系数f弹性系数
Em黏性系数
ηm弹性系数
Ek黏性系数
ηk瞬时应变量 负相关 不相关 不相关 不相关 负相关 起始应变量 不相关 不相关 不相关 不相关 不相关 起始应变率 不相关 负相关 不相关 负相关 不相关 稳定应变量 负相关 负相关 负相关 不相关 不相关 稳定应变率 不相关 不相关 不相关 不相关 不相关 瞬时恢复量 负相关 不相关 不相关 不相关 不相关 弹性后效回复率 不相关 不相关 负相关 负相关 不相关 残余应变量 不相关 负相关 不相关 不相关 不相关 表 4 研究海域海底底质物理力学性质范围
Table 4 Physical and mechanical properties of seabed sediments in the study area
物理性质 结果范围 力学性质 结果范围 天然含水率/% 88.13~137.79 粘聚力
/kPa3.7~6.9 天然密度/(g/cm3) 1.3~1.5 内摩擦角
/(°)2.4~3.9 孔隙比 2.46~3.85 压缩系数/MPa−1 1.86~3.73 液性指数 0.96~1.97 压缩模量/MPa 1.26~2.13 塑性指数 37.2~57.8 贯入阻力
/N0.19~1.32 表 5 三轴试验方案
Table 5 Triaxial test scheme
试验次数 试样尺寸
(直径/mm×高度/mm)剪切类型 围压
/kPa加载速率
/
(mm/min)1 38×76 固结排水 100 0.008 2 38×76 固结排水 150 0.008 3 38×76 固结排水 200 0.008 表 6 模型微观力学参数
Table 6 Micromechanical parameters of the model
参数类型 Maxwell体 Kelvin体 滑动摩擦
系数
f弹性系数
Em/MPa黏性系数
ηm/(MP·s)弹性系数
Ek/MPa黏性系数
ηk/MP·s参数取值 0.2 10 1.5 0.03 0.1 表 7 不同工况下的模型微观力学参数
Table 7 The micromechanical parameters of the model under different working conditions
工况 Maxwell体 Kelvin体 滑动摩擦
系数f弹性系数
Em/MPa黏性系数
ηm/(MP·s)弹性系数
Ek/MPa黏性系数
ηk/(MP·s)1 0.2 10 1.5 0.03 0.1 2 0.1 10 1.5 0.03 0.1 3 0.3 10 1.5 0.03 0.1 4 0.2 5 1.5 0.03 0.1 5 0.2 15 1.5 0.03 0.1 表 8 不同工况下采矿车各行驶速度的下陷数据
Table 8 Subsidence data of a mining car at different driving speeds under different working conditions
参数 工况 行驶速度 0 0.5/(m/s) 1.0/(m/s) 1.5/(m/s) 2.0/(m/s) 稳定
时间
/s1 9.5 7.4 6.8 5.8 3.5 2 9.6 7.8 6.5 6.1 3.8 3 9.5 8.0 6.3 5.6 3.2 4 9.7 8.9 7.9 7.0 6.3 5 8.2 6.3 5.6 4.9 2.6 稳定
下陷量/m1 0.331 0.183 0.162 0.132 0.114 2 0.372 0.226 0.184 0.165 0.130 3 0.281 0.163 0.137 0.116 0.091 4 0.335 0.186 0.163 0.130 0.115 5 0.331 0.179 0.159 0.129 0.109 -
[1] 王贤觉, 陈毓蔚, 吴明清. 铁锰结核的稀土和微量元素地球化学及其成因[J]. 海洋与湖沼, 1984, 15(6):501-514 WANG Xianjue, CHEN Yuwei, WU Mingqing. Geochemistry of RE and trace elements in ferromanganese nodules and their genesis[J]. Oceanologia et Limnologia Sinica, 1984, 15(6):501-514.]
[2] 仲义, 陈忠, 莫爱彬, 等. 南海北部铁锰结核成因及元素的赋存状态[J]. 热带海洋学报, 2017, 36(2):48-59 ZHONG Yi, CHEN Zhong, MO Aibin, et al. Genetic types and elemental occurrence phases of ferromanganese nodules in the northern South China Sea[J]. Journal of Tropical Oceanography, 2017, 36(2):48-59.]
[3] 戴瑜, 刘少军. 深海采矿机器人研究: 现状与发展[J]. 机器人, 2013, 35(3):363-375 doi: 10.3724/SP.J.1218.2013.00363 DAI Yu, LIU Shaojun. Researches on deep ocean mining robots: status and development[J]. Robot, 2013, 35(3):363-375.] doi: 10.3724/SP.J.1218.2013.00363
[4] 陈小玲. 大洋多金属结核矿区表层沉积物的物理性质[J]. 东海海洋, 2004, 22(1):28-33 CHEN Xiaoling. The physical properties of surface sediments in oceanic polymetallic nodule[J]. Donghai Marine Science, 2004, 22(1):28-33.]
[5] Deepak C R, Shajahan M A, Atmanand M A, et al. Development tests on the underwater mining system using flexible riser concept[C]//Proceedings of the 4th Ocean Mining Symposium. Cupertino, USA: ISOPE, 2001: 94-98.
[6] 殷征欣, 王海峰, 韩金生, 等. 南海边缘海多金属结核与大洋多金属结核对比[J]. 吉林大学学报:地球科学版, 2019, 49(1):261-277 YIN Zhengxin, WANG Haifeng, HAN Jinsheng, et al. Comparison between the marginal-sea polymetallic nodules in South China Sea and Ocean Polymetallic nodules[J]. Journal of Jilin University:Earth Science Edition, 2019, 49(1):261-277.]
[7] 吴绍渊. 南海海底稀土元素研究进展[J]. 海洋科学, 2014, 38(3):116-121 doi: 10.11759/hykx20130918004 WU Shaoyuan. Research progress of rare earth elements of the seafloor in the South China Sea[J]. Marine Sciences, 2014, 38(3):116-121.] doi: 10.11759/hykx20130918004
[8] 马雯波, 饶秋华, 吴鸿云, 等. 深海稀软底质土宏观性能与显微结构分析[J]. 岩土力学, 2014, 35(6):1641-1646 MA Wenbo, RAO Qiuhua, WU Hongyun, et al. Macroscopic properties and microstructure analyses of deep-sea sediment[J]. Rock and Soil Mechanics, 2014, 35(6):1641-1646.]
[9] 饶秋华, 刘泽霖, 许锋, 等. 深海稀软底质特性及采矿车行走性能研究进展[J]. 中国有色金属学报, 2021, 31(10):2795-2816 doi: 10.11817/j.ysxb.1004.0609.2021-42120 RAO Qiuhua, LIU Zelin, XU Feng, et al. Research progress on characteristics of deep-sea soft sediment and walking performance of mining vehicle[J]. The Chinese Journal of Nonferrous Metals, 2021, 31(10):2795-2816.] doi: 10.11817/j.ysxb.1004.0609.2021-42120
[10] Leng D X, Shao S, Xie Y C, et al. A brief review of recent progress on deep sea mining vehicle[J]. Ocean Engineering, 2021, 228:108565. doi: 10.1016/j.oceaneng.2020.108565
[11] Schulte E, Handschuh R, Schwarz W. Transferability of soil mechanical parameters to traction potential calculation of a tracked vehicle[C]//Proceedings of the Fifth ISOPE Ocean Mining Symposium. Tsukuba, Japan: The International Society of Offshore and Polar Engineers, 2003: 123-131.
[12] Hong S, Kim H W, Choi J S. Transient dynamic analysis of tracked vehicles on extremely soft cohesive soil[C]//Proceedings of the 5th ISOPE Pacific/Asia Offshore Mechanics Symposium. Daejeon, Korea: The International Society of Offshore and Polar Engineers, 2002: 100-107.
[13] 戴瑜, 刘核, 张滔, 等. 海底履带式采矿车行走牵引通过性能研究[J]. 中国科技论文, 2015, 10(10):1203-1208 DAI Yu, LIU He, ZHANG Tao, et al. Research on traction trafficability of seafloor tracked mining vehicle[J]. China Sciencepaper, 2015, 10(10):1203-1208.]
[14] 陈峰, 桂卫华, 王随平, 等. 深海底履带机器车建模及仿真研究[J]. 机器人, 2004, 26(6):510-514 doi: 10.3321/j.issn:1002-0446.2004.06.007 CHEN Feng, GUI Weihua, WANG Suiping, et al. Modeling and simulation of a deep seabed tracked vehicle[J]. Robot, 2004, 26(6):510-514.] doi: 10.3321/j.issn:1002-0446.2004.06.007
[15] Qi C L, Rao Q H, Liu Q, et al. Traction rheological properties of simulative soil for deep-sea sediment[J]. Journal of Oceanology and Limnology, 2019, 37(1):62-71. doi: 10.1007/s00343-018-7258-7
[16] 房后国, 刘娉慧, 袁志刚. 海积软土固结过程中微观结构变化特征分析[J]. 水文地质工程地质, 2007, 34(2):49-52,56 doi: 10.3969/j.issn.1000-3665.2007.02.012 FANG Houguo, LIU Pinghui, YUAN Zhigang. Analysis on characteristics of microstructure change during marine soft soil consolidation[J]. Hydrogeology and Engineering Geology, 2007, 34(2):49-52,56.] doi: 10.3969/j.issn.1000-3665.2007.02.012
[17] 石崇, 张强, 王盛年. 颗粒流(PFC5.0)数值模拟技术及应用[M]. 北京: 中国建筑工业出版社, 2018 SHI Chong, ZHANG Qiang, WANG Shengnian. Numerical Simulation Technology and Application with Particle Flow Code (PFC5.0)[M]. Beijing: China Architecture & Building Press, 2018.]
[18] 杨振伟, 金爱兵, 周喻, 等. 伯格斯模型参数调试与岩石蠕变特性颗粒流分析[J]. 岩土力学, 2015, 36(1):240-248 doi: 10.16285/j.rsm.2015.01.033 YANG Zhenwei, JIN Aibing, ZHOU Yu, et al. Parametric analysis of Burgers model and creep properties of rock with particle flow code[J]. Rock and Soil Mechanics, 2015, 36(1):240-248.] doi: 10.16285/j.rsm.2015.01.033
[19] Itasca Consulting Group. PFC3D (Particle flow code in 3 dimensions) theory and background[R]. Minneapolis, USA: Itasca Consulting Group, 2008: 33-39.
[20] 中华人民共和国住房和城乡建设部. GB/T 50123-2019 土工试验方法标准[S]. 北京: 中国计划出版社, 2019 Ministry of Housing and Urban Rural Development of the People's Republic of China. GB/T 50123-2019 Standard for geotechnical test method[S]. Beijing: China Planning Press, 2019.]
[21] Huang X, Hanley K J, O'Sullivan C, et al. Effect of sample size on the response of DEM samples with a realistic grading[J]. Particuology, 2014, 15:107-115. doi: 10.1016/j.partic.2013.07.006
[22] 石崇, 徐卫亚. 颗粒流数值模拟技巧与实践[M]. 北京: 中国建筑工业出版社, 2015 SHI Chong, XU Weiya. Numerical Simulation Technique and Practice with Particle Flow Code[M]. Beijing: China Architecture & Building Press, 2015.]
[23] Shafipour R, Soroush A. Fluid coupled-DEM modelling of undrained behavior of granular media[J]. Computers and Geotechnics, 2008, 35(5):673-685. doi: 10.1016/j.compgeo.2007.12.003
[24] Shimizu H, Murata S, Ishida T. The distinct element analysis for hydraulic fracturing in hard rock considering fluid viscosity and particle size distribution[J]. International Journal of Rock Mechanics and Mining Sciences, 2011, 48(5):712-727. doi: 10.1016/j.ijrmms.2011.04.013
[25] Tomac I, Gutierrez M. Fluid lubrication effects on particle flow and transport in a channel[J]. International Journal of Multiphase Flow, 2014, 65:143-156. doi: 10.1016/j.ijmultiphaseflow.2014.04.007
[26] Catalano E, Chareyre B, Barthélémy E. Pore-scale modeling of fluid-particles interaction and emerging poromechanical effects[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2014, 38(1):51-71. doi: 10.1002/nag.2198
[27] 吴静. 岩溶隧道突水通道扩展机理、最小防突厚度及逃生路线优化研究[D]. 山东大学博士学位论文, 2018 WU Jing. Expansion of water inrush channel, minimum rock thickness and escape routes optimization of karst tunnel[D]. Doctor Dissertation of Shandong University, 2018.]
[28] Bekker M G. Introduction to Terrain-vehicle Systems[M]. Ann Arbor: University of Michigan Press, 1969.
[29] Xu F, Rao Q H, Ma W B. Predicting the sinkage of a moving tracked mining vehicle using a new rheological formulation for soft deep-sea sediment[J]. Journal of Oceanology and Limnology, 2018, 36(2):230-237. doi: 10.1007/s00343-018-6344-1
-
期刊类型引用(1)
1. 高懿鹏,金国庆,邹丽,于宗冰,程阳锐. 深海采矿车履齿结构沉陷特性数值研究. 中国造船. 2025(02): 146-157 . 百度学术
其他类型引用(0)