|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
程序如下：
2 n: u h$ D2 ^6 A& {function out=mult_syg7(z,x)2 y4 M6 Z- u: A
Pi=3.14159;5 l8 ~7 b9 Y/ @$ D2 p
r0= 1.24*10^-3; % Outer diameter of the membrane, [m]6 s9 L" N/ z, p# X2 g
ri= 9.4*10^-4;% Inner diameter of the membrane, [m]
0 x0 |& Z/ m0 @- m* s& hr=(r0-ri)/log(r0/ri);% equivalent diameter ,[m]
& ^) k0 P! g* Q- FL=0.05;% membrane length,[m]
: Z# C2 R6 m& _3 a6 ^, M% D2 l# E% N& ^Ac=Pi*r*L;% membrane area [m^2]
; J) n9 X9 ?' L: LR=8.314;4 M0 Y% X5 y- Y& R2 R; J
T0=1000;2 m* }4 h: j/ [/ s. H: B
T=T0+273;
8 O2 b. C Y/ [" Y/ f j- FP=101.325;%kpa7 U6 I: ? |! L4 `2 F! `
F=2.678*10^-2*2;
% k9 v* E8 |! Z/ h7 A. k# km=3;& `- A1 @: Z, ?1 l
n=1;$ F# b. n3 s% B
S=(1+m+n+2*x(1))/P;: i v: ? c2 c4 Z% M4 h
PCH4=(1-x(1))/S;
+ V& ^( L K' y# vPH2O=(m-x(1)-x(2))/S;
' H& s7 c3 |( ?) Mout= ...2 M6 ?+ z/ m5 T4 O3 \7 ]0 d) A
[Ac/F*1.5*10^-3*exp(-117*10^3/(R*T))*PCH4.^0.97.*PH2O.^-0.08. @/ O. ]& J* ^
(3*x(1)+x(2)).*x(2)./(1-x(2))./(m-x(1)-x(2))-exp(4400/T-4.063)];
3 w! M' R f! F( W: w# |+ d& N* u! l% f2 h6 h
调用如下:5 d- e# Z8 `$ y! G7 J
>> x0 = [0; 0];
i+ ]+ B* @' T1 B5 vtspan = [0:0.1:1];
. o; I0 w( c0 v! w PM = [1 0 ; 0 0];! A. ^3 u- Q+ O! J: b) N
options = odeset('Mass',M,'RelTol',1e-4,'AbsTol',[1e-6 1e-10 ]);
- ]% c/ T/ }1 X4 W2 o. k5 K[t,x] = ode15s(@mult_syg7,tspan,x0,options)% w5 E0 L1 h, b" M
7 ]8 d- a: i! }( m8 M7 j$ ?% x报错如下：
8 L8 ~" o# Q& p* _ n4 O9 C错误使用 daeic12 (line 76), x; G4 Z+ @9 t* U* A
此 DAE 的索引大于 1。
2 E: E7 @( A# y7 f( y% d! V9 C$ I# z, S) Y" @# J& [
出错 ode15s (line 310)& I6 ] C9 X. h: q n9 j
[y,yp,f0,dfdy,nFE,nPD,Jfac] = daeic12(odeFcn,odeArgs,t,ICtype,Mt,y,yp0,f0,...+ E" d8 ^. ?4 h" O
|
|
/1 
发表于 2021-3-19 14:48
收藏
淘帖
支持！
反对！