| ²é¿´: 560 | »Ø¸´: 2 | ||
zhaoyubosdͳæ (СÓÐÃûÆø)
|
[ÇóÖú]
ÇóÖúC++±à³Ì£¬²»Ê¤¸Ð¼¤
|
|
ÔËÐÐC++³ÌÐò£¬ÏÔʾÔËÐгɹ¦£¬ÀïÃæÓÐÕâÑùÒ»¾ä»°£¬ÇëÎÊÊÇʲôÒâ˼£¿The thread 0x2cb0 has exited with code 0 (0x0).ÈçºÎ½â¾ö¡£»¹ÓУ¬ºÚÉ«ÔËÐнçÃæÉÏΪʲôûÓÐÔËÐнá¹û£¬Çë¸÷λ´óÉñÖ¸½Ì¡£¸½ÉϳÌÐò´úÂë¡£ #include <iostream> #include <fstream> #include <iomanip> #include <cmath> #include <string> using namespace std; void fibermodel (double X[], int n0, double U[], int n1, double V[], int n2, double P[], int n3, double R[],int n4, double r, double l, double mu, double dp, double den, double converge, int n, double &dX); int main () { //Hollow fiber double X[301],t[301], U[301], V[301], Pressure[301], R[301]; double r, l, T,Rm, mu, dp, den, dX,converge,cp; double ResultP[51][51],ResultU[51][51],ResultV[51][51],ResultR[51][51]; int n,m,IncrementT,IncrementX; string OutputFileName; //Information about the software cout << "HOLLOW FIBER MEMBRANE SIMULATION SOFTWARE Version 1.0" << "\n" << "developed by Dr. Lianfa Song, College of Engineering, Texas Tech University" << "\n" << "APRIL 2015" << "\n" << "\n" << "\n"; // Open input file and read parameter values ifstream input_file("numeric.txt" ; input_file >> r >> l >> T >> Rm >> mu >> dp >> den >> converge >> n >> m >> IncrementX >> IncrementT >> cp >> OutputFileName; cout << "The Name of Output File = " << OutputFileName << "\n" << "\n"; // Start calculation // Call fibermodel int j,i; t[0]=0.0; for (j=1;j<=m;j++) {t[j]=t[j-1]+T/m;} for (j=0;j<=m;j++) { if (j=0) { for(i=0;i<=n;i++) { R=Rm;} } fibermodel (X, 301, U, 301, V, 301, Pressure, 301, R, 301, r, l, mu, dp, den, converge, n, dX); if (j%6==0) { for(i=0;i<=n;i=i+6) { ResultP[j]=Pressure;ResultU[j]=U;ResultV[j]=V;ResultR[j]=R;} } for(i=0;i<=n;i++) { R=R+cp*V*T/m;} } //Output results // Open output file ofstream output_file(OutputFileName); output_file << "CONDITIONS" << "\n" ; output_file << fixed; output_file.precision(0); output_file << setw(13) << "Space-Step=" << setw(13) << n << "\n"; output_file << scientific; output_file.precision(4); output_file << setw(13) << "InnerRadius=" << setw(13) << r << "\n"; output_file << setw(13) << "FiberLength=" << setw(13) << l << "\n"; output_file << setw(13) << "MemResisten=" << setw(13) << Rm << "\n"; output_file << setw(13) << "WaterViscos=" << setw(13) << mu << "\n"; output_file << setw(13) << " Pressure=" << setw(13) << dp << "\n"; output_file << setw(13) << "ConvergeCre=" << setw(13) << converge << "\n" << "\n"; output_file << fixed << "RESULTS1" << "\n"; //PµÄ·Ö²¼Çé¿ö output_file << fixed << "Pressure" << "\n"; output_file << scientific; output_file.precision(3); output_file<< "X"; for (int j=0; j<=m; j=j+IncrementT) {output_file<< setw(13) << "t="<< t[j] << setw(13);} output_file<< endl; for (int i=0; i<=n; i=i+IncrementX) { output_file<< X ; for (int j=0; j<=n; j=j+IncrementT) { output_file<< setw(13) << ResultP[j] << setw(13); } output_file<< endl; } output_file<< endl; //uµÄ·Ö²¼Çé¿ö output_file << fixed << "u" << "\n"; output_file << scientific; output_file.precision(3); output_file<< "X"; for (int j=0; j<=m; j=j+IncrementT) {output_file<< setw(13) << "t="<< t[j] << setw(13);} output_file<< endl; for (int i=0; i<=n; i=i+IncrementX) { output_file<< X ; for (int j=0; j<=n; j=j+IncrementT) { output_file<< setw(13) << ResultU[j] << setw(13); } output_file<< endl; } output_file<< endl; //vµÄ·Ö²¼Çé¿ö output_file << fixed << "v" << "\n"; output_file << scientific; output_file.precision(3); output_file<< "X"; for (int j=0; j<=m; j=j+IncrementT) {output_file<< setw(13) << "t="<< t[j] << setw(13);} output_file<< endl; for (int i=0; i<=n; i=i+IncrementX) { output_file<< X ; for (int j=0; j<=n; j=j+IncrementT) { output_file<< setw(13) << ResultV[j] << setw(13); } output_file<< endl; } output_file<< endl; //RµÄ·Ö²¼Çé¿ö output_file << fixed << "R" << "\n"; output_file << scientific; output_file.precision(3); output_file<< "X"; for (int j=0; j<=m; j=j+IncrementT) {output_file<< setw(13) << "t="<< t[j] << setw(13);} output_file<< endl; for (int i=0; i<=n; i=i+IncrementX) { output_file<< X ; for (int j=0; j<=n; j=j+IncrementT) { output_file<< setw(13) << ResultR[j] << setw(13); } output_file<< endl; } output_file<< endl; system("PAUSE" ;return 0; } // fibermodel defination void fibermodel (double X[], int n0, double U[], int n1, double V[], int n2, double P[], int n3,double R[], int n4, double r, double l, double mu, double dp, double den, double converge, int n, double &dX) { //Preparation double UP, LP, K, pi = 3.1415926; K=dX; dX = l/double(n); X[0]=0.0; for (int i=1; i<=n; i++) {X=X[i-1]+dX;} //Start Calculation int iter=0; UP=dp; LP=0.0; // loop: iter = iter +1; P[n] = 0.5*(UP+LP); U[n] =0.0; for (int i=n; i>=1; i--) { V = P/R; U[i-1] =U+2.0*dX*V/r; P[i-1]=P+8.0*mu*U[i-1]*dX/(r*r)+K*3.0*den*U[i-1]*V*dX/r; if (P[i-1] > dp) {UP=P[n]; goto loop;} } if (abs((P[0]-dp)/dp) > converge) {LP=P[n]; goto loop;} V[0] = dp/R[0]; return; } |
» ²ÂÄãϲ»¶
Ò»Ö¾Ô¸211 0703»¯Ñ§ 346·ÖÇóµ÷¼Á
ÒѾÓÐ26È˻ظ´
ÉúÎïѧµ÷¼Á
ÒѾÓÐ13È˻ظ´
293Çóµ÷¼Á
ÒѾÓÐ15È˻ظ´
335Çóµ÷¼Á
ÒѾÓÐ17È˻ظ´
È˹¤ÖÇÄÜ320µ÷¼Á08¹¤À໹Óлú»áÂð
ÒѾÓÐ14È˻ظ´
271Çóµ÷¼Á
ÒѾÓÐ27È˻ظ´
¿¼ÑÐÇóµ÷¼Á
ÒѾÓÐ7È˻ظ´
085408¹âµçÐÅÏ¢¹¤³Ìר˶355Ò»Ö¾Ô¸³¤´º¹â»úËùµ÷¼Á
ÒѾÓÐ3È˻ظ´
071000ÉúÎïѧ£¬Ò»Ö¾Ô¸ÉîÛÚ´óѧ296·Ö£¬Çóµ÷¼Á
ÒѾÓÐ15È˻ظ´
302Çóµ÷¼Á
ÒѾÓÐ7È˻ظ´
zzyµÎºÅ
½ð³æ (СÓÐÃûÆø)
- Ó¦Öú: 18 (СѧÉú)
- ½ð±Ò: 5449.9
- É¢½ð: 367
- ºì»¨: 1
- Ìû×Ó: 282
- ÔÚÏß: 895.6Сʱ
- ³æºÅ: 1540568
- ×¢²á: 2011-12-17
- ÐÔ±ð: GG
- רҵ: ÀíÂۺͼÆË㻯ѧ

2Â¥2015-10-21 22:58:08
leetcoder
гæ (³õÈëÎÄ̳)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 119.6
- Ìû×Ó: 28
- ÔÚÏß: 6.1Сʱ
- ³æºÅ: 4132819
- ×¢²á: 2015-10-10
- רҵ: ¼ÆËã»úÓ¦Óü¼Êõ
3Â¥2015-10-22 09:14:39













;
»Ø¸´´ËÂ¥
5