| 查看: 1058 | 回复: 0 | ||
snipers2004金虫 (初入文坛)
|
[求助]
udf读取txt文件的数据错误求助
|
|
气体的物性参数太多了，需用udf导入相应的txt文件，如“氩气输运系数.txt”。但是我编写了一个程序，却无法将数据导入。 #include "udf.h" #define pi 3.14159265358979 #define mu0 (4*pi*1.0e-7) //真空磁导率 #define DataMax 298 //共导入298个数据 #define Tmin 300 //允许的最小温度 #define Tmax 30000 //允许的最大温度 real Press[298],TT[298],vis[298],thc[298],elc[298],den[298],sph[298]; DEFINE_EXECUTE_ON_LOADING(on_loading, libname) { int i,j; FILE *fPRO; fPRO=fopen("氩气输运系数.txt","r"  ;// Open and read the Gas Transport and Propertiesfor(i=0;i<298;i++) { fscanf(fPRO,"%le %le %le %le %le\n",&vis,&thc,&elc,&den,&sph); //fscanf(fPRO,"%e %e %e %e %e",&vis,&thc,&elc,&den,&sph); //fscanf(fPRO,"%e %e %e %e %e",&(vis),&(thc),&(elc),&(den),&(sph)); //fscanf(fPRO,"%le %le %le %le %le",&(vis),&(thc),&(elc),&(den),&(sph)); //Message("%le %le %le %le %le\n",vis,thc,elc,den,sph); } fclose(fPRO); } real sigma(cell_t c,Thread *t) { int ip; real ele_con; if (C_T(c,t)<Tmin) C_T(c,t)=Tmin; else if (C_T(c,t)>Tmax) C_T(c,t)=Tmax; ip=floor((C_T(c,t)-200)/100); ele_con=elc[ip]+(C_T(c,t)-100*ip-200)/100*(elc[ip+1]-elc[ip]); C_UDMI(c,t,0)=ele_con; return ele_con; } DEFINE_PROPERTY(density,c,t)/*密度*/ { int i1; real rho; real temp1=C_T(c,t); i1=floor((temp1-200)/100); rho=den[i1]+(temp1-100*i1-200)/100*(den[i1+1]-den[i1]); C_UDMI(c,t,1)=rho; return rho; } DEFINE_PROPERTY(viscosity,c,t)/*粘滞率*/ { int i2; real mu; real temp2=C_T(c,t); i2=floor((temp2-200)/100); mu=vis[i2]+(temp2-100*i2-200)/100*(vis[i2+1]-vis[i2]); C_UDMI(c,t,5)=mu; return mu; } DEFINE_PROPERTY(thercond,c,t)/*热导率*/ { int i3; real kp; real temp3=C_T(c,t); i3=floor((temp3-200)/100); kp=thc[i3]+(temp3-100*i3-200)/100*(thc[i3+1]-thc[i3]); C_UDMI(c,t,8)=kp; return kp; } 请问各位，上面的代码有什么错误？ |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有74人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
-
招聘 金属平磨液，抛光液研发工程师
已经有3人回复
25