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
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
|
/////////////////////////////////////////////////////////////
// //
// Copyright (c) 2003-2014 by The University of Queensland //
// Centre for Geoscience Computing //
// http://earth.uq.edu.au/centre-geoscience-computing //
// //
// Primary Business: Brisbane, Queensland, Australia //
// Licensed under the Open Software License version 3.0 //
// http://www.apache.org/licenses/LICENSE-2.0 //
// //
/////////////////////////////////////////////////////////////
#ifndef __ADHESIVEFRICTION_H
#define __ADHESIVEFRICTION_H
#include "Model/FrictionInteraction.h"
/*!
\brief Interaction parameters for adhesive frictional interaction
*/
class CAdhesiveFrictionIGP : public AIGParam
{
public:
virtual std::string getTypeString() const {return "AdhesiveFriction";}
void setTimeStepSize(double timeStepSize)
{
this->dt = timeStepSize;
}
double k;
double mu;
double k_s;
double dt;
double r_cut;
};
/*!
\brief Frictional+adhesive Elastic interaction between particles
*/
class CAdhesiveFriction : public CFrictionInteraction
{
public: // types
typedef CAdhesiveFrictionIGP ParameterType;
typedef double (CAdhesiveFriction::* ScalarFieldFunction)() const;
typedef pair<bool,double> (CAdhesiveFriction::* CheckedScalarFieldFunction)() const;
typedef Vec3 (CAdhesiveFriction::* VectorFieldFunction)() const;
protected:
double m_r_cut;
double m_r_cut_h;
public:
CAdhesiveFriction();
CAdhesiveFriction(CParticle*,CParticle*,const CAdhesiveFrictionIGP&);
virtual ~CAdhesiveFriction();
static string getType() {return "AdhesiveFriction";};
virtual void calcForces();
friend class TML_PackedMessageInterface;
};
#endif // __ADHESIVEFRICTION_H
|
