TC18 钛合金双圆锥台试样热变形有限元模拟研究-sup-①-_sup-_毛颖.pdf

返回 相似 举报
TC18 钛合金双圆锥台试样热变形有限元模拟研究-sup-①-_sup-_毛颖.pdf_第1页
第1页 / 共6页
TC18 钛合金双圆锥台试样热变形有限元模拟研究-sup-①-_sup-_毛颖.pdf_第2页
第2页 / 共6页
TC18 钛合金双圆锥台试样热变形有限元模拟研究-sup-①-_sup-_毛颖.pdf_第3页
第3页 / 共6页
TC18 钛合金双圆锥台试样热变形有限元模拟研究-sup-①-_sup-_毛颖.pdf_第4页
第4页 / 共6页
TC18 钛合金双圆锥台试样热变形有限元模拟研究-sup-①-_sup-_毛颖.pdf_第5页
第5页 / 共6页
点击查看更多>>
资源描述:
TC18 钛合金双圆锥台试样热变形有限元模拟研究 ① 毛 颖1, 雷 雨2, 丁 旭1, 张晓泳3, 汪冰峰1,3, 樊 凯4, 冯抗屯5, 谢 静6, 王海鹏7, 雷家峰8 (1.中南大学 材料科学与工程学院,湖南 长沙 410083; 2.中铝洛阳铜业有限公司,河南 洛阳 471039; 3.中南大学 粉末冶金国家重点实验室,湖南 长 沙 410083; 4.湖南金天钛业科技有限公司,湖南 常德 413000; 5.中航飞机起落架有限公司,湖南 长沙 410200; 6.中国第二重型机械集团德阳万航模 锻有限责任公司,四川 德阳 618000; 7.西安三角防务有限公司,陕西 西安 710089; 8.中国科学院金属研究所,辽宁 沈阳 110016) 摘 要 以新型高强韧 TC18 钛合金为例,研究了双圆锥台试样热变形规律。 基于热模拟试验,获得 TC18 钛合金的高温流变行为, 采用 Arrhenius 双曲正弦函数建立 TC18 钛合金的本构方程。 对双圆锥台试样热变形进行有限元模拟,获得了压缩后变形试样纵向 截面的应变梯度分布。 对比双圆锥台试样不同应变区域与相同变形条件下 GLEEBLE 热模拟试样的 α 相平均晶粒尺寸,结果表明 一个双圆锥台试样上可以快速获取多个不同应变条件的热变形组织。 关键词 钛合金; TC18; 双圆锥台; 热变形; 应变梯度分布; 有限元模拟 中图分类号 TG146文献标识码 Adoi10.3969/ j.issn.0253-6099.2020.05.032 文章编号 0253-6099(2020)05-0128-06 Finite Element Simulation for Thermal Deformation of TC18 Titanium Alloy Specimen on a Double Cone MAO Ying1, LEI Yu2, DING Xu1, ZHANG Xiao⁃yong3, WANG Bing⁃feng1,3, FAN Kai4, FENG Kang⁃tun5, XIE Jing6, WANG Hai⁃peng7, LEI Jia⁃feng8 (1.School of Material Science and Engineering, Central South University, Changsha 410083, Hunan, China; 2.Chinalco Luoyang Copper Industry Co Ltd, Luoyang 471039, Henan, China; 3. State Key Laboratory of Powder Metallurgy, Central South University, Changsha 410083, Hunan, China; 4.Hunan Gold Sky Titanium Technology Co Ltd, Changde 413000, Hunan, China; 5.AVIC Landing Gear Advanced Manufacturing Corp, Changsha 410200, Hunan, China; 6.China National Erzhong Group Deyang Wanhang Die Forging Co Ltd, Deyang 618000, Sichuan, China; 7.Xi′an Triangle Defense Incorporated Company, Xi′an 710089, Shaanxi, China; 8.Institute of Metal Research, Chinese Academy of Sciences, Shenyang 110016, Liaoning, China) Abstract A new type of TC18 titanium alloy with high strength and toughness was taken as a sample to study the thermal deformation law of the specimen on a double cone. Based on the thermal simulation experiments, the high⁃temperature rheological behavior of TC18 titanium alloy was obtained, and a constitutive equation was established for the TC18 titanium alloy with the Arrhenius hyperbolic⁃sine function. The thermal deformation of the double cone specimen was simulated by using finite⁃element method, and the strain gradient distribution at longitudinal section of the deformed sample after compression was obtained. Based on the comparison of the mean grain size of α⁃phase in the double cone specimen in different strain region and GLEEBLE thermal simulation specimen under the same deformation conditions, it is found that several thermal deformation structures under different strain conditions can be quickly obtained from the double cone specimen. Key words titanium alloy; TC18; double cone; thermal deformation; strain gradient distribution; finite element simulation TC18 钛合金是一种高强韧近 β 型钛合金,具有强 度高、塑韧性好以及耐高温、耐腐蚀等优点,被广泛用 于航天航空和海洋军事等领域[1-3]。 该合金由前苏联 航空材料所研制,是退火状态下强度最高的钛合金。 近年来,随着高通量设计概念的提出,不少学者基 于这一概念提出了各种试验设计方法,双圆锥台试样 ①收稿日期 2020-04-22 基金项目 粉末冶金国家重点实验室自主课题(621021721);湖南省创新科技项目(2017GK2292);国家自然科学基金面上项目(51771231) 作者简介 毛 颖(1997-),女,湖北荆门人,硕士研究生,主要研究方向为钛合金高通量变形技术与微观结构。 通讯作者 汪冰峰(1978-),男,湖南岳阳人,副教授,博士,主要从事材料动态行为研究工作。 第 40 卷第 5 期 2020 年 10 月 矿矿 冶冶 工工 程程 MINING AND METALLURGICAL ENGINEERING Vol.40 №5 October 2020 变形技术是当前材料高通量测试的热点技术,也是一 种快速获取材料变形参数的很有前途的试验方法。 Tan 等人[4]设计了一种高温镍基合金的双圆锥台试 样,通过有限元模拟获得试样梯度应变分布,并探讨了 不同应变状态下合金的组织情况和热力学行为。 本文 选用 TC18 合金高强韧钛合金为研究对象,采用高通 量双圆锥台试样设计这一理念,通过圆柱形和双圆锥 台两种形状试样的热压缩试验,对比验证双圆锥台试 样热变形试验可行性,并试图在解决热变形过程中变 形参数与组织演变关系的同时减少试验次数,提高试 验效率。 1 试验材料与过程 试验材料为热锻 TC18 钛合金,其名义成分为 Ti⁃5Al⁃5Mo⁃5V⁃1Cr⁃1Fe,合金相变点在 870 ℃左右[5-6]。 热模拟压缩试验的试样为 Φ8 mm 12 mm 的圆 柱试样,压缩试验在 Gleeble3800 热模拟试验机上进 行。 将试样加热到目标温度保温 3 min,随后以一定的 应变速率将试样压缩至设定变形量。 变形温度分别为 790,820,850,880 和 900 ℃,应变速率为 0.01,0.1,1.0 和 10 s -1 ,设定应变分别为 0.2,0.3 和 0.4。 压缩后,水 冷淬火保留原始金相组织。 淬火后试样沿压缩方向线 切割,经研磨、抛光和 Kroll 试剂(2 mL HF+4 mL HNO3 +100 mL H2O)腐蚀,在光学显微镜下观察微观组织并 采用 NanoMeasurer 软件测量其中 α 相晶粒尺寸。 双圆锥台试样热压缩试验的试样加工尺寸及外形 如图 1 所示,压缩试验在 MTS 322T 型工作台试验机 上进行。 将双圆锥台试样加热到目标温度保温 15 min, 随后以一定的应变速率将试样压缩至设定变形量。 目 标温度为 880 ℃,应变速率为 0.01 s -1 ,设定变形量为 8 mm。 压缩后,水冷淬火。 沿压缩方向切开,随后进 行镶嵌、金相砂纸打磨、抛光和腐蚀,腐蚀液配比为 5 mL HF+10 mL HNO3+85 mL H2O,在光学显微镜下 观察不同区域的微观组织并采用 NanoMeasurer 软件 测量其中 α 相晶粒尺寸。 图 1 双圆锥台试样加工尺寸及外形图 (a) 加工尺寸; (b) 试样外形 2 试验结果与讨论 2.1 建立 TC18 钛合金的本构方程 在不同变形条件下,TC18 钛合金经热模拟压缩试 验所得的真应力⁃真应变曲线如图 2 所示。 可见真应 力⁃真应变曲线存在明显峰值应力,试样在热压缩过程 中发生了动态再结晶[7-10]。 真应变 -400 -300 -200 -100 0 -0.10.0-0.2-0.3-0.4-0.5 真应力/ MPa 真应变 -300 -250 -200 -150 -100 -50 0 -0.10.0-0.2-0.3-0.4-0.5 真应力/ MPa 真应变 -180 -150 -120 -90 -60 -30 0 -0.10.0-0.2-0.3-0.4-0.5 真应力/ MPa 真应变 -100 -80 -60 -40 -20 0 -0.10.0-0.2-0.3-0.4-0.5 真应力/ MPa 790 ℃ 820 ℃ 850 ℃ 880 ℃ 900 ℃ ● ● ■ ■ ▲ ▲ ▲ ▲ ◆ 790 ℃ 820 ℃ 850 ℃ 880 ℃ 900 ℃ ● ■ ▲ ▲ ◆ ◆ ● ■ ▲ ▲ ◆ a b 790 ℃ 820 ℃ 850 ℃ 880 ℃ 900 ℃ ● ● ■ ■ ▲ ▲ ▲ ▲ ◆ ◆ ● ■ ▲ ▲ ◆ c 790 ℃ 820 ℃ 850 ℃ 880 ℃ 900 ℃ ● ■ ▲ ▲ ◆ d 图 2 不同变形条件下 TC18 钛合金的真应力⁃真应变曲线 (a) ̇ ε=0.01 s -1 ; (b) ̇ ε=0.1 s -1 ; (c) ̇ ε=1.0 s -1 ; (d) ̇ ε=10 s -1 921第 5 期毛 颖等 TC18 钛合金双圆锥台试样热变形有限元模拟研究 采用文献[11]基于指数型本构关系 ̇ ε=A1σn1和幂 指数型本构关系 ̇ ε=A2exp(βσ)提出的 Arrhenius 双曲 正弦函数来建立 TC18 钛合金的本构方程 ̇ ε = A [sinh(ασ)] nexp[ - Q/ (RT)] (1) 式中 ̇ ε 为应变速率,s -1 ;Q 为变形激活能,J/ mol;σ 为 流变应力,MPa;T 为绝对温度,K;R 为摩尔气体常数, R=8.314 J/ (molK);A、A1、A2、n、n1、α 和 β 均为与材 料有关的常数。 分别对 ln ̇ ε⁃lnσ 和 ln ̇ ε⁃σ 进行一元线性回归处理, 拟合结果见图 3。 可得 n1=5.191 94,β=0.033 08 MPa -1 , 因此 α=β/ n1=0.006371 MPa -1 。 lnσ/MPa 3 0 -3 -6 4.03.54.55.05.56.0 lnε/s-1 . 790 ℃ 820 ℃ 850 ℃ 880 ℃ 900 ℃ ● ■ ▲ ▲ ◆ ● ● ●■ ■ ▲ ▲ ▲ ▲ ▲ ▲ ◆ ◆ ◆ ●■ ■ ▲ ▲ ◆ σ/MPa 3 0 -3 -6 1000200300400 lnε/s-1 . 790 ℃ 820 ℃ 850 ℃ 880 ℃ 900 ℃ ● ■ ▲ ▲ ◆ ● ● ● ■ ■ ▲ ▲ ▲ ▲ ▲ ▲ ◆ ◆ ◆ ●■ ■ ▲ ▲ ◆ a b 图3 流变应力与应变速率的关系 (a) lṅ ε⁃lnσ; (b) lṅ ε⁃σ 作出 lṅ ε⁃ln[sinh(ασ)]和 ln[sinh(ασ)]⁃T -1 的线 性回归曲线如图 4 所示。 分别取两条曲线斜率的平均 值,算得变形激活能Q=303.366 9 kJ/ mol。 根据文献[12],材料在高温下发生塑性变形时应 变速率与温度之间的关系可用参数 Z 表示 Z = ̇ εexp Q RT = A[sinh(ασ)] n (2) 式中 Z 为温度补偿应变速率因子。 对 lnZ⁃ln[sinh(ασ)] 进行一元线性回归处理,如图 5 所示。 得到 n=3.767 81 和 A=3.055 991013。 将上述求得的 A,α,n 和 Q 代入 Arrhenius 函数模 型,得到 TC18 钛合金热压缩变形时的本构方程为 ̇ ε = 3.056 1013 [sinh(6.371 10 -3 )]3.768 exp[ - 303.367 103/ (RT)] (3) ln[sinhασ] 3 0 -3 -6 -1.50.00.51.0-1.0-0.51.52.0 lnε/s-1 . 790 ℃ 820 ℃ 850 ℃ 880 ℃ 900 ℃ ● ■ ▲ ▲ ◆ ● ● ● ■ ■▲ ▲ ▲ ▲ ▲ ▲ ◆ ◆ ◆ ●■ ■ ▲ ▲ ◆ T -1/10-3 K-1 2.0 1.0 0.0 -1.0 -2.0 0.840.860.880.900.920.94 ln[sinhασ] 0.01 s-1 0.1 s-1 1.0 s-1 10 s-1 ● ■ ▲ ▲ ● ● ● ● ● ▲ ▲ ▲ ▲▲ ▲ ▲ ▲ ▲ ▲ ■ ■ ■ ■ ■ a b 图 4 流变应力与应变速率、温度的关系 (a) lṅ ε⁃ln[sinh(ασ)]; (b) ln[sinh(ασ)]⁃T -1 ln[sinhασ] 38 36 34 32 30 28 26 -1.51.51.00.50.0-0.5-1.02.0 lnZ ■ ■ ■ ■■ ■■ ■ ■ ■ ■■ ■ ■ ■ ■ ■ ■ ■ ■ 图 5 lnZ⁃ln[sinh(ασ)]拟合直线 2.2 TC18 钛合金双圆锥台试样热压缩试验有限元模拟 将热模拟压缩试验所得本构方程导入有限元前处 理过程[13-14],通过 Deform 3D 有限元模拟计算,得到温 度 880 ℃、应变速率 0.01 s -1 热压缩条件下试样截面的 等效应变云图如图 6 所示。 试样中心与试样端面应变 量存在一定的梯度,从中心部位到边缘应变量逐渐降 低,热压缩后试样心部到边缘的应变值与其相对位置 的关系曲线如图 7 所示。 图 8 为 880 ℃、0.01 s -1 条件下试验与 Deform 模拟 所得力⁃位移曲线对比图。 可见模拟值与试验值的最 031矿 冶 工 程第 40 卷 大误差约为 12.5%,平均误差在 10%以内。 图 6 等效应变云图 心部至边缘距离/mm ■ ■ ■ ■ ■ ■ ■ ■ ■ 0.8 0.6 0.4 0.2 0.0 4081216 真应变 图 7 双圆锥台试样截面水平轴线上应变分布曲线 位移/mm ■ -35 -30 -25 -20 -15 -10 -5 0 -20-4-6-8 力/kN ● 模拟曲线 试验曲线 ■ ● 图 8 880 ℃、0.01 s -1 下力⁃位移曲线试验与有限元模拟对比 图 9 为 TC18 在温度 880 ℃、应变速率 0.01 s -1 条 件下热压缩至不同应变的金相图。 由图可知,不同应 变程度下,试样中 α 相含量不同,TC18 钛合金在热压 缩过程中发生了相变。 图 10 为双圆锥台试样纵截面选取的金相观察区 域,图 11 为双圆锥台形 TC18 试样在 880 ℃、0.01 s -1 压缩条件下试样截面不同区域的金相组织图,其中 (a) ~(e)分别对应图 10 中的相应区域(下同)。 由图 可知,从边缘到心部随着应变量增加,α 相含量不同。 压缩过程中,β 相发生了明显的动态再结晶现象。 采用 NanoMeasurer 软件测量图 9 和图 11 金相组 织中 α 相晶粒尺寸,并使用 Origin2017 软件绘制得到 相应 α 相晶粒尺寸统计图,如图 12~13 所示。 图 9 880 ℃、0.01 s -1 不同应变下热模拟金相组织 (a) ε=0.2; (b) ε=0.3; (c) ε=0.4 图 10 双圆锥台试样纵截面金相观察区域 图 11 880 ℃、0.01 s -1 热压缩双圆锥台试样截面不同区域金相 131第 5 期毛 颖等 TC18 钛合金双圆锥台试样热变形有限元模拟研究 晶粒尺寸/μm 0.4 0.3 0.2 0.1 0.0 0.01.02.03.04.0 比例 a 晶粒尺寸/μm 0.4 0.3 0.2 0.1 0.0 0.01.02.03.04.0 比例 b 晶粒尺寸/μm 0.4 0.3 0.2 0.1 0.0 0.01.02.03.04.0 比例 c 图 12 880 ℃、0.01 s -1 不同应变下热模拟组织中α相晶粒尺 寸统计图 (a) ε=0.2; (b) ε=0.3; (c) ε=0.4 由图 12~13 可知,圆柱形试样热模拟压缩后,在 ε=0.4,0.3 和 0.2 应变下分别对应的 α 相平均晶粒尺 寸约为 2.52,2.36 和1.89 μm。 双圆锥台试样热压缩 后,截面(a)、(b)、(c)、(d)、(e)对应的 α 相平均晶粒 尺寸约为 2.65,2.67,2.20,2.02 和 1.70 μm。 对比可 知,ε=0.4 时圆柱形试样热压缩组织与双圆锥台试样 截面区域(a)、(b)的组织较为相似,区域(a)、(b)模 拟对应的应变范围为 0.42~0.75。 ε= 0.3 时圆柱形试 样热压缩组织与双圆锥台试样截面区域(c)的组织较 为相似,区域(c)模拟对应的应变范围为 0.24~0.35。 ε=0.2 时圆柱形试样热压缩组织与双圆锥台试样区域 (d)、(e)的组织较为相似,区域(d)、(e)模拟对应的 应变范围为 0.07~0.18。 总的来说,双圆锥台试样热 压缩试验的组织结果与相同变形条件下的圆柱形热模 拟试样组织结果较为一致。 晶粒尺寸/μm 0.4 0.3 0.2 0.1 0.0 0.01.02.03.04.0 比例 a 晶粒尺寸/μm 0.4 0.3 0.2 0.1 0.0 0.01.02.03.04.0 比例 b 晶粒尺寸/μm 0.45 0.30 0.15 0.00 0.01.02.03.04.0 比例 c 晶粒尺寸/μm 0.4 0.3 0.2 0.1 0.0 0.01.02.03.04.0 比例 d 晶粒尺寸/μm 0.5 0.4 0.3 0.2 0.1 0.0 0.01.02.03.04.0 比例 e 图 13 880 ℃、0.01 s -1 热压缩双圆锥台试样截面不同区域α相 晶粒尺寸统计图 231矿 冶 工 程第 40 卷 3 结 论 1) 通过热模拟压缩试验获得了 TC18 钛合金在 790~900 ℃下的真应力⁃真应变曲线。 曲线有明显峰 值应力,试样在热压缩过程中发生了动态再结晶,采用 Arrhenius 函数模型建立了 TC18 钛合金的热变形本构 方程。 2) 在 880 ℃、̇ ε=0.01 s -1 变形条件下,有限元模拟 得到双圆锥台试样截面应变呈梯度分布。 3) 双圆锥台试样不同应变区域 α 相晶粒尺寸与 相同应变条件下圆柱形热模拟试样 α 相晶粒尺寸较 为一致,说明双圆锥台试样热变形有限元模拟所得应 变分布结果较为准确,一个双圆锥台试样上可以快速 获取多个不同应变条件的热变形组织。 参考文献 [1] 莱茵斯皮特尔斯. 钛与钛合金[M]. 陈振华,译. 北京化学工 业出版社, 2005. [2] Lin Y C, Huang J, He D G, et al. Phase transformation and dynamic recrystallization behaviors in a Ti55511 titanium alloy during hot com⁃ pression[J]. Journal of Alloys and Compounds, 2019,795471-482. [3] Zherebtsov S V, Murzinova M A, Klimova M V, et al. Microstructure e⁃ volution during warm working of Ti⁃5Al⁃5Mo⁃5V⁃1Cr⁃1Fe at 600 and 800 ℃[J]. Materials Science and Engineering A, 2013,563168-176. [4] Tan L M, Huang Z W, Liu F, et al. Effects of strain amount and strain rate on grain structure of a novel high Co nickel⁃based polycrys⁃ talline superalloy[J]. Materials & Design, 2017,13160-68. [5] 刘岚逸,汪冰峰,张晓泳. TC18 钛合金中剪切局域化及其微观结 构研究[J]. 矿冶工程, 2018,38(4)149-151. [6] Luo J, Wang L F, Liu S F, et al. The correlation between the flow behavior and the microstructure evolution during hot working of TC18 alloy[J]. Materials Science and Engineering A, 2016,654213-220. [7] 尹 勇,谭 真,肖振宇. 微量 Ca、Y 元素对 AZ91 镁合金动态再 结晶及力学性能的影响[J]. 矿冶工程, 2019,39(1)114-117. [8] 郭俊成,肖振宇,杨续跃. 两步锻压提升铸态 AZ80 镁合金延伸率 的研究[J]. 矿冶工程, 2018,38(4)144-148. [9] Sakai T, Belyakov A, Kaibyshev R, et al. Dynamic and post⁃dynamic recrystallization under hot, cold and severe plastic deformation condi⁃ tions[J]. Progress in Materials Science, 2014,60(1)130-207. [10] Jia B H, Song W D, Tang H P, et al. Hot deformation behavior and constitutive model of TC18 alloy during compression[J]. Rare Met⁃ als, 2014,33(4)383-389. [11] Sellars C M, Mctegart W J. On the mechanism of hot deformation[J]. Acta Metallurgica, 1966,14(9)1136-1138. [12] Zener C, Hollomon J H. Effect of strain rate upon plastic flow of steel[J]. Journal of Applied Physics, 1944,15(1)22-32. [13] 彭柳锋,郭道强,魏玉勇,等. 2A12 铝合金锻件成形的有限元模 拟[J]. 矿冶工程, 2017,37(5)99-104. [14] 张 莉,李升军. DEFORM 在金属塑性成形中的应用[M]. 北 京机械工业出版社, 2009. 引用本文 毛 颖,雷 雨,丁 旭,等. TC18 钛合金双圆锥台试样热变 形有限元模拟研究[J]. 矿冶工程, 2020,40(5)128-133. (上接第 123 页) [8] Lv K, Li J, Qing X, et al. Synthesis and photo⁃degradation applica⁃ tion of WO3/ TiO2hollow spheres[J]. Journal of Hazardous Materials, 2011,189(1-2)329-335. [9] 曾爱香,熊惟浩,王采芳,等. 溶胶⁃凝胶法制备空心微珠表面钡铁 氧体包覆层的研究[J]. 材料保护, 2004,37(9)19-20. [10] 张 姣,江学良,余 露,等. 水热法制备二氧化铈纳米空心球及 其吸附性能研究[J]. 材料研究学报, 2016,30(5)365-371. [11] 黄 彬,江 奇,谢德钰,等. 模板法制备聚苯胺空心微球材料及 其电化学性能[J]. 功能材料, 2011,42(1)144-147. [12] 张 珂,汪太生,李自力,等. 细乳液聚合法制备聚合物空心球[J]. 高分子学报, 2018(4)475-481. [13] 叶敏伟. 新型粉末喷涂装备[J]. 电镀与涂饰, 2011,30(10)73-75. [14] 吴子健. 热喷涂技术与应用[M]. 北京 机械工业出版社, 2006. 引用本文 张胜全,何亮亮,王 胜,等. 铜渣尾矿制备空心陶瓷微球及 其表征[J]. 矿冶工程, 2020,40(5)120-123. (上接第 127 页) [14] Fehir J, Richard J. Differential polarization of spin and charge densi⁃ ty in substituted phenoxy radicals[J]. Journal of Physical Chemistry A, 2009,113(32)9249-9260. [15] Vogt P, Padova P D, Quaresima C, et al. Silicene compelling ex⁃ perimental evidence for graphenelike two⁃dimensional silicon[J]. Physical Review Letters, 2012,108(15)155501. [16] Mizoguchi H, Okunaka M, Kitano M, et al. Hydride⁃based electride material, LnH2( Ln = La, Ce, or Y) [ J]. Inorganic Chemistry, 2016,55(17)8833-8838. 引用本文 台俊杰,肖 璇. 桥位吸附硼原子硅烯的结构与电子性质[J]. 矿冶工程, 2020,40(5)124-127. 331第 5 期毛 颖等 TC18 钛合金双圆锥台试样热变形有限元模拟研究
展开阅读全文

资源标签

最新标签

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

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

矿业文库合伙人QQ群 30735420