Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F81102824
igfem_enrichment.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
Wed, Sep 4, 05:28
Size
3 KB
Mime Type
text/x-c
Expires
Fri, Sep 6, 05:28 (2 d)
Engine
blob
Format
Raw Data
Handle
20489897
Attached To
rAKA akantu
igfem_enrichment.cc
View Options
/**
* @file igfem_enrichment.hh
*
* @author Aurelia Isabel Cuba Ramos <aurelia.cubaramos@epfl.ch>
*
*
* @brief Implementation of IGFEM enrichment
*
*
* Copyright (©) 2010-2012, 2014 EPFL (Ecole Polytechnique Fédérale de Lausanne)
* Laboratory (LSMS - Laboratoire de Simulation en Mécanique des Solides)
*
*/
/* -------------------------------------------------------------------------- */
#include "igfem_enrichment.hh"
namespace
akantu
{
/* -------------------------------------------------------------------------- */
IGFEMEnrichment
::
IGFEMEnrichment
(
Mesh
&
mesh
)
:
intersector_sphere
(
mesh
)
{}
/* -------------------------------------------------------------------------- */
void
IGFEMEnrichment
::
initialize
()
{
intersector_sphere
.
init
();
}
/* -------------------------------------------------------------------------- */
void
IGFEMEnrichment
::
update
(
ID
domain
)
{
if
(
domain
==
""
)
domain
=
default_geometry
;
Geometry
&
geometry
=
getGeometry
(
domain
);
intersector_sphere
.
buildIGFEMMeshFromSpheres
(
geometry
);
}
/* -------------------------------------------------------------------------- */
void
IGFEMEnrichment
::
unRegisterGeometryObject
(
const
ID
&
domain
)
{
GeometryMap
::
iterator
it
=
geometries
.
find
(
domain
);
AKANTU_DEBUG_ASSERT
(
it
!=
geometries
.
end
(),
"Geometry object with domain "
<<
domain
<<
" was not found"
);
geometries
.
erase
(
it
);
if
(
!
geometries
.
empty
())
default_geometry
=
(
*
geometries
.
begin
()).
first
;
}
/* -------------------------------------------------------------------------- */
void
IGFEMEnrichment
::
registerGeometryObject
(
Geometry
&
geometry
,
const
ID
&
domain
)
{
if
(
geometries
.
size
()
==
0
)
default_geometry
=
domain
;
#ifndef AKANTU_NDEBUG
GeometryMap
::
iterator
it
=
geometries
.
find
(
domain
);
AKANTU_DEBUG_ASSERT
(
it
==
geometries
.
end
(),
"Geometry object with domain "
<<
domain
<<
" was already created"
);
#endif
std
::
stringstream
sstr
;
sstr
<<
"geometry:"
<<
domain
;
geometries
[
domain
]
=
&
geometry
;
}
/* -------------------------------------------------------------------------- */
IGFEMEnrichment
::
Geometry
&
IGFEMEnrichment
::
getGeometry
(
ID
&
domain
)
const
{
AKANTU_DEBUG_IN
();
if
(
domain
==
""
)
domain
=
default_geometry
;
GeometryMap
::
const_iterator
it
=
geometries
.
find
(
domain
);
AKANTU_DEBUG_ASSERT
(
it
!=
geometries
.
end
(),
"The geometry "
<<
domain
<<
" is not registered"
);
AKANTU_DEBUG_OUT
();
return
*
(
it
->
second
);
}
/* -------------------------------------------------------------------------- */
void
IGFEMEnrichment
::
moveInterface
(
Real
new_position
,
ID
domain
)
{
if
(
domain
==
""
)
domain
=
default_geometry
;
Geometry
&
geometry
=
getGeometry
(
domain
);
/// for this type of IGFEM enrichment the geometry consists of a list of
/// spheres
/// -> need to loop over spheres and change their radius,
/// which specifies the position of interfaces
Geometry
::
const_iterator
query_it
=
geometry
.
begin
();
Geometry
sphere_list
;
for
(;
query_it
!=
geometry
.
end
();
++
query_it
)
{
SK
::
Sphere_3
sphere
(
query_it
->
center
(),
new_position
*
new_position
);
sphere_list
.
push_back
(
sphere
);
}
geometry
.
zero
();
geometry
=
sphere_list
;
this
->
update
(
domain
);
}
}
// namespace akantu
Event Timeline
Log In to Comment