專利名稱：多重?zé)嵫h(huán)下渦輪葉片熱障涂層危險(xiǎn)區(qū)域的預(yù)測(cè)方法
技術(shù)領(lǐng)域：
本發(fā)明屬于高性能航空發(fā)動(dòng)機(jī)隔熱防護(hù)涂層系統(tǒng)可靠性預(yù)測(cè)技術(shù)領(lǐng)域，尤其涉及
一種多重?zé)嵫h(huán)下渦輪葉片熱障涂層危險(xiǎn)區(qū)域的預(yù)測(cè)方法。
背景技術(shù)：
熱障涂層(Thermal barrier coatings,簡(jiǎn)稱TBCs)是一種陶瓷材料，它具有熔點(diǎn) 高、熱傳導(dǎo)率低、蒸汽壓低、輻射率低和反射率高等特點(diǎn)，它是將陶瓷粉末噴涂或者沉積在 高溫合金熱端部件(尤其是渦輪葉片)表面，用以降低高溫部件的工作溫度，使其免受高溫 腐蝕和高溫氧化，從而達(dá)到延長(zhǎng)高溫部件的使用壽命的目的，滿足現(xiàn)代航空燃?xì)鉁u輪發(fā)動(dòng) 機(jī)內(nèi)高溫合金部件在高于其熔點(diǎn)溫度的服役環(huán)境中工作的要求，提高了航空發(fā)動(dòng)機(jī)燃?xì)鉁?度和熱效率。因此，它被廣泛應(yīng)用于航空航天、化工、冶金和能源等領(lǐng)域。
然而在實(shí)際應(yīng)用中，由于各層材料參數(shù)不匹配、高溫蠕變、高溫界面氧化和陶瓷材 料高溫相變等因素共同導(dǎo)致熱障涂層受到熱應(yīng)力和殘余壓縮應(yīng)力的交替作用，而且隨著應(yīng) 用時(shí)間的增加，陶瓷層內(nèi)受到越來(lái)越大的殘余壓縮應(yīng)力作用，同時(shí)還伴隨著涂層界面孔洞 或界面裂紋的不斷成核、擴(kuò)展和裂紋連接。隨著熱循環(huán)次數(shù)的增加，逐漸增大的殘余壓縮應(yīng) 力和界面裂紋共同導(dǎo)致了陶瓷涂層以屈曲和剝落形式與金屬基體相脫離而破壞。 一旦涂層 發(fā)生剝落，熱端金屬部件將直接暴露在高溫惡劣環(huán)境下，其后果是十分嚴(yán)重的。因此國(guó)內(nèi)外 許多研究人員采用理論分析、實(shí)驗(yàn)研究和模擬方法等多種手段來(lái)研究在特定工作環(huán)境下， 陶瓷涂層與金屬基底之間的破壞過(guò)程和破壞機(jī)理，從而預(yù)測(cè)熱障涂層系統(tǒng)的工作壽命或服 役時(shí)間，提高其應(yīng)用的可靠性。但是對(duì)于復(fù)雜結(jié)構(gòu)的熱障涂層系統(tǒng)(例如渦輪葉片、導(dǎo)向葉 片)，一般很難采用理論解析解來(lái)進(jìn)行相關(guān)研究，必須依靠實(shí)驗(yàn)測(cè)試和有限元模擬方法來(lái)實(shí) 現(xiàn)。其中，反復(fù)的實(shí)驗(yàn)測(cè)試需要花費(fèi)大量的人力、物力和財(cái)力。這就使得有限元模擬方法成 為一種主流的研究方法，該方法不僅可以降低試驗(yàn)成本，而且能夠減少設(shè)計(jì)和研制周期。
目前在熱障涂層領(lǐng)域中，已經(jīng)開(kāi)展了比較廣泛的有限元模擬研究，預(yù)測(cè)熱障涂層 系統(tǒng)在某工作環(huán)境下的溫度場(chǎng)、位移場(chǎng)、應(yīng)力場(chǎng)以及破壞過(guò)程。但是，大部分研究工作還主 要集中在基礎(chǔ)研究，其研究對(duì)象主要是簡(jiǎn)單的幾何構(gòu)型(平板模型、半圓形模型等)的熱障 涂層系統(tǒng)(部件)。對(duì)于復(fù)雜曲面實(shí)際形狀的渦輪葉片熱障涂層系統(tǒng)的有限元模擬，我們已 率先提出了一個(gè)實(shí)用的建模方法(公開(kāi)號(hào)CN101567025)。但是，當(dāng)需要模擬的熱循環(huán)次數(shù) (N)非常大時(shí)(例如N > 100次，甚至上萬(wàn)次)，前期處理的工作量變得非常大(在ABAQUS/ CAE中，500個(gè)循環(huán)的設(shè)置需要約為50 200小時(shí)的工作量，且極易出錯(cuò))，并且操作繁瑣，容 易出錯(cuò)。因此迫切需要發(fā)展相關(guān)有限元特殊前期處理手段或程序來(lái)實(shí)現(xiàn)。

發(fā)明內(nèi)容
本發(fā)明的目的在于，提供一種多重?zé)嵫h(huán)下渦輪葉片熱障涂層危險(xiǎn)區(qū)域的預(yù)測(cè)方 法，用以解決在使用ABAQUS/CAE進(jìn)行渦輪葉片熱障涂層系統(tǒng)大數(shù)目熱循環(huán)(N > 100)時(shí)， 前期處理的工作繁瑣復(fù)雜，過(guò)度消耗人力和時(shí)間的問(wèn)題。
3
技術(shù)方案是， 一種多重?zé)嵫h(huán)下渦輪葉片熱障涂層危險(xiǎn)區(qū)域的預(yù)測(cè)方法，在渦輪 葉片熱障涂層系統(tǒng)的有限元分析模型中，賦予材料參數(shù)，設(shè)置一個(gè)熱循環(huán)的分析步以及邊 界條件，對(duì)模型進(jìn)行網(wǎng)格劃分后，生成一個(gè)熱循環(huán)分析步的inp文件，其特征是所述方法還 包括 —對(duì)生成的所述一個(gè)熱循環(huán)分析步的inp文件進(jìn)行處理，去掉的所述一個(gè)熱循環(huán) 分析步的inp文件中的分隔符； —從處理后的所述一個(gè)熱循環(huán)分析步的inp文件中，提取出一個(gè)熱循環(huán)所對(duì)應(yīng)的 歷史數(shù)據(jù)； —將所述一個(gè)熱循環(huán)所對(duì)應(yīng)的歷史數(shù)據(jù)，進(jìn)行設(shè)定次數(shù)的熱循環(huán)，生成設(shè)定次數(shù) 的熱循環(huán)分析步數(shù)據(jù)； —修改設(shè)定次數(shù)的熱循環(huán)分析步的注釋名和名稱； —將經(jīng)過(guò)修改的所述設(shè)定次數(shù)的熱循環(huán)分析步數(shù)據(jù)覆蓋所述一個(gè)熱循環(huán)分析步 的inp文件中的數(shù)據(jù)，修改inp文件名并保存； —分析修改后的inp文件的數(shù)據(jù)，獲得渦輪葉片熱障涂層系統(tǒng)危險(xiǎn)區(qū)域的預(yù)測(cè)結(jié) 果。 所述一個(gè)熱循環(huán)分析步包括升溫過(guò)程、保濕過(guò)程及降溫過(guò)程。 所述設(shè)置邊界條件包括設(shè)置力學(xué)邊界條件和設(shè)置熱學(xué)邊界條件。所述渦輪葉片熱障涂層系統(tǒng)的有限元分析模型，通過(guò)CATIA軟件和ABAQUS軟件共
同獲得。 所述分析修改后的inp文件的數(shù)據(jù)，獲得渦輪葉片熱障涂層系統(tǒng)危險(xiǎn)區(qū)域的預(yù)測(cè) 結(jié)果，具體是采用ABAQUS軟件實(shí)現(xiàn)。 本發(fā)明僅需要設(shè)置一個(gè)熱循環(huán)，其余熱循環(huán)所需要的設(shè)置可以自動(dòng)完成，極大地 降低了前處理成本，降低人工操作的出錯(cuò)率，提高預(yù)測(cè)結(jié)果的準(zhǔn)確性，特別適用于大數(shù)目熱 循環(huán)的渦輪葉片熱障涂層系統(tǒng)危險(xiǎn)區(qū)域的預(yù)測(cè)。


圖1是本發(fā)明提供的多重?zé)嵫h(huán)下渦輪葉片熱障涂層危險(xiǎn)區(qū)域的預(yù)測(cè)方法流程 圖； 圖2是單個(gè)渦輪葉片熱障涂層系統(tǒng)的整體模型圖；其中，1是熱障涂層，2是氧化 層，3是粘結(jié)層，4是基底，5是底座； 圖3是熱障涂層系統(tǒng)承受的熱循環(huán)實(shí)驗(yàn)方式示意圖； 圖4是單個(gè)渦輪葉片熱障涂層系統(tǒng)在經(jīng)歷500個(gè)熱循環(huán)后，陶瓷層內(nèi)應(yīng)力云圖的 分布圖； 圖5是表示應(yīng)力集中區(qū)域每個(gè)熱循環(huán)最大應(yīng)力隨熱循環(huán)的演化示意圖。
具體實(shí)施例方式
下面結(jié)合附圖，對(duì)優(yōu)選實(shí)施例作詳細(xì)說(shuō)明。應(yīng)該強(qiáng)調(diào)的是，下述說(shuō)明僅僅是示例性 的，而不是為了限制本發(fā)明的范圍及其應(yīng)用。 選取單個(gè)渦輪葉片熱障涂層系統(tǒng)作為實(shí)施對(duì)象，進(jìn)行其在設(shè)定次數(shù)的熱循環(huán)(本實(shí)施例設(shè)定為500個(gè)熱循環(huán))作用下的有限元模擬，最終預(yù)測(cè)出陶瓷層內(nèi)危險(xiǎn)區(qū)域的分布。 圖1是本發(fā)明提供的多重?zé)嵫h(huán)下渦輪葉片熱障涂層危險(xiǎn)區(qū)域的預(yù)測(cè)方法流程
圖，圖l中，本發(fā)明包括以下步驟 步驟101 :建立有限元分析模型。 通過(guò)CATIA軟件與ABAQUS軟件，先建立渦輪葉片熱障涂層系統(tǒng)的有限元分析模 型，如圖2(a)所示。我們選取的各層厚度為he = 0. 40mm，ht = 0. 10mm，hb = 0. 20mm，hs = 1.80mm，如圖2(b)所示，渦輪葉片由外到內(nèi)依次為熱障涂層1、氧化層2、粘結(jié)層3、基底4。
為了方便、簡(jiǎn)化起見(jiàn)，本發(fā)明作如下假設(shè)1)各層涂層、基底及底座材料均勻且均 各向同性；2)各層涂層厚度均勻；3)底座簡(jiǎn)化為長(zhǎng)方體；4)采用理想彈塑性模型，且僅考 慮TBC層的高溫蠕變；5)僅僅考慮陶瓷層內(nèi)的破壞，不考慮界面之間的破壞；6)本實(shí)施例 只描了單個(gè)渦輪葉片熱障涂層系統(tǒng)的有限元仿真。
步驟102 :賦予材料參數(shù)，設(shè)置熱循環(huán)分析部和邊界條件。 分別對(duì)TBC層、TGO層、BC層、葉片基底定義材料參數(shù)，熱障涂層系統(tǒng)各層材料的關(guān) 鍵材料參數(shù)均考慮隨溫度變化而變化，具體數(shù)值如表1-7所示。 力學(xué)邊界條件設(shè)置1)左下邊u = w = 0 ;2)右下邊w = 0。如圖2 (a)所示。
熱學(xué)邊界條件是根據(jù)熱循環(huán)來(lái)設(shè)置TBC層外表面與冷卻通道的溫度，熱障涂層系 統(tǒng)所經(jīng)受的熱循環(huán)如圖4所示，主要包括2個(gè)過(guò)程 l)TBCs系統(tǒng)制備工序完成后，整個(gè)系統(tǒng)從626. 85°C (900K)冷卻至20°C ;
2)設(shè)定的熱循環(huán)方式是在IO分鐘內(nèi)陶瓷層表面升溫至112rC，渦輪葉片(基 底)內(nèi)冷卻通道升溫至70(TC;然后保持這種溫度梯度持續(xù)時(shí)間為40分鐘；再在10分鐘內(nèi) 整個(gè)系統(tǒng)均冷卻至2(TC。 分析步的設(shè)置采用隱式的熱_力耦合類型的分析步。在本步驟中，僅需要設(shè)置一
個(gè)熱循環(huán)的分析步即可，包括一個(gè)熱循環(huán)的各個(gè)過(guò)程，如升溫過(guò)程、保濕過(guò)程及降溫過(guò)程等。 熱循環(huán)的分析步命名規(guī)則如下
St印-lc-l表示第1個(gè)循環(huán)的升溫過(guò)程；
St印-lc-2表示第1個(gè)循環(huán)的保溫過(guò)程；
St印-lc-3表示第1個(gè)循環(huán)的降溫過(guò)程。 步驟103 :對(duì)模型進(jìn)行網(wǎng)格劃分，并生成只含一個(gè)熱循環(huán)的inp文件。 網(wǎng)格類型的選擇葉身與底座連接處選為C3D4T單元，其它各處選取C3D8RHT。這
樣可以既保證一定的計(jì)算精度，又能控制計(jì)算規(guī)模。總共使用了 49495個(gè)單元。 網(wǎng)格劃分完成之后，保存為Patents-02-cuboid-L15-900-lc. cae文件，然后輸出
Job-Patents-02-cuboid-L15-900-lc. inp文件(只含一個(gè)熱循環(huán))。其中，文件名中的"lc"
表示該inp文件僅含一個(gè)熱循環(huán)。 步驟104 :從上述inp文件中，提取出一個(gè)熱循環(huán)所對(duì)應(yīng)的歷史數(shù)據(jù)。 本發(fā)明使用由Visual C++6. 0自主開(kāi)發(fā)的"基于ABAQUS的渦輪葉片熱障涂層系統(tǒng)
大數(shù)目熱循環(huán)仿真系統(tǒng)"進(jìn)行inp文件處理。 把Job-Patents-02-cuboid-L15-900-lc. inp導(dǎo)入程序"基于ABAQUS的渦輪葉片 熱障涂層系統(tǒng)大數(shù)目熱循環(huán)仿真系統(tǒng)"，找到inp文件的歷史數(shù)據(jù)，把第一個(gè)循環(huán)所對(duì)應(yīng)的歷史數(shù)據(jù)提取出來(lái)，即下述部分 氺氺---------------------------- ** **STEP :St印-lc-l
** 氺St印，name = Step_lc_l， nlgeom ......
YES
**STEP :St印-lc-3氺氺
氺St印，name = Step_lc_3， nlgeom
YES
氺End Step
步驟105 :將上述一個(gè)熱循環(huán)所對(duì)應(yīng)的歷史數(shù)據(jù)，進(jìn)行設(shè)定次數(shù)的熱循環(huán)，生成設(shè)
定次數(shù)的熱循環(huán)分析步數(shù)據(jù)。設(shè)定次數(shù)的熱循環(huán)，在本實(shí)施例中為500次熱循環(huán)。設(shè)定的次數(shù)在"基于ABAQUS的渦輪葉片熱障涂層系統(tǒng)大數(shù)目熱循環(huán)仿真系統(tǒng)"中輸入。
步驟106 :修改設(shè)定次數(shù)的熱循環(huán)分析步的注釋名和名稱。 把生成的設(shè)定次數(shù)的熱循環(huán)分析步數(shù)據(jù)的注釋名與名字分別修改為對(duì)應(yīng)的名字，其它各數(shù)據(jù)保持不變，如 把第2個(gè)循環(huán)的升溫過(guò)程對(duì)應(yīng)的分析步改名為St印-2c-l
把第2個(gè)循環(huán)的保溫過(guò)程對(duì)應(yīng)的分析步改名為St印-2c-2
把第2個(gè)循環(huán)的降溫過(guò)程對(duì)應(yīng)的分析步改名為St印-2c-3
把第3個(gè)循環(huán)的升溫過(guò)程對(duì)應(yīng)的分析步改名為St印-3c-l
把第3個(gè)循環(huán)的保溫過(guò)程對(duì)應(yīng)的分析步改名為St印-3c-2
把第3個(gè)循環(huán)的降溫過(guò)程對(duì)應(yīng)的分析步改名為St印-3c-3
……其余依此類推。 步驟107 :將經(jīng)過(guò)修改的所述設(shè)定次數(shù)的熱循環(huán)分析步數(shù)據(jù)覆蓋所述一個(gè)熱循環(huán)分析步的inp文件中的數(shù)據(jù)，修改inp文件名并保存 在進(jìn)行500個(gè)熱循環(huán)后，將生成的數(shù)據(jù)覆蓋前述的一個(gè)熱循環(huán)分析步的i即文件中的數(shù)據(jù)，將文件名修改為Job-Patents-02-cuboid-L15-900-500c. i即。其中，文件名中的"500c"表示該inp文件含500個(gè)熱循環(huán)。這樣，我們就得到了模擬渦輪葉片熱障涂層系統(tǒng)500個(gè)熱循環(huán)所需要的inp文件。 步驟108 :將Job-Patents-02-cuboid-L15-900-500c. inp文件提交ABAQUS軟件進(jìn)行分析，獲得渦輪葉片熱障涂層系統(tǒng)危險(xiǎn)區(qū)域的預(yù)測(cè)結(jié)果。 模擬的結(jié)果可得到渦輪葉片熱障涂層系統(tǒng)的溫度場(chǎng)、位移場(chǎng)、應(yīng)力應(yīng)變場(chǎng)等關(guān)鍵物理量。由于陶瓷層的主要成分為氧化鋯，屬于脆性材料，我們主要利用第一強(qiáng)度理論來(lái)評(píng)價(jià)該層的應(yīng)力分布情況。因此，我們主要觀察最大主應(yīng)力(Max. Principal stress)的分布，最大主應(yīng)力集中之處就是可能最先發(fā)生失效的地方，即危險(xiǎn)區(qū)域，危險(xiǎn)區(qū)域之外的區(qū)域的應(yīng)力水平一般都相對(duì)較低， 一般不會(huì)發(fā)生破壞，因此為安全區(qū)域。 本發(fā)明分析處理過(guò)程基于ABAQUS/CAE軟件完成，操作方便，不需要利用其它編譯
6軟件額外編程實(shí)現(xiàn)；而且，在ABAQUS/CAE中僅需要設(shè)置一個(gè)熱循環(huán)即可，其余熱循環(huán)所需要的設(shè)置可以自動(dòng)完成，極大地降低了前處理成本(在ABAQUS/CAE中，500個(gè)循環(huán)的設(shè)置約為50 200小時(shí)的工作量，且極易出錯(cuò))，使用本發(fā)明則避免了人工出錯(cuò)，準(zhǔn)確率高，特別適用于大數(shù)目熱循環(huán)(循環(huán)數(shù)目N> 100)的處理。另外，本方法可以很方便地得到每個(gè)熱循環(huán)關(guān)鍵時(shí)間點(diǎn)的物理量信息，如升溫結(jié)束時(shí)刻、保溫結(jié)束時(shí)刻、冷卻結(jié)束時(shí)刻的應(yīng)力值等物理量信息。這樣就可以使用更長(zhǎng)的分析步長(zhǎng)，減少計(jì)算時(shí)間，降低計(jì)算成本。最后，本方法同樣適用于無(wú)涂層的渦輪葉片(發(fā)動(dòng)機(jī))或理想形狀的渦輪葉片熱障涂層系統(tǒng)的大數(shù)目熱循環(huán)的模擬。 以上所述，僅為本發(fā)明較佳的具體實(shí)施方式
，但本發(fā)明的保護(hù)范圍并不局限于此，任何熟悉本技術(shù)領(lǐng)域的技術(shù)人員在本發(fā)明揭露的技術(shù)范圍內(nèi)，可輕易想到的變化或替換，都應(yīng)涵蓋在本發(fā)明的保護(hù)范圍之內(nèi)。因此，本發(fā)明的保護(hù)范圍應(yīng)該以權(quán)利要求的保護(hù)范圍為準(zhǔn)。
權(quán)利要求
一種多重?zé)嵫h(huán)下渦輪葉片熱障涂層危險(xiǎn)區(qū)域的預(yù)測(cè)方法，在渦輪葉片熱障涂層系統(tǒng)的有限元分析模型中，賦予材料參數(shù)，設(shè)置一個(gè)熱循環(huán)的分析步以及邊界條件，對(duì)模型進(jìn)行網(wǎng)格劃分后，生成一個(gè)熱循環(huán)分析步的inp文件，其特征是所述方法還包括--對(duì)生成的所述一個(gè)熱循環(huán)分析步的inp文件進(jìn)行處理，去掉的所述一個(gè)熱循環(huán)分析步的inp文件中的分隔符；--從處理后的所述一個(gè)熱循環(huán)分析步的inp文件中，提取出一個(gè)熱循環(huán)所對(duì)應(yīng)的歷史數(shù)據(jù)；--將所述一個(gè)熱循環(huán)所對(duì)應(yīng)的歷史數(shù)據(jù)，進(jìn)行設(shè)定次數(shù)的熱循環(huán)，生成設(shè)定次數(shù)的熱循環(huán)分析步數(shù)據(jù)；--修改設(shè)定次數(shù)的熱循環(huán)分析步的注釋名和名稱；--將經(jīng)過(guò)修改的所述設(shè)定次數(shù)的熱循環(huán)分析步數(shù)據(jù)覆蓋所述一個(gè)熱循環(huán)分析步的inp文件中的數(shù)據(jù)，修改inp文件名并保存；--分析修改后的inp文件的數(shù)據(jù)，獲得渦輪葉片熱障涂層系統(tǒng)危險(xiǎn)區(qū)域的預(yù)測(cè)結(jié)果。
2. 根據(jù)權(quán)利要求1所述的一種多重?zé)嵫h(huán)下渦輪葉片熱障涂層危險(xiǎn)區(qū)域的預(yù)測(cè)方法， 其特征是所述一個(gè)熱循環(huán)分析步包括升溫過(guò)程、保濕過(guò)程及降溫過(guò)程。
3. 根據(jù)權(quán)利要求1所述的一種多重?zé)嵫h(huán)下渦輪葉片熱障涂層危險(xiǎn)區(qū)域的預(yù)測(cè)方法， 其特征是所述設(shè)置邊界條件包括設(shè)置力學(xué)邊界條件和設(shè)置熱學(xué)邊界條件。
4. 根據(jù)權(quán)利要求1所述的一種多重?zé)嵫h(huán)下渦輪葉片熱障涂層危險(xiǎn)區(qū)域的 預(yù)測(cè) 方法，其特征是所述渦輪葉片熱障涂層系統(tǒng)的有限元分析模型，通過(guò)CATIA軟件和ABAQUS 軟件共同獲得。
5. 根據(jù)權(quán)利要求1所述的一種多重?zé)嵫h(huán)下渦輪葉片熱障涂層危險(xiǎn)區(qū)域的預(yù)測(cè)方法， 其特征是所述分析修改后的inp文件的數(shù)據(jù)，獲得渦輪葉片熱障涂層系統(tǒng)危險(xiǎn)區(qū)域的預(yù)測(cè) 結(jié)果，具體是采用ABAQUS軟件實(shí)現(xiàn)。
全文摘要
本發(fā)明公開(kāi)了高性能航空發(fā)動(dòng)機(jī)隔熱防護(hù)涂層系統(tǒng)可靠性預(yù)測(cè)技術(shù)領(lǐng)域中的一種多重?zé)嵫h(huán)下渦輪葉片熱障涂層危險(xiǎn)區(qū)域的預(yù)測(cè)方法。包括對(duì)生成的所述一個(gè)熱循環(huán)分析步的inp文件進(jìn)行處理；從處理后的的inp文件中，提取出一個(gè)熱循環(huán)所對(duì)應(yīng)的歷史數(shù)據(jù)；進(jìn)行設(shè)定次數(shù)的熱循環(huán)，生成設(shè)定次數(shù)的熱循環(huán)分析步數(shù)據(jù)；修改設(shè)定次數(shù)的熱循環(huán)分析步的注釋名和名稱；將經(jīng)過(guò)修改的所述設(shè)定次數(shù)的熱循環(huán)分析步數(shù)據(jù)覆蓋所述一個(gè)熱循環(huán)分析步的inp文件中的數(shù)據(jù)，修改inp文件名并保存；分析修改后的inp文件的數(shù)據(jù)，獲得渦輪葉片熱障涂層系統(tǒng)危險(xiǎn)區(qū)域的預(yù)測(cè)結(jié)果。本發(fā)明降低了前處理成本和人工操作出錯(cuò)率，提高預(yù)測(cè)結(jié)果的準(zhǔn)確性。
文檔編號(hào)F01D5/28GK101777086SQ201010034048
公開(kāi)日2010年7月14日 申請(qǐng)日期2010年1月12日 優(yōu)先權(quán)日2010年1月12日
發(fā)明者劉奇星, 周益春, 毛衛(wèi)國(guó) 申請(qǐng)人:湘潭大學(xué)