Homec4science
Diffusion akantu (bugfixes/bug-fix-in-contact-and-interface-extension) Browse

rAKA/src/model/solid_mechanics/materials/material_embedded/material_reinforcement_template_tmpl.hh
ecf274c1ee81bugfixes/bug-fix-in…

Last Change
Raw File

material_reinforcement_template_tmpl.hh
View Options

/**
* @file material_reinforcement_template_tmpl.hh
*
* @author Lucas Frérot <lucas.frerot@epfl.ch>
*
* @date creation: Mon Mar 16 2015
* @date last modification: Mon Mar 16 2015
*
* @brief Reinforcement material templated with constitutive law
*
* @section LICENSE
*
* Copyright (©) 2010-2012, 2014 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/>.
*
*/
// /!\ no namespace here !
/* -------------------------------------------------------------------------- */
template<UInt dim, class ConstLaw>
MaterialReinforcementTemplate<dim, ConstLaw>::MaterialReinforcementTemplate(SolidMechanicsModel & model,
const ID & id):
Material(model, id),
MaterialReinforcement<dim>(model, id),
ConstLaw(model, id) {
this->element_filter.free();
this->model->getInterfaceMesh().initElementTypeMapArray(this->element_filter, 1, 1,
false, _ek_regular);
}
/* -------------------------------------------------------------------------- */
template<UInt dim, class ConstLaw>
MaterialReinforcementTemplate<dim, ConstLaw>::~MaterialReinforcementTemplate()
{}
/* -------------------------------------------------------------------------- */
/// TODO this function is super ugly. Find better way to initialize internals
template<UInt dim, class ConstLaw>
void MaterialReinforcementTemplate<dim, ConstLaw>::initMaterial() {
MaterialReinforcement<dim>::initMaterial();
// Gotta change the dimension from here onwards
this->ConstLaw::spatial_dimension = 1;
ConstLaw::initMaterial();
this->ConstLaw::gradu.free();
this->ConstLaw::stress.free();
this->ConstLaw::delta_T.free();
this->ConstLaw::sigma_th.free();
this->ConstLaw::potential_energy.free();
Mesh::type_iterator
type_it = this->MaterialReinforcement<dim>::model->getInterfaceMesh().firstType(),
type_end = this->MaterialReinforcement<dim>::model->getInterfaceMesh().lastType();
// Reshape the ConstLaw internal fields
for (; type_it != type_end ; ++type_it) {
UInt nb_elements = this->MaterialReinforcement<dim>::element_filter(*type_it).getSize();
FEEngine & interface_engine =
this->MaterialReinforcement<dim>::model->getFEEngine("EmbeddedInterfaceFEEngine");
UInt nb_quad = interface_engine.getNbQuadraturePoints(*type_it);
this->ConstLaw::gradu.alloc(nb_elements * nb_quad, 1, *type_it, _not_ghost);
this->ConstLaw::stress.alloc(nb_elements * nb_quad, 1, *type_it, _not_ghost);
this->ConstLaw::stress.initialize(1);
this->ConstLaw::delta_T.alloc(nb_elements * nb_quad, 1, *type_it, _not_ghost);
this->ConstLaw::sigma_th.alloc(nb_elements * nb_quad, 1, *type_it, _not_ghost);
this->ConstLaw::potential_energy.alloc(nb_elements * nb_quad, 1, *type_it, _not_ghost);
if (this->ConstLaw::use_previous_stress) {
this->ConstLaw::stress.initializeHistory();
this->ConstLaw::stress.previous().free();
this->ConstLaw::stress.previous().alloc(nb_elements * nb_quad, 1, *type_it, _not_ghost);
}
if (this->ConstLaw::use_previous_gradu) {
this->ConstLaw::gradu.initializeHistory();
this->ConstLaw::gradu.previous().free();
this->ConstLaw::gradu.previous().alloc(nb_elements * nb_quad, 1, *type_it, _not_ghost);
}
if (this->ConstLaw::use_previous_stress_thermal) {
this->ConstLaw::sigma_th.initializeHistory();
this->ConstLaw::sigma_th.previous().free();
this->ConstLaw::sigma_th.previous().alloc(nb_elements * nb_quad, 1, *type_it, _not_ghost);
}
this->initPlastic(*this);
}
}
/* -------------------------------------------------------------------------- */
template<UInt dim, class ConstLaw>
void MaterialReinforcementTemplate<dim, ConstLaw>::computeGradU(const ElementType & el_type,
GhostType ghost_type) {
AKANTU_DEBUG_IN();
MaterialReinforcement<dim>::computeGradU(el_type, ghost_type);
const UInt voigt_size = Material::getTangentStiffnessVoigtSize(dim);
Array<Real>::matrix_iterator full_gradu_it =
this->MaterialReinforcement<dim>::gradu(el_type, ghost_type).begin(dim, dim);
Array<Real>::matrix_iterator full_gradu_end =
this->MaterialReinforcement<dim>::gradu(el_type, ghost_type).end(dim, dim);
Array<Real>::scalar_iterator gradu_it =
this->ConstLaw::gradu(el_type, ghost_type).begin();
Array<Real>::matrix_iterator cosines_it =
this->directing_cosines(el_type, ghost_type).begin(voigt_size, voigt_size);
for (; full_gradu_it != full_gradu_end ; ++full_gradu_it,
++gradu_it,
++cosines_it) {
Matrix<Real> & full_gradu = *full_gradu_it;
Real & gradu = *gradu_it;
Matrix<Real> & C = *cosines_it;
computeInterfaceGradUOnQuad(full_gradu, gradu, C);
}
AKANTU_DEBUG_OUT();
}
/* -------------------------------------------------------------------------- */
template<UInt dim, class ConstLaw>
void MaterialReinforcementTemplate<dim, ConstLaw>::computeInterfaceGradUOnQuad(const Matrix<Real> & full_gradu,
Real & gradu,
const Matrix<Real> & C) {
const UInt voigt_size = Material::getTangentStiffnessVoigtSize(dim);
Matrix<Real> epsilon(dim, dim);
Vector<Real> e_voigt(voigt_size);
Vector<Real> e_interface_voigt(voigt_size);
epsilon = 0.5 * (full_gradu + full_gradu.transpose());
MaterialReinforcement<dim>::strainTensorToVoigtVector(epsilon, e_voigt);
e_interface_voigt.mul<false>(C, e_voigt);
gradu = e_interface_voigt(0);
}
/* -------------------------------------------------------------------------- */
template<UInt dim, class ConstLaw>
void MaterialReinforcementTemplate<dim, ConstLaw>::computeTangentModuli(const ElementType & el_type,
Array<Real> & tangent,
GhostType ghost_type) {
AKANTU_DEBUG_IN();
AKANTU_DEBUG_ASSERT(tangent.getNbComponent() == 1, "Reinforcements only work in 1D");
ConstLaw::computeTangentModuli(el_type, tangent, ghost_type);
AKANTU_DEBUG_OUT();
}
/* -------------------------------------------------------------------------- */
template<UInt dim, class ConstLaw>
void MaterialReinforcementTemplate<dim, ConstLaw>::computeStress(ElementType type,
GhostType ghost_type) {
AKANTU_DEBUG_IN();
ConstLaw::computeStress(type, ghost_type);
Array<Real>::matrix_iterator full_sigma_it =
this->MaterialReinforcement<dim>::stress(type, ghost_type).begin(dim, dim);
Array<Real>::matrix_iterator full_sigma_end =
this->MaterialReinforcement<dim>::stress(type, ghost_type).end(dim, dim);
Array<Real>::scalar_iterator sigma_it =
this->ConstLaw::stress(type, ghost_type).begin();
Array<Real>::scalar_iterator pre_stress_it =
this->MaterialReinforcement<dim>::pre_stress(type, ghost_type).begin();
for (; full_sigma_it != full_sigma_end ; ++full_sigma_it, ++sigma_it, ++pre_stress_it) {
Matrix<Real> & sigma = *full_sigma_it;
sigma(0, 0) = *sigma_it + *pre_stress_it;
}
AKANTU_DEBUG_OUT();
}
/* -------------------------------------------------------------------------- */
template<UInt dim, class ConstLaw>
Real MaterialReinforcementTemplate<dim, ConstLaw>::getEnergy(std::string id) {
return MaterialReinforcement<dim>::getEnergy(id);
}
/* -------------------------------------------------------------------------- */
template <UInt dim, class ConstLaw>
void MaterialReinforcementTemplate<dim, ConstLaw>::computePotentialEnergy(ElementType type,
GhostType ghost_type) {
ConstLaw::computePotentialEnergy(type, ghost_type);
}
/* -------------------------------------------------------------------------- */
template <UInt dim, class ConstLaw>
void MaterialReinforcementTemplate<dim, ConstLaw>::flattenInternal(const std::string & field_id,
ElementTypeMapArray<Real> & internal_flat,
const GhostType ghost_type,
ElementKind element_kind) {
if (field_id == "stress_embedded" || field_id == "inelastic_strain")
MaterialReinforcement<dim>::flattenInternal(field_id, internal_flat, ghost_type, element_kind);
}
template<UInt dim, class ConstLaw>
void MaterialReinforcementTemplate<dim, ConstLaw>::initPlastic(
MaterialReinforcementTemplate<dim, MaterialElastic<1> > & _this) {
}
template<UInt dim, class ConstLaw>
void MaterialReinforcementTemplate<dim, ConstLaw>::initPlastic(
MaterialReinforcementTemplate<dim, MaterialLinearIsotropicHardening<1> > & _this) {
_this.ConstLaw::iso_hardening.free();
_this.ConstLaw::plastic_energy.free();
_this.ConstLaw::d_plastic_energy.free();
_this.ConstLaw::inelastic_strain.free();
Mesh::type_iterator
type_it = _this.MaterialReinforcement<dim>::model->getInterfaceMesh().firstType(),
type_end = _this.MaterialReinforcement<dim>::model->getInterfaceMesh().lastType();
for (; type_it != type_end ; ++type_it) {
UInt nb_elements = _this.MaterialReinforcement<dim>::element_filter(*type_it).getSize();
FEEngine & interface_engine =
_this.MaterialReinforcement<dim>::model->getFEEngine("EmbeddedInterfaceFEEngine");
UInt nb_quad = interface_engine.getNbQuadraturePoints(*type_it);
_this.ConstLaw::iso_hardening.alloc(nb_elements * nb_quad, 1, *type_it, _not_ghost);
_this.ConstLaw::plastic_energy.alloc(nb_elements * nb_quad, 1, *type_it, _not_ghost);
_this.ConstLaw::d_plastic_energy.alloc(nb_elements * nb_quad, 1, *type_it, _not_ghost);
_this.ConstLaw::inelastic_strain.alloc(nb_elements * nb_quad, 1, *type_it, _not_ghost);
_this.plastic_energy.initialize(1);
_this.d_plastic_energy.initialize(1);
_this.iso_hardening.initialize(1);
_this.inelastic_strain.initialize(1);
_this.ConstLaw::iso_hardening.initializeHistory();
_this.ConstLaw::iso_hardening.previous().free();
_this.ConstLaw::iso_hardening.previous().alloc(nb_elements * nb_quad, 1, *type_it, _not_ghost);
_this.ConstLaw::plastic_energy.initializeHistory();
_this.ConstLaw::plastic_energy.previous().free();
_this.ConstLaw::plastic_energy.previous().alloc(nb_elements * nb_quad, 1, *type_it, _not_ghost);
_this.ConstLaw::d_plastic_energy.initializeHistory();
_this.ConstLaw::d_plastic_energy.previous().free();
_this.ConstLaw::d_plastic_energy.previous().alloc(nb_elements * nb_quad, 1, *type_it, _not_ghost);
_this.ConstLaw::inelastic_strain.initializeHistory();
_this.ConstLaw::inelastic_strain.previous().free();
_this.ConstLaw::inelastic_strain.previous().alloc(nb_elements * nb_quad, 1, *type_it, _not_ghost);
}
}
template<UInt dim, class ConstLaw>
void MaterialReinforcementTemplate<dim, ConstLaw>::savePreviousState() {
ConstLaw::savePreviousState();
}
c4science · Help