文章快速检索  
  高级检索
钛合金薄壁件高速超声椭圆振动铣削机理和试验
张明亮1,2, 张德远1,3, 刘佳佳1, 高泽1, 韩雄4     
1. 北京航空航天大学 机械工程及自动化学院, 北京 100083;
2. 北京航空航天大学 沈元荣誉学院, 北京 100083;
3. 北京航空航天大学 生物医学工程高精尖创新中心, 北京 100083;
4. 航空工业成都飞机工业(集团)有限责任公司, 成都 610073
摘要: 针对航空领域中钛合金薄壁件在铣削过程中存在加工精度差、加工效率低等问题,提出了钛合金薄壁件高速超声椭圆振动铣削的方法。首先,在铣削加工中引入高速超声振动切削理念,使切削刀具刀尖附加椭圆振动进行超声频断续切削,可改善加工质量且突破了超声振动加工对临界速度的限制。然后,分析了该方法的分离原理,采用独立研制的超声椭圆振动铣削刀柄装置针对钛合金薄壁件进行了切削试验。试验结果显示,高速超声椭圆振动铣削相较于普通铣削,切削力降低20%~30%,且已加工表面让刀量降幅20%~30%,表面缺陷减少,表面粗糙度降低。
关键词: 高速超声椭圆振动     铣削     钛合金薄壁件     切削力     让刀量    
Mechanism and experiment of high-speed ultrasonic elliptical vibration milling of thin-walled titanium alloy parts
ZHANG Mingliang1,2, ZHANG Deyuan1,3, LIU Jiajia1, GAO Ze1, HAN Xiong4     
1. School of Mechanical Engineering and Automation, Beihang University, Beijing 100083, China;
2. Shenyuan Honors College, Beihang University, Beijing 100083, China;
3. Beijing Advanced Innovation Center for Biomedical Engineering, Beihang University Beijing 100083, China;
4. AVIC Chengdu Aircraft Industrial(Group) Co., Ltd., Chengdu 610073, China
Received: 2018-12-02; Accepted: 2019-01-18; Published online: 2019-03-18 16:05
Foundation item: National Natural Science Foundation of China (51475031)
Corresponding author. ZHANG Deyuan, E-mail: zhangdy@buaa.edu.cn
Abstract: A high-speed ultrasonic elliptical vibration milling method is proposed for poor processing quality and low processing efficiency in milling of thin-walled titanium alloy parts in aerospace. Firstly, high-speed ultrasonic elliptic vibration milling combines high-speed ultrasonic vibration cutting concept with milling. A high-frequency intermittent cutting where the cutting tool tip moves in accordance with the elliptical trajectory can effectively improve the quality of machining and break through the limitations of the ultrasonic vibration processing on the critical speed. Then, the separation principle of the process is analyzed, and a self-developed ultrasonic vibration milling device is used in the experiment with thin-walled titanium alloy parts. The experimental results show that, compared to the common milling processing, the cutting force reduction and the relieving amount reduction for high-speed ultrasonic elliptical vibration milling method are both about 20% to 30%, and meanwhile the surface defects and roughness are cut down.
Keywords: high-speed ultrasonic elliptical vibration     milling     thin-walled titanium alloy parts     cutting force     relieving amount    

钛合金具有较好的耐热性和耐腐蚀性,以及较高的比强度,尤其是在高温(300~400℃)条件下可保持这种特性,因此钛合金得到了广泛应用,已成为航空航天工程领域中较为重要的材料之一[1-5]

为了保证零件的可靠性及使用性能, 钛合金薄壁件对已加工表面质量、加工精度及加工效率要求较高。而钛合金薄壁件刚度较低,在加工中易导致变形大。目前,薄壁零件加工变形的研究集中于在线数控补偿技术、适当分配加工余量等辅助措施[6-7],或是改良走刀路径、切削参数、刀具参数和改良刀具结构等[8-9]

超声椭圆振动切削技术通过在加工工件或切削刀具上附加超声椭圆振动,使刀具刀尖相对于工件按椭圆运动轨迹进行振动,从而实现超声频断续振动切削[10-11]。与普通切削相比,超声椭圆振动切削技术的“分离特性”可使工件与刀尖分离,打开切削区,降低平均切削力,此时切削液可进入内部切削区,从而对切削区进行冷却,降低切削温度,改善刀具磨损并有效延长刀具寿命[12-16]。超声振动轨迹分离时,存在临界切削速度[17-18]。高速超声振动切削的提出打破了临界切削速度的限制,有助于提升加工效率[19]

本文将高速超声椭圆振动应用于铣削加工过程中,来解决钛合金薄壁件加工中的质量与效率问题,分析了其分离原理,并设计试验验证其在薄壁加工中的效果。与普通铣削加工对比,在一定工艺参数下,高速超声振动铣削能够有效降低钛合金薄壁件让刀量和已加工表面粗糙度。

1 高速超声椭圆振动铣削机理分析

超声椭圆振动切削的特点是切削刀具刀尖相对工件按照椭圆运动轨迹,对工件进行超声频断续切削。超声椭圆振动铣削加工过程中,铣刀随着主轴产生进给运动和回转运动的同时,亦在进行超声椭圆振动[20-21];即铣刀的刀尖运动轨迹是由以上3种运动合成,如图 1所示。以刀具回转轴为z轴,xy所在平面垂直于z轴且三坐标轴相互垂直建立直角坐标系。

图 1 试验用工件坐标轴 Fig. 1 Workpiece axis for experiment

图 1中坐标系为基准建立刀尖轨迹方程。工件如图 1所示,工件高度H远大于其厚度b。超声椭圆振动铣削运动过程中,铣刀刀尖的运动轨迹是由进给运动、主轴回转运动和超声椭圆振动这3种运动进行合成,所以超声椭圆振动铣削单个刀尖上任意一点P的运动轨迹为

(1)

m个刀尖上P点的运动轨迹为

(2)

式中:D为铣刀的直径;Vf为进给速度;ωz为椭圆振动的角速度,ωz=2πff为超声振动频率;α为铣刀转过的角度,α=ωitωi为主轴的回转角速度;AB分别为超声椭圆振动轨迹在椭圆长轴和短轴2个方向上的振幅。

对于超声椭圆振动切削,有两种切削情况:一种是低速超声椭圆振动切削,即切削速度小于临界切削速度(VL=2πAf),每个刀尖切削过程中进行高频椭圆运动,实现高频断续切削;另一种即高速超声椭圆振动切削,即切削速度大于临界切削速度VL,刀尖的椭圆运动轨迹被拉开,此切削过程要实现分离需要借助每齿之间运动轨迹的相位差、振幅以及每齿进给量等因素。

对于传统超声加工来说,当切削速度大于临界切削速度时,超声振动切削由断续切削变为连续切削,要实现刀具与工件的高频断续切削的条件更为苛刻,对于低速椭圆超声来说,只需满足每齿进给fzA+B,即可实现刀尖与工件之间的高频断续分离切削,较为容易实现;而对于高速超声而言,要实现刀尖与工件之间的高频断续分离切削的条件不仅要考虑到每齿进给量与振幅的关系,还要考虑每齿之间的相位差的关系,换言之,高速超声断续铣削的实现,是依靠每齿进给、振幅以及相邻齿之间相位差的关系来实现的。

低速情况下每个齿的切削自身都带有分离,而高速椭圆振动铣削是依靠相邻齿之间的表面轮廓的叠加实现分离的,如图 2所示,选取单个刀齿切削弧上一部分进行分析。工件表面轮廓是4个齿的刀尖轨迹形成的包络线,后3个齿的工件表面去除量都不同,也即每个刀齿的空切时间都不同,最终工件表面形成如黄色曲线的表面形貌。

图 2 刀尖轨迹分离情况 Fig. 2 Tool tip trajectory separation

为了便于观察分析第2个刀齿的切削机理(侧铣),如图 2所示,对于第2齿单个椭圆振动切削周期(A1A2)可分为4个阶段:①A1B阶段为空切阶段,此阶段切削力为零,刀尖逐渐向工件靠近,在B点处开始切入工件。②BC阶段为切入阶段,该阶段从B点开始逐渐切入工件,刀具对工件开始进行切削,同时切削由厚变薄,相当于低速椭圆振动铣削的切入阶段,并且符合顺铣的特点。③CD阶段为切出阶段,切削过程中切屑与前刀面产生的摩擦力方向发生反转,可以促进切屑的流出,前刀面摩擦力方向反转特性,使前刀面上产生积屑瘤的可能性大大降低,保证了刀具的锋利程度,降低了切削力,提高了刀具耐用度。④DA2阶段为分离阶段,该阶段铣刀刀尖与工件分离,切削力为零;该阶段切削区完全打开,切削液能够顺利进入切削区,使铣刀刀刃和工件得到充分冷却,切削温度得到大幅度降低,加工精度提高,有助于提升刀具切削性能。

对于m齿铣刀,当主轴转速为n(r/min)时,后一齿落后前一齿时间为,则由相位差造成的振动周期差值为

(3)

则由进给造成的振动周期差值为

(4)

式中:V为刀尖切削速度。

相邻两齿的振动周期差值为

(5)

图 2易知,只要相邻两齿相差的周期数ω=ω1+ω2+ω3即可,常数ω=(0,1),都可以在加工过程中产生分离。

同时,在刀具转过任意α后,刀尖N点法向切削厚度h=fzcos α < A+B,本试验中A+B=16.4 μm,fz < A+B可保证刀具旋转一周中都有分离。

2 试验 2.1 试验平台设计

为验证高速超声椭圆振动铣削在薄壁侧铣加工中的效果,设计了如图 3所示的试验系统。该系统主要由加工中心、工件、超声振动系统三部分组成。试验用加工中心为BV100立式加工中心。试验工件为钛合金(Ti-6Al-4V)薄壁件,尺寸为70 mm(长)×3.5 mm(厚)×60 mm(高)。工件固定在Kistler 9254测力仪上方,试验过程中可以直接采集力数据,测力仪固定于加工中心工作台。超声振动系统由超声电源、铣削换能器组成,换能器变幅杆末端是直径为12 mm的标准四刃硬质合金铣刀,采用微量润滑进行冷却,在铣削过程中机床开启主轴中心气冷,超声铣削换能器在气冷条件下能长时间稳定工作。

图 3 钛合金薄壁件高速超声椭圆振动铣削试验系统 Fig. 3 High-speed ultrasonic elliptical vibration milling experiment system for thin-walled titanium alloy workpiece
2.2 试验方案设计

本文所提的高速超声椭圆振动铣削方法主要针对精密切削加工,在进行高速超声椭圆振动铣削试验时,超声椭圆振动铣削刀柄的具体振动参数如表 1所示。

表 1 超声椭圆振动铣削实现振动参数 Table 1 Vibration parameters of ultrasonic elliptical vibration milling
参数 数值
频率f/Hz 17 880
振幅A/μm 8.6
振幅B/μm 7.8

试验方案设计如下:铣刀加工位置为薄壁件侧壁,为了准确测量2种切削方法的平均切削力并考虑减小让刀等因素,普通切削和高速超声椭圆振动铣削的切削力在一次切削走刀中完成测量。在试验过程中先进行普通铣削的切削力测量,等采集稳定信号之后再开启超声电源,转换为高速超声椭圆振动铣削,直到采集到稳定的切削力信号。

目前, 常用钛合金铣削切削速度为60~80 m/min。超声椭圆振动铣削的临界切削速度为:VL=2πAf=2π×8.6×17 880≈58 m/min,本试验高速超声椭圆振动铣削使用更高切削速度,具体参数如表 2所示。其中切削速度V=nπD,即试验中切削速度的变化通过改变机床转速。

表 2 加工参数 Table 2 Processing parameters
参数 数值
切削速度/(m·min-1) 100, 120, 140, 160, 180
每齿进给量/mm 0.015
切深/mm 5
切宽/mm 0.2

2.3 试验结果测量

本试验主要观测加工后让刀量和表面粗糙度,试验过程中测量切削力。让刀量通过螺旋测微仪测量加工后表面剩余厚度来衡量,螺旋测微仪型号为日本三丰量具mitutoyo,测量精度为0.001 mm,测量区域为直径4 mm的圆形,对高速超声椭圆振动铣削入刀位置中心和普通铣削出刀位置中心进行测量,两位置关于整个薄壁对称,以便减少其他因素干扰。

表面粗糙度使用Tylor Hobso粗糙度仪进行测量,测量精度为0.001 μm,测量3次取平均值,作为试验最后测量结果。为完整表征表面形貌,使用扫描电子显微镜(SEM),型号为JSM 6010,分别对高速超声椭圆振动铣削与普通铣削的表面形貌作对比。

试验过程中的切削力通过Kistler 9254测力仪进行测量,测量范围为0~20 kN,精度为0.01 N。测力仪采样频率设定为50 000 Hz,可对xyz及其复合的多方向力数据进行采集。

3 试验结果与讨论 3.1 让刀量

钛合金薄壁件精铣加工中,影响铣削精度的一个主要因素是在铣削过程中工件产生让刀现象,这种现象会导致实际切削深度小于理论值。在切宽0.2 mm,切深5 mm, 每齿进给0.015 mm时,已加工表面让刀量随切削速度变化如图 4所示。

图 4 让刀量随切削速度变化曲线 Fig. 4 Curves of relieving amount varying with cutting speed

随着切削速度的增大,高速超声椭圆振动铣削与普通铣削2种形式的加工让刀量都呈现出逐渐减小的趋势,但是在此过程中,超声椭圆振动铣削的让刀量始终小于普通铣削的让刀量,降幅在20%~30%之间,这很大程度上提高了钛合金薄壁件铣削的加工精度。

钛合金薄壁件相对于硬质合金铣刀来讲属于弱刚性零件,所以本文在模拟实际加工中薄壁件的让刀时采用如图 5所示的模型。普通铣侧面时,刀具整体受力可以简化为切削反力Fy和切深抗力Fx,由图可知影响让刀的切削力主要是Fx,即在切深抗力Fx作用下,产生让刀误差。

图 5 薄壁件切削让刀模型 Fig. 5 Model of relieving amount for milling thin-walled workpiece

高速超声椭圆振动铣削能够消除让刀量、提高铣削精度的原因主要有瞬间局部高能切削和大幅度降低铣削力两方面。

随着切削速度的增大,切削力成逐渐减小的趋势,相比普通铣削,高速超声椭圆振动铣削在xy两个方向的切削力都有明显下降;针对钛合金侧铣,主切削力是x方向的FxFx大小直接影响刀具寿命和薄壁件铣削精度。由图 6可以看出,随着切削速度V的增大,切削力逐渐减小,在此过程中,高速超声椭圆振动铣削方式在Fx方向的切削力始终低于普通铣削方式在Fx方向的切削力,下降20%~30%,这对于钛合金薄壁件的高速铣削具备优势。

图 6 切削力随切削速度变化曲线 Fig. 6 Curves of cutting force varying with cutting speed

在薄壁件加工中,让刀量主要受切削力影响,高速超声椭圆振动铣削分离机制带来的断续切削可以有效降低加工过程中的平均力,进而降低让刀量。

3.2 表面粗糙度

在切宽0.2 mm,切深5 mm,每齿进给0.015 mm时,已加工表面粗糙度值随切削速度变化如图 7所示。

图 7 表面粗糙度随切削速度变化曲线 Fig. 7 Curves of surface roughness varying with cutting speed

在固定其他参数不变时,切削速度逐步提高,由于钛合金薄壁件刚性不足且材料回弹大,高速小进给铣削过程中易导致颤振,普通加工与高速超声椭圆振动加工的表面粗糙度也都随之有小幅提高,在0.59~0.83 μm中变化,从图 7中可以看出,高速超声椭圆振动铣削在试验参数下能够稳定实现表面粗糙度值下降。

采用SEM拍摄已加工表面形貌如图 8所示,可以观察到普通铣削表面每一齿切过的痕迹清晰明显,整个表面由多个单齿加工形成的圆弧组成,而高速超声椭圆振动铣削表面,每齿加工过程中都形成了均匀致密的更为微细的结构。

图 8 试验工件表面形貌 Fig. 8 Surface topography of workpiece

进一步放大观察,如图 9所示,在普通铣削表面会有鳞刺缺陷,这些缺陷导致粗糙度值得增大。而高速超声椭圆振动铣削表面留有超声频冲击产生的均一致密的纹理,无明显表面缺陷,从微观上验证了高速超声椭圆振动铣削的精密性。

图 9 试验工件件微观形貌 Fig. 9 Micro-topography of workpiece

通过已加工表面形貌可以发现:高速超声椭圆振动铣削依靠相邻齿之间表面轮廓的叠加实现分离,使得每次刀尖切削钛合金过程中去除量更小,会形成规律致密的表面纹理,且均匀分布,没有明显的缺陷,一定程度上使表面粗糙度值降低。

4 结论

本文应用自行研制的超声椭圆振动铣削刀柄装置进行了钛合金薄壁件高速超声椭圆振动铣削加工精度与表面加工质量的试验研究,结合高速超声椭圆振动铣削原理和试验结果得到以下主要结论:

1) 高速超声椭圆振动铣削加工突破传统超声振动铣削加工对临界速度的限制,扩大了超声振动铣削的应用范围。铣削是机加工中应用最为广泛的加工方式,将高速超声理念与铣削结合,扩宽了高速超声理念的应用领域。

2) 相比于普通铣削,高速超声椭圆振动铣削过程中薄壁件所受到的切削力更小,降幅20%~30%,同时让刀量下降20%~30%,有助于提升薄壁件的加工精度。

3) 相比于普通铣削,高速超声椭圆振动铣削能够改善加工表面缺陷,降低加工表面粗糙度。

参考文献
[1]
裴东, 王波. 钛合金的切屑加工工艺综述[J]. 黑龙江科技信息, 2011(28): 95.
PEI D, WANG B. Review on chip processing technology of titanium alloy[J]. Heilongjiang Science and Technology Information, 2011(28): 95. DOI:10.3969/j.issn.1673-1328.2011.28.110 (in Chinese)
[2]
周子同, 陈志同, 熊曦耀, 等. 钛合金TB6侧铣表面完整性实验[J]. 北京航空航天大学学报, 2014, 40(6): 849-854.
ZHOU Z T, CHEN Z T, XIONG X Y, et al. Experiment on the surface integrity of side milling titanium TB6[J]. Journal of Beijing University of Aeronautics and Astronautics, 2014, 40(6): 849-854. (in Chinese)
[3]
LEYENS C, PETERS M. Titanium and titanium alloys[M]. Weinheim: Wiley-VCH, 2003.
[4]
刘彬, 刘延斌, 杨鑫, 等. TITANIUM 2008:国际钛工业、制备技术与应用的发展现状[J]. 粉末冶金材料科学与工程, 2009, 14(2): 67-73.
LIU B, LIU Y B, YANG X, et al. TITANIUM 2008:Development of international titanium industry, preparation technology and applications[J]. Materials Science and Engineering of Powder Matallurgy, 2009, 14(2): 67-73. DOI:10.3969/j.issn.1673-0224.2009.02.001 (in Chinese)
[5]
PETERS M, KUMPFERT J, WARD C H, et al. Titanium alloys for aerospace applications[J]. Advanced Engineering Materials, 2003, 5(6): 419-427. DOI:10.1002/adem.200310095
[6]
梅中义, 王运巧, 范玉青. 飞机结构件数控加工变形控制研究与仿真[J]. 航空学报, 2005, 26(2): 234-239.
MEI Z Y, WANG Y Q, FAN Y Q. Researching and simulating deformation of aerocraft structure part in NC machining[J]. Acta Aeronautica et Astronautica Sinica, 2005, 26(2): 234-239. DOI:10.3321/j.issn:1000-6893.2005.02.022 (in Chinese)
[7]
王运巧, 梅中义, 范玉青. 航空薄壁结构件数控加工变形控制研究[J]. 现代制造工程, 2005(1): 31-33.
WANG Y Q, MEI Z Y, FAN Y Q. Research on deformation in NC machining of aerospace thin-walled structure part[J]. Modern Manufacturing Engineering, 2005(1): 31-33. DOI:10.3969/j.issn.1671-3133.2005.01.013 (in Chinese)
[8]
TLUSTY J. High-speed machining[J]. CIRP Annals, 1993, 42(2): 733-738. DOI:10.1016/S0007-8506(07)62536-0
[9]
TLUSTY J, SMITH S, WINFOUGH W R. Techniques for the use of long slender end mills in high-speed milling[J]. CIRP Annals, 1996, 45(1): 393-396. DOI:10.1016/S0007-8506(07)63088-1
[10]
MORIWAKI T, SHAMOTO E. Ultrasonic elliptical vibration cutting[J]. CIRP Annals, 1995, 44(7): 31-34.
[11]
季远, 李勋, 张德远. 超声椭圆振动精密切削[J]. 航空制造技术, 2005(4): 92-95.
JI Y, LI X, ZHANG D Y. Ultrasonic elliptic vibration precision cutting[J]. Journal of Aviation Manufacturing Technology, 2005(4): 92-95. DOI:10.3969/j.issn.1671-833X.2005.04.023 (in Chinese)
[12]
SHEN X H, ZHANG J, XING D X, et al. A study of surface roughness variation in ultrasonic vibration-assisted milling[J]. International Journal of Advanced Manufacturing Technology, 2012, 58(5-8): 553-561. DOI:10.1007/s00170-011-3399-y
[13]
ZHANG C, ZHANG J, FENG P. Mathematical model for cutting force in rotary ultrasonic face milling of brittle materials[J]. International Journal of Advanced Manufacturing Technology, 2013, 69(1-4): 161-170. DOI:10.1007/s00170-013-5004-z
[14]
GRAŽEVIČIŪTE J, SKIEDRAITE I, JŪRENAS V, et al. Applications of high frequency vibrations for surface milling[J]. Mechanika, 2015, 69(1): 46-49.
[15]
CHERN G L, CHANG Y C. Using two-dimensional vibration cutting for micro-milling[J]. International Journal of Machine Tools & Manufacture, 2006, 46(6): 659-666.
[16]
SHEN X H, ZHANG J H, LI H, et al. Ultrasonic vibration-assisted milling of aluminum alloy[J]. International Journal of Advanced Manufacturing Technology, 2012, 63(1-4): 41-49. DOI:10.1007/s00170-011-3882-5
[17]
房善想, 赵慧玲, 张勤俭. 超声加工技术的应用现状及其发展趋势[J]. 机械工程学报, 2017, 53(19): 22-32.
FANG S X, ZHAO H L, ZHANG Q J. The application status and development trends of ultrasonic machining technology[J]. Journal of Mechanical Engineering, 2017, 53(19): 22-32. (in Chinese)
[18]
ABDULLAH A. Investigation of the effect of cutting speed and vibration amplitude on cutting forces in ultrasonic-assisted milling[J]. Proceedings of the Institution of Mechanical Engineers Part B:Journal of Engineering Manufacture, 2011, 226(7): 1185-1191.
[19]
SUI H, ZHANG X Y, ZHANG D Y, et al. Feasibility study of high-speed ultrasonic vibration cutting titanium alloy[J]. Journal of Materials Processing Technology, 2017, 247: 111-120. DOI:10.1016/j.jmatprotec.2017.03.017
[20]
张明亮, 姜兴刚, 刘佳佳, 等. 钛合金超声椭圆振动铣削参数对切削力的影响[J]. 电加工与模具, 2017(6): 39-41.
ZHANG M L, JIANG X G, LIU J J, et al. Influence of ultrasonic elliptical vibration milling parameters with titanium alloy on the cutting force[J]. Electromachining & Mould, 2017(6): 39-41. DOI:10.3969/j.issn.1009-279X.2017.06.010 (in Chinese)
[21]
辛文龙, 刘佳佳, 姜兴刚, 等.钛合金薄壁件腹板超声椭圆振动铣削精度研究[C]//第16届全国特种加工学术会议.苏州: 中国机械工程学会特种加工分会, 2015: 406-410.
XIN W L, LIU J J, JIANG X G, et al.Research on accuracy of ultrasonic elliptic vibration milling of thin walled titanium alloy webs[C]//The 16th National Conference on Non-Traditional Machining.Suzhou: Non-Traditional Machining Institution of the Chinese Mechanical Engineering Society, 2015: 406-410(in Chinese). http://cpfd.cnki.com.cn/Article/CPFDTOTAL-ZGVE201510003060.htm
http://dx.doi.org/10.13700/j.bh.1001-5965.2018.0712
北京航空航天大学主办。
0

文章信息

张明亮, 张德远, 刘佳佳, 高泽, 韩雄
ZHANG Mingliang, ZHANG Deyuan, LIU Jiajia, GAO Ze, HAN Xiong
钛合金薄壁件高速超声椭圆振动铣削机理和试验
Mechanism and experiment of high-speed ultrasonic elliptical vibration milling of thin-walled titanium alloy parts
北京航空航天大学学报, 2019, 45(8): 1606-1612
Journal of Beijing University of Aeronautics and Astronsutics, 2019, 45(8): 1606-1612
http://dx.doi.org/10.13700/j.bh.1001-5965.2018.0712

文章历史

收稿日期: 2018-12-02
录用日期: 2019-01-18
网络出版时间: 2019-03-18 16:05

相关文章

工作空间