Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F75802206
needleman.cc
No One
Temporary
Actions
Download File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Subscribers
None
File Metadata
Details
File Info
Storage
Attached
Created
Sun, Aug 4, 09:13
Size
6 KB
Mime Type
text/x-c++
Expires
Tue, Aug 6, 09:13 (2 d)
Engine
blob
Format
Raw Data
Handle
19610911
Attached To
rLIBMULTISCALE LibMultiScale
needleman.cc
View Options
/**
* @file needleman.cc
*
* @author Guillaume Anciaux <guillaume.anciaux@epfl.ch>
* @author Till Junge <till.junge@epfl.ch>
*
* @date Wed Jan 15 17:00:43 2014
*
* @brief Needleman style discrete dislocation coupler between DD and FEM
*
* @section LICENSE
*
* Copyright (©) 2010-2011 EPFL (Ecole Polytechnique Fédérale de Lausanne)
* Laboratory (LSMS - Laboratoire de Simulation en Mécanique des Solides)
*
* LibMultiScale 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.
*
* LibMultiScale 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 LibMultiScale. If not, see <http://www.gnu.org/licenses/>.
*
*/
#include "needleman.hh"
#include "lib_bridging.hh"
#include "lm_common.hh"
/* -------------------------------------------------------------------------- */
__BEGIN_LIBMULTISCALE__
Needleman::Needleman(const std::string &name)
: LMObject(name), CouplingDDContinuum(name),
bridging_nodes("bridging:" + name)
//,
// bridging_sources(this->domA.getContainerNodes(),this->domC.getContainer())
{}
/* -------------------------------------------------------------------------- */
Needleman::~Needleman() { clearAll(); }
/* -------------------------------------------------------------------------- */
template <typename DomainDD, typename DomainC>
void Needleman::init(DomainDD &domDD, DomainC &domC) {
// initialize the bridging_zone object
bridging_nodes.init(domDD, domC);
// //forward name
// sprintf(temp,"%s-bridging",Parent::name.c_str());
// bridging_sources.setName(temp);
// // initialize the bridging_zone object
// bridging_sources.init();
}
/* -------------------------------------------------------------------------- */
template <typename DomainDD, typename DomainC>
void Needleman::coupling(CouplingStage stage, DomainDD &domDD, DomainC &domC) {
if (stage != COUPLING_STEP3)
return;
// since dislocations move I must recompute their positions and associations
// with elements
recomputeBridgingInformationOnNodes();
LM_TOIMPLEMENT;
// applyStressesOnNodes(domDD, domC);
applyStressesOnSources();
// apply boundary conditions to the FEM zone from image stresses of the
// dislocation stress fields
applyBoundaryConditionToFEM(domC, domDD);
}
/* -------------------------------------------------------------------------- */
void Needleman::recomputeBridgingInformationOnNodes() {
LM_FATAL("to adapt");
// construct sublist of coupled nodes
// bridging_nodes.buildPointList();
// // build positions of atom vector
// bridging_nodes.buildPositions();
// // build shape matrix
// bridging_nodes.buildShapeMatrix(bridging_nodes->pointList->nbElem());
// //filter unmatched atomic DOFs
// bridging_nodes.filterPointListForUnmatched();
// //filter assoc vector
// bridging_nodes.filterAssoc();
}
/* -------------------------------------------------------------------------- */
template <typename ContC>
void Needleman::applyStressesOnNodes(ContC &meshList) {
LM_TOIMPLEMENT;
// assign the stress to nodes
// auto &elems = meshList.getContainerElems();
// for (auto &&[nd, assoc] :
// zip(bridging_nodes.pointList, bridging_nodes.getAssoc())) {
// // assoc return the element index in which point i resides
// // elementList is a container of elements
// auto &&el = elems.get(assoc);
// // conversion from 2D to 3D (temporary ?)
// nd.stress(0) = el.stress(0);
// nd.stress(1) = el.stress(1);
// nd.stress(2) = 0;
// nd.stress(3) = el.stress(2);
// nd.stress(4) = 0;
// nd.stress(5) = 0;
// }
}
/* -------------------------------------------------------------------------- */
void Needleman::applyStressesOnSources() {
// // assign the stress to sources
// IteratorDDSourcesSubset itsources =
// bridging_sources->pointList->getIterator(); RefDDSource nd =
// itsources.getFirst(); UInt i = 0; RefElem el; while (!itsources.end()){
// // assoc return the element index in which poUInt i resides
// // elementList is a conatiner of elements
// el = bridging_sources->elementList.get(bridging_sources->getAssoc(i));
// // conversion from 2D to 3D (temporary ?)
// nd.stress(0) = el.stress(0);
// nd.stress(1) = el.stress(1);
// nd.stress(2) = 0;
// nd.stress(3) = el.stress(2);
// nd.stress(4) = 0;
// nd.stress(5) = 0;
// ++i;
// nd = itsources.getNext();
// }
}
/* -------------------------------------------------------------------------- */
template <typename DomC, typename DomA>
void Needleman::applyBoundaryConditionToFEM(DomC &domC,
DomA &domA [[gnu::unused]]) {
auto &contElem = domC.getContainerElems();
LM_FATAL("what was that ?");
// contElem.boundarycontainerSwitch(1);
auto contNodes = domC.getContainerNodes();
for (auto &&el [[gnu::unused]] : contElem) {
std::vector<UInt> nodes;
LM_TOIMPLEMENT;
// el.globalIndexes(nodes);
// Real nx, ny, nz;
// el.normal(nx, ny, nz);
// Real normals[3];
// normals[0] = nx;
// normals[1] = ny;
// normals[2] = 0;
UInt nnodes = nodes.size();
Real quadrature_point[2] = {0., 0.};
LM_TOIMPLEMENT;
// Real h = el.length();
for (UInt n = 0; n < nnodes; ++n) {
RefNode nd = contNodes.get(n);
LM_TOIMPLEMENT;
// quadrature_point[0] += nd.position0(0);
// quadrature_point[1] += nd.position0(1);
}
quadrature_point[0] /= 2;
quadrature_point[1] /= 2;
LM_TOIMPLEMENT;
// Real tractions[3];
// domA.imageForce(quadrature_point, normals, tractions);
for (UInt n = 0; n < nnodes; ++n) {
RefNode nd = contNodes.get(n);
LM_TOIMPLEMENT;
// Real fx = nd.force()[0] * h;
// Real fy = nd.force()[1] * h;
LM_TOIMPLEMENT;
// nd.setForce(0, fx + 0.5 * tractions[0]);
// nd.setForce(0, fy + 0.5 * tractions[1]);
}
}
LM_FATAL("what was that");
// contElem.boundarycontainerSwitch(0);
}
/* -------------------------------------------------------------------------- */
/* LMDESC NEEDLEMAN
Needleman style discrete dislocation coupler between DD and FEM
*/
/* LMEXAMPLE
Will be added
*/
void Needleman::declareParams() {}
/* -------------------------------------------------------------------------- */
DECLARE_COUPLER_INIT_MAKE_CALL(Needleman, domDD, domCONTINUUM)
DECLARE_COUPLING_MAKE_CALL(Needleman, domDD, domCONTINUUM)
/* ------------------------------------------------------------------------ */
__END_LIBMULTISCALE__
Event Timeline
Log In to Comment