24СʱÈÈÃŰæ¿éÅÅÐаñ    

Znn3bq.jpeg
²é¿´: 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;
        }
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

zzyµÎºÅ

½ð³æ (СÓÐÃûÆø)

´úÂëÌ«»ìÂÒÁË£¬¶¼²»Äܳɹ¦±àÒë¹ýÈ¥¡£Êý×éÓë±äÁ¿Ö±½ÓÏàµÈ£¬²»Í¬´óСµÄÊý×éÖ±½ÓÓÃ=²Ù×÷£¬»¹ÓÐ10¼¸¸ö²ÎÊýµÄº¯Êý¡£»¹ÊÇÈ·ÈÏÒ»ÏÂÕâÊDz»ÊÇÕýÈ·µÄ´úÂë°É¡£
ÌìµÀ³êÇÚ
2Â¥2015-10-21 22:58:08
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

leetcoder

гæ (³õÈëÎÄ̳)

ÏÈ˵ÄãÕâ¸ö³ÌÐòÓÃÀ´¸ÉÂï¡ú_¡ú

·¢×ÔСľ³æAndroid¿Í»§¶Ë
3Â¥2015-10-22 09:14:39
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ zhaoyubosd µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] »¯¹¤Ñ§Ë¶ 285Çóµ÷¼Á +28 Wisjxn 2026-04-07 28/1400 2026-04-13 10:56 by ѧzh
[¿¼ÑÐ] 297Çóµ÷¼Á +18 ORCHID1 2026-04-10 19/950 2026-04-13 09:21 by lhj2009
[¿¼ÑÐ] 327Çóµ÷¼Á +14 Xxjc1107. 2026-04-06 14/700 2026-04-13 08:31 by Delta2012
[¿¼ÑÐ] 297¹¤¿Æ£¬Çóµ÷¼Á? +11 ºÓÄÏũҵ´óѧ-ÄÜ 2026-04-12 11/550 2026-04-13 08:29 by wp06
[¿¼ÑÐ] 277Çóµ÷¼Á +23 Äß½¨Éè 2026-04-06 23/1150 2026-04-13 00:40 by Íõͯ×ÓÁú
[¿¼ÑÐ] Ò»Ö¾Ô¸Õã´óÉúÎï325·ÖÇóµ÷¼Á +9 zysheng 2026-04-12 9/450 2026-04-12 22:31 by yuyin1233
[¿¼ÑÐ] µ÷¼ÁÇóÊÕÁô +29 ¹ûÈ»ÓÐÎÒ 2026-04-10 30/1500 2026-04-12 22:14 by zxcwyt
[¿¼ÑÐ] 086000µ÷¼Á +6 Ê®Æßsa 2026-04-07 6/300 2026-04-12 11:05 by ´óÁ¦Ë®ÊÖÁ¦´óÎÞÇ
[¿¼ÑÐ] Çóµ÷¼Á£¬262»úеר˶ +8 àÅyyl 2026-04-08 8/400 2026-04-12 02:31 by Çï¶¹²ËÑ¿
[¿¼ÑÐ] 305Çóµ÷¼Á +6 77Qi 2026-04-07 6/300 2026-04-12 02:30 by Çï¶¹²ËÑ¿
[ÕÒ¹¤×÷] ɽ¶«¸ßУ½Ìʦ¿¼ºË³¬¼¶ÎÞµ×Ïߣ¬Ô±¹¤¹ý²»ÏÂÈ¥À² +4 qut2026 2026-04-09 9/450 2026-04-12 00:54 by qut2026
[¿¼ÑÐ] 280Çóµ÷¼Á +7 ÙâÙâÒ¹Ò¹ 2026-04-09 10/500 2026-04-12 00:33 by À¶ÔÆË¼Óê
[¿¼ÑÐ] Çóµ÷¼Á +6 СÄô°®Ñ§Ï° 2026-04-11 9/450 2026-04-11 21:20 by À¶ÔÆË¼Óê
[¿¼ÑÐ] 085400 328·Ö Çóµ÷¼Á +10 ιÄãÒ»¸ö´ó³È×Ó 2026-04-09 14/700 2026-04-11 19:53 by lqspecial
[¿¼ÑÐ] 283Çóµ÷¼Á +22 ÄǸöàà×Ó 2026-04-09 22/1100 2026-04-11 10:41 by ÄæË®³Ë·ç
[¿¼ÑÐ] 087100³õÊÔ311Çóµ÷¼Á +4 ÈÎÑÅÇÙ 2026-04-09 4/200 2026-04-11 10:33 by zhq0425
[¿¼ÑÐ] »¹Óл¯¹¤¶þÂÖµ÷¼ÁµÄѧУÂð 5+14 »¯¹¤ÈË999 2026-04-09 48/2400 2026-04-11 10:27 by 89436494
[¿¼ÑÐ] »úеר368 ÓÐÈ¥´¦Âð +4 ÖÖ´óÊ÷ 2026-04-10 4/200 2026-04-10 15:31 by jiajinhpu
[¿¼ÑÐ] 297Çóµ÷¼Á +27 GENJIOW 2026-04-07 30/1500 2026-04-09 23:20 by wolf97
[¿¼ÑÐ] ÇóÖú +3 ¿¨¿¨¶«88 2026-04-06 4/200 2026-04-06 15:28 by going home
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û