Homec4science
Diffusion akantu (features/updating_coverage_config) Browse

rAKA/src/model/solid_mechanics/materials/material_thermal.hh
461babc7c047features/updating_cover…

Last Change
Raw File

material_thermal.hh
View Options

/**
* Copyright (©) 2010-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 "material.hh"
/* -------------------------------------------------------------------------- */
#ifndef AKANTU_MATERIAL_THERMAL_HH_
#define AKANTU_MATERIAL_THERMAL_HH_
namespace akantu {
template <Int dim> class MaterialThermal : public Material {
/* ------------------------------------------------------------------------ */
/* Constructors/Destructors */
/* ------------------------------------------------------------------------ */
public:
MaterialThermal(SolidMechanicsModel & model, const ID & id = "",
const ID & fe_engine_id = "");
/* ------------------------------------------------------------------------ */
/* Methods */
/* ------------------------------------------------------------------------ */
public:
/// constitutive law for all element of a type
void computeStress(ElementType el_type, GhostType ghost_type) override;
/// local computation of thermal stress
template <class Args> inline void computeStressOnQuad(Args && args);
/* ------------------------------------------------------------------------ */
template <Int dim_ = dim>
decltype(auto) getArguments(ElementType el_type, GhostType ghost_type) {
return zip_append(Material::getArguments<dim_>(el_type, ghost_type),
"delta_T"_n = delta_T(el_type, ghost_type),
"sigma_th"_n = sigma_th(el_type, ghost_type),
"previous_sigma_th"_n =
sigma_th.previous(el_type, ghost_type));
}
template <Int dim_ = dim>
decltype(auto) getArgumentsTangent(Array<Real> & tangent_matrices,
ElementType el_type,
GhostType ghost_type) {
return Material::getArgumentsTangent<dim_>(tangent_matrices, el_type,
ghost_type);
}
/* ------------------------------------------------------------------------ */
/* Class Members */
/* ------------------------------------------------------------------------ */
protected:
/// Young modulus
Real E{0.};
/// Poisson ratio
Real nu{1. / 2.};
/// Thermal expansion coefficient
/// TODO : implement alpha as a matrix
Real alpha{0.};
/// Temperature field
InternalField<Real> & delta_T;
/// Current thermal stress
InternalField<Real> & sigma_th;
};
/* ------------------------------------------------------------------------ */
/* Inline impl */
/* ------------------------------------------------------------------------ */
template <Int dim>
template <class Args>
inline void MaterialThermal<dim>::computeStressOnQuad(Args && args) {
auto && sigma = args["sigma_th"_n];
auto && deltaT = args["delta_T"_n];
sigma = -this->E / (1. - 2. * this->nu) * this->alpha * deltaT;
}
template <>
template <class Args>
inline void MaterialThermal<1>::computeStressOnQuad(Args && args) {
auto && sigma = args["sigma_th"_n];
auto && deltaT = args["delta_T"_n];
sigma = -this->E * this->alpha * deltaT;
}
} // namespace akantu
#endif /* AKANTU_MATERIAL_THERMAL_HH_ */
c4science · Help