src/engine/engineTopoChangerMesh/pistonSliding/setPistonSlidingConstraint.H

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