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