Homec4science
Diffusion akantu (bugfix/69-segmentation-fault-caused-by-model-applybc-aka-fromtraction-trac-traction) Browse

rAKA/src/model/phase_field/phasefield.hh
bf637b29f909bugfix/69-segmentation…

Last Change
Raw File

phasefield.hh
View Options

/**
* Copyright (©) 2020-2023 EPFL (Ecole Polytechnique Fédérale de Lausanne)
* Laboratory (LSMS - Laboratoire de Simulation en Mécanique des Solides)
*
* This file is part of Akantu
*
* 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 "constitutive_law.hh"
/* -------------------------------------------------------------------------- */
#ifndef AKANTU_PHASEFIELD_HH_
#define AKANTU_PHASEFIELD_HH_
/* -------------------------------------------------------------------------- */
namespace akantu {
class PhaseFieldModel;
class PhaseField;
} // namespace akantu
namespace akantu {
using PhaseFieldFactory =
Factory<PhaseField, ID, Int, const ID &, PhaseFieldModel &, const ID &>;
class PhaseField : public ConstitutiveLaw<PhaseFieldModel> {
using Parent = ConstitutiveLaw<PhaseFieldModel>;
/* ------------------------------------------------------------------------ */
/* Constructors/Destructors */
/* ------------------------------------------------------------------------ */
public:
PhaseField(PhaseFieldModel & model, const ID & id = "",
const ID & fe_engine_id = "");
/* ------------------------------------------------------------------------ */
/* Methods */
/* ------------------------------------------------------------------------ */
public:
/// initialize the phasefield computed parameter
virtual void initPhaseField() { Parent::initConstitutiveLaw(); }
void initConstitutiveLaw() override { this->initPhaseField(); }
///
virtual void beforeSolveStep();
/// assemble the residual for this phasefield
virtual void assembleInternalForces(GhostType ghost_type);
/// assemble the stiffness matrix for this phasefield
virtual void assembleStiffnessMatrix(GhostType ghost_type);
/// compute the driving force for this phasefield
virtual void computeAllDrivingForces(GhostType ghost_type = _not_ghost);
protected:
/// compute the dissipated energy by element
void computeDissipatedEnergyByElements();
/// function called to updatet the internal parameters when the
/// modifiable parameters are modified
void updateInternalParameters() override;
// constitutive law for driving force
virtual void computeDrivingForce(ElementType /* el_type */,
GhostType /* ghost_type */ = _not_ghost) {
AKANTU_TO_IMPLEMENT();
}
/// compute the dissiapted energy
virtual void computeDissipatedEnergy(ElementType el_type);
/// compute the dissipated energy for an element
virtual void
computeDissipatedEnergyByElement(const Element & /*element*/,
Vector<Real> & /*edis_on_quad_points*/) {
AKANTU_TO_IMPLEMENT();
}
/// compute the dissipated energy for an element
virtual void
computeDissipatedEnergyByElement(ElementType /*type*/, Idx /*index*/,
Vector<Real> & /*edis_on_quad_points*/) {
AKANTU_TO_IMPLEMENT();
}
/* ------------------------------------------------------------------------ */
/* Accessors */
/* ------------------------------------------------------------------------ */
public:
/// static method to reteive the material factory
static PhaseFieldFactory & getFactory();
/// return the damage energyfor the subset of elements contained
/// by the phasefield
[[nodiscard]] Real getEnergy(const ID & energy_id) override;
/// Compute dissipated energy for an individual element
[[nodiscard]] Real getEnergy(const ID & energy_id,
const Element & element) override;
AKANTU_GET_MACRO_BY_ELEMENT_TYPE_CONST(Strain, strain, Real);
AKANTU_GET_MACRO_AUTO(Strain, strain);
AKANTU_GET_MACRO_AUTO_NOT_CONST(Strain, strain);
AKANTU_GET_MACRO_BY_ELEMENT_TYPE_CONST(Damage, damage_on_qpoints, Real);
AKANTU_GET_MACRO_AUTO_NOT_CONST(Damage, damage_on_qpoints);
AKANTU_GET_MACRO_AUTO(Damage, damage_on_qpoints);
/* ------------------------------------------------------------------------ */
/* Class Members */
/* ------------------------------------------------------------------------ */
protected:
/// length scale parameter
Real l0{0.};
/// Young's modulus
Real E{0.};
/// Poisson ratio
Real nu{0.};
/// Isotropic formulation
bool isotropic{true};
/// Lame's first parameter
Real lambda{0.};
/// Lame's second paramter
Real mu{0.};
/// critical energy release rate
// Real g_c;
DefaultRandomInternalField<Real> & g_c;
/// damage arrays ordered by element types
InternalField<Real> & damage_on_qpoints;
/// grad_d arrays ordered by element types
InternalField<Real> & gradd;
/// phi arrays ordered by element types
InternalField<Real> & phi;
/// strain arrays ordered by element types
InternalField<Real> & strain;
/// driving force ordered by element types
InternalField<Real> & driving_force;
/// driving energy ordered by element types
InternalField<Real> & driving_energy;
/// damage energy ordered by element types
InternalField<Real> & damage_energy;
/// damage energy density ordered by element types
InternalField<Real> & damage_energy_density;
/// dissipated energy by element
InternalField<Real> & dissipated_energy;
};
} // namespace akantu
#include "phasefield_inline_impl.hh"
#include "internal_field_tmpl.hh"
#include "random_internal_field_tmpl.hh"
namespace akantu {
namespace {
template <template <Int> class PF> bool instantiatePhaseField(const ID & id) {
return PhaseFieldFactory::getInstance().registerAllocator(
id, [](Int dim, const ID &, PhaseFieldModel & model, const ID & id) {
return tuple_dispatch<AllSpatialDimensions>(
[&](auto && _) -> std::unique_ptr<PhaseField> {
constexpr auto && dim_ = aka::decay_v<decltype(_)>;
return std::make_unique<PF<dim_>>(model, id);
},
dim);
});
}
} // namespace
} // namespace akantu
#endif
c4science · Help