Homec4science
Diffusion akantu (bugfixes/65-Making-the-integrators-available-in-Python) Browse

rAKA/src/model/contact_mechanics/resolution.hh
271e73ac5baabugfixes/65-Making-the…

Last Change
Raw File

resolution.hh
View Options

/**
* @file resolution.hh
*
* @author Mohit Pundir <mohit.pundir@epfl.ch>
*
* @date creation: Mon Jan 7 2019
* @date last modification: Mon Jan 7 2019
*
* @brief Mother class for all contact resolutions
*
* @section LICENSE
*
* Copyright (©) 2010-2018 EPFL (Ecole Polytechnique Fédérale de Lausanne)
* Laboratory (LSMS - Laboratoire de Simulation en Mécanique des Solides)
*
* Akantu is free software: you can redistribute it and/or modify it under the
* terms of the GNU Lesser General Public License as published by the Free
* Software Foundation, either version 3 of the License, or (at your option) any
* later version.
*
* Akantu is distributed in the hope that it will be useful, but WITHOUT ANY
* WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
* A PARTICULAR PURPOSE. See the GNU Lesser General Public License for more
* details.
*
* You should have received a copy of the GNU Lesser General Public License
* along with Akantu. If not, see <http://www.gnu.org/licenses/>.
*
*/
/* -------------------------------------------------------------------------- */
#include "aka_factory.hh"
#include "aka_memory.hh"
#include "parsable.hh"
#include "parser.hh"
#include "fe_engine.hh"
#include "contact_element.hh"
#include "resolution_utils.hh"
#include "geometry_utils.hh"
/* -------------------------------------------------------------------------- */
#ifndef __AKANTU_RESOLUTION_HH__
#define __AKANTU_RESOLUTION_HH__
/* -------------------------------------------------------------------------- */
namespace akantu {
class Model;
class ContactMechanicsModel;
} // namespace akantu
namespace akantu {
/**
* Interface of all contact resolutions
* Prerequisites for a new resolution
* - inherit from this class
* - implement the following methods:
* \code
*
* virtual void computeNormalForce();
* virtual void computeFrictionForce();
*
* \endcode
*
*/
class Resolution : public Memory,
public Parsable {
/* ------------------------------------------------------------------------ */
/* Constructor/Destructor */
/* ------------------------------------------------------------------------ */
public:
/// instantiate contact resolution with defaults
Resolution(ContactMechanicsModel & model, const ID & id = "");
/// Destructor
~Resolution() override;
protected:
void initialize();
/// computes coordinates of a given element
void computeCoordinates(const Element & , Matrix<Real> &);
/* ------------------------------------------------------------------------ */
/* Functions that resolutions should reimplement for force */
/* ------------------------------------------------------------------------ */
public:
/// computes the force vector due to normal traction
virtual void computeNormalForce(__attribute__((unused)) const ContactElement &,
__attribute__((unused)) Vector<Real> &) {
AKANTU_TO_IMPLEMENT();
}
/// computes the tangential force vector due to frictional traction
virtual void computeTangentialForce(__attribute__((unused)) const ContactElement &,
__attribute__((unused)) Vector<Real> &) {
AKANTU_TO_IMPLEMENT();
}
protected:
/// local computation of trial tangential traction due to friction
virtual void computeTrialTangentialTraction(__attribute__((unused)) const ContactElement &,
__attribute__((unused)) const Matrix<Real> &,
__attribute__((unused)) Vector<Real> &) {
AKANTU_TO_IMPLEMENT();
}
/// local computation of tangential traction due to stick
virtual void computeStickTangentialTraction(__attribute__((unused)) const ContactElement &,
__attribute__((unused)) Vector<Real> &,
__attribute__((unused)) Vector<Real> &) {
AKANTU_TO_IMPLEMENT();
}
/// local computation of tangential traction due to slip
virtual void computeSlipTangentialTraction(__attribute__((unused)) const ContactElement &,
__attribute__((unused)) const Matrix<Real> &,
__attribute__((unused)) Vector<Real> &,
__attribute__((unused)) Vector<Real> &) {
AKANTU_TO_IMPLEMENT();
}
/* ------------------------------------------------------------------------ */
/* Functions that resolutions should reimplement for stiffness */
/* ------------------------------------------------------------------------ */
public:
/// compute the normal moduli due to normal traction
/*virtual void computeNormalModuli(__attribute__((unused)) const ContactElement &,
__attribute__((unused)) Matrix<Real> &,
__attribute__((unused)) Vector<Real> &,
__attribute__((unused)) Matrix<Real> & ) {
AKANTU_TO_IMPLEMENT();
}*/
/// compute the normal moduli due to normal traction
virtual void computeNormalModuli(__attribute__((unused)) const ContactElement &,
__attribute__((unused)) Matrix<Real> & ) {
AKANTU_TO_IMPLEMENT();
}
/// compute the tangent moduli due to tangential traction
virtual void computeTangentialModuli(__attribute__((unused)) const ContactElement &,
__attribute__((unused)) Matrix<Real> &,
__attribute__((unused)) Vector<Real> &,
__attribute__((unused)) Matrix<Real> & ) {
AKANTU_TO_IMPLEMENT();
}
/// compute the tangent moduli due to frictional traction
virtual void computeFrictionalModuli(__attribute__((unused)) Matrix<Real> &,
__attribute__((unused)) Array<Real> &,
__attribute__((unused)) Array<Real> &,
__attribute__((unused)) Array<Real> &,
__attribute__((unused)) Array<Real> &,
__attribute__((unused)) Matrix<Real> &,
__attribute__((unused)) Vector<Real> &,
__attribute__((unused)) ContactElement &) {
AKANTU_TO_IMPLEMENT();
}
/* ------------------------------------------------------------------------ */
/* Methods */
/* ------------------------------------------------------------------------ */
public:
/// assemble the residual for this resolution
void assembleInternalForces(GhostType ghost_type);
/// assemble the stiffness matrix for this resolution
void assembleStiffnessMatrix(GhostType ghost_type);
private:
/// assemble the residual for this resolution
void assembleInternalForces();
/// assemble the local array to global array for a contact element
void assembleLocalToGlobalArray(const ContactElement & , Vector<Real> & , Array<Real> & );
/// assemble the local stiffness to global stiffness for a contact element
void assembleLocalToGlobalMatrix(const ContactElement &, const Matrix<Real> &, SparseMatrix &);
public:
virtual void beforeSolveStep();
virtual void afterSolveStep(bool converged = true);
public:
/// function to print the contain of the class
void printself(std::ostream & stream, int indent = 0) const override;
/* ------------------------------------------------------------------------ */
/* Class Members */
/* ------------------------------------------------------------------------ */
protected:
/// friction coefficient : mu
Real mu;
/// spatial dimension
UInt spatial_dimension;
/// is master surface deformable
bool is_master_deformable;
/// Link to the fem object in the model
FEEngine & fem;
/// resolution name
std::string name;
/// model to which the resolution belong
ContactMechanicsModel & model;
};
/// standard output stream operator
inline std::ostream & operator<<(std::ostream & stream,
const Resolution & _this) {
_this.printself(stream);
return stream;
}
} // namespace akantu
/* -------------------------------------------------------------------------- */
namespace akantu {
using ResolutionFactory =
Factory<Resolution, ID, UInt, const ID &, ContactMechanicsModel &, const ID &>;
/// macaulay bracket to convert positive gap to zero
template <typename T>
T macaulay(T var) {return var < 0 ? 0 : var; }
template <typename T>
T heaviside(T var) {return var < 0 ? 0 : 1.0; }
} // namespace akantu
#define INSTANTIATE_RESOLUTION_ONLY(res_name) \
class res_name
#define RESOLUTION_DEFAULT_PER_DIM_ALLOCATOR(id, res_name) \
[](UInt dim, const ID &, ContactMechanicsModel & model, \
const ID & id) -> std::unique_ptr<Resolution> { \
switch (dim) { \
case 1: \
return std::make_unique<res_name>(model, id); \
case 2: \
return std::make_unique<res_name>(model, id); \
case 3: \
return std::make_unique<res_name>(model, id); \
default: \
AKANTU_EXCEPTION("The dimension " \
<< dim << "is not a valid dimension for the contact resolution " \
<< #id); \
} \
}
#define INSTANTIATE_RESOLUTION(id, res_name) \
INSTANTIATE_RESOLUTION_ONLY(res_name); \
static bool resolution_is_alocated_##id[[gnu::unused]] = \
ResolutionFactory::getInstance().registerAllocator( \
#id, RESOLUTION_DEFAULT_PER_DIM_ALLOCATOR(id, res_name))
#endif /* __AKANTU_RESOLUTION_HH__ */
c4science · Help