contact::AnalyticalPlaneWallContact Class Reference

#include <wallContact.h>

Inheritance diagram for contact::AnalyticalPlaneWallContact:

Collaboration diagram for contact::AnalyticalPlaneWallContact:

Public Member Functions
void	apply (data::ModelData &data, PairForce *pair, bool use_node_damping) const override
bool	skipsMeshedGrainWall () const override
Public Member Functions inherited from contact::WallContact
virtual	~WallContact ()=default

Detailed Description

Definition at line 49 of file wallContact.h.

Member Function Documentation

apply()

void contact::AnalyticalPlaneWallContact::apply ( data::ModelData &  data, PairForce *  pair, bool  use_node_damping  ) const

overridevirtual

Implements contact::WallContact.

Definition at line 29 of file wallContact.cpp.

                                                                    {
  if (data.d_particlesListTypeWall.empty())
    return;
 
  const unsigned n_workers = util::parallel::getNThreads();
  tf::Executor executor(n_workers);
  tf::Taskflow taskflow;
 
  taskflow.for_each_index(
      (std::size_t)0, data.d_fContCompNodes.size(), (std::size_t)1,
      [&data, use_node_damping](std::size_t II) {
        const auto i = data.d_fContCompNodes[II];
        const auto &ptIdi = data.getPtId(i);
        auto &pi = data.getParticleFromAllList(ptIdi);
        if (pi->isWall())
          return;
 
        const auto &yi = data.d_x[i];
        const auto &vi = data.d_v[i];
        const double voli = data.d_vol[i];
        util::Point force_i;
 
        for (auto *wall : data.d_particlesListTypeWall) {
          if (!wall || !wall->d_geom_p)
            continue;
 
          geom::WallContactHit hit;
          if (!wall->d_geom_p->wallContactQuery(yi, hit) || !hit.active)
            continue;
 
          const auto &contact =
              data.d_particleDeck_p->d_contactDeck.getContact(
                  pi->getGroupId("contact_id"),
                  wall->getGroupId("contact_id"));
 
          // Signed gap from geom (positive = free space). Contact when gap < Rc.
          const double R = hit.signed_gap;
          if (!util::isLess(R, contact.d_contactR))
            continue;
 
          // en points from grain toward wall (into the material).
          const util::Point en = -1. * hit.outward_n;
          auto scalar_f = contact.d_Kn * (R - contact.d_contactR) * voli;
          if (scalar_f > 0.)
            scalar_f = 0.;
 
          util::Point f = scalar_f * en;
 
          if (contact.d_frictionOn) {
            const util::Point vji = -1. * vi; // wall velocity ~ 0
            const double vn = vji * en;
            util::Point et = vji - vn * en;
            if (util::isGreater(et.length(), 0.))
              et = et / et.length();
            else
              et = util::Point();
            f += contact.d_mu * scalar_f * et;
          }
 
          if (use_node_damping && contact.d_dampingOn && voli > 0. &&
              contact.d_K > 0. && contact.d_contactR > 0.) {
            const util::Point vji = -1. * vi;
            const double vn = vji * en;
            if (util::isLess(vn, 0.)) {
              const double meq =
                  util::equivalentMass(pi->getDensity() * voli,
                                       wall->getDensity() * voli);
              const double beta_n =
                  contact.d_betan *
                  std::sqrt(contact.d_K * contact.d_contactR * meq);
              f += (beta_n * vn / voli) * en;
            }
          }
 
          force_i += f;
        }
 
        data.d_f[i] += force_i;
      });
 
  executor.run(taskflow).get();
}

References util::parallel::getNThreads().

Here is the call graph for this function:

skipsMeshedGrainWall()

bool contact::AnalyticalPlaneWallContact::skipsMeshedGrainWall ( ) const

inlineoverridevirtual

True when meshed grain–wall pairs must be skipped in Contact::computeForces.

Reimplemented from contact::WallContact.

Definition at line 53 of file wallContact.h.

{ return true; }

---

The documentation for this class was generated from the following files:

• /home/user/project/src/contact/wallContact.h
• /home/user/project/src/contact/wallContact.cpp