myGully.com - Ausgabe Datei erstellen, Frage

Ja ich das Program macht immer noch das selbe, es speichert immer noch jede Zeile. Ich hab schon nach 1sek. 9MB :D und das Program läuft par Minuten bis es fertig ist.

Ja ich hab Zeahlerdateizeile mit double definiert. Aber ich was es genau ist nicht, nur das es = 0 ist.

hier mal das ganze Program:
PHP-Code:

		
			
#include <stdio.h>

#include <stdlib.h>

#include <math.h>



double M_PI = 3.1415926535897932384626433832795028841971693993751058209749445923078164062862089986280348253421170679;

    FILE *ausgabe;

 

  

  double Zaehlerdateizeile = 1;

  double I;                                                                                // Grundlage: I = U/Rt

  double U;

  double Rt;

 

  // BERECHNUNG DER SPANNUNG IN ABHAENGIGKEIT VON DER ZEIT:  U(t) = 230 * sqrt(2) * sin (t * 2 * Math.PI * f)

  

  double t = 0;                                                                            // Zeit in Sekunden

  int FREQUENZ = 50;                                                                    // Frequenz = 50 Hz, genormte Netzfrequenz                                   

 

  // BERECHNUNG DES WIDERSTANDES IN ABHAENGIGKEIT VON DER TEMPERATUR: Rt = Rel * (1 + alpha * (T1-T0)), wobei (T1-T0) = delta_T ist



  double ALPHA = 0.0048;                                                                // linearer Temperaturkoeffizient 4,8 * 10 ^-3 K^-1 bei 20°C

  double delta_T;                                                                        // T1-T0 (T1 => Drahttemperatur, T0 => Raumtemperatur)

  double Tkelv;                                                                            // Temperatur des Drahtes in °Kelvin   

  double Tkelv_g;                                                                        // Temperatur des Glases in °Kelvin

  double raumTemp = 293.15;                                                                // Raumtemperatur in Kelvin                                          

  double Rel;                                                                            // Elektrisches Widerstand            

  double roSpez = 0.0528;                                                                // Spezifischer Widerstand von Wolfram in Ohm^2

  double laenge = 0.62;                                                                    // Laenge in m 0.062

  double A_d;                                                                            // Querschnitt in m^2 des Drahtes A_d = r ^2 * Math.PI

  double durchmesser = 0.000015;                                                        // Durchmesser in m                                    

  



  // BERECHNUNG DER TEMPERATURDIFFERENZ (delta_T)

  

  // C = DeltaQ/DeltaT  <=> DeltaT = DeltaQ/C

  // C = c * m

  

  double C_w;                                                                            // Waermekapazitaet in                            J/K (Wolfram)

  double C_g;                                                                            // Wearmekapazitaet in                            J/K (Glass)

  int c_w = 134;                                                                        // spezifische Waermekapazitaet von Wolfram        in J/(kg*K)  

  int c_g = 700;                                                                        // spezifische Waermekapazitaet von Glas        in J/(kg*K) 

  double m_d;                                                                            // Masse des Drahtes in kg                        m = V_d * Rho_W

  double m_g;                                                                            // Masse Glaskolben in kg

  double V_d;                                                                            // Volumen des Drahtes in m^3

  double V_g = 0.12;                                                                    // Volumen des Glaskolbens in m^3 ergibbt sich aus Subration 2er Kugeln

  int Rho_W = 19270;                                                                    // Dichte von Wolfram in kg/m^3

  int Rho_G = 2500;                                                                        // Dichte von Glass in kg/m^3

    

  // BERECHNUNG DER THERMISCHEN ENERGIE



  double delta_Q;                                                                        // thermische Energie in J DeltaQ = Pw * delta_t

  double Pw;                                                                            // Waermeleistung in Watt

  

  // BERECHNEN DER WAERMELEISTUNG

  

  double Pel;                                                                            // elektrische Leistung in Watt  Pel = Pw + Pstr <=> Pw = Pel - Pstr

  double Pstr;                                                                            // Strahlungsleistung in Watt

  

  

  // STRAHLUNGSLEISTUNG BERECHNEN

  

  double BK = 0.0000000567;                                                                // Stefan-Boltzmann-Konstante J/K Pstr = T^4 * Bk * Ob * Epsilon(T)

  double Ob_d;                                                                            // Oberflaeche m^2

  double Ob_g;

  double eps_w;                                                                            // Emissionsgrad Wolfram

  double eps_g;                                                                            // Emissionsgrad Glass

  

  // ZUR RECHNUNG BENÖTIGTEN WERTE



  double INTERVALL = 0.000001;            //0.000001                                                // entspricht Delta t in Sekunden

  double Terwaerm_draht;                                                                // Erwaermung des Drahtes

  double Terwaerm_glas;                                                                    // Erwearmung des Glases

  double ausdehn_l;                                                                        // Laenge nach Ausdehnung

  double ausdehn_l_glas;                                                                // Leange nach Ausdehnung Glas

  double ausdehn_d;                                                                        // Durchmesser nach Ausdehnung

  double ausdehn_d_glas;                                                                // Durchmesser Ausdehnung glas

  double ausdehnungskoeff_w = 0.0000045;                                                // Ausdehnungskoeffizient Wolfram

  double ausdehnungskoeff_g    = 0.0000076;                                                // Ausdehnungskoeffizient Glas

  double dauerSimulation = 1.0;                                                            // Simulationsdauer in Sekunden       

  double Kelvin = 273.15;                                   

 

 

 // BERECHNUNG DER KONSTANEN WERTE DIE FÜR DEN ERSTEN DURCHLAUF BENÖTIGT WERDEN



 void erzeugeKonstante() {

      

      A_d = pow((durchmesser/2),2) * M_PI;

      Rel = 58.5;                                                                        // Rel = roSpez * laenge / A

      Tkelv = raumTemp;

      Tkelv_g = raumTemp;

      Ob_d = 2 * (durchmesser/2) * M_PI * laenge;

      m_d = A_d * laenge * Rho_W;                                                        // Berechnung der Masse des Drahtes

      m_g = Rho_G * V_g;                                                                // Berechnung des Masse des Glaskolebsn (Volumen ergibt sich auch Subtraktion zweier Kugeln, 55mm und 54mm, daraus resultiert die Wandstaerke)

      C_w = c_w * m_d;                                                                    // Berechnung der Waermekapazitaet (Wolfram)

      C_g = c_g * m_g;                                                                    // Berechnung der Waermekapazitaet (Glas)



}

 // BERECHNUNG DER WERTE DIE SICH AUFGRUND DER AUSDEHNUNG AENDERN  



void berechneVariable() {



        delta_T = Tkelv - Tkelv_g;

        ausdehn_l = laenge * ( 1 + ausdehnungskoeff_w * delta_T );

        ausdehn_d = durchmesser * ( 1 + ausdehnungskoeff_w * delta_T);

        Ob_d = 2 * (ausdehn_d/2) * M_PI * ausdehn_l;

        Ob_g = M_PI * pow((ausdehn_d_glas),2);

}



// BERECHNUNG EINES INTERVALLS



void durchlaufEinesIntervalls() {

        

        berechneVariable();                                                                // BERECHNUNG DER VERAENDERTEN VARIABLEN

         

        U = 230 * sqrt(2.0) * sin(t * 2 * M_PI * FREQUENZ);

        Rt = Rel * ( 1 + ALPHA * delta_T);

        I = U / Rt;

        Pel = (U*U)/Rt;

        eps_w = 0.39;                                                                    // Emissionsgrad Wolfram

        eps_g = 0.94;                                                                    // Emissionsgrad  Glas

        Pstr = BK * ((Ob_d * ((pow(Tkelv,4))-(pow(Tkelv_g,4))) * eps_w) - (Ob_g * ((pow(Tkelv_g,4))-(pow(raumTemp,4))) * eps_g)) ; // Waermeabgabe von Draht zu Glaskolen und von Glaskolben zur Umgebungsluft bzw. Raumptemperatur

        Pw = Pel - Pstr;

        Terwaerm_draht = (Pw * INTERVALL) / C_w;

        Terwaerm_glas = (Pw * INTERVALL) / C_g;

        Tkelv = Tkelv + Terwaerm_draht;

        Tkelv_g = Tkelv_g + Terwaerm_glas;

        t = t + INTERVALL; 

}



// AUSGABE DER WERTE IN EINE DATEI

void gebeAus() {                                                                        //Datei oeffnen

if(Zaehlerdateizeile%100)

{

    Zaehlerdateizeile++;

    return;

}

    ausgabe = fopen("Ergebnis.txt","a");

if(ausgabe != NULL) {



     fprintf(ausgabe,                                                                    // Die Ausgabe

     "%2.6lf s %15.6lf A %15.6lf V %15.6lf Ohm %10.6lf W (Pw) %15.6lf W (Pel) %15.6lf W (Pstr) %20.10lf Kelvin (Temp. d. Drahtes) %20.10lf Kelvin (Temp. d. Glaskolbens) \n",\

     t,I,U,Rt,Pw,Pel,Pstr,Tkelv,Tkelv_g);  

     Zaehlerdateizeile +=1;

    }

    fclose(ausgabe);                                                                    // Datei schließen

}



// DIE SIMULATION DES EINSCHALTENS

void simuliere() {

    erzeugeKonstante();                                                                    // Schleifendurchlaeufe

    

double i = 0;



for ( i; i <= dauerSimulation; i += INTERVALL) {

      printf("%2.6f / %2.6f\r", i , dauerSimulation);

      durchlaufEinesIntervalls();

      gebeAus();   

    }

}



int main()

{

    printf("SIMULATION DES EINSCHALTENS EINER GLUELAMPE.\n\nBITTE WARTEN BIS DIE SIMULATION BEENDET IST, DANKE\n");

    

    simuliere();                                                                         // Beginn der Simulation

        

    printf("\n\n\nDIE SIMULATION IST BEENDET. SIE KOENNEN NUN DIE ERGEBNIS \nDATEI OEFFNEN. DIESE FINDEN SIE IM MICROSOFT VISUAL C++ PROJEKT ORDNER\n");

    system("PAUSE");    

    return 0;