applications/solvers/solidMechanics/icoFsiElasticNonLinULSolidFoam/createStressPointMesh.H

   1
   2 // Create point stress mesh for interpolation
   3 pointMesh pStressMesh(stressMesh);
   4
   5 pointPatchInterpolation patchPointInterpolator(stressMesh);
   6
   7
   8 // Create point displacement field (check fixedValue patches)
   9
  10 wordList types
  11 (
  12     pStressMesh.boundary().size(),
  13     calculatedFvPatchVectorField::typeName
  14 );
  15
  16 // wordList types = DU.boundaryField().types();
  17
  18 forAll(DU.boundaryField().types(), patchI)
  19 {
  20     if
  21     (
  22         DU.boundaryField().types()[patchI]
  23      == fixedValueFvPatchVectorField::typeName
  24     )
  25     {
  26         types[patchI] = fixedValueFvPatchVectorField::typeName;
  27     }
  28 }
  29
  30 pointVectorField pointDU
  31 (
  32     IOobject
  33     (
  34         "pointDU",
  35         runTime.timeName(),
  36         stressMesh
  37     ),
  38     pStressMesh,
  39     dimensionedVector("zero", dimLength, vector::zero),
  40     types
  41 );