Homec4science
Diffusion akantu (master) Browse

rAKA/model/solid_mechanics/solid_mechanics_model_material.cc
f5eb8635db83master

Last Change
Raw File

solid_mechanics_model_material.cc
View Options

/**
* @file solid_mechanics_model_material.cc
* @author Guillaume ANCIAUX <guillaume.anciaux@epfl.ch>
* @date Thu Nov 25 10:48:53 2010
*
* @brief
*
* @section LICENSE
*
* Copyright (©) 2010-2011 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 "solid_mechanics_model.hh"
#include "material.hh"
#include "aka_math.hh"
/* -------------------------------------------------------------------------- */
__BEGIN_AKANTU__
/* -------------------------------------------------------------------------- */
void SolidMechanicsModel::readMaterials(const std::string & filename) {
MaterialParser parser;
parser.open(filename);
std::string mat_type = parser.getNextMaterialType();
UInt mat_count = 0;
while (mat_type != ""){
std::stringstream sstr_mat; sstr_mat << id << ":" << mat_count++ << ":" << mat_type;
Material * material;
MaterialID mat_id = sstr_mat.str();
/// read the material properties
if(mat_type == "elastic") material = parser.readMaterialObject<MaterialElastic>(*this,mat_id);
else if(mat_type == "damage") material = parser.readMaterialObject<MaterialDamage>(*this,mat_id);
else if(mat_type == "mazars") material = parser.readMaterialObject<MaterialMazars>(*this,mat_id);
else AKANTU_DEBUG_ERROR("Malformed material file : unknown material type "
<< mat_type);
materials.push_back(material);
mat_type = parser.getNextMaterialType();
}
}
/* -------------------------------------------------------------------------- */
void SolidMechanicsModel::initMaterials() {
AKANTU_DEBUG_ASSERT(materials.size() != 0, "No material to initialize !");
Material ** mat_val = &(materials.at(0));
/// fill the element filters of the materials using the element_material arrays
const Mesh::ConnectivityTypeList & type_list = getFEM().getMesh().getConnectivityTypeList(_not_ghost);
Mesh::ConnectivityTypeList::const_iterator it;
for(it = type_list.begin(); it != type_list.end(); ++it) {
if(Mesh::getSpatialDimension(*it) != spatial_dimension) continue;
UInt nb_element = getFEM().getMesh().getNbElement(*it);
UInt * elem_mat_val = element_material[*it]->values;
for (UInt el = 0; el < nb_element; ++el) {
mat_val[elem_mat_val[el]]->addElement(*it, el);
}
}
/// @todo synchronize element material
/// fill the element filters of the materials using the element_material arrays
const Mesh::ConnectivityTypeList & ghost_type_list =
getFEM().getMesh().getConnectivityTypeList(_ghost);
for(it = ghost_type_list.begin(); it != ghost_type_list.end(); ++it) {
if(Mesh::getSpatialDimension(*it) != spatial_dimension) continue;
UInt nb_element = getFEM().getMesh().getNbGhostElement(*it);
UInt * elem_mat_val = ghost_element_material[*it]->values;
for (UInt el = 0; el < nb_element; ++el) {
mat_val[elem_mat_val[el]]->addGhostElement(*it, el);
}
}
std::vector<Material *>::iterator mat_it;
for(mat_it = materials.begin(); mat_it != materials.end(); ++mat_it) {
/// init internals properties
(*mat_it)->initMaterial();
}
}
__END_AKANTU__
c4science · Help