已知某二元理想物系的A B两组分纯物质的蒸汽压数据如附表所示 现将25℃含A的质量百分数为30％的混
已知某二元理想物系的A、B两组分纯物质的蒸汽压数据如附表所示，现将25℃含A的质量百分数为30％的混合液，以每小时20000kg的流量加入常压操作的精馏塔中，要求塔底产品的流量为15000kg/h。塔顶全凝器冷凝物每10kg就有7kg回流入塔，实际回流比为最小回流比的1.8倍。（1）用M.T.法求理论板数及理论进料位置。若全塔总板效率为45％，求该精馏塔的实际板数及实际进料板位置；（2）设各实际板以汽相组成表示及以液相组成表示的单板效率均为0.7，求塔顶第二实际塔板下降的液相浓度和塔底第二块实际塔板上升汽相的浓度。已知A、B的摩尔质量分别是78kg/kmol和92kg/kmol，在平均操作温度范围内A、B组分的比热容为1.84kJ/（kg·℃)，组分A、B的汽化潜热分别为390kJ/kg和360kJ/kg。 附表 80.1 85 90 95 100 105 110.6 pAS kPa 101.3 116.9 135.5 155.7 179.2 204.2 240.0 pAS kPa 40
回答
共1条回复-
1．计算原理(化学原理和计算方法) 某溶液由两种物质组成，c A 为物质A的质量浓度/g·L -1 ，c B 为物质B的质量浓度/g· -1 ，μ为物质的黏度/cP。设此溶液黏度的数学模型为 μ=a 0 +a 1 c A +a 2 c B 将表中所列实验数据代入多元线性回归子程序(其中调用高斯消去法子程序解正规方程组)，通过计算确定方程中的各个系数a 0 、a 1 、a 2 值和复相关系数。 2．程序框图 3．源程序 programmain DOUBLE PRECISION XV,S,A,X,B0,QT,U,QE,SE,RF dilnension XV(11),S(11) COMMON/A×/A(10,11),X(50,10) OPEN(5,FILE='2006TE1.TXT',STATUS='OLD') OPEN(6,FILE='2006TE2.TXT',STATUS='UNKNOWN') READ(5,*) M,N M1=M+1 READ(5,*) ((X(I,J),J=1,M1),I=1,N) CALL MLR(M,N,XV,S,BO,QT,U,QE,SE,R,F) WRITE(6,20) M,N 20 FORMAT(/4X,'INPUT DATA:','Xn1,'1X,'Xn2,…',1X,'Xnm,Yn(m=1,'& ,i2,',n=1,',i2,')'/1X,66(1H-)/6X,'X1',8X,'X2',8X, & 'Y',7X,'Ycalc',7X,'DIFF'/1X,66 (1H-)) DO 40 I=1,N YCALC=BO DO 30 J=1,M 30 YCALC=YCALC+A(J,M1) * X(I,J) DIFF=YCALC-X(I,M1) 40 WRITE(6,50) (X(I,J),J=1,M1),YCALC,DIFF 50 FORMAT(3X,10G10.4) WRITE(6,60) (XV(I),I=1,M1) 60 FORMAT(/1X,'Xv',10G10.4) WRITE(6,70) (S(I),I=1,M1) 70 FORMAT(/1X,'Sv',10G10.4) WRITE(6,80) 80 FORMAT (/5X,'Regression equahon is:',/1x,60(1H-)) WRITE(6,90)B0,(A(I,M1),I=1,M) 90 FORMAT(1X,'μ=',G10.4,'+',G10.4,' * CA','+', & G10.4,' * CB') WRITE(6,100) QT,U,QE,SE,U/QT,R,F 100 FORMAT(/5X,'Qt=',G18.6,/5X,'U=',G18.6,/5X,'Qe=',G18.6,/5X, & G18 6) END SUBROUTINE MLR(M,N,XV,S,BO,QT,U,QE,SE,R,F) DOUBLE PRECISION XV,S,A,X,B0,QT,U,QE,SE,R F DIMENSION XV(M+1),S(M+1),ALI0(50) COMMON /AX/A(10,11),X(50,10) M1=M+1 DO 20 I=1,M1 XV(I)=0.D0 DO 10 K=1,N 10 XV(I)=XV(I)+X(K,I) 20 XV(I)=XV(I)/DBLE(N) DO 40 I=1,M1 S(I)=0 DO 30 K=1,N 30 S(I)=S(I)+(X(K,I)-XV(I)) ** 2 40 S(I)=DSQRT(S(I)) DO 50 I=1,M ALI0(I)=0 DO 50 K=1,N 50 ALI0(I)=ALI0(I)+(X(K,I)-XV(I)) * (X(K,M1)-XV(M1)) DO 60 I=1,M1 DO 60 K=1,N 60 X(K,I)=(X(K,I)-XV(I))/S(I) DO 70 I=1,M A(I,I)=1.0 DO 70 J=I+1,M1 A(I,J)=0 DO 70 K=1,N 70 A(I,J)=A(I,J)+X(K,I) * X(K,J) DO 80 I=1,M-1 DO 80 J=I+1,M 80 A(J,I)=A(I,J) WRITE(6,11) M 11 FORMAT(4X,'Normal equahon A(m,m+1)',2x,'m=',i2/1x 60(1H-)] DO 90 I=1,M 90 WRITE(6,33) (A(I,J),J=1,M1) CALL GS(A,M,M1,1.D-06) WRITE(6,22) WRITE(6,33) (A(I,M1),I=1,M) B0=0.D0 U=0.D0 DO 100 I=1,M A(I,Ml)=A(I,M1) * S(M1)/S(I) B0=B0+A(I,M1) * XV(I) 100 U=U+A(I,M1)*ALI0(I) B0=XV(M1)-B0 QT=S(M1) * S(M1) QE=QT-U SE=DSQRT(QE/(N-M-1)) R=DSQRT(1.D0-QE/QT) F=(N-M-1) * U/QE/M 22 FORMAT(/4X,'Roots of mormal equahon','A(i,m+l),i=1,m'& /lx,60(1H-)) 33 FORMAT(1X,10G12.4) RETURN END SUBROUTINE GS(A,N,M,EPS) DOUBLE PRECISION EPS,A DIMENSION A(10,11) DO 50 K=1,N P=0.0 DO 10 I=K,N IF(P GE ABS(A(I,K)))GOTO 10 P=ABS(A(I,K)) L=I 10 CONTINUE IF(P.LT.EPS)STOP'NO SOLUTION' IF(L.EQ.K)GOTO 30 DO 20 J=K,M TMP=A(L,J) A(L,J)=A(K,J) A(K,J)=TMP 20 CONTINUE 30 P=10/A(K,K) DO 40 J=K+1,M A(K,J)=A(K,J)*P DO 40 I=K+1,N A(I,J)=A(I,J)-A(I,K) * A(K,J) 40 CONTINUE 50 CONTINUE DO 60 K=2,N I=M-K DO 60 J=I+1,N A(I,M)=A(I,M)-A(I,J) * A(J,M) 60 CONTINUE RETURN END 4．运行结果 Regression equation is: μ=-27.43 + .2327 *CA+ .4095 * CB Qt= 2428.38 U= 1355.86 Qe= 1072.52 Se= 9.45390 N= .558341 R= .747222 F= .758515