RK 方程展开型体积根的求解 - 针对具体物质的 ExcelVBA 特殊函数求解法-1753286009191.pptVIP

单元2气体的热力学性质计算*针对具体物质的RK方程特殊函数求解法*RK方程展开型体积根的求解

——针对具体物质的ExcelVBA特殊函数求解法*针对具体物质的RK方程特殊函数求解法*步骤1：编写求解一元3次非线性方程的通用函数Newton3ExcelVBA通用函数代码如下：PublicFunctionNewton3(V0,ESP,A,B,C,D)AsDoubleDimy0,Y1,X,x1AsDoublex1=V0DoX=x1y0=A*X^3+B*X^2+C*X+DY1=3*A*X^2+2*B*X+Cx1=X-y0/Y1LoopUntilAbs((x1-X)/X)ESPNewton3=x1EndFunction*针对具体物质的RK方程特殊函数求解法*步骤2：调用Newton3函数，编写针对R-K方程体积根求解的通用函数，代码如下：PublicFunctionRKEqV_pTpcTcESP(p,T,pc,Tc,ESP)AsDoubleDimV0,R,a,b,A1,A2,A3,A4AsDoubleR=8.314V0=R*T/pa=R^2*Tc^2.5/9/(2^(1/3)-1)/pcb=(2^(1/3)-1)*R*Tc/3/pcA1=1A2=-R*T/pA3=(a/T^0.5-b*R*T-p*b^2)/pA4=-a*b/p/T^0.5RKEqV_pTpcTcESP=Newton3(V0,ESP,A1,A2,A3,A4)EndFunction步骤3：编写查询临界参数的自定义函数Tc_C3H6()，Pc_C3H6()，代码如下：PublicFunctionTc_C3H6Tc_C3H6=365EndFunctionPublicFunctionPc_C3H6Pc_C3H6=4.62EndFunction*针对具体物质的RK方程特殊函数求解法*步骤4：编写针对具体物质的R-K方程特殊型函数,代码如下:PublicFunctionRKV_C3H6(n,Ct,MPa)AsDoubleDimTK,P,TC,PcAsDoubleTK=Ct+273.15P=MPa*1000000#Tc=Tc_C3H6()‘调用丙烯临界温度自定义函数Pc=1000000#*Pc_C3H6()‘调用丙烯临界压力自定义函数RKV_C3H6=n*RKEqV_pTpcTcESP(P,TK,Pc,TC,0.0001)EndFunction步骤5：针对原始条件n=1mol、t=125℃、p=10MPa,直接调用自定义特殊函数RKV_C3H6(n,t,p)，在单元格B26中输入：=RKV_C3H6(1,125,10)，按“Enter”键,即得V值,如下图所示:注意:采用Excel单变量求解体积时可能得不到正确解!为什么?*注意:采用Excel单变量求解体积时可能得不到正确解!为什么?*注意:采用Excel单变量求解体积时可能得不到正确解!为什么?*单元2气体的热力学性质计算注意:采用Excel单变量求解体积时可能得不到正确解!为什么?*注意:采用Excel单变量求解体积时可能得不到正确解!为什么?*注意:采用Excel单变量求解体积时可能得不到正确解!为什么?*