src/engine/engineTopoChangerMesh/pistonSliding/setPistonSlidingConstraint.H

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