applications/solvers/solidMechanics/elasticSolidFoam/writeFields.H

   1 if (runTime.outputTime())
   2   {
   3     volScalarField epsilonEq
   4       (
   5        IOobject
   6        (
   7     "epsilonEq",
   8     runTime.timeName(),
   9     mesh,
  10     IOobject::NO_READ,
  11     IOobject::AUTO_WRITE
  12     ),
  13        sqrt((2.0/3.0)*magSqr(dev(epsilon)))
  14        );
  15
  16     Info<< "Max epsilonEq = " << max(epsilonEq).value()
  17     << endl;
  18
  19     volScalarField sigmaEq
  20       (
  21        IOobject
  22        (
  23     "sigmaEq",
  24     runTime.timeName(),
  25     mesh,
  26     IOobject::NO_READ,
  27     IOobject::AUTO_WRITE
  28     ),
  29        sqrt((3.0/2.0)*magSqr(dev(sigma)))
  30        );
  31
  32     Info<< "Max sigmaEq = " << max(sigmaEq).value()
  33     << endl;
  34
  35     //- boundary traction
  36     // volVectorField traction
  37     //   (
  38     //    IOobject
  39     //    (
  40     //     "traction",
  41     //     runTime.timeName(),
  42     //     mesh,
  43     //     IOobject::NO_READ,
  44     //     IOobject::AUTO_WRITE
  45     //     ),
  46     //    mesh,
  47     //    dimensionedVector("zero", dimForce/dimArea, vector::zero)
  48     //    );
  49     // forAll(traction.boundaryField(), patchi)
  50     //   {
  51     //     traction.boundaryField()[patchi] =
  52     //       n.boundaryField()[patchi] & sigma.boundaryField()[patchi];
  53     //   }
  54
  55     // contact fields
  56     // if(contactBoundaries)
  57     //   {
  58     //     volScalarField stickSlipFaces
  59     //       (
  60     //        IOobject
  61     //        (
  62     //         "stickSlipFaces",
  63     //         runTime.timeName(),
  64     //         mesh,
  65     //         IOobject::NO_READ,
  66     //         IOobject::AUTO_WRITE
  67     //         ),
  68     //        mesh,
  69     //        dimensionedScalar("zero", dimless, 0.0)
  70     //        );
  71     //     forAll(stickSlipFaces.boundaryField(), patchi)
  72     //       {
  73     //         if(U.boundaryField()[patchi].type()
  74     //            ==
  75     //            solidContactFvPatchVectorField::typeName)
  76     //           {
  77     //         const solidContactFvPatchVectorField& Upatch =
  78     //           refCast<const solidContactFvPatchVectorField>
  79     //           (U.boundaryField()[patchi]);
  80
  81     //         if(!Upatch.master())
  82     //           {
  83     //             stickSlipFaces.boundaryField()[patchi] =
  84     //               Upatch.frictionContactModelPtr()->stickSlipFaces();
  85     //           }
  86     //           }
  87     //       }
  88     //     stickSlipFaces.write();
  89     //   }
  90
  91
  92     //- boundary forces
  93 //     Info << nl;
  94 //     forAll(mesh.boundary(), patchi)
  95 //       {
  96 //     Info << "Patch " << mesh.boundary()[patchi].name() << endl;
  97 //     vectorField totalForce = mesh.Sf().boundaryField()[patchi] & sigma.boundaryField()[patchi];
  98
  99 //     vector force = sum( totalForce );
 100 //     Info << "\ttotal force is " << force << " N" << endl;
 101
 102 //     const vectorField& nb = n.boundaryField()[patchi];
 103 //     scalar normalForce = sum( nb & totalForce );
 104 //     Info << "\tnormal force is " << normalForce << " N" << endl;
 105 //     scalar shearForce = mag(sum( (I - sqr(nb)) & totalForce ));
 106 //     Info << "\tshear force is " << shearForce << " N" << endl;
 107 //     }
 108
 109     runTime.write();
 110   }