基于有限元法的疊層薄壁工件自動化制孔壓緊力預(yù)測

劉雪鋒，朱偉東，楊國榮，黨曉娟

?

基于有限元法的疊層薄壁工件自動化制孔壓緊力預(yù)測

劉雪鋒1，朱偉東1，楊國榮2，黨曉娟2

(1. 浙江大學(xué) 機(jī)械工程學(xué)院，浙江 杭州，310027；2. 西安飛機(jī)工業(yè)(集團(tuán))有限責(zé)任公司，陜西 西安，710089)

針對實(shí)驗(yàn)法選擇合適壓緊力時間周期長、成本高的問題，分析疊層薄壁工件層間毛刺的產(chǎn)生原理，利用簡化的梁模型研究壓緊力對層間間隙的影響。根據(jù)層間毛刺和層間間隙的關(guān)系，融合有限元仿真和函數(shù)擬合，提出基于有限元的壓緊力理論預(yù)測方法，以實(shí)現(xiàn)制孔壓緊力的快速準(zhǔn)確預(yù)測。采用理論預(yù)測的壓緊力進(jìn)行2024疊層鋁板制孔實(shí)驗(yàn)。研究結(jié)果表明：構(gòu)建的層間毛刺高度模型預(yù)測誤差在10%以內(nèi)，實(shí)驗(yàn)孔的層間毛刺高度小于 0.1 mm，達(dá)到“無毛刺”的制孔要求，為自動化制孔壓緊力的選擇提供了理論參考。

自動化制孔；壓緊力；有限元；層間毛刺；層間間隙

在飛機(jī)制造領(lǐng)域，以鉚釘連接和螺栓連接為主的機(jī)械連接是飛機(jī)裝配的主要連接方式[1]，對連接孔的質(zhì)量要求較高，傳統(tǒng)的手工制孔逐漸被機(jī)器人自動化制孔所取代。飛機(jī)蒙皮、長珩、角片等疊層薄壁工件的自動化制孔，疊層工件之間會形成毛刺，嚴(yán)重影響裝配質(zhì)量[2?3]。工業(yè)上普遍通過壓緊裝置施加單側(cè)壓緊力，以抑制層間毛刺，提高制孔質(zhì)量。然而，單側(cè)壓緊力的存在，同樣會造成疊層工件局部受壓變形，降低制孔的精度。目前，單側(cè)壓緊力的選擇依靠經(jīng)驗(yàn)，針對不同的疊層結(jié)構(gòu)，需經(jīng)多次實(shí)驗(yàn)確定合適的壓緊力，缺乏有效的理論預(yù)測。國內(nèi)外學(xué)者對單層工件制孔的研究較多，研究主要針對剛度較大的材料，通過數(shù)理統(tǒng)計(jì)和實(shí)驗(yàn)相結(jié)合的方法，研究工藝參數(shù)、刀具材質(zhì)和幾何角度對制孔毛刺的影響[4?9]。近年來，部分學(xué)者利用有限元仿真進(jìn)行鉆削研究。SAUNDERS[10]建立二維有限元模型預(yù)測工件的溫度分布和應(yīng)力狀態(tài)。ISBILIR等[11]建立三維有限元鉆削模型，從軸向力、毛刺、應(yīng)力3個方面分析鉆削過程。對于飛機(jī)裝配領(lǐng)域的弱剛度疊層薄壁工件制孔問題，一直缺乏研究，相關(guān)文獻(xiàn)較少。LIANG[12]提出疊層弱剛度工件制孔的簡化模型，通過實(shí)驗(yàn)驗(yàn)證壓緊力對層間毛刺具有抑制作用。CHOI等[13]對疊層制孔過程進(jìn)行有限元仿真，分析層間間隙的形成以及層間間隙對制孔質(zhì)量的影響。王珉等[14]借助有限元分析層間貼合面間隙隨壓緊力的變化情況，得到最佳壓緊力的取值范圍。但對壓緊力的準(zhǔn)確預(yù)測研究較少。綜上所述，壓緊力對疊層薄壁工件制孔的質(zhì)量有著重要影響，但對壓緊力的精確選擇缺乏研究，尚無理論方法進(jìn)行準(zhǔn)確預(yù)測。目前，工業(yè)上采用實(shí)驗(yàn)法選擇合適的壓緊力，一旦制孔工件結(jié)構(gòu)發(fā)生變化，需要重新通過實(shí)驗(yàn)獲取最佳壓緊力，時間周期長、成本高。本文作者分析層間毛刺的生成過程，研究壓緊力對層間間隙的抑制作用，根據(jù)“無毛刺”的要求，提出基于有限元的壓緊力預(yù)測方法，構(gòu)建層間毛刺高度預(yù)測模型，能方便快捷地對壓緊力的取值進(jìn)行準(zhǔn)確預(yù)測，提高制孔質(zhì)量。

1 理論分析

1.1 毛刺的產(chǎn)生

在機(jī)械鉆削加工中，由于工件材料存在塑性變形，毛刺的產(chǎn)生無法避免。毛刺一般產(chǎn)生于刀具的入口和出口處，如圖1所示。毛刺的產(chǎn)生降低了工件表面加工質(zhì)量，產(chǎn)生應(yīng)力集中，降低了結(jié)構(gòu)件的疲勞強(qiáng)度，嚴(yán)重影響機(jī)械連接的使用壽命[15]。

層間毛刺在疊層工件制孔中產(chǎn)生，層間毛刺的存在，影響疊層工件的緊密貼合。相比單層工件的入口和出口毛刺，層間毛刺不易去除。傳統(tǒng)的制孔工藝需將制孔后的疊層薄壁工件拆開后再去除毛刺，然后重新裝配。據(jù)估計(jì)，毛刺的去除成本可占產(chǎn)品總成本的30%[16]，因而，抑制層間毛刺的生成是疊層薄壁工件制孔研究的重點(diǎn)。

圖1 毛刺的產(chǎn)生部位

層間毛刺的產(chǎn)生與層間間隙緊密相關(guān)。刀具鉆穿上層工件并鉆入下層工件，下層工件在軸向鉆削力和上層工件的共同作用下發(fā)生彈性變形，制孔附近區(qū)域形成層間間隙，給毛刺的生成提供了空間。上層工件的出口毛刺和下層工件的入口毛刺被擠壓進(jìn)入層間間隙，形成層間毛刺。

1.2 層間毛刺和層間間隙

層間毛刺的高度是指上層工件出口毛刺和下層工件入口毛刺的高度總和，假設(shè)上層工件出口毛刺高度為1，下層工件入口毛刺高度為2。根據(jù)層間毛刺生成時，層間間隙對層間毛刺的生長是否起到抑制作用，可將層間毛刺高度和層間間隙的關(guān)系分為2種情況，如圖2所示。

(a) 層間間隙抑制層間毛刺；(b) 層間間隙無法抑制層間毛刺

當(dāng)層間間隙較小時，層間間隙對層間毛刺的生成有抑制作用，層間出口毛刺和入口毛刺緊密貼合， 此時，

1+2=(1)

當(dāng)層間間隙較大時，層間間隙無法抑制層間毛刺的生成，層間出口毛刺和層間入口毛刺之間留有空隙，此時，

1+2＜(2)

疊層薄壁工件制孔一般要求層間緊密貼合，抑制層間間隙的產(chǎn)生，因而，制孔時層間間隙較小。本文研究層間間隙較小的情況。

1.3 壓緊力和層間間隙

式中：X，Y和C均為與工件材料有關(guān)的系數(shù)；為刀具直徑；為進(jìn)給速度；K為修正系數(shù)。

由式(3)可知：當(dāng)工件的材料和制孔刀具確定后，制孔時產(chǎn)生的軸向鉆削力只與進(jìn)給速度有關(guān)，與壓緊力無關(guān)。

圖3 簡化的疊層薄壁工件制孔模型

當(dāng)?shù)毒咩@穿上層工件，鉆入下層工件時，軸向鉆削力作用在下層工件上，此時疊層工件層間間隙最大。假設(shè)初始的層間間隙為0 mm，不斷增大壓緊力，制孔過程中的最大層間間隙的變化如圖4所示。圖4(a)所示為沒有壓緊力作用的情況，即=0 N，下層工件在鉆削力1的作用下發(fā)生變形，形成層間間隙，層間間隙較大。圖4(b)所示為通過壓腳施加單側(cè)壓緊力的情況，即＞0 N，上層工件和下層工件均發(fā)生彈性變形，層間間隙減小。圖4(c)所示為壓緊力繼續(xù)增大的情況，上層工件和下層工件固支端附近開始緊密貼合，層間間隙繼續(xù)減小。當(dāng)貼合面延伸至壓腳兩側(cè)時，壓緊力達(dá)到臨界值，繼續(xù)增加壓緊力幾乎對層間間隙無任何影響，如圖4(d)所示。疊層工件在制孔過程中，由于鉆削力1的存在，導(dǎo)致層間間隙無法完全消除，但可以通過施加單側(cè)壓緊力來有效地減小層間間隙。

(a) F=0 N；(b) F增大；(c) F繼續(xù)增大；(d) F達(dá)到臨界值

2 有限元仿真

疊層薄壁工件的種類繁多，結(jié)構(gòu)各異，通過梁模型可以定性地分析施加壓緊力對疊層間隙的影響，但無法進(jìn)行精確的理論計(jì)算。對于不同結(jié)構(gòu)的疊層薄壁工件，均可采用有限元法對壓緊力和層間間隙的關(guān)系進(jìn)行進(jìn)一步研究。

疊層工件的鉆削仿真涉及熱力耦合[18]，需要考慮工件之間的接觸和摩擦，仿真過程復(fù)雜，且花費(fèi)時間較長。為達(dá)到快速預(yù)測的目的，提出等效鉆削仿真的概念，即忽略鉆削過程中工件的振蕩，通過施加靜態(tài)鉆削力模擬毛刺成型時刻工件的層間間隙。

在Solidworks中建立簡化的等效鉆削模型，并導(dǎo)入Abaqus軟件進(jìn)行有限元仿真，如圖5所示。模型由支座、疊層鋁板、壓腳構(gòu)成，通過約束疊層鋁板4個螺栓孔內(nèi)表面的自由度，模擬螺栓連接，將疊層鋁板固定在支座上。表1所示為疊層鋁板參數(shù)[19]，疊層鋁板長×寬×高為200 mm×140 mm×2 mm，壓腳內(nèi)徑1=17.5 mm，外徑2=22.5 mm，厚度Y=2 mm。預(yù)先通過鉆削實(shí)驗(yàn)，利用Kistler9257B型測力儀測得相同實(shí)驗(yàn)條件下，當(dāng)主軸轉(zhuǎn)速=3 000 r/min，進(jìn)給速度=180 mm/min時，鉆削力1=140 N。

設(shè)置支座和壓腳為剛體，上、下鋁板為可變形三維實(shí)體，網(wǎng)格類型為C3D8R。有限元仿真采用準(zhǔn)靜態(tài)仿真，分為2個分析步，第1步在壓腳上表面施加壓緊力，第2步在下層鋁板的制孔位置施加靜態(tài)鉆削力。

首先，希望大家可以進(jìn)一步凝聚共識、明確目標(biāo)。牢固樹立共建共治共享的理念，著眼于促進(jìn)西江經(jīng)濟(jì)帶發(fā)展和保護(hù)西江母親河的共同愿景，充分發(fā)揮共建各方的優(yōu)勢，形成強(qiáng)大監(jiān)管合力，進(jìn)一步夯實(shí)西江治理能力和治理體系現(xiàn)代化管理基礎(chǔ)，規(guī)范西江航運(yùn)市場秩序，合力打擊違法違規(guī)經(jīng)營行為，提高安全生產(chǎn)水平，構(gòu)建“安全命運(yùn)共同體”。

表1 疊層薄壁工件參數(shù)

(a) 視圖1；(b) 視圖2

在上層鋁板的下表面和下層鋁板的上表面對應(yīng)的制孔邊緣位置分別設(shè)置參考點(diǎn)，通過觀察參考點(diǎn)在壓緊力方向上的位移差，研究層間間隙的變化。依次設(shè)置0，100，200，300，400，500和600 N的壓緊力，進(jìn)行有限元仿真計(jì)算。在Abaqus后處理中，測量得不同壓緊力作用下的層間間隙，如表2所示。由表2可見：當(dāng)壓緊力小于200 N時，層間間隙高度隨壓緊力的增大迅速減??；當(dāng)壓緊力在200~600 N時，層間間隙的高度隨壓緊力的增大緩慢減小，層間間隙趨于穩(wěn)定，壓緊力對層間間隙的影響速率逐漸減弱，與理論分析結(jié)果相吻合。

表2 壓緊力作用下的層間間隙

3 壓緊力預(yù)測

根據(jù)理論分析，當(dāng)疊層薄壁工件層間間隙較小時，層間間隙對層間毛刺的生成具有抑制作用，此時，層間間隙高度可代表層間毛刺高度。基于層間間隙和層間毛刺的關(guān)系，提出基于有限元的壓緊力預(yù)測方法。即通過有限元仿真獲取壓緊力和層間間隙的關(guān)系，并進(jìn)行函數(shù)擬合，建立層間毛刺高度預(yù)測模型，結(jié)合制孔質(zhì)量要求進(jìn)行壓緊力預(yù)測。

針對擬合函數(shù)，常用判定系數(shù)2和均方根誤差RMS(root mean square error)進(jìn)行適應(yīng)性檢驗(yàn)，驗(yàn)證擬合效果。判定系數(shù)2和均方根誤差RMS的求解 公式[20]：

判定系數(shù)2取值為0~1，2越大，RMS越小，代表函數(shù)擬合效果越好，對模型的解釋能力越強(qiáng)。根據(jù)不同壓緊力作用下的有限元仿真結(jié)果，對層間間隙和壓緊力的關(guān)系進(jìn)行函數(shù)擬合，建立層間毛刺高度預(yù)測模型，結(jié)果如下式所示：

式中：的取值范圍為0~600 N。擬合函數(shù)的判定系數(shù)2=0.982，均方根誤差RMS=9.927，表明擬合效果較好。圖6所示為有限元仿真結(jié)果擬合曲線。

由圖6可知：壓緊力對層間間隙具有抑制作用，從而降低層間毛刺高度，提高制孔質(zhì)量。但壓緊力的增大會使局部變形增大，導(dǎo)致制孔精度下降。所以，在滿足毛刺高度要求的情況下，盡可能選擇較小的壓緊力。在飛機(jī)制造領(lǐng)域，根據(jù)供應(yīng)商的不同，“無毛刺”的要求為毛刺高度控制在0.10~0.15 mm[1]?，F(xiàn)取“無毛刺”的高度閥值為0.10 mm，結(jié)合擬合函數(shù)式(6)，令

1—仿真結(jié)果；2—擬合曲線。

根據(jù)式(7)，可求得壓緊力的最小值=123 N。由于加工、裝夾等因素存在誤差，為保證層間毛刺高度低于0.10 mm，故選取安全系數(shù)=1.1，在實(shí)際加工中，壓緊力exp為

在實(shí)際制孔中，選擇壓緊力exp135 N，可滿足工程制孔要求。

4 實(shí)驗(yàn)驗(yàn)證

為驗(yàn)證基于有限元的壓緊力理論預(yù)測方法的可行性，采用2024疊層鋁板進(jìn)行鉆削制孔實(shí)驗(yàn)，疊層鋁板長×寬×高為200 mm×140 mm×2 mm，刀具采用硬質(zhì)合金鉆頭，直徑=5.1 mm，頂角=120°，螺旋升角=30°。疊層鋁板通過螺栓連接固定在支座上，和有限元模型保持一致。支座后端連接Kistler9257B型測力儀，實(shí)驗(yàn)在自動化制孔平臺上進(jìn)行，如圖7所示。手動調(diào)節(jié)氣壓閥，調(diào)節(jié)壓緊力=135 N，設(shè)置轉(zhuǎn)速= 3 000 r/min，進(jìn)給速度=180 mm/min。保持壓緊力不變，在疊層鋁板的中心位置進(jìn)行鉆削實(shí)驗(yàn)，實(shí)驗(yàn)共鉆削4個實(shí)驗(yàn)孔。

利用Dino-Lite顯微鏡觀察實(shí)驗(yàn)孔的層間毛刺情況，結(jié)果如圖8所示。由圖8可見：層間毛刺集中在層間出口處，層間入口處的毛刺很小，肉眼無法直接辨識。

圖7 自動化制孔實(shí)驗(yàn)平臺

利用百分表測量孔1~孔4的層間毛刺高度，測得實(shí)驗(yàn)孔的毛刺高度為0.088~0.098 mm，均小于“無毛刺”的閥值0.10 mm，如圖9所示。通過層間毛刺高度預(yù)測模型，求得當(dāng)=135 N時，毛刺的理論預(yù)測高度=0.091 mm。毛刺的理論預(yù)測高度和實(shí)際測量高度基本一致，相對誤差在10%以內(nèi)，表明毛刺高度預(yù)測模型具有精確的預(yù)測效果。

實(shí)驗(yàn)采用計(jì)算得到的壓緊力，能夠控制毛刺高度低于0.10 mm，達(dá)到“無毛刺”的制孔要求，滿足工程制孔要求，表明基于有限元的壓緊力預(yù)測方法具有可行性。

(a) 孔1出口毛刺；(b) 孔1入口毛刺；(c) 孔3出口毛刺；(d) 孔3入口毛刺

1—實(shí)驗(yàn)值；2—理論值。

5 結(jié)論

1) 壓緊力對工件層間間隙具有抑制作用，隨著壓緊力的增大，抑制速率逐漸減弱，層間間隙逐漸減小，最終趨于穩(wěn)定。施加壓緊力能有效地減小層間毛刺，但無法完全消除層間毛刺。

2) 建立的層間毛刺高度預(yù)測模型具有準(zhǔn)確的預(yù)測效果，滿足“無毛刺”的制孔要求，表明基于有限元的壓緊力預(yù)測方法具有可行性，為自動化制孔壓緊力的精確選擇提供了理論參考。

[1] 王珉, 薛少丁, 陳文亮, 等. 面向飛機(jī)自動化裝配的單向壓緊制孔毛刺控制技術(shù)[J]. 航空制造技術(shù), 2011, 381(9): 26?29. WANG Min, XUE Shaoding, CHEN Wenliang, et al. Burr control technology of one-side clamping force for aircraft automatic assembly[J]. Aeronautical Manufacturing Technology, 2011, 381(9): 26?29.

[2] AURICH J C, DORNFELD D, ARRAZOLA P J, et al. Burrs-analysis, control and removal[J]. CIRP Annals Manufacturing Technology, 2009, 58(2): 519?542.

[3] TIAN W, HU J, LIAO W, et al. Formation of interlayer gap and control of interlayer burr in dry drilling of stacked aluminum alloy plates[J]. Chinese Journal of Aeronautics, 2016, 29(1): 283?291.

[4] KARNIK S R, GAITONDE V N, DAVIM J P. A comparative study of the ANN and RSM modeling approaches for predicting burr size in drilling[J]. International Journal of Advanced Manufacturing Technology, 2008, 38(9/10): 868?883.

[5] SAUNDERS L K L, MAUCH C A. An exit burr model for drilling of metals[J]. Journal of Manufacturing Science and Engineering, 2001, 123(4): 562?566.

[6] BU Y, LIAO W H, TIAN W, et al. An analytical model for exit burrs in drilling of aluminum materials[J]. The International Journal of Advanced Manufacturing Technology, 2016, 85(9/10/11/12): 2783?2796.

[7] GAITONDE V N, KARNIK S R. Minimizing burr size in drilling using artificial neural network (ANN)-particle swarm optimization (PSO) approach[J]. Journal of Intelligent Manufacturing, 2012, 23(5): 1783?1793.

[8] KILICKAP E. Modeling and optimization of burr height in drilling of Al-7075 using Taguchi method and response surface methodology[J]. The International Journal of Advanced Manufacturing Technology, 2010, 49(9/10/11/12): 911?923.

[9] MELKOTE S N, NEWTON T R, HELLSTERN C, et al. Interfacial burr formation in drilling of stacked aerospace materials[C]//Burrs-Analysis, Control and Removal. Berlin: Springer, 2010: 89?98.

[10] SAUNDERS L K L. A finite element model of exit burrs for drilling of metals[J]. Finite Elements in Analysis and Design, 2003, 40(2): 139?158.

[11] ISBILIR O, GHASSEMIEH E. Finite element analysis of drilling of titanium alloy[J]. Procedia Engineering, 2011, 10(7): 1877?1882.

[12] LIANG Jie. The formation and effect of interlayer gap in dry drilling of stacked metal materials[J]. The International Journal of Advanced Manufacturing Technology, 2013, 69(5/6/7/8): 1263?1272.

[13] CHOI J, MIN S, DORNFELD D, et al. Modeling of inter-layer gap formation in drilling of a multi-layered material[C]// Proceedings of the 6th CIRP International Workshop on Modeling of Machining. Hamilton: McMaster University, 2003: 19?20.

[14] 王珉, 薛少丁, 蔣宏宇, 等. 飛機(jī)大部件對接自動化制孔單向壓緊力分析[J]. 南京航空航天大學(xué)學(xué)報, 2012, 44(4): 553?558. WANG Min, XUE Shaoding, JIANG Hongyu, et al. One-side pressure-force analysis of automatic drilling of aircraft fuselage section-joint assembly[J]. Journal of Nanjing University of Aeronautics & Astronautics, 2012, 44(4): 553?558.

[15] 卜泳, 許國康, 肖慶東. 飛機(jī)結(jié)構(gòu)件的自動化精密制孔技術(shù)[J]. 航空制造技術(shù), 2009(24): 61?64. BU Yong, XU Guokang, XIAO Qingdong. Automatic precision drilling technology of aircraft structural part[J]. Aeronautical Manufacturing Technology, 2009(24): 61?64.

[16] GILLESPIE L K. Deburring precision miniature parts[J]. Precision Engineering, 1979, 1(4): 189?198.

[17] 洪華舟, 韋紅余, 陳文亮, 等. 航空薄壁件制孔毛刺生長控制工藝研究[J]. 中國機(jī)械工程, 2012, 23(19): 2312?2316. HONG Huazhou, WEI Hongyu, CHEN Wenliang, et al. Study on the burr growth control technology of aeronautical thin-walled workpiece[J]. China Mechanical Engineering, 2012, 23(19): 2312?2316.

[18] 莊茁, 由小川, 廖劍暉, 等. 基于ABAQUS的有限元分析和應(yīng)用[M]. 北京: 清華大學(xué)出版社, 2009: 287?350. ZHUANG Zhuo, YOU Xiaochuan, LIAO Jianhui, et al. The finite element analysis and application based on ABAQUS[M]. Beijing: TsingHua University Press, 2009: 287?350.

[19] BUYUK M, KAN S, LOIKKANEN M J. Explicit finite-element analysis of 2024-T3/T351 aluminum material under impact loading for airplane engine containment and fragment shielding[J]. Journal of Aerospace Engineering, 2009, 22(3): 287?295.

[20] MONTGOMERY D C. 試驗(yàn)設(shè)計(jì)與分析[M]. 汪仁官, 陳榮昭, 譯. 北京: 中國統(tǒng)計(jì)出版社, 1998: 563?564. MONTGOMERY D C. Design and analysis of experiments[M]. WANG Renguan, CHEN Rongzhao, trans. Beijing: China Statistics Press, 1998: 563?564.

(編輯 劉錦偉)

Clamping force prediction for automatic drilling of stacked thin-walled workpieces based on finite element method

LIU Xuefeng1, ZHU Weidong1, YANG Guorong2, DANG Xiaojuan2

(1. Department of Mechanical Engineering, Zhejiang University, Hangzhou 310027, China; 2. Xi’an Aircraft Industry(Group) Ltd. Co., Xi’an 710089, China)

To solve the problem of the long time and high cost of clamping force selection by experimental method, the interlayer burr forming theory of thin-walled workpieces was analyzed and the influences of clamping force on the interlayer gap by a simplified beam model were explained. According to the relationship between interlayer burrs and the interlayer gap, a clamping force prediction method was proposed, which was based on the finite element simulation and function fitting to achieve a fast and accurate prediction for clamping force. The predicted clamping force was used in 2024 laminated aluminum plate drilling experiments. The results show that the error of the prediction model of the conducted interlayer burr height is less than 10%, and the interlayer burr height of experimental hole is less than 0.1 mm, which can meet the quality requirements of drilled holes and provide a theoretical reference for clamping force selection.

automatic drilling; clamping force; finite element; interlayer burr; interlayer gap

10.11817/j.issn.1672?7207.2018.02.011

V262.4

A

1672?7207(2018)02?0339?06

2017?03?16；

2017?04?27

國家自然科學(xué)基金資助項(xiàng)目(51205352)(Project(51205352) supported by the National Natural Science Foundation of China)

朱偉東，博士(后)，副教授，從事機(jī)器人制孔和航空復(fù)材研究；E-mail：wdzhu@zju.edu.cn