applications/solvers/solidMechanics/elasticOrthoNonLinULSolidFoam/writeFields.H

   1 if (runTime.outputTime())
   2 {
   3     //C.write();
   4
   5     volScalarField epsilonEq
   6     (
   7         IOobject
   8         (
   9             "epsilonEq",
  10             runTime.timeName(),
  11             mesh,
  12             IOobject::NO_READ,
  13             IOobject::AUTO_WRITE
  14         ),
  15         sqrt((2.0/3.0)*magSqr(dev(epsilon)))
  16     );
  17
  18     Info<< "Max epsilonEq = " << max(epsilonEq).value() << endl;
  19
  20     volScalarField sigmaEq
  21     (
  22         IOobject
  23         (
  24             "sigmaEq",
  25             runTime.timeName(),
  26             mesh,
  27             IOobject::NO_READ,
  28             IOobject::AUTO_WRITE
  29         ),
  30         sqrt((3.0/2.0)*magSqr(dev(sigma)))
  31     );
  32
  33     Info<< "Max sigmaEq = " << max(sigmaEq).value() << endl;
  34
  35 // volVectorField traction
  36 // (
  37 //     IOobject
  38 //     (
  39 //         "traction",
  40 //         runTime.timeName(),
  41 //         mesh,
  42 //         IOobject::NO_READ,
  43 //         IOobject::AUTO_WRITE
  44 //     ),
  45 //     mesh,
  46 //     dimensionedVector("zero", dimForce/dimArea, vector::zero)
  47 // );
  48 // forAll(mesh.boundary(), patchi)
  49 // {
  50 //     traction.boundaryField()[patchi] =
  51 //         n.boundaryField()[patchi] & sigma.boundaryField()[patchi];
  52 // }
  53
  54 //     //- patch forces
  55 //     forAll(mesh.boundary(), patchi)
  56 //     {
  57 //         Info << "Patch " << mesh.boundary()[patchi].name() << endl;
  58 //         vectorField totalForce = mesh.Sf().boundaryField()[patchi] & sigma.boundaryField()[patchi];
  59 //         vector force = sum( totalForce );
  60 //         Info << "\ttotal force is " << force << " N" << endl;
  61 //     tensorField F = I + gradDU.boundaryField()[patchi];
  62 //         tensorField Finv = inv(F);
  63 //         scalar normalForce = sum( n.boundaryField()[patchi] & totalForce );
  64 //         Info << "\tnormal force is " << normalForce << " N" << endl;
  65 //         scalar shearForce = mag(sum( (I - sqr(n.boundaryField()[patchi])) & totalForce ));
  66 //         Info << "\tshear force is " << shearForce << " N" << endl;
  67 //     }
  68
  69     runTime.write();
  70 }