src/engine/engineTopoChangerMesh/accordionEngineMesh/setAccordionEngineMeshConstraintsNoDeformation.H

   1 {
   2
   3     const faceZoneMesh& fZones = faceZones();
   4     const pointZoneMesh& pZones = pointZones();
   5
   6     label constraintSize = 0;
   7
   8     forAll(valves_, valveI)
   9     {
  10
  11         vector valveVel =
  12             valves_[valveI].curVelocity()*valves_[valveI].cs().axis();
  13
  14         {
  15
  16             label movingPtsIndex = pZones.findZoneID("movingPointsV"+Foam::name(valveI + 1));
  17             const labelList& movingPointsV = pZones[movingPtsIndex];
  18
  19             forAll(movingPointsV, mpI)
  20             {
  21                 constrainedPoints.append(movingPointsV[mpI]);
  22                 constrainedVelocity.append(valveVel);
  23                 constraintSize++;
  24             }
  25
  26         }
  27
  28         {
  29
  30             label staticPtsIndex = pZones.findZoneID("staticPointsV"+Foam::name(valveI + 1));
  31             const labelList& staticPointsV = pZones[staticPtsIndex];
  32
  33             forAll(staticPointsV, spI)
  34             {
  35                 constrainedPoints.append(staticPointsV[spI]);
  36                 constrainedVelocity.append(vector::zero);
  37                 constraintSize++;
  38             }
  39
  40         }
  41
  42         {
  43
  44             label movingFacesIndex = fZones.findZoneID("movingFacesZoneV"+Foam::name(valveI + 1));
  45             const labelList& movingFacesV = pZones[movingFacesIndex];
  46
  47             forAll(movingFacesV, mfI)
  48             {
  49                 constrainedPoints.append(movingFacesV[mfI]+points().size());
  50                 constrainedVelocity.append(valveVel);
  51                 constraintSize++;
  52             }
  53
  54         }
  55
  56         {
  57
  58             label staticFacesIndex = fZones.findZoneID("staticFacesZoneV"+Foam::name(valveI + 1));
  59             const labelList& staticFacesV = pZones[staticFacesIndex];
  60
  61             forAll(staticFacesV, sfI)
  62             {
  63                 constrainedPoints.append(staticFacesV[sfI]+points().size());
  64                 constrainedVelocity.append(vector::zero);
  65                 constraintSize++;
  66             }
  67
  68         }
  69
  70     }
  71
  72 }