Influence of Sand Sampling Method on Bearing Capacity Calculation of Shallow Foundation in Discrete Element Method
-
摘要:
离散元数值模拟中，不同的制样方法会导致土体孔隙比和均匀性存在差异，进而对浅基础承载力的模拟计算结果产生影响，因此需要分析不同制样对浅基础承载力影响的问题. 本文分别使用粒径放大法、Distribute法、GM （grid method）法和欠层压实法对无黏性砂土进行制样，且试样在10
g 的重力场下进行地应力平衡；利用测量圆对不同位置土体孔隙比、水平应力和竖直应力进行监测，得到试样平均孔隙比e 和小于1的侧向土压力系数K 0值；通过在试样表面放置刚性墙体并以相同的速度加载来模拟浅基础承载力试验，研究不同制样方法对浅基础承载力的影响. 研究结果表明：GM法与欠层压实法生成的试样，其孔隙比接近最初设置的目标孔隙比，误差约为3.5%；而粒径放大法与Distribute法生成的试样，其孔隙比会小于目标孔隙比，误差为20.0%左右；在试样整体均匀性方面，GM法得到的试样均匀性最好，随后依次是欠层压实法、Distribute法和粒径放大法；由于不同制样方法所得的试样孔隙比和K 0不同，在浅基础承载力模拟计算中不同制样方法得到的承载力关系为：GM法 < 欠层压实法 < 粒径放大法 < Distribute法.Abstract:In the discrete element numerical simulation, different sample preparation methods will lead to differences in soil void ratio and uniformity, which will affect the simulation results of bearing capacity of shallow foundation. Therefore, it is necessary to analyze the influence of different sample preparation on bearing capacity of shallow foundation. Four methods (e.g., particle amplification method, distribute method, grid method, and under compaction method) were used to prepare the samples of cohesionless sand, and the samples were balanced under the gravity field of 10
g . The void ratio, horizontal stress and vertical stress of soil at different positions were monitored by measuring circle, and the average void ratioe and the lateral earth pressure coefficientK 0 value less than 1 were obtained. The influence of different sample preparation methods on the bearing capacity of shallow foundation was studied by placing rigid wall on the surface of the sample and loading at the same speed to simulate the bearing capacity test of shallow foundation. The results show that the porosity ratios of samples generated by the GM and under compaction method are closer to the original target porosity ratio with an error of about 3.5%. In comparison, the porosity ratios generated by particle amplification method and distribute method are smaller than the target porosity ratio, with an error of about 20.0%. Additionally, GM presents the most homogeneous sand samples, followed by the under compaction method, distribute method and particle amplification method, respectively. Due to the varying porosity ratios andK 0 of samples, the obtained bearing capacity of shallow foundation also changes. The relationship of bearing capacity obtained by different sample preparation methods in the simulation of bearing capacity of shallow foundation is : GM < underlayer compaction method < particle size amplification method < distribute method. -
兰新高铁地处中国西北地区，东起兰州，西至乌鲁木齐，全长1776 km，是世界上第一条一次性建成运营里程最长的高速铁路. 线路穿越我国新疆地区的百里、三十里等四大风区，风区内常年盛行七八级大风，曾发生过列车被大风吹翻的重大铁路安全事故[1]. 为保证列车安全顺利地通过风区，铁路建设部门在风区线路轨旁修建了高度为3.5～4.0 m的挡风墙. 但挡风墙“防车不防网”，气流经过挡风墙之后出现了明显的加速效应，使得正馈线发生剧烈舞动，导致线间放电、金具磨损加剧、线索疲劳断股，危及行车安全. 由于铁路运输对安全性要求很高，能够采用的防舞措施有限，低风压正馈线不改变接触网结构，有较好的适应性. 因此，设计一种新型低风压正馈线，用于降低正馈线的舞动幅值，保障牵引变电系统的安全可靠运行十分必要.
作用于架空输电导线的风载荷占整个线路所受风载荷的50%以上[2]，降低正馈线所受风压不仅可以减小线索舞动幅值，还可以降低因线间放电而引起的跳闸风险. 20世纪70年代，日本关西电气与住友电气等单位对多种不同表面形状的导线进行了风洞试验，研究结果表明，导线风阻力系数与导线表面形状有较大关系[3]. 20世纪90年代起，关西电气和住友电气对低风压导线的运行机理进行了持续研究，研究结果表明，低风压导线负压区面积小于普通导线，风阻力系数较小[4-5]. 我国在低风压导线领域起步较晚，但发展较快. 近年来，上海电缆研究所、无锡华能电缆有限公司、江苏中天科技股份有限公司等单位研发了不同表面结构的低风压导线产品，并申请了专利[6].
以6种不同表面形状的低风压正馈线和常规正馈线模型为研究对象，对各型正馈线在12～24 m/s （间隔3 m/s）风载荷作用下的舞动情况进行仿真，监测记录不同型号的低风压正馈线和常规正馈线在不同风载荷下的风阻力系数和跨距中点位移值，分析低风压正馈线表面结构参数对风阻力系数和跨距中点位移值变化的影响规律；对不同风载荷作用下风阻力系数较小的低风压正馈线建立三维有限元模型，并施加轴向拉力，分析其在轴向拉力作用下的形变及应力情况，研究成果为低风压正馈线的制造、选型及现场维护提供理论依据.
1. 模型建立
以现场架设的常规正馈线作为依据，建立常规正馈线二维模型，其二维截面结构如图1所示. 图中，R = 11.88 mm，为正馈线半径. 正馈线是由2层钢股和3层铝股相互绞合而成，其中，钢股直径为2.22 mm，铝股直径为2.85 mm.
图 1 常规正馈线二维截面示意Figure 1. Two-dimensional cross-sectional diagram of conventional catenary positive feeder本文以低风压导线运行机理为基础，设计凹槽数为8、不同凹槽半径的6种低风压正馈线模型，其中，一种低风压正馈线模型截面结构如图2所示，线径均为23.76 mm. 与常规正馈线相比，低风压正馈线模型的最外层由8根中间带有凹槽的铝股线构成，凹槽的小圆弧半径r 与常规正馈线半径R的比值为0.10～0.15，其余层股线结构参数与常规正馈线相同.
图 2 r/R = 0.14型低风压正馈线截面示意Figure 2. Cross-sectional diagram of low-wind-pressure catenary positive feeder with r/R = 0.142. 仿真计算与分析
2.1 仿真计算
对r/R = 0.10～0.15的6种低风压正馈线和常规正馈线模型进行12、15、18、21、24 m/s风载荷下的气动力特性仿真. 首先，设置计算域边界条件，将计算域左侧边界设置为速度入口，导线表面采用无滑移壁面边界，右侧边界设置为压力出口. 迭代收敛残差值取1 × 10−5，时间步长为0.005 s，计算1000步，以r/R = 0.14型低风压正馈线为例，仿真计算得到常规正馈线和低风压正馈线在18 m/s风载荷下的阻力系数CD时程图，如图3所示. r/R = 0.14型低风压正馈线在18 m/s风载荷下的阻力系数时程曲线幅值明显小于常规正馈线，为各低风压正馈线中阻力系数幅值最小的一种. 对45 m跨距正馈线在只受重力作用下线索悬垂状态找形成功后，将阻力系数以线性插值的方式添加，得到常规正馈线和r/R = 0.14型低风压正馈线跨距中点在18 m/s风载荷作用下的垂向位移时程图，如图4所示. 在同一坐标系下，r/R = 0.14型低风压正馈线垂向位移曲线幅值明显小于常规正馈线，即r/R = 0.14型低风压正馈线能够降低舞动幅值.
2.2 气动力参数分析
对6种低风压正馈线和常规正馈线在不同风载荷下的气动力特性进行仿真计算，得到不同线型在不同风载荷下的阻力系数，如图5所示. 统计6种低风压正馈线在18 m/s风载荷作用下跨距中点横向和垂向位移最大值，如图6所示.
图 5 不同低风压正馈线模型阻力系数曲线Figure 5. Drag coefficient curves for different types ofcatenary positive feeder models由图5可知，常规正馈线和6种低风压正馈线风阻力系数基本上随风速的增大呈减小趋势，且在整个测试风速范围内，低风压正馈线阻力系数均小于常规正馈线，说明本文设计的低风压正馈线模型具有防舞效果. 其中：常规正馈线气动力参数在整个测试风速范围内随风速的增大而减小， r/R = 0.10～ 0.11型低风压正馈线阻力系数交替出现下降和上升；r/R = 0.12～0.15型低风压正馈线阻力系数先上升之后一直下降；r/R = 0.14型低风压正馈线除在24 m/s风速时阻力系数略大于r/R = 0.15外，在其余测试风速时阻力系数均小于其他测试对象，为测试对象中降阻效果最好的线型.
由图6可知，常规正馈线跨距中点横向位移最大值为1.338 m，垂向位移最大值为1.240 m，低风压正馈线跨距中点位移最大值均较常规正馈线小，且最大位移值随r/R值的增大整体上呈现下降的趋势，但r/R = 0.15型低风压正馈线位移最大值较r/R = 0.14型低风压正馈线大，说明低风压正馈线设计的关键在于找到恰当的r/R值；r/R = 0.14型低风压正馈线位移幅值最小，其横向位移最大值为0.990 m，垂向位移最大值为0.910 m，较常规正馈线分别下降26.0%和26.6%，与图5风阻力系数相对应，说明了本文仿真方法的合理性以及仿真结果的准确性.
3. 低风压正馈线受力特性分析
正馈线在舞动过程中线夹出口处会产生大小不断变化的集中应力，会加剧连接金具的磨损以及线索疲劳断股. 为研究低风压正馈线在舞动时的受力特性，根据国际大电网会议的建议，以线夹出口89 mm处导线的动弯应变来评估低风压正馈线的运行状况[7]. 对图5、6各正馈线仿真结果分析发现，r/R = 0.12，0.13，0.14型3种低风压正馈线防舞效果较佳. 因此，建立3种低风压正馈线89 mm长度三维模型，在有限元软件中对自由端面施加一定拉伸载荷，模拟正馈线舞动时的受力情况，分析低风压正馈线发生舞动时的形变量及应力特性[8-9].
3.1 建立有限元模型
兰新高铁现场架设的常规正馈线结构参数如表1所示，该正馈线铝的杨氏模量为59 GPa，泊松比为0.30，钢的杨氏模量为190 GPa，泊松比为0.28，额定拉断力为83.42 kN. 本文以常规正馈线结构参数为依据，建立89 mm长度低风压正馈线三维模型，在固定端建立三维坐标轴O-XYZ. 其中，r/R = 0.14型低风压正馈线三维模型如图7所示，最外层铝股线中央带有凹槽，股线数为8股，其他层股线结构参数与常规正馈线相同.
表 1 常规正馈线结构参数Table 1. Structural parameters of conventional positive feeder材料 层数 股数/股 直径/mm 节径/mm 节距/mm 绞向 钢 最内层 1 2.22 次内层 6 2.22 21 139.86 左 铝 次外层 10 2.85 13 160.68 右 邻外层 16 2.85 12 216.72 左 最外层 22 2.85 11 261.36 右 
图 7 r/R = 0.14型低风压正馈线三维模型Figure 7. Three-dimensional model of positive feeder with r/R = 0.143.2 仿真方法验证
1） 仿真设置. 在有限元软件中设置钢铝股线的材料属性，将每根股线的中心节点等效为一个整体，在相邻股线之间建立接触对. 在设置边界条件时，将正馈线模型位于线夹一端完全固定约束，即在模型的固定端端面上约束X、Y、Z 3个方向的自由度. 为防止产生端部效应，在模型的自由端端面上建立一个刚域点，将该端面上所有自由度与该刚域点耦合形成一个刚域面，在刚域点上施加运行张力. 最后，采用扫掠法对模型进行正六面体网格划分[10-12]，得到的模型网格图如图8所示.
2） 方法验证. 我国一般将导线的运行张力设定为15%～25% 额定拉断力[13]（rated tensile strength，RTS），兰新高铁考虑正馈线新线系数后正馈线的最大许用张力为32137 N，为额定拉断力的38.5%，但正馈线长期在最大许用张力下运行，容易导致疲劳断股. 本文对常规正馈线施加25% RTS，即20.855 kN拉力，仿真得出各层股线轴向张力，将仿真结果与式（1）正馈线各层股线理论张力计算式结果进行对比[14]，结果见表2.
表 2 常规正馈线各层股线轴向张力Table 2. Axial tension of each layer of conventional positive feederkN 材料 层数 理论值 仿真值 误差 钢 最内层 5.963 5.921 0.242 次内层 7.839 7.818 0.221 铝 次外层 3.956 3.651 0.305 邻外层 5.045 4.889 0.156 最外层 3.383 3.112 0.271 Fn=πd2n4znEncosαn×cos(αn√(1+εi)2+(tanαn−μεi)2−1)， (1) 式中：Fn为第n层股线轴向张力；dn为股线直径；zn为股线数；En为股线弹性模量；αn为拉伸前股线捻角；εi为股线轴向伸长率；μ为该层股线的泊松比.
由表2可知，常规正馈线在25% RTS作用下，2层钢芯承担了54.1%的拉力，3层铝股线承担了45.9%的拉力，这与正馈线的设计初衷和现场运行情况相符. 仿真值与理论值存在一定误差，是因为理论计算时未考虑股线之间的挤压和摩擦对正馈线应力的影响. 但误差仍处于合理范围内，故本文所采用的仿真方法能够恰当地模拟正馈线在受到轴向拉力时的应力-应变特性.
3.3 结果分析
3.3.1 低风压正馈线受力形变及应力分析
在有限元软件中对r/R = 0.12，0.13，0.14型低风压正馈线三维模型自由端面施加25% RTS，分别比较3种正馈线的形变位移及应力变化情况[15-16]. 其中，r/R = 0.14型低风压正馈线在施加25% RTS后，轴向整体形变及各层股线轴向形变如图9所示.
图 9 r/R = 0.14型低风压正馈线形变云图Figure 9. Deformation nephogram of positive feeder with r/R = 0.14由图9（a）可知，正馈线在受到轴向拉力的作用下会发生一定轴向形变，固定端形变量最小，越往自由端，形变量越大，这与现实情况相吻合. 由图9（b）可知，各层股线在受到轴向拉力的情况下，形变量是不同的，轴向形变量从钢层到铝层逐渐增大，至最外层铝股线达到最大，即在现实中正馈线舞动时，最外层铝股线最容易发生金属疲劳断股，这是由于钢的杨氏模量大于铝的杨氏模量，铝的延展性较钢强[17-18]，在同一拉力作用下，铝的形变量大于钢. 其次，正馈线在受到轴向拉力时，内层股线受到外层股线的挤压，使得内层股线的形变量小于外层股线.
对r/R = 0.14型低风压正馈线施加25% RTS，等效应力云图见图10，其纵向中心截面应变及形变云图见图11.
图 10 r/R = 0.14型低风压正馈线应力Figure 10. Stress nephogram of positive feeder with r/R = 0.14
图 11 r/R = 0.14型低风压正馈线轴向截面云图Figure 11. Axial section nephogram of positive feeder with r/R = 0.14由图10（a）可知，次内层钢股线在正馈线振动时承担了大部分轴向应力，且除最内层直钢股应力分布均匀外，其他层股线应力分布不均，但具有一定的规律，即应力极值点的位置与正馈线每层股线的绞向相同，这是因为正馈线在轴向拉力作用下，同层相邻股线在绞合方向上相互挤压，形成应力集中点. 从图10（b）可知，钢股等效应力远大于铝股，即在正馈线舞动时钢股承受了很大一部分应力，且应力在直钢芯两侧对称分布，应力从固定端和自由端到中间的变化趋势相似，保证了正馈线受力对称.
由图11（a）可知，股线轴向应变沿最内层钢股两侧对称分布，中心钢股承受了最大拉变，这符合正馈线的设计初衷，邻外层钢股出现了最大压变，这是因为该层股线受到相邻两层股线的挤压. 从图11（b）可知，由固定端到自由端整体形变量逐渐增大，自由端形变量远大于固定端；由于中心钢股非螺旋结构，且钢的杨氏模量大，延展性较铝差，抗拉强度大，所以形变量最小；次内层钢股和邻外层铝股在绞线中间位置附近受到相邻股线的严重挤压，出现了最大形变量.
对r/R = 0.14型低风压正馈线施加25% RTS，沿轴向坐标轴Z轴截取固定端Z = 0至自由端端面，间隔8.9 mm，共11个轴向应力截面，分析其轴向应力，其中4个应力截面如图12所示.
图 12 25% RTS作用下r/R = 0.14型低风压正馈线横向截面应力云图Figure 12. Stress nephogram of horizontal cross-section of positive feeder with r/R = 0.14 under 25% RTS由图12可知，正馈线在轴向载荷作用下会发生一定程度的扭转. 正馈线的不同截面轴向应力分布不同，中心直钢芯上应力分布较均匀，螺旋钢层和铝层同一股线应力呈阶梯状分布，且应力极值点位置与正馈线层绞合方向相关；从固定端端面到自由端端面的不同截面上，应力极值点位置从与绞合方向相反到与绞合方向相同，但均位于股线接触处，说明在股线接触位置容易出现应力集中点，易引起股线磨损断股；中间位置截面轴向应力呈现中心对称，且在提取的截面中应力最小，轴向应力整体上呈现两端大中间小的情况.
3.3.2 3种低风压正馈线受力特性比较
为更准确地反映3种低风压正馈线在25% RTS作用下的受力特性，对状态变化参数进行比较，得到3种低风压正馈线25% RTS作用下状态参数变化曲线，如图13所示.
图 13 3种低风压正馈线25% RTS作用下状态参数变化曲线Figure 13. Variation of state parameter for the three types of positive feeders under 25% RTS由图13（a）可知，2层钢股承受了绝大部分应力，邻外层铝股在3层铝股中承受的应力最大，不同型号低风压正馈线在25% RTS作用下，从内向外第1、4层股线承受的应力基本相同，最外层股线最大应力随r/R值的增大而增大. 由图13（b）可知，在轴向拉力的作用下，股线位移量从内层到外层依次增大，表明最外层铝线最容易疲劳断股. 由于中心钢股没有螺旋，抗拉强度大，所以r/R = 0.12，0.13，0.14型低风压正馈线中心钢股位移量基本相等，r/R = 0.13，0.14型低风压正馈线第2、4层股线位移量基本相等；r/R = 0.12型低风压正馈线第2、5层股线位移量均小于其他低风压正馈线，第5层股线位移量随r/R值的增大而增大.
4. 结 论
通过对设计的低风压正馈线进行防舞有效性仿真分析，发现低风压正馈线达到了预期效果，并对其中3种防舞效果较佳的低风压正馈线进行受力特性分析，得出以下结论，并对低风压正馈线的制造、选型及现场维护给出以下建议：
1） 低风压正馈线在受到轴向拉力时，铝股线的轴向形变量大于钢股线，为了平衡正馈线抗拉强度和导电性能，在制造正馈线时可以考虑将钢股线和铝股线交替绞制.
2） 在轴向载荷作用下，低风压正馈线在线夹处受到很大的应力，加剧了正馈线与线夹的磨损，容易导致股线断裂. 在正馈线的日常维护中应该加强巡视，确保线夹转动灵活，减小线夹处正馈线的静应力.
3） 低风压正馈线在轴向载荷作用下会发生一定程度的扭转，在股线接触处出现应力集中点，应力集中点位置与正馈线绞合方向相关. 在制造正馈线时可以考虑在股线表面覆缓冲层，减缓股线之间的振荡冲击，延长正馈线使用寿命.
4） 对3种防舞效果较为理想的低风压正馈线进行受力分析，发现在同一轴向载荷作用下，最外层铝股线层的位移量与r/R比值成正比例关系. 而r/R值越大，最外层铝股线横截面积越小，在正馈线舞动时越容易疲劳断股，因此，在低风压正馈线选型时应该综合考虑，平衡防舞有效性与使用寿命.
致谢：兰州交通大学天佑创新团队计划（TY202010）资助.
-
图 4 位置1～3土体平均水平应力与竖向应力分布
Figure 4. Distribution of average horizontal and vertical stresses of soil in location 1−3
表 1 模型细观参数
Table 1. Microscopic parameters of the model
参数 值 粒径/mm 0.8～1.2 颗粒密度/（kg•m−3） 2650 初始孔隙比 0.2 颗粒法向刚度 kn/ （N•m−1） 5 × 107 颗粒切向刚度 ks/ （N•m−1） 2.5 × 107 墙体法向刚 knw/ （N·m−1） 1 × 109 墙体切向刚 ksw/ （N•m−1） 1 × 109 初始摩擦系数 μ1 8.7 × 10−4 最终摩擦系数 μ2 0.839 墙体与颗粒摩擦系数 μ3 0 
下载: 导出CSV
表 2 试样的极限承载力
Table 2. Ultimate bearing capacity of samples
制样方法 试样 e 均匀性 极限承载力/kPa GM 法 G-S 0.194 1 12.00 欠层压实法 Q-S 0.193 2 12.50 Distribute 法 D-S 0.157 3 18.19 粒径放大法 L-S 0.164 4 14.30
下载: 导出CSV
-
[1] CUNDALL P A, STRACK O D L. A discrete numerical model for granular assemblies[J]. Géotechnique, 1979, 29(1): 47-65. [2] JIANG M J, JIANG T, CROSTA G B, et al. Modeling failure of jointed rock slope with two main joint sets using a novel DEM bond contact model[J]. Engineering Geology, 2015, 193: 79-96. doi: 10.1016/j.enggeo.2015.04.013 [3] BYM T, MARKETOS G, BURLAND J B, et al. Use of a two-dimensional discrete-element line-sink model to gain insight into tunnelling-induced deformations[J]. Géotechnique, 2013, 63(9): 791-795. [4] ZHENG G, YANG X Y, ZHOU H Z, et al. A simplified prediction method for evaluating tunnel displacement induced by laterally adjacent excavations[J]. Computers and Geotechnics, 2018, 95: 119-128. doi: 10.1016/j.compgeo.2017.10.006 [5] LIU Z N, KOYI H A. Kinematics and internal deformation of granular slopes: insights from discrete element modeling[J]. Landslides, 2013, 10(2): 139-160. doi: 10.1007/s10346-012-0318-8 [6] XU W J, DONG X Y. Simulation and verification of landslide tsunamis using a 3D SPH-DEM coupling method[J]. Computers and Geotechnics, 2021, 129: 103803.1-103803.16. [7] CHANG C S, CHAO S J. Discrete element method for bearing capacity analysis[J]. Computers and Geotechnics, 1991, 12(4): 273-288. doi: 10.1016/0266-352X(91)90026-C [8] Itasca. Particle flow code in two dimensions[M]. Minnesota: Itasca Consulting Group, Inc. 2004. [9] 戴轩,郑刚,程雪松,等. 基于DEM-CFD方法的基坑工程漏水漏砂引发地层运移规律的数值模拟[J]. 岩石力学与工程学报,2019,38(2): 396-408.DAI Xuan, ZHENG Gang, CHENG Xuesong, et al. Numerical simulation of ground movement induced by leakage of groundwater and sand in excavations based on the DEM-CFD method[J]. Chinese Journal of Rock Mechanics and Engineering, 2019, 38(2): 396-408. [10] 徐琨,周伟,马刚,等. 基于离散元法的颗粒破碎模拟研究进展[J]. 岩土工程学报,2018,40(5): 880-889. doi: 10.11779/CJGE201805013XU Kun, ZHOU Wei, MA Gang, et al. Review of particle breakage simulation based on DEM[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(5): 880-889. doi: 10.11779/CJGE201805013 [11] 金磊,曾亚武. 土石混合体宏细观力学特性和变形破坏机制的三维离散元精细模拟[J]. 岩石力学与工程学报,2018,37(6): 1540-1550.JIN Lei, ZENG Yawu. Refined simulation for macro-and meso-mechanical properties and failure mechanism of soil-rock mixture by 3D DEM[J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(6): 1540-1550. [12] 景路,郭颂怡,赵涛. 基于流体动力学-离散单元耦合算法的海底滑坡动力学分析[J]. 岩土力学,2019,40(1): 388-394.JING Lu, GUO Songyi, ZHAO Tao. Understanding dynamics of submarine landslide with coupled CFD-DEM[J]. Rock and Soil Mechanics, 2019, 40(1): 388-394. [13] 徐文杰. 滑坡涌浪流-固耦合分析方法与应用[J]. 岩石力学与工程学报,2020,39(7): 1420-1433.XU Wenjie. Fluid-solid coupling method of landslide tsunamis and its application[J]. Chinese Journal of Rock Mechanics and Engineering, 2020, 39(7): 1420-1433. [14] 马华庆,赵永志. 喷动流化床中杆状颗粒混合特性的CFD-DEM模拟[J]. 浙江大学学报(工学版),2020,54(7): 1347-1354. doi: 10.3785/j.issn.1008-973X.2020.07.013MA Huaqing, ZHAO Yongzhi. CFD-DEM investigation on mixing of rod-like particles in spout-fluid bed[J]. Journal of Zhejiang University (Engineering Science), 2020, 54(7): 1347-1354. doi: 10.3785/j.issn.1008-973X.2020.07.013 [15] CHO S E, LEE S R. Instability of unsaturated soil slopes due to infiltration[J]. Computers and Geotechnics, 2001, 28(3): 185-208. doi: 10.1016/S0266-352X(00)00027-6 [16] CHEN X Y, ZHANG L L, ZHANG L M, et al. Modelling rainfall-induced landslides from initiation of instability to post-failure[J]. Computers and Geotechnics, 2021, 129: 103877.1-103877.16. [17] YANG S C, CI W, LESHCHINSKY B, et al. Scale effects on the ultimate bearing capacity of rectangular footings placed on slopes[J]. Computers and Geotechnics, 2021, 137: 104254.1-104254.7. [18] YANG S C, LESHCHINSKY B, CUI K, et al. Influence of failure mechanism on seismic bearing capacity factors for shallow foundations near slopes[J]. Géotechnique, 2021, 71(7): 594-607. [19] YANG S C, LESHCHINSKY B, CUI K, et al. Unified approach toward evaluating bearing capacity of shallow foundations near slopes[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2019, 145(12): 04019110.1-04019110.16. [20] NG T T, MEYERS R. Side resistance of drilled shafts in granular soils investigated by DEM[J]. Computers and Geotechnics, 2015, 68: 161-168. doi: 10.1016/j.compgeo.2015.04.009 [21] WANG Z J, JACOBS F, ZIEGLER M, et al. Visualisation and quantification of geogrid reinforcing effects under strip footing loads using discrete element method[J]. Geotextiles and Geomembranes, 2020, 48(1): 62-70. doi: 10.1016/j.geotexmem.2019.103505 [22] FU Z Z, CHEN S S, LIU S H. Discrete element simulations of shallow plate-load tests[J]. International Journal of Geomechanics, 2016, 16(3): 1-12. [23] XU R, LIU E L, JIANG X Q, et al. Analysis on evolution of mesostructure of cohesionless soil ground upon loading[J]. Powder Technology, 2020, 368: 1-17. doi: 10.1016/j.powtec.2020.04.048 [24] CHA M S, SANTAMARINA J C. Effect of dissolution on the load-settlement behavior of shallow foundations[J]. Canadian Geotechnical Journal, 2016, 53(8): 1353-1357. doi: 10.1139/cgj-2014-0370 [25] 蒋明镜,李磊,孙超. 地基承载力尺寸效应离散元分析[J]. 同济大学学报(自然科学版),2015,43(5): 649-656.JIANG Mingjing, LI Lei, SUN Chao. Analysis on scale effect of bearing capacity of foundation by discrete element method[J]. Journal of Tongji University (Natural Science), 2015, 43(5): 649-656. [26] ZHONG W H, LIU H L, WANG Q, et al. Investigation of the penetration characteristics of snake skin-inspired pile using DEM[J]. Acta Geotechnica, 2021, 16(6): 1849-1865. doi: 10.1007/s11440-020-01132-2 [27] DUAN N, CHENG Y P. A modified method of generating specimens for a 2D DEM centrifuge model[C]//Geo-Chicago 2016: Sustainability, Energy, and the Geoenvironment. Chicago: American Society of Civil Engineers, 2016: 610-620. [28] JIANG M J, KONRAD J M, LEROUEIL S. An efficient technique for generating homogeneous specimens for DEM studies[J]. Computers and Geotechnics, 2003, 30(7): 579-597. doi: 10.1016/S0266-352X(03)00064-8 [29] GOURVENEC S M, MANA D S K. Undrained vertical bearing capacity factors for shallow foundations[J]. Géotechnique Letters, 2011, 1(4): 101-108. [30] KUSAKABE O. Geotechnical centrifuge technology[M]. London: Blackie Academic & Professional. 1995: 118-167. [31] LAU C K, BOLTON M D. The bearing capacity of footings on granular soils Ⅱ: experimental evidence[J]. Géotechnique, 2011, 61(8): 639-650. [32] CERATO A B, LUTENEGGER A J. Bearing capacity of square and circular footings on a finite layer of granular soil underlain by a rigid base[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2006, 132(11): 1496-1501. doi: 10.1061/(ASCE)1090-0241(2006)132:11(1496) [33] DAVIS R O, SELVADURAI A P S. Plasticity and Geomechanics[M]. Cambridge: Cambridge University Press, 2002. [34] 申志福,蒋明镜,朱方园,等. 离散元微观参数对砂土宏观参数的影响[J]. 西北地震学报,2011,33(增1): 160-165.SHEN Zhifu, JIANG Mingjing, ZHU Fangyuan, et al. Influence of the micro parameters of discrete element on the macro parameters of sands[J]. Northwestern Seismological Journal, 2011, 33(S1): 160-165. [35] TRAN V D H, MEGUID M A, CHOUINARD L E. A finite-discrete element framework for the 3D modeling of geogrid-soil interaction under pullout loading conditions[J]. Geotextiles and Geomembranes, 2013, 37: 1-9. doi: 10.1016/j.geotexmem.2013.01.003 [36] KIM B S, PARK S W, KATO S. DEM simulation of collapse behaviours of unsaturated granular materials under general stress states[J]. Computers and Geotechnics, 2012, 42: 52-61. doi: 10.1016/j.compgeo.2011.12.010 [37] NGO N T, INDRARATNA B, RUJIKIATKAMJORN C. DEM simulation of the behaviour of geogrid stabilised ballast fouled with coal[J]. Computers and Geotechnics, 2014, 55: 224-231. doi: 10.1016/j.compgeo.2013.09.008 [38] BHANDARI A, HAN J. DEM study of a shallow foundation under vertical loading[C]//2009 International Foundation Congress and Equipment. Florida: American Society of Civil Engineers, 2009: 465-472. [39] LI L, JIANG M J, LI T, et al. DEM analyses of the whole failure process of shallow foundation in plate load test on dense sand[C]//IOP Conference Series: Earth and Environmental Science. Warwick: Earth and Environmental Science, 2015: 012039.1-012039.6. 期刊类型引用(1)
1. 詹宝容，庾锡昌. 联合收割机裂纹转子与滚动轴承故障诊断系统研究——基于卷积神经网络. 农机化研究. 2024(05): 187-191 . 
百度学术其他类型引用(4)
-
下载:
下载:
百度学术
点击查看大图
计量
- 文章访问数: 365
- HTML全文浏览量: 183
- PDF下载量: 43
- 被引次数: 5