汽车工艺与材料

元素	C	Mn	Si	Al	Nb	B	V	Cr	Ti
质量分数	0.30~0.25	1.30~1.80	0.60~1.00	0.03~0.05	0.03~0.06	0.002~0.003	0.03~0.06	0.80~1.20	0.02~0.04

元素	C	Mn	Si	Al	Nb	B	V	Cr	Ti
质量分数	0.30~0.25	1.30~1.80	0.60~1.00	0.03~0.05	0.03~0.06	0.002~0.003	0.03~0.06	0.80~1.20	0.02~0.04

材质牌号	冲击吸能值（峰值载荷时）/J	冲击吸能值（侵入量为80 mm时）/J
微合金1.8 GPa热成形钢	771.76	1 249.2
普通1.8 GPa热成形钢	616.81	1 075.2
1.5 GPa高强钢	449.7	1 013.4

材质牌号	冲击吸能值（峰值载荷时）/J	冲击吸能值（侵入量为80 mm时）/J
微合金1.8 GPa热成形钢	771.76	1 249.2
普通1.8 GPa热成形钢	616.81	1 075.2
1.5 GPa高强钢	449.7	1 013.4

微合金化超高强热成形车门防撞梁性能验证

PDF下载

冯毅 ¹^,² , 黄光杰 ¹ , 路洪洲 ³ , 张龙柱 ⁴ , 蒲霞 ⁵

汽车工艺与材料 | 中国汽车工程学会汽车材料分会第23届学术年会优秀论文 2023,(8): 29-37

收起

汽车工艺与材料 | 中国汽车工程学会汽车材料分会第23届学术年会优秀论文 2023, (8): 29-37

微合金化超高强热成形车门防撞梁性能验证

全屏

冯毅¹^,², 黄光杰¹, 路洪洲³, 张龙柱⁴, 蒲霞⁵

作者信息

¹ 重庆大学，重庆 400044

² 中国汽车工程研究院股份有限公司，重庆 401122

³ 中信金属股份有限公司，北京 100004

⁴ 河钢集团邯钢公司，邯郸 056001

⁵ 重庆长安汽车股份有限公司，重庆 400023

冯毅（1980—），男，正高级工程师，博士研究生，研究方向为汽车轻量化、汽车新材料、新工艺。

Performance Verification of Microalloyed Ultra High Strength for Hot Stamped Door Crash Beam

Yi Feng¹^,², Guangjie Huang¹, Hongzhou Lu³, Longzhu Zhang⁴, Xia Pu⁵

Affiliations

¹ Chongqing University, Chongqing 400044

² China Automotive Engineering Research Institute Co., Ltd., Chongqing 401122

³ CITIC Metal Co., Ltd., Beijing 100004

⁴ Handan Iron and Steel Company of Hegang group, Handan 056001

⁵ Chongqing Changan Automobile Co., Ltd., Chongqing 400023

出版时间: 2023-08-20 doi: 10.19710/J.cnki.1003-8817.20230176

文章导航

摘要

收起

将某新开发的微合金化1.8 GPa热成形钢牌号与传统钢牌号进行了从材料到零件级的安全性能测评。结果表明：相比于普通钢种，微合金化1.8 GPa热成形钢基于组织细化、第二相析出、残余奥氏体三大关键因素，具有更加明显的安全性优势。建立了2种材料的动态断裂模型，微合金钢在相同应力状态下具有更高的极限断裂应变，显示了更强的断裂抗力。对2种1.8 GPa热成形及1.5 GPa高强钢车门防撞梁进行了落锤冲击试验，微合金钢1.8 GPa热成形车门防撞梁有更加优异的抗碰撞侵入及碰撞吸能性能。

关键词

热成形钢 / 车门防撞梁 / 强韧性 / 微合金化 / 安全性能

Abstract

收起

The safety performance of a newly developed microalloyed 1.8 GPa hot stamped steel brand and the traditional brand was evaluated from material to part level. The results show that microalloyed 1.8 GPa hot stamped steel has more obvious safety advantages than ordinary steel based on 3 key factors: microstructure refinement, second phase precipitation and retained austenite. The dynamic fracture models of the two materials were established. The microalloyed steel has a higher ultimate fracture strain under the same stress state, showing a stronger fracture resistance. Drop weight impact tests were carried out on 2 kinds of 1.8 GPa hot stamped and 1.5 GPa high-strength steel door anti-collision beams. The micro alloy steel 1.8 GPa hot stamped door anti-collision beam has better anti-collision intrusion and collision energy absorption properties.

Key words

Hot stamped steel / Door crash beam / Strength and toughness / Microalloying / Safety performance

引用本文

冯毅, 黄光杰, 路洪洲, 张龙柱, 蒲霞. 微合金化超高强热成形车门防撞梁性能验证. 汽车工艺与材料, 2023 , (8) : 29 -37 . DOI: 10.19710/J.cnki.1003-8817.20230176

Yi Feng, Guangjie Huang, Hongzhou Lu, Longzhu Zhang, Xia Pu. Performance Verification of Microalloyed Ultra High Strength for Hot Stamped Door Crash Beam[J]. Automobile Technology & Material, 2023 , (8) : 29 -37 . DOI: 10.19710/J.cnki.1003-8817.20230176

正文

收起

1 前言

收起

热成形作为目前几乎唯一能够较好解决钢的强度与成形之间矛盾且最具产业化发展前景的超高强度汽车零件制造技术，近年来发展迅速。当前国内外热成形产线总数已超过200条，国内最早从2000年左右开始发展热成形技术，截至目前全国共计拥有约140条热成形生产线。当下国内热成形技术发展如火如荼，各汽车重点产业地区纷纷建起了热成形产线，热成形技术在乘用车领域上的应用已非常普遍^[1-2]。

随着近年来国内外汽车轻量化技术持续发展和汽车碰撞安全法规日趋严格，对于汽车典型安全构件用钢强度级别要求不断提升，而热成形钢顺应行业发展趋势，其强度已经迈向1.8~2.0 GPa强度级^[3-4]。当前1.8~2.0 GPa强度级热成形钢在国内乘用车领域的应用尚处于起步阶段，距离成熟推广尚有一定差距，亟待解决以下关键技术问题。

当前1 500 MPa热成形零件应用已经成熟，但是此强度条件下，延伸率多为6%~8%，塑性及韧性明显偏低（材料强塑积偏低），导致零件极限尖冷弯性能低下（根据VDA 238-100（07/2020） Plate Bending Test for metallic materials，极限尖冷弯角度达到60°以上才能有良好的吸能效果，而目前国内钢企生产的热成形钢均难以达到宝马（60°）公司和奔驰（65°）公司认证要求），具体反映在零件的碰撞吸能能力不足，在发生安全事故时冲击载荷无法充分被零件本身所吸纳而向乘员舱内传递，对车内乘员造成伤害。当零件强度进一步提升至1 800~2 000 MPa时，其延伸率更难保证，导致目前1 800 MPa强度级热成形零件一直未在国内外实现规模化应用。

热成形零件在汽车上一般是作为最为重要的安全构件，在实际碰撞过程中，热成形零件起到抗载、吸能作用。因此了解热成形钢材的动态特性，对模拟零件在冲击载荷下的变形和断裂，评价零件碰撞性能，改进零件设计和加工，提高汽车的安全性至关重要。当前，热成形钢材动态塑性失效准则的建立及其动态性能与热成形零件碰撞安全性能间的对应关系研究已成为国内外汽车行业重点关注的技术领域。其次，对零件成形、碰撞过程的模拟除需本构模型预测动态变形，还需要知道动态塑性失效准则，以预测零件的失效情况。目前应用较为普遍的动态断裂失效模型有Johnson-Cook、Xue-Wierzbicki、Damage Initiation and Evolution Criteria (DIEC)、Gissmo及最大剪切应力模型^[5-8]。

此外，最终还需要通过对热成形后的零件进行碰撞测试，直接测量零件的侵入量、承载冲击力和吸收能等指标来直接评价零件的碰撞性能。评价的标准主要分为两类：一类是以碰撞过程中部件的侵入量、冲击力、吸能量等指标为准，欧洲NCAP将碰撞速度下部件侵入量分为Good、Acceptable、Marginal和Poor，国内相应标准还在制定中。另一类标准是评价碰撞过程中驾驶员的头颈部、躯干、腿部等器官的损伤程度，这类标准还涉及医学领域，评价难度较大。总体而言，测评方法复杂而且成本高，一般作为零件最后定型后的测试检测。在此之前需要通过大量的数值模拟对碰撞过程进行仿真以评价和改进零件设计。未来发展趋势是对材料基础力学性能与零件使用性能相结合进行评价，通过建立材料强韧化组织、性能和零件最终服役性能间的纽带关系，更好地对整车轻量化及安全设计起到指导意义。

材料是零件性能基础，近年来微合金技术在热成形领域应用日趋普遍化^[9]。本文基于某微合金化的1.8 GPa热成形用钢，开出了某电动车型车门防撞梁零件。基于测试和仿真手段，建立了该热成形钢的材料级断裂本构模型，并其他材质车门防撞梁零件的动态安全性能（落锤冲击）测评对比，验证了微合金热成形钢对于提升车门防撞梁零件安全性能的技术优势。

2 试验材料

收起

本文车门防撞梁开发所选材料共计2种，第1种为某钢企1.8 GPa微合金热成形钢，其成分如表1所示。第2种材质为非微合金化处理的1.8 GPa微合金热成形钢。微合金钢成分设计思路如下：通过调整C、Si、Mn等主要合金元素含量，并添加微量的强韧化元素，全面提升热成形零件的综合服役性能水平。基于当前国内外热成形主流强韧化组织调控思路，合金成分的优化聚焦组织细化、第二相粒子析出及残余奥氏体协同增韧^[10]。

基于马氏体强度与其碳含量关系和对应强度要求，设计C质量分数为0.300%~0.035%。

基于钢种综合服役性能优化Ti、Al、Cr、B、Mn、Si等元素。其中，Ti与C、N结合形成化合物，可细化晶粒，提高钢的强度和韧性。Ti对钢的力学性能影响取决于它的存在形态和Ti/C质量分数比，微量Ti（0.03%~0.10%）使钢的屈服强度有所提高，但 Ti/C质量分数比超过4时，其强度和韧性急剧下降。Ti对钢的低温冲击韧性改善作用微弱，大颗粒、带棱角的TiN在形成过程中会降低钢机体的连续性，从而导致微裂纹的产生。钢中加入少量的Al，可细化晶粒，提高冲击韧性。Al还具有抗氧化性和抗腐蚀性能，与Cr、Si合用，可显著提高钢的高温不起皮性和耐高温腐蚀能力。Cr在奥氏体中扩散速度比较小，并阻碍C在奥氏体中的扩散，对贝氏体的细化有一定的效果，同时可提高钢的抗氧化性和耐蚀性。B含量一般控制在0.001%~0.005%，少量B使钢淬透性显著增加。B、Nb复合加入极大改善钢的显微组织和力学性能，随着B、Nb的增加，钢中珠光体逐渐降低，马氏体和贝氏体含量增加。B的质量分数超过0.005%时会出现连续网状B化物，降低钢的强韧性。Mn可扩大奥氏体相区，可通过固溶强化提高铁机体强度，有利于提高低合金钢强度。Mn还会阻碍富碳奥氏体进一步分解，在相变过程中抑制铁素体+珠光体的形成，促进粒状贝氏体的形成。此外，Mn适当增加可提高钢的强度，对韧性影响不大。Si的添加可以减缓C原子的扩散，阻碍C化物析出，抑制珠光体的形成，使钢中残余奥氏体数量增加。Si有利于提高钢的强度，钢中添加质量分数为1.5%~2%的Si，屈服强度可提高15%~20%。但是，Si过高会降低轧材表面质量，并在一定程度上恶化锻造、焊接和表面涂覆能力。

添加微量的Nb和V（质量分数为0.03%~0.06%）。Nb具有细化晶粒、析出强化的作用，添加一定量的Nb元素有利于提高钢的强度、韧性和改善钢的焊接性能，还可以提高钢的耐大气腐蚀能力，防止晶间腐蚀，改善焊接性能。少量的V可细化奥氏体晶粒，提高韧性，对钢的低温冲击韧性有利。但是，V过高会导致碳化物聚集长大，会降低钢的强度，若碳化物析出在晶界上还会降低钢的室温韧性。Nb、V复合添加时钢具有较好的强韧配合，原因是Nb与V的沉淀强化和 Nb的晶粒细化可以在不降低钢强度的同时改善韧性。

微合金1.8 GPa热成形钢在淬火态下的组织检测结果如图1所示。可以看出，通过复合微合金化处理使得该钢种具有非常细小的淬火马氏体板条组织，基体中析出了一定量的Nb、V等第二相颗粒，且形成了一定量的残余奥氏体组织，这些组织特征均对提升材料（零件）基体的韧性起到促进作用。

3 材料级断裂模型建立

收起

断裂模型建立流程如下：首先，进行2种1.8 GPa热成形钢在淬火态条件下的准静态单向拉伸试验，获得弹塑性力学性能数据，用于标定本构模型参数。通过R5缺口拉伸试验、R10缺口拉伸试验、纯剪切拉伸试验、拉剪切拉伸试验、中心孔拉伸试验及杯突试验，结合数字图像相关（DIC）设备，获得材料在不同应力状态下的断裂试验数据，用于标定LS_DYNA中的MMC断裂失效模型参数。利用标定后的材料本构模型及断裂失效模型参数，精确表征不同材料的变形及断裂特性。

试验试样如图2所示，每组试验重复3次。高速拉伸试验在ZWICK HTM 5020高速拉伸试验机上进行，其余均在CMT5305电子万能试验机上进行。选取LS_DYNA中的MAT_24号材料模型，基于准静态单向拉伸试验数据，结合试验数据和有限元的方法进行材料本构模型参数标定。断裂分析中，基于纯剪试验、拉剪试验、R5缺口拉伸试验、R10缺口拉伸试验、中心孔拉伸试验及杯突试验数据，选取LS_DYNA中MAT_ADD_ERISION内的MMC断裂失效模型进行断裂失效模型参数标定。

准静态单向拉伸试验、纯剪试验、拉剪试验、R5缺口拉伸试验、R10缺口拉伸试验、中心孔拉伸试验拉伸速率分别为3 mm/min、0.6 mm/min、0.6 mm/min、0.6 mm/min、0.6 mm/min、0.6 mm/min。杯突试验也在CMT5305电子万能试验机上进行，冲头速率为2 mm/min。高速拉伸试验应变率设为0.1 /s、1 /s、10 /s、100 /s和500 /s。试验过程中，准静态单向拉伸试验、纯剪试验、拉剪试验、R5缺口拉伸试验、R10缺口拉伸试验及中心孔拉伸试验均选用50 mm引伸计进行相对变形测量。高速拉伸试验选用8 mm虚拟引伸计进行测量。

如上所述，材料断裂模型的建立，包括材料本构模型以及断裂失效模型、断裂准则的建立。首先通过建立准静态单向拉伸试验数值模型，基于LS_DYNA中的MAT_24号材料本构模型，通过不同本构方程拟合得到不同真应力-塑性应变曲线，挑选拟合效果最好的2种模型进行归一化处理，选择不同归一化系数，通过仿真对比得到最好归一化参数的本构模型。其次，MMC断裂失效模型是目前国内外应用最为普遍的断裂模型类别，MMC断裂失效模型中，断裂时刻等效塑形应变是应力三轴度及归一化洛德角参数的函数^[11]。基于纯剪试验、拉剪试验、R5缺口拉伸试验、R10缺口拉伸试验、中心孔拉伸试验及杯突试验的试样尺寸及试验条件，建立各试验的数值模型。在LS_DYNA的MAT_24材料模型中输入标定后的真应力-塑性应变曲线，在不输入断裂失效模型参数的情况下进行模拟仿真。从仿真结果中选取各试验下断裂单元的应力三轴度。基于提取到的信息以及不同应力状态断裂试验中采用DIC测试得到的临界断裂应变，选取MMC断裂失效模型及GISSMO断裂准则进行参数标定及优化，使MMC断裂失效模型及GISSMO断裂准则能准确表征材料在复杂受力状态下的断裂失效行为。

2种1.8 GPa热成形钢的断裂性能对比如图3所示。可以看出，明显的微合金钢的断裂曲线或曲面总体均位于普通钢种之上，体现出了更加优异的断裂抗性，微合金钢热成形零件具有更加优异的碰撞安全性能优势。

车门防撞梁作为车身侧面碰撞中的重要承力零件，起到抵抗车门向内弯曲，保护乘员舱完整性的作用，是侧面碰撞中吸能的主要部件，其侵入量、侵入速度和变形模式决定乘员的安全。因此，开展了基于不同材质的目标车型车门防撞梁的动态落锤对比试验。

大能量动态落锤试验系统如图4所示。该试验设备主要用于测试材料和结构在动态冲击载荷下的变形行为和力学性能，是系统测试和研究材料结构在高应变率下的能量吸收行为的重要设备之一。设备参数如下：试验最大高度15 m；最大冲击速度17.3 m/s；桥压5 V，应变仪增益300 dB；红外线光栅触发。落锤试验机数据采集系统（图4）由NI-PXI采集卡（例如PXI6358）、接线盒、信号调理器以及动态应变放大器组成，主要用于测试位移、速度、振动、定时以及温度等信号。设备操作软件为NI SignalExpress2013。此外还配备了高速相机采集系统，用于记录结构变形过程。设备操作软件为Vic-Snap8，采样频率最高可达900 000帧/s。图5为落锤冲击样品外观，所配备的工装、及试验装夹过程。

图6~图8为2种成分的1.8 GPa热成形车门防撞梁及其对应的原替代1.5 GPa强度级的车门防撞梁，共计3种材质零件的落锤冲击过程记录。冲击过程中可分为接触、变形（最大动态变形位移值表征）、冲头卸载样品回弹（最大静态变形位移值表征）。试验结果表明，微合金1.8 GPa热成形车门防撞梁最大动态变形位移约为82.0 mm，最大静态变形位移为17.6 mm。普通1.8 GPa热成形车门防撞梁最大动态变形位移约为108.3 mm，最大静态变形位移为48.0 mm。1.5 GPa强度级的车门防撞梁的最大动态变形位移约为89.2 mm，最大静态变形位移为58.0 mm。图9为3种材质车门防撞梁经落锤试验后的外观。分析认为微合金1.8 GPa热成形车门防撞梁在冲击过程中的最大动态位移、最终静态位移均低于另外2种材质的零件，说明其具有更加优异的抗碰撞侵入性能，从这一方面验证了其安全性能优势。

3种材质车门防撞件的碰撞性能检测曲线如图10所示。首先，如图10a所示，碰撞过程中在同等时间条件下，微合金钢车门防撞梁的位移量始终最低；如图10b所示，随着碰撞过程进行，在同等时间条件下，微合金钢车门防撞梁的运动速度也始终处于最低状态。如图10c和图10d所示，不论随位移或时间变化，微合金钢车门防撞梁的峰值载荷力始终最高。综上显示微合金钢车门防撞梁具有最优异的抗碰撞变形能力。

图10e为3种材质防撞梁的碰撞吸能曲线示例，在整个碰撞过程中，微合金钢车门防撞梁的吸能量也是最高的，从另一角度体现了微合金钢基于其优异的韧性实现了零件更为优异的碰撞吸能能力。表2为侵入量为80 mm时3种材质零件冲击吸能值对比，可以看出微合金钢车门防撞梁的吸能量高于其余2种材质零件。上述落锤冲击试验结果验证了基于微合金热成形钢的性能优势和热成形车门防撞梁这一构件的的安全防护功能优势（注：图10中P1为微合金1.8 GPa热成形车门防撞梁；B1为普通1.8 GPa热成形车门防撞梁；Y1为原被替代的1.5 GPa高强钢防撞管）。

4 结束语

收起

本文将某新开发的微合金化1.8 GPa热成形钢牌号与传统牌号进行了从材料到零件级的安全性能测评。结果表明：微合金化1.8 GPa热成形钢基于组织细化、第二相析出、残余奥氏体三大关键因素，相比于普通钢种，具有更加明显的安全性优势。建立了2种材料的动态断裂模型，微合金钢在相同应力状态下具有更高的极限断裂应变，显示了更强的断裂抗力。对2种1.8 GPa热成形及1.5 GPa高强钢车门防撞梁进行了落锤冲击试验，微合金钢1.8 GPa热成形车门防撞梁有更加优异的抗碰撞侵入及碰撞吸能性能。

需要指出的是本文开展2种PHS1800钢种的断裂模型对比，构建断裂模型对应的测试样品是基于钢板平板模淬火处理得到的，其制备工况和实际的零件热成形工况存在差异(实际零件热成形后其各部位基体组织性能均存在一定的差异性，而钢板淬火后各部位组织性能具有较强的一致性)，在零件上取样经过测试和仿真处理构建断裂模型更有工程意义。因此，未来应致力于研究考虑零件取样组织性能差异的更高精度断裂模型的构建流程方法，并将其应用于热成形零件的安全性能开发。

此外，本文所涉微合金化的PHS1800钢种，其合金元素较多，有必要研究与其相适应的生产及其零件的热成形工艺控制方法，以充分发挥微合金化的技术优势，提升材料和零件的性能水平。还需要基于进一步的钢种及其的生产工艺优化研究成果，不断优化现有成分方案，以实现性能和成本更完美的契合。

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

[1]

金学军, 龚煜, 韩先洪, 等. 先进热成形汽车钢制造与使用的研究现状与展望[J]. 金属学报, 2020, 56(4): 411-428.

[2]

赵征志, 陈伟健, 高鹏飞, 等. 先进高强度汽车用钢研究进展及展望[J]. 钢铁研究学报, 2020, 32(12): 1059-1076.

[3]

梁江涛. 2 000 MPa级热成形钢的强韧化机制及应用技术研究[D]. 北京: 北京科技大学, 2019.

[4]

胡宽辉. 2 000 MPa级高强塑积热成形钢的研究[D]. 武汉: 武汉科技大学, 2019.

[5]

周佳, 梁宾, 赵岩, 等. 复杂应力状态下车用高强钢断裂失效行为表征与应用研究[J]. 塑性工程学报, 2021, 28(3): 153-163.

[6]

穆磊. 面向先进高强钢的韧性断裂预测模型研究与应用[D]. 北京: 北京科技大学, 2018.

[7]

BASARAN

. Stress State Dependent Damage Modeling with A Focus on The Lode Angle Influence[D]. Aachen: Rheinisch-Westf Lischen Technischen Hochschule Aachen, 2011.

[8]

BORRVALL

. A General Damage Initiation and Evolution Model (DIEM) in LS-DYNA[C]// European LS-DYNA Conference. Manchester, 2013: 1-13.

[9]

路洪洲, 赵岩, 冯毅, 等. 微合金化热成形钢开发应用进展及展望[J]. 机械工程材料, 2020, 44(12): 1-10.

[10]

梁江涛, 赵征志, 路洪洲. Nb微合金化技术在2 000 MPa级热成形钢中的作用[C]// 中国第二届汽车EVI及高强度钢氢致延迟断裂大会论文集, 2019.

[11]

BAI

, WIERZBICKI

. A New Model of Metal Plasticity and Fracture with Pressure and Lode Dependence[J]. International Journal of Plasticity, 2008, 24(6): 1071-1096.

2023年第卷第8期

PDF下载

184

引用本文

BibTeX

文章信息

doi: 10.19710/J.cnki.1003-8817.20230176

首发时间：2026-01-12
出版时间：2023-08-20

补充材料

相关文章

文章信息

作者

出版历史

基金

作者信息

¹ 重庆大学，重庆 400044

² 中国汽车工程研究院股份有限公司，重庆 401122

³ 中信金属股份有限公司，北京 100004

⁴ 河钢集团邯钢公司，邯郸 056001

⁵ 重庆长安汽车股份有限公司，重庆 400023

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/qcgyycl/CN/10.19710/J.cnki.1003-8817.20230176

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT