C 課程設計報告精餾塔理論塔板數的計算

c++課程設計報告

精餾塔理論塔板數的計算

化學化工學院

高階工程人才實驗班（化學）

趙靜冰學號1507110114

目錄（1）設計報告3

（2）摘要4

（3）程式設計思路5

（4）流程框圖6

（5）c 語言原始程式15

（6）程式執行結果17

（7）後記18

（8）參考文獻19

二、設計報告

1、撰寫要求

（1）c語言課程設計任務書；

（2）目錄；

（3）摘要；

（4）程式設計思路和流程框圖；

（5）c語言原始程式；

（6）程式執行結果；

（7）後記（程式除錯過程出現問題的討論）。

2、印裝和上交要求

（1）設計報告要求採用a4紙雙面列印，訂書針裝訂。

（2）將源程式按統一檔名100a-bc形式發給輔導老師

課程設計報告

摘要：此次c++程式課程設計，是在給定的設計條件下，利用梯級**法計算苯—甲苯精餾塔理論塔板數。第乙個程式使用了直接計算法解方程組，計算出塔頂d和塔底w，使用最小二乘擬合法和高斯消元法求得相平衡方程，同時使用牛頓迭代的法求得q線方程與相平衡方程的交點以及q線方程與精餾段操作線方程等的交點，迴圈使用牛頓迭代法最後求得精餾塔理論塔板的數目和進料板的位置；在設計的過程中所需要的c++語言方法有：

迴圈、選擇、陣列、繪圖等的綜合使用，最後完成整個設計過程。

程式設計思路：由於沒有學習過化工原理方面的知識，所以初開始，對於陌生的專有名詞感覺這次的設計很不好做，不知從何下手。後來老師給出了公式等推導原理，並在老師講解後我才大概懂得題目的要求，參考了之前學過的課本和一些範例，再詢問其他合作同學，我終於有了設計思路。

我將整個程式設計過程分為兩步，計算部分和作圖部分，而計算又分為兩個部分，程式一過程使用了直接計算法解方程組，計算出塔頂d和塔底w，程式二過程先利用高斯消元法求得相平衡方程，再利用牛頓迭代法求得精餾塔理論塔板數的過程。

1、高斯函式：利用高斯主元素消去法可以求解平衡線方程中的係數及常數a0,a1,a2,a3。

2、矩陣板塊：我們在求解平衡線方程中的係數及常數a0,a1,a2,a3時，有乙個問題，那就是矩陣的行數和列數不相同，不能完全照搬書上所說的高斯主元素，而是需要利用資料的擬合。我們可以用矩陣的轉置與矩陣的乘法（參考實際程式），達到此目的。

3、牛頓迭代法：通過翻閱資料，運用牛頓迭代法是求方程根，在所，先定義一下牛頓迭代法，然後再在所需要的地方呼叫一下，

4、迴圈板塊：陣列可以解決資料龐大的問題，迴圈的過程中，輸出其角標即可。採用的是do….while語句。

c++設計流程圖：

c++原始程式：

程式一 : //求解方程組

#include<>

#define f 46.61 //f為處理量

#define xf 0.45 //xf為苯-甲苯混合液中苯的含量

#define xd 0.996 //xd為塔頂產品摩爾分率

#define xw 0.0118 //xw為塔底產品摩爾分率

int main()

程式二 ://求精餾段與提鎦段方程，並求塔板數

#include <>

double diedai(double a,double b,double c,double d迭代定義

while(fabs(xn1-xn)>=1e-6);

return xn1;

}void main()

,y=,

a,b,c,da,b,c,d為函式的係數

for(i=0;i<25;i++)

}for(i=0;i<3;i求解

for(k=1;k<4-i;k++)

}for(i=3;i>=0;i求三次函式的係數a，b,c,d

cout<<"三次函式的係數分別為："< cout<<"a="< cout<<"b="< cout<<"c="< cout<<"d="< cout< i=0;j=0;

int q;

double m=21.4048,g=46.61,k=25.

2052,xf,xw,xd,ye,xe,rmin,t,a0=a[0],a1=a[1],a2=a[2],a3=a[3],r,xn,xn1,f,f1,g,g1,xq,yq,xm,xm1,x[50],p,y1,y2,x0,x1,x2,e;

double q=1; //進料熱狀態q的值為1

cout< cout< cout<<"請輸入xf: ";

cin>>xf;

cout< cout<<"請輸入xd: ";

cin>>xd;

cout< cout<<"請輸入xw: ";

cin>>xw;

cout< xe=xf;

ye=a0+a1*xe+a2*pow(xe,2)+a3*pow(xe,3); //平衡方程

rmin=(xd-ye)/(ye-xe);

cout<<"輸入最小回流比前的係數t(t=1.1~2.0): ";

cin>>t;

cout< r=t*rmin;

xq=xf;

yq=r*xq/(r+1)+xd/(r+1);

cout<<"xe="< cout< cout<<"最小回流比為："< cout< cout<<"回流比為："< cout< cout<<"精餾操作線方程為：

"< cout<<"提餾操作線方程為："< cout< cout<<"xq="< cout< cout< cout< x[1]=diedai(a0-xd,a1,a2,a3); //牛頓迭代法求解

i=1;

do{i++;

C 課程設計報告精餾塔理論塔板數的計算

篩板精餾塔課程設計

化工原理篩板精餾塔課程設計

課程設計指導書板式精餾塔

C 課程設計報告精餾塔理論塔板數的計算

篩板精餾塔課程設計

化工原理篩板精餾塔課程設計

課程設計指導書 板式精餾塔

相關推薦

課程設計指導書板式精餾塔