DOI:10.19476/j.ysxb.1004.0609.2019.04.12
正弦振动引起的BGA焊点Sn-Cu金属间化合物失效机理
胡 丹,沈 骏,陈 旭,翟大军,高润华
(重庆大学 材料科学与工程学院,重庆 400044)
摘 要:研究不同回流焊工艺下得到的球栅阵列(BGA)焊点在正弦振动疲劳试验中的失效行为。借助扫描电镜观察在不同加热因子下形成的Sn-Cu金属间化合物的形态和厚度,运用有限元模拟分析球栅阵列焊点在正弦振动下的应力集中和分布。结果表明:当金属间化合物层的形态和厚度不同时,裂纹的萌生和扩展机理不同;随着金属间化合物层厚度的增加,焊点的振动疲劳寿命先是缓慢提高,随后急剧下降;当金属间化合物层的厚度为1.5~3.0 μm时,接头的振动疲劳寿命达到最大值。
关键词:金属间化合物;球栅阵列;振动;疲劳寿命
文章编号:1004-0609(2019)-04-0756-08 中图分类号:TG 425.1 文献标志码:A
集成电路I/O密度的提高推动着现代电子封装产业倒装芯片技术的发展[1-2]。在倒装芯片技术中,焊点在起连接作用的同时,还起着导电、散热的作用,且单一焊点在服役过程中承受的载荷越来越大。而整个封装系统的功能依赖于焊点的完整性,因此焊点的完整性决定了整个系统的可靠性[3]。热循环引起的低周疲劳是表面贴装技术可靠性中最关心的问题,而振动产生的高周疲劳也严重影响封装的可靠性[4]。美国空军评估:约有20%航空电子的失效都是由振动和冲击引起[5]。电子元件在生产加工过程、运输和服役过程中总会受到振动载荷的作用,尤其是航空电子系统和汽车电子系统在服役过程中会经受大量的振动。因 此,由振动引起的失效在可靠性问题中应引起高度重视[6-7]。
目前研究人员对焊点热疲劳的研究较多,对振动疲劳的研究较少,因此许多研究者都很关注对处于振动环境下的电子元件疲劳寿命的预测[8-9]。CHOI等[10]通过分析温度对Sn3Ag0.5Cu(SAC305)印刷线路组件在振动环境下耐用性的影响,发现薄型BGA元件在反复的随机振动刺激下,其高温耐用性相比于室温耐用性有所降低。YU等[6]基于振动测试和有限元分析提出了一种预测电子元件在随机振动载荷下疲劳寿命的评估方法,基于雨流计数法、应力-寿命曲线和Miner法则计算出的疲劳寿命与实际测试结果一致。CHE等[11]运用振动疲劳测试和分析方法对倒装芯片的焊点进行了疲劳寿命评估,在这项研究中,累计破坏指标因子被用于预测疲劳寿命。CHEN等[12]结合振动疲劳测试、有限元分析和理论公式,预测了电子元件处于正弦振动载荷下的疲劳寿命,研究人员发现最大应力产生在塑料焊球阵列封装(Plastic ball grid array package,PBGA)元件边缘处的焊球上。YANG等[3]用循环扫描的方式表征了PBGA的动态性能,同样发现疲劳失效总是发生在PBGA元件中边缘处的焊球上。MAIO等[13-14]研究了不同振动频率下的振动可靠性,研究表明断口位置非常靠近IMC和焊点的界面处。YONG等[15]研究了多种无铅钎料在振动和不同温度环境下的可靠性,研究表明在所有的钎料中,裂纹都萌生于IMC和钎料基体的界面处,受振动的影响,裂纹沿着钎料基体缓慢扩展,但是,KO等没有考虑到IMC层形态的影响。TU等[16]研究了IMC对BGA焊点振动失效的影响,研究表明IMC层越厚,焊点的疲劳寿命越短。
上述文献中没有提及BGA焊点中Sn与Cu形成的IMC层受正弦振动引起的失效相关的机理,因此,本文研究了不同温度回流焊工艺下BGA焊点受正弦振动引起的失效行为,通过扫描电子显微镜观察了IMC层的形态和厚度,借助有限元模拟分析了BGA焊点在正弦振动下的应力集中。
1 实验
本实验所用试样的封装形式均为倒装球栅阵列(Flip chip ball grid array, FCBGA),FCBGA尺寸为35 mm×35 mm×3.63 mm,焊料为Sn3.0Ag0.5Cu,直径为0.6 mm,焊球间距为1 mm。PCB由FR4材料制成,尺寸为214 mm×214 mm×1.5 mm。将锡膏印刷到PCB上的铜焊盘之后,利用回流炉(HELLER 1900EXL)将FCBGA焊接在PCB上铜焊盘相应位置,在回流焊过程中设置了6组不同的温度工艺参数。所测液相温度对时间的积分近似等于t×ΔT,这个积分就是加热因子,是回流焊工艺的特征工艺参数[16]。回流焊工艺的具体加热因子、峰值温度和液相线以上停留时间如表1所列。
表1 试样的回流参数
Table 1 Reflow parameters of samples
多次回流焊后,将组装好的PCB放在振动台上,并用5个螺母柱进行固定。分别用两个传感器来测试试样的动态响应,其中一个位于振动台上,另一个位于PCB上,由此可确定测试设备的传递系数。通过扫频试验可以确定测试样品的固有频率,随后试样在一阶固有频率下进行定频正弦振动。用探测器实时监测振动过程中菊花链电阻,并记录失效时间。一旦焊点在振动过程中出现裂纹,菊花链的电阻就会增大。
振动测试后,用环氧树脂镶嵌试样,再按照标准金相步骤(粗磨、精磨、抛光、并用5% HNO3+92% C2H5OH+3% HCl(体积分数)腐蚀5 s)制备焊点截面金相试样。为了得到IMC层准确厚度,通过扫描电子显微镜(TESCAN VEGA 3 LMH SEM)观察焊点的形态并确定IMC层的厚度。IMC层的厚度由截面金相试样的扫描电镜图片结合图像分析得出:1) 在恰当的倍数下得到扫描电镜图片;2) 通过Image Tool 3.0软件计算IMC层的厚度;3) IMC层的平均等效厚度由式(1)计算得出:
(1)
式中:HIMC是IMC层的平均等效厚度;S是IMC层的面积;L是IMC层的长度。
由于焊球太小而不能直接测量其应力,故用ANSYS 13.0软件进行有限元分析,模拟焊点在拉伸载荷下的Mises应力分布。为了简化模型,焊点视为直径和高度均为0.6 mm的圆柱体。SOLID45 用于模拟元件,网格生成后的模型如图1所示。研究表明,最大应力位于PBGA元件中边缘的焊点上[3, 5-6]。包含焊球、Cu基板、Cu6Sn5在内的材料性能如表2所列。
图1 采用ANSYS 13.0构建的振动模型示意图
Fig. 1 Schematic diagrams of model constructed with ANSYS 13.0 for vibration test
表2 有限元分析模型中所用材料的性能[13, 15]
Table 2 Material properties used in FEA model[13, 15]
2 结果与分析
2.1 IMC微观结构演变
BGA焊点的微观结构演变过程如图2所示。IMC层的厚度记录于表3中。LAURILA等[17-19]的研究结果表明IMC层主要是由Cu6Sn5相构成,该相是在回流焊过程中由Sn3.0Ag0.5Cu与Cu形成的。当加热因子为30 s·℃时,IMC层平滑且较薄(厚度约为0.74 μm)(见图2(a));当加热因子为160 s·℃时,形成厚度为1.57 μm连续的波纹状Cu6Sn5层(见图2(b));当加热因子为390 s·℃时,形成厚度约为1.95 μm的圆滑的扇贝状Cu6Sn5层(见图2(c))。由图2(d)和2(f)可知,当加热因子增大时,扇贝状的Cu6Sn5层逐渐增厚,当加热因子达到1680 s·℃时,Cu6Sn5层整体上呈扇贝状,但是边缘有棱角(见图2(f))。
表3 IMC层厚度
Table 3 Thickness of IMC layer
在Sn与Cu作用初期,Sn、Cu原子之间相互反应,扩散距离较短,所以Cu6Sn5层较薄,Cu6Sn5层的生长行为遵从线性动力学[20-21]。随着加热因子的增大,形成了连续的波纹状IMC层,这导致在扩散过程中体扩散占主导,Cu6Sn5层的生长速度减慢。随后形成圆滑的波纹状的Cu6Sn5层,在这个阶段中,同时进行着IMC层晶界扩散、晶界开槽、晶粒粗化、溶解进熔化的钎料。当晶界扩散和晶界开槽在扩散过程中占主导时,Cu6Sn5层的生长速率增大。当加热因子较大时,IMC层虽然呈波纹状,但是边缘有棱角,Cu6Sn5层晶粒粗化降低了晶界强度,使凹槽作为扩散路径,因此Cu6Sn5层生长速度变慢。
2.2 焊点在振动载荷下的应力分布
有限元分析模型中检测了最大应力,焊点在振动载荷下的应力分布如图3所示。为简化计算,模拟时假设材料为各向同性弹性材料,由于对称性,可用四分之一截面进行计算。由图3可知,振动载荷作用在多相结构上会导致两种不同材料边界的尖角处产生应力集中,这是由于边界处相不连续,异种材料的力学性能不同,导致振动载荷下变形不协调造成的。为了研究Cu6Sn5层对应力集中的影响,保持钎料有限元模型的总高度不变,Cu6Sn5层的厚度分别设为5 μm、10 μm、15 μm。由此得到:在1 μN拉伸力下模拟时,焊点中的最大应力分别为2697.9 Pa,2733.7 Pa和2760.7 Pa,因此IMC层厚度的增加促进了边角的应力集中。实际上,在IMC生成时,体积减小,这也导致边角处应力集中。3种材料Cu6Sn5、铜和钎料的硬度分别为3.70 GPa、160HV和12~20HV[22],且最大应力出现在焊点中最硬的地方(Cu6Sn5层中)。
图2 加热因子不同时回流焊后IMC层的形貌SEM像
Fig. 2 SEM images of IMC layers after reflow at different heating factors
图3 焊点在拉伸载荷下的应力分布
Fig. 3 Von Misses stress distribution inside solder joint under tensile loading
2.3 IMC对焊点振动断裂机理的影响
图4 IMC厚度在1~2 μm左右时焊点振动失效断裂形貌
Fig. 4 Fracture morphologies of solder joint subjected to vibration load with IMC thickness of about 1-2 μm
图5 IMC厚度在2~3 μm左右时焊点振动失效断裂形貌
Fig. 5 Fracture morphologies of solder joint subjected to vibration load with interfacial IMC thickness of about 2-3 μm
图4所示为IMC厚度在1~2 μm左右时焊点振动失效断裂形貌。从图4(a)可见,IMC层的平均厚度为1~2 μm,IMC层界面薄且平坦,焊点瓶颈处出现穿晶断裂,随后裂纹在钎料基体中沿某角度扩展。如图4(b)所示,断口呈韧窝状,这是钎料留在表面的体现。图5所示为IMC厚度在2~3 μm左右时焊点振动失效断裂形貌。由图5(a)可知,当加热因子增大时,IMC层厚度增加到2~3 μm。裂纹同样起源于焊点的瓶颈处,但是在沿着钎料与Cu6Sn5层界面扩展。扇贝状的Cu6Sn5晶粒在断面图片中清晰可见(见图5(b)),因此这是典型的沿晶断裂。
裂纹起源于焊点的瓶颈处,由焊点结构及异种材料之间的热膨胀系数不匹配导致。根据SHEN等[23]的研究表明,Cu6Sn5的热膨胀系数(16.3×10-6 K-1)与Sn3.0Ag0.5Cu的热膨胀系数(25×10-6 K-1)有很大差距。基板被钎料固定而被完全约束,然而PCB在X轴方向和Y轴方向仍保有自由度,因此当温度改变时,PCB在X轴方向和Y轴方向会产生较大变形。这便导致钎料中产生较大的剪切力,使裂纹在瓶颈处萌生。当IMC层厚度增加到2~3 μm时,裂纹沿着钎料与Cu6Sn5层界面扩展,因为最大应力出现在IMC层中,且根据以上有限元分析可知最大应力随IMC层厚度的增加而增加,因此裂纹萌生的起始位置移至Cu6Sn5层。此时,不会发生断裂,裂纹也不在IMC层中扩展,因为IMC层的强度比铜和钎料的强度高[24],因此,裂纹倾向于沿着钎料与Cu6Sn5层界面不连续处扩展。
当加热因子继续增大时,断裂机理有所改变。图6所示为IMC厚度3~4 μm时焊点振动失效断裂形貌。由图6可知,当IMC厚度达到3~4 μm时,裂纹在突出的Cu6Sn5层尖端腰部处萌生,并且沿着Cu6Sn5层扩展。作为解理断裂,断面形态呈解理面,这是破裂的Cu6Sn5附在表面的体现(见图6(b))。当加热因子继续增大,IMC层的总厚度超过4 μm(见图7),裂纹开始萌生的位置移至Cu6Sn5层的内部,并且沿着铜和Cu6Sn5界面扩展(见图7(b)),同样呈解理断裂。
图6 IMC厚度3~4 μm时焊点振动失效断裂形貌
Fig. 6 Fracture morphologies of solder joint subjected to vibration load with IMC thickness about 3-4 μm
图7 IMC厚度大于4 μm时焊点振动失效断裂形貌
Fig. 7 Fracture morphologies of solder joint subjected to vibration load with IMC thickness of more than 4 μm
一般认为如果界面IMC层生长形式是平坦的,那么纵向形变应力为0。然而在实际中,界面IMC层在生长过程中倾向于变得粗糙,结果如图3所示。非均匀生长导致纵向形变应力积累,在Sn-Cu二元合金系中形成Cu6Sn5时体积减小[25-26]。IMC层的收缩导致 应力与周围材料不匹配,存在大量残余应力。FIELDS等[27]的研究表明,Cu6Sn5、Cu和Sn3.0Ag0.5Cu的断裂韧性分别是1.4、10和42[27]。对于低KIC的Cu6Sn5而言,一旦在突出的Cu6Sn5层表面形成裂纹或沟槽,裂纹或沟槽很快便会沿着尖端撕裂。值得注意的是:随着IMC层的生长,IMC层断裂机理由塑形断裂转变为脆性断裂。
2.4 IMC对焊点振动疲劳寿命的影响
采用6种不同温度工艺进行回流焊的试样的正弦振动循环实验的结果如表4所列。焊点寿命分布可用N50%来表征,N50%表示菊花链电阻增加到初始值50%的疲劳寿命。N50%与IMC层厚度的关系如图8所示。由图8可知:随着IMC厚度的增加,疲劳寿命先提高再降低。当IMC层厚度约为2 μm时,疲劳寿命最长。
表4 焊点的振动疲劳寿命
Table 4 Vibration fatigue life of solder joints
图8 焊点振动疲劳寿命和IMC层厚度的关系
Fig. 8 Relationship between statistical fatigue lifetime and thickness of IMC layer of solder joints
当IMC层非常薄且Cu6Sn5与钎料界面较光滑时,钎料不能稳固连接,极易从铜上剥离。当IMC层呈波纹状,且厚度增加到1.5~3.0 μm时,裂纹在钎料基体中萌生,并沿着某一角度扩展或沿着Cu6Sn5与钎料界面扩展。当裂纹沿Cu6Sn5与钎料界面扩展时,裂纹会受到突出的Cu6Sn5层的阻碍。由于IMC层的脆性断裂机理,IMC层太厚会提高焊点失效的可能性[28]。
3 结论
1) IMC厚度增加大幅促进两种不同材料共同边界边缘处应力集中,改变应力分布,由此使裂纹起源发生变化。
2) 随着IMC厚度的增加,断裂机理由塑形断裂转变为脆性断裂。当IMC层厚度为1~2 μm,界面IMC层薄且平坦时,裂纹在焊点瓶颈处萌生,并在钎料内沿着与钎料和Cu6Sn5界面约呈45°夹角扩展;当IMC厚度增加到2~3 μm时,裂纹扩展路径沿着钎料与Cu6Sn5界面;当IMC厚度达到3~4 μm时,裂纹在突出的Cu6Sn5尖端腰部处萌生,并沿着Cu6Sn5扩展;当IMC层厚度超过4 μm时,裂纹萌生的位置转移至Cu6Sn5内部,并沿着铜与Cu6Sn5界面扩展,实验效果符合模拟预期。
3) IMC厚度约为2 μm时,焊点的振动疲劳寿命达到最大值。
REFFERENCES
[1] HUANG Yue,LIN Chun,YE Zhen-hua,DING Rui-jun. Reflow flip-chip bonding technology for infrared detectors[J]. Journal of Micromechanics and Microengineering, 2015, 25(8): 1-6.
[2] TONG H M, LAI Y S, WONG C P. Advanced flip chip package [M]. USA: Springer, 2013: 31-45.
[3] 杨淼森, 孙凤莲, 孔祥霞, 周云芳. Sn-Ag-Cu无铅球栅阵列焊点塑性表征[J]. 中国有色金属学报, 2015, 25(11): 3119-3126.
YANG Miao-sen, SUN Feng-lian, KONG Xiang-xia, ZHOU Yun-fang. Plastic characterization and performance of Sn-Ag-Cu lead-free BGA solder joint[J]. The Chinese Journal of Nonferrous Metals, 2015, 25(11): 3119-3126.
[4] KIM J M, WOO S W, CHANG Y S, KIM Y J, CHOI J B, JI K Y. Impact reliability estimation of lead free solder joint with IMC layer[J]. Thin Solid Films, 2009, 517: 4255-4259.
[5] LIU F, MENG G. Random vibration reliability of BGA lead-free solder joint[J]. Microelectronics Reliability, 2014, 54(1): 226-232.
[6] YU D, AL-YAFAWI A, NGUYEN T T, PARK S, CHUNG S. High-cycle fatigue life prediction for Pb-free BGA under random vibration loading[J]. Microelectronics Reliability, 2011, 51(3): 649-656.
[7] 李雪梅, 孙凤莲, 张 浩, 辛 瞳. 微焊点Cu/SAC305/C固-液界面反应及电迁移行为[J]. 焊接学报, 2016, 37(9): 61-64.
LI Xue-mei, SUN Feng-lian, ZHANG Hao, XIN Tong. Micro-solder joints Cu/SAC305/Cu solid-liquid interfacial reaction and electromigration behavior[J]. Transactions of the China Welding Institution, 2016, 37(9) : 61-64.
[8] LIU F, LU Y, WANG Z, ZHANG Z. Numerical simulation and fatigue life estimation of BGA packages under random vibration loading[J]. Microelectronics Reliability, 2015, 55(12): 2777-2785.
[9] XIE X J. Fatigue life prediction of board level solder joints for special equipment under vibration loads[J]. Electronic Components and Materials, 2016, 35(6): 98-102.
[10] CHOI C, DASGUPTA A. Effect of temperature on vibration durability of SAC305 printed wiring assemblies[C]// Thermal and Thermomechanical Phenomena in Electronic Systems. IEEE, 2012: 745-752.
[11] CHE F, PANG J H. Vibration reliability test and finite element analysis for flip chip solder joints[J]. Microelectronics Reliability, 2009, 49(7): 754-760.
[12] CHEN Y S, WANG C S, YANG Y J. Combining vibration test with finite element analysis for the fatigue life estimation of PBGA components[J]. Microelectronics Reliability, 2008, 48(4): 638-644.
[13] MAIO D D, HUNT C P. High-frequency vibration tests of Sn-Pb and lead-free solder joints[C]// Electronics System-Integration Technology Conference. IEEE, 2008: 819-824.
[14] SONA M, PRABHU K N. Review on microstructure evolution in SAC solders and its effect on mechanical integrity of solder joints[J]. Journal of Materials Science (Materials in Electronics), 2013, 24( 9): 3149-3169.
[15] YONG H K, YOO S, CHANG W L. Evaluation of lead-free solder reliability under vibration at elevated temperature[J]. International Journal of Materials & Structural Integrity, 2014, 8(1/2/3): 53.
[16] TU P L, CHAN Y C, HUNG K C, LAI J K L. Comparative study of micro-BGA reliability under bending stress[J]. IEEE Transactions on Advanced Packaging, 2000, 23(4): 750-756.
[17] LAURILA T, VUORINEN V, KIVILAHTI J K. Interfacial reactions between lead-free solders and common base materials[J]. Cheminform, 2006, 37(30): 1-60.
[18] 鞠国魁, 韦习成, 孙 鹏, 刘建影. 界面金属间化合物对铜基Sn-3.0Ag-0.5Cu焊点拉伸断裂性能的影响[J]. 中国有色金属学报, 2007, 17(12): 1936-1942.
JU Guo-kui, WEI Xi-cheng, SUN Peng, LIU Jian-ying. Effects of interfacial IMC on tensile fracture behavior ofSn-3.0Ag-0.5Cu solder joints on copper substrates[J]. The Chinese Journal of Nonferrous Metals, 2007, 17(12): 1936-1942.
[19] KONG Yi-gang, KONG Zhi-gang, SHI Feng-min. Microstructure and mechanical property of Sn-Ag-Cu solder material[J]. Rare Metals, 2017, 36(3): 193-197.
[20] SHEN J, CHAN, Y C, LIU S. Growth mechanism of Ni3Sn4 in a Sn/Ni liquid/solid interfacial reaction[J]. Acta Materialia, 2009, 57(17): 5196-5206.
[21] PAWELKIEWICZ M, DANIELEWSKI M, JANCZAK- RUSCH J. Intermetallic layer growth kinetics in Sn-Ag-Cu system using diffusion multiple and reflow techniques[J]. Advanced Engineering Materials, 2015, 17(4): 512-522.
[22] LEE H T, CHEN M H, JAO H M, LIAO T L. Influence of interfacial intermetallic compound on fracture behavior of solder joints[J]. Materials Science and Engineering A, 2003, 358(1): 134-141.
[23] SHEN J, ZHAI D, CAO Z, ZHAO M, PU Y. Fracture behaviors of Sn-Cu intermetallic compound layer in ball grid array induced by thermal shock[J]. Journal of Electronic Materials, 2014, 43(2): 567-578.
[24] SHANG J K, ZENG Q L, ZHANG L, ZHU Q S. Mechanical fatigue of Sn-rich Pb-free solder alloys[J]. Journal of Materials Science Materials in Electronics, 2007, 18(1/3): 211-227.
[25] GONG J, LIU C, CONWAY P P, SILBERSCHMIDT V V. Evolution of CuSn intermetallics between molten SnAgCu solder and Cu substrate[J]. Acta Materialia, 2008, 56(16): 4291-4297.
[26] YANG Y, LI Y, LU H, CHEN J. Interdiffusion at the interface between Sn-based solders and Cu substrate[J]. Microelectronics Reliability, 2013, 53(2): 327-333.
[27] FIELDS R J, LOW S R Ⅲ, LUCEY G K. The metal science of joining[M]. Warrendale, PA: TMS, 1992: 165-173.
[28] 姚 健, 卫国强, 石永华, 谷 丰. 电迁移极性效应及其对 Sn-3.0Ag-0.5Cu无铅焊点拉伸性能的影响[J]. 中国有色金属学报, 2011, 21(12): 3094-3099.
YAO Jian, WEI Guo-qiang, SHI Yong-hua, GU Feng. Polarity effect of electromigration and its influence on tensile properties of Sn-3.0Ag-0.5Cu lead-free solder joint[J]. The Chinese Journal of Nonferrous Metals, 2011, 21(12): 3094-3099.
Fracture mechanism of Sn-Cu intermetallic compound layer in BGA solder joint induced by sinusoidal vibration
HU Dan, SHEN Jun, CHEN Xu, ZHAI Da-jun, GAO Run-hua
(School of Materials Science and Engineering, Chongqing University, Chongqing 400044, China)
Abstract: The sinusoidal vibration failure behavior of ball grid array(BGA) solder joints reflowed with different temperature profiles was studied. The morphologies and thickness of Sn-Cu intermetallic compound (IMC) formed under diverse heating factor were observed by scanning electron microscopy. The stress concentration and distribution of BGA solder joints under the sinusoidal vibration were analyzed by finite element simulation. The results show that the origin and propagation of cracks mechanism differs when the morphologies and thickness of IMC layer are different. The vibration fatigue lifetime of solder joints increases gradually at first and then drops rapidly with the increase of the thickness of IMC layer, reaching the maximum when the thickness of IMC layer is about 1.5-3.0 μm.
Key words: intermetallic compound; ball grid array; vibration; fatigue lifetime
Foundation item: Project (XJ201608) supported by the Special Project for Technological Innovation of Key Industries in Xiangcheng District, Suzhou, China; Project (2017) supported by the Industrial Technology Innovation Project and Industrial Joint Innovation Project in Suzhou, China; Project (2016-27) supported by the Science and Technology Project in Beibei District, Chongqing, China; Project (SKLMT-ZZKT- 2017M15) supported by the Key Laboratory Project of Mechanical Transmission in Chongqing University, China
Received date: 2018-03-19; Accepted date: 2018-05-31
Corresponding author: SHEN Jun; Tel: +86-13883111150; E-mail: shenjun@cqu.edu.cn
(编辑 龙怀中)
基金项目:苏州市相城区重点产业技术创新专项(XJ201608);苏州市产业技术创新专项-产业联合创新项目(2017);重庆市北碚区科技专项(2016-27);重庆大学机械传动重点实验室项目(SKLMT-ZZKT-2017M15)
收稿日期:2018-03-19;修订日期:2018-05-31
通信作者:沈 骏,教授,博士;电话:13883111150;E-mail:shenjun@cqu.edu.cn