Fatigue properties of AZ31 magnesium alloy welded joint by double-sided friction stir welding
-
摘要: 针对搅拌摩擦单面焊两侧热输入不均匀性导致疲劳强度低的问题，采用双面对称搅拌摩擦焊方法对10 mm厚的AZ31镁合金板材进行焊接，并研究其疲劳性能. 结果表明，双面对称搅拌摩擦焊接头的屈服强度为130 MPa，与单面焊的屈服强度123 MPa相比提高了5%；其疲劳极限为88 MPa，比单面焊接头的50 MPa提高了76%；双面对称接头疲劳裂纹萌生在上/下侧的前进侧位置，并跨越上/下侧焊缝交界面，最终在下/上侧焊缝的后退侧RS区域瞬断，其疲劳断口均为以解理特征为主的脆性断裂. 双面对称焊接头其中一面应变范围与单面搅拌摩擦焊的应变较高的后退侧接近. 通过双面搅拌摩擦焊接的镁合金接头疲劳强度得到了大幅度提升，疲劳寿命得到了延长.Abstract: In order to solve the problem of the low fatigue strength of single-sided welding (SSFSW) caused by the uneven heat input, double-sided friction stir welding (DSFSW) was used to weld AZ31 magnesium alloy with a thickness of 10 mm. The results show that the yield strength of DSFSW joint is 130 MPa, which is 5% higher than that of the SSFSW joint. The fatigue limit of the DSFSW joint is 88 MPa, which is 76% higher than that of the SSFSW joint. The fatigue crack of DSFSW joint starts at the advancing side of the upper/lower side, crossing the interface of the upper/lower weld seam and breaking in the retreating side (RS) of the one weld seam finally. The fatigue fracture surface of DSFSW is brittle fracture mainly composed of cleavage. The strain range of one side of DSFSW joint is close to the RS strain of SSFSW. The fatigue strength and life of magnesium alloy joint welded by DSFSW have been greatly increased.
-
0. 序言
激光焊接中高能束能量传输是通过狭长的激光匙孔实现的，此时匙孔开口直径不足1 mm，而焊接熔透特征区域位于匙孔底部，且比匙孔开口直径还要小1个数量级属于介观尺度范畴[1]. 因此，采用直接法测试的熔透特征信号易被周围干扰信号掩盖. 同时，施焊中匙孔上方还存在大量等离子体、金属飞沫、金属蒸气焰及粒子团簇等喷射物质，具有较强的光辐射、光线遮挡及散射等干扰作用，难以实现对匙孔底部熔透行为的直接有效检测. 因此，业内一般采用熔池曲面形态、焊缝热场分布、孔壁声学共振特征等间接测试方法[2-6]. 虽然这些方法能在一定程度上获取激光焊接过程熔透状态的某些特征信息，但是间接测量受电、磁场、气流及温度等环境影响因素较多、不确定性大，测试可靠性难以保证. 文中介绍一种针对Q235钢激光焊接熔透识别的直接检测方法. 该方法利用不同发光体的谱段特性获取匙孔内部受激辐射行为透视实像，通过研究熔透特征区域在匙孔成像中的位置及介观尺度下的寻位测试方法，有效提高测试敏感度. 并利用统计学方法实现检测信号的趋势性特征与熔透状态的关联识别，实现熔透特征抽象化分离，并获得明显的试验效果. 文中拟通过介观尺度信号萃取与统计学特征分析两项新技术的有机结合，寻求一种激光焊接熔透特征可靠在线识别的新途径.
1. 试验装置及试验条件
1.1 试验装置
激光焊接熔透检测系统由特效光路单元、图像采集单元/光电检测单元、数据采集与处理单元及微距监控单元组成. 其中，特效光路单元与激光焊接头同轴监控口连接，其内部光路是一组由变倍镜头、变焦镜头、物镜镜头组成的具有增长工作距离、加大景深、汇聚视角功能的光学镜头组，可将焊接头同轴监控口处的光信号汇聚成放大的光学实像. 图像采集单元选取FastCam Ultima 512高速摄像机，用于观察同轴影像特征. 光电检测单元由三维微调机构、高敏光感元件组成，可通过针孔摄取及介观寻位方式将感应芯片高精度定位至熔透特征区域，实现熔透有效识别信号的最大化增强. 数据采集与处理单元由检测工控机、集成化控制模块、预处理芯片及高速数据采集卡组成，用于对光电传感器的传输信号进行数据采集、处理，并具有计算获得焊接熔透状态等功能. 微距监控单元由微距镜头、微型相机、图像处理板卡及显示器组成，用于介观尺度下的辅助寻位与调节. 检测装置实物及内部组成如图1所示.
![]()
图 1 熔透检测装置及内部组成Figure 1. Penetration experimental apparatus and its internal composition. (a) penetration experimental setup; (b) internal composition1.2 试验条件
试验针对Q235钢板，板厚3 ~ 10 mm. 选用6 kW通快叠式激光器，焊接参数为激光功率2 ~ 4 kW，离焦量为0 mm，焊接倾角13°，氩气保护流量为20 L / min，焊接速度为2 m / min. 焊接激光光束入射至碳钢板形成激光匙孔，用配有970 nm ± 20 nm红外滤光片的光学镜头组观察激光枪同轴的匙孔图像，激光匙孔形态由FastCam Ultima 512高速摄像机拍摄获得，频率选取在(128 × 128) pixels下使用200 fps采样频率. 熔透识别信号由高敏光感元件获取，采样频率10 kHz.
2. 试验原理
2.1 熔透特征信号提取
覆盖在匙孔上方的等离子体、金属蒸气焰等干扰信号的连续谱几乎覆盖了整个紫外和可见光谱段，但是当进入到近红外谱段后(λ > 900 nm)就已经得到了大幅衰减，而此时匙孔内的一种新的辐射光信号就显现出来. 该信号是激光匙孔内壁在剧烈的能量交互作用下被激发出的微弱荧光[7]，其辐射波的增强范围大致在900 ~ 1 200 nm. 由于该辐射信号产生于匙孔内壁表面，且与匙孔形态变化有直接关联，因此可作为反映匙孔内壁变化的直接检测信号. 考虑到试验所采用的摄像机感应极限在950 nm附近，所以试验选取滤光片的截止范围是970 nm ± 20 nm.
焊接头同轴光路信号在进入光学镜头组前先经过红外滤光，利用不同发光体的谱段特性，有效分离并抑制覆盖匙孔上方的干扰信号，使匙孔内壁荧光信号得到有效增强，然后通过光学镜头组汇聚成放大实像，借助高速摄像可以从施焊点上端捕捉到匙孔内壁受激光能流激发的荧光辐射信号. 通过对比4 kW激光入射至3 mm厚钢板发生熔透焊接和入射至8 mm厚钢板非熔透焊接时激光匙孔内的热辐射映像差异，可知在熔透时匙孔热辐射照片中存在一个明显的低灰度暗区，如图2所示. 其产生的原因为：激光熔透焊接时会存在激光贯穿母材时刻，此时激光束会在母材的上下表面产生两个开口，上端开口较大，下端开口较小，由于匙孔底部缺口区域失去了荧光效应，因此在匙孔上方观测到的热辐射信号就会明显衰减，而产生了图中的低灰度暗区. 同时，该区域存在周期性特征，并随着焊缝熔透情况的加大停留时间和面积会明显增加，因此，可称为熔透特征区域. 由此，可通过检测匙孔下端熔透特征区域产生的幅值、时间及频率情况，来定性及定量监测激光焊接熔透情况.
![]()
图 2 非熔透和熔透状态下的匙孔内部辐射情况Figure 2. Internal radiation of the key holes in non-penetration and penetration state. (a) non-penetration state; (b) penetration state在功能实现上，试验将激光焊接头同轴提取的辐射信号，经由具有增长工作距离、加大景深及汇聚视角功能的特效光学镜头组汇聚成放大的光学实像. 高速摄像可以拍到匙孔热辐射图像，是因为镜头将激光同轴的匙孔实像打在高速摄像的CCD成像面上，所以熔透特征信号提取也需要在镜头的成像面上完成. 另一方面，由于激光焊接头引导光红光的投射位置几乎与激光匙孔荧光映像位置完全重合，所以红光可作为调节检测点位置的有效依据. 在保持光路不变的情况下，移除高速摄像的CCD，并在其成像面上放置一块毛玻璃. 此时，在垂直上方的微距监控单元拍摄的毛玻璃背面中，可以观察到红光在镜头成像面上的位置及大小. 通过红光确定匙孔位置，并根据熔透特征区域与匙孔的相对位置关系，可在监控显示器上标出熔透特征区域的位置及大小，如图3所示.
![]()
图 3 微距监控系统拍摄的毛玻璃成像Figure 3. Imaging on frosted glass taken by micro-distance monitoring system取下毛玻璃，并将底端带有针孔的传感器定位座固定在镜头组成像面上. 此时，在微距监控系统的显示器中会清晰看到针孔的位置和之前标出的熔透特征区域的位置，如图4a所示. 通过三维微调机构将针孔逐渐移至标记位置，调节误差精度应在0.02 mm以内，如图4b所示.
![]()
图 4 调节针孔位置Figure 4. Adjust pinhole position. (a) before adjustment position; (b) after adjustment position将高敏度传感器装入定位座内，由此可实现传感器光感面的准确定位. 利用传感器对介观尺度下荧光辐射强度的高分辨能力，实现对熔透有效识别信号的最大化增强.
2.2 熔透特征信号处理
由于激光焊接熔透特征区域的产生具有周期性特征，同时检测信号会伴随焊接状态及匙孔波动产生急剧变化，因此检测数据中存在着复杂的时域及频域变化. 但前期研究发现，如果在焊接条件不变的情况下，一段时间内选取的所有检测信号其整体宏观分布会趋于一致[7-8]，而且所收集的信号数量越多分布特征越明显. 这说明决定信号分布特征的主要因素不在于检测信号随时域及频域的变化过程，而取决于焊接热反应本质，即实际焊接条件及焊接状态[9-10]. 因此，通过统计学方法对信号产生的规律性点阵进行特征提取与分析，可以实现对焊接熔透行为的关联识别[11-12].
试验设计4 kW激光功率焊接3 mm厚和8 mm厚Q235钢板，通过介观萃取方法获取二者在1 s内的光感检测信号分布信息，分别如图5所示.
![]()
图 5 检测信号空间分布Figure 5. Spatial distribution of detection signal. (a) signals in penetration characteristic region; (b) signals in non-penetration characteristic region从图5不难发现二者在空间分布上均为有规律的点阵，不同检测信号强度区域内所获得的信号点数量存在明显差异. 说明激光焊接熔透状态与熔透特征区域辐射行为存在规律性变化特征，因此试验拟通过数学统计方法获得信号在空间分布的相对密度. 具体运算方法为：利用VB软件将检测数据按照时间顺序分解成一系列的连续分析区间，每个区间取100个采样数据点(间隔约10 ms的检测数据量)，然后将每个区间内的数据点按升序排列，并从最小的数据点开始依次获取其后的第20个点与该点的差值，比较得到差值最小的值域范围即为信号分布的最密集区域. 同时，该区域的中值可视为信号最集中出现的强度数值. 计算公式为
$$\gamma = \min ({D_{n + 20}} - {D_n}),n \in (1,2,3 \cdots 80)$$ (1) 式中：γ为信号分布的最密集的值域长度；Dn为区间内的第n个点，Dn + 20为Dn后的第20个点. 通过计算及比对可获得每个分析区间的γ值和对应的n值.
$$\beta = {D_n} + \frac{1}{2} \cdot \gamma $$ (2) 式中：β为信号最集中出现的强度数值.
由于β数值附近出现的信号频率最多，概率最大，作用最集中，因此可作为与焊接熔透关联的一项重要准稳态特征参数. 依此方法获得每个分析区间内的特征参数动态连接，可获得激光焊接熔透行为检测曲线. 由于该方法属于一种趋势性捕捉方法，所以能够有效屏蔽掉焊接时信号波动及随机干扰信号对检测结果产生的不利影响，大幅简化复杂光学信号的分析难度，实现对焊接熔透信息的有效识别.
3. 试验结果与分析
3.1 不等厚试板熔透检测试验
为了验证采用信号介观萃取和统计分析方法对相同焊接参数下激光焊接熔透识别效果，试验选取10 mm厚Q235钢板作为焊接对象，焊前在钢板背面加工了2个同样为7.5 mm深的矩形凹槽. 采用2 kW激光功率，以2 m/min焊接速度对试验板进行表面堆焊. 试验板无凹槽处获得的是未熔透焊缝，而在矩形槽处获得的是适度熔透的焊缝，通过信号介观萃取和统计分析方法的检测结果如图6所示.
![]()
图 6 不等厚试板熔透检测结果Figure 6. Penetration test result for non-equal thickness test plate. (a) welding plate; (b) signals in penetration characteristic region; (c) data analysis of penetration detection图6a为加工的不等厚试板实物照片，图6b是采用信号介观萃取方法获得的熔透特征区域检测信号，图6c为采用信号统计分析方法获得的激光焊接熔透行为检测曲线. 对比试板实际熔透情况不难发现，经由信号介观萃取方法获取的检测信号对激光焊接熔透行为具有较好的感知效果. 该信号在经过信号统计分析方法处理后可以获得明确的熔透识别结果. 同时，该检测结果与焊接实际熔透情况完全一致.
3.2 变功率熔透检测试验
为了验证在相同板厚情况下采用文中方法的熔透识别可靠性，试验选取3 mm厚Q235钢试板作为焊接对象，采用激光功率2 300 W – 2 900 W – 2 300 W – 2 900 W方波式功率输出堆焊形式，获得了熔透与未熔透交替出现的焊缝. 此时由于激光功率变化和熔透特征均较小，所以信号介观萃取方法获得的检测信号中肉眼已经很难分辨出信号的波动情况. 但是采用信号统计分析方法处理后的激光焊接熔透行为检测曲线却能清晰地区分出熔透与未熔透区域，而且检测到的结果与实际焊缝熔透情况完全一致，如图7所示.
![]()
图 7 变功率熔透试验结果Figure 7. Penetration test results for variable laser power. (a) signals in penetration characteristic region; (b) data analysis of penetration detection; (c) penetration result on the back of the welding plate4. 结论
(1)利用不同发光体的谱段特性分离并抑制干扰信号，使匙孔内壁荧光信号得到有效增强后，可作为激光焊接熔透识别的直接检测信号.
(2)通过采用针孔摄取及介观寻位方式将感应芯片高精度定位至熔透特征区域直接获取信号的方法可实现熔透有效识别信号的最大化增强.
(3)将检测信号通过统计学概率特征识别，可实现对焊接熔透特征信息的趋势性特征抽象分离，获得可靠的熔透检测结果.
-
![]()
图 1 双面FSW焊接过程及接头截面区域分布
Figure 1. DSFSW welding process and joint cross section area distribution. (a) DSFSW welding process; (b) cross section area distribution of joint
![]()
图 2 DSFSW接头横截面不同区域微观组织
Figure 2. Microstructure of various zones on the cross-section of DSFSW joint. (a) HAZ on the AS of the upper side weld; (b) HAZ on the RS of the lower side weld; (c) TMAZ on the AS of the upper side weld; (d) TMAZ on the RS of the lower side weld; (e) NZ
![]()
图 3 AZ31镁合金拉伸裂纹示意图
Figure 3. Schematic diagram of tensile crack of AZ31 magnesium alloy. (a) SSFSW joint; (b) DSFSW joint
![]()
图 6 SSFSW接头疲劳断口形貌及疲劳裂纹
Figure 6. Fatigue fracture morphology and fatigue crack of SSFSW joint. (a) fracture surface at crack initiation stage; (b) partial enlargement of Fig.6a; (c) fatigue striation; (d) cleavage plane in fatigue fracture; (e) fracture apperance of instantaneous fracture area
![]()
图 8 DSFSW接头疲劳断口形貌及疲劳裂纹
Figure 8. Fatigue fracture morphology and fatigue crack of DSFSW joint. (a) enlargement of region A; (b) enlargement of region B; (c) enlargement of region C
![]()
图 9 AZ31镁合金接头疲劳应力-应变曲线
Figure 9. Fatigue stress-strain curve of AZ31 Mg alloy joint. (a) SSFSW; (b) DSFSW
![]()
图 10 AZ31镁合金接头疲劳过程受力示意图
Figure 10. Stress diagram of AZ31 Mg alloy joint during the fatigue process. (a) sample loading; (b) SSFSW; (c) DSFSW
![]()
图 11 AZ31镁合金接头AS HAZ的EBSD结果分析
Figure 11. EBSD analysis results of AS HAZ on AZ31 Mg alloy joint. (a) Schmid factor distribution; (b) microstructure
![]()
图 12 AZ31镁合金接头RS HAZ的EBSD结果分析
Figure 12. EBSD analysis results of RS HAZ on AZ31 Mg alloy joint. (a) Schmid factor distribution; (b) microstructure
表 1 AZ31镁合金化学成分(质量分数，%)
Table 1 Chemical compositions of AZ31 Mg alloy
Al Mn Zn Ca Si Cu Ni Mg 3.02 0.6 1.10 0. 04 0.1 0.01 0.001 余量 
下载: 导出CSV
表 2 AZ31镁合金的力学性能
Table 2 Mechanical properties of AZ31 Mg alloy
抗拉强度
Rm/MPa屈服强度
Rp0.2/MPa断后伸长率
A(%)断面收缩率
Z(%)242 140 19.0 14.5
下载: 导出CSV
表 3 焊接设备参数
Table 3 Welding equipment parameters
焊接
方法搅拌头轴肩直径
d0/mm搅拌针
d1/mm搅拌针最大直径
dmax/mm单面焊 20 9.85 10 双面焊 20 4.85 5
下载: 导出CSV
表 4 优化的焊接参数
Table 4 Optimized welding parameters
焊接方法 搅拌速度w/( r·min−1) 焊接速度v/(mm·min−1) 单面焊 750 300 双面焊 1 200 300
下载: 导出CSV
-
[1] 许楠, 冯若男, 宋亓宁, 等. 微观组织对镁合金FSW焊缝应变硬化行为的影响[J]. 焊接学报, 2020, 41(11): 7−12. Xu Nan, Feng Ruonan, Song Qining, et al. Effects of microstructure on strain hardening behavior offriction stir welded magnesium alloy[J]. Transactions of the China Welding Institution, 2020, 41(11): 7−12.
[2] 郭少飞, 刘雪松, 张红霞, 等. 基于能量耗散的AZ31B镁合金接头疲劳极限快速评估[J]. 焊接学报, 2020, 41(12): 38 − 43. doi: 10.12073/j.hjxb.20200919002 Guo Shaofei, Liu Xunsong, Zhang Hongxia, et al. Rapid evaluation of fatigue limit of AZ31B magnesium alloy joints based on energy dissipation[J]. Transactions of the China Welding Institution, 2020, 41(12): 38 − 43. doi: 10.12073/j.hjxb.20200919002
[3] Kamal M, Rahman M M. Advances in fatigue life modeling: A review[J]. Renewable and Sustainable Energy Reviews, 2018, 82(1): 940 − 949.
[4] 魏巍, 孙屹博, 杨光, 等. 基于能量耗散的Q460焊接接头疲劳强度评估[J]. 焊接学报, 2021, 42(4): 49 − 55. doi: 10.12073/j.hjxb.20200907001 Wei Wei, Sun Yibo, Yang Guang, et al. Fatigue strength evaluation of Q460 weld joints based on energy dissipation[J]. Transactions of the China Welding Institution, 2021, 42(4): 49 − 55. doi: 10.12073/j.hjxb.20200907001
[5] Meng X C, Huang Y X, Cao J, et al. Recent progress on control strategies for inherent issues in friction stir welding[J]. Progress in Materials Science, 2020, 115: 100706.
[6] 秦丰, 周军, 侯振国, 等. 6082铝合金双面搅拌摩擦焊接头组织与性能[J]. 焊接学报, 2021, 42(2): 75 − 80. doi: 10.12073/j.hjxb.20201231001 Qin Feng, Zhou Jun, Hou Zhenguo, et al. Research on microstructure and properties of double-sided friction stir welding joint of 6082 aluminum alloy[J]. Transactions of the China Welding Institution, 2021, 42(2): 75 − 80. doi: 10.12073/j.hjxb.20201231001
[7] Kulwant S, Gurbhinder S, Harmeet S. Review on friction stir welding of magnesium alloys[J]. Journal of Magnesium and Alloys, 2018, 6(4): 399 − 416. doi: 10.1016/j.jma.2018.06.001
[8] Haghshenas M, Gerlich A P. Joining of automotive sheet materials by friction-based welding methods: A review[J]. Engineering Science and Technology, an International Journal, 2018, 21(1): 130 − 148. doi: 10.1016/j.jestch.2018.02.008
[9] Liu X C, Zhen Y Q, Sun Y F, et al. Local inhomogeneity of mechanical properties in stir zone of friction stir welded AA1050 aluminum alloy[J]. Transactions of Nonferrous Metals Society of China, 2020, 30(9): 2369 − 2380. doi: 10.1016/S1003-6326(20)65385-7
[10] Mabuwa S, Msomi V. Fatigue behaviour of the multi-pass friction stir processed AA8011-H14 and AA6082-T651 dissimilar joints[J]. Engineering Failure Analysis, 2020, 118: 104876. doi: 10.1016/j.engfailanal.2020.104876
[11] 陈洪胜, 王文先, 陈伟, 等. 镁/铝层合板FSW接头微观组织及力学性能[J]. 焊接学报, 2020, 41(3): 38 − 44. Chen Hongsheng, Wang Wenxian, Chen Wei, et al. Microstructure and mechanical properties of FSW joint of Mg/Al clad sheets[J]. Transactions of the China Welding Institution, 2020, 41(3): 38 − 44.
[12] Darmadi D B, Talice M. Improving the strength of friction stir welded joint by double side friction welding and varying pin geometry[J]. Engineering Science and Technology an International Journal, 2021, 24(3): 637 − 647. doi: 10.1016/j.jestch.2020.11.001
[13] Thakur A, Sharma V, Bhadauria S S. Effect of tool tilt angle on weld joint strength and microstructural characterization of double-sided friction stir welding of AZ31B magnesium alloy[J]. CIRP Journal of Manufacturing Science and Technology, 2021, 35: 132 − 145. doi: 10.1016/j.cirpj.2021.05.009
[14] Wang X P, Yoshiaki M. Interface development and microstructure evolution during double-sided friction stir spot welding of magne- sium alloy by adjustable probes and their effects on mechanical properties of the joint[J]. Journal of Materials Processing Technology, 2021, 296: 117104.
[15] 马宗义, 商乔, 倪丁瑞, 等. 镁合金搅拌摩擦焊接的研究现状与展望[J]. 金属学报, 2018, 54(11): 1597 − 1617. doi: 10.11900/0412.1961.2018.00392 Ma Zongyi, Shang Qiao, Ni Dingrui, et al. Friction stir welding of magnesium alloys: a review[J]. Acta Metallurgica Sinica, 2018, 54(11): 1597 − 1617. doi: 10.11900/0412.1961.2018.00392
[16] 周利, 张仁晓, 舒凤远, 等. Q235钢搅拌摩擦焊接头微观组织与力学性能分析[J]. 焊接学报, 2019, 40(3): 80 − 84. doi: 10.12073/j.hjxb.2019400076 Zhou Li, Zhang Renxiao, Shu Fengyuan, et al. Microstructure and mechanical properties of friction stir welded joint of Q235 steel[J]. Transactions of the China Welding Institution, 2019, 40(3): 80 − 84. doi: 10.12073/j.hjxb.2019400076
[17] Ashtiani H, Shayanpoor A A. New constitutive equation utilizing grain size for modeling of hot deformation behavior of AA1070 aluminum[J]. Transactions of Nonferrous Metals Society of China, 2021, 31(2): 345 − 357. doi: 10.1016/S1003-6326(21)65500-0
[18] Li L, Yang J, Yang Z Y, et al. Towards revealing the relationship between deformation twin and fatigue crack initiation in a rolled magnesium alloy[J]. Materials Characterization, 2021, 179: 111362. doi: 10.1016/j.matchar.2021.111362
[19] Yan Zhifeng , Wang Denghui, He Xiuli, et al. Deformation behaviors and cyclic strength assessment of AZ31B magnesium alloy based on steady ratcheting effect[J]. Materials Science & Engineering:A, 2018, 723(18): 212 − 220.
-
期刊类型引用(2)
1. 孙谦，黄瑞生，徐富家，曹浩，李林，宋扬，马强. 介观领域中的激光熔透信号物理特性研究. 焊接学报. 2024(07): 27-33+40 . 
本站查看
2. 陈月峰，王树昂，赵飞，陈金强，程定富，孙云翔. 自动焊图像处理技术的发展及应用研究. 焊接技术. 2021(03): 1-5 . 百度学术
其他类型引用(1)
计量
- 文章访问数: 331
- HTML全文浏览量: 41
- PDF下载量: 52
- 被引次数: 3
