亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        采用每股流體必進(jìn)化策略的換熱網(wǎng)絡(luò)優(yōu)化方法

        2020-08-07 09:54:36崔國(guó)民陳家星
        關(guān)鍵詞:算例流體個(gè)體

        曹 美, 崔國(guó)民, 陳家星, 李 劍

        (上海理工大學(xué) 新能源科學(xué)與工程研究所, 上海 200093)

        1 前 言

        換熱網(wǎng)絡(luò)綜合已成為化工系統(tǒng)研究的難點(diǎn),國(guó)內(nèi)外相關(guān)研究從優(yōu)化方法的角度主要分為夾點(diǎn)分析和數(shù)學(xué)規(guī)劃法[1]。夾點(diǎn)分析[2]最早由LINNHOFF 等提出,以熱力學(xué)為基礎(chǔ),分析過程系統(tǒng)中的溫度分布,選擇合適的夾點(diǎn)溫度及組合匹配方式以實(shí)現(xiàn)能量的最大回收利用。

        數(shù)學(xué)規(guī)劃法可分為確定性方法和啟發(fā)式算法[3-5]。確定性方法如牛頓法[6-8]、最速下降法[9]、分支定界法[10]等,這些方法在求解具有非凸特性的混合整數(shù)非線性規(guī)劃問題時(shí),求解效率較低,不能滿足大規(guī)模系統(tǒng)計(jì)算需求。而啟發(fā)式算法如模擬退火算法(simulated annealing, SA)[11-12]、禁忌算法(Tabu search, TS)[13]、微分進(jìn)化算法(differential evolution, DE)[14]、遺傳算法(genetic algorithm, GA)[15]及粒子群算法(particle swarm optimization, PSO)[16]等,能夠利用已有的經(jīng)驗(yàn)規(guī)則,隨機(jī)搜索求解域內(nèi)的全局最優(yōu)解。啟發(fā)式算法可操作性強(qiáng),計(jì)算效率高,因此在換熱網(wǎng)絡(luò)最優(yōu)化問題中得到越來越多的應(yīng)用。但是在種群進(jìn)化后期,由于種群多樣性的缺失,啟發(fā)式算法容易出現(xiàn)早熟收斂的問題。同時(shí)大規(guī)模系統(tǒng)的求解難度和復(fù)雜程度不斷增加,影響優(yōu)化效率和優(yōu)化精度。

        為了提高啟發(fā)式算法應(yīng)用于大規(guī)模系統(tǒng)的優(yōu)化質(zhì)量,提升算法的全局優(yōu)化性能。肖媛等[17]總結(jié)了啟發(fā)式算法應(yīng)用于換熱網(wǎng)絡(luò)優(yōu)化的特點(diǎn),揭示了其早熟收斂現(xiàn)象的本質(zhì),并在此基礎(chǔ)上提出了一種新的強(qiáng)制進(jìn)化隨機(jī)游走算法(random walk algorithm with compulsive evolution, RWCE)。相比于其他啟發(fā)式算法,RWCE 具有程序簡(jiǎn)單、算法適應(yīng)性強(qiáng)的特點(diǎn),具備跳出局部最優(yōu)解的能力,從而提升了換熱網(wǎng)絡(luò)的優(yōu)化質(zhì)量。

        采用RWCE 優(yōu)化換熱網(wǎng)絡(luò),分析其優(yōu)化進(jìn)程。發(fā)現(xiàn)個(gè)體在每次進(jìn)化時(shí)都是以一定的幾率隨機(jī)參與游走,因此在優(yōu)化過程中不一定每股流體上都有換熱單元被選中參與進(jìn)化。對(duì)于沒有參與進(jìn)化的流股,該次迭代便沒有意義。如要保證每股流體都進(jìn)化一次需要增加大量的迭代次數(shù),使得換熱網(wǎng)絡(luò)的進(jìn)化效率低下。故本文提出以熱流體為基準(zhǔn)每股流體必進(jìn)化策略,即每隔一定周期強(qiáng)制使得每股流體上都有換熱單元被選中參與進(jìn)化,并每隔一定步數(shù)強(qiáng)制接受每股流體都游走的結(jié)構(gòu),使其跳出局部極值,提升算法的優(yōu)化效率和優(yōu)化質(zhì)量。

        2 換熱網(wǎng)絡(luò)優(yōu)化模型

        換熱網(wǎng)絡(luò)優(yōu)化問題描述為:以年度總費(fèi)用(total annual cost,TAC)最低為目標(biāo),通過在NH股熱流體與NC股冷流體之間設(shè)置換熱單元及冷、熱公用工程,使冷熱流體達(dá)到目標(biāo)溫度,實(shí)現(xiàn)能量回收。本文采用換熱網(wǎng)絡(luò)節(jié)點(diǎn)非結(jié)構(gòu)模型[18],以2 股熱流體和2 股冷流體為例,非結(jié)構(gòu)模型如圖 1 所示。其中,H1、H2 表示兩股熱流體,C1、C2 表示兩股冷流體,公用工程布置在冷熱流股末端,每股冷熱流體上分別預(yù)設(shè)NGDC、NGDH個(gè)產(chǎn)生換熱單元的位置。

        本文基于換熱網(wǎng)絡(luò)非結(jié)構(gòu)模型,以TAC 最小為目標(biāo),數(shù)學(xué)描述為:

        該模型包含NH股熱流體、NC股冷流體,熱流體上共有NKH個(gè)節(jié)點(diǎn),冷流體上共有NKC個(gè)節(jié)點(diǎn),NKH= NH× NGDH,NKC= NC× NGDC。TAC 如式(1)和(2)所示。

        圖1非結(jié)構(gòu)模型結(jié)構(gòu)圖Fig.1 Non-structural model structure

        式中,TAC 表示年度綜合費(fèi)用;FG表示固定投資費(fèi)用;CA表示換熱單元的面積費(fèi)用系數(shù);A 表示設(shè)備面積,m2;β 表示面積費(fèi)用指數(shù);Q 表示公用工程熱負(fù)荷,kW;C 表示公用工程費(fèi)用系數(shù);下角標(biāo)CU 表示冷公用工程,HU 表示熱公用工程;z, i, j 分別表示換熱單元數(shù)、熱流體數(shù)和冷流體數(shù),z∈Nh,c,i∈NH,j∈NC;Nh,c表示可能存在的所有換熱單元的個(gè)數(shù),N′h,c實(shí)際存在的換熱單元的個(gè)數(shù),C’表示換熱單元是否存在的0-1 變量。

        換熱單元的熱負(fù)荷及面積的熱平衡關(guān)系式如式(3)和(4)所示。

        3 改進(jìn)RWCE 優(yōu)化換熱網(wǎng)絡(luò)

        3.1 RWCE 原理

        RWCE 的基本思想[17]是:隨機(jī)產(chǎn)生初始換熱網(wǎng)絡(luò)結(jié)構(gòu),以目標(biāo)函數(shù)TAC 減小為強(qiáng)制進(jìn)化方向,隨機(jī)改變換熱單元的面積實(shí)現(xiàn)連續(xù)變量和整型變量的同步優(yōu)化;當(dāng)游走到局部最優(yōu)點(diǎn)時(shí),以一定幾率隨機(jī)接受差解,使其跳出局部極值,尋找全局最優(yōu)。主要操作步驟為:

        (1) 種群初始化。隨機(jī)產(chǎn)生初始種群M,其中包含n 個(gè)個(gè)體Zi,i = 1,…,n,每個(gè)個(gè)體對(duì)應(yīng)一個(gè)換熱網(wǎng)絡(luò)結(jié)構(gòu)。初始種群表示如式(7)~(9)所示。

        其中,η 為保留系數(shù),指可保留的最小熱負(fù)荷與最大游走步長(zhǎng)ΔL 的比值。

        b. 新個(gè)體產(chǎn)生階段。隨機(jī)選擇冷熱流股上的節(jié)點(diǎn),判斷該節(jié)點(diǎn)處是否存在換熱單元,若不存在,則在該位置以一定概率隨機(jī)生成一定熱負(fù)荷的換熱單元。如式(15)~(17)所示。

        其中,γ,γ1,γ2為0 到1 之間均勻分布的隨機(jī)數(shù);u 為熱流體上節(jié)點(diǎn),v 為冷流體上節(jié)點(diǎn),Qp,q為在u,v節(jié)點(diǎn)處新生成換熱單元的熱負(fù)荷。

        (3) 個(gè)體選擇階段。迭代后,若TAC 比前一次迭代值更低,則接受該次迭代的換熱網(wǎng)絡(luò)結(jié)構(gòu),并將該結(jié)構(gòu)對(duì)應(yīng)的熱負(fù)荷及節(jié)點(diǎn)編號(hào)作為下一次迭代的初始結(jié)構(gòu);若迭代后,TAC 不低于上一次迭代的值,則以一定概率δ 接受該次迭代的換熱網(wǎng)絡(luò)結(jié)構(gòu)作為下一次迭代的初始結(jié)構(gòu);否則不更新該個(gè)體。如式(18)所示。

        其中,r 為0 到1 之間均勻分布的隨機(jī)數(shù)。

        (4) 終止條件。記錄迭代過程中TAC 最低的換熱網(wǎng)絡(luò)結(jié)構(gòu),當(dāng)滿足終止條件時(shí),迭代終止。

        3.2 以熱流體為基準(zhǔn)每股流體必游走策略

        相較于其他優(yōu)化方法,RWCE 優(yōu)化換熱網(wǎng)絡(luò)能得到較低的TAC,但是需要較長(zhǎng)的計(jì)算時(shí)間,優(yōu)化效率仍有待提高。研究發(fā)現(xiàn)個(gè)體在每次進(jìn)化時(shí)都以一定幾率隨機(jī)選擇已存在的換熱單元參與進(jìn)化,且進(jìn)化幾率值較低,因此不能確定每次進(jìn)化時(shí)每股流體上都有換熱單元被選中參與進(jìn)化。以13 股熱流體7 股冷流體的算例為例,進(jìn)化過程中的節(jié)點(diǎn)選中幾率設(shè)為0.2,分析每次進(jìn)化時(shí)流體不被選中的幾率。

        由圖2 可知,進(jìn)化前期,換熱單元數(shù)目存在較大波動(dòng)時(shí),流股被選中幾率也有較大變化;進(jìn)化后期,換熱單元數(shù)基本固定在某一數(shù)值附近,此時(shí)每股流體不被選中的幾率基本固定在0.74 左右。即在進(jìn)化后期,13 股熱流體中平均每次迭代只能選中3 股左右參與進(jìn)化。對(duì)于當(dāng)前次進(jìn)化沒有被選中的流股,該次迭代就沒有意義;要保證每股流體都進(jìn)化一次需要消耗大量的迭代次數(shù),對(duì)換熱網(wǎng)絡(luò)的進(jìn)化產(chǎn)生不利影響。針對(duì)此特點(diǎn),采取以熱流體為基準(zhǔn)每股流體必游走策略,即采取措施使得每隔一定周期每股流體上都至少有一個(gè)換熱單元被選中參與進(jìn)化,從而均勻地增加每股流體的進(jìn)化幾率,以提升換熱網(wǎng)絡(luò)的優(yōu)化效率。

        主要策略是在存在的個(gè)體進(jìn)化處采取措施。在個(gè)體進(jìn)化之前,預(yù)先判斷每股流體上已經(jīng)存在的換熱單元數(shù),重新編號(hào),并與原來的編號(hào)一一對(duì)應(yīng)。每隔 itG1次迭代,采取一次每股流體上都有一個(gè)換熱單元被選中參與進(jìn)化的措施;其余迭代步數(shù)以原有的進(jìn)化方式隨機(jī)參與進(jìn)化。其中,itG1表示采取必游走策略的迭代周期。每隔一定周期采取一次每股流體必游走策略,既能在個(gè)體進(jìn)化到一定程度之后給種群一定的方向指導(dǎo),又能在其他迭代期間使種群保持原有的隨機(jī)進(jìn)化方向,保留其隨機(jī)性特點(diǎn)。

        以13 股熱流體7 股冷流體的算例[19](算例1)為例,流股參數(shù)如表1 所示。采用RWCE 優(yōu)化換熱網(wǎng)絡(luò)并采取每股流體必游走策略,當(dāng)參數(shù)設(shè)為:NGDH= 10,NGDC= 10,ΔL = 100,η = 0.2,δ = 0.01,Qmax= 100,itG1= 1 000 時(shí),TAC 變化曲線如圖 3 示。

        圖2 流股未被選中幾率Fig.2 The probability of the stream unselected in each iteration

        由圖 3 可知,采取必游走策略的RWCE 能比原始 RWCE 以更快的速度尋找到具有更低 TAC 的結(jié)構(gòu),有效提高了算法的優(yōu)化質(zhì)量。

        3.3 以熱流體為基準(zhǔn)每股流體必進(jìn)化策略

        采用加入每股流體必游走策略的RWCE 優(yōu)化換熱網(wǎng)絡(luò),分析發(fā)現(xiàn),只在存在的個(gè)體進(jìn)化處采取措施,不能保證每次進(jìn)化之后的個(gè)體一定會(huì)被接受。采取必游走策略之后,每股流體上都有換熱單元的熱負(fù)荷發(fā)生變化,故該次迭代后的費(fèi)用很大幾率上比前一次迭代的費(fèi)用高。此時(shí),進(jìn)化后較高費(fèi)用的結(jié)構(gòu)便不會(huì)被接受。因此該策略只能在一定程度上改變進(jìn)化方向,不能從根本上跳出局部極值。為充分發(fā)揮必游走策略的方向指導(dǎo)作用,進(jìn)一步在個(gè)體選擇位置采取措施,每隔 itG2次迭代強(qiáng)制接受一次每股流體必游走之后的結(jié)構(gòu),即每股流體必進(jìn)化策略,其中 itG2表示強(qiáng)制接受每股流體必游走之后的結(jié)構(gòu)的迭代周期。

        表1 算例1 參數(shù)數(shù)據(jù)Table 1 Flow parameters of case 1

        圖3 必游走TAC 變化曲線圖Fig.3 TAC curve for the participatory evolution strategy for streams

        表2 算例2 參數(shù)數(shù)據(jù)Table 2 Flow parameters of case 2

        每隔一定周期采取一次每股流體必進(jìn)化策略,使得換熱網(wǎng)絡(luò)的優(yōu)化既能夠保留原有進(jìn)化方法的隨機(jī)性特性,又能給該隨機(jī)進(jìn)化方法一定的方向指導(dǎo),更快速地找到更優(yōu)的結(jié)構(gòu)。以13 股熱流體7 股冷流體和8 股熱流體7 股冷流體的算例[26](算例2)為例,采用RWCE 優(yōu)化換熱網(wǎng)絡(luò)并采取每股流體必進(jìn)化策略,8 股熱流體7 股冷流體的算例流股參數(shù)如表2 所示。當(dāng)將13 股熱流體7 股冷流體的算例參數(shù)設(shè)為:NGDH= 10,NGDC= 10,ΔL=100,η = 0.2,δ = 0.01,Qmax= 100,itG1= 1 000,itG2= 5.5×106時(shí),采取必進(jìn)化策略的TAC 與原程序的TAC 變化曲線如圖4、圖5 示。

        圖4 必進(jìn)化實(shí)時(shí)TAC 變化曲線圖 Fig.4 TAC curve for compulsive participatory evolutionary strategy in the real time

        圖5 必進(jìn)化最優(yōu)TAC 變化曲線圖 Fig.5 Lowest TAC curve for compulsive participatory evolutionary strategy

        當(dāng)將8 股熱流體7 股冷流體的算例參數(shù)設(shè)為:NGDH= 10,NGDC= 10,ΔL = 110,η = 0.2,δ = 0.01,Qmax= 100,ITG1= 1 000,ITG2= 5.5×106時(shí),采取必進(jìn)化策略的TAC 與原程序的TAC 變化曲線如圖6、圖7 所示。

        圖6 必進(jìn)化實(shí)時(shí)TAC 變化曲線圖 Fig.6 TAC curve for the compulsive participatory evolutionary strategy in the real time

        圖7 必進(jìn)化最優(yōu)TAC 變化曲線圖 Fig.7 Lowest TAC curve for compulsive participatory evolutionary strategy

        由圖可知,采取每股流體必進(jìn)化策略的RWCE,能以更高的速率找到更優(yōu)解,顯著地提高了RWCE的求解速率及求解質(zhì)量。在每股流體必進(jìn)化策略中,每一次強(qiáng)制接受每股流體必游走的結(jié)構(gòu)之后,費(fèi)用會(huì)急劇增加,但每次的費(fèi)用增加,都會(huì)打破原有的進(jìn)化方向,跳出局部極值,從而增強(qiáng)RWCE 的局部搜索能力。

        3.4 改進(jìn)RWCE 步驟

        將結(jié)合上述策略的改進(jìn)RWCE 應(yīng)用于換熱網(wǎng)絡(luò)優(yōu)化,存在的個(gè)體進(jìn)化部分的改進(jìn)步驟如下:

        Step 1:判斷已經(jīng)存在換熱單元的位置,將已存在的換熱單元重新編號(hào),并與原來的節(jié)點(diǎn)編號(hào)一一對(duì)應(yīng);

        Step 2:對(duì)種群內(nèi)部的所有個(gè)體進(jìn)行隨機(jī)游走操作:每隔一定迭代次數(shù)采取一次每股流體上至少有一個(gè)換熱單元被選中參與游走的措施,其余迭代步數(shù)仍按照原來的進(jìn)化幾率進(jìn)行隨機(jī)游走操作;

        個(gè)體選擇部分的改進(jìn)步驟為:

        Step 1:每隔一定迭代次數(shù)后強(qiáng)制接受一次當(dāng)前迭代后的結(jié)構(gòu):經(jīng)過it 次迭代,無論TAC 是否比上一次迭代更低,都強(qiáng)制接受該換熱網(wǎng)絡(luò)結(jié)構(gòu)及每個(gè)節(jié)點(diǎn)的熱負(fù)荷,并將該結(jié)構(gòu)的熱負(fù)荷作為第(it+1)次迭代的初始熱負(fù)荷。如式(19)所示。

        Step 2:其余迭代次數(shù)仍以得到更低的TAC 為目標(biāo),進(jìn)行隨機(jī)游走操作。經(jīng)過it 次迭代,若TAC 比上一次迭代更低,則接受換熱網(wǎng)絡(luò)結(jié)構(gòu)及每個(gè)節(jié)點(diǎn)的熱負(fù)荷,并將該結(jié)構(gòu)的熱負(fù)荷作為第(it+1)次迭代的初始熱負(fù)荷,否則不更新該個(gè)體,保持第it 次迭代的換熱網(wǎng)絡(luò)結(jié)構(gòu)及熱負(fù)荷。如式(20)所示。

        Step 3:變異操作。若經(jīng)過it 次迭代,TAC 不低于上一次迭代的費(fèi)用,則以一定概率δ 接受該次迭代的換熱網(wǎng)絡(luò)結(jié)構(gòu)及熱負(fù)荷,作為第(it+1)次迭代的初始結(jié)構(gòu)及初始熱負(fù)荷。如式(21)所示。

        其中,r 為0 到1 之間均勻分布的隨機(jī)數(shù),δ 為當(dāng)(it+1)次迭代的TAC 大于it 次迭代時(shí),更新?lián)Q熱網(wǎng)絡(luò)結(jié)構(gòu),接受差解的概率。其算法流程圖如圖8,9,10 所示。

        圖8 存在的個(gè)體進(jìn)化Fig.8 Random walk evolution for each individualin the population

        圖9 個(gè)體選擇階段Fig.9 Individual selections

        圖10 必進(jìn)化策略流程圖Fig.10 Flow chart of the weighted evolutionary strategy

        4 算例驗(yàn)證

        4.1 算例1

        文獻(xiàn)中的優(yōu)化方法大多基于YEE 等[20-22]提出的分級(jí)超結(jié)構(gòu)模型進(jìn)行的。 ESCOBAR 等[23]采用等溫混合假設(shè)法,得到算例1 的TAC 為1 537 086 $·a-1。Pav?o 等[24]采用改進(jìn)的粒子群算法,得到TAC為 1 516 482 $·a-1;BAO 等[25]采用保護(hù)最優(yōu)解策略的 RWCE,得到 TAC 為 1 462 363 $·a-1。本文采用每股流體必游走策略的RWCE,得到的TAC 為1 432 904 $·a-1,費(fèi)用結(jié)構(gòu)圖如圖11 所示;采用每股流體必進(jìn)化策略的RWCE,得到的TAC 為1 411 131 $·a-1,費(fèi)用結(jié)構(gòu)圖如圖12 所示;算例1 的結(jié)果對(duì)比如表3 所示。

        圖11 采用每股流體必游走策略的RWCE 優(yōu)化結(jié)果Fig.11 RWCE optimization results using the participatory evolutionary strategy for streams

        圖12 采用每股流體必進(jìn)化策略的RWCE 優(yōu)化結(jié)果Fig.12 RWCE optimization results using the compulsive participatory evolutionary strategy

        4.2 算例2

        BJORK 等[26]采用數(shù)學(xué)規(guī)劃法改造換熱網(wǎng)絡(luò)得到算例 2 的 TAC為 1 530 063 $·a-1;PENG 等[27]采用改進(jìn)模擬退火算法得到的TAC 為1 527 240 $·a-1;LEANDRO 等[28]采用模擬退火算法和粒子群算法相結(jié)合的雙層優(yōu)化方法,得到的TAC 為1 525 394 $·a-1;本文采用每股流體必游走策略的RWCE,得到的TAC 為1 519 791 $·a-1,費(fèi)用結(jié)構(gòu)圖如圖13 所示;采用每股流體必進(jìn)化策略的RWCE,得到的TAC 為1 518 969 $·a-1,費(fèi)用結(jié)構(gòu)圖如圖14 所示;算例2 的結(jié)果對(duì)比如表4 所示。

        表3 算例1 優(yōu)化結(jié)果對(duì)比Table 3 Comparison on optimized results of case 1

        圖13 采用每股流體必游走策略的RWCE 優(yōu)化結(jié)果Fig.13 RWCE optimization results using the participatory evolutionary strategy for streams

        圖14 采用每股流體必進(jìn)化策略的RWCE 優(yōu)化結(jié)果Fig.14 RWCE optimization results using the compulsive participatory evolutionary strategy

        由算例 1 及算例 2的優(yōu)化結(jié)果及對(duì)比可知,在相同的目標(biāo)函數(shù)、約束條件和費(fèi)用系數(shù)等條件下,RWCE 算法能比其他優(yōu)化方法以更快的速度得到更優(yōu)的結(jié)果,尤其是改進(jìn)的RWCE 算法。采用以熱流體為基準(zhǔn)每股流體必進(jìn)化策略的RWCE算法,有效提高了換熱網(wǎng)絡(luò)的優(yōu)化效率及優(yōu)化質(zhì)量。

        表4 算例2 優(yōu)化結(jié)果對(duì)比Table 4 Comparison on optimized results of case 2

        5 結(jié) 論

        本文提出的以熱流體為基準(zhǔn)每股流體必進(jìn)化策略的RWCE 算法,能夠在保留非結(jié)構(gòu)模型一定自由性的基礎(chǔ)上,給非結(jié)構(gòu)模型一定的方向指導(dǎo)。既不破壞非結(jié)構(gòu)模型的自由性與靈活性,也能夠使非結(jié)構(gòu)模型的優(yōu)化目標(biāo)更趨向于全局最優(yōu)解,從而提高換熱網(wǎng)絡(luò)的優(yōu)化質(zhì)量。

        通過兩個(gè)不同規(guī)模的算例驗(yàn)證提出的策略,得到的結(jié)果分別為1 411 131 和1 518 969 $·a-1。結(jié)果表明,采用每股流體必進(jìn)化策略的換熱網(wǎng)絡(luò)優(yōu)化方法能夠有效提高非結(jié)構(gòu)模型的優(yōu)化質(zhì)量。

        但是每股流體必進(jìn)化策略仍存在一定的不足。在換熱網(wǎng)絡(luò)優(yōu)化后期,出現(xiàn)不存在公用工程的流股,此時(shí)采用每股流體必進(jìn)化策略的效果便不如預(yù)期。因此每股流體必進(jìn)化策略的優(yōu)化方法仍需繼續(xù)深入研究,以改進(jìn)該策略在換熱網(wǎng)絡(luò)優(yōu)化后期的不足,進(jìn)一步提升該策略在換熱網(wǎng)絡(luò)非結(jié)構(gòu)模型中的優(yōu)化質(zhì)量。

        猜你喜歡
        算例流體個(gè)體
        流體壓強(qiáng)知多少
        山雨欲來風(fēng)滿樓之流體壓強(qiáng)與流速
        關(guān)注個(gè)體防護(hù)裝備
        等效流體體積模量直接反演的流體識(shí)別方法
        基于振蕩能量的低頻振蕩分析與振蕩源定位(二)振蕩源定位方法與算例
        互補(bǔ)問題算例分析
        個(gè)體反思機(jī)制的缺失與救贖
        How Cats See the World
        基于CYMDIST的配電網(wǎng)運(yùn)行優(yōu)化技術(shù)及算例分析
        燃煤PM10湍流聚并GDE方程算法及算例分析
        亚洲精品中文字幕一区二区| 亚洲性无码av在线| 亚洲中文字幕永久网站| 成人av综合资源在线| 色婷婷亚洲一区二区三区| 国产va免费精品观看| 欧美激情国产一区在线不卡| 亚洲av成人波多野一区二区| 激情内射人妻1区2区3区| 久久精品娱乐亚洲领先| 中文字幕人成乱码熟女| 国产成人精品一区二区三区免费| 玩弄放荡人妻一区二区三区| 国产极品喷水视频| 少妇人妻无奈的跪趴翘起| 国内熟女啪啪自拍| 香蕉人妻av久久久久天天| 日韩人妻无码精品二专区| 免费在线视频亚洲色图| 人妻中文字幕无码系列| 亚洲国产成人久久综合一区77| 99精品免费视频| 国产大片在线观看91| 日本中文字幕一区二区有码在线| 射死你天天日| 色综合999| 亚洲av第一区国产精品| 99久久婷婷国产综合精品青草免费| 亚洲视频在线看| 久久久国产视频久久久| 无遮挡很爽很污很黄的女同| 日韩欧群交p片内射中文| 色综合久久久久久久久五月| 人妻丰满熟妇av一区二区| av天堂中文亚洲官网| 97久人人做人人妻人人玩精品| 国产午夜精品一区二区三区| 中文字幕不卡在线播放| 亚洲一区二区三区av天堂| 亚洲精品国偷拍自产在线| 狠狠躁夜夜躁人人爽天天天天97|