截面温度不均匀对受约束波纹腹板梁火灾下悬链线效应的影响.pdf

返回 相似 举报
截面温度不均匀对受约束波纹腹板梁火灾下悬链线效应的影响.pdf_第1页
第1页 / 共7页
截面温度不均匀对受约束波纹腹板梁火灾下悬链线效应的影响.pdf_第2页
第2页 / 共7页
截面温度不均匀对受约束波纹腹板梁火灾下悬链线效应的影响.pdf_第3页
第3页 / 共7页
亲,该文档总共7页,到这儿已超出免费预览范围,如果喜欢就下载吧!
资源描述:
第33卷 第6期 2016年11月 建筑科学与工程学报 Journal of Architecture and Civil Engineering Vol. 33 No. 6 Nov.2016 文章编号1673 - 2049201606 - 0053 - 07 收稿日期2016 - 05 - 13 基金项目 国家自然科学基金项目51578322 ; 山东省自然科学基金项目ZR2015EM041 ; 山东省重点研发计划项目2016GGX103013 作者简介 王培军1976 - , 男, 山东菏泽人, 副教授, 工学博士,E - mail pjwangsdu. edu. cn。 截面温度不均匀对受约束波纹腹板梁火灾下 悬链线效应的影响 王培军1, 刘昌斌1, 杨振民2 1.山东大学 土建与水利学院, 山东 济南 250061;2.中通钢构股份有限公司, 山东 聊城 252000 摘要 以受约束波纹腹板梁为研究对象, 利用有限元模型进行了截面温度均匀和不均匀分布对轴向 约束波纹腹板梁悬链线效应影响的对比分析。结果表明 由于不均匀温度下产生的热弯曲效应, 约 束钢梁在受压阶段的竖向挠度更大而梁内轴力更小; 均匀温度受约束钢梁的临界温度比不均匀温 度钢梁高6. 2; 不均匀温度钢梁的温度比越小, 钢梁的热弯曲效应越明显, 相应的临界温度越低。 关键词 波纹腹板梁; 轴向约束钢梁; 火灾; 不均匀温度分布; 悬链线效应 中图分类号TU375. 1 文献标志码 A Effects of Non - uni Temperature Distribution Across Section on Catenary Action of Restrained Corrugated Web Beam in Fire WANG Pei - jun1,LIU Chang - bin1,YANG Zhen - min2 1.School of Civil Engineering,Shandong University,Jinan 250061,Shandong,China; 2.Zhongtong Steel Structure,Co.,Ltd,Liaocheng 252000,Shandong,China AbstractThe restrained corrugated web beam was taken as research object,and the influences of uni and non - uni temperature distribution across section on catenary actions of axially restrained corrugated web beam were analyzed by using the finite element model.The results show that due to the thermal bending effect produced by the non - uni temperature,the vertical deflection of the restrained steel beam is larger at the compression stage,and the axial force in beam is smaller.The failure temperature of the steel beam with uni temperature is 6. 2 higher than that of the steel beam with non - uni temperature.The smaller the temperature ratio is,the more obvious the thermal bending effect of steel beam is,and the lower the failure temperature is. Key wordscorrugated web beam;axially restrained steel beam;fire;non - uni temperature distribution;catenary action 0 引 言 波纹腹板梁广泛应用于大跨结构和桥梁中, 腹 板的波折提高了其平面外刚度和抗剪屈曲承载能 力, 并避免设置横向加劲肋。李立峰等[ 1]通过 4根 波形钢腹板H型钢梁的屈曲加载试验, 对弹性屈曲 强度的计算方法进行了分析, 给出了弹性屈曲强度 简化计算公式和适用范围。张哲等[ 2]通过对 4根试 万方数据 件进行抗弯试验, 研究了波纹腹板H型钢组合梁的 受弯性能。 实际结构中钢梁受到周围梁柱的约束, 火灾下 大变形阶段时, 钢梁轴向拉力与竖向挠度形成的悬 链线效应可以协助钢梁继续承载, 提高钢梁的耐火 极限。Liu等[ 3 - 4]采用试验和有限元方法参数化研 究了火灾下钢梁的悬链线效应。Wu等[ 5]给出了高 温下钢筋混凝土梁悬链线效应的分析计算程序。波 纹腹板梁轴向刚度小于相应的平腹板梁, 其悬链线 效应 与 平 腹 板 梁 有 较 大 差 别。张 哲 等[ 6 - 7]采 用 ABAQUS有限元模型对火灾下轴向约束波纹腹板 梁的悬链线效应进行了参数分析, 并且对比研究了 波纹腹板梁和平腹板梁在悬链线效应方面的异同。 实际结构中的钢梁上翼缘受到混凝土楼板的保 护, 处于三面受火状态, 梁截面温度分布不均匀, 上 翼缘的温度比下翼缘和腹板的低。王培军等[ 8]通过 对9根受约束钢柱的抗火试验发现截面温度不均匀 导致钢柱在绕截面对称轴弯曲时同时发生扭转。 Dwaikat等 [9]通过对三面受火钢梁的有限元分析, 给出了考虑不均匀温度影响时P - M曲线 其中,P 为钢梁的轴向承载力,M为钢梁的抗弯承载力 的 简化修正方法。梁截面温度的不均匀分布导致钢梁 悬链线效应与均匀温度分布时有较大差别。目前尚 无对火灾下截面温度不均匀分布时轴向约束波纹腹 板梁悬链线效应的研究。 本文采用ABAQUS有限元软件对均匀温度和 不均匀温度分布下轴向约束波纹腹板梁悬链线效应 阶段的轴力、 挠度和悬链线弯矩的发展进行对比分 析, 研究不均匀温度分布中高度比和温度比对钢梁 悬链线效应的影响。 1 波纹腹板梁及不均匀温度分布 1. 1 波纹腹板梁参数 本文拟研究的波纹腹板梁的腹板形状采用 波 纹腹板钢结构技术规程 [10]推荐的波形, 其尺寸如 图1所示。梁的跨度为9 600 mm, 沿跨度方向上有 40个完整的波形。 1. 2 不均匀温度分布参数 钢梁截面温度简化如图2所示 其中, κ,ω分别 为不均匀温度分布中高度比和温度比 。以温度比 为0. 5, 高度比为1/4为例, 上翼缘的温度是下翼缘 温度的50, 温度的变化范围为腹板高度的1/4。 Yin等 [11]在分析火灾下钢梁的弯扭屈曲时采用了 同样的温度分布形式。Liu等[ 3]的火灾试验结果也 图1波纹腹板梁尺寸 单位mm Fig. 1Dimensions of Corrugated Web BeamUnitmm 图2梁截面均匀温度和不均匀温度分布 Fig. 2Uni and Non - uni Temperature Distribution Across Section of Beam 验证了这种温度分布形式。本文不均匀温度分布下 钢梁的温度均指下翼缘的温度。 2 有限元模型及模型验证 2. 1 有限元模型 轴向约束波纹腹板梁关于平面 yOz 对称, 有限 元模型中仅模拟半跨钢梁, 如图3所示。钢材的泊 松比和热膨胀系数分别为0. 3,1. 410-5℃-1。常 温下 钢 材 的 屈 服 强 度 为345 MPa, 弹 性 模 量 为 2. 05105MPa。火灾高温对钢材的强度和弹性模 量有显著影响, 随着温度升高, 钢材的弹性模量和屈 服强度均逐渐降低。本文高温下材料模型采用欧洲 规范[ 12]的规定。单元类型为 S4R壳单元, 通过网格 敏感性分析确定单元宽度为20 mm。 在钢梁端部建立MPC约束ABAQUS中的多 节点约束 , 并限制MPC控制点沿y轴和z轴的位 移及绕z轴的转动。采用线弹性弹簧Spring2单元 连接MPC控制点, 模拟梁端弹性轴向约束。梁跨 中采用对称的边界条件。为防止钢梁发生弯扭屈 曲, 限制梁上翼缘的平面外位移。实际工程中, 钢梁 上翼缘设置有混凝土楼板, 可有效限制钢梁上翼缘 的侧移。有限元分析过程采用恒载升温方式, 第1 步在梁上翼缘施加均布荷载, 第2步使整个模型升 温, 梁截面的不均匀温度分布采用ABAQUS中的 45 建筑科学与工程学报 2016年 万方数据 图3波纹腹板梁有限元模型 Fig. 3Finite Element Model of Corrugated Web Beam 自定义函数实现。为避免火灾下钢梁大变形分析时 的不收敛问题, 引入人工阻尼系数α, 并取能量耗散 系数为尽可能小值[ 4]。 2. 2 模型验证 采用Liu等[ 3]轴向约束刚度为 62 kNmm-1 时火灾下平腹板梁的试验结果来验证有限元模型的 准确性。荷载比 常温下简支梁最大弯矩与塑性抗 弯承载力的比值 分别为0. 5,0. 7。有限元模型与 试验结果对比如图4所示, 其中正值代表拉力。有 限元分析采用了ABAQUS默认的收敛准则, 即迭 代的残差小于相应荷载或位移增量的0. 5。由图 4可以看出, 本文有限元分析结果与Liu等 [3]的试 验结果以及Yin等[ 4]的有限元模拟结果吻合较好, 不同人工阻尼系数对分析结果影响较小, 下文分析 采用α110-10。验证结果说明了本文有限元模 型的适用性。 3 钢梁悬链线效应对比分析 在钢梁上翼缘施加荷载比为0. 5的均布荷载, 大小为143. 6 kPa。梁 端 的 轴 向 约 束 刚 度KC 0. 1KB, 其中KB为钢梁常温下的轴向刚度, 则KC 的大小为17 083 kNm-1。 3. 1 轴力对比分析 图5为2种温度分布下波纹腹板梁的轴力变化 规律。在受压阶段, 由于不均匀温度分布引起的热 图4有限元模型验证 Fig. 4Verification of Finite Element Model 弯曲提前释放了钢梁的膨胀变形, 且不均匀温度分 布时钢梁截面的平均温度较低, 不均匀温度分布下 约束钢梁轴力随温度的增长速度比均匀温度分布时 的低, 钢梁在更高温度下达到了较低的最大轴向压 力。在悬链线阶段 轴力为拉力时 , 不均匀温度分 布下约束钢梁的上翼缘温度较低, 上翼缘的轴向拉 力较大, 钢梁整体的最大轴向拉力比均匀温度分布 时的高约10。高温下均匀温度分布时钢梁的轴 55 第6期 王培军, 等 截面温度不均匀对受约束波纹腹板梁火灾下悬链线效应的影响 万方数据 图5约束钢梁的轴力变化 Fig. 5Changes of Axial Force in Restrained Steel Beam 向拉力达到其抗拉屈服强度, 并随着温度的升高而 降低。不均匀温度分布下钢梁的上翼缘温度较低, 且在悬链线效应阶段其拉应力远低于屈服强度。 悬链线温度为Tc, 该温度下钢梁的轴力为0, 标 志着悬链线效应的开始,Tc与无约束钢梁的临界温 度相同; 临界温度为Tf, 该温度下钢梁失去抗弯平 衡能力, 是钢梁考虑悬链线效应时的临界温度。2 个临界温度之间的温差ΔT反映了考虑悬链线效应 后钢梁临界温度的提高值。因此, 温差ΔT反映了 悬链线作用。 约束钢梁在2种温度分布下悬链线温度差距较 小, 如表1所示。不均匀温度分布下钢梁悬链线温 度为608 ℃, 均匀温度分布时为591 ℃, 不考虑悬链 线效应时不均匀温度分布下钢梁的临界温度略高, 这是因为不均匀温度分布下钢梁的轴力较小, 抗弯 承载力较大, 其受力变形的发展速度略小于均匀温 度分布时。均匀温度分布时钢梁挠度的变化如图6 所示。均匀温度分布时约束钢梁高温下的悬链线效 应显著, ΔT比不均匀温度分布时高66 ℃。考虑悬 链线效应时均匀温度分布下钢梁的临界温度比不均 匀温度分布时高约6. 249 ℃ 。 表1约束钢梁的临界温度 Tab. 1Critical Temperatures of Restrained Steel Beam 温度分布形式 TcTfΔT 不均匀温度分布/℃ 608865257 均匀温度分布/℃ 591914323 3. 2 挠度对比分析 2种温度分布下波纹腹板梁的挠度变化如图6 所示。由图6可知 温度低于460 ℃时, 由于不均匀 温度分布下的热弯曲效应, 钢梁的挠度比均匀温度 分布时的大; 温度高于750 ℃时, 由于均匀温度分布 下钢梁轴向拉力减小, 钢梁的挠度迅速增加, 约束钢 梁破坏时均匀温度下钢梁的挠度明显大于不均匀温 度分布时的挠度。 3. 3 悬链线弯矩对比分析 约束钢梁的悬链线弯矩Mc等于钢梁梁端轴力 图6约束钢梁的挠度变化 Fig. 6Changes of Deflection in Restrained Steel Beam Ft和跨中挠度δ的乘积, 反映了悬链线效应对钢梁 抗弯承载力的贡献。 在受压阶段, 相比于均匀温度分布钢梁, 不均匀 温度分布约束钢梁在更高温度下达到其最大负悬链 线弯矩, 且弯矩值比均匀温度分布时的大34. 6 27 kNm , 如图7所示。进入悬链线阶段后,2 种温度分布下钢梁的悬链线弯矩几乎重合。悬链线 阶段后期, 均匀温度分布下钢梁的悬链线弯矩等于 外荷载产生的弯矩, 独自承担外荷载。不均匀温度 分布时, 高温下上翼缘的轴向拉力比下翼缘大, 导致 钢梁截面内产生负弯矩, 悬链线弯矩必须大于外荷 载弯矩才能满足平衡方程。 图7约束钢梁的悬链线弯矩变化 Fig. 7Changes of Catenary Bending Moment of Restrained Steel Beam 3. 4 轴力对比分析 2种温度分布下波纹腹板梁上下翼缘及腹板上 的轴力变化如图8所示。2种温度分布中下翼缘温 度是相同的, 高温时下翼缘的轴向拉力都达到其抗 拉屈服强度, 并随着温度的升高而降低。上翼缘则 不同, 均匀温度分布下上翼缘的轴向压力在高温下 达到其抗压屈服强度, 并随着温度升高而降低。不 均匀温度分布下上翼缘的轴向压力远没有达到其屈 服强度, 由于梁截面受力平衡, 该轴向压力随着相应 下翼缘轴向拉力的降低而降低。 悬链线阶段后期, 均匀温度分布下钢梁上下翼 缘都达到了抗拉屈服强度, 钢梁变成真正的悬索构 65 建筑科学与工程学报 2016年 万方数据 图8约束钢梁的轴力变化 Fig. 8Changes of Axial Force of Restrained Steel Beam 件。不均匀温度分布下钢梁上翼缘的轴向拉力比下 翼缘大, 这也是钢梁内负弯矩产生的原因。无论温 度分布如何, 波纹腹板上产生的轴力远小于翼缘上 的轴力, 其对轴力的贡献可以忽略不计。 3. 5 弯矩贡献对比分析 火灾高温下受约束钢梁的弯矩平衡如图9所示 其中, q为梁上翼缘作用的均布荷载,L为梁跨度, Fy为梁端支座的竖向支反力 。高温下钢材软化使 钢梁的挠度迅速增加, 弯曲变形释放了热膨胀变形, 钢梁轴力由压力变成拉力, 成为拉弯构件。钢梁轴 向拉力Ft与跨中挠度δ形成悬链线弯矩Mc, 与跨 中截面弯矩Mu共同承担外荷载产生的弯矩Mp。 图9火灾下约束钢梁的弯矩平衡 Fig. 9Bending Moment Equilibrium of Restrained Steel Beam in Fire 火灾高温下受约束钢梁的弯矩平衡方程为 MuFtδMp 1 McFtδ 2 Mp 1 8 qL 2 3 2种温度分布下波纹腹板梁上下翼缘及腹板上 的弯矩变化如图10所示。2种温度分布下由波纹 腹板提供的弯矩远小于翼缘上的弯矩, 腹板作用可 图10约束钢梁的弯矩变化 Fig. 10Changes of Bending Moment of Restrained Steel Beam 以忽略不计。悬链线弯矩Mc和整体钢梁内部的弯 矩曲线对称, 两弯矩之和等于外荷载产生的弯矩 Mp, 满足约束钢梁的弯矩平衡公式1 。在整个升 温过程中, 外荷载产生的弯矩不变, 因此压力阶段的 负悬链线弯矩增大了钢梁的弯矩荷载; 悬链线阶段 的正悬链线弯矩协助钢梁分担弯矩荷载; 后期钢梁 内部产生负弯矩, 导致悬链线弯矩大于外荷载产生 的弯矩。 4 不均匀温度分布影响分析 4. 1 高度比的影响 约束钢梁的参数设置如表2所示。研究6组高 度比对约束钢梁悬链线效应的影响, 高度比分别为 0,1/8,1/4,1/2,3/4,1, 其中高度比为0代表梁截面 只有上翼缘的温度低, 高度比为1代表不均匀温度 沿整个梁截面线性分布。不均匀温度分布的温度比 为0. 5, 分析结果如图11所示。 高度比不同仅导致钢梁腹板的温度分布变化, 梁截面上下翼缘温度不变。由于火灾下波纹腹板对 钢梁悬链线效应的影响可以忽略, 不同高度比时波 纹腹板梁的各种反应几乎完全重合。 4. 2 温度比的影响 分别研究温度比为0. 1,0. 3,0. 5,0. 7,0. 9,1. 0 时钢梁悬链线效应的变化, 其中温度比为1. 0代表 沿梁截面温度均匀分布。 温度比越小, 上下翼缘的 表2约束钢梁的参数设置 Tab. 2Parameters Setting of Restrained Steel Beam κω 荷载比轴向约束刚度KC跨度/mm 0,1/8,1/4,1/2,3/4,10. 50. 50. 1KB9 600 1/40. 1,0. 3,0. 5,0. 7,0. 9,1. 00. 50. 1KB9 600 75 第6期 王培军, 等 截面温度不均匀对受约束波纹腹板梁火灾下悬链线效应的影响 万方数据 图11高度比的影响 Fig. 11Influences of Height Ratio 温度差越大。不均匀温度分布的高度比为1/4, 分 析结果如图12所示, 钢梁的临界温度见表3。 温度比越小, 上下翼缘的温差越大, 约束钢梁在 升温阶段的热弯曲效应越明显, 压力阶段钢梁的挠 度越大, 轴向压力越小。压力阶段约束钢梁的最大 悬链线弯矩随着温度比的增大先增大后减小。温度 比越小, 上翼缘的温度越低, 高温下钢梁的轴向拉力 越大, 主要依靠悬链线弯矩抵抗外荷载, 依据钢梁的 抗弯平衡公式1 , 当需求的悬链线弯矩相同时, 轴 向拉力越大, 维持弯矩平衡所需的挠度越小。值得 注意的是, 温度比为0. 9和1. 0 均匀温度分布 仅 相差0. 1, 但温度比为1. 0时钢梁在高温下的挠度 已显著增大。 由表3可知, 对于不均匀温度分布, 即温度比小 于1. 0时, 悬链线温度接近定值, 几乎不受温度比的 影响, 但约束钢梁的ΔT随着温度比的增大而增大, 温度比较大时钢梁上下翼缘温差较小, 悬链线效应 发展更充分。考虑悬链线效应时的临界温度随着温 度比的增大而增大, 当温度比从0. 1增大到0. 9时, 钢梁的临界温度从841 ℃升高到888 ℃, 温度上升 了约5. 647 ℃ 。 5 结语 1 不均匀温度分布引起的热弯曲释放了钢梁 的膨胀变形, 受压阶段约束钢梁的挠度更大, 轴向压 力更小。悬链线阶段不均匀温度分布时钢梁的最大 图12温度比的影响 Fig. 12Influences of Temperature Ratio 表3不同温度比时约束钢梁的临界温度 Tab. 3Critical Temperatures of Restraint Steel Beam with Different Temperature Ratios 温度比 0. 10. 30. 50. 70. 91. 0 Tc/℃602610608607605591 Tf/℃841854865870888914 ΔT/℃239244257263283323 轴向拉力比均匀温度分布时的高约10。高温下 梁截面的不均匀温度分布使梁内产生负弯矩。由于 均匀温度分布时钢梁挠度的充分发展, 其ΔT比不 均匀温度分布时高66 ℃。均匀温度分布时钢梁临 界温度比不均匀温度分布时高约6. 249 ℃ 。 2 受约束钢梁达到临界温度时, 截面温度均匀 分布的受约束波纹腹板钢梁上下翼缘均达到了受拉 屈服, 但对于截面温度不均匀分布的受约束波纹腹 板钢梁, 下翼缘达到了受拉屈服。由于上翼缘温度 较低, 其拉应力大于下翼缘拉应力, 且拉应力小于钢 材在对应温度下的屈服强度。 3 不同温度分布形式中高度比对波纹腹板梁 85 建筑科学与工程学报 2016年 万方数据 的悬链线效应几乎没有影响。温度比越小, 上下翼 缘的温差越大, 约束钢梁的热弯曲效应越明显, 压力 阶段钢梁的挠度越大, 轴向压力越小。钢梁临界温 度随着温度比的增大而增大, 当温度比从0. 1增大 到0. 9时, 钢梁的临界温度从841 ℃升高到888 ℃, 温度上升了约5. 647 ℃ 。 参考文献 References [1] 李立峰, 侯立超, 孙君翠.波形钢腹板抗剪性能的研究 [J].湖南大学学报 自然科学版,2015,4211 56 - 63. LI Li - feng,HOU Li - chao,SUN Jun - cui. Research on Shear Mechanical Property of Corrugated Steel Webs [J]. Journal of Hunan UniversityNatural Sciences, 2015,4211 56 - 63. [2] 张 哲, 李国强, 孙飞飞.波纹腹板H型钢组合梁抗 弯性能[J].土木工程学报,2014, 475 34 - 39. ZHANG Zhe,LI Guo - qiang,SUN Fei - fei.Flexural Behavior of the Composite Beam with Corrugated Webs[J]. China Civil Engineering Journal,2014,47 5 34 - 39. [3]LIU T C H,FAHAD M K,DAVIES J M. Experimen - tal Investigation of Behaviour of Axially Restrained Steel Beams in Fire[J]. Journal of Constructional Steel Research,2002,589 1211 - 1230. [4]YIN Y Z,WANG Y C. A Numerical Study of Large Deflection Behavior of Restrained Steel Beams at Ele - vated Temperatures[J].Journal of Constructional Steel Research,2004,607 1029 - 1047. [5]WU B,LU J Z. A Numerical Study of the Behaviour of Restrained RC Beams at Elevated Temperatures [J]. Fire Safety Journal,2009,444 522 - 531. [6] 张 哲, 刘昌斌, 王培军.火灾下轴向约束波纹腹板梁 悬链线效应分析[J].建筑钢结构进展,2015,175 58 - 64. ZHANG Zhe,LIU Chang - bin,WANG Pei - jun. Cate - nary Action Analysis of the Axially Restrained Corru - gated Web Steel Beams in a Fire[J]. Progress in Steel Building Structures,2015,175 58 - 64. [7] 张 哲, 刘昌斌, 王培军.火灾高温下轴向受约束波纹 腹板梁受力性能参数分析[J].钢结构,2016,311 94 - 99. ZHANG Zhe,LIU Chang - bin,WANG Pei - jun. Pa - rameter Analysis of the Mechanics Perance of Axially Restrained Corrugated Web Beams at Elevat - ed Temperatures in a Fire[J]. Steel Construction, 2016,311 94 - 99. [8] 王培军, 王 晓.截面温度不均匀受约束钢柱抗火性 能试验研究与有限元分析[J].建筑结构学报,2013, 343 28 - 34. WANG Pei - jun,WANG Xiao. Experimental and Nu - merical Studies on Restrained Steel Columns with Non - uni Temperature Distribution Across Sec - tion in Fire[J]. Journal of Building Structures,2013, 343 28 - 34. [9]DWAIKAT M,KODUR V. A Simplified Approach for uating Plastic Axial and Moment Capacity Curves for Beam - columns with Non - uni Thermal Gradients[J]. Engineering Structures,2010,325 1423 - 1436. [10]CECS 2912011, 波纹腹板钢结构技术规程[S]. CECS 2912011,TechnicalSpecificationforSteel Structures with Corrugated Webs[S]. [11]YIN Y Z,WANG Y C. Numerical Simulations of the Effects of Non - uni Temperature Distributions on Lateral TorsionalBucklingResistance of SteelI - beams[J]. Journal of Constructional Steel Research, 2003,598 1009 - 1033. [12]EN 1993 - 1 - 22005,Eurocode 3Design of Steel Struc - turesPart 1. 2General RulesStructural Fire Design[S]. 95 第6期 王培军, 等 截面温度不均匀对受约束波纹腹板梁火灾下悬链线效应的影响 万方数据
展开阅读全文

资源标签

最新标签

长按识别或保存二维码,关注学链未来公众号

copyright@ 2019-2020“矿业文库”网

矿业文库合伙人QQ群 30735420