applications/solvers/solidMechanics/elasticIncrSolidFoam/calculateNetForces.H

   1 //- net forces on boundary patches
   2
   3 vectorField netForces(mesh.boundary().size(), vector::zero);
   4
   5 Info << nl;
   6 forAll(netForces, patchI)
   7 {
   8   netForces[patchI] = gSum(mesh.Sf().boundaryField()[patchI] & sigma.boundaryField()[patchI]);
   9
  10   Info << "patch\t" << mesh.boundary()[patchI].name() << "\t\tnet force is\t"
  11        << netForces[patchI] << " N" << endl;
  12 }