1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
|
#include <pkg/dem/CapillaryPhys.hpp>
#include <pkg/dem/ScGeom.hpp>
namespace yade { // Cannot have #include directive inside.
YADE_PLUGIN((CapillaryPhys)(Ip2_FrictMat_FrictMat_CapillaryPhys));
// The following code was moved from Ip2_FrictMat_FrictMat_CapillaryPhys.cpp
void Ip2_FrictMat_FrictMat_CapillaryPhys::go(
const shared_ptr<Material>& b1 //FrictMat
,
const shared_ptr<Material>& b2 // FrictMat
,
const shared_ptr<Interaction>& interaction)
{
ScGeom* geom = YADE_CAST<ScGeom*>(interaction->geom.get());
if (geom) {
if (!interaction->phys) {
const shared_ptr<FrictMat>& sdec1 = YADE_PTR_CAST<FrictMat>(b1);
const shared_ptr<FrictMat>& sdec2 = YADE_PTR_CAST<FrictMat>(b2);
if (!interaction->phys) interaction->phys = shared_ptr<CapillaryPhys>(new CapillaryPhys());
const shared_ptr<CapillaryPhys>& contactPhysics = YADE_PTR_CAST<CapillaryPhys>(interaction->phys);
Real Ea = sdec1->young;
Real Eb = sdec2->young;
Real Va = sdec1->poisson;
Real Vb = sdec2->poisson;
Real Da = geom->radius1; // FIXME - multiply by factor of sphere interaction distance (so sphere interacts at bigger range that its geometrical size)
Real Db = geom->radius2; // FIXME - as above
Real fa = sdec1->frictionAngle;
Real fb = sdec2->frictionAngle;
Real Kn = 2 * Ea * Da * Eb * Db / (Ea * Da + Eb * Db); //harmonic average of two stiffnesses
Real Ks = 2 * Ea * Da * Va * Eb * Db * Vb
/ (Ea * Da * Va + Eb * Db * Va); //harmonic average of two stiffnesses with ks=V*kn for each sphere
contactPhysics->tangensOfFrictionAngle = math::tan(math::min(fa, fb));
contactPhysics->kn = Kn;
contactPhysics->ks = Ks;
}
}
};
} // namespace yade
|
