DreamBlooms/HansPub-15000-to-45000 huggingface dataset:
system: [SYSTEM]

user: 5. 结论与建议 通过分析，BIM在建筑领域处于引入期向增长期过渡，相容性方面，我们可以看出相比较传统CAD，BIM在建筑全寿命周期都具有一定的相对优势，且BIM与当下的建筑理念和需求是一致的，且BIM可以和传统CAD协同作业，更好的为建筑业服务。同时通过易懂性分析可以看出BIM理念是较易理解，但软件操作稍难，同时也说明BIM作为新兴技术还是存在一定的问题，使得难以推广。这些问题包括BIM的不确定性的风险，经营管理层对BIM新模式转变困难是不确定性风险的主要因素。但是总体而言，不论是设计人员还是施工人员都有较高采用BIM技术的意愿。针对本次研究，本着可操作的原则，提出如下建议。 5.1. 以BIM使用者而言 这里的使用者不仅仅指个人，也指企业、单位甚至是国家。BIM使用者不仅仅是作为BIM产品带来好处的享受者，更应当具有推广BIM技术的责任，目前BIM人才本就缺乏，作为新技术的受益者将BIM技术更好地运用在日常的工作中的同时还应当发掘新的创新点，这个创新点包括BIM技术的应用点和发展可能采纳BIM技术的人员。使用者内部也可以有相应的激励奖惩制度，为培养人才提供有力保障。 5.2. 以BIM创造者而言 当然这里的创造者也不仅仅指个人。所谓的BIM创造者不是指创造BIM的人，而是指对BIM软硬件进行开发(包括二次开发)的人员和编制BIM标准的人员。在专心研究使BIM技术更加成熟的同时，可以拓宽眼界融合其他技术、学科、领域发展BIM。 5.3. 以BIM推广者而言 BIM的推广者指无论是否从事BIM相关工作，有意愿使用BIM为建筑领域添砖加瓦的人员。在推广BIM应寻找合适、合理、有效的办法。研究熟悉市面上各类BIM软件及BIM方案的特性及优缺点，并针对项目需求，为企业提出合理建议，使其自愿学习利用BIM，同时自愿加入BIM的推广者中。
system: [SYSTEM]

user: 摘 要 针对定边采油厂某注水项目区域内主采层位为三叠系长8层酸化常规土酸酸化反应速度很快、作业距离短、二次、三次伤害严重、过度溶蚀、腐蚀严重等问题，开发了一种包含盐酸、氢氟酸、有机酸、氟硼酸以及多种添加剂，具有低腐蚀、低伤害、低反应速度、配伍性良好性能的新型酸液，并对新型酸液及添加剂性能进行了评价，结果显示：相比常规土酸，新型酸液缓速效果明显，适合深度酸化；同等添加量时，H03型缓蚀剂的缓蚀效果明显优于与H01、H02型缓蚀剂，且相同缓蚀效果条件下，H03型缓蚀剂的用量最少；螯合剂添加量1500 ppm时，稳铁能力达到最大值；表面活性剂HW-3用量500~2000 ppm，新型酸液与原油表面张力维持在18 mN/m左右，为维持一定的安全余量，3%新型酸液中添加1000 ppm的HW-3型表面活性剂。 关键词 :降压增注，深度酸化，缓蚀剂，铁离子稳定剂，表面活性剂 Copyright © 2020 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 工区概况 定边采油厂某注水项目区域内主采层位为三叠系长8层，油藏埋深2500~2800 m，主要为特低孔超低渗岩性油藏。该区块目前开发中存在问题，主要有以下几点：初期产量高，递减快，油井产量低；地层物性差，低孔低渗，注水压力高，注入难度大；地层致密，注入压力高，达不到配注要求，没有建立有效水驱关系，注水开发效果差；常规土酸酸化解堵、小规模压力措施效果不大，注入压力仍居高不下，存在主要问题包括 [ 1 ] [ 2 ] [ 3 ] [ 4 ] [ 5 ]：①常规土酸酸液与砂岩油藏中粘土的初始反应速度很快，导致酸液作用距离过短，地层中深部的堵塞得不到解除，对地层的改造不彻底。②常规土酸酸液与砂岩油藏中各种成分反应副产物氟硅化物和氟铝化物发生二、三次反应，形成一系列不溶性氟硅酸盐沉淀物，影响常规土酸酸化效果。③常规土酸液中都是以强酸为主，对与结构较疏松的砂岩油藏易过度溶蚀，造成近井地带井壁坍塌。④常规土酸酸化体系还易对施工设备、井下管柱、井筒等造成腐蚀。 因此针对该区域下组合深部油藏高压注水急需展开系统分析，优选合适酸液，达到降压增注目的，提高注水效果。针对定边采油厂某注水项目区酸化需要，研发了一种低腐蚀、低伤害、低反应速度、配伍性良好的新型酸液体系，并对酸液体系性能进行了评价。
system: [SYSTEM]

user: 2. 酸液选择依据 酸化主要是消除近井地带的污染，所选择的酸液要求在有效消除储层原伤害的同时，不会对储层造成二次伤害。 2.1. 砂岩类型分析 为探究定边采油厂某注水项目区域内层系岩石特点，开展了2组岩心薄片分析，着重了解碎屑岩的物质成分和结构，鉴定岩石中矿物成分，估计各矿物的百分含量，为选择合适的酸液提供理论依据。薄片分析正交偏光显微示图如图1、图2所示。 图1. 试样1岩心薄片正交偏光显微示图 图2. 试样2岩心薄片正交偏光显微示图 从图1、图2可以看出，岩石为砂状碎屑结构，颗粒支撑，孔隙式胶结，碎屑成分主要为石英砂岩，部分长石砂岩，少量岩屑。填隙物主要为粘土矿物。 2.2. 伤害类型 储层的伤害主要包括：钻井液固相颗粒和滤液浸入、粘土膨胀和微粒运移、地层结垢(CaCO 3 、CaSO 4 、BaSO 4 等)、有机沉积物、乳化堵塞、润湿性改变、细菌、水锁等 [ 6 ] [ 7 ] [ 8 ] [ 9 ]，对于不同伤害类型的储层，要选择相应的酸液体系。 将盐酸、氢氟酸、少量氟硼酸、有机酸等按一定比例及配伍性添加剂混合而成新型酸液。当盐酸足量时有机酸几乎不参与反应；当盐酸与储层矿物发生反应，有效浓度低时，有机酸才与储层矿物发生缓慢反应，延缓了氢氟酸的消耗，随着氢氟酸的消耗量不断增加，氟硼酸水解反应速度增加，生成氢氟酸，维持较低的pH值，增加了酸液穿透半径，从而达到深度酸化目的。一般乙酸或甲酸的浓度为0%~9%，盐酸的浓度为5%~15%，氢氟酸的浓度为0%~4%。氟硼酸的浓度为0%~3%。
system: [SYSTEM]

user: 4. 结论 新型酸液酸化是注水井降压增注的有效措施，相比常规土酸，新型酸液缓速效果明显，适合深度酸化；同等添加量时，H03型缓蚀剂的缓蚀效果明显优于H01、H02型缓蚀剂，且相同缓蚀效果条件下，H03型缓蚀剂的用量最少；螯合剂添加量1500 ppm时，稳铁能力达到最大值；表面活性剂HW-3用量500~2000 ppm，新型酸液与原油表面张力维持在18 mN/m左右，为维持一定的安全余量，3%新型酸液中添加1000 ppm的HW-3型表面活性剂。
system: [SYSTEM]

user: 摘 要 本文以长沙中建梅溪一号(3B地块)基坑工程为例，通过前期的地质调查、多方案论证，确定了初步的土钉墙支护结构和放坡挂网素喷的基坑支护方案，应用FLAC3D对支护方案模拟，模拟分析了采用土钉墙支护条件下，基坑开挖过程中基坑土体的应力场分布规律、支护结构受力状况，模拟结果显示：1) 从土钉受力情况来看，土钉轴力基本为拉力。同时，随着开挖工程的开展，土钉受力情况是随之而发生变化的，数值模拟结果显示，开挖过程中，土钉所承受的最大轴力为5.2 t。2) 土钉在长度方向上轴力分布遵循“中性点”原则，即中间数值大、两边数值小。最大轴力位置在潜在滑移面处。靠近坡面土钉段的轴力大于土钉深部段轴力这说明，土钉除了承受土体变形后施加的动力外，还受到喷射混凝土面层施加在土钉上的抗拔力作用。3) 随着开挖的进行，土钉受力逐渐增大。这主要是由于开挖深度增加，基坑边坡土体变形增加，将土体的位移力和土体膨胀力传递给土钉所致。 关键词 :基坑支护，土钉墙支护，应力状态 Copyright © 2020 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 土钉是一种基于新奥隧道法原理发展起的一种基坑支护技术 [ 1 ] - [ 8 ]。新奥隧道法(The New Austrian Tunneling Method)最早见于六十年代初期，到六十年代中期，随着经验施工观测资料的积累与施工经验的丰富，陆续将此法从硬质岩层推广到含软岩层、粉土、砾石和砂石的隧道施工当中 [ 9 ] [ 10 ]。随着计算机技术的发展及数学理论的不断提升，有关土钉支护的研究方法主要有以下四种：1) 试验方法 [ 11 ]；2) 极限平衡法 [ 12 ] [ 13 ] [ 14 ]；3) 有限元法 [ 15 ]；4) 有限差分法 [ 16 ] [ 17 ] [ 18 ]。试验方法的开展大量的资金投入，对于一般的科研机构由于资金不足而不能够系统的展开试验研究，尤其是大型足尺试验的开展更为不现实，而对于生产单位则为了追求经济效益，一般也不愿意将资金投入到此方面。极限平衡分析法在土钉有关变形信息与内力精确分布两大方面存在一定的不足，从而难以全面评估土钉支护的工作性能及安全性。计算机技术的发展，既节约了大量的资金，又能够通过计算获得土钉支护内力、变形等较全面详细的数据分析，因此通过计算机数值模拟可以为基坑工程的设计和施工提供指导性依据 [ 19 ] [ 20 ] [ 21 ] [ 22 ] [ 23 ]。而利用有限元开发的大型通用软件ANSYS和利用有限差分原理开发的FLAC3D软件相比，后者运算速度快，是专门的岩土模拟程序，针对性强 [ 24 ] [ 25 ] [ 26 ]。 本文长沙中建梅溪一号(3B地块)基坑支护方案为例，利用FLAC3D程序对其位移场及稳定性进行数值模拟分析，并与现场实测数据结果进行比较，以期得到一些有参考意义的结论，为以后本地区及相似工程地质条件的基坑支护研究工作提供一定的参考。
system: [SYSTEM]

user: 2. 地质条件概述 2.1. 场地地层条件 场地内共涉及11层土，分别为第1层：杂填土(Q 4-3 m1 )；第2层：新近沉积粉土(Q 4-3 a1 )；第3层：新近沉积粉质粘土夹粉土(Q 4-3 a1 )；第4层：新近沉积粉质粘土(Q 4-3 a1 )；第5层(Q 4-2 1 )：粉土；第6层：粉质粘土(Q 4-2 1 )；第7层：粉土(Q 4-2 1 )；第8层：粉质粘土(Q 4-2 1 )；第8夹层：粉土(Q 4-2 1 )；第9层：粉土(Q 4-2 1 )；10-1层：粉砂(Q 4-1 a1 )；第10层：细砂(Q 4-1 a1 )；第11层：粉质粘土(Q 3 a1 )。 2.2. 场地水文地质条件 前期工程勘查资料表明，勘探深度范围内场地的地下水可分为两层，上层为潜水，下部为承压水，两者被相对隔水层第8层粉质粘土隔开。场地内自然水位约为地表下1.5 m。根据规范要求以及周边工程施工经验，建议上部潜水含水层的综合渗透系数取K = 0.5 m/d，下部承压水含水层的综合渗透系数取K = 8.0 m/d。 2.3. 主要土层物理力学指标 场地内各主要土层的物理力学指标见表1。 Table 1 Physical and mechanical indicators of major soil layer 1 2 3 4 5 6 7 8 9 土层名称 杂填土 沉积粉土 粉质粘土 沉积粉质粘土 粉土 粉质粘土 粉土 粉质粘土 粉土 厚度(m) 0.3~2.35 2.43 1.84 2.46 3.39 1.69 1.83 3.12 1.44 体积模量(MPa) 4.3 4.04 5.59 6.44 8.05 8.33 6.35 9.32 5.59 剪切模量(MPa) 1 1.1 2.29 2.63 3.29 3.20 2.44 4.3 2.29 质量密度 (×10 3 Kg/m 3 ) 1.7 1.9 1.9 2.0 2.0 2.0 1.9 2.0 1.9 粘聚力(KPa) 15 20 30 27 15 32 12 31 30 内摩擦角(˚) 15 20 20 21 21 20 20 20 20 剪胀角 0 0 0 0 0 0 0 0 0 表1. 主要土层物理力学指标
system: [SYSTEM]

user: 3. 基坑支护设计 本工程的基坑支护对三个1-1’、2-2’、3-3’剖面(图1，图2，图3)采用土钉支护的方案，对4-4’剖面采用放坡挂网素喷的支护方式，土钉支护方案的支护参数见表2，表3，表4；素喷坡面的设计参数见图4。
system: [SYSTEM]

user: 4. 应力–应变规律 当各个三角形内的应变张量确定以后，然后初始设定的本构规律计算相应的应力。在此之前，需要对原有的应力进行两个方面的修正：1) 转动应力作修正；2) 扣除孔隙压力。应变与应力间的本构关系给成增量形式： Δ σ i j = f ( Δ l i j ， σ i j ， ⋯ ) (1) 式中： Δ σ i j ——应力增量张量， i ， j = 1 ， 2 ； f ( ) ——包含有增量应变，现存总应力状态，材料常数等的本构方程。 图1. 剖面1-1’示意图 图2. 剖面2-2’示意图 图3. 剖面3-3’示意图 Table 2 Soil nailing support parameters of 1-1’ profil 排号 埋置深度(m) 土钉长度(m) 钢筋配置 水平间距(m) 倾角(˚) 成孔直径(mm) 1 1.3 8.0 1φ18 1.3 5~15 100 2 2.8 10.0 1φ18 1.3 5~15 100 3 4.3 9.0 1φ18 1.3 5~15 100 4 5.8 9.0 1φ18 1.3 5~15 100 5 7.3 6.0 1φ18 1.3 5~15 100 表2. 剖面1-1’土钉支护参数 Table 3 Soil nailing support parameters of 2-2’ profil 排号 埋置深度(m) 土钉长度(m) 钢筋配置 水平间距(m) 倾角(˚) 成孔直径(mm) 1 1.3 8.0 1φ18 1.3 5~15 100 2 2.7 9.0 1φ18 1.3 5~15 100 3 4.1 8.0 1φ18 1.3 5~15 100 4 5.5 6.0 1φ18 1.3 5~15 100 表3. 剖面2-2’土钉支护参数 Table 4 Soil nailing support parameters of 3-3’ profil 排号 埋置深度(m) 土钉长度(m) 钢筋配置 水平间距(m) 倾角(˚) 成孔直径(mm) 1 1.3 6.0 1φ18 1.3 5~15 100 2 2.8 9.0 1φ18 1.3 5~15 100 3 4.3 8.0 1φ18 1.3 5~15 100 4 5.8 6.0 1φ18 1.3 5~15 100 5 7.3 4.5 1φ18 1.3 5~15 100 表4. 剖面3-3’土钉支护参数 4.1. 网点处不平衡力的确定 一旦四边形域内的应力张量已经求得，便可以用它来计算各网点处的不平衡力。重写运动方程式 [ 24 ] [ 27 ]： p [ ∂ u ˙ i ∂ τ ] = ∂ σ i j ∂ x j + p g i (2) 其结果为 p [ ∂ u ˙ i ∂ τ ] = 1 A ∑ σ i j ε j k Δ x k + p g i (3) 或 ∂ u ˙ i ∂ τ = F i m + g i 式中： F i = ∑ σ i j ε j k Δ x k ，     m = p A 因此网点力变为： F i = σ i j a ε j k ( x k ( 1 ) − x k ( 4 ) ) + σ i j b ε j k ( x k ( 2 ) − x k ( 1 ) ) + σ i j c ε j k ( x k ( 3 ) − x k ( 2 ) ) + σ i j d ε j k ( x k ( 4 ) − x k ( 3 ) ) (4) 这个力在所有应力和坐标，在时刻t，t + Δt等时都是已知的。在下一个半时间步长时的网点速度为： u ˙ i ( t + Δ t 2 ) = u ˙ i ( t − Δ t 2 ) + F i t m Δ t + g i (5) 应该指出，如果域中出现任何一个错误，那么出现错误的域中的各项将由式(4)中直接省略掉。 如果模拟过程中，选择尺度较大的应变模型，则把修正以后的网格坐标用于式(5)的，求得各新速度： x i ( t + Δ t ) = x i ( t ) + u ˙ i ( t + Δ t 2 ) Δ t (6) 4.2. 应力转动修正项 当利用代码中的大应变逻辑时，由于物体的转动，必须修正应力。应力的正规转动方程式为： σ ¯ i j = σ α β J i α J j β (7) J i j = [ cos θ sin θ − sin θ cos θ ] (二维) 式中： σ ¯ i j ——大应变转动时修正了的应力张量。 当角度 Δ θ 非常小的条件下，因为 cos θ ≈ 1 ， sin θ ≈ Δ θ ，则 J i j = [ 1 Δ θ − Δ θ 1 ] (二维) 因此， σ ¯ 11 = σ 11 J 11 J 11 + σ 21 J 12 J 11 + σ 12 J 11 J 12 + σ 22 J 12 J 12           = σ 11 + σ 21 Δ θ + σ 12 Δ θ + σ 22 Δ θ 2 Δ σ 11 = 2 σ 12 Δ θ ( 因 为 σ 21 = σ 12 ， Δ θ 很 小 时 ， Δ θ 2 = 0 ) (8) 特殊条件下， σ 11 的修正项为 Δ σ 11 = − 2 σ 12 Δ θ (9) 同理，有 σ ¯ 22 = σ 11 J 21 J 21 + σ 21 J 22 J 21 + σ 12 J 21 J 22 + σ 22 J 22 J 22 = σ 11 Δ θ 2 − σ 21 Δ θ − σ 12 Δ θ + σ 22 2 = σ 22 − 2 σ 12 Δ θ + σ 11 Δ θ 2 (10) 反号并略去 Δ θ 2 项后，有 Δ σ 22 = − 2 σ 12 Δ θ (11) 还有 σ ¯ 12 = σ 11 J 11 J 21 + σ 21 J 12 J 21 + σ 12 J 11 J 22 + σ 22 J 12 J 22 = − σ 11 Δ θ − σ 21 Δ θ 2 + σ 12 + σ 22 Δ θ (12) Δ σ 12 = ( σ 11 − σ 22 ) Δ θ (13) 当采用大应变逻辑时，把这些修正项加到各应力中。
system: [SYSTEM]

user: 5. 基坑模型及参数设置 5.1. 模拟边界设置 本基坑位于长沙市岳麓区梅溪湖片区境内，南邻梅溪湖路，东邻E2路，场地交通便利，地理位置优越。实际基坑为长方形，基坑外轴线面积约102 × 371.9 m 2 ，南北侧开挖深度为地表下6 m，中部为8 m。根据场地工程地质条件，基坑支护方案拟采用以下两种形式：1) 土钉墙支护；2) 放坡挂网素喷支护方案。据工程经验与现场实际情况，本次数值模拟模型尺寸设置为30 m × 20 m × 20 m。 5.2. 模型参数设置 模型中岩土物理力学参数按基坑支护方案设计中参数取值。 5.3. 计算模型 本次共建立三个力学模型，分别是工况一基坑开挖6 m (设计方案剖面2，位于基坑北侧)、工况二基坑开挖深8 m (设计方案剖面1，位于基坑东侧)、工况三基坑开挖深8 m (设计方案剖面3，位于基坑西侧)。各剖面数值模拟模型图分别见图4~6。 图4. 剖面二模型图 图5. 剖面一模型图 图6. 剖面三模型图 5.4. 支护结构的模拟 土钉采用FLAC3D中的锚索结构(cable)，喷射混凝土用壳结构(shell-type)。土钉及混凝土层面的参数如表5和表6所示。 Table 5 Soil nail material parameter 弹性模量/GP 截面积/m 2 屈服强度/KPa 水泥黏结刚度/N/M 水泥黏结强度/KPa 20 2.54 × 10 −4 800 4 × 10 9 89 表5. 土钉材料参数 Table 6 Concrete parameter 弹性模量/GP 泊松比 2.1 0.28 表6. 混凝土参数
system: [SYSTEM]

user: 6. 数值模拟结果及分析 6.1. 土层初始应力状态模拟 在长期自重应力作用下，基坑在开挖施工前，施工区的地层土体已经处于稳定状态。所以在模型过程中，需要对模型进行应力初始化。开挖计算前，位移场和速度场需要清零，以避免出现误差。见图7所示。 从模拟结果来看，研究区基坑深度范围内土体垂直应力呈现由上至下呈递增分布的规律，在地表附近约为−0.05 MPa (负号表示为压力)，直至基坑深处时，其数值大约为−0.38 MPa。模拟结果与现场土层的应力分布状态基本吻合。 6.2. 土钉轴力分析 1) 剖面1土钉轴力(图8~13) 图7. 应力初始化云图 图8. 开挖1.5 m土钉轴力图 图9. 开挖3 m土钉轴力图 图10. 开挖4.5 m土钉轴力图 图11. 开挖6 m土钉轴力图 图12. 开挖8 m土钉轴力图 图13. 开挖稳定后土钉轴力曲线 2) 剖面2土钉轴力(图14~18) 图14. 开挖1.5 m土钉轴力图 图15. 开挖3 m土钉轴力图 图16. 开挖4.5 m土钉轴力图 图17. 开挖6 m土钉轴力图 图18. 开挖稳定后土钉轴力曲线 3) 剖面3土钉轴力(图19~24) 图19. 开挖1.5 m土钉轴力图 图20. 开挖3 m土钉轴力图 图21. 开挖4.5 m土钉轴力图 图22. 开挖6 m土钉轴力图 图23. 开挖8 m土钉轴力图 图24. 开挖稳定后土钉轴力曲线 对基坑开挖过程中土钉受力进行了实时监测如图8~24所示，得出如下结论： 1) 从土钉受力情况来看，土钉轴力基本为拉力。同时，随着开挖工程的开展，土钉受力情况是随之而发生变化的，数值模拟结果显示，开挖过程中，土钉所承受的最大轴力为5.2 t。 2) 土钉在长度方向上轴力分布遵循“中性点”原则，即中间数值大、两边数值小。最大轴力位置在潜在滑移面处。靠近坡面土钉段的轴力大于土钉深部段轴力这说明，土钉除了承受土体变形后施加的动力外，还受到喷射混凝土面层施加在土钉上的抗拔力作用。 3) 随着开挖的进行，土钉受力逐渐增大。这主要是由于开挖深度增加，基坑边坡土体变形增加，将土体的位移力和土体膨胀力传递给土钉所致。
system: [SYSTEM]

user: 7. 结论 1) 从土钉受力情况来看，土钉轴力基本为拉力。数值模拟结果显示，开挖过程中土钉所承受的最大轴力为5.2 t。 2) 土钉在长度方向上轴力分布遵循“中性点”原则，即中间数值大、两边数值小。 3) 随着开挖的进行，土钉受力逐渐增大。
system: [SYSTEM]

user: 摘 要 本文将三维数字散斑相关方法(Digital Image correlation， DIC)引入纳米混凝土变形测量当中，基于数字散斑技术对纳米混凝土在单轴受压过程中的全场位移和应变进行采集，分析混凝土在不同纳米材料掺量、不同受力阶段下应变云图的变化特征和变形参数的变化规律，并对比不同纳米掺量混凝土力与位移曲线。结果表明：纳米含量为1.5%的混凝土于纳米含量为1%的混凝土相比，其强度降低了15.7%。其降低原因可能是由于随着纳米含量的增加，纳米二氧化硅在分散时出现团聚的可能性增加，从而导致混凝土强度降低。数字散斑相关方法可有效观测纳米混凝土试件在变形前后应变演化趋势，从而对试件裂纹的扩展进行计算。混凝土受压过程中，发生的应变突变现象是因为混凝土局部破坏导致。混凝土破坏是由横向裂缝和纵向裂缝共同作用。 关键词 :数字散斑，纳米混凝土，裂缝扩展，力与位移曲线 Copyright © 2020 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 随着经济的飞速发展，混凝土在现代生活中所占比例越来越重。近年来对混凝土研究地不断深入，人们发现混凝土内部的孔隙是影响混凝土强度，导致混凝土破坏的主要因素 [ 1 ]。所以很多学者选择在混凝土中掺加改性材料，以此来提高混凝土强度 [ 2 ]。二十世纪八十年代 [ 3 ] 以来随着光电子技术、数字图像处理和计算机等新技术的诞生，一种光学测力新方法——数字散斑相关方法(digital image correlation， DIC)形成发展起来。数字散斑相关方法是由光学原理结合力学理论，并加以计算机辅助，以实验的手段对混凝土的各项力学性能进行测定。数字散斑相关方法具有测量方便性、非接触性、测量高精度性 [ 4 ] 等特点，使得这项技术在工程材料领域 [ 5 ] [ 6 ]、医学、生物、国防等 [ 7 ] 各个领域获得了广泛地应用。Jennifer D’Anna等 [ 8 ] 研究了DIC技术在玄武岩纤维增强胶凝材料(BFRCM)条拉伸性能中的应用。Yinyin Li、Wei Jiang [ 9 ] 利用数字图像相关方法和积分理论，对不同质量的纳米碳化钨掺入量进行了激光修复，评价了其断裂性。洪哲、杨树 [ 10 ] 利用DIC技术研究了沥青混凝土断裂特性，结果表明裂缝尖端区域具有不规则、不唯一、不确定性。黄兴震等 [ 11 ] 利用电子散斑干涉技术(ESPI)和数字散斑技术(DIC)对普通混凝土和高强混凝土的断裂特性进行研究，结果表明ESPI和DIC均可代替传统测量仪器；当试件变形较大时应用DIC技术测量的误差较小且较稳定。李东洋 [ 12 ] 在混凝土断裂性能测试及描述方法研究中，表明DIC技术确定混凝土中裂纹尖端的方法是有效且可行的。时金娜等 [ 13 ]，利用DIC技术对不同玄武岩纤维掺量的混凝土进行单轴压缩实验进行分析，结果表明掺入玄武岩纤维可减小混凝土的变形。本文采用数字散斑技术(DIC)，对纳米二氧化硅改性混凝土压缩变形过程进行实时监测分析，研究试件在不同力作用下的变形量、试件表面应变场变化以及裂纹开展情况。
system: [SYSTEM]

user: 2. 实验概况 2.1. 试件制作 实验采用100 mm × 100 mm × 100 mm混凝土试块，纳米含量为1%、1.5%的二氧化硅纳米混凝土试块，其试验配合见表1： Table 1 Test mix rati 编号 W/C C/kg 天然骨料 纳米材料/g W/kg S/kg G/kg 15 nm SiO 2 NS1.0 0.5 8.70 14.28 26.54 87.84 4.39 NS1.5 0.5 8.65 14.28 26.54 131.76 4.39 表1. 试验配合比 注：表中的符号简写，W/C：水灰比，C：水泥，S：砂子，G：石子，W：水，NS：纳米SiO 2 混凝土。 2.2. 散斑制作 在本实验中采用人工散斑，首先在试件表面先喷洒一层白漆，再用装有黑色哑光漆的喷雾器在试件表面撒布黑色斑点。喷涂时应将黑色喷雾器提高一定高度水平喷射，使黑色雾状油气在自身重力作用下在试件表面形成随机且均匀的分布散斑点。散斑点的大小与灰度的对比度是评价人工散斑质量的主要指标，人工散斑的制作对灰度计算中的相关系数值及数字散斑相关测量的精确度具有一定影响。 2.3. 实验设备及装置 实验试件为混凝土试件，实验器材：MTS公司生产的万能伺服试验机、冷白光源、高精度低速摄像机、图像卡、计算机以及软件。试验装置和光路如下图1所示： 图1. 试验装置和光路 2.4. 原理 如图1所示，光源照射试件表面，根据试件表面喷涂的颗粒，形成随机分布散斑。用低速摄像机记录散斑图，试件表面散斑随着试件的变形而发生变化，分析变形前后散斑的位置变化，得到散斑沿着横轴和纵轴发生的相对位移。变形前后试件的散斑图有一定的相关性，若试件变形不大，其表面的散斑场灰度变形可忽略。设(x， y)为变形前点，(x * ， y * )为变形后点。二者关系如下： { x ∗ = x + u + Δ x ∂ u ∂ x + Δ y ∂ u ∂ y y ∗ = y + v + Δ x ∂ v ∂ x + Δ y ∂ v ∂ y (1) P ( x i ∗ ， y j ∗ ) 表示变形后对应点的灰度值，则二者关系为： C = ∑ i = 1 m s ∑ j = 1 m s [ f ( x i ， y j ) − f − ] [ g ( x i ∗ ， y j ∗ ) − g ¯ ] ∑ i = 1 m s ∑ j = 1 m s [ f ( x i ， y j ) − f ¯ ] 2 ∑ i = 1 m s ∑ j = 1 m s [ g ( x i ∗ ， y j ∗ ) − g ¯ ] 2 (2) 式中 f ¯ 、 g ¯ 分别为参考子目和目标子目的平均灰度值，当C = 1时，完全相关；C = 0时，无关。 2.5. 试验过程 图像采集系统由两台高分辨像机、图像采集板、计算机以及控制系统组成。利用相机进行图像采集时，摄像机同步到以恒定时间间隔捕捉被调查试样的同步图像，以保证荷载同变形场的演化与时间轴对应关系。系统中采用磁性表座将两台低速像机固定于试验机下台面，并使相机呈一定角度，保证低速相机与试件之间相对稳定性，相机轴线垂直于试件被测面，与试件表面交于形心位置，试验中采用轴向压缩的加载方式，设计加载速率为0.5 kN/s，加载至试件破坏。
system: [SYSTEM]

user: 3. 试验结果分析 3.1. 力与位移曲线分析 单轴受压过程中力与位移曲线，如图2所示。实验结果发现，纳米二氧化硅含量为1%时混凝土强度为36.9 MPa；纳米二氧化硅含量为1.5%时混凝土强度为31.1 MPa，与纳米含量为1%的混凝土相比其强度降低了15.7%。其降低原因可能是由于随着纳米含量的增加，纳米二氧化硅在分散时出现团聚的可能性增加，从而导致混凝土强度降低。 图2. 力与位移曲线 3.2. 应变分析 图3~6为不同纳米材料掺量在不同受压阶段的横向、纵向应变云图。分别为混凝土受压前期(a)、中期(b) (c)、后期(d)~(g)、试件破坏(h)四个阶段。 混凝土受压前期，如图3~6中图(a)所示，横、纵向云图中均有密压区出现，且两种含量的试件左下部都有应变集中区域产生，由此可知此时试件内部存在的原始孔隙已基本被挤压密实。 混凝土受压中期，如图3~6中(b) (c)所示，不同掺量的混凝土在受压中期时均出现了横向应变集中区域，由图中标尺可读出，含量为1%的纳米混凝土受压前期左下部的集中应变区域已经被破坏，随之延伸出两条应变加深带；含量为1.5%的纳米混凝土应变集中区域增大，从上下方向各延伸出一条应变加深带；因为混凝土本身就是带裂缝工作，所以试验继续进行。混凝土受压中期局部最大应变见表2，通过对不同含量纳米混凝土在中期局部最大应变值比较，结果表明应力较低时，含量为1%的混凝土对变形的抵抗能力较弱。 图3. 纳米掺量为1%的混凝土的横向应变云图 图4. 纳米掺量为1.5%的混凝土的横向应变云图 图5. 纳米掺量为1%的混凝土的纵向应变云图 图6. 纳米掺量为1.5%的混凝土的纵向应变云图 Table 2 Local maximum strain value in the middle stage of concrete compressio 纳米含量 横向局部最大应变 纵向局部最大应变 1% 0.5235 × 10 −2 至1.706 × 10 −2 0.1266 × 10 −2 至0.1738 × 10 −2 1.5% 2.353 × 10 −2 至2.868 × 10 −2 0.6388 × 10 −2 至0.6454 × 10 −2 表2. 混凝土受压中期局部最大应变值 混凝土受压后期，如图3~6中(d)~(g)所示，局部最大应变继续增加，由图中标尺可读出，混凝土受压后期局部最大应变见表3，由此可得，混凝土受压后期纳米含量为1%的混凝土较纳米含量为1.5%的混凝土局部最大应变变化更加明显，且与中期比较1%纳米混凝土应变剧增，从应变云图中可以看出是因为含量为1%的混凝土左下角破坏导致。混凝土在受压过程中除纵向压应变外，还有横向膨胀应变，所以结合横、纵向集中应变云图及三维应变云图(图7，图8)可更加明显地看出，纳米含量为1%的混凝土产生了一条近乎贯穿的横向裂缝，和试件左右两条近乎贯穿的纵向裂缝；而纳米含量为1.5%的混凝土仅有试件左右两侧出现近乎贯穿的裂缝。 Table 3 Local maximum strain value in the later stage of concrete compressio 纳米含量 横向局部最大应变 纵向局部最大应变 1% 5.356 × 10 − 2 至9.489 × 10 − 2 1.988 × 10 -2 至6.239 × 10 −2 1.5% 3.002 × 10 − 2 至5.427 × 10 − 2 0.6479 × 10 -2 至0.7704 × 10 −2 表3. 混凝土受压后期局部最大应变值 图7. 纳米掺量为1%的混凝土的横向应变云图 图8.纳米掺量为1.5%的混凝土的横向应变云图 试件破坏，如图3~6中(h)所示，在试件破坏时，各含量的纳米混凝土左右两侧均出现贯穿竖向裂缝，纳米含量为1%的混凝土出现横向贯穿裂缝，纳米含量为1.5%的混凝土也出现了很多较小的横向裂缝。说明混凝土破坏是由横向裂缝和纵向裂缝共同形成。
system: [SYSTEM]

user: 4. 结论 本实验基于数字散斑技术对纳米混凝土在单轴受压过程中的全场位移和应变进行了采集，通过对比不同掺量纳米混凝土荷载与位移曲线，分析了混凝土在不同纳米材料掺量、不同受力阶段下应变云图的变化特征和变形参数的变化规律。 1) 纳米含量为1.5%的混凝土于纳米含量为1%的混凝土相比，其强度降低了15.7%。其降低原因可能是由于随着纳米含量的增加，纳米二氧化硅在分散时出现团聚的可能性增加，从而导致混凝土强度降低。 2) 数字散斑相关方法可有效观测纳米混凝土试件在变形前后应变演化趋势，从而对试件裂纹的扩展进行计算。 3) 混凝土受压过程中，发生的应变突变现象是因为混凝土局部破坏导致。 4) 混凝土最终破坏模式是由横向裂缝和纵向裂缝共同形成。
system: [SYSTEM]

user: 摘 要 随着互联网的日益盛行，网上商店遍地开花，网上购物正潜移默化地影响我们的生活方式，融入到新生代生活的点点滴滴。为了更好地了解新生代消费者对网购的多样化需求与对网购的相关评价，展开了此次调查。文章在研究了新生代购物行为和新生代网购满意度后，得出了新生代对网购的总体态度，也发现了质量、信誉、服务和购物环境等问题，并给出了相应建议。 关键词 :新生代，网络购物，满意度 Copyright © 2020 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 研究背景 近年来，网络购物凭借着实惠的价格、便捷的模式成为新的消费趋势，并且网络购物呈现了年轻化的特点，新生代用户的规模不断增加。作为新生代，他们到底对于网络购物是持什么样的态度呢？ 为了更好地了解新生代消费者对网购的看法与需求，本次研究从新生代消费者对卖家服务、产品信息、物流、售后满意度等角度展开调查，了解该年龄段的网购满意度与购物体验、不同的产品需求以及对卖家服务的要求，探索其网购特点，对于存在的问题寻求解决办法，以提高新生代消费者网络购物满意度，给新生代带来更好的网络购物体验。
system: [SYSTEM]

user: 2. 研究过程与方法 (一) 研究的对象 通过对新生代消费者网络购物行为进行研究，通过对新生代物流满意度、客服满意度、商品性价比满意度等方面展开调查新生代消费者网络购物总体情况，并由此对当前网络购物存在的问题提出解决方法和措施。本次调查共回收144份有效问卷，包含90后65份，00后79份。 (二) 研究方法 1、文献分析 目前对于线上购物满意度的研究中，多数文献对研究对象的选择多为大学生或是某一特定职业群体，而选择的研究方向常从研究对象的特征、线上购物的核心优势、物流配送等角度进行研究 [ 1 ]，同时有少数学者如史晓丹从提供线上购物服务的网站质量和安全性方面进行研究 [ 2 ] [ 3 ]。 2、问卷调查 本次制定了两份问卷——90后与00后对线上购物的满意度调查问卷。调研的对象都为00后90后人群。本问卷的前四题为基本情况调查，了解新生代线上购物的月购买次数月消费和偏好等基本情况。五到八题为对线上购物各项指标满意度情况调查，为本调查的重点。九到十题为对是否愿意退货及给差评情况原因的调查，以此了解新生代面对网购不好结果时的心态。十二题为网购与实体购物满意度对比，从而分析线上于线下购物的发展情况。十四题与十五题为填空题，了解新生代选择线上购物的原因及对线上购物的建议。 3、访谈法 事先准备关于网络购物满意度的题项清单，按照事先设计好的问题来询问调查中的新生代消费者，了解该消费者对网络购物的真实态度，并且对网络购物中的各个方面进行询问，以期能够了解最真实的数据，提高网络用户购物体验。
system: [SYSTEM]

user: 4. 研究的结论与建议 (一) 结论 基于电商行业蓬勃发展以及新生代消费者独特的购物特征，本文对新生代消费者网购满意度进行具体研究。发放网络问卷的同时运用文献分析法、访谈法等研究方法，深入研究了新生代消费者对网购满意度的相关情况。具体从新生代消费者的特征、影响新生代消费者消费行为的因素以及新生代消费者对线上购物全程各主体的满意程度三大角度进行具体研究。从调查结果来看，新生代消费者对于线上购物整体的满意程度较高，从线上购物的全过程来看，无论是在线上购物前期对客服的满意程度或是中期对物流服务的满意程度还是后期对于售后服务的满意程度，新生代消费者都是持较为满意的态度，同时对于线上购物全程所遇见的问题，也都是持包容理解的态度。 对新生代消费者特征的研究。新生代消费者群体女性占比大，未来，不仅女性购物潜力巨大而且伴随网购行业的日益成熟，男性线上购物群体也会进一步扩大，同时整体上新生代消费者每月的购物次数有待进一步提高，并且由于新生代消费者思想的开放、对网购的接受程度较高以及熟悉网购的具体操作和流程等各方面的原因，网购行业的发展前景是十分美好的。 对影响新生代消费者消费行为因素的研究。网上购物方便快捷、价格低廉等核心优势对新生代消费者的影响程度较大，但是也存在商品质量差、平台商家失信和购物环境恶劣等问题，影响甚至阻碍新生代消费者线上购物的欲望，因此要想网购行业能进一步健康、良性发展，就需尽快对这些问题进行针对性的纠正。 最后，对新生代消费者对线上购物全程各主体的满意程度的研究。可见新生代消费者对于物流的时效性、客服的服务态度等各方面都会有较高的要求，这也迫使相应主体进行完善，促进网购行业可持续发展。 (二) 建议 鉴于研究结果与分析，虽然整体上网购行业的发展前景是十分美好的，但是也存在一些不可回避且阻碍网购行业发展的问题，因此针对几大突出问题，提出策略建议： 1、关于商品质量的问题 现阶段各网购平台都存在不同程度的商品质量差、以次充好等商品质量问题。因此，从国家角度出发，相关监察机构应加强对商品质量的检查力度，同时根据网购业发展特点创新检查方式，依法打击制售的不法行为，以提高产品质量，维护消费者和产品生产企业的合法权益。从网购服务平台角度出发，各平台首先应加强平台自身管理，建立严格的商品筛选机制。使商品质量的好坏与商家的等级高低挂钩，逐步减少平台内劣质商品数量以达到提高平台内商品质量整体提升的目的。从商家角度出发，商家自身首先要提高自身思想道德素质，从源头拒绝售卖假冒伪劣产品，同时，商家也要服从、配合平台的管理。平台和商家形成配合，共同促进电子商务健康有序发展。 2、关于购物环境的问题 存在用户评价、交易数据失信，即好评返利、差评骚扰、找人刷单等许多严重影响购物环境的不良行为。这些行为具有普遍性的特点，同时这些行为表面上看是会促进商家的发展，实际上却是会导致商家之间恶性竞争，增加商家的负担，同时也影响消费者评论的真实性和公正性，进而导致商家无法客观地了解消费者的满意程度，无法改善产品的质量和服务水平。鉴于这些问题，提出如下建议：社会层面：积极完善诚信环境，培育诚信文化，促进消费者自身诚信水平的提高，从而减少或避免用户评论失信的情况出现；网购平台层面：平台可以严格消费者身份认证，同时也可以通过建立第三方信用评价体系或者与政府相关部门合作共通信用体系进一步完善信用评价体系。同时平台也可以拓展提高技术水平，加大对平台不良竞争行为的监管力度；商家层面：商家自身应明确这些行为的不良影响，杜绝不良竞争行为，避免恶性竞争。 3、关于售后服务的问题 售后服务也属于购物的一部分，售后服务的好坏将影响平台服务水平，商家的服务水平甚至品牌的名誉。但是，平台内各种商品的售后服务水平参差不齐，同时也都存在售后客服不负责、不专业的情况。针对这些问题，我们给出如下建议：国家层面：建立或完善针对网购的统一售后服务制度；网购平台层面：规范平台内商家售后服务，同时也可提高售后服务在店铺评级中所占比重；商家层面：商家可对售后服务人员进行培训，同时有能力的商家也可以将线上线下的售后服务进行结合，尽可能的给消费者创造便利，实现真正的顾客至上。 4、关于线上店铺和实体店铺经营 第一、对于线上店铺，可以推出一些实景试用模式，例如现在许多美妆店铺推出的线上试妆等方式，提高消费者的体验感，从而帮助消费者可以更好地做出选择。第二、对于实体店来说，店家可以注册一些公众号，及时的上传一些店铺新款信息也可以通过公众号了解顾客的爱好，需求和建议。再者可以提高实体店设置的美观度，以迎合新生代消费者对美的需求，增强其好感度。
system: [SYSTEM]

user: 摘 要 随着对地下矿产的开采，采空区边坡地质灾害突发现象日益突出，尤其是崩塌和滑坡造成了重大生命财产损失，近年来引起了广泛的关注。本文以2017年8月28日贵州纳雍县张家湾镇普洒村崩塌滑坡为研究对象，通过现场勘查，对崩滑区的概况有一个初步掌握，深入分析其成因机制，并结合PFC数值模拟，对边坡的崩滑发展过程进行模拟。分析研究认为，该边坡在自然风化作用以及开采作用下，岩体节理裂隙开始发育，由于坡体中部陡缓交界处应力集中，促进坡体局部软岩崩滑，使坡顶岩体开始出现拉裂缝，随着开采的继续，裂隙逐渐贯通，并沿陡倾结构面发展，随着裂隙的发展，坡体顶部危险岩体以崩塌为主，中下部危险岩体以边破碎边滑动的方式沿坡面往下运动，整个边坡破坏由崩塌破坏为主变为滑移破坏为主。研究贵州纳雍县张家湾镇普洒村崩塌滑坡成因机制，对我国山区存在类似条件的灾害隐患点的研究以及未来的减灾防灾工作，具有重要的指导意义。 关键词 :采空区边坡，崩滑，成因机理，PFC 2D Copyright © 2020 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 从工程力学角度看，矿产资源原属于地质体的一部分，自然形成的地质体在长久的地质作用下，内部应力趋于平衡，然而工业建设的发展需要大量的矿产资源，自然资源的开发使得地应力重分布，岩体结构劣化，诱发了地面塌陷、山体崩塌、山体滑坡、地表移动等一系列地质灾害问题。近年来，地下开采引起的地质灾害，以及地表环境变化十分严重，主要表现为采空区塌陷和山体开裂，从而导致塌陷、滑坡等地质灾害，造成重大生命财产损失 [ 1 ]。如贵州省大方县德兴煤矿崩滑，贵州凯里市龙场镇崩塌等 [ 2 ]。 本文以贵州纳雍县张家湾镇普洒村崩塌滑坡为例，通过现场调查测绘，研究其变形和破坏情况。根据采空区滑坡的工程地质条件，采用数值模拟方法对地下开采诱发滑坡的形成机制进行了详细研究，并分析煤层地下开采过程中滑坡的采空区稳定性及应力应变变化。为我国山区存在类似条件灾害隐患点的研究以及未来的减灾防灾工作，具有重要的指导意义。
system: [SYSTEM]

user: 2. 概况 2017年8月28日10时30分左右，贵州省纳雍县张家湾镇普洒村山体发生高位崩塌地质灾害，造成重大生命财产损失，伤亡人员高达43人，引发广泛关注 [ 3 ]。 根据现场调查，崩塌山体位于煤矿采空区上部。山体崩塌失稳后，岩体迅速破碎、崩落、铲刮，并沿坡体高速滚落滑动，由于坡角处，村庄建筑物的阻挡，崩滑体停止运动。在此过程中，坡面原有的松散堆积体被铲刮并随崩滑体一起运动，形成了“滚雪球”效应。沿坍塌方向最长地层约610 m，平均长度约575 m，中部宽度360~380 m，平均厚度约4 m，体积约82.3 × 104 m 3 (见图1)，是典型的由地下矿区开采而引发的大型岩质滑坡 [ 4 ] [ 5 ]。 图1. 现场崩塌概况图
system: [SYSTEM]

user: 3. 成因机制分析 为进一步研究采空区边坡崩滑成因机理，此次将该滑坡成因机制分为五个阶段：1) 冻融循环、干湿循环造成岩体节理裂隙的形成；2) 坡体地形地貌变化差异大，砂泥岩互层，造成局部应力集中；3) 随着煤矿区的继续开采，采空区域不断扩大，岩体应力重分布，结构面不断劣化；4) 量变导致质变，当岩体中部裂隙发育到一定程度时，上部岩体在自重作用下挤压下部岩体并引起破坏；5) 局部崩塌带动整体失稳 [ 6 ] [ 7 ] [ 8 ]。具体分析如下： 1) 冻融循环、干湿循环等自然风化作用导致岩质边坡发生风化，加之矿区内爆破震动及开采等工程活动，促进坡体内外岩体节理裂隙的发育和贯通(图2)。 图2. 裂隙发育阶段 图3. 坡顶拉裂阶段 2) 坡体下部坡面较缓，而上部坡面较陡，在中部区域交汇，此处岩体由于坡面倾角变化较大，在振动荷载下容易形成应力集中。同时，坡面中部中厚层泥岩和泥质砂岩等软岩风化破碎速度较快，应力集中促进坡体局部软岩崩滑，使上部风化较慢的灰岩成悬空状态，坡顶岩体开始出现拉裂缝，同时为边坡岩体变形提供空间条件(图3)。 3) 随着煤矿区的继续开采，采空区域不断扩大，加剧边坡岩体应力重分布，结构面不断劣化，节理裂隙不断贯通，坡体后缘竖向裂隙向下发展。裂隙往临空面方向的不稳定岩体在重力作用下产生压缩下部岩体。由于顶部岩体背靠在后墙上，而趾部岩体没有坍塌，因而岩体中部由于不平衡力矩引起弯曲并逐渐断裂(图4)。 图4. 裂隙贯通阶段 4) 量变导致质变。当岩体中部裂隙发育到一定程度时，上部岩体在自重作用下挤压下部岩体，导致下部岩体不断破坏，结构面不断劣化，并引起破坏。此时，岩体下部的岩体首先会因失稳而下沉，最后会呈现出整体坐落的特征。当节理裂隙缓慢发育的岩体中应力达到其极限强度时，难以承受上部岩体的自重压力，发生突发性压碎，进一步，上部不稳定岩体沿控制性节理裂隙发生突发性局部大规模崩滑(图5)。 图5. 裂隙扩展、局部崩滑阶段 5) 崩滑开始时，坡体顶部危险岩体以崩塌为主，由于岩体节理裂隙十分发育，在崩塌的过程中，中下部危险岩体并没有发生整体外翻，而是以边破碎边滑动的方式沿坡面往下运动，主要堆积于坡脚附近，即由崩塌破坏为主变为滑移破坏为主(图6)。 图6. 整体崩滑阶段 根据崩滑后的地形反推坡面的安全角度约为30˚±5˚，破碎岩体边坡综合内摩擦角约为30˚±10˚ (取30˚作为平均值)。崩滑前后地形线存在一个交点，交点以上的岩体均是不稳定的，在不平衡力的作用下，节理裂隙逐渐发育和局部贯通。量变产生质变，控制性节理裂隙贯通后，不稳定岩体突然启动崩滑模式。
system: [SYSTEM]

user: 4. 数值模拟 基于上述分析，选取主剖面(图7)构建计算模型，采用二维离散元软件PFC 2D (5.0)对边坡的崩滑发展过程进行模拟 [ 9 ] [ 10 ]。 图7. 工程地质剖面图 计算模型采用崩滑前的地形线，为了减少计算量，对原始剖面进行简化是有必要和可行的，因为高程在83.22 m以下的坡体对整个坡体的崩滑过程没有实质性影响，因此建模时可以忽略这部分坡体，顶点与底部的相对高度差232.33 m，计算模型长度464.89 m。 计算模型由10，949个半径为1 m的球充填而成，球之间呈正三角形排列，采用平行粘结接触模型将球体粘在一起，粘结后的颗粒集合体具有硬岩(灰岩、白云岩)的强度和弹性模量，球体和接触模型的计算参数见表1。 这种规则排列构成的虚拟模型与实际坡体的主要区别在于：1) 实际岩体岩性不同，而虚拟岩体采用统一细观参数，岩性相同；2) 组成岩体的颗粒组成分不同，而虚拟岩体由统一粒径的颗粒构成。上述两方面的区别可归结为虚拟岩体试样没有考虑实际岩体各向异性的特性。因此，各向异性对岩体节理裂隙发育的影响无法通过这种模型直接获取。然而，根据第二节的成因机制分析，可以通过控制节理裂隙的发育，模拟坡体崩塌滑坡启动和运动过程。 Table 1 Calculation parameters of parallel bond mode 模型 Bond Gap Emod (linear) Kratio (linear) Emod (pb) Kratio (pb) Pb_ten Pb_coh Pb_fa Dp_nratio Fric 平行粘结 0.5e-4 1.0e9 1.0 1.0e9 1.0 10e6 50e6 0.0 0.5 0.577 表1. 平行粘结模型计算参数 统计岩体中存在的节理裂隙是个比较困难和危险的工作，然而在PFC中，可以采用Discrete Fracture Network (DFN)功能模块，随机的在坡体内生成节理裂隙。节理裂隙长度为1~3倍的颗粒直径，倾斜角度服从高斯分布，平均值80˚，标准差10˚。 节理裂隙的数量根据截面面积确定。首先，现场收集试样，最小尺寸大于0.3 m，然后将试样沿中间切成2半，取一半，磨平切面，在切面上用记号笔画一个半径为0.1 m的标准圆，采用放大镜，寻找出切面上的节理裂隙，1 cm以下的裂隙不考虑，用记号笔做好标记，然后统计标准圆内节理裂隙的数量N。重复随机取样3次，统计3次节理裂隙数量并求取平均值，最终计算单位面积内节理裂隙的数量N/A (A为半径0.1 m标准圆的面积)约为0.5%。根据虚拟模型的面积和统计的单位面积内节理裂隙的数量得到虚拟模型中节理数量为164条。 节理裂隙采用光滑节理接触模型，计算参数见表2。默认的接触模型采用Hertz接触模型，相关参数见表3。如图8所示，在自重应力下，坡体上的颗粒发生的位移很小，处于稳定状态。 Table 2 Calculation parameters of smooth joint mode 模型 Kn/N/m Ks/N/m Fric Da/(˚) Ten/Pa Shear/Pa Cohe/Pa Fa Large 光滑节理 2e9 2e9 0.577 0.0 0.0 0.0 0.0 0.0 1 表2. 光滑节理模型计算参数 Table 3 Hertz contact model calculation parameters and sphere parameter 类别 Hertz_shear Hertz_poiss Fric Damp Dp_nratio Hertz接触 2.45e10 0.3 0.35 - 0.5 球 2.45e10 0.3 0.35 0.0 - 表3. Hertz接触模型计算参数及球体参数 根据纳雍崩滑资料，坡体中上部的软岩先被压碎，紧接着，上部硬岩发生崩滑。说明，该处软岩是上部硬岩保持稳定的一个支点，推测硬岩内部发育有基本贯通的控制性裂隙，使硬岩脱离母岩，成为危险性岩体。因此模拟坡顶崩滑的启动过程时，在软岩处，植入不断发育的控制性裂隙4条，如图9所示，倾角分别为0˚，120˚，135˚和−45˚，点A是它们的发育始点。每条裂隙按照既定的方向以1 m/s的速度增长，即每秒增长1 m，且各自在增长方向上存在偏差，均为±10˚。随着裂隙的发育，A点下部出现局部岩体崩滑，使A点上部岩体成悬空状态。 裂隙的发育导致，带来局部岩体崩滑的同时，也导致岩体内应力不均衡，在裂隙两端或坡面转角处发生岩体破碎，触发新的局部崩滑，如图10所示。 随着软岩的破碎和崩滑，失去支点的顶部硬岩最终启动崩滑模式，如图11，硬岩以破碎的方式沿坡面滑动，同时，后缘裂隙开口变大。 启动崩滑后的岩体沿坡面运动，大块岩体不断破碎成小块岩体，最终堆积成新的坡面，如图12所示。 图8. 纳雍边坡计算模型(自重下的位移量) 图9. 局部软岩崩滑 图10. 岩体局部崩滑 图11. 硬岩启动崩滑模式 图12. 崩滑后新的坡面
system: [SYSTEM]

user: 5. 结语 1) 通过现场勘查，研究分析了贵州纳雍县张家湾镇普洒村崩塌滑坡的特征及成因机理。 2) 在掌握该采空区边坡崩塌的特征及成因机理的基础上，运用二维离散元软件PFC 2D (5.0)对边坡进行数值建模，对崩滑发展过程进行模拟，基本再现了采空区边坡崩滑过程。 3) 该边坡在自然风化作用以及开采作用下，岩体节理裂隙开始发育，由于坡体中部陡缓交界处应力集中，中厚层泥岩和泥质砂岩等软岩风化破碎速度较快，促进坡体局部软岩崩滑，使上部风化较慢的灰岩成悬空状态，坡顶岩体开始出现拉裂缝，随着开采的继续，裂隙逐渐贯通，坡顶拉裂缝向下发展，坡体顶部危险岩体以崩塌为主，下部危险岩体以边破碎边滑动的方式沿坡面往下运动，整个边坡破坏由崩塌破坏为主变为滑移破坏为主。整个破坏过程可以概括为“坡体中部关键性块体失稳–坡顶拉裂缝发育–坡体下部岩体压碎”。
system: [SYSTEM]

user: 摘 要 医药冷链配送模式需求量日益增加，北京作为国家首都，人口密集，药品的需求量远高于其他城市。药品的冷链配送对于物流行业来说一直备受关注，面对高需求量、高标准的药品冷链配送，我们要做到严格要求，确保全程冷链。本文通过蚁群算法以及MATLAB仿真软件，对A医药集团在北京地区的配送路径进行优化，并对优化结果从环境、企业、患者三个角度进行分析。 关键词 :医药物流，冷链配送，路径优化，蚁群算法 Copyright © 2020 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 随着我国物流业的发展，物流配送水平正在不断提高，冷链业务已渗透到各行各业，药品安全问题关乎民生，超温的情况下会使药品失效，患者使用失去药效的药品会拖延病情，增加治疗难度，降低治愈率，不断弥补药品冷链配送行业的不足，对于现在的医药冷链物流业的发展非常重要 [ 1 ]。为了使药品在配送过程中的药效得到保障，全程冷链配送模式逐渐走进医药圈。全程冷链成本非常高，对配送路径进行优化是最直观、有效的降低成本方式，不断缩短配送时间，为企业带来最大收益，为患者争取治疗时间 [ 2 ]。 本文基于A医药集团在北京地区某一天的冷链配送情况，选取28家固定配送对象，并针对原有的配送路径进行合理的路径优化。
system: [SYSTEM]

user: 2. 研究过程 2.1. 研究步骤 1) 阅读配送路径优化方面的文献，学习路径优化算法，并记录学习结果。 2) 对A医药集团在北京地区医药冷链物流情况进行简单介绍，并发现配送路径存在的问题。 3) 针对于A医药集团在北京地区的配送情况，采用蚁群算法对配送路径进行优化，并利用MATLAB进行仿真实验并分析结果。 2.2. 研究方法 1) 文献学习法。通过阅读大量关于医药冷链配送、路径优化方面的文献，对该领域有了充分的了解，并将二者结合进行分析。 2) 数据统计法。利用现有的公开数据进行分析，了解企现状。对具体配送对象和配送中心的经纬度坐标进行搜集统计，为仿真实验做准备。 3) 建模法。利用蚁群算法对已有数据进行建模分析，通过MATLAB工具软件进行算法实现。