土工格栅与尾矿界面摩擦特性的拉拔数值试验研究.pdf

第7 2 卷第4 期有 色 金属 矿山部分 2 0 2 0 年7 月 D o i l O ．3 9 6 9 ／j ．i s s n ．1 6 7 l ‘4 1 7 2 ．2 0 2 0 ．0 4 ．0 2 3 土工格栅与尾矿界面摩擦特性的拉拔数值 试验研究 高太1 ，王汝杰1 ，高锦财1 ，杜常博2 1 ．鞍山钢铁集团有限公司齐大山选矿厂，辽宁 鞍山1 1 4 0 4 3 ； 2 ．辽宁工程技术大学土木工程学院，辽宁阜新1 2 3 0 0 0 摘要研究土工格栅加筋尾矿坝的首要问题是格栅一尾矿界面摩擦特性，采用了F L A C “’数值模拟分析了土 工格栅在尾矿填料中的拉拔过程，探究含水率与上覆压力对于格栅一尾矿界面摩擦特性的影响规律，并且将其与 室内拉拔试验的数据进行比较分析。研究结果表明土工格栅与尾矿数值的计算结果与试验得到的数据基本符 合；认为在土工格栅与尾矿数值计算过程当中，F L A c ”数值模拟可以很好响应实验室测得的界面摩擦参数 似黏 聚力和似摩擦角 ；在数值模拟过程中，界面摩擦参数的选择需要进行适度的调节，其中似黏聚力数值可取1 ．1 倍 试验值，似摩擦角数值可取o ．7 倍试验值。研究结果可为土工格栅加筋尾矿坝设计提供一定参考。 关键词F L A C3 0 数值模拟；土工格栅；尾矿；界面摩擦特性 中图分类号T D 9 2 6 ．4 1文献标志码A文章编号1 6 7 1 4 1 7 2 2 0 2 0 0 4 一0 1 2 3 一0 6 P u l l o u tn u m e r i c a le x p e r i m e n t a t i o ns t u d yo ni n t e r f a c ef r i c t i O n c h a r a c t e r i s t i c sb e t w e e ng e O g r i da n dt a i l i n g s G A OT a i l ，W A N GR u j i e l ，G A OJ i n c a i l ，D UC h a n g b 0 2 1 ．Q i d a s h a nC o n c e n t r a t o ro fA n s h a nI r o na n dS t e e lG r o u pC o ．，L t d ．，A n s h a nL i a o n i n g1 1 4 0 4 3 ，C h i n a ； 2 ．C o l l e g eo fC i v 订E n g i n e e r i n g ，L i a o n i n gT e c h n i c a lU n i v e r s i t y ，F u x i nL i a o n i n g12 3 0 0 0 ，C h i n a A b s t r a c t I n t e r f a c ef r i c t i o nc h a r a c t e r i s t i c so fg e o g r i d t a i l i n g sa r et h ek e yq u e s t i o ni nt h es t u d yo fg e o g r i d r e i n f o r c e dt a i l i n g sd a m ．F L A C 3 0n u m e r i c a ls i m u l a t i o nw a su s e dt oa n a l y z et h ep u l 卜o u tp r o c e s so fg e o g r i di n t a i l i n g s ．T h ei n f l u e n c eo fo v e r l y i n gp r e s s u r ea n dw a t e rc o n t e n to nt h ei n t e r f a c ef r i c t i o nc h a r a c t e r i s t i c so fg e o g r i d t a 订i n g sw a sd i s c u s s e da n dc o m p a r e dw i t ht h ei n d o o rp u l l 一o u tt e s t ．T h er e s u l t si n d i c a t e dt h a tt h en u m e r i c a lr e s u l t s b e t w e e ng e o g r i d sa n dt a i l i n g sw e r eb a s i c a l l yc o n s i s t e n tw i t ht h ee x p e r i m e n t a ld a t a ．1 ti sc o n s i d e r e dt h a tF L A C 3 D n u m e r i c a Is i m u I a t i o nc a nw e l lr e s p o n dt ot h ei n t e r f a c ef r i c t i o np a r a m e t e r s p s e u d oc o h e s i o na n dp s e u d of r i c t i o na n 9 1 e m e a s u r e di n1 a b o r a t o r y ．I nt h ep r o c e s so fn u m e r i c a ls i m u l a t i o n ，t h ei n t e r f a c ef r i c t i o np a r a m e t e r sn e e dt ob ea d ju s t e d a p p r o p r i a t e l y ，i nw h i c ht h ep s e u d oc o h e s i o nv a l u e sc a nb ec h o s e na s1 ．1t i m e so ft h ee x p e r i m e n t a lv a l u e ，a n dt h e p s e u d of r i c t i o na n 9 1 ev a l u e sc a nb ec h o s e na sO ．7t i m e so ft h ee x p e r i m e n t a lv a l u e s ．T h er e s e a r c hr e s u l t sc a np r o v i d e s u p p o r tf o rt h ed e s i g no fg e o g r i dr e i n f o r c e dt a i l i n g sd a m ． K e yw o r d s F L A C 3 Dn u m e “c a ls i m u l a t i o n ；g e o g “d ；t a i l i n g s ；i n t e r f a c ef r i c t i o nc h a r a c t e r i s t i c s 在土工格栅加筋填料的受力过程当中，由于填 料和格栅两种不同材料的差异性，使他们存有相对 滑动趋势，但是格栅和填料之间存在的摩阻力会阻 拦这种趋势，这使格栅受到相对来说较大的拉力，这 基金项目鞍钢科研项目 2 0 1 8 ～科A 1 9 作者简介高太 1 9 6 0 一 ，男，高级工程师，矿山工程专业。 通信作者杜常博 1 9 7 1 一 。男，博士研究生，矿物加工工程专业。 基本相当于提升了加筋土体的抗拉强度，同时，土工 格栅会对加筋土体造成约束，这样也使加筋土体的 抗剪强度有了相比之下的提高，所以筋土这种复合 体的加固机制和破环形式是有探究的必要的，近几 年针对筋土界面的摩擦特性展开的国内外研究主要 在试验和数值模拟两个方向。 数值模拟研究是对试验研究的一个有效补充和 万方数据 1 2 4 有 色金 属 矿山部分第7 2 卷 验证，可以得到传统方法无法得到的数据和规律，许 多学者使用有限元方法、离散元方法及有限差分方 法研究筋土界面的摩擦特性[ 卜引万亮Ⅲ提出了一个 A N S Y S 计算模型用来模拟筋土拉拔的实验过程， 模拟结果与原本实验结果基本相同，更好地验证了 模型的合理性；李首龙[ 5 ] 利用A B A Q U S 建立了土 工格栅拉拔有限元模型，改变网格尺寸、拉拔速率等 相关的参数，并且将计算值与试验值进行比较分析； 高文华等[ 6 ] 基于F L A C 软件进行了加筋边坡稳定 性计算，计算结果与试验结果吻合良好；杜勇立[ 7 ] 采 用F L A C 3 D 分析了筋土界面特性，与试验结果较为 符合；C H E N 等[ 8 ] 采用离散元方法讨论了周期荷载 下约束区范围、深度、栅孔对颗粒材料的约束作用； 赖汉江等‘9 3 和周健等‘1 0 3 均采用离散元软件P F C 2 D 模拟拉拔试验，从细观角度分析筋土界面特性；苗 晨曦等[ 3 和郑俊杰等[ 1 2 3 均采用“c l u m p ”方法建立 三维离散元模型，较真实地模拟了筋土界面细观 结构。 本文采用F L A C 3 D 数值模拟，对在尾矿填料中 的土工格栅的拉拔过程进行分析，探究含水率、上覆 压力对土工格栅和尾矿界面摩擦特性的影响规律， 并且将其与在室内开展的拉拔试验进行对比并分 析，互相监测数值计算结果和室内试验结果，以便更 加准确和全面地对格栅一尾矿的界面摩擦特性进行 分析。 1 拉拔数值计算分析 1 ．1 数值计算假定 进行数值模拟计算时假定尾矿填料各向同性、 连续、均质；土工格栅是线弹性材料；试验箱各面刚 度无限大 除加载面 。 1 ．2 拉拔试验模拟方案 由于尾矿采用干堆，其天然含水率在8 ％～1 0 ％， 由于将含水率的变化纳人考虑，试验模拟采用 1 ．5 ％、9 ．7 ％、1 9 ．3 ％三种含水率进行模拟试验，具 体拉拔模拟试验方案设计见表1 。 表l 试验方案 T 曲I el1 b tp h n 1 ．3 数值计算模型 拉拔试验的数值计算模型如图1 所示，土工格栅 置于尾矿填料中间，试验箱的内尺寸为5 0 0m m 4 0 0m m 4 0 0m m ，筋材尺寸为5 0 0m m 1 8 0m m ， 试验模拟中尾矿填料的本构模型采用M o h r C o u l o m b 模型。模型的底面固定，4 个侧面限制其 X 、y 方向的位移，模型顶面为应力边界，在其上面 施加法向荷载。 图l 数值计算模型 F i g ．1 N u m e r i c a Ic a I c u I a t j 蛐n - O d e I 数值计算选用的材料参数如下 1 尾矿采用实体单元，M o h r - C o u l o m b 模型输 入的参数有弹性体积模量 b u l k 、弹性切变模量 s h e a r 、黏聚力 c o h e s i o n 、内摩擦角 f r i c t i o n 、密度 d e n s i t y 等。不同含水率下尾矿的建模参数见表2 。 2 土工格栅选用G e o g r i d 单元，输入的参数 有弹性模量 E 、泊松比 v 、厚度 t h i c k 、耦合 弹簧单位面积上的刚度 c S _ s k 、耦合弹簧黏聚力 c S s c o h 、耦合弹簧摩擦角 c S j r i c 。不同含水率 下格栅一尾矿的界面摩擦参数见表3 。 表2 不同含水率的尾矿建模参数 7 r a b I e2P a r a m e t e 体o ft a i l i n 笋u n d e rd i f f b 他n tw a t e r ∞n t e n t s 万方数据 第4 期高太等土工格栅与尾矿界面摩擦特性的拉拔数值试验研究1 2 5 衰3 不同含水率下筋土一尾矿界面摩擦参数 T a b l e3I n t e r f a ∞f r j c t i 蚰a Ip a r 绷e t e 隋o ft a i l i n 擎 u n d e rd i f f e r e n tw a t e rc o n t e n t s 数值计算结果及分析 经过数值计算得出不同上覆压力的作用下土工 格栅拉拔位移与拉拔力之间的关系 以含水率、 1 ．5 ％上覆压力1 ．0 2k P a 为例 ，如图2 所示，具体 数值见表4 。通过对数值的分析可以得知，土工格 栅拉拔力随拉拔位移的变化规律都是前面先表现出 线性增长的趋势，然后增长幅度逐渐变小，直至减小 至o ，拉拔力达到最大值，稳定在最大拉拔力。因为 数值模拟的有限性 界面模型和材料单元 ，拉拔力 达到最大后不再发生变化。 在三个互不相同的尾矿含水率的前提之下，不 同的上覆压力下相应的土工格栅最大拉拔力的数值 如图3 所示。由图3 中可知，在不相同的上覆压力 作用下，土工格栅最大拉拔力跟着上覆压力增加已 呈现线性增加；在相同的上覆压力作用下，随着尾矿 含水率这一指标的增加最大拉拔力先增大后减小， 如上覆压力为2 ．0 4k P a 时，含水率从1 ．5 ％增大到 9 ．7 ％时，最大拉拔力增幅6 ．2 1 ％，而含水率从 9 ．7 ％增大到1 9 ．3 ％时，最大拉拔力降幅5 ．9 8 ％，这 是由于尾矿含水率在增加，尾矿颗粒间以及其与格 栅间的相互作用会由胶凝作用逐渐向润滑作用 改变。 图2 拉拔位移与拉拔力变化关系 含水宰1 ．5 ％．上覆压力1 ．0 2k P a F i 昏2R e l a t i o 璐h i pb e t w ∞np u I 卜伽t 曲p I a o e 呲n ta n dp u I 卜伽tf o r ∞ 万方数据 有 色金 属 矿山部分 第7 2 卷 Z 芒 R 蜓 d _ K 嘣 上覆压力／k P a 图3 最大拉拔力随上覆压力变化的数值模拟结果 F i g ．3 N u m e r i c a ls i m u l a t i o nr e s u I t so ft h em a x i m u mp u I l o u t f o r c eV a r y i n gw i t ho V e r l y i n gp r e s s u r e 3 数值计算结果与室内试验数据对比 分析 3 ．1 室内拉拔试验结果分析 不同含水率下最大拉拔力随上覆压力变化关系 的室内试验结果[ 13 。1 4 ] ，如图4 所示。根据图4 我们 可以知道，土工格栅承受的最大拉拔力在不一样的 含水率现状下随着上覆压力变大而增大，表现出线 性增长的趋势，在同样的上覆压力条件下，最大拉拔 力随着尾矿含水率的增大也先增大后减小，如上覆 压力为2 ．0 4k P a 时，含水率从1 ．5 ％增大到9 ．7 ％、 1 9 ．3 ％，土工格栅最大拉拔力先增大9 ．0 6 ％后减 小3 ．2 5 ％。 进行适当的调整、试算，将试验值与通过反算得到的 模拟值进行对比分析。如图5 所示为土工格栅最大 剪应力的试验值与模拟值对比情况。从图5 可知， 试验测得数据和土工格栅与尾矿填料的数值计算结 果规律基本吻合，说明数值模拟与实验室测得的结 果可以相互验证分析。 5 8 5 ．6 54 謇52 巷5o 稼 蕊4 - 8 4 ．6 4 4 4 ．2 6 0 58 56 謇5 ．4 謇，， 密 蔷s ．o 4 ．8 4 6 4 4 6 0 5 ．8 5 ．6 生54 型52 楹 蕊s ．o 4 ．8 4 ．6 上覆压力／l P a r a 、w 1 ．5 ％ 上覆压力，k P a r } 1 w 9 ．7 ％ 3 ．2 数值计算结果与室内试验数据对比分析 数值计算过程中对格栅一尾矿的界面摩擦参数 1 ．0 22 ．0 4 上覆压力，k P a c l w 1 9 ．3 ％ 图5 最大剪应力的模拟值与试验值对比 F i g ．5C o m p a r i s O nb e t w e e ns i m u l a t e dV a l u e sa n d e x p e r i m e n t a lV a l u e so ft h em a x i m u ms h e a rs t r e s s 万方数据 第4 期高太等土工格栅与尾矿界面摩擦特性的拉拔数值试验研究 1 2 7 如图6 所示为不同含水率情况下法向应力与最 大剪应力变化关系的试验结果与数值计算结果对比 情况，法向应力与最大剪应力的拟合结果呈线性关 系，该拟合直线符合莫尔一库伦定律，由此可得出拉 拔模拟试验的界面摩擦参数似黏聚力和似摩擦角。 法向应力，k P a a 含水率15 ％ 法向应力，I 【P a ㈤含水率97 ％ 法向应力，l P a c 含水率1 9 ．3 ％ 图6不同含水率下法向应力与最大剪应力变化关系 F i g ．6R e l a t i o n s h i pb e t w e e nn o r m a ls t r e s sa n d t h em a x i m u ms h e a rs t r e s su n d e rd i f f e r e n tw a t e rc o n t e n t 土工格栅与尾矿界面摩擦参数似黏聚力和似摩 擦角的试验值和模拟值的比较关系，如图7 所示。 从图7 中能够看出，格栅尾矿界面摩擦参数的变 动规律都和室内试验数值变动规律相同，其中似黏 聚力的试验值与模拟值差距比较小，试验值比模拟 值小8 ％左右，似摩擦角的试验值与模拟值差距比 较大，试验值比模拟值大3 0 ％左右，这种情况出现 的原因认定为是室内试验误差及数值模拟边界条件 为理想边界条件导致，所以在土工格栅与尾矿的数 值计算过程当中，界面摩擦参数似黏聚力和似摩擦 角要适度的进行调整，似黏聚力数值可取1 ．1 倍试 验值，似摩擦角数值可取o ．7 倍试验值。 a 似黏聚力 含水率，％ 1 ， 似摩擦角 图7界面摩擦参数模拟值与试验值的对比关系 F i g ．7C o m p a r i s o nb e t w e e ns i m u I a t e dV a I u e sa n d e x p e r i m e n t a lv a l u e so fi n t e r f a c ef r i c t i o np a r a m e t e r s 4结论 1 试验测得数据和土工格栅与尾矿数值计算 结果基本相符，考虑到无法避开的室内试验误差 及数值模拟参数取值上的偏差，认为室内试验和 数值模拟都能够很好地反映出格栅一尾矿的界 面摩擦特性，试验结果和数值计算结果能够互相 检验分析。 万方数据 1 2 8 有 色金 属 矿山部分 第7 2 卷 2 认为F L A C 3 D 数值计算可以适度的响应室内 试验测得的界面摩擦参数，但是界面摩擦参数需进 行一定的调节，在土工格栅与尾矿数值计算过程当 中，其中似黏聚力数值可取1 ．1 倍试验值，似摩擦角 数值可取o ．7 倍试验值。 参 考 文 献 [ 1 ] [ 2 ] [ 3 ] [ 4 ] [ 5 ] [ 6 ] [ 7 ] 何友芳．土工格栅与细粒尾矿界面作用特性试验研究及数值 分析[ D ] ．重庆重庆大学，2 0 0 5 ． H EY o u f a n g ．T h ee x p e “m e n t a ls t u d yo ni n t e r f a c ep r o p e r t i e s b e t w e e ng e o g r i da n dg r a i nt a i l i n ga n dn u m e “c a la n a l y s i s [ D ] ． C h o n g q i n g C h o n g q i n gU n i v e r s i t y ，2 0 0 5 ． 陈建峰，李辉利，柳军修，等．土工格栅与砂土的细观界面特性 研究[ J ] ．岩土力学，2 0 1 1 ，3 2 增刊1 6 6 7 1 ． C H E NJ i a n f e n g ，L IH u i l i ，L 1 U 】u n x i u ，e ta 1 ．M e s o s c o p i cs t u d y o “n t e r f a c ep r o p e r t i e so fg e o g r i d r e i n f o r c e ds o i l [ J ] ．R o c ka n d S o i lM e c h a n i c s ，2 0 1 1 ，3 2 S 1 6 6 7 1 ． 王家全，王宇帆，黄世斌．等．循环荷载作用下土工格栅剪切特 性的颗粒流细观分析[ J ] ．水利学报，2 0 1 4 ，4 5 9 1 0 8 2 1 0 9 0 ． W A N GJ i a q u a n ，W A N GY u f a n ，H U A N GS h i b i n ，e ta 1 ．T h e p a r t i c l ef l o wm e s o s c o p i ca n a l y s i so fg e o g “ds h e a rp r o p e r t i e s u n d e rc y c l i cl o a d i n g [ J ] ．J o u r n a lo fH y d r a u l i cE n g i n e er i n g ， 2 0 1 4 ，4 5 9 1 0 8 2 1 0 9 0 ． 万亮．土工格栅与膨胀土相互作用的拉拔试验研究[ D ] ．长沙 长沙理工大学，2 0 0 7 ． W A NI ，i a n g ．R e s e a r c ho nt h eI n t e r a c t i o nb e t w e e ne x p a n s i v 。 s o i la n dg e o g r i d sb yp u l l o u tt e s t s [ D ] ．c h a n g s h a c h a n g s h a U n i v e r s i t yo fS c i e n c ea n dT e c h n o l o g y ，2 0 0 7 ． 李首龙．黏性土中土工格栅拉拔特性试验及数值模拟分析研 究[ D ] ．大连大连理工大学，2 0 0 9 ． L IS h o u l o n g ．E x p e r i m e n t a la n dn u m e r i c a ls i m u l a t i o ns t u d yo n p u l l o u tb e h a v i o ro fg e o g r i du n d e rc l a y [ D ] ．D a l i a n D a “a n U n i v e r s i t yo fT e c h n o l o g y ．2 0 0 9 ． 高文华，R I C H A R DJB ．条形荷载下加筋土边坡破坏机制的数 值模拟口] ．土木工程学报，2 0 0 7 6 5 4 5 8 ． G A 0W e n h u a ，R I C H A R DJB ．N u m e r i c a ls i m u l a t i o no ff a i l u r e m e c h a n i s mf o rr e j n f o r c e da n du n r e i n f o r c e ds l o p e su n d e rs t r i p f o o t i n g 口] ．c h i n ac i v i lE n g i n e e r i n gJ o u r n a l ，2 0 0 7 6 5 4 5 8 ． 杜勇立．土工格栅加筋红砂岩粗粒土挡墙试验研究和数值分 析[ D ] ．长沙湖南大学，2 0 1 0 ． D UY o n g l i ．T e s ta n a l y s i sa n dn u m er i c a ls i m u l a t i o no f r e i n f o r c e dr e ds a n d s t o n ee a r t hr e t a i n i n gw a l lb yg e o g r i d [ 【 ] ． C h a n g s h a H u n a nU n i v e r s i t y ，2 0 1 0 ． [ 8 ]c H E Nc ，M c D 0 w E L LGR ，T H O MNH ．D i s c r e t ee l e m e n t m o d e l l i n go fc y c l i cl o a d so fg e o g r i d r e i n f o r c e db a I l a s tu n d e r c o n f i n e da n du n c o n f i n e dc o n d i t i o n s [ J ] ．G e o t e x t 订e sa n d G e o m e m b r a n e s ，2 0 1 2 ，3 5 7 6 8 6 ． [ 9 ] 赖汉江，郑俊杰，甘甜，等．土工格栅加筋拉拔试验界面特性细 观分析[ J ] ．土木工程与管理学报，2 0 1 3 ，2 9 4 4 5 4 9 ． L A IH a n j i a n g ，Z H E N GJ u n j i e ，G A NT i a n ，e ta 1 ．M i c r o s c o p i c a n a l y s i so fi n t e r f a c ec h a r a c t e r i s t i c so fg e o g r i da n ds o i lt h r o u g h p u l l o u tt e s t [ J ] ．J o u m a lo fc i v 订E n g i n e e r i n ga n dM a n a g e m e n t ， 2 0 1 3 ，2 9 4 4 5 4 9 ． [ 1 0 ] 周健，} L 祥利，鞠庆海，等．土工合成材料与土界面的细观研 究口] ．岩石力学与工程学报，2 0 0 7 ，2 6 增刊1 3 1 9 6 3 2 0 2 ． Z H O UJ i a n ，K O N GX i a n g l i ，J UQ i n g h a i ，e ta 1 ．M e s o s c o p i c a l s t u d yo ni n t e r f a c eb e t w e e ng e o s y n t h e t i c sa n ds o i l [ J ] ．c h i n e s e J o u r n a lo fR o c kM e c h a n i c sa n dE n g i n e er i n g ，2 0 0 7 ，2 6 S 1 3 1 9 6 3 2 0 2 ． [ 1 1 ] 苗晨曦，郑俊杰，崔明娟，等．三向土工格栅筋土界面及摩擦特 性的离散元模拟[ J ] ．岩土力学，2 0 1 4 增刊1 4 2 3 4 3 0 ． M I A OC h e n x i ，Z H E N GJ u n j i e ，C U IM i n 西u a n ，e ta 1 ．D i s c r e t e e l e m e n ts i m u l a t i o no fi n t e r f a c ec h a r a c t e r i s t i ca n df “c t i o n a l p r o p e r t i e so ft “a x i a lg e o g r ．d [ J ] ．R o c ka n ds o i lM e c h a n i c s ， 2 0 1 4 S 1 4 2 3 4 3 0 ． [ 1 2 ] 郑俊杰，苗晨曦，谢明星，等．界面特性及填料粒径影响格栅加 筋性能的离散元研究[ J ] ．岩土工程学报，2 0 1 3 ，3 5 8 1 4 2 3 1 4 2 8 ． Z H E N GJ u n j i e ，M l A OC h e n x i 。X I EM i n g x i n g ，e ta 1 ．I n t e r f a c 8 p r o p e r t i e s a n di n f l u e n c eo f p a r t i c l e s i z eo n g e o g r i d r e i n f o r c e m e n tp e r f o r m a n c eb yD E M [ J ] ．C h i n e s eJ o u r n a lo f G e o t e c h n i c a lE n g i n e e r i n g ，2 0 1 3 ，3 5 8 1 4 2 3 1 4 2 8 ． [ 1 3 ] 杜镀。张向阳，王肖辉，等．运用非线性M 。h r - C o u l o m b 破坏准则分 析尾矿坝稳定性[ J ] ．有色金属 矿山部分 ，2 0 1 9 ，7 1 5 8 1 8 4 ． D UD u ，Z H A N GX i a n g y a n g ，W A N GX i a o h u i ，e ta 1 ．S t a b i l i t y a n a l y s i so ft a i l i n g sd a mu s i n gt h en o n l i n e a rM o h r C o u l o m b f a i l u r ec “t e r i o n [ J ] ．N o n f e r r o u sM e t a l s M i n i n gs e c t i o n ，2 0 1 9 ， 7 1 5 8 1 8 4 ． [ 1 4 ] 谢宝珊，张向东，杜常博．不同含水率下尾矿与土工格栅界面 特性的试验研究[ J ] ．实验力学．2 0 】8 。3 3 1 1 2 7 1 3 3 ． X I EB a o j i n ，Z H A N GX i a n g d o n g ，D UC h a n g b o ．E x p er i m e n t a l s t u d yo ft h ei n t e r f a c ec h a r a c t e “s t i c sb e t w e e ng e o g “da n d t a i l i n g s w i t hd i f f e r e n tw a t e rc o n t e n t [ J ] ． J o u r n a l o f E x D e r i m e n t a lM e c h a n i c s ，2 0 1 8 ，3 3 1 1 2 7 1 3 3 ． 万方数据