蘇俊省,王君杰,郭 進,彭運動(.同濟大學土木工程學院,上海0009;.石家莊鐵道大學土木工程學院,河北石家莊05004;.中交公路規(guī)劃設計院有限公司,北京00088)
?
基于鋼筋低周疲勞的橋墩地震易損性分析
蘇俊省1,王君杰1,郭 進2,彭運動3
(1.同濟大學土木工程學院,上海200092;2.石家莊鐵道大學土木工程學院,河北石家莊050043;3.中交公路規(guī)劃設計院有限公司,北京100088)
基于Coffin- Manson鋼筋疲勞損傷模型,采用非線性纖維梁柱單元,對縱筋采用HRB335和HRB500E的圓柱形橋墩擬靜力試驗進行數值模擬,研究擬靜力作用下低周疲勞對鋼筋和構件承載力退化的影響.為了進一步研究鋼筋低周疲勞對試件累積損傷的影響,將Takemura和Kunnath提出的兩組不同加載模式下的擬靜力試驗進行數值模擬.結果表明:基于Coffin- Manson模型,采用纖維單元,在材料層面上考慮鋼筋的低周疲勞,可以較好地模擬試件在不同加載模式下的累積損傷和承載力退化.
鋼筋混凝土;橋墩;低周疲勞;累積損傷;承載力退化;纖維單元
準確預測結構在地震作用下的損傷程度是結構基于性能抗震設計的關鍵.地震作用具有持時長、往復振動次數多的特點,因此地震尤其強震作用下,結構會產生較大的累計損傷,導致承載力降低,試件破壞,進而引發(fā)整體結構倒塌破壞.因此,為了更好地模擬結構構件在地震作用下的響應,需要在結構非線性動力分析時考慮累積損傷導致的承載力降低.
反復荷載作用下,鋼筋低周疲勞、混凝土保護層剝落和黏結滑移的劣化引起鋼筋混凝土構件承載力退化.Park等[1]指出,鋼筋混凝土構件的累積損傷與試件的最大位移和累積滯回耗能相關.文獻[2-4]提出了基于最大位移的承載力退化模型,文獻[5-7]提出了基于滯回耗能的承載力退化模型.郭進等[8]在彎矩-曲率層面提出了以損傷值為變量的強度退化模型.曲哲等[9]基于力-位移層面的鋼筋混凝土構件承載力退化模型,提出了在纖維模型桿單元中考慮鋼筋混凝土構件的承載力退化,并將構件承載力的退化全部計入鋼筋纖維的滯回本構中.
地震損傷破壞主要包括材料、構件和結構整體[10]3個層面,其中材料損傷是結構損傷的本質,而現有的承載力退化模型均為構件或結構層面上的承載力退化模型.因此,基于材料損傷層面建立強度退化模型,通過材料層面的強度退化反映結構構件上的強度退化,能夠更直觀地反映試件損傷破壞的本質.
Kunnath等[11]和劉伯權等[12]的鋼筋混凝土柱低周疲勞試驗表明,縱筋的低周疲勞損傷是導致鋼筋混凝土試件累積損傷的重要因素.劉鳴等[13]通過對鋼筋混凝土柱低周疲勞試驗研究表明,鋼筋混凝土試件的強度退化與構件的損傷過程發(fā)展一致.因此,鋼筋的低周疲勞損傷與試件的強度退化直接相關,在鋼筋本構模型中必須考慮低周疲勞,建立低周疲勞損傷與強度退化之間的關系,在材料層面上解決由累積損傷引起的強度退化問題.
結構彈塑性分析模擬地震作用下結構動力響應時,纖維單元模型與梁柱單元和實體單元相比,能以較低的計算成本獲得較高的求解精度[14-15].本文采用纖維模型,考慮鋼筋材料的低周疲勞性能,對2個圓柱形橋墩進行數值模擬,研究鋼筋的低周疲勞性能對試件抗震性能的影響.通過模擬Takemura和Kunnath 2組采用不同加載模式的試驗,進一步研究了在材料層面上考慮鋼筋低周疲勞時,材料損傷與強度退化之間的關系,可以較好地模擬出不同加載模式下,累積損傷對試件強度退化的影響.
地震作用下,鋼筋的低周疲勞對試件的抗震性能有較大影響.鋼筋的低周疲勞模型主要包括以下3種:基于材料應變幅值和加載周數的Coffin -Manson疲勞模型[16-17],考慮鋼筋應力、應變和加載周數的Smith-Watson-Topper(SWT)疲勞模型[18]和修正的SWT疲勞模型[19-20],基于耗能的鋼筋疲勞模型[21-23].由于Coffin- Manson模型形式簡單,經過大量試驗驗證能較好地模擬鋼筋在大應變下的低周疲勞性能,因此本文采用Coffin- Manson疲勞模型.
Coffin[16]和Manson[17]首先提出了基于材料應變幅值和加載周數的疲勞壽命方程,如下所示;
式中:εp為塑性應變;Cf為疲勞試驗測定的材料疲勞常數;α為疲勞指數;2 Nf為按εp進行等幅值加載直到鋼筋斷裂時可完成的半周循環(huán)加載次數.εp根據下式計算可得:
式中:σt和εt分別為當前加載循環(huán)中材料的總應力和總應變(見圖1a);Es為鋼筋的彈性模量.
Brown等[24]根據鋼筋的ε-2 Nf低周疲勞公式提出,當鋼筋應變幅值為εpi時,半周循環(huán)引起的鋼筋損傷Di為加載周數的倒數(見式(3)),累積損傷D為Di之和,計算公式如下所示:D=0時表示鋼筋無損傷,D=1時表示鋼筋斷裂.鋼筋斷裂后,強度迅速退化為零.
圖1 鋼筋疲勞損傷模型Fig.1 Fatigue damage model for steel bars
試件循環(huán)加載過程中,鋼筋疲勞損傷導致鋼筋的強度降低,如圖1b所示.其中鋼筋強度降低系數φSR與損傷指數D線性相關,如下所示:
其中,Cd為強度降低常數.因此,鋼筋強度降低系數φSR可以簡化為
2.1 纖維單元
本文采用纖維模型進行數值模擬,有限單元剛度法在梁柱單元非線性分析中程序較為簡單,收斂性較好,因此本文采用剛度法非線性梁柱單元.每個試件沿高度劃分為9個纖維梁柱單元,每個單元取5個積分點.梁柱單元截面采用纖維單元,各構件均將截面離散為核心混凝土、保護層混凝土和鋼筋.試件的有限元模型如圖2所示.數值模擬時采用OpenSees有限元軟件,核心區(qū)混凝土采用Kent-Park約束混凝土模型,保護層混凝土采用Kent-Park非約束混凝土模型[25](由Concrete02 Material模型模擬).鋼筋模型采用修正的Chang- Mander鋼筋模型[26],并考慮鋼筋的屈曲和低周疲勞[27-28](由ReinforcingSteelMaterial模型模擬).
圖2 鋼筋混凝土橋墩有限元模型Fig.2 Finite element model of reinforced concrete model
2.2 鋼筋材料
Reinforcing Steel Material鋼筋本構[29]模型是基于Chang- Mander鋼筋本構[26]發(fā)展而來的,鋼筋的應力應變曲線如圖3所示.Chang- Mander鋼筋模型能較好地模擬鋼筋的骨架曲線形狀,包括彈性階段、屈服平臺段、應變硬化段和應變軟化段.為了更好地模擬鋼筋的受壓屈曲和低周疲勞特性,Kunnath等[28]對該模型進行了改進,考慮了鋼筋的屈曲和低周疲勞,并將該改進后的Reinforcing Steel Material本構模型加入到OpenSees有限元軟件中[29].
鋼筋本構模型如圖3所示.圖3a表示鋼筋的骨架曲線.其中,fy為鋼筋屈服應力,fu為鋼筋極限應力;εsu為鋼筋極限應力對應的應變,εsh為鋼筋開始發(fā)生硬化時對應的應變;Es為鋼筋初始彈性模量,Esh為鋼筋硬化彈性模量.圖3b表示鋼筋的滯回規(guī)則.其中,①、②表示骨架段,③、④表示反向段,⑤、⑥表示返回段,⑦、⑧表示第1過渡段,⑨、⑩表示第2過渡段.
圖3 OpenSees Reinforcing Steel Material模型應力應變關系Fig.3 Stress-strain relationship for OpenSeesReinforcing Steel Material model
Concrete02Material本構應力應變關系如圖4所示,該模型是基于Kent-Park本構發(fā)展而來,混凝土拉壓側軟化段均為線性下降.其中,fpc和εpc分別為混凝土最大壓應力和對應的應變;fpu和εpu分別為混凝土壓碎時的應力和對應的應變;ft和εt分別為混凝土極限拉應力和對應的應變;E0為混凝土初始彈性模量,計算如式(8)所示;λ表示混凝土為壓應變εpu時卸載斜率與初始剛度的比值;Et為混凝土受拉卸載剛度.
2.3 混凝土材料
圖4 OpenSees Concrete02Material模型應力應變關系Fig.4 Stress-strain relationship for OpenSeesConcrete02Material model
3.1 試驗介紹
本文進行了2個圓柱形橋墩的試驗,試件的尺寸與配筋構造完全相同,所采用鋼筋類型不同(見圖5).試件C -LM-C40和C -MM-C40縱筋分別采用HRB335和HRB500E鋼筋,箍筋均采用HRB500E鋼筋,鋼筋的力學性能見表1.試件為鋼筋混凝土懸臂柱,截面直徑600mm,保護層厚度25mm,截面有效高度550mm,加載中心至墩底高度2 650mm,剪跨比4.82;采用C40混凝土,混凝土立方體抗壓強度標準值為42.9MPa;縱筋采用D20鋼筋,配筋率1.56%,箍筋采用D10鋼筋,箍筋間距60mm,體積配箍率為0.8%.試驗時柱頂施加760kN的軸壓力,軸壓比為0.083.保持軸力恒定,然后施加水平荷載,水平荷載通過位移控制,加載位移0~10mm范圍內位移增量為2mm,10~50mm位移加載段位移增量為5mm,50mm之后位移增量為10mm;每級位移幅值下進行3次循環(huán)加載,每級加載后,進行1次上一級位移幅值的加載.試驗加載規(guī)則見圖6所示.
圖5 圓柱形試件尺寸和配筋圖(單位:mm)Fig.5 Dimension and reinforcements of circularspecimens(unit:mm)
表1 鋼筋材料性能指標Tab.1_Material properties of reinforcement
圖6 加載規(guī)則曲線Fig.6 Loading time history
3.2 數值模擬
鋼筋混凝土柱鋼筋和混凝土材料分別采用上文中介紹的Reinforcing Steel Material和Concrete02 Material模型進行模擬,材料的力學性能如表2和3所示.
表2 數值模擬中混凝土力學特性Tab.2 Mechanical properties of concrete material in numerical analysis
試件采用不同強度的鋼筋,因此在考慮鋼筋低周疲勞時,低周疲勞參數略有差別.試驗中得到的混凝土強度為立方體抗壓強度標準值,Concrete02 Material模型中為圓柱體抗壓強度,因此根據文獻[30]中規(guī)定,采用立方體抗壓強度乘以0.85獲得圓柱體抗壓強度.
表3 數值模擬中鋼筋力學特性Tab.3 Mechanical properties of reinforcing steel material in numeric alanalysis
為了研究縱筋低周疲勞對試件數值模擬結果的影響,采用了考慮和不考慮縱筋低周疲勞2種數值模擬方法,數值模擬與試驗滯回曲線比較如圖7所示.為了研究試驗過程中鋼筋的低周疲勞損傷性能,提取了數值模擬時圓形截面最外層縱筋的應力應變曲線,如圖8所示.根據式(4)計算得到的最外層鋼筋損傷指數,如圖9所示.
圖7 試驗與數值模擬結果比較Fig.7 Comparison between experimental and analysis results
由圖9可知,圓柱形試件2與圓柱形試件1相比,隨著縱筋強度由HRB335變?yōu)镠RB500E,鋼筋低周疲勞性能提高,疲勞常數Cf由0.6變?yōu)?.8,加載時鋼筋損傷值增長變慢,試件的延性性能提高.
加載初期,鋼筋疲勞損傷值較小,試件損傷由表層混凝土開裂和脫落引起;加載后期,隨著鋼筋應變增大,鋼筋低周疲勞上升逐漸增大,最后導致鋼筋損傷斷裂,試件破壞.
圖8 數值模擬鋼筋應力應變曲線Fig.8 Stress-strain curves of reinforcing steel based on numerical analysis
圖9 數值模擬鋼筋損傷Fig.9 Damage of reinforcing steel based on numerical analysis
由圖8可知,考慮鋼筋低周疲勞時,鋼筋應力應變曲線隨著疲勞指數增大出現明顯的強度降低,當疲勞指數為1時,鋼筋強度退化為零.結合圖7中力-位移骨架曲線,當考慮鋼筋低周疲勞時,鋼筋強度明顯降低,從而可以較好地模擬鋼筋混凝土橋墩強度退化.考慮鋼筋低周疲勞時,采用纖維模型可以較好地模擬出鋼筋混凝土橋墩的力-位移滯回曲線.
為了進一步研究鋼筋低周疲勞對試件累積損傷的影響,本文對Takemura等[31]試驗的1組采用不同加載模式的6個鋼筋混凝土柱進行數值模擬.6個試件的尺寸與配筋構造均完全相同,采用了6種不同的加載制度,對應的構件編號分別為TP001~TP006.
圖10 數值模擬結果與Takemura試驗結果對比Fig.10 Comparison between numerical analysis and experimental results of Takemura test
數值模擬時采用了考慮鋼筋低周疲勞和不考慮鋼筋低周疲勞2種方式,結果如圖10所示.通過考慮與不考慮鋼筋低周疲勞數值模擬結果對比發(fā)現,考慮鋼筋低周疲勞時,數值模擬可以較好地模擬不同加載模式下鋼筋混凝土柱的累積損傷導致的強度退化,從而較好地模擬鋼筋混凝土柱的力-位移滯回曲線.
Kunnath等[28]完成的1組試驗,包括以等幅值等加載模式模擬地震波.通過Kunnath試驗的數值模擬,進一步驗證鋼筋的低周疲勞損傷模型,以及考慮鋼筋低周疲勞損傷時,采用纖維單元對不同加載模式下累積損傷的模擬效果.12個試件的尺寸與配筋構造均完全相同,僅加載制度不同,12種不同加載模式對應的構件編號分別為A01~A12.
圖11 數值模擬結果與Kunnath試驗結果對比(試樣A01~A08)Fig.11 Comparison between numerical analysis and experimental results of Kunnath test(A01~A08)
試驗數值模擬采用了考慮和不考慮鋼筋低周疲勞2種方式,結果如圖11和12所示.由圖11d~f可知,在相同加載幅值時,采用Coffin- Manson鋼筋疲勞模型,可以較好地模擬鋼筋混凝土試件的低周疲勞損傷破壞.由圖11g~h、圖12a~d可知,考慮鋼筋低周疲勞時,纖維模型能更好地模擬試件累積損傷導致的承載力退化.因此,采用纖維單元時,可以在材料層面考慮鋼筋的損傷,并且反映到試件的累積損傷,較準確地模擬鋼筋混凝土柱的力-位移滯回曲線.
圖12 數值模擬結果與Kunnath試驗結果對比(試樣A09~A12)Fig.12 Comparison between numerical analysis and experimental results of Kunnath test(A09~A12)
(1)基于纖維單元,采用Coffin- Manson鋼筋疲勞模型,可以較好地模擬鋼筋的強度退化以及由鋼筋強度退化導致的鋼筋混凝土試件強度退化.
(2)試件加載初期,鋼筋損傷值較小,鋼筋混凝土橋墩損傷主要由表層混凝土開裂和脫落引起;隨著加載位移增大,鋼筋疲勞損傷增大,試件承載力降低;最后鋼筋損傷斷裂,試件破壞.
(3)通過縱筋采用HRB500E與HRB335的鋼筋混凝土柱對比可知,高強度鋼筋的低周疲勞性能更好,因此采用高強度鋼筋的試件具有更好的承載力和變形能力.
(4)通過對Takemura和Kunnath 2組采用不同加載模式的試驗進行數值模擬可知,采用纖維單元,鋼筋低周疲勞可以較好地模擬試件在等幅值和變幅值等不同加載模式下的強度退化和試件破壞.因此,通過考慮鋼筋的低周疲勞,可以建立起材料層面上的損傷模型,從而基于材料損傷模擬鋼筋混凝土試件的累積損傷和承載力退化.
本文僅考慮了鋼筋材料的低周疲勞損傷,未考慮混凝土的損傷,結果與結論適用于低軸壓比下以鋼筋低周疲勞破壞為主的鋼筋混凝土試件,對高軸壓比下混凝土試件有待進一步驗證.
[1]Park Y J,Ang A H S.Mechanistic seismic damage model for reinforced concrete[J].Journal of Structural Engineering,ASCE,1985,111(4):722.
[2]Roufaiel M S L,Meyer C.Analytical modeling of hysteretic behavior of R/C frames[J].Journal of Structural Engineering,ASCE,1987,113(3):429.
[3]Chung Y S,Meyer C,Masanobu S.Modeling of concrete damage[J].ACI Structural Journal,1989,86(3):259.
[4]Lai S S,Will G T,Otani S.Model for inelastic biaxial bending of concrete members[J].Journal of Structural Engineering,ASCE,1984,110(11):2563.
[5]Youssef M,Ghobarah A.Strength deterioration due to bond slip and concrete crushing in modeling of reinforced concrete members[J].ACI Structural Journal,1999,96(6):956.
[6]Kunnath S K,Reinhorn A M,Park Y J.Analytical modeling of inelastic seismic response of R/C structures[J].Journal of Structural Engineering,ASCE,1990,116(4):996.
[7]Mork K J.Response analysis of reinforced concrete structures under seismic excitation[J].Earthquake Engineering and Structural Dynamics,1991,23(1):33.
[8]郭進,王君杰.基于累積損傷的彎矩-曲率強度退化滯回模型[J].同濟大學學報:自然科學版,2013,41(12):1767.GUO Jin,WANG Junjie.Reinforced concrete strength degradation hysteretic moment-curvature model based on cumulative damage[J].Journal of Tongji University:Natural Sceince,2013,41(12):1767.
[9]曲哲,葉列平.基于有效累積滯回耗能的鋼筋混凝土構件承載力退化模型[J].工程力學,2011,28(6):45.QU Zhe,YE Lieping.Strength deterioration model based on effective hysteretic energy dissipation for RC members undercyclic loading[J].Engineering Mechanics,2011,28(6):45.
[10]連業(yè)達,王憲杰,張洵安,等.結構地震累積損傷指標研究進展[J].工業(yè)建筑,2012,42(4):118.LIAN Yeda,WANG Xianjie,ZHANG Xunan,et al.Research advances of structural seismic cumulative damage index[J].Industrial Construction,2012,42(4):118.
[11]Kunnath S K,El-bahy A,Taylor A W,et al.Cumulative seismic damage of reinforced concrete bridge piers[R].Buffalo:National Center for Earthquake Engineering Research,1997.
[12]劉伯權,白紹良,徐云中,等.鋼筋混凝土柱低周疲勞性能的試驗研究[J].地震工程與工程振動,1998,18(4):82.LIU Boquan,BAI Shaoliang,XU Yunzhong,et al.Experimental study of low-cyclic behavior of concrete columns [J].Earthquake Engineering and Engineering Vibration,1998,18(4):82.
[13]劉鳴,劉伯權.鋼筋混凝土柱低周疲勞力學性能分析[J].長安大學學報:自然科學版,2012,32(5):65.LIU Ming,LIU Boquan.Low-cycle fatigue damage for shear behavior of reinforced concrete columns.[J].Journal of Chang’an University:Natural Science Edition,2012,32(5):65.
[14]Spacone E,Filippou F C,Taucer F F.Fibre beam-column model for non-linear analysis of R/C frames.Part I:formulation[J].Earthquake Engineering and Structural Dynamics,1996,25(7):711.
[15]Spacone E,Filippou F C,Taucer F F.Fibre beam-column model for non-linear analysis of R/C frames.Part II:applications[J].Earthquake Engineering and Structural Dynamics,1996,25(7):727.
[16]Coffin L F,Jr.A study of the effects of cyclic thermal stresses on a ductile metal[J].American Society of Mechanical Engineers,1954,76:931.
[17]Manson S S.Behavior of materials under conditions of thermal stress[R].Washington D C:National Advisory Committee for Aeronautics,1953.
[18]Smith K,Topper T,Watson P.A stress-strain function for the fatigue of metals:stress-strain function for metal fatigue including mean stress effect[J].Journal of Materials,1970,5:767.
[19]Lorenzo F,Laird C.A new approach to predicting fatigue life behavior under the action of mean stresses[J].Materials Science and Engineering,1984,62(2):205.
[20]Koh S,Stephens R.Mean stress effects on low cycle fatigue for a high strength steel[J].Fatigue &Fracture of Engineering Materials &Structures,1991,14(4):413.
[21]Tong X Y,Wang D J,Xu H.Investigation of cyclic hysteresis energy in fatigue failure process[J].International Journal of Fatigue,1989,11(5):353.
[22]Sugiura K,Chang K,Lee G.Evaluation of low-cycle fatigue strength of structural metals[J].Journal of Engineering Mechanics,1991,117(10):2373.
[23]Mander J,Panthaki F,Kasalanati A.Low-cycle fatigue behavior of reinforcing steel[J].Journal of Materials in Civil Engineering,1994,6(4):453.
[24]Brown J,Kunnath S K.Low-cycle fatigue behavior of longitudinal reinforcement in reinforced concrete bridge columns[R].Buffalo:Multidisciplinary Center for Earthquake Engineering Research,2000.
[25]Scott B,Park R,Priestley M.Stress-strain behavior of concrete confined by overlapping hoops at low and high strain rates[J].ACI Journal Proceedings,1982,79(2):13.
[26]Chang G A,Mander J B.Seismic energy based fatigue damage analysis of bridge columns part 1 Evaluation of Seismic Capacity [R].Buffalo:National Center for Earthquake Engineering Research,1994.
[27]Gomes A,Appleton J.Nonlinear cyclic stress-strain relationship of reinforcing bars including buckling[J].Engineering Structures,1997,19(10):822.
[28]Kunnath S K,Heo Y,Mohle J F.Nonlinear uniaxial material model for reinforcing steel bars[J].Journal of Structural Engineering,ASCE,2009,135(4):335.
[29]Mazzoni S,Mckenna F,Scott M H,et al.OpenSees command language manual[R].Berkeley:Pacific Earthquake Engineering Research Center,2007.
[30]范立礎,卓衛(wèi)東.橋梁延性抗震設計[M].北京:人民交通出版社,2001.FAN Lichu,ZHUO Weidong.Ductility seismic design of bridge [M].Beijing:China Communications Press,2001.
[31]Takemura H,Kawashima K.Effect of loading hysteresis on ductility capacity of reinforced concrete bridge piers[J].Journal of Structural Engineering,1997,43A:849.
Seismic Fragility Analysis of Bridge Piers Based on Low-cycle Fatigue Life of Reinfrocing Steel
SU Junsheng1,WANG Junjie1,GUO Jin2,PENG Yundong3
(1.College of Civil Engineering,Tongji University,Shanghai 200092,China;2.College of Civil Engineering,Shijiazhuang Tiedao University,Shijiazhuang 050043,China;3.China Highway Planning and Design Institute Co.,Ltd.,Beijing 100088,China)
A nonlinear fiber beam-column element was used to simulate quasi-static tests of circular concrete bridge piers reinforced with HRB335 and HRB500E longitudinal bars based on Coffin-Manson low-cycle fatigue damage model for reinforcing steel.The effects of low-cycle fatigue damage on the strength degradation of reinforcing steel and reinforced concrete specimens were studied.In order to further verify the effect of low-cycle fatigue of reinforcing steel on the cumulative damage of reinforced concrete specimens,two quasi-static tests proposed by Takemura and Kunnath with different loading modes were conducted.The results indicate that,the fiber beam-column element could simulate the cumulative damage and strength degradation of reinforced concrete columns under different loading modes when considering the low-cycle fatigue of reinforcing steel on material level by using Coffin-Manson model.
reinforced concrete;bridge piers;low-cycle fatigue;cumulative damage;strength degradation;fiber beam-column element
TU375.3
A
0253-374X(2016)01-0029-08
10.11908/j.issn.0253-374x.2016.01.005
2014-12-26
貴州省“十二五”重大科技專項([2011]6014);國家"九七三"重點基礎研究發(fā)展計劃(2013CB036305);國家自然科學基金(51278373,51438010,51508351)
蘇俊?。?990—),男,博士生,主要研究方向為橋梁抗震和高性能材料.E-mail:sujunsheng@yeah.net
王君杰(1962—),男,教授,博士生導師,工學博士,主要研究方向為橋梁抗震和橋梁船撞.E-mail:jjwang@#edu.cn