-
大小: 371KB文件類型: .rar金幣: 2下載: 0 次發(fā)布日期: 2021-06-08
- 語言: C/C++
- 標(biāo)簽: SPH無網(wǎng)格??
資源簡介
SPH無網(wǎng)格程序,可以調(diào)試使用的工程文件,C++和Fortran語音混合編程!

代碼片段和文件信息
/**************************************************************************************************/
void?export_vart(int?mint?ntdouble?t)
{
int?inreglength;
char?filename[10]filename0[10];
FILE?*fp2;
double?c0u0p0d0e0;
itoa(mfilename010);//把整數(shù)轉(zhuǎn)換為字符串。
length=strlen(filename0);//求字符串長度。
if(length==1)?strcpy(filename“0000“);//將一個字符串復(fù)制到另一字符串中。
if(length==2)?strcpy(filename“000“);
if(length==3)?strcpy(filename“00“);
if(length==4)?strcpy(filename“0“);
if(length==5)?strcpy(filename““);
strcat(filenamefilename0);//將一個字符串附加到另一個字符串上。
strcat(filename“.plt“);
fp2=fopen(filename“w“);//打開流!
fprintf(fp2“%s\n““VARIABLES?=?X?U?P?D?E?H?S?A?U0?P0?D0?E0“);
fprintf(fp2“%s\n““TEXTCS=frameC=REDS=LOCALX=60.0Y=2.0“);
fprintf(fp2“%s\n““HU=POINTLS=1AN=LEFTBX=FILLEDBXM=20LT=0.4“);
fprintf(fp2“%s\n““BXO=WHITEBXF=WHITEF=TIMES-BOLDH=20A=0MFC=\“\““);
fprintf(fp2“%s%12.5le%s\n““T=\“Time?=?“t/1.0e-6“?m s\““);
double?u1=0.0;//Input?Velocity1!
double?u2=0.0;//Input?Velocity2!
double?c1=1290.0;
double?c2=6519.97863934307;
double?p1=1.0e6;//Input?Pressure1!
double?p2=3.5e10;//Input?Pressure2!
double?p3=2.570826e10;//2.461595124e10;//Input?Pressure3!
double?d1=15750.0;//15.75;//Input?Density1!
double?d2=2470.0;//Input?Density2!
double?e1=15.8730158730159;
double?e2=7.08502024291498e6;
double?k1=5.0k2=3.0;
double?a2=2.32261477582487*1.0e9/pow(1000.0k2);//Input?a2!
double?d3=d1*(p1*(k1-1.0)+p3*(k1+1.0)+2.0*c1*c1*d1)/(p3*(k1-1.0)+p1*(k1+1.0)+2.0*c1*c1*d1);
double?d4=pow(p3/a21.0/k2);
double?u3=-sqrt((p3-p1)*(1.0/d1-1.0/d3));
double?e3=(p3-c1*c1*(d3-d1))/(k1-1.0)/d3;
double?e4=p3/d4/(k2-1.0);
double?D=d3*u3/(d3-d1);
double?XJ=u3*t;
double?XD=D*t;
double?XRH=c2*t;
double?XRT=(c2+(k2+1.0)/2.0*u3-(k2-1.0)/2.0*u2)*t;
for(nreg=1;nreg<=mreg;nreg++)
{
if(openreg[nreg]==0)?continue;
fprintf(fp2“%s%s%s%3d%s\n““ZONE?T=“namereg[nreg]“?I=“iip[nreg]“?F=POINT“);
for(i=1;i<=iip[nreg];i++)
{
if(rp[nm[i][nreg]][nt].x<=XD?&&?nreg==1)
{
u0=u1;
p0=p1;
d0=d1;
e0=e1;
}
if(rp[nm[i][nreg]][nt].x>XD?&&?nreg==1)
{
u0=u3;
p0=p3;
d0=d3;
e0=e3;
}
if(rp[nm[i][nreg]][nt].x<=XRT?&&?nreg==2)
{
u0=u3;
p0=p3;
d0=d4;
e0=e4;
}
if(rp[nm[i][nreg]][nt].x>XRT?&&?rp[nm[i][nreg]][nt].x {
c0=(k2-1.0)/(k2+1.0)*(rp[nm[i][nreg]][nt].x/t-u2+2.0/(k2-1.0)*c2);
u0=(k2-1.0)/(k2+1.0)*(2.0/(k2-1.0)*rp[nm[i][nreg]][nt].x/t+u2-2.0/(k2-1.0)*c2);
p0=p2*pow(c0/c22.0*k2/(k2-1.0));
d0=d2*pow(p0/p21.0/k2);
e0=p0/d0/(k2-1.0);
}
if(rp[nm[i][nreg]][nt].x>=XRH?&&?nreg==2)
{
u0=u2;
p0=p2;
d0=d2;
e0=e2;
}
fprintf(fp2“%15.5le%15.5le%15.5le%15.5le%15.5le%15.5le%5d%5d%15.5le%15.5le%15.5le%15.5le\n“
????rp[nm[i][nreg]][nt].xvp[nm[i][nreg]][nt].x
????????pp[nm[i][nreg]][nt]??dp[nm[i][nreg]][nt]??
ep[nm[i][nreg]][nt]??hp[nm[i][nreg]][nt]??
sp[nm[i][nreg]]??????ap[nm[i][nreg]]??????u0p0d0e0);
}
}
fclose(fp2);//關(guān)閉流!
}
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????2454??2008-12-03?15:31??SPH_YJW\artificial.for
?????文件????????592??2008-12-03?15:31??SPH_YJW\density.for
?????文件????????341??2008-12-03?15:31??SPH_YJW\eos.for
?????文件???????1239??2008-12-03?14:36??SPH_YJW\exact.for
?????文件???????1184??2009-03-05?19:21??SPH_YJW\hupdate.for
?????文件???????5890??2008-12-07?10:31??SPH_YJW\intergral.for
?????文件???????1514??2008-06-15?09:43??SPH_YJW\internal.for
?????文件???????1264??2008-02-28?17:29??SPH_YJW\kernel.for
?????文件???????3051??2009-10-14?15:44??SPH_YJW\M4_exact.cpp
?????文件???????3011??2009-10-14?15:44??SPH_YJW\M5_exact.cpp
?????文件???????2745??2008-12-07?08:59??SPH_YJW\model.for
?????文件???????1192??2009-10-15?14:10??SPH_YJW\neighbor.for
?????文件???????2141??2008-12-03?15:59??SPH_YJW\output.for
?????文件????????777??2009-10-14?08:46??SPH_YJW\param.inc
?????文件???????2454??2008-12-03?15:31??SPH_YJW\R4\artificial.for
?????文件???????3432??2009-10-15?13:11??SPH_YJW\R4\Debug\BuildLog.htm
?????文件???????2048??2009-10-15?13:11??SPH_YJW\R4\Debug\SPH.exe.em
?????文件????????592??2008-12-03?15:31??SPH_YJW\R4\density.for
?????文件????????489??2009-10-14?19:27??SPH_YJW\R4\eos.for
?????文件???????2184??2009-10-15?14:03??SPH_YJW\R4\exact.for
?????文件???????1184??2009-03-05?19:21??SPH_YJW\R4\hupdate.for
?????文件???????4315??2009-10-15?14:12??SPH_YJW\R4\intergral.for
?????文件???????1514??2008-06-15?09:43??SPH_YJW\R4\internal.for
?????文件???????1264??2008-02-28?17:29??SPH_YJW\R4\kernel.for
?????文件???????3051??2009-10-14?15:44??SPH_YJW\R4\M4_exact.cpp
?????文件???????3011??2009-10-14?15:44??SPH_YJW\R4\M5_exact.cpp
?????文件???????2511??2009-10-15?13:09??SPH_YJW\R4\model.for
?????文件???????1155??2009-10-15?14:07??SPH_YJW\R4\neighbor.for
?????文件???????5040??2009-10-14?22:31??SPH_YJW\R4\output.for
?????文件????????777??2009-10-14?08:46??SPH_YJW\R4\param.inc
............此處省略100個文件信息
評論
共有 條評論