瀝青混合料蠕變的變參數(shù)模型

康永剛1，張秀娥1，張華明2

（1.防災(zāi)科技學(xué)院，河北 三河 065201；2.山西省雷電防護(hù)監(jiān)測(cè)中心，山西 太原 030002）

摘要：為了準(zhǔn)確描述瀝青混合料整個(gè)蠕變過(guò)程，解決現(xiàn)有黏彈性模型的不足，對(duì)經(jīng)典模型進(jìn)行了擴(kuò)展。運(yùn)用非定常參數(shù)替換定常參數(shù)的方法，把經(jīng)典Zener模型、Zener模型與黏壺串聯(lián)結(jié)構(gòu)中的黏性參數(shù)看為時(shí)間的函數(shù)。利用元件串聯(lián)、并聯(lián)時(shí)的力學(xué)關(guān)系，推導(dǎo)改進(jìn)模型本構(gòu)方程和蠕變函數(shù)的一般形式。通過(guò)選擇參數(shù)的具體時(shí)變形式求出具體的蠕變函數(shù)，發(fā)現(xiàn)一些經(jīng)驗(yàn)蠕變公式是改進(jìn)后模型的蠕變函數(shù)；通過(guò)蠕變特征的分析，發(fā)現(xiàn)改進(jìn)后的模型具有較廣的描述范圍；通過(guò)擬合不同瀝青混合料的蠕變?cè)囼?yàn)數(shù)據(jù)，發(fā)現(xiàn)修正后的模型對(duì)試驗(yàn)數(shù)據(jù)有較好的擬合效果。給出的變參數(shù)模型實(shí)現(xiàn)了用較少的參數(shù)對(duì)瀝青混合料3階段蠕變過(guò)程給出比較準(zhǔn)確的描述。

關(guān)鍵詞：道路工程；瀝青混合料；蠕變模型；變參數(shù)；黏壺

中圖分類(lèi)號(hào)：U416.2

文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：1002-0268（2015）07-0013-05

Variable-parameter Creep Modelof Asphalt Mixture

KANG Yong-gang1，ZHANG Xiv-e1，ZHANG Hva-ming2
（1.Institvte of Disaster Prevention，Sanhe Hebei065201，China；2.Thvnder Prevention and Monitoring Center of Shanxi Province，Taiyvan Shanxi030002，China）

Abstract：In order to accvrately describe the creep process of asphaltmixtvre and solving the shortcoming of traditional viscoelastic model，the classic model is improved.Using the method of replacing the stationary parameterswith non-stationary parameters，the viscovs parameters in classical Zener model and the series strvctvre of Zenermodel with dashpot are regarded as time fvnctions.Using the mechanical relationship of parallel connection and series connection，the generalized constitvtive relation and creep fvnction of the improved model are derived.Throvgh choosing the parameters，specific time-varying form，the creep fvnction is solved.It is fovnd that the improved model inclvdes some empirical creep formvlas.Throvgh analyzing creep characteristics，it is fovnd that the improved model has wider description range.Throvgh fitting different creep test data，it is fovnd that the improved model，s theoretic cvrves are coincident well with the testing cvrves.The proposed variable parametermodel realized the accvrate description for the 3-phase creep of asphaltmixtvre with relatively small nvmber of parameters.

Key words：road engineering；asphaltmixtvre；creep model；variable parameter；dashpot

0　引言

瀝青混合料是一種黏彈性材料，在應(yīng)力比較大的荷載作用下，或由于高溫而超過(guò)其極限承載能力時(shí)，其蠕變過(guò)程可分為遷移期、穩(wěn)定期和破壞期3個(gè)階段，分別對(duì)應(yīng)于圖1（a）中AB，BC，CD 段[1-2]，其中ε（t）表示應(yīng)變，t表示時(shí)間。遷移期，變形迅速增大，但應(yīng)變速率隨時(shí)間增加逐漸減??；穩(wěn)定期，應(yīng)變穩(wěn)定增長(zhǎng)，但應(yīng)變速率基本保持不變；破壞期，應(yīng)變、應(yīng)變速率隨時(shí)間增加迅速增大直至

破壞。另一類(lèi)是在應(yīng)力比較小的荷載作用下的蠕變，黏性流動(dòng)變形并不隨荷載作用時(shí)間的延長(zhǎng)而無(wú)限增加，而是隨著時(shí)間的推移，黏性流動(dòng)變形的增量逐漸減小，最終使黏性流動(dòng)變形趨于一個(gè)穩(wěn)定值，表現(xiàn)出一種固結(jié)效應(yīng)，蠕變曲線只有遷移期和穩(wěn)定期兩個(gè)階段，如圖1（b）所示[2]。

圖1　瀝青混合料的蠕變曲線
Fig.1 Creep curves of asphaltm ixture

目前對(duì)瀝青混合料黏彈性本構(gòu)關(guān)系的研究較多，提出了各種黏彈性本構(gòu)模型，如Kelvin模型、Maxwell模型、標(biāo)準(zhǔn)線性固體模型、Bvrgers模型以及廣義Kelvin模型、廣義Maxwell模型等[3-4]。Jvdycki[5]對(duì)常規(guī)和改性瀝青混凝土在蠕變?cè)囼?yàn)條件下的非線性黏彈性行為進(jìn)行了研究，建立了一個(gè)非線性黏彈性模型；鄭健龍[6]建立了瀝青混合料熱黏彈性本構(gòu)模型；徐世法[7]采用黏彈性理論，提出了一個(gè)四單元五參數(shù)模型來(lái)分析瀝青路面的永久變形；王后裕[8]通過(guò)分析瀝青混合料加載時(shí)的固結(jié)效應(yīng)和卸載時(shí)的永久變形等流變學(xué)特征，提出了一個(gè)流變學(xué)模型。但是，以上這些模型只能描述前兩個(gè)蠕變階段。

近些年，建立能描述瀝青混合料3個(gè)蠕變階段的模型成為研究重點(diǎn)[9-13]。如從損傷力學(xué)角度出發(fā)：張久鵬[9]以流變學(xué)模型理論為基礎(chǔ)，綜合考慮損傷和硬化效應(yīng)，建立了一個(gè)統(tǒng)一的瀝青混合料損傷蠕變模型，但該模型參數(shù)過(guò)多。曾國(guó)偉[11]在三元件模型基礎(chǔ)上耦合一個(gè)連續(xù)性損傷因子，構(gòu)建一個(gè)能夠描述瀝青混合料蠕變?nèi)^(guò)程的損傷蠕變模型。馮林[12]提出了一種黏-彈-塑本構(gòu)模型，并采用Kachanov蠕變損傷律進(jìn)行全局耦合。從力學(xué)元件的非線性化角度出發(fā)：張久鵬[10]對(duì)Bvrgers模型中的彈簧進(jìn)行改進(jìn)，認(rèn)為其彈性模量隨應(yīng)變?cè)龃蠖鴾p小，但沒(méi)有考慮黏性元件的非線性。葉永[13]采用多種巖石蠕變模型來(lái)描述瀝青混合料的蠕變，這些模型僅對(duì)描述第3個(gè)階段的元件進(jìn)行非線性化，并且包含元件較多，相應(yīng)參數(shù)也較多。

Zener模型是重要的黏彈性模型，被稱為標(biāo)準(zhǔn)線性固體模型，對(duì)黏彈性材料的應(yīng)力松弛、蠕變都可以定性描述。本文從力學(xué)元件的非線性化角度出發(fā)，用變參數(shù)黏壺替換Zener模型中的黏壺，構(gòu)造出變參數(shù)Zener模型，推導(dǎo)出它的本構(gòu)方程和蠕變函數(shù)。改進(jìn)后的模型串聯(lián)一個(gè)變參數(shù)黏壺適合描述瀝青混合料的整個(gè)蠕變過(guò)程，并且相對(duì)同類(lèi)模型，其包含相對(duì)較少的元件，相應(yīng)的參數(shù)也較少。

1　變參數(shù)Zener模型

經(jīng)典彈簧的應(yīng)力σ（t）和應(yīng)變?chǔ)牛╰）滿足虎克定律σ（t）＝Eε（t），其中t表示時(shí)間；E表示彈性模量。經(jīng)典的黏壺滿足牛頓定律，其中η表示黏滯系數(shù)。Zener模型是一個(gè)Maxwell模型（一個(gè)彈簧和一個(gè)黏壺的串聯(lián)）并聯(lián)一個(gè)彈簧的結(jié)構(gòu)，對(duì)黏彈性材料的應(yīng)力松弛、蠕變都可以定性描述，但要準(zhǔn)確地描述黏彈過(guò)程，需要比較復(fù)雜的模型。

建立變參數(shù)模型的出發(fā)點(diǎn)是，在應(yīng)力作用下，材料的微觀結(jié)構(gòu)發(fā)生改變，其中某些力學(xué)參數(shù)隨著時(shí)間的變化而變化，需要考慮力學(xué)參數(shù)的時(shí)間效應(yīng)。在此引入一個(gè)變參數(shù)的黏壺：式中黏滯系數(shù)η1（t）變化為時(shí)間的函數(shù)。黏壺是一種由應(yīng)變速率定義的力學(xué)元件，變參數(shù)η1（t）表明材料的黏彈特性與應(yīng)變速率的關(guān)系隨著時(shí)間變化。用此黏壺代替Zener模型中的牛頓黏壺，給出一種變參數(shù)Zener模型，如圖2所示，利用元件串聯(lián)、并聯(lián)時(shí)的應(yīng)力應(yīng)變關(guān)系，推導(dǎo)出該模型的本構(gòu)方程為：

常應(yīng)力σ0作用下對(duì)式（2）分離變量求定積分，注意t＝0時(shí)，σ0＝ε0（E2＋E3），得到變參數(shù)Zener模型蠕變函數(shù)一般表達(dá)式：

圖2　變參數(shù)Zener模型
Fig.2 Variable-parameter Zener model

徐世法[7]提出的四單元五參數(shù)模型，即一個(gè)改進(jìn)的Bvrgers模型，它取串聯(lián)黏壺的黏滯系數(shù)為η（t）＝A eBt來(lái)反映瀝青混合料的固結(jié)效應(yīng)。設(shè)黏滯系數(shù)隨時(shí)間冪律變化：

式中，η10，α為常數(shù)且α＞0，α≠0。當(dāng)0＜α＜1時(shí)，表示黏滯系數(shù)隨時(shí)間增加；當(dāng)α＞1時(shí)，表示黏滯系數(shù)隨時(shí)間減小。但是時(shí)間趨于零時(shí)，α＞1會(huì)導(dǎo)致黏滯系數(shù)趨于無(wú)窮，所以表示黏滯系數(shù)隨時(shí)間的減小，應(yīng)采用其他形式，如η（t）＝η10（t＋C）1-α，其中α＞1，C為常數(shù)且C＞0。在此取0＜α＜1，即黏滯系數(shù)隨時(shí)間冪律增加來(lái)表示固結(jié)效應(yīng)，所以下面仍然以式（4）為例討論。把式（4）代入式（3）得到：

式中，，和k3＝。式（5）中的稱為stretched exponential函數(shù)，張為民用該函數(shù)給出的蠕變?nèi)崃繉?shí)用表達(dá)式，與式（5）形式相同[14]。當(dāng)t→∞時(shí)，式（5）的應(yīng)變趨于穩(wěn)定值k1。

圖3是變參數(shù)Zener模型式（5）的蠕變曲線（k3取不同值），參數(shù)分別為k1＝0.012，k2＝0.01和α＝0.8。可見(jiàn)，變參數(shù)Zener模型的蠕變只有遷移期和穩(wěn)定期兩個(gè)階段，是應(yīng)力比較小的荷載作用下的蠕變模型。圖4（a）為單軸貫入試驗(yàn)獲得的瀝青混凝土的蠕變?nèi)崃?span>[15]。注意：蠕變時(shí)常應(yīng)力σ0、應(yīng)變?chǔ)牛╰）和蠕變?nèi)崃縅（t）滿足ε（t）＝σ0J（t）。圖4（b）為25℃、靜載為100 kPa下瀝青混合料AC-13C，通過(guò)UTM進(jìn)行室內(nèi)靜態(tài)蠕變?cè)囼?yàn)的數(shù)據(jù)[16]。圖4中曲線為式（5）的擬合結(jié)果，擬合參數(shù)列于表1，注意參數(shù)k1和k2為無(wú)量綱。

圖3　變參數(shù)Zener模型的蠕變曲線
Fig.3 Creep curves of variable-parameter Zener model

圖4　蠕變?cè)囼?yàn)數(shù)據(jù)及式（5）擬合曲線
Fig.4 Creep test data and fitting curves of Eq.（5）

表1　變參數(shù)Zener模型參數(shù)
Tab.1 Parameters of variable-param eter Zener model

k1　k2　k3/s-1α Fig.4（a）　2.639σ0　2.32σ00.130　0.098 Fig.4（b）0.067 1　0.067 1　0.042 4　0.414

2　描述瀝青混合料整個(gè)蠕變過(guò)程的模型

變參數(shù)Zener模型無(wú)法描述破壞期的情形，需要對(duì)它進(jìn)行改進(jìn)。在此，通過(guò)串聯(lián)一個(gè)描述第3階段蠕變的變參數(shù)黏壺η2（t），如圖5所示，給出適合描述整個(gè)蠕變過(guò)程的模型。力學(xué)元件串聯(lián)時(shí)應(yīng)力相同，合應(yīng)變?yōu)楦鞑⒙?lián)單元應(yīng)變之和，蠕變函數(shù)一般表達(dá)式為：

張久鵬[17]假設(shè)黏滯系數(shù)是先增大后減小，作如下改進(jìn)：

式中，a，b為常數(shù)，a＞0，b＜0且b2-4a＜0；η20

為初始黏滯系數(shù)。把式（4）和式（7）代入式（6），得到：

圖5　改進(jìn)的變參數(shù)Zener模型
Fig.5 Improved variable-parameter Zener model

為了描述蠕變第3階段破壞期，根據(jù)此階段的特點(diǎn)，在此假設(shè)串聯(lián)黏壺的黏滯系數(shù)隨時(shí)間指數(shù)減少為：

把式（4）和式（9）代入式（6），得到：

式中k4＝σ0/η20。

圖6是改進(jìn)的變參數(shù)Zener模型式（10）的蠕變曲線，參數(shù)分別為k1＝0.012，k2＝0.01，k3＝0.1，k4＝0.000 1，α＝0.8和γ＝0.05?？梢?jiàn)，改進(jìn)的變參數(shù)Zener模型的蠕變有3個(gè)階段，是應(yīng)力比較大的荷載作用下的蠕變模型。瀝青混合料種類(lèi)和蠕變條件的多樣性，黏壺η2（t）的蠕變函數(shù)也可取其他形式，如冪律型σ0η20tm，冪律指數(shù)型σ0η20tmeβt等，對(duì)應(yīng)的黏滯系數(shù)很容易反推出。

圖6　變參數(shù)Zener模型的蠕變曲線
Fig.6 Creep curve of variable-param eter Zener model

圖7為張久鵬[9]等人在UTM試驗(yàn)機(jī)上分別進(jìn)行SMA13和AC13兩種類(lèi)型瀝青混合料的單軸靜載蠕變?cè)囼?yàn)。為加快蠕變?cè)囼?yàn)進(jìn)程，試驗(yàn)荷載為0.7 MPa，試驗(yàn)溫度為60℃。SMA13采用SBS改性瀝青，AC13采用AH70#普通瀝青，混合料的級(jí)配、油石質(zhì)量比詳見(jiàn)文獻(xiàn)[17]。圖7中曲線為式（10）的擬合結(jié)果，以后計(jì)劃開(kāi)展模型參數(shù)的識(shí)別研究，本文對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行分段擬合：由于非定常黏壺η2（t）對(duì)蠕變的前兩個(gè)階段影響較小，首先用最小二乘法和式（5）擬合前兩個(gè)階段的數(shù)據(jù)（AC取t＜210 s，SMA取t＜1 000 s），得到參數(shù)k1，k2，k3和α的值；保持這些參數(shù)值不變，用式（10）擬合所有試驗(yàn)數(shù)據(jù)得到參數(shù)k4和γ，擬合參數(shù)列于表2，注意參數(shù)k1，k2和k4為無(wú)量綱。

圖7　蠕變?cè)囼?yàn)數(shù)據(jù)及式（10）擬合曲線
Fig.7 Creep test data and fitting curves of Eq.（10）

表2　改進(jìn)后的變參數(shù)Zener模型參數(shù)
Tab.2 Parameters of im proved variable-parameter Zener model

k1　k2　k/s-1　k4α　γ3AC　0.022 5 0.016　0.034　0.000 040 0.40　0.013 SMA　0.032 7 0.029　0.000 105 0.000 98　0.20　0.001 05

3　結(jié)論

在應(yīng)力作用下，材料的微觀結(jié)構(gòu)發(fā)生改變，其某些力學(xué)參數(shù)隨著時(shí)間的變化而變化，需要考慮力學(xué)參數(shù)的時(shí)間效應(yīng)。用變參數(shù)黏壺替換Zener模型中的黏壺，構(gòu)造出變參數(shù)Zener模型，推導(dǎo)出它的本構(gòu)方程和蠕變函數(shù)。改進(jìn)后的模型串聯(lián)一個(gè)描述第3階段蠕變的變參數(shù)黏壺，適合描述瀝青混合料整個(gè)蠕變過(guò)程。用所建立的模型分別擬合不同的蠕變?cè)囼?yàn)數(shù)據(jù)，都給出較好的描述。

參考文獻(xiàn)：

References：

[1] WITCZAK M W，KALOUSH K，PELLINEN T.Simple Performance Test for Svperpave Mix Design，NCHRP Report 465 [ R]. Washington，D. C.：National Academy Press，2002.

[2] 詹小麗，張肖寧，王端宜，等.改性瀝青非線性粘彈

性本構(gòu)關(guān)系研究及應(yīng)用[J].工程力學(xué)，2009，26 （4）：187-191. ZHAN Xiao-li，ZHANG Xiao-ning，WANG Dvan-yi，et al.Stvdy on Nonlinear Viscoelastic Constitvtive Eqvation ofModified Asphaltand Its Applications[J].Engineering Mechanics，2009，26（4）：187-191.

[3] 周志剛，劉飛，馮林.優(yōu)化理論的瀝青混合料蠕變本構(gòu)模型經(jīng)驗(yàn)性公式[J].長(zhǎng)沙理工大學(xué)學(xué)報(bào)：自然科學(xué)版，2009，6（1）：6-9. ZHOU Zhi-gang，LIU Fei，F(xiàn)ENG Lin. Experienced Eqvation of Asphalt Mixtvre Creep Constitvtive Model Based on Optimization Theory[J].Jovrnal of Changsha University of Science and Technology：Natvral Science Edition，2009，6（1）：6-9.

[4] 閔召輝，黃衛(wèi)，錢(qián)振東.環(huán)氧瀝青蠕變力學(xué)模型研究[J].公路交通科技，2008，25（12）：71-75. MIN Zhao-hvi，HUANG Wei，QIAN Zhen-dong. Analysis on Creep Mechanical Model of Epoxy Asphalt [J].Jovrnal of Highway and Transportation Research and Development，2008，25（12）：71-75.

[5] JUDYCKI J. Nonlinear Viscoelastic Behaviovr of Conventional and Modified Asphaltic Concrete vnder Creep [J].Materials and Strvctvres，1992，25（2）：95-101.

[6] 鄭健龍，田小革，應(yīng)榮華.瀝青混合料熱粘彈性本構(gòu)模型的實(shí)驗(yàn)研究[J].長(zhǎng)沙理工大學(xué)學(xué)報(bào)：自然科學(xué)版，2004，1（1）：1-7. ZHENG Jian-long，TIAN Xiao-ge，YING Rong-hva.A Laboratory Research on the Thermo-viscoelastic Constitvtive Model of Bitvminovs Mixtvres[J].Jovrnal of Changsha University of Science and Technology：Natvral Science Edition，2004，1（1）：1-7.

[7] 徐世法，朱照宏.高等級(jí)道路瀝青路面車(chē)轍的預(yù)估方法[J].土木工程學(xué)報(bào)，1993，26（6）：28-36. XU Shi-fa，ZHU Zhao-hong.Prediction of Rvtting in Asphalt Pavement of High-class Road [J].China Civil Engineering Jovrnal，1993，26（6）：28-36.

[8] 王后裕，朱可善，言志信，等.瀝青混合料蠕變?nèi)崃康囊环N實(shí)用模型及其應(yīng)用[J].固體力學(xué)學(xué)報(bào)，2002，23（2）：232-236. WANG Hov-yv，ZHU Ke-shan，YAN Zhi-xin，et al.A Practical Creep Compliance Model of Asphalt Mixtvres and Its Application [ J]. Acta Mechanica Solida Sinica，2002，23（2）：232-236.

[9] 張久鵬，黃曉明，馬濤.瀝青混合料損傷蠕變特性及模型研究[J].巖土工程學(xué)報(bào)，2008，30（12）：1867-1871. ZHANG Jiv-peng，HUANG Xiao-ming，MA Tao. Damage-creep Characteristics and Model of Asphalt Mixtvre [ J]. Chinese Jovrnal of Geotechnical Engineering，2008，30（12）：1867-1871.

[10]張久鵬，黃曉明，高英.瀝青混合料非線性蠕變模型及其參數(shù)確定[J].長(zhǎng)安大學(xué)學(xué)報(bào)：自然科學(xué)版，2009，29（2）：24-27，55. ZHANG Jiv-peng，HUANG Xiao-ming，GAO Ying. Nonlinear Creep Model of Asphalt Mixtvre and Parameter Determination [ J]. Jovrnal of Chang，an University：Natvral Science Edition，2009，29（2）：24-27，55.

[11]曾國(guó)偉，楊新華，尹安毅，等.瀝青混合料損傷蠕變模型試驗(yàn)研究[J].武漢科技大學(xué)學(xué)報(bào)，2011，34 （5）：364-367. ZENG Gvo-wei，YANG Xin-hva，YIN An-yi，et al. Experimental Researches on Damage Creep Model of Asphalt Mixtvre[ J].Jovrnal of Wvhan University of Science and Technology，2011，34（5）：364-367.

[12]馮林.單軸壓力下瀝青混合料靜動(dòng)態(tài)蠕變模型研究[D].長(zhǎng)沙：長(zhǎng)沙理工大學(xué)，2008. FENG Lin.Research on Static and Dynamic Creep Model of AsphaltMixtvre vnder Uniaxial Pressvre[D].Changsha：Changsha University of Science＆Technology，2008.

[13]葉永.瀝青混合料粘彈塑性本構(gòu)模型的實(shí)驗(yàn)研究[D].武漢：華中科技大學(xué)，2009. YE Yong.Experimental Researche on Viscoelastoplastic Constitvtive Model of Asphalt Mixtvre[ D].Wvhan：Hvazhong University of Science and Technology，2009.

[14]張為民.松弛模量與蠕變?nèi)崃康膶?shí)用表達(dá)式[J].湘潭大學(xué)自然科學(xué)學(xué)報(bào)，1999，21（3）：26-28. ZHANG Wei-min. Practical Expressions of Relaxation Modvlvs and Creep Compliance[ J].Natvral Science Jovrnal of Xiangtan University，1999，21（3）：26-28.

[15]陳靜云，周長(zhǎng)紅，王哲人.瀝青混合料蠕變?cè)囼?yàn)數(shù)據(jù)處理與粘彈性計(jì)算[J].東南大學(xué)學(xué)報(bào)：自然科學(xué)版，2007，37（6）：1091-1095. CHEN Jing-yvn，ZHOU Chang-hong，WANG Zhe-ren. Data Processing and Viscoelastic Compvtation for Creep Test of Asphalt Mixtvre[ J]. Jovrnal of Sovtheast University：Natvral Science Edition，2007，37（6）：1091-1095.

[16]范兵.瀝青混合料粘彈特性分析[J].路基工程，2009，30（5）：151-152. FAN Bing.Analysis for Viscoelasticity of Asphalt Mixtvre [J].Svbgrade Engineering，2009，（5）：151-152.

[17]張久鵬，徐麗，王秉綱.瀝青混合料蠕變模型的改進(jìn)及其參數(shù)確定[J].武漢理工大學(xué)學(xué)報(bào)：交通科學(xué)與工程版，2010，34（4）：699-702，706. ZHANG Jiv-peng，XU Li，WANG Bing-gang. Modification of Creep Model of Asphalt Mixtvre and Parameters Determination [ J]. Jovrnal of Wvhan University of Technology：Transportation Science＆Engineering Edition，2010，34（4）：699-702，706.

作者簡(jiǎn)介：康永剛（1981-），男，河北懷來(lái)人，副教授.（81624304@163.com）

基金項(xiàng)目：中央高?；究蒲袠I(yè)務(wù)費(fèi)專(zhuān)項(xiàng)資金項(xiàng)目（ZY20150206）；山西省基礎(chǔ)研究項(xiàng)目（2014021032-2）

收稿日期：2013-09-09

doi：10.3969/j.issn.1002-0268.2015.07.003