賈立哲，陳德珅

汶川地震鋼混框排架主廠房反應(yīng)特性分析

賈立哲，陳德珅

（哈爾濱工業(yè)大學(xué)土木工程學(xué)院，150090哈爾濱）

在對比El Centro波、Taft波、汶川-江油波等不同強震主廠房反應(yīng)特性的基礎(chǔ)上，直接采用不同峰值的汶川-江油波開展動力彈塑性時程分析，探究主廠房這類典型結(jié)構(gòu)在汶川這類特殊地震動下的響應(yīng)特性和震損機理.結(jié)果表明：汶川-江油波作用下結(jié)構(gòu)的響應(yīng)偏小，但變化更不規(guī)律；主廠房能夠抵御遠超設(shè)防烈度的汶川-江油地震作用而不發(fā)生嚴(yán)重的失效破壞；結(jié)構(gòu)扭轉(zhuǎn)效應(yīng)顯著，煤斗層始終為結(jié)構(gòu)薄弱層；框架部分剛度衰退快于排架部分，排架部分承擔(dān)總地震作用比例先提高后不變，總體比例處于較低水平，無法成為結(jié)構(gòu)的第二道抗震防線；雖汽機房網(wǎng)架屋面不會因桿件強度破壞而垮塌，但由于支座處承受剪力巨大且變形過大，支座的率先破壞將引發(fā)屋蓋系統(tǒng)的失穩(wěn)破壞；主廠房率先進入塑性的構(gòu)件主要集中在底層和煤斗層，表現(xiàn)為明顯的“強梁弱柱”特點.

汶川地震；框排架；主廠房；屋蓋系統(tǒng)；地震反應(yīng)

汶川地震震級大、震源淺、烈度高，釋放的能量特別大，破壞性極強.與一般點源地震動圓形衰減規(guī)律不同，其衰減特征為瘦長形衰減［1-2］，其加速度反應(yīng)譜中長周期成分較多，持時特別長［3］，波及范圍特別廣，造成災(zāi)難性的破壞.大型火力發(fā)電廠是重要的生命線工程，其核心主廠房為了適應(yīng)生產(chǎn)工藝的要求，通常采用鋼筋混凝土框排架結(jié)構(gòu)體系.由于結(jié)構(gòu)的不規(guī)則性，在汶川地震中遭受了不同程度破壞，直接和間接損失慘重［4-5］.因此，研究鋼混框排架主廠房這類典型結(jié)構(gòu)在汶川地震這類特殊地震動作用下的響應(yīng)具有重要現(xiàn)實意義.現(xiàn)有研究中，吳濤等［6-7］對于鋼混框排架結(jié)構(gòu)的設(shè)計方法、單雙向地震作用效應(yīng)、子結(jié)構(gòu)擬動力試驗、結(jié)構(gòu)體系創(chuàng)新等方面進行了系統(tǒng)深入的研究，但未涉及該類結(jié)構(gòu)在汶川強震作用下的反應(yīng)特征.

本文基于ANSYS軟件首先通過對比不同強震作用下主廠房的變形和內(nèi)力響應(yīng)，揭示鋼混框排架主廠房這類典型結(jié)構(gòu)在汶川地震動下的響應(yīng)特性.然后直接采用江油地震臺記錄的地震動進行彈塑性時程分析，以期模擬主廠房在汶川地震中的真實受力狀態(tài)，通過對其位移反應(yīng)、內(nèi)力響應(yīng)、框排架協(xié)同工作性、汽機房屋面網(wǎng)架的破壞情況及結(jié)構(gòu)塑性發(fā)展情況的分析，以探究該類結(jié)構(gòu)的震損原因.

圖1 原型主廠房結(jié)構(gòu)布置（mm）

1 框排架主廠房結(jié)構(gòu)模型

以單機容量600 MW的某已投產(chǎn)火電廠鋼混框排架主廠房為原型［7］，結(jié)構(gòu)布置見圖1，橫向為框排架結(jié)構(gòu)，縱向為框架-剪力墻和框架-支撐體系，工程按7度設(shè)防，8度抗震構(gòu)造要求，汽機房屋面為變厚度雙面起坡網(wǎng)架，面板為現(xiàn)澆鋼筋混凝土板.選用BEAM188單元模擬框架梁、柱，按配筋面積相等的原則將鋼筋均勻分配在截面邊緣來模擬鋼筋混凝土組合截面；選用LINK8單元模擬汽機房網(wǎng)架，對于僅考慮受拉的支撐選用LINK10單元模擬；選用SHELL181單元模擬樓板、屋面板及剪力墻；選用MASS21單元模擬非結(jié)構(gòu)構(gòu)件的質(zhì)量，將其分布到相應(yīng)樓層上.混凝土和鋼材的應(yīng)力-應(yīng)變關(guān)系分別采用GB50010—2010《混凝土結(jié)構(gòu)設(shè)計規(guī)范》［8］中分段式模型和理想彈塑性模型.

2 不同地震作用下主廠房地震響應(yīng)

考慮到框排架結(jié)構(gòu)抗震性能分析的普遍適用性，地震動選擇考慮Ⅱ、Ⅲ類場地條件，取El Centro波、Taft波、汶川-江油波三條實際地震動進行計算，地震動持時為20 s，按抗震規(guī)范調(diào)整有效峰值加速度至設(shè)防烈度7度（0.15 g）罕遇地震所對應(yīng)的310 Gal.每條地震動水平加速度記錄中峰值較大的沿結(jié)構(gòu)橫向、峰值較小的沿結(jié)構(gòu)縱向，豎向加速度記錄沿結(jié)構(gòu)高度方向同時輸入［9-10］.

2.1 位移反應(yīng)

結(jié)構(gòu)頂層縱橫向位移時程曲線見圖2，位移反應(yīng)對比統(tǒng)計見圖3.可看出，結(jié)構(gòu)位移反應(yīng)相差較大，頂層縱橫向位移在不同時刻達到最大值，波形曲線形狀變化較大；El Centro波作用下結(jié)構(gòu)橫向位移反應(yīng)最為突出，Taft波作用下結(jié)構(gòu)縱向位移反應(yīng)最為突出，汶川-江油波作用下結(jié)構(gòu)縱橫向位移反應(yīng)相對較小，但排架柱橫向變形更不均勻.

圖2 不同地震動作用頂層位移時程曲線

圖3 不同強震作用結(jié)構(gòu)位移反應(yīng)對比統(tǒng)計

2.2 內(nèi)力響應(yīng)

框架部分層間剪力及排架柱剪力對比統(tǒng)計見圖4，總體上汶川-江油波作用下結(jié)構(gòu)內(nèi)力響應(yīng)較其他波小，但框架底層橫向?qū)娱g剪力大于其他強震作用的情況，結(jié)構(gòu)均在第四層煤斗層發(fā)生突變，煤斗層始終為結(jié)構(gòu)薄弱層；排架柱剪力始終為Taft波作用時最大，El Centro波作用時次之，汶川-江油波作用時最小，頂層與底層排架柱剪力差別相比中間位置大，由于汽機房采用面內(nèi)剛度較大網(wǎng)架上覆鋼筋混凝土板的形式，水平地震作用通過支座直接傳遞到柱頂，故排架柱頂層剪力較中間部位大.

圖4 不同地震動作用結(jié)構(gòu)各層剪力對比統(tǒng)計

3 汶川地震主廠房地震反應(yīng)分析

3.1 位移反應(yīng)

以汶川-江油實際地震動為激勵，對比分析峰值加速度為小震55 Gal、大震310 Gal和實際地震511 Gal時，結(jié)構(gòu)框架部分縱橫向的層位移、層間位移角及排架柱柱頂縱橫向位移見圖5.從圖中可看出:結(jié)構(gòu)縱橫向位移反應(yīng)相差明顯，橫向較縱向位移反應(yīng)大得多.由于結(jié)構(gòu)橫向的不規(guī)則性問題較縱向突出，地震作用引起結(jié)構(gòu)薄弱部位變形集中，率先進入塑性受力狀態(tài)，甚至提前發(fā)生破壞.中國抗震規(guī)范［11］要求鋼筋混凝土框架-抗震墻結(jié)構(gòu)彈性及彈塑性最大層間位移角應(yīng)分別小于1／800和1／100，橫向?qū)娱g位移角曲線在第四層煤斗層位置位移角有明顯突變，彈性階段結(jié)構(gòu)最大層間位移角為1／1 299，彈塑性階段結(jié)構(gòu)最大層間位移角為1／214，均能滿足規(guī)范要求.隨著地震作用的增大，排架柱橫向變形不均勻程度增加，這種趨勢加劇了屋蓋系統(tǒng)扭轉(zhuǎn)效應(yīng)，對屋蓋系統(tǒng)支座受力極為不利.

3.2 內(nèi)力響應(yīng)

框排架主廠房模型框架部分層間剪力分布、排架柱各層剪力分布見圖6.層間剪力最大值都出現(xiàn)在煤倉間第四層煤斗層，由于煤斗大梁為大尺寸深梁，局部剛度、質(zhì)量過分懸殊，在層間剛度分配中容易吸收較大地震作用，在設(shè)計驗算中應(yīng)予以足夠重視.由于排架柱縱向靠連梁和柱間支撐連成一體，剛度較排架橫向平面外剛度大得多，排架柱縱向剪力明顯大于橫向剪力；屋蓋系統(tǒng)水平地震作用由支座直接傳遞到排架柱柱頂，故排架柱上部剪力較中間部位大，如果支座強度或構(gòu)造不滿足要求將率先破壞，危及屋蓋整體安全.

圖5 不同強度地震動作用結(jié)構(gòu)位移反應(yīng)對比統(tǒng)計

圖6 不同強度地震動作用結(jié)構(gòu)各層剪力對比統(tǒng)計

3.3 框排架協(xié)同工作性能

鋼混框排架結(jié)構(gòu)在遭受地震作用時，框架部分與排架部分按剛度比和變形協(xié)調(diào)條件分配地震作用.在端部設(shè)有剪力墻的鋼混框排架結(jié)構(gòu)中，剪力墻能夠吸收縱向結(jié)構(gòu)體系90%以上的地震剪力，分析中只考慮框排架柱的基底剪力.不同峰值加速度的汶川-江油波作用下，排架承擔(dān)的地震剪力與基底總剪力之間統(tǒng)計關(guān)系見圖7（a）、（b）.在結(jié)構(gòu)彈性階段，各模型中排架承擔(dān)的縱、橫向地震剪力分別占結(jié)構(gòu)基底縱、橫向總剪力的18.99%和4.02%；隨著地震動峰值加速度的提高，框架部分剛度退化快于排架部分，在罕遇地震作用下，排架部分承擔(dān)的縱橫向地震作用占總地震作用的比例較多遇地震時均有提高，分別提高到24.07%和4.63%，但進入塑性階段后比例基本不變，分別為24.19%和4.69%.各排架柱承擔(dān)縱橫向剪力與所在軸線基底總剪力比見圖7（c）、（d），由圖中結(jié)果可知，排架柱縱向剪力與相應(yīng)軸線基底剪力的比例明顯大于橫向，①、②和⑦、⑧軸線排架柱縱向剪力占所在軸線總剪力比例低于其他軸線，主要由于該處設(shè)有剪力墻，與剪力墻相連的框架柱剪力較大，①、②和⑦、⑧軸線總剪力明顯大于其他軸線.結(jié)構(gòu)在塑性階段，框架與排架塑性內(nèi)力重分布不明顯，總體上排架部分仍處于較低水平，幾乎無法起到繼框架-剪力墻之后第二道抗震防線作用.

3.4 汽機房屋面網(wǎng)架破壞分析

通過構(gòu)件強度驗算和整體變形驗算來評估屋蓋系統(tǒng)的地震安全性.由于罕遇地震屬于偶然荷載，構(gòu)件截面強度驗算及結(jié)構(gòu)變形驗算均采用重力荷載和罕遇地震作用組合的標(biāo)準(zhǔn)值.構(gòu)件強度驗算以應(yīng)力比為衡量指標(biāo)，構(gòu)件驗算應(yīng)力值與材料設(shè)計應(yīng)力之比為構(gòu)件應(yīng)力比，其表達式為

式中:N為構(gòu)件軸力，φ為軸心受壓構(gòu)件穩(wěn)定系數(shù)，其值按GB50017—2003《鋼結(jié)構(gòu)設(shè)計規(guī)范》要求計算（受拉構(gòu)件φ=1.0），A為構(gòu)件截面面積，f為鋼材屈服強度設(shè)計值，屋架結(jié)構(gòu)鋼材為Q235級，屈服強度設(shè)計值取215 MPa.

圖7 不同強度地震動作用框排架承擔(dān)縱橫向剪力情況統(tǒng)計

網(wǎng)架屋面超限桿件（應(yīng)力比大于1.0）應(yīng)力比云圖及其統(tǒng)計圖見圖8（a）、（b）.應(yīng)力比云圖為正放四角錐網(wǎng)架水平投影圖，水平和豎直桿件為上、下弦桿，傾斜桿件均為腹桿.從圖8（a）可看出，網(wǎng)架超限桿件總數(shù)為238根，占桿件總數(shù)比例達到8.85%，超限桿件大部分為腹桿，少量弦桿超限破壞，主要受力桿件應(yīng)力水平仍在安全范圍內(nèi)，網(wǎng)架屋面體系因桿件強度破壞而導(dǎo)致整體垮塌的可能性較小.但由于屋蓋系統(tǒng)采用無檁體系，鋼筋混凝土板直接鋪設(shè)在支撐系統(tǒng)上弦，對上弦桿件起到有利保護作用.在極罕遇地震作用下，如果屋面板由于連接失效等原因破壞后，結(jié)構(gòu)上弦桿件應(yīng)力水平將迅速升高，還將有部分上弦桿件發(fā)生強度破壞.

網(wǎng)架屋面總位移極值云圖見圖8（c）.圖中上部深色位移值最小，圖右側(cè)灰色位移值最大.由計算結(jié)果可知:屋蓋整體位移極值發(fā)生于邊榀排架柱支座處，其值為0.164 m，橫向位移極值亦為該處值為0.151 m，排架柱支承側(cè)縱向位移明顯大于框架柱支承側(cè)，其最大值為0.063 m，扭轉(zhuǎn)變形趨勢相當(dāng)明顯.豎向位移極值發(fā)生于邊跨跨中，其值為0.078 m.網(wǎng)架屋面排架柱頂支座附近位移極為不均勻，兩端支座位移大于中間支座，支座平面扭轉(zhuǎn)效應(yīng)顯著.網(wǎng)架支座為平板壓力支座，連接螺栓為4M24（8.8級普通螺栓），支座抗剪承載力為451.8 kN.本文計算得出網(wǎng)架支座承擔(dān)最大地震剪力高達975 kN，超出其抗剪承載力一倍以上，因此強度或構(gòu)造不滿足要求的支座必然率先發(fā)生破壞，改變了網(wǎng)架原有的受力狀態(tài)，最終導(dǎo)致了屋蓋系統(tǒng)失去穩(wěn)定直至垮塌.這在一定程度上證明了汶川地震中江油電廠汽機房網(wǎng)架屋面的破壞原因.

圖8 網(wǎng)架屋面強度校核和整體變形

3.5 結(jié)構(gòu)塑性發(fā)展情況

以時間間隔0.02 s的汶川-江油實際地震動作為激勵，定性研究主廠房結(jié)構(gòu)構(gòu)件應(yīng)力變化發(fā)展順序，借此判斷結(jié)構(gòu)中的薄弱構(gòu)件.由于構(gòu)件處于復(fù)雜應(yīng)力狀態(tài)，取Von Mises應(yīng)力作為衡量指標(biāo).

不同時刻結(jié)構(gòu)應(yīng)力分布情況見圖9.在瞬態(tài)分析初期所有構(gòu)件均未出現(xiàn)較大的應(yīng)力區(qū)域；隨著加載的持續(xù)進行，1.76 s時刻底層柱底、煤倉間大梁梁柱節(jié)點、除氧間大梁梁柱節(jié)點逐步出現(xiàn)應(yīng)力集中的現(xiàn)象；在3.66 s時刻，結(jié)構(gòu)底層及關(guān)鍵梁柱節(jié)點區(qū)應(yīng)力區(qū)域逐步擴大，剪力墻、煤斗層柱端及排架柱柱底應(yīng)力逐步提高；在8.08 s時刻，與剪力墻相連的框架柱底層柱端及煤斗層梁柱節(jié)點區(qū)應(yīng)力水平提高較快，結(jié)構(gòu)兩端邊榀構(gòu)件應(yīng)力明顯大于相應(yīng)中間榀構(gòu)件，柱尤其是梁柱節(jié)點部位應(yīng)力水平明顯高于梁；隨后加載至12.44 s時，煤斗層柱端、除氧器大梁梁柱節(jié)點、排架柱柱底應(yīng)力水平顯著提高，剪力墻與框架交接處應(yīng)力水平高于墻身內(nèi)部；隨著時間推移到19.62 s時刻，原有高應(yīng)力區(qū)域逐步擴展.

圖9 不同時刻結(jié)構(gòu)應(yīng)力分布云圖（未顯示樓板與屋面板）

端部設(shè)有剪力墻的鋼混框排架結(jié)構(gòu)體系在汶川-江油實際地震動作用下，主體結(jié)構(gòu)安全性是能夠得到保證的.由于結(jié)構(gòu)錯層、局部縮進等特殊構(gòu)造，各層間豎向承載力變化較大，各層破壞程度差異明顯，底層和煤斗層構(gòu)件應(yīng)力明顯高于其他部位.由此上述分析可定性判斷此類結(jié)構(gòu)總體上的屈服順序:與剪力墻相連的框架柱底端-煤倉間大梁、除氧間大梁梁柱節(jié)點區(qū)-煤倉間大梁所在的第四層柱端-框架主梁、底層其他柱柱底-剪力墻與框架交接處，表現(xiàn)為明顯的“強梁弱柱”特點.為滿足煤斗層的特殊工藝要求，煤斗大梁寬700 mm，高3 600 mm，跨度12 m，跨高比小于5，屬深受彎構(gòu)件，在層間剛度分配中，由于剛度過大在地震過程中承擔(dān)了絕大部分的層間地震作用，并通過梁柱節(jié)點傳遞給相鄰的框架柱，“塑性鉸”難以在梁上形成，梁柱節(jié)點及柱端位置成為薄弱環(huán)節(jié).

中國建筑抗震設(shè)計通常強調(diào)“強柱弱梁、強節(jié)點弱構(gòu)件”［11］，保證結(jié)構(gòu)的屈服順序為“先梁后柱再節(jié)點”，使結(jié)構(gòu)在地震作用下能夠反復(fù)吸收、耗散地震能量，不至于在短時間內(nèi)形成破壞機構(gòu).而鋼混框排架結(jié)構(gòu)的屈服順序顯然不能滿足要求，結(jié)構(gòu)薄弱層為結(jié)構(gòu)底層和煤斗層，但由于整體結(jié)構(gòu)剛度大、安全儲備極高，在極罕遇地震中主體結(jié)構(gòu)仍能保證安全.汶川地震中江油電廠震害也證明了這一點:主體結(jié)構(gòu)基本安全，除汽機房屋面垮塌和填充墻大量破壞外，結(jié)構(gòu)底層柱混凝土剝落鋼筋外露，部分支承梁的牛腿局壓破壞，部分梁柱節(jié)點開裂嚴(yán)重.

4 結(jié) 論

1）與El Centro波、Taft波相比，汶川-江油波作用下結(jié)構(gòu)的位移響應(yīng)和內(nèi)力響應(yīng)偏小，但其變化趨勢更不規(guī)律.

2）框排架主廠房能夠抵御遠超設(shè)防烈度的實際汶川地震作用而不發(fā)生嚴(yán)重的失效破壞.結(jié)構(gòu)扭轉(zhuǎn)效應(yīng)顯著，煤斗層始終為結(jié)構(gòu)薄弱層.隨地震作用增大排架柱橫向變形不均勻程度增加，加劇了屋蓋系統(tǒng)扭轉(zhuǎn)效應(yīng)，對屋蓋支座受力極為不利.

3）隨地震動峰值提高，框架部分剛度衰退快于排架部分，排架部分承擔(dān)總地震作用比例提高.但結(jié)構(gòu)進入塑性階段后，排架部分承擔(dān)總體地震作用的比例幾乎不變，總體比例處于較低水平，無法成為繼框架-剪力墻之后的第二道抗震防線.

4）汽機房屋面網(wǎng)架強度超限桿件比例為8.85%，大部分為次要受力桿件腹桿，不會因大量桿件強度破壞而導(dǎo)致屋蓋整體垮塌.但網(wǎng)架支座處承受剪力最大值超出其抗剪承載力一倍以上，加之屋蓋扭轉(zhuǎn)導(dǎo)致支座變形過大，強度或構(gòu)造不滿足要求的支座率先破壞，將引發(fā)屋蓋系統(tǒng)的失穩(wěn)破壞.從數(shù)值模擬的角度給出了汶川地震中江油電廠汽機房網(wǎng)架屋面的垮塌原因.

5）主廠房結(jié)構(gòu)各層破壞程度差異明顯，底層和煤斗層構(gòu)件應(yīng)力明顯高于其他部位.定性分析此類結(jié)構(gòu)的屈服順序為:與剪力墻相連的框架柱底端-煤倉間大梁、除氧間大梁梁柱節(jié)點區(qū)-煤倉間大梁所在的第四層柱端-框架主梁、底層其他柱柱底-剪力墻與框架交接處，表現(xiàn)為明顯的“強梁弱柱”特點.

［1］劉浪，李小軍.汶川8.0級地震地震動峰值加速度衰減特性分析［J］.北京工業(yè)大學(xué)大學(xué)學(xué)報，2012，38（2）:173-179.

［2］LI Xiaojun.Strong motions and engineering structure performancesinrecentmajorearthquakes［J］. Earthquake Science，2010，23（1）:1-3.

［3］WEN Zengping，XIE Junju.Near-source strong ground motion characteristics of the 2008 Wenchuan earthquake［J］.Bulletin of the Seismological Society of America，2008，100（5B）:2425-2439.

［4］KRAUSMANN E，CRUZ A M，AFFELTRANGER B. The impact of the 12 may 2008 Wenchuan earthquake on industrial facilities［J］.Journal of Loss Prevention in the Process Industries，2010，23（2）:242-248.

［5］姜云甫，黃佑驗.江油電廠主廠房震害分析［J］.武漢大學(xué)學(xué)報:工學(xué)版，2009，42:172-176.

［6］吳濤.大型火力發(fā)電廠鋼筋混凝土框排架結(jié)構(gòu)抗震性能及設(shè)計方法研究［D］.西安:西安建筑科技大學(xué)，2003.

［7］WU Tao，BAI Guoliang，LIU Boquan.Experimental study on seismic behavior of RC frame-bent structure of a large-scale power plant［J］.Journal of Building Structures，2007，28（3）:46-52.

［8］GB50010—2010混凝土結(jié)構(gòu)設(shè)計規(guī)范［S］.北京:中國建筑工業(yè)出版社，2010.

［9］CHOPRA A K.Dynamics of structures:theory and applications to earthquake engineering［M］.3rd ed.New York:Pearson Education Inc，2007.

［10］NAKAMURA N，AKITA S，SUZUKI T，et al.Study of ultimate seismic response and fragility evaluation of nuclear power building using nonlinear three-dimensional finite element model［J］.Nuclear Engineering and Design，2010，240（1）:166-180.

［11］GB50011—2010建筑抗震設(shè)計規(guī)范［S］.北京:中國建筑工業(yè)出版社，2010.

（編輯 趙麗瑩）

Seismic response characteristics for R.C.frame-bent main building under Wenchuan earthquake

JIA Lizhe，CHEN Deshen
（School of Civil Engineering，Harbin Institute of Technology，150090 Harbin，China）

The responses of main building under different strong earthquakes，such as El Centro，Taft，Wenchuan-Jiangyou were compared，and dynamic time history analysis based on Wenchuan-Jiangyou record with different peak values was used to reveal the damage mechanism of this typical structure.The results indicate that the response of Wenchuan-Jiangyou record is the smallest，but variability is the biggest. Wenchuan earthquake，which is more than design intensity，does not cause serious damage of main building. The torsional effect is remarkable，and coal scuttle layer is always weak layer.The stiffness degradation of frame is faster than the bent part，and the proportion of bent part within total earthquake action is enhanced firstly and then remains same.But the overall proportion is still in low level and bent part is unable to become the second resistance component.Although space grid roof system in turbine hall does not collapse because of strength failure of components，the bearing carries over huge shearing force with oversized distortion，and the failure of bearing will lead to collapse of roof system.Plastic deformations mainly concentrate in the first floor and coal scuttle layer，and“strong beam and weak column”are exposed.

wenchuan earthquake；frame-bent structure；main building；roof system；seismic response

TU375

A

0367-6234（2014）06-0008-06

2013-07-13.

國家自然科學(xué)基金資助項目（51008100）；

中國博士后科學(xué)基金特別資助項目（201104413）.

賈立哲（1978—），男，副教授，碩士生導(dǎo)師.

賈立哲，jializhe＠hit.edu.cn.