applications/solvers/solidMechanics/elasticAcpSolidFoam/writeFields.H

   1 if (runTime.outputTime() || topoChange)
   2 {
   3     volScalarField sigmaEq
   4     (
   5         IOobject
   6         (
   7             "sigmaEq",
   8             runTime.timeName(),
   9             mesh,
  10             IOobject::NO_READ,
  11             IOobject::AUTO_WRITE
  12         ),
  13         sqrt((3.0/2.0)*magSqr(dev(sigma)))
  14     );
  15
  16     Info<< "Max sigmaEq = " << max(sigmaEq).value()
  17         << endl;
  18
  19     volScalarField epsilonEq
  20     (
  21         IOobject
  22         (
  23             "epsilonEq",
  24             runTime.timeName(),
  25             mesh,
  26             IOobject::NO_READ,
  27             IOobject::AUTO_WRITE
  28         ),
  29         sqrt((2.0/3.0)*magSqr(dev(epsilon)))
  30     );
  31
  32     Info<< "Max epsilonEq = " << max(epsilonEq).value()
  33         << endl;
  34
  35 //     Info << "\nCalculate maximal principal stress ..." << flush;
  36 //     // Principal stresses
  37 //     volVectorField sigmaMax
  38 //     (
  39 //         IOobject
  40 //         (
  41 //             "sigmaMax",
  42 //             runTime.timeName(),
  43 //             mesh,
  44 //             IOobject::NO_READ,
  45 //             IOobject::AUTO_WRITE
  46 //         ),
  47 //         mesh,
  48 //         dimensionedVector("sigmaMax", dimPressure, vector::zero)
  49 //     );
  50 //     vectorField& sigmaMaxI = sigmaMax.internalField();
  51
  52 //     forAll (sigmaMaxI, cellI)
  53 //     {
  54 //         vector eValues = eigenValues(sigma.internalField()[cellI]);
  55 //         tensor eVectors = eigenVectors(sigma.internalField()[cellI]);
  56
  57 //         scalar maxEValue = 0;
  58 //         label iMax = -1;
  59 //         forAll(eValues, i)
  60 //         {
  61 //             if (eValues[i] > maxEValue)
  62 //             {
  63 //                 maxEValue = eValues[i];
  64 //                 iMax = i;
  65 //             }
  66 //         }
  67
  68 //         if (iMax != -1)
  69 //         {
  70 //             if (iMax == 0)
  71 //             {
  72 //                 sigmaMaxI[cellI] = eVectors.x()*eValues.x();
  73 //             }
  74 //             else if (iMax == 1)
  75 //             {
  76 //                 sigmaMaxI[cellI] = eVectors.y()*eValues.y();
  77 //             }
  78 //             else if (iMax == 2)
  79 //             {
  80 //                 sigmaMaxI[cellI] = eVectors.z()*eValues.z();
  81 //             }
  82 //         }
  83 //     }
  84
  85
  86     //- cohesive damage and cracking, and GII and GII
  87     volScalarField damageAndCracks
  88       (
  89        IOobject
  90        (
  91         "damageAndCracks",
  92         runTime.timeName(),
  93         mesh,
  94         IOobject::NO_READ,
  95         IOobject::AUTO_WRITE
  96         ),
  97        mesh,
  98        dimensionedScalar("zero", dimless, 0.0),
  99        calculatedFvPatchVectorField::typeName
 100        );
 101     volScalarField GI
 102       (
 103        IOobject
 104        (
 105         "GI",
 106         runTime.timeName(),
 107         mesh,
 108         IOobject::NO_READ,
 109         IOobject::AUTO_WRITE
 110         ),
 111        mesh,
 112        dimensionedScalar("zero", dimless, 0.0),
 113        calculatedFvPatchVectorField::typeName
 114        );
 115     volScalarField GII
 116       (
 117        IOobject
 118        (
 119         "GII",
 120         runTime.timeName(),
 121         mesh,
 122         IOobject::NO_READ,
 123         IOobject::AUTO_WRITE
 124         ),
 125        mesh,
 126        dimensionedScalar("zero", dimless, 0.0),
 127        calculatedFvPatchVectorField::typeName
 128        );
 129     forAll(U.boundaryField(), patchi)
 130       {
 131         // if(U.boundaryField()[patchi].type() == cohesiveLawMultiMatFvPatchVectorField::typeName)
 132         if(U.boundaryField()[patchi].type() == solidCohesiveFvPatchVectorField::typeName)
 133           {
 134             // cohesiveLawMultiMatFvPatchVectorField& Upatch =
 135             //   refCast<cohesiveLawMultiMatFvPatchVectorField>(U.boundaryField()[patchi]);
 136             solidCohesiveFvPatchVectorField& Upatch =
 137               refCast<solidCohesiveFvPatchVectorField>(U.boundaryField()[patchi]);
 138
 139             GI.boundaryField()[patchi] = Upatch.GI();
 140             GII.boundaryField()[patchi] = Upatch.GII();
 141             damageAndCracks.boundaryField()[patchi] = Upatch.crackingAndDamage();
 142           }
 143       }
 144     volScalarField GTotal("GTotal", GI + GII);
 145     GTotal.write();
 146
 147
 148     runTime.writeNow();
 149 }