embedded_interface_model.cc
No OneTemporary
Actions

Subscribers

None

File Metadata

Created: Sun, Nov 10, 09:09

embedded_interface_model.cc
View Options

	/**
	* @file embedded_interface_model.cc
	*
	* @author Lucas Frérot <lucas.frerot@epfl.ch>
	*
	* @date creation: Mon Mar 9 2015
	* @date last modification: Mon Mar 9 2015
	*
	* @brief Model of Solid Mechanics with embedded interfaces
	*
	* @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/>.
	*
	*/

	/* -------------------------------------------------------------------------- */

	#include "aka_common.hh"

	#include "embedded_interface_model.hh"
	#include "material_reinforcement.hh"

	#ifdef AKANTU_USE_IOHELPER
	# include "dumper_paraview.hh"
	# include "dumpable_inline_impl.hh"
	#endif

	/* -------------------------------------------------------------------------- */

	__BEGIN_AKANTU__

	const EmbeddedInterfaceModelOptions
	default_embedded_interface_model_options(_explicit_lumped_mass, false, false);

	/* -------------------------------------------------------------------------- */
	EmbeddedInterfaceModel::EmbeddedInterfaceModel(Mesh & mesh,
	Mesh & primitive_mesh,
	UInt spatial_dimension,
	const ID & id,
	const MemoryID & memory_id) :
	SolidMechanicsModel(mesh, spatial_dimension, id, memory_id),
	intersector(mesh, primitive_mesh),
	interface_mesh(NULL),
	primitive_mesh(primitive_mesh),
	interface_material_selector(NULL)
	{
	// This pointer should be deleted by ~SolidMechanicsModel()
	MaterialSelector * mat_sel_pointer =
	new MeshDataMaterialSelector<std::string>("physical_names", *this);

	this->setMaterialSelector(*mat_sel_pointer);

	interface_mesh = &(intersector.getInterfaceMesh());

	// Create 1D FEEngine on the interface mesh
	registerFEEngineObject<MyFEEngineType>("EmbeddedInterfaceFEEngine", *interface_mesh, 1);
	}

	/* -------------------------------------------------------------------------- */
	EmbeddedInterfaceModel::~EmbeddedInterfaceModel() {
	delete interface_material_selector;
	}

	/* -------------------------------------------------------------------------- */
	void EmbeddedInterfaceModel::initFull(const ModelOptions & options) {
	const EmbeddedInterfaceModelOptions & eim_options =
	dynamic_cast<const EmbeddedInterfaceModelOptions &>(options);

	// We don't want to initiate materials before shape functions are initialized
	SolidMechanicsModelOptions dummy_options(eim_options.analysis_method, true);

	// Do no initialize interface_mesh if told so
	if (!eim_options.no_init_intersections)
	intersector.constructData();

	// Initialize interface FEEngine
	FEEngine & engine = getFEEngine("EmbeddedInterfaceFEEngine");
	engine.initShapeFunctions(_not_ghost);
	engine.initShapeFunctions(_ghost);

	SolidMechanicsModel::initFull(dummy_options);

	// This will call SolidMechanicsModel::initMaterials() last
	this->initMaterials();

	#if defined(AKANTU_USE_IOHELPER)
	this->mesh.registerDumper<DumperParaview>("reinforcement", id);
	this->mesh.addDumpMeshToDumper("reinforcement", *interface_mesh,
	1, _not_ghost, _ek_regular);
	#endif
	}

	/* -------------------------------------------------------------------------- */
	// This function is very similar to SolidMechanicsModel's
	void EmbeddedInterfaceModel::initMaterials() {
	Element element;

	delete interface_material_selector;
	interface_material_selector =
	new InterfaceMeshDataMaterialSelector<std::string>("physical_names", *this);

	for (ghost_type_t::iterator gt = ghost_type_t::begin(); gt != ghost_type_t::end(); ++gt) {
	element.ghost_type = *gt;

	Mesh::type_iterator it = interface_mesh->firstType(1, *gt);
	Mesh::type_iterator end = interface_mesh->lastType(1, *gt);

	for (; it != end ; ++it) {
	UInt nb_element = interface_mesh->getNbElement(it, gt);

	element.type = *it;

	Array<UInt> & mat_indexes = material_index.alloc(nb_element, 1, it, gt);

	for (UInt el = 0 ; el < nb_element ; el++) {
	element.element = el;
	UInt mat_index = (*interface_material_selector)(element);

	AKANTU_DEBUG_ASSERT(mat_index < materials.size(),
	"The material selector returned an index that does not exist");
	mat_indexes(element.element) = mat_index;
	materials.at(mat_index)->addElement(it, el, gt);
	}
	}
	}

	SolidMechanicsModel::initMaterials();
	}

	/**
	* DO NOT REMOVE - This prevents the material reinforcement to register
	* their number of components. Problems arise with AvgHomogenizingFunctor
	* if the material reinforcement gives its number of components for a
	* field. Since AvgHomogenizingFunctor verifies that all the fields
	* have the same number of components, an exception is raised.
	*/
	ElementTypeMap<UInt> EmbeddedInterfaceModel::getInternalDataPerElem(const std::string & field_name,
	const ElementKind & kind) {
	if (!(this->isInternal(field_name,kind))) AKANTU_EXCEPTION("unknown internal " << field_name);

	for (UInt m = 0; m < materials.size() ; ++m) {
	if (materials[m]->isInternal(field_name, kind)) {
	Material * mat = NULL;

	switch(this->spatial_dimension) {
	case 1:
	mat = dynamic_cast<MaterialReinforcement<1> *>(materials[m]);
	break;

	case 2:
	mat = dynamic_cast<MaterialReinforcement<2> *>(materials[m]);
	break;

	case 3:
	mat = dynamic_cast<MaterialReinforcement<3> *>(materials[m]);
	break;
	}

	if (mat == NULL && field_name != "stress_embedded")
	return materials[m]->getInternalDataPerElem(field_name,kind);
	else if (mat != NULL && field_name == "stress_embedded")
	return mat->getInternalDataPerElem(field_name, kind, "EmbeddedInterfaceFEEngine");
	}
	}

	return ElementTypeMap<UInt>();
	}

	/* -------------------------------------------------------------------------- */
	void EmbeddedInterfaceModel::addDumpGroupFieldToDumper(const std::string & dumper_name,
	const std::string & field_id,
	const std::string & group_name,
	const ElementKind & element_kind,
	bool padding_flag) {
	#ifdef AKANTU_USE_IOHELPER
	dumper::Field * field = NULL;

	// If dumper is reinforcement, create a 1D elemental field
	if (dumper_name == "reinforcement")
	field = this->createElementalField(field_id, group_name, padding_flag, 1, element_kind);
	else {
	try {
	SolidMechanicsModel::addDumpGroupFieldToDumper(dumper_name, field_id, group_name, element_kind, padding_flag);
	} catch (...) {}
	}
	if (field) {
	DumperIOHelper & dumper = mesh.getGroupDumper(dumper_name,group_name);
	Model::addDumpGroupFieldToDumper(field_id,field,dumper);
	}

	#endif
	}

	__END_AKANTU__

embedded_interface_model.ccNo OneTemporaryActions

File Metadata

embedded_interface_model.ccView Options

Event Timeline

embedded_interface_model.cc
No OneTemporary
Actions

embedded_interface_model.cc
View Options