MercuryDPM  Trunk
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
BondedInteraction.cc
Go to the documentation of this file.
1 //Copyright (c) 2013-2020, 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 
27 #include "BondedInteraction.h"
29 #include "Particles/BaseParticle.h"
30 #include "InteractionHandler.h"
31 #include <iomanip>
32 #include <fstream>
33 
40  : BaseInteraction(P, I, timeStamp)
41 {
42  bonded_ = false;
43 #ifdef DEBUG_CONSTRUCTOR
44  std::cout<<"BondedInteraction::BondedInteraction() finished"<<std::endl;
45 #endif
46 }
47 
52  : BaseInteraction(p)
53 {
55  bonded_ = p.bonded_;
56 #ifdef DEBUG_CONSTRUCTOR
57  std::cout<<"BondedInteraction::BondedInteraction(const BondedInteraction &p finished"<<std::endl;
58 #endif
59 }
60 
62 {
63 #ifdef MERCURY_USE_MPI
64  logger(FATAL,"ChargedBondedInteractions are currently not implemented in parallel MercuryDPM");
65 #endif
66 }
67 
72 {
73 #ifdef DEBUG_DESTRUCTOR
74  std::cout<<"BondedInteraction::~BondedInteraction() finished"<<std::endl;
75 #endif
76 }
77 
81 void BondedInteraction::write(std::ostream& os) const
82 {
83  os << " bonded " << bonded_;
84 }
85 
89 void BondedInteraction::read(std::istream& is)
90 {
91  std::string dummy;
92  is >> dummy >> bonded_;
93 }
94 
99 {
100  const BondedSpecies* species = getSpecies();
101  if (bonded_ && getOverlap() >= 0)
102  {
103  addForce(getNormal() * (-species->getBondForceMax()
105  }
106 }
107 
112 {
113  return static_cast<const BondedSpecies*>(getBaseSpecies()->getAdhesiveForce());
114 ;
115 }
116 
121 {
122  return "Bonded";
123 }
124 
131 {
132  if (!bonded_)
133  return 0.0;
134  else
135  return -getSpecies()->getBondForceMax() * getOverlap();
136 }
137 
139 {
140  return bonded_;
141 }
142 
144 {
145  bonded_ = bonded;
146 }
147 
148 
150 {
151  bonded_ = true;
152 }
153 
155 {
156  bonded_ = false;
157 }
Mdouble getBondDissipation() const
Allows the spring constant to be accessed.
BondedSpecies contains the parameters used to describe a linear irreversible short-range force...
Definition: BondedSpecies.h:38
void write(std::ostream &os) const override
Interaction print function, which accepts an std::ostream as input.
~BondedInteraction() override
Destructor.
Logger< MERCURY_LOGLEVEL > logger("MercuryKernel")
double Mdouble
Definition: GeneralDefine.h:34
BaseAdhesiveForce * getAdhesiveForce() const
Definition: BaseSpecies.h:152
const BondedSpecies * getSpecies() const
Returns a pointer to the adhesive force species BondedSpecies.
Stores information about interactions between two interactable objects; often particles but could be ...
const Vec3D & getNormal() const
Gets the normal vector between the two interacting objects.
Mdouble getNormalRelativeVelocity() const
Returns a double which is the norm (length) of the relative velocity vector.
const BaseSpecies * getBaseSpecies() const
Return a constant point to BaseSpecies of the interaction.
Mdouble getElasticEnergy() const override
Returns the elastic energy stored in the adhesive spring.
Mdouble getOverlap() const
Returns a Mdouble with the current overlap between the two interacting objects.
bool getBonded() const
Get function.
void addForce(Vec3D force)
add an force increment to the total force.
void setBonded(bool bonded)
Set function.
std::string getBaseName() const
Returns the name of the interaction, see Interaction.h.
void computeAdhesionForce()
Computes the Adhesive force.
Defines the basic properties that a interactable object can have.
Mdouble getBondForceMax() const
Allows bondForceMax_ to be accessed.
void read(std::istream &is) override
Interaction read function, which accepts an std::istream as input.
bool bonded_
A history parameter to store if the particles were in contact or not. Useful to compute adhesive forc...