MercuryDPM  Trunk
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
NurbsWall.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 #include "NurbsWall.h"
27 #include "array"
28 #include "InteractionHandler.h"
29 #include "WallHandler.h"
30 #include "DPMBase.h"
31 #include "Math/ExtendedMath.h"
32 #include "Particles/BaseParticle.h"
33 using mathsFunc::square;
40 {
41  logger(DEBUG, "NurbsWall() constructor finished.");
42 }
43 
48  : BaseWall(other), nurbsSurface_(other.nurbsSurface_)
49 {
50  logger(DEBUG, "NurbsWall(const NurbsWall&) copy constructor finished.");
51 }
52 
53 NurbsWall::NurbsWall(const NurbsSurface& nurbsSurface)
54  : nurbsSurface_(nurbsSurface)
55 {
56  logger(DEBUG, "NurbsWall(const NurbsSurface&) constructor finished.");
57 }
58 
60 {
61  logger(DEBUG, "~NurbsWall() finished, destroyed the wall.");
62 }
63 
65 {
66  return new NurbsWall(*this);
67 }
68 
69 void NurbsWall::set(const NurbsSurface& nurbsSurface) {
70  nurbsSurface_ = nurbsSurface;
71 }
72 
73 
74 bool NurbsWall::getDistanceAndNormal(const BaseParticle& p, Mdouble& distance, Vec3D& normal_return) const
75 {
76  //transform coordinates into position-orientation frame
77  Vec3D position = p.getPosition() - getPosition();
78  getOrientation().rotateBack(position);
80  if (nurbsSurface_.getDistance(position, p.getRadius()+s->getInteractionDistance(), distance, normal_return)) {
81  getOrientation().rotate(normal_return);
82  return true;
83  } else {
84  return false;
85  }
86 }
87 
91 void NurbsWall::read(std::istream& is)
92 {
93  BaseWall::read(is);
94  std::string dummy;
95  is >> dummy;
96  is >> nurbsSurface_;
97 }
98 
102 void NurbsWall::write(std::ostream& os) const
103 {
104  BaseWall::write(os);
105  os << " NurbsSurface " << nurbsSurface_;
106 }
107 
111 std::string NurbsWall::getName() const
112 {
113  return "NurbsWall";
114 }
115 
117 {
118  unsigned nu = 30;
119  unsigned nv = 30;
120 
121  //create points
122  size_t nPoints = vtk.points.size();
123  for (double u=0; u<=nu; u++) {
124  for (double v=0; v<=nv; v++) {
125  Vec3D p = nurbsSurface_.evaluate(u/nu,v/nu);
126  getOrientation().rotate(p);
127  p += getPosition();
128  vtk.points.push_back(p);
129  }
130  }
131 
132  //create connectivity matrix
133  //vtk.triangleStrips.reserve(vtk.triangleStrips.size()+nu);
134  for (unsigned i=0; i<nu; ++i) {
135  std::vector<double> cell;
136  cell.reserve(2*nv);
137  for (unsigned j=0; j<nv; ++j) {
138  cell.push_back(nPoints + j + i * (nv+1));
139  cell.push_back(nPoints + j + (i+1) * (nv+1));
140  }
141  vtk.triangleStrips.push_back(cell);
142  }
143 }
bool getDistance(Vec3D P, double radius, double &distance, Vec3D &normal) const
Find projection onto surface, return distance (and contactPoint)
std::enable_if<!std::is_pointer< typename U::MixedSpeciesType >::value, typename U::MixedSpeciesType * >::type getMixedObject(const U *S, const U *T)
const Vec3D & getPosition() const
Returns the position of this BaseInteractable.
BaseSpecies is the class from which all other species are derived.
Definition: BaseSpecies.h:49
This function defines a wall via a NurbsSurface.
Definition: NurbsWall.h:35
void rotateBack(Vec3D &position) const
Applies the inverse rotation to a position.
Definition: Quaternion.cc:592
Logger< MERCURY_LOGLEVEL > logger("MercuryKernel")
double Mdouble
Definition: GeneralDefine.h:34
const ParticleSpecies * getSpecies() const
Returns a pointer to the species of this BaseInteractable.
const std::complex< Mdouble > i
Definition: ExtendedMath.h:50
NurbsWall()
Default constructor: make a wall with default parameters.
Definition: NurbsWall.cc:39
std::string getName() const final
Returns the name of the object, here the string "Screw".
Definition: NurbsWall.cc:111
void write(std::ostream &os) const override
Function that writes a BaseWall to an output stream, usually a restart file.
Definition: BaseWall.cc:102
void rotate(Vec3D &position) const
Applies the rotation to a position.
Definition: Quaternion.cc:563
Mdouble getInteractionDistance() const
returns the largest separation distance at which adhesive short-range forces can occur.
Definition: BaseSpecies.h:146
void write(std::ostream &os) const override
Writes this wall to an output stream, for example a restart file.
Definition: NurbsWall.cc:102
WallHandler * getHandler() const
A function which returns the WallHandler that handles this BaseWall.
Definition: BaseWall.cc:136
std::vector< std::vector< double > > triangleStrips
Definition: BaseWall.h:39
Mdouble getRadius() const
Returns the particle's radius.
Definition: BaseParticle.h:345
Basic class for walls.
Definition: BaseWall.h:47
SpeciesHandler speciesHandler
A handler to that stores the species type i.e. LinearViscoelasticSpecies, etc.
Definition: DPMBase.h:1319
void read(std::istream &is) override
Function that reads a BaseWall from an input stream, usually a restart file.
Definition: BaseWall.cc:80
~NurbsWall()
Default destructor.
Definition: NurbsWall.cc:59
NurbsWall * copy() const final
Copy this wall and return a pointer to the copy.
Definition: NurbsWall.cc:64
std::vector< Vec3D > points
Definition: BaseWall.h:38
void writeVTK(VTKContainer &vtk) const override
Definition: NurbsWall.cc:116
Vec3D evaluate(double u, double v) const
Evaluate point on a nonrational NURBS surface.
void read(std::istream &is) override
Reads this wall from an input stream, for example a restart file.
Definition: NurbsWall.cc:91
bool getDistanceAndNormal(const BaseParticle &P, Mdouble &distance, Vec3D &normal_return) const final
Compute the distance from the Screw for a given BaseParticle and return if there is a collision...
Definition: NurbsWall.cc:74
Definition: Vector.h:49
T square(const T val)
squares a number
Definition: ExtendedMath.h:104
DPMBase * getDPMBase()
Gets the problem that is solved using this handler.
Definition: BaseHandler.h:725
const Quaternion & getOrientation() const
Returns the orientation of this BaseInteractable.
void set(const NurbsSurface &nurbsSurface)
Defines a standard wall, given an outward normal vector s.t. normal*x=normal*point for all x of the w...
Definition: NurbsWall.cc:69
NurbsSurface nurbsSurface_
The centre of the lower end of the screw.
Definition: NurbsWall.h:95