FreeFOAM The Cross-Platform CFD Toolkit
Hosted by SourceForge:
Get FreeFOAM at SourceForge.net.
            Fast, secure and Free Open Source software downloads

createFields.H

Go to the documentation of this file.
00001 Info<< "Reading displacement field D\n" << endl;
00002 volVectorField D
00003 (
00004     IOobject
00005     (
00006         "D",
00007         runTime.timeName(),
00008         mesh,
00009         IOobject::MUST_READ,
00010         IOobject::AUTO_WRITE
00011     ),
00012     mesh
00013 );
00014 
00015 Info<< "Reading/setting displacement correction field Dcorr\n" << endl;
00016 volVectorField Dcorr
00017 (
00018     IOobject
00019     (
00020         "Dcorr",
00021         runTime.timeName(),
00022         mesh
00023     ),
00024     D
00025 );
00026 Dcorr *= 0.0;
00027 
00028 Info<< "Calculating stress field sigmaD\n" << endl;
00029 volSymmTensorField sigmaD
00030 (
00031     IOobject
00032     (
00033         "sigmaD",
00034         runTime.timeName(),
00035         mesh,
00036         IOobject::NO_READ,
00037         IOobject::NO_WRITE
00038     ),
00039     mu*twoSymm(fvc::grad(D)) + (lambda*I)*tr(fvc::grad(D))
00040 );
00041 
00042 volTensorField sigmaExp
00043 (
00044     IOobject
00045     (
00046         "sigmaExp",
00047         runTime.timeName(),
00048         mesh,
00049         IOobject::NO_READ,
00050         IOobject::NO_WRITE
00051     ),
00052     (lambda - mu)*fvc::grad(Dcorr) + mu*fvc::grad(Dcorr)().T()
00053   + (lambda*I)*tr(fvc::grad(Dcorr))
00054 );
00055 
00056 // ************************ vim: set sw=4 sts=4 et: ************************ //
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines