一種高速?zèng)_壓裝備關(guān)鍵部件熱態(tài)特性的跨尺度仿真方法
【專利摘要】本發(fā)明公開了一種高速?zèng)_壓裝備關(guān)鍵部件熱態(tài)特性的跨尺度仿真方法。該方法首先確定高速?zèng)_壓裝備關(guān)鍵部件在工作過程中的主要熱源及其發(fā)熱量。然后確定影響關(guān)鍵部件溫度場(chǎng)分布的熱傳遞途徑與熱傳遞方式，計(jì)算對(duì)流換熱系數(shù)。獲取關(guān)鍵部件固體粗糙表面的分形維數(shù)和特征長(zhǎng)度，模擬固體粗糙表面，綜合考慮彈性、塑性、彈塑性三種不同變形機(jī)制，計(jì)算其接觸熱阻，從而得到各類固體接觸面的接觸傳熱系數(shù)。最后對(duì)高速?zèng)_壓裝備關(guān)鍵部件進(jìn)行熱態(tài)特性有限元仿真分析。本發(fā)明考慮多變形機(jī)制，在考慮動(dòng)摩擦因素情況下，將變形分為三種情況分別計(jì)算基體熱阻，使得仿真方法對(duì)高速?zèng)_壓裝備關(guān)鍵部件熱態(tài)特性的分析更加全面、精確，更加符合實(shí)際情況。
【專利說明】
一種高速?zèng)_壓裝備關(guān)鍵部件熱態(tài)特性的跨尺度仿真方法
技術(shù)領(lǐng)域
[0001] 本發(fā)明涉及高速?zèng)_壓技術(shù)領(lǐng)域，尤其涉及一種高速?zèng)_壓裝備關(guān)鍵部件熱態(tài)特性的 跨尺度仿真方法。
【背景技術(shù)】
[0002] 高速?zèng)_壓裝備關(guān)鍵部件熱態(tài)特性分析主要有試驗(yàn)法和仿真方法。試驗(yàn)法需進(jìn)行大 量實(shí)驗(yàn)來完成整機(jī)溫度場(chǎng)的測(cè)量，溫度信息采集時(shí)，無法整機(jī)都布置傳感器，而通常是在沖 壓裝備某些部位進(jìn)行布置，傳感器的位置直接影響溫度信息的準(zhǔn)確釆集，最終影響熱分析 模型的精度。此外試驗(yàn)法需要的設(shè)備較多，成本較昂貴，不具備通用性。仿真方法借助現(xiàn)有 的有限元軟件對(duì)高速?zèng)_壓裝備關(guān)鍵部件進(jìn)行仿真建模與分析計(jì)算，能夠高效低成本地獲得 高速?zèng)_壓裝備關(guān)鍵部件內(nèi)部溫度場(chǎng)的變化規(guī)律與瞬態(tài)熱平衡時(shí)間，在工程中得到廣泛應(yīng) 用。
[0003] 要保證高速?zèng)_壓裝備關(guān)鍵部件熱態(tài)特性仿真分析結(jié)果的有效性，其關(guān)鍵在于準(zhǔn)確 計(jì)算熱態(tài)特性仿真分析時(shí)各類熱環(huán)境參數(shù)，構(gòu)建與實(shí)際相符的熱態(tài)特性仿真分析模型，其 中，固體接觸面接觸傳熱參數(shù)即接觸熱阻的準(zhǔn)確計(jì)算是難點(diǎn)。接觸熱阻作為決定固體接觸 面間熱傳遞的一個(gè)重要參數(shù)，對(duì)高速?zèng)_壓裝備關(guān)鍵部件的熱態(tài)特性有著直接而重要的影 響，因此，高速?zèng)_壓裝備關(guān)鍵部件的熱態(tài)特性仿真分析時(shí)必須充分考慮固體接觸面接觸熱 阻的影響。
[0004] 高速?zèng)_壓裝備運(yùn)行過程中，滾動(dòng)軸承、滑動(dòng)軸承連接處作為主要熱源部分，其接觸 熱阻對(duì)溫度場(chǎng)的截流作用更為明顯，需充分考慮其接觸面的微觀形貌特征和彈性、塑性、彈 塑性三種形變方式，建立與實(shí)際相符的固體接觸面的接觸熱阻計(jì)算模型，準(zhǔn)確計(jì)算各接觸 面的接觸傳熱參數(shù)，從而保證高速?zèng)_壓裝備關(guān)鍵部件熱態(tài)特性仿真分析結(jié)果的準(zhǔn)確性。

【發(fā)明內(nèi)容】

[0005] 本發(fā)明的目的在于針對(duì)現(xiàn)有技術(shù)的不足，提供一種高速?zèng)_壓裝備關(guān)鍵部件熱態(tài)特 性的跨尺度仿真方法。首先根據(jù)高速?zèng)_壓裝備關(guān)鍵部件的具體運(yùn)行情況確定該關(guān)鍵部件在 工作過程中的主要熱源，并計(jì)算其發(fā)熱量。然后根據(jù)各部件接觸面的具體接觸類型確定熱 傳遞方式，計(jì)算強(qiáng)制對(duì)流換熱系數(shù)和自然對(duì)流換熱系數(shù)。在考慮固體接觸面真實(shí)微觀形貌 與多邊形機(jī)制的情況下，給出各類固體接觸面的接觸熱阻計(jì)算公式，從而得到相應(yīng)的接觸 傳熱系數(shù)。最后利用Pro/E建立關(guān)鍵部件三維模型，導(dǎo)入ANSYS workbench中，劃分網(wǎng)格，設(shè) 置材料屬性和計(jì)算得到的熱環(huán)境參數(shù)，對(duì)高速?zèng)_壓裝備關(guān)鍵部件的熱態(tài)特性進(jìn)行有限元仿 真分析。
[0006] 為實(shí)現(xiàn)上述目的，本發(fā)明采用的技術(shù)方案是:一種高速?zèng)_壓裝備關(guān)鍵部件熱態(tài)特 性的跨尺度仿真方法，該方法包括以下步驟：
[0007] 1)確定影響高速?zèng)_壓裝備關(guān)鍵部件溫度場(chǎng)分布的各主要熱源，并計(jì)算各主要熱源 的發(fā)熱量。
[0008] 2)確定影響高速?zèng)_壓裝備關(guān)鍵部件溫度場(chǎng)分布的熱傳遞途徑與熱傳遞方式，計(jì)算 出對(duì)流換熱系數(shù)。
[0009] 3)細(xì)觀尺度下考慮多變形機(jī)制的固體接觸粗糙表面接觸傳熱系數(shù)計(jì)算，具體包括 以下子步驟：
[0010] 3.1)利用功率譜法與W-M分形函數(shù)獲取高速?zèng)_壓裝備關(guān)鍵部件固體接觸粗糙表面 的分形維數(shù)D和特征長(zhǎng)度G，利用分形-蒙特卡羅法模擬沖壓裝備關(guān)鍵部件的固體接觸粗糙 表面;所述粗糙表面由微凸峰不規(guī)則排列構(gòu)成，所述微凸峰由微凸體經(jīng)自仿射變換構(gòu)造新 微凸體堆疊，再由新生成微凸體經(jīng)自仿射變換構(gòu)造新微凸體堆疊，如此往復(fù)，直到新生成的 微凸體的微接觸直徑小于設(shè)定值；
[0011] 3.2)計(jì)算關(guān)鍵部件固體接觸粗糙表面上構(gòu)成微凸峰的微凸體的收縮熱阻和該微 凸峰的不同面積接觸點(diǎn)間熱阻；
[0012] 3.3)將關(guān)鍵部件固體接觸粗糙表面上構(gòu)成微凸峰的微凸體，根據(jù)其微接觸半徑的 大小，分成彈性變形微凸體、彈塑性變形微凸體和塑性變形微凸體;具體分類方法如下:計(jì) 算彈性變形的臨界微接觸半徑和塑性變形的臨界微接觸半徑;將兩個(gè)臨界微接觸半徑作為 分段點(diǎn)，若微凸體的微接觸半徑小于彈性變形的臨界微接觸半徑，則將其歸類為彈性變形 微凸體;若微凸體的微接觸半徑大于塑性變形的臨界微接觸半徑，則將其歸類為塑性變形 微凸體;其余歸為彈塑性變形微凸體；
[0013] 3.4)針對(duì)發(fā)生不同形變的微凸體，分別計(jì)算其基體熱阻；
[0014] 3.5)將固體接觸粗糙表面上構(gòu)成微凸峰的微凸體的基體熱阻、收縮熱阻和該微凸 峰的不同面積接觸點(diǎn)間熱阻串聯(lián)，得到高速?zèng)_壓裝備關(guān)鍵部件固體接觸粗糙表面上各微凸 峰的接觸熱阻；
[0015] 3.6)將固體接觸粗糙表面上所有的微凸峰的接觸熱阻并聯(lián)，計(jì)算得到關(guān)鍵部件固 體接觸粗糙表面的接觸熱阻。
[0016] 3.7)將關(guān)鍵部件固體接觸粗糙表面的接觸熱阻計(jì)算結(jié)果取倒數(shù)，得到固體接觸粗 糙表面的接觸傳熱系數(shù)。
[0017] 4)宏觀尺度下高速?zèng)_壓裝備關(guān)鍵部件熱態(tài)特性仿真分析:利用Pro/E構(gòu)建關(guān)鍵部 件的三維模型，導(dǎo)入ANSYS workbench，并劃分網(wǎng)格，設(shè)置部件中各零件的材料屬性。根據(jù)步 驟1)、2)、3)的計(jì)算結(jié)果，設(shè)置關(guān)鍵部件熱態(tài)特性仿真所需的發(fā)熱量、對(duì)流換熱系數(shù)和接觸 傳熱系數(shù)，在宏觀尺度下對(duì)高速?zèng)_壓裝備關(guān)鍵部件進(jìn)行熱態(tài)特性仿真分析。
[0018] 進(jìn)一步地，基于多變形機(jī)制的彈性形變的臨界微接觸半徑心和塑性變形的臨界微 接觸半徑R2的計(jì)算公式如下：
[0021 ]式中：〇y為材料屈服強(qiáng)度;K是兩接觸材料各自硬度與屈服強(qiáng)度的比值的較小值；1 為微凸體變形前截面橫向?qū)挾?E為有效彈性模量;G為特征長(zhǎng)度;D為分形維數(shù)。
[0022] 本發(fā)明的有益效果是:在考慮動(dòng)摩擦因素的情況下，綜合考慮了微凸體的彈性、彈 塑性、完全塑性三種變形情況下基體熱阻應(yīng)有不同的計(jì)算方式。同時(shí)，接觸面上的熱阻網(wǎng)絡(luò) 模型綜合考慮了基體熱阻、收縮熱阻以及不同面積接觸點(diǎn)之間的熱阻，使計(jì)算模型更全面， 使得所提出仿真方法對(duì)高速?zèng)_壓裝備關(guān)鍵部件熱態(tài)特性的分析更加全面、精確，更加符合 實(shí)際情況。
【附圖說明】
[0023] 圖1為高速?zèng)_壓裝備關(guān)鍵部件熱態(tài)特性的跨尺度仿真流程圖；
[0024]圖2為高速?zèng)_壓裝備施力機(jī)構(gòu)1/2簡(jiǎn)化對(duì)稱結(jié)構(gòu)圖；
[0025] 圖3為用分形-蒙特卡羅法產(chǎn)生粗糙表面的流程圖；
[0026] 圖4為跨尺度仿真方法所得的高速?zèng)_壓裝備施力機(jī)構(gòu)的溫度場(chǎng)變化趨勢(shì)圖。
【具體實(shí)施方式】
[0027] 以下結(jié)合附圖和具體實(shí)施例對(duì)本發(fā)明作進(jìn)一步說明。
[0028] 高速?zèng)_壓裝備關(guān)鍵部件熱態(tài)特性的跨尺度仿真流程如圖1所示。以高速?zèng)_壓裝備 施力機(jī)構(gòu)熱態(tài)特性仿真分析為例，說明高速?zèng)_壓裝備關(guān)鍵部件熱態(tài)特性跨尺度仿真的主要 步驟。高速?zèng)_壓裝備施力機(jī)構(gòu)主要由主軸、連桿、上橫梁、滑塊、軸承等構(gòu)成。圖2是本發(fā)明相 應(yīng)的高速?zèng)_壓裝備施力機(jī)構(gòu)1/2簡(jiǎn)化對(duì)稱結(jié)構(gòu)圖。仿真方法具體如下：
[0029] 1)高速?zèng)_壓裝備施力機(jī)構(gòu)熱源分析與發(fā)熱量計(jì)算；
[0030] 根據(jù)高速?zèng)_壓裝備的具體運(yùn)行情況確定影響高速?zèng)_壓裝備施力機(jī)構(gòu)溫度場(chǎng)分布 的主要熱源來自滑動(dòng)軸承發(fā)熱量和滾動(dòng)軸承發(fā)熱量。
[0031] 上橫梁與主軸的軸瓦、連桿與主軸的軸瓦實(shí)際均是滑動(dòng)軸承，與滾動(dòng)軸承相似，滑 動(dòng)軸承工作時(shí)，軸承的熱量由摩擦功轉(zhuǎn)變而來。軸承每秒產(chǎn)生的熱量可利用軸承的摩擦因 數(shù)，軸承所受載荷和軸承圓周線速度計(jì)算求得。
[0032] 滾動(dòng)軸承發(fā)熱由外加載荷摩擦力矩引起的滾動(dòng)體與滾道的摩擦生熱，由潤(rùn)滑劑粘 度力矩引起的克服滾動(dòng)體粘性阻力而產(chǎn)生的粘性摩擦生熱和由滾動(dòng)體在與內(nèi)外圈滾道接 觸區(qū)上發(fā)生自旋產(chǎn)生的自旋摩擦生熱所構(gòu)成，可分別計(jì)算各摩擦熱從而得到滾動(dòng)軸承的發(fā) 熱量。
[0033] 高速?zèng)_壓裝備施力機(jī)構(gòu)各熱源發(fā)熱量計(jì)算結(jié)果如表1所示，W為單位一一瓦特。
[0034] 表 1
[0035]
[0036] 2)高速?zèng)_壓裝備施力機(jī)構(gòu)熱傳遞分析與對(duì)流換熱系數(shù)計(jì)算；
[0037] 2.1)高速?zèng)_壓裝備施力機(jī)構(gòu)熱傳遞分析
[0038]確定影響高速?zèng)_壓裝備施力機(jī)構(gòu)溫度場(chǎng)分布的熱傳遞途徑與熱傳遞方式，計(jì)算出 對(duì)流換熱系數(shù)。
[0039]熱量的傳遞是由溫差引起的，在溫差作用下，熱量自發(fā)地從高溫物體傳向低溫物 體，實(shí)現(xiàn)熱量的傳遞。高速?zèng)_壓裝備施力機(jī)構(gòu)主要熱源的熱傳遞方式如表2所示。
[0042] 2.2)高速?zèng)_壓裝備施力機(jī)構(gòu)對(duì)流換熱系數(shù)計(jì)算
[0043] 高速?zèng)_壓裝備施力機(jī)構(gòu)對(duì)流換熱主要包括:主軸上軸承接觸部分模擬潤(rùn)滑油的強(qiáng) 制熱對(duì)流、主軸與空氣的強(qiáng)制熱對(duì)流、高速?zèng)_壓裝備外殼與空氣的自然熱對(duì)流。熱對(duì)流遵循 牛頓冷卻公式，其大小由對(duì)流換熱系數(shù)表示。由雷諾數(shù)判定冷卻液流動(dòng)狀態(tài)，并依據(jù)努塞爾 數(shù)可計(jì)算出主軸不同部件的對(duì)流換熱系數(shù)。
[0044] 強(qiáng)制流換熱系數(shù)可以通過部件周圍流體的導(dǎo)熱系數(shù)、努塞爾數(shù)、軸承所在軸頸的 直徑或者主軸當(dāng)量直徑進(jìn)行求解。
[0045] 自然流換熱系數(shù)可以通過格拉曉夫數(shù)和普朗特?cái)?shù)進(jìn)行求解。
[0046] 高速?zèng)_壓裝備施力機(jī)構(gòu)對(duì)流換熱系數(shù)計(jì)算結(jié)果如表3所示，W、m和K為單位一一瓦 特、米和開爾文。
[0047] 表 3
[0049] 3)細(xì)觀尺度下考慮多變形機(jī)制的固體接觸粗糙表面接觸傳熱系數(shù)計(jì)算，具體包括 以下子步驟：
[0050] 3.1)利用功率譜法與W-M分形函數(shù)獲取高速?zèng)_壓裝備關(guān)鍵部件固體接觸粗糙表面 的分形維數(shù)D和特征長(zhǎng)度G，利用分形-蒙特卡羅法模擬沖壓裝備關(guān)鍵部件的固體接觸粗糙 表面;所述粗糙表面由微凸峰不規(guī)則排列構(gòu)成，所述微凸峰由微凸體經(jīng)自仿射變換構(gòu)造新 微凸體堆疊，再由新生成微凸體經(jīng)自仿射變換構(gòu)造新微凸體堆疊，如此往復(fù)，直到新生成的 微凸體的微接觸直徑小于設(shè)定值；
[0051] 3.2)計(jì)算關(guān)鍵部件固體接觸粗糙表面上構(gòu)成微凸峰的微凸體的收縮熱阻和該微 凸峰的不同面積接觸點(diǎn)間熱阻；
[0052] 3.3)將關(guān)鍵部件固體接觸粗糙表面上構(gòu)成微凸峰的微凸體，根據(jù)其微接觸半徑的 大小，分成彈性變形微凸體、彈塑性變形微凸體和塑性變形微凸體;具體分類方法如下:計(jì) 算彈性變形的臨界微接觸半徑心和塑性變形的臨界微接觸半徑R2,計(jì)算公式如下：
[0055] 式中：〇y為材料屈服強(qiáng)度;K是兩接觸材料各自硬度與屈服強(qiáng)度的比值的較小值；1 為微凸體變形前截面橫向?qū)挾?Ε為有效彈性模量;G為特征長(zhǎng)度;D為分形維數(shù)。
[0056] 將兩個(gè)臨界微接觸半徑作為分段點(diǎn)，若微凸體的微接觸半徑小于彈性變形的臨界 微接觸半徑，則將其歸類為彈性變形微凸體;若微凸體的微接觸半徑大于塑性變形的臨界 微接觸半徑，則將其歸類為塑性變形微凸體;其余歸為彈塑性變形微凸體；
[0057] 3.4)針對(duì)發(fā)生不同形變的微凸體，分別計(jì)算其基體熱阻；
[0058] 3.5)將固體接觸粗糙表面上構(gòu)成微凸峰的微凸體的基體熱阻、收縮熱阻和該微凸 峰的不同面積接觸點(diǎn)間熱阻串聯(lián)，得到高速?zèng)_壓裝備關(guān)鍵部件固體接觸粗糙表面上各微凸 峰的接觸熱阻；
[0059] 3.6)將固體接觸粗糙表面上所有的微凸峰的接觸熱阻并聯(lián)，計(jì)算得到關(guān)鍵部件固 體接觸粗糙表面的接觸熱阻。
[0060] 3.7)將關(guān)鍵部件固體接觸粗糙表面的接觸熱阻計(jì)算結(jié)果取倒數(shù)，得到固體接觸粗 糙表面的接觸傳熱系數(shù)，計(jì)算結(jié)果如表4所示，W、m和Κ為單位--瓦特、米和開爾文。
[0061] 表 4
[0062]
[0063] 4)宏觀尺度下高速?zèng)_壓裝備施力機(jī)構(gòu)熱態(tài)特性仿真分析:利用Pro/E構(gòu)建高速?zèng)_ 壓裝備施力機(jī)構(gòu)的三維模型，將其導(dǎo)入ANSYS workbench中，根據(jù)需要確定網(wǎng)格大小和劃分 方式。網(wǎng)格劃分后，設(shè)置材料屬性，并根據(jù)步驟2)和3)計(jì)算結(jié)果，設(shè)置該熱態(tài)特性仿真模型 各熱源的發(fā)熱量、對(duì)流換熱系數(shù)和接觸傳熱系數(shù)，對(duì)高速?zèng)_壓裝備施力機(jī)構(gòu)熱態(tài)特性進(jìn)行 仿真分析，仿真結(jié)果如圖4所示。
【主權(quán)項(xiàng)】
1. 一種高速?zèng)_壓裝備關(guān)鍵部件熱態(tài)特性的跨尺度仿真方法，其特征在于，該方法包括 以下步驟： 1) 確定影響高速?zèng)_壓裝備關(guān)鍵部件溫度場(chǎng)分布的各主要熱源，并計(jì)算各主要熱源的發(fā) 熱量。 2) 確定影響高速?zèng)_壓裝備關(guān)鍵部件溫度場(chǎng)分布的熱傳遞途徑與熱傳遞方式，計(jì)算出對(duì) 流換熱系數(shù)。 3) 細(xì)觀尺度下考慮多變形機(jī)制的固體接觸粗糙表面接觸傳熱系數(shù)計(jì)算，具體包括以下 子步驟： 3.1) 利用功率譜法與W-M分形函數(shù)獲取高速?zèng)_壓裝備關(guān)鍵部件固體接觸粗糙表面的分 形維數(shù)D和特征長(zhǎng)度G，利用分形-蒙特卡羅法模擬沖壓裝備關(guān)鍵部件的固體接觸粗糙表面； 所述粗糙表面由微凸峰不規(guī)則排列構(gòu)成，所述微凸峰由微凸體經(jīng)自仿射變換構(gòu)造新微凸體 堆疊，再由新生成微凸體經(jīng)自仿射變換構(gòu)造新微凸體堆疊，如此往復(fù)，直到新生成的微凸體 的微接觸直徑小于設(shè)定值； 3.2) 計(jì)算關(guān)鍵部件固體接觸粗糙表面上構(gòu)成微凸峰的微凸體的收縮熱阻和該微凸峰 的不同面積接觸點(diǎn)間熱阻； 3.3) 將關(guān)鍵部件固體接觸粗糙表面上構(gòu)成微凸峰的微凸體，根據(jù)其微接觸半徑的大 小，分成彈性變形微凸體、彈塑性變形微凸體和塑性變形微凸體;具體分類方法如下:計(jì)算 彈性變形的臨界微接觸半徑和塑性變形的臨界微接觸半徑;將兩個(gè)臨界微接觸半徑作為分 段點(diǎn)，若微凸體的微接觸半徑小于彈性變形的臨界微接觸半徑，則將其歸類為彈性變形微 凸體;若微凸體的微接觸半徑大于塑性變形的臨界微接觸半徑，則將其歸類為塑性變形微 凸體;其余歸為彈塑性變形微凸體； 3.4) 針對(duì)發(fā)生不同形變的微凸體，分別計(jì)算其基體熱阻； 3.5) 將固體接觸粗糙表面上構(gòu)成微凸峰的微凸體的基體熱阻、收縮熱阻和該微凸峰的 不同面積接觸點(diǎn)間熱阻串聯(lián)，得到高速?zèng)_壓裝備關(guān)鍵部件固體接觸粗糙表面上各微凸峰的 接觸熱阻； 3.6) 將固體接觸粗糙表面上所有的微凸峰的接觸熱阻并聯(lián)，計(jì)算得到關(guān)鍵部件固體接 觸粗糙表面的接觸熱阻。 3.7) 將關(guān)鍵部件固體接觸粗糙表面的接觸熱阻計(jì)算結(jié)果取倒數(shù)，得到固體接觸粗糙表 面的接觸傳熱系數(shù)。 4) 宏觀尺度下高速?zèng)_壓裝備關(guān)鍵部件熱態(tài)特性仿真分析:利用Pro/E構(gòu)建關(guān)鍵部件的 三維模型，導(dǎo)入ANSYS workbench，并劃分網(wǎng)格，設(shè)置部件中各零件的材料屬性。根據(jù)步驟 1)、2)、3)的計(jì)算結(jié)果，設(shè)置關(guān)鍵部件熱態(tài)特性仿真所需的發(fā)熱量、對(duì)流換熱系數(shù)和接觸傳 熱系數(shù)，在宏觀尺度下對(duì)高速?zèng)_壓裝備關(guān)鍵部件進(jìn)行熱態(tài)特性仿真分析。2. 根據(jù)權(quán)利要求1所述的高速?zèng)_壓裝備關(guān)鍵部件熱態(tài)特性的跨尺度仿真方法，其特征 在于，基于多變形機(jī)制的彈性形變的臨界微接觸半徑心和塑性變形的臨界微接觸半徑R 2的 計(jì)算公式如下：式中：〇y為材料屈服強(qiáng)度;K為兩接觸材料各自硬度與屈服強(qiáng)度的比值的較小值；1為微 凸體變形前截面橫向?qū)挾?Ε為有效彈性模量;G為特征長(zhǎng)度;D為分形維數(shù)。
【文檔編號(hào)】G06F17/50GK105868446SQ201610172556
【公開日】2016年8月17日
【申請(qǐng)日】2016年3月24日
【發(fā)明人】程錦, 周振棟, 林志強(qiáng), 馮毅雄, 劉振宇, 譚建榮
【申請(qǐng)人】浙江大學(xué)