Beiträge: 272
Registriert: 02.10.2008

erstellt am: 27. Jul. 2010 21:34      <-- editieren / zitieren -->      Unities abgeben: 

entschuldige, habe deine PM erst sehr spät gelesen.

Es freut mich, dass es jetzt bei dir funktioniert. Über die Variable staune ich auch, dass deine Version das fordert. Kann aber auch an den Solver liegen. Die Variable "true" war ja nur geschätzt und es scheint ja auch die richtige Richtung gewesen zu sein. 

Beiträge: 13
Registriert: 30.04.2011

erstellt am: 25. Jul. 2011 18:45      <-- editieren / zitieren -->      Unities abgeben: 

..und die hier erwähnte Ausgabe von Kräften hat auch nicht funktioniert. Gibt es eine bekannte und auch funktionierende Möglichkeit das Drehmoment zu erhalten?

Grüße
Christian

[Diese Nachricht wurde von ChristianE36 am 25. Jul. 2011 editiert.]

Beiträge: 272
Registriert: 02.10.2008

erstellt am: 25. Jul. 2011 20:03      <-- editieren / zitieren -->      Unities abgeben: 

Zitat:

---

.und die hier erwähnte Ausgabe von Kräften hat auch nicht funktioniert.

---

Die Funktion habe ich schon lange nicht mehr verwendet, aber sollte trotzdem funktionieren und dir die Momente und Kräfte für jeden Zeitschritt ausgeben. Was genau funktioniert nicht?

Zitat:

---

--calculateTorque: ließ sich nicht kompilieren.

---

Zitat:

---

-- computeTorque:Fehlermeldung; "no patch named rotor"

---

Beiträge: 13
Registriert: 30.04.2011

erstellt am: 25. Jul. 2011 21:59      <-- editieren / zitieren -->      Unities abgeben: 

1.  http://www.cfd-online.com/OpenFOAM_Discus /messages/1/__Client_C__Documents_and_Settings_c563640_Desktop_calculateTorque_tar-10038.unk
Das ist der direkte Link zum Download. das ganze wird wie üblich mit wmake kompiliert. Das habe ich heute einige stunden erfolglos probiert.
Die entsprechenden Pfadangaben stimmen nicht mehr. (und ich hab versagt beim "korrigieren". bin zu doof zum programmieren  )
Wenn du das schaffst, erkläre mir bitte wie du es gemacht hast.
Ich hab das schon bei einer Solverkompilierung nicht verstanden.

2. "computeTorque"
Wurde mal für 1.4.1 geschrieben. Das ließ sich sogar anstandslos kompilieren. aber beim benutzen hat es dann gehapert. Der Link: http://openfoamwiki.net/index.php/ComputeTorque
Erst wird der Code gezeigt, darunter dann der Link zum Download (1).

3. Das was in diesem Thread zuvor diskutiert wurde hatte ich probiert.
Dann werde ich es mir nochmal anschauen, was ich genau in die controlDict reinschreiben muss. Ich habe es so verstanden, dass es auf den "libraries" für das Post-Processing basiert.(habe ich UserGuide gelesen) Aus einer bereits erfolgreichen Berechnung habe ich die controlDict wie folgt verändert:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.5                                   |
|   \\  /    A nd           | Web:      http://www.OpenFOAM.org                |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
libs ("libOpenFoamTurbo.so")
libs ("libNewCylidricallInletBCs.so")
application     MRFSimpleFoam;

startFrom       startTime;

startTime       0;

stopAt          endTime;

endTime         5000;

deltaT          1;

writeControl    timeStep;

writeInterval   1000;

purgeWrite      0;

writeFormat     ascii;

writePrecision  8;

writeCompression compressed;

timeFormat      general;

timePrecision   6;

runTimeModifiable yes;

functions
(

forces
{
type forces;
functionObjectLibs ("libforces.so"); //Lib to load
patches (IMPELLER); // change to your patch name
rhoName rhoInf;
rhoInf 1000; //Reference density for fluid
CofR (2 0 0); //Origin for moment calculations
        outputControl   timeStep;
        outputInterval  1;
}
forceCoeffs
{
// rhoInf - reference density
// CofR - Centre of rotation
// dragDir - Direction of drag coefficient
// liftDir - Direction of lift coefficient
// pitchAxis - Pitching moment axis
// magUinf - free stream velocity magnitude
// lRef - reference length
// Aref - reference area
type forceCoeffs;
functionObjectLibs ("libforces.so");
patches (IMPELLER);
rhoName rhoInf;
rhoInf 1000;
CofR (2 0 0);
liftDir (0 1 0);
dragDir (1 0 0);
pitchAxis (0 0 1);
magUInf 1e-7;
lRef 0.06; // sphere diameter
Aref 0.0014137; //1/2 * projected area = pi*r²/2

        outputControl   timeStep;
        outputInterval  1;
}

);

//libs ("libOpenFoamTurbo.so");
// ************************************************************************* //
Die Rechnung läuft zwar auch so, aber es wird trotzdem nichts ausgegeben.
Insgesamt bin ich mir nicht sicher, welcher der drei Ansätze der "richtige" ist.
Gruß, Christian

Beiträge: 13
Registriert: 30.04.2011

Beiträge: 272
Registriert: 02.10.2008

erstellt am: 26. Jul. 2011 11:22      <-- editieren / zitieren -->      Unities abgeben: 

calculateTorque-Bin.zip calculateTorque.zip

hier hast ersteinmal die Version "calculateTorque" für OF-2.0.x.

Gruß Thomas

Beiträge: 13
Registriert: 30.04.2011

erstellt am: 26. Jul. 2011 11:26      <-- editieren / zitieren -->      Unities abgeben: 

In der Kräfte-datei stehen: forces(pressure, viscous) und moment(pressure, viscous). Dabei wird für jede dieser Angaben ein Vektor(?) rausgeschrieben. ich vermute es handelt sich dabei um x,y,z.

Meine unbestätigte Vermutung ist: bei dem moment(pressure) handelt es prinzipiell um den Wert, den Ansys mir ausgibt, wenn ich im Function calculator mir das Drehmoment bezüglich einer Achse ausrechnen lasse.
ist das richtig?
Um nun das korrekte Moment zu erhalten, müsste ich doch alle Patches angeben, also alle Flächen orthogonal zur Rotationsachse?(in meinem Fall z)

Für eine kompilierbare Version von "CalculateTorque" wäre ich sehr dankbar 
Gruß
Christian
Edit: Das ging ja schnell! Während ich schrieb hast du schon gepostet! Das ist großartig! Danke! Werde es mir zu gemüte führen!

[Diese Nachricht wurde von ChristianE36 am 26. Jul. 2011 editiert.]

Beiträge: 13
Registriert: 30.04.2011

Beiträge: 15
Registriert: 30.06.2009

Beiträge: 272
Registriert: 02.10.2008

erstellt am: 02. Aug. 2011 15:57      <-- editieren / zitieren -->      Unities abgeben: 

ein Blick im Quellcode:

pDyn = 1/2 * rhoRef * magUInf^2

Cm = pitchMoment/( Aref * lRef * pDyn )

Code:

---

// Free-stream conditions
//- Velocity magnitude

scalar magUInf_;

// Reference scales
//- Length
scalar lRef_;

//- Area
scalar Aref_;

---

Gruß Thomas

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP