資源描述:
《計(jì)算機(jī)在化學(xué)化工中的應(yīng)用資料》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在應(yīng)用文檔-天天文庫(kù)。
1、YibinUniversity《計(jì)算機(jī)在化學(xué)化工中的應(yīng)用》期末總結(jié)報(bào)告題目用計(jì)算機(jī)解決《化工原理》中的計(jì)算問(wèn)題專業(yè)應(yīng)用化學(xué)學(xué)生姓名XXXXXX學(xué)號(hào)XXXXXXXXX年級(jí)2014級(jí)指導(dǎo)教師XXXXXXXXX 化學(xué)與化工學(xué)院用計(jì)算機(jī)解決《化工原理》中的計(jì)算問(wèn)題1引言隨著現(xiàn)代科學(xué)技術(shù)的發(fā)展和計(jì)算機(jī)的廣泛應(yīng)用,各學(xué)科對(duì)計(jì)算機(jī)的依賴程度越來(lái)越高,化學(xué)化工領(lǐng)域也不例外。例如,實(shí)驗(yàn)數(shù)據(jù)的處理及擬合、模型參數(shù)的確定、非線性方程組的求解、化工過(guò)程模擬,都離不開(kāi)計(jì)算機(jī)的幫助。2問(wèn)題的提出現(xiàn)在,對(duì)于現(xiàn)階段的我們來(lái)說(shuō),需要用到計(jì)
2、算機(jī)來(lái)解決的問(wèn)題,主要是用計(jì)算機(jī)解決一些復(fù)雜的計(jì)算求解。本次總結(jié)報(bào)告,我將著重介紹用計(jì)算機(jī)解決下列問(wèn)題的過(guò)程及結(jié)果:?jiǎn)栴}一、本學(xué)期化工原理課程設(shè)計(jì)中“年處理40000噸乙醇-水溶液連續(xù)精餾塔的設(shè)計(jì)”中,精餾塔塔板數(shù)的計(jì)算和進(jìn)料板位置的確定。已知:進(jìn)料組成xF(摩爾分?jǐn)?shù),下同)為0.174,相對(duì)揮發(fā)度α為3.04,泡點(diǎn)進(jìn)料q=1,最小回流比Rmin=1.8,操作回流比R=1.8Rmin=3.24,要求塔頂產(chǎn)品組成xD為0.7789,塔底產(chǎn)品組成xW為0.00196,產(chǎn)品氣化率e=D/F=0.2214。問(wèn)題二
3、、化工原理(下冊(cè))第八章吸收,課后習(xí)題21題(4)小問(wèn)中,吸收因子法:NOG=AA-1lnA-1Ay1-mx2y2-mx2+1A已知:NOG、y1、y2、m、x2,求解吸收因子A。3用計(jì)算機(jī)解決問(wèn)題3.1問(wèn)題一的解決針對(duì)問(wèn)題一:“精餾塔塔板數(shù)的計(jì)算和進(jìn)料板位置的確定?!笨紤]用圖解法和逐板計(jì)算兩種方法來(lái)解決。當(dāng)然,這兩種方法均可以不用計(jì)算機(jī),但不用計(jì)算機(jī)就會(huì)有圖解法不夠精確、誤差大和逐板計(jì)算運(yùn)算量過(guò)大的問(wèn)題,而用計(jì)算機(jī)就不會(huì)有這些問(wèn)題。3.1.1圖解法考慮到精確度的問(wèn)題,解決這個(gè)問(wèn)題不考慮用Office作圖
4、,而考慮用Origin7.5漢化版。圖解法的具體操作步驟如下:①根據(jù)給定體系的平衡數(shù)據(jù)繪制平衡曲線y-x圖;也可有相對(duì)揮發(fā)度α作圖。②根據(jù)精餾段操作線方程,過(guò)點(diǎn)(xD,xD)作精餾段操作曲線。③根據(jù)q線方程作q線;當(dāng)q>0時(shí),q線過(guò)點(diǎn)(xF,xF)及點(diǎn)(xF/q,0)。當(dāng)q≤0時(shí),q線過(guò)點(diǎn)(xF,xF)及點(diǎn)(0,xF/(1-q))。當(dāng)q=1時(shí),q線過(guò)點(diǎn)(xF,xF)及點(diǎn)(xF,0)。④根據(jù)提餾段操作線方程,過(guò)點(diǎn)(xW,xW)作提餾段操作曲線。⑤在圖中,從(xD,xD)出發(fā),向左作水平線交相平衡線于1點(diǎn),過(guò)
5、1點(diǎn)向下作垂線交精餾段操作線于1'點(diǎn),完成一個(gè)梯級(jí);以此類推,當(dāng)水平線與平衡線交點(diǎn)橫坐標(biāo)小于xF時(shí),垂線改交于提餾段操作曲線,繼續(xù)作提餾段操作線與平衡線之間的梯級(jí),直到水平線與平衡線交點(diǎn)橫坐標(biāo)小于xW為止。根據(jù)梯級(jí)圖確定塔板數(shù)和進(jìn)料板位置。由問(wèn)題一知:平衡線方程、精餾段操作線方程、提餾段操作線方程、q線方程分別為:y=3.04x1+2.04x、y=0.764x+0.1837、y=1.8294x-0.001626、x=0.174。最終繪制梯級(jí)圖如下:[鍵入文檔的引述或關(guān)注點(diǎn)的摘要。您可將文本框放置在文檔中的
6、任何位置??墒褂谩拔谋究蚬ぞ摺边x項(xiàng)卡更改重要引述文本框的格式。]放大由圖可知:塔板數(shù)為13;進(jìn)料板位置為第4塊板。3.1.2主板計(jì)算法主板計(jì)算法步驟:①令y1=xD,并代入平衡線方程,解出x1。②將x1代入精餾段操作線方程,解出y2,再將y2代入平衡線方程,解出x2。以此類推,當(dāng)xm≤xq(q=1時(shí),xq=xF)時(shí),該板為進(jìn)料板(即為第m塊板)。③將xm代入提餾段操作線方程,解出ym+1。再將ym+1代入平衡線方程,解出xm+1。以此類推,當(dāng)xm+n≤xw時(shí),結(jié)束計(jì)算。理論塔板數(shù)為:NT=m+n。用VB編
7、寫程序代碼如下:PrivateSubCommand1_Click()Dima,x,y,xd,xw,xf,r,k,yy,xx,ea=Val(InputBox("相對(duì)揮發(fā)度a"))r=Val(InputBox("回流比r"))xd=Val(InputBox("塔頂產(chǎn)品的摩爾分?jǐn)?shù)xd"))xw=Val(InputBox("塔底產(chǎn)品的摩爾分?jǐn)?shù)xw"))xf=Val(InputBox("進(jìn)料的摩爾分?jǐn)?shù)xf"))e=Val(InputBox("氣化率e"))y=xdk=1x=y/(a-(a-1)*y)100yy=((
8、r*x)/(r+1))+(xd/(r+1))xx=yy/(a-(a-1)*yy)k=k+1Ifxx>xfThenx=xxGoTo100ElseEndIfPrint"xx=";xx;"進(jìn)料板位置k=";kn=0x=xx200y=(((r+(1/e))*x)/(r+1))-((((1/e)-1)/(r+1))*xw)xx=y/(a-(a-1)*y)n=n+1Ifxx>xwThenx=xxGoTo200ElseEndIfq=k+nP