pd::ReferenceGapSelfContact Class Reference

#include <selfContact.h>

Inheritance diagram for pd::ReferenceGapSelfContact:

Collaboration diagram for pd::ReferenceGapSelfContact:

Public Member Functions
util::Point	force (const util::Point &yji, double volj, double Kn, double Rc, double r0) const override
Public Member Functions inherited from pd::SelfContact
virtual	~SelfContact ()=default

Detailed Description

Bhattacharya–Lipton 2023 §4.4: act only when current separation < Rc. Natural length is r0 when r0 > 0, else Rc (short-range: only when R < Rc).

Definition at line 49 of file selfContact.h.

Member Function Documentation

force()

util::Point pd::ReferenceGapSelfContact::force ( const util::Point &  yji, double  volj, double  Kn, double  Rc, double  r0  ) const

overridevirtual

Force contribution on node i from neighbor j (broken / unbound pair).

Parameters

yji	Current relative position xj+uj - (xi+ui)
volj	Neighbor nodal volume
Kn	Contact stiffness
Rc	Contact radius for this particle (used by broken_bond_kn)
r0	Reference distance |xj - xi| (used by reference_gap)

Implements pd::SelfContact.

Definition at line 48 of file selfContact.cpp.

                                                            {
  // Short-range only: act when current separation < Rc. Natural length is the
  // unbroken reference distance r0 when available (often r0 > Rc for Mode-I
  // crack-crossing bonds); otherwise Rc. Using min(r0, Rc) collapsed to Rc
  // whenever r0 > Rc and made this law identical to broken_bond_kn.
  const double Rji = yji.length();
  if (!(Rji > 0.) || !(Rc > 0.) || !(Rji < Rc))
    return util::Point();
  const double natural = (r0 > 0.) ? r0 : Rc;
  return cappedRepulsive(yji, volj, Kn, natural);
}

References util::Point::length().

Here is the call graph for this function:

---

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

• /home/user/project/src/pd/selfContact.h
• /home/user/project/src/pd/selfContact.cpp