MercuryDPM  Alpha
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
ParhamiMcMeekingSinterInteraction.cc
Go to the documentation of this file.
1 //Copyright (c) 2013-2014, 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 
29 #include "Particles/BaseParticle.h"
30 #include "InteractionHandler.h"
31 #include <iomanip>
32 #include <fstream>
33 
40  : BaseInteraction(P, I, timeStamp)
41 {
42 #ifdef DEBUG_CONSTRUCTOR
43  std::cout<<"ParhamiMcMeekingSinterInteraction::ParhamiMcMeekingSinterInteraction() finished"<<std::endl;
44 #endif
45 }
46 
51  : BaseInteraction(p)
52 {
53 #ifdef DEBUG_CONSTRUCTOR
54  std::cout<<"ParhamiMcMeekingSinterInteraction::ParhamiMcMeekingSinterInteraction(const ParhamiMcMeekingSinterInteraction &p finished"<<std::endl;
55 #endif
56 }
57 
62 {
63 #ifdef DEBUG_DESTRUCTOR
64  std::cout<<"ParhamiMcMeekingSinterInteraction::~ParhamiMcMeekingSinterInteraction() finished"<<std::endl;
65 #endif
66 }
67 
72 {}
73 
78 {}
83 {
84  //std::cout << "ParhamiMcMeekingSinterInteraction::computeAdhesionForce" << std::endl;
85  const SpeciesType* species = getSpecies();
86  Mdouble effectiveDiameter = 2.0*getEffectiveRadius(); // effectiveRadius = (r1*r2)/(r1+r2)
87  Mdouble contactRadiusSquared = 2.0 * effectiveDiameter * getOverlap();
88  Vec3D tangentialRelativeVelocity = getRelativeVelocity() - getNormal() * getNormalRelativeVelocity();
89 
90  //viscous force is viscosityCoefficient_*contactRadius^4*normalRelativeVelocity
91  //adhesion force is adhesionCoefficient_*radius
92  Mdouble normalForce = - species->getAdhesionCoefficient() * effectiveDiameter
93  - species->getViscosityCoefficient()*contactRadiusSquared*contactRadiusSquared * getNormalRelativeVelocity();
94  //tangential force is slidingFrictionCoefficient_*contactRadius^2*radius*tangentialRelativeVelocity
95  Vec3D tangentialForce = - species->getSlidingFrictionCoefficient()*contactRadiusSquared*effectiveDiameter*tangentialRelativeVelocity;
96  //std::cout << "P" << species->getAdhesionCoefficient() << species->getViscosityCoefficient() << species->getSlidingFrictionCoefficient() << std::endl;
97  Mdouble attractiveForce = - species->getAdhesionCoefficient() * effectiveDiameter;
98  //std::cout << effectiveDiameter << "Fs=" << species->getAdhesionCoefficient() * effectiveDiameter
99  //<< " Fv/del/del/del=" << 4.0*species->getViscosityCoefficient()*effectiveDiameter*effectiveDiameter << std::endl;
100  //<< species->getViscosityCoefficient() << " d" <<contactRadiusSquared<< " d" <<getNormalRelativeVelocity() << std::endl;
101  addForce(getNormal() * normalForce + tangentialForce);
102 }
103 
109 {
110  return 0.0;
111 }
116 {
117  return dynamic_cast<const SpeciesType*> (getBaseSpecies()); //downcast
118 }
123 {
124  return "ParhamiMcMeekingSinter";
125 }
Mdouble getEffectiveRadius() const
Returns a Mdouble to the effective radius of the interaction. (Not corrected for the overlap) ...
ParhamiMcMeekingSinterInteraction(BaseInteractable *P, BaseInteractable *I, Mdouble timeStamp)
Constructor.
void write(std::ostream &os) const
Interaction print function, which accepts an std::ostream as input.
const Vec3D & getRelativeVelocity() const
Returns a constant reference to a vector of relative velocity.
double Mdouble
Stores information about interactions between two interactable objects; often particles but could be ...
ParhamiMcMeekingSinterSpecies contains the parameters used to describe a linear reversible short-rang...
Mdouble getNormalRelativeVelocity() const
Returns a double which is the norm (length) of the relative velocity vector.
std::string getBaseName() const
Returns the name of the interaction, see Interaction.h.
const BaseSpecies * getBaseSpecies() const
Return a constant point to BaseSpecies of the interaction.
const Vec3D & getNormal() const
Gets the normal vector between the two interacting objects.
#define UNUSED
Definition: GeneralDefine.h:39
Mdouble getOverlap() const
Returns a Mdouble with the current overlap between the two interacting objects.
void addForce(Vec3D force)
add an force increment to the total force.
const SpeciesType * getSpecies() const
A dynamic_cast of BaseSpecies pointer type to a pointer to an object of type ReversibleAdhesiveSpecie...
void computeAdhesionForce()
Computes the adhesive forces.
Defines the basic properties that a interactable object can have.
Implementation of a 3D vector (by Vitaliy).
Definition: Vector.h:45
void read(std::istream &is)
Interaction read function, which accepts an std::istream as input.
Mdouble getElasticEnergy() const
Returns the amount of Elastic energy involved in an interaction. Basically used in case you want to w...