MindlinRollingTorsionInteraction.h
Go to the documentation of this file.
1 //Copyright (c) 2013-2023, The MercuryDPM Developers Team. All rights reserved.
2 //For the list of developers, see <http://www.MercuryDPM.org/Team>.
3 //
4 //Redistribution and use in source and binary forms, with or without
5 //modification, are permitted provided that the following conditions are met:
6 // * Redistributions of source code must retain the above copyright
7 // notice, this list of conditions and the following disclaimer.
8 // * Redistributions in binary form must reproduce the above copyright
9 // notice, this list of conditions and the following disclaimer in the
10 // documentation and/or other materials provided with the distribution.
11 // * Neither the name MercuryDPM nor the
12 // names of its contributors may be used to endorse or promote products
13 // derived from this software without specific prior written permission.
14 //
15 //THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND
16 //ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
17 //WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
18 //DISCLAIMED. IN NO EVENT SHALL THE MERCURYDPM DEVELOPERS TEAM BE LIABLE FOR ANY
19 //DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES
20 //(INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
21 //LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND
22 //ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
23 //(INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
24 //SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
25 
26 #ifndef MindlinRollingTorsionInteraction_H
27 #define MindlinRollingTorsionInteraction_H
28 
29 #include "MindlinInteraction.h"
30 #include "Math/Vector.h"
31 
32 class BaseParticle;
33 
35 
36 class BaseInteractable;
37 
45 {
46 public:
51 
56 
57  //used for mpi
59 
64 
69 
74  void computeFrictionForce();
75 
79  void read(std::istream& is) override;
80 
84  void write(std::ostream& os) const override;
85 
90  void integrate(Mdouble timeStep) override;
91 
95  Mdouble getElasticEnergy() const override;
96 
100  std::string getBaseName() const;
101 
106 
112  void reverseHistory() override;
113 
114  void rotateHistory(Matrix3D& rotationMatrix) override;
115 
116  Vec3D getRollingSpring() const;
117 
118  Vec3D getTorsionSpring() const;
119 
120  void setRollingSpring(Vec3D rollingSpring);
121 
122  void setTorsionSpring(Vec3D torsionSpring);
123 
124 private:
145 };
146 
147 #endif
Defines the basic properties that a interactable object can have.
Definition: BaseInteractable.h:55
Definition: BaseParticle.h:54
Implementation of a 3D matrix.
Definition: Matrix.h:38
Computes the forces corresponding to sliding friction.
Definition: MindlinInteraction.h:43
This class allows one to take all three types of frictional interactions into account....
Definition: MindlinRollingTorsionInteraction.h:45
~MindlinRollingTorsionInteraction() override
Destructor.
Definition: MindlinRollingTorsionInteraction.cc:78
const MindlinRollingTorsionSpecies * getSpecies() const
Returns a const pointer of type FrictionSpecies*.
Definition: MindlinRollingTorsionInteraction.cc:233
Vec3D getRollingSpring() const
Definition: MindlinRollingTorsionInteraction.cc:271
MindlinRollingTorsionInteraction()
Definition: MindlinRollingTorsionInteraction.cc:52
Vec3D rollingSpringVelocity_
Stores the rate at which the rolling spring compresses or relaxes. Set in computeFrictionForce(),...
Definition: MindlinRollingTorsionInteraction.h:134
void integrate(Mdouble timeStep) override
Computes the amount of compression in all the springs, i.e., increments the rollingSpring_,...
Definition: MindlinRollingTorsionInteraction.cc:213
void read(std::istream &is) override
Interaction read function, which accepts an std::istream as input.
Definition: MindlinRollingTorsionInteraction.cc:98
void rotateHistory(Matrix3D &rotationMatrix) override
When periodic particles are used, some interactions need certain history properties rotated (e....
Definition: MindlinRollingTorsionInteraction.cc:262
void reverseHistory() override
A useful feature if one wants to return to the initial state of the springs. However,...
Definition: MindlinRollingTorsionInteraction.cc:250
void computeFrictionForce()
Computes the forces arising due to all three types of friction, i.e., sliding, rolling and torsional.
Definition: MindlinRollingTorsionInteraction.cc:109
Mdouble getElasticEnergy() const override
Returns the global amount of energy stored in all the springs (rolling, sliding and torsional).
Definition: MindlinRollingTorsionInteraction.cc:223
std::string getBaseName() const
Returns interaction name/type.
Definition: MindlinRollingTorsionInteraction.cc:242
void setTorsionSpring(Vec3D torsionSpring)
Definition: MindlinRollingTorsionInteraction.cc:286
void write(std::ostream &os) const override
Interaction print function, which accepts an std::ostream as input.
Definition: MindlinRollingTorsionInteraction.cc:88
Vec3D getTorsionSpring() const
Definition: MindlinRollingTorsionInteraction.cc:276
Vec3D rollingSpring_
Stores the amount of rolling spring compression. Set in integrate(), used in computing frictional for...
Definition: MindlinRollingTorsionInteraction.h:129
Vec3D torsionSpring_
Stores the amount of torsional spring compression. Set in integrate(), used in computing frictional f...
Definition: MindlinRollingTorsionInteraction.h:139
MindlinRollingTorsionSpecies SpeciesType
An alias for FrictionSpecies.
Definition: MindlinRollingTorsionInteraction.h:50
void setRollingSpring(Vec3D rollingSpring)
Definition: MindlinRollingTorsionInteraction.cc:281
Vec3D torsionSpringVelocity_
Stores the rate at which the torsional spring compresses or relaxes. Set in computeFrictionForce(),...
Definition: MindlinRollingTorsionInteraction.h:144
MindlinRollingTorsionSpecies contains the parameters used to describe sliding, rolling and torsional ...
Definition: MindlinRollingTorsionSpecies.h:39
Definition: Vector.h:51
double P
Uniform pressure.
Definition: TwenteMeshGluing.cpp:73