File Metadata

Created: Sun, Nov 3, 23:14

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

	#include <CGAL/Cartesian.h>

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

	__BEGIN_AKANTU__

	typedef CGAL::Cartesian<Real> K;

	EmbeddedInterfaceModel::EmbeddedInterfaceModel(Mesh & mesh,
	UInt spatial_dimension,
	const ID & id,
	const MemoryID & memory_id) :
	SolidMechanicsModel(mesh, spatial_dimension, id, memory_id),
	interface_mesh(NULL),
	interface_material_selector(NULL),
	interface_container(mesh)
	{}

	EmbeddedInterfaceModel::~EmbeddedInterfaceModel() {
	delete interface_material_selector;
	}

	void EmbeddedInterfaceModel::initModel() {
	interface_container.constructData();
	instanciateInterfaces();

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

	SolidMechanicsModel::initModel();
	}

	void EmbeddedInterfaceModel::initMaterials() {
	delete interface_material_selector;
	interface_material_selector = new EmbeddedInterfaceMaterialSelector<std::string>("material", *this);
	this->setMaterialSelector(*interface_material_selector);

	Element element;
	//Material ** mat_val = &(materials.at(0));

	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();
	}

	void EmbeddedInterfaceModel::instanciateInterfaces() {
	const std::pair<Parser::const_section_iterator, Parser::const_section_iterator> &
	sub_sect = this->parser->getSubSections(_st_embedded_interface);

	Parser::const_section_iterator section_it = sub_sect.first;

	std::list<std::pair<K::Segment_3, std::string> > interface_list;

	for (; section_it != sub_sect.second ; ++section_it) {
	EmbeddedInterface interface(this->spatial_dimension);
	interface.parseSection(*section_it);

	interface_list.push_back(std::make_pair(interface.getPrimitive(), interface.getMaterialName()));
	}

	initInterface(interface_list);
	}

	void EmbeddedInterfaceModel::initInterface(const std::list<std::pair<K::Segment_3, std::string> > & interface_list) {
	AKANTU_DEBUG_IN();

	interface_mesh = &(interface_container.meshOfLinearInterfaces(interface_list));

	registerFEEngineObject<MyFEEngineType>("EmbeddedInterfaceFEEngine", *interface_mesh, 1);
	FEEngine & engine = getFEEngine("EmbeddedInterfaceFEEngine");
	engine.initShapeFunctions(_not_ghost);
	engine.initShapeFunctions(_ghost);

	AKANTU_DEBUG_OUT();
	}

	void EmbeddedInterfaceModel::assembleStiffnessMatrix() {
	SolidMechanicsModel::assembleStiffnessMatrix();
	}

	dumper::Field * EmbeddedInterfaceModel::createElementalField(const std::string & field_name,
	const std::string & group_name,
	bool padding_flag,
	const ElementKind & kind,
	const std::string & fe_engine_id) {

	if (field_name == "stress_embedded") {
	return SolidMechanicsModel::createElementalField(field_name,
	group_name,
	padding_flag,
	kind,
	"EmbeddedInterfaceFEEngine");
	} else {
	return SolidMechanicsModel::createElementalField(field_name,
	group_name,
	padding_flag,
	kind,
	fe_engine_id);
	}
	}

	__END_AKANTU__

embedded_interface_model.cc
No OneTemporary
Actions

File Metadata

embedded_interface_model.cc
View Options

Event Timeline

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