Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F79779151
solid_mechanics_model_cohesive.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, Aug 28, 04:59
Size
24 KB
Mime Type
text/x-c
Expires
Fri, Aug 30, 04:59 (2 d)
Engine
blob
Format
Raw Data
Handle
20225891
Attached To
rAKA akantu
solid_mechanics_model_cohesive.cc
View Options
/**
* @file solid_mechanics_model_cohesive.cc
*
* @author Mauro Corrado <mauro.corrado@epfl.ch>
* @author Nicolas Richart <nicolas.richart@epfl.ch>
* @author Marco Vocialta <marco.vocialta@epfl.ch>
*
* @date creation: Tue May 08 2012
* @date last modification: Wed Jan 13 2016
*
* @brief Solid mechanics model for cohesive elements
*
* @section LICENSE
*
* Copyright (©) 2010-2012, 2014, 2015 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_cohesive.hh"
#include "aka_iterators.hh"
#include "cohesive_element_inserter.hh"
#include "element_synchronizer.hh"
#include "facet_synchronizer.hh"
#include "fe_engine_template.hh"
#include "integrator_gauss.hh"
#include "material_cohesive.hh"
#include "parser.hh"
#include "shape_cohesive.hh"
#include "dumpable_inline_impl.hh"
#ifdef AKANTU_USE_IOHELPER
#include "dumper_iohelper_paraview.hh"
#endif
/* -------------------------------------------------------------------------- */
#include <algorithm>
/* -------------------------------------------------------------------------- */
namespace
akantu
{
const
SolidMechanicsModelCohesiveOptions
default_solid_mechanics_model_cohesive_options
(
_explicit_lumped_mass
,
false
);
/* -------------------------------------------------------------------------- */
SolidMechanicsModelCohesive
::
SolidMechanicsModelCohesive
(
Mesh
&
mesh
,
UInt
dim
,
const
ID
&
id
,
const
MemoryID
&
memory_id
)
:
SolidMechanicsModel
(
mesh
,
dim
,
id
,
memory_id
,
ModelType
::
_solid_mechanics_model_cohesive
),
tangents
(
"tangents"
,
id
),
facet_stress
(
"facet_stress"
,
id
),
facet_material
(
"facet_material"
,
id
)
{
AKANTU_DEBUG_IN
();
auto
&&
tmp_material_selector
=
std
::
make_shared
<
DefaultMaterialCohesiveSelector
>
(
*
this
);
tmp_material_selector
->
setFallback
(
this
->
material_selector
);
this
->
material_selector
=
tmp_material_selector
;
#if defined(AKANTU_USE_IOHELPER)
this
->
mesh
.
registerDumper
<
DumperParaview
>
(
"cohesive elements"
,
id
);
this
->
mesh
.
addDumpMeshToDumper
(
"cohesive elements"
,
mesh
,
Model
::
spatial_dimension
,
_not_ghost
,
_ek_cohesive
);
#endif
if
(
this
->
mesh
.
isDistributed
())
{
/// create the distributed synchronizer for cohesive elements
this
->
cohesive_synchronizer
=
std
::
make_unique
<
ElementSynchronizer
>
(
mesh
,
"cohesive_distributed_synchronizer"
);
auto
&
synchronizer
=
mesh
.
getElementSynchronizer
();
this
->
cohesive_synchronizer
->
split
(
synchronizer
,
[](
auto
&&
el
)
{
return
Mesh
::
getKind
(
el
.
type
)
==
_ek_cohesive
;
});
this
->
registerSynchronizer
(
*
cohesive_synchronizer
,
_gst_material_id
);
this
->
registerSynchronizer
(
*
cohesive_synchronizer
,
_gst_smm_stress
);
this
->
registerSynchronizer
(
*
cohesive_synchronizer
,
_gst_smm_boundary
);
}
this
->
inserter
=
std
::
make_unique
<
CohesiveElementInserter
>
(
this
->
mesh
,
id
+
":cohesive_element_inserter"
);
AKANTU_DEBUG_OUT
();
}
/* -------------------------------------------------------------------------- */
SolidMechanicsModelCohesive
::~
SolidMechanicsModelCohesive
()
=
default
;
/* -------------------------------------------------------------------------- */
void
SolidMechanicsModelCohesive
::
setTimeStep
(
Real
time_step
,
const
ID
&
solver_id
)
{
SolidMechanicsModel
::
setTimeStep
(
time_step
,
solver_id
);
#if defined(AKANTU_USE_IOHELPER)
this
->
mesh
.
getDumper
(
"cohesive elements"
).
setTimeStep
(
time_step
);
#endif
}
/* -------------------------------------------------------------------------- */
void
SolidMechanicsModelCohesive
::
initFullImpl
(
const
ModelOptions
&
options
)
{
AKANTU_DEBUG_IN
();
const
auto
&
smmc_options
=
dynamic_cast
<
const
SolidMechanicsModelCohesiveOptions
&>
(
options
);
this
->
is_extrinsic
=
smmc_options
.
extrinsic
;
inserter
->
setIsExtrinsic
(
is_extrinsic
);
if
(
mesh
.
isDistributed
())
{
auto
&
mesh_facets
=
inserter
->
getMeshFacets
();
auto
&
synchronizer
=
dynamic_cast
<
FacetSynchronizer
&>
(
mesh_facets
.
getElementSynchronizer
());
this
->
registerSynchronizer
(
synchronizer
,
_gst_smmc_facets
);
this
->
registerSynchronizer
(
synchronizer
,
_gst_smmc_facets_conn
);
synchronizeGhostFacetsConnectivity
();
/// create the facet synchronizer for extrinsic simulations
if
(
is_extrinsic
)
{
facet_stress_synchronizer
=
std
::
make_unique
<
ElementSynchronizer
>
(
mesh_facets
);
facet_stress_synchronizer
->
updateSchemes
(
[
&
](
auto
&
scheme
,
auto
&
proc
,
auto
&
direction
)
{
scheme
.
copy
(
const_cast
<
const
FacetSynchronizer
&>
(
synchronizer
)
.
getCommunications
()
.
getScheme
(
proc
,
direction
));
});
this
->
registerSynchronizer
(
*
facet_stress_synchronizer
,
_gst_smmc_facets_stress
);
}
inserter
->
initParallel
(
*
cohesive_synchronizer
);
}
ParserSection
section
;
bool
is_empty
;
std
::
tie
(
section
,
is_empty
)
=
this
->
getParserSection
();
if
(
not
is_empty
)
{
auto
inserter_section
=
section
.
getSubSections
(
ParserType
::
_cohesive_inserter
);
if
(
inserter_section
.
begin
()
!=
inserter_section
.
end
())
{
inserter
->
parseSection
(
*
inserter_section
.
begin
());
}
}
SolidMechanicsModel
::
initFullImpl
(
options
);
AKANTU_DEBUG_OUT
();
}
// namespace akantu
/* -------------------------------------------------------------------------- */
void
SolidMechanicsModelCohesive
::
initMaterials
()
{
AKANTU_DEBUG_IN
();
// make sure the material are instantiated
if
(
!
are_materials_instantiated
)
instantiateMaterials
();
/// find the first cohesive material
UInt
cohesive_index
=
UInt
(
-
1
);
for
(
auto
&&
material
:
enumerate
(
materials
))
{
if
(
dynamic_cast
<
MaterialCohesive
*>
(
std
::
get
<
1
>
(
material
).
get
()))
{
cohesive_index
=
std
::
get
<
0
>
(
material
);
break
;
}
}
if
(
cohesive_index
==
UInt
(
-
1
))
AKANTU_EXCEPTION
(
"No cohesive materials in the material input file"
);
material_selector
->
setFallback
(
cohesive_index
);
// set the facet information in the material in case of dynamic insertion
// to know what material to call for stress checks
if
(
is_extrinsic
)
{
this
->
initExtrinsicMaterials
();
}
else
{
this
->
initIntrinsicMaterials
();
}
AKANTU_DEBUG_OUT
();
}
// namespace akantu
/* -------------------------------------------------------------------------- */
void
SolidMechanicsModelCohesive
::
initExtrinsicMaterials
()
{
const
Mesh
&
mesh_facets
=
inserter
->
getMeshFacets
();
facet_material
.
initialize
(
mesh_facets
,
_spatial_dimension
=
spatial_dimension
-
1
,
_with_nb_element
=
true
,
_default_value
=
material_selector
->
getFallbackValue
());
for_each_elements
(
mesh_facets
,
[
&
](
auto
&&
element
)
{
auto
mat_index
=
(
*
material_selector
)(
element
);
auto
&
mat
=
dynamic_cast
<
MaterialCohesive
&>
(
*
materials
[
mat_index
]);
facet_material
(
element
)
=
mat_index
;
mat
.
addFacet
(
element
);
},
_spatial_dimension
=
spatial_dimension
-
1
);
SolidMechanicsModel
::
initMaterials
();
this
->
initAutomaticInsertion
();
}
/* -------------------------------------------------------------------------- */
#if 0
void SolidMechanicsModelCohesive::initIntrinsicCohesiveMaterials(
const std::string & cohesive_surfaces) {
AKANTU_DEBUG_IN();
#if defined(AKANTU_PARALLEL_COHESIVE_ELEMENT)
if (facet_synchronizer != nullptr)
inserter->initParallel(facet_synchronizer, cohesive_element_synchronizer);
// inserter->initParallel(facet_synchronizer, synch_parallel);
#endif
std::istringstream split(cohesive_surfaces);
std::string physname;
while (std::getline(split, physname, ',')) {
AKANTU_DEBUG_INFO(
"Pre-inserting cohesive elements along facets from physical surface: "
<< physname);
insertElementsFromMeshData(physname);
}
synchronizeInsertionData();
SolidMechanicsModel::initMaterials();
auto && tmp = std::make_shared<MeshDataMaterialCohesiveSelector>(*this);
tmp->setFallback(material_selector->getFallbackValue());
tmp->setFallback(material_selector->getFallbackSelector());
material_selector = tmp;
// UInt nb_new_elements =
inserter->insertElements();
// if (nb_new_elements > 0) {
// this->reinitializeSolver();
// }
AKANTU_DEBUG_OUT();
}
#endif
/* -------------------------------------------------------------------------- */
void
SolidMechanicsModelCohesive
::
synchronizeInsertionData
()
{
if
(
inserter
->
getMeshFacets
().
isDistributed
())
{
inserter
->
getMeshFacets
().
getElementSynchronizer
().
synchronizeOnce
(
*
inserter
,
_gst_ce_groups
);
}
}
/* -------------------------------------------------------------------------- */
void
SolidMechanicsModelCohesive
::
initIntrinsicMaterials
()
{
AKANTU_DEBUG_IN
();
SolidMechanicsModel
::
initMaterials
();
this
->
insertIntrinsicElements
();
AKANTU_DEBUG_OUT
();
}
/* -------------------------------------------------------------------------- */
/**
* Initialize the model,basically it pre-compute the shapes, shapes derivatives
* and jacobian
*/
void
SolidMechanicsModelCohesive
::
initModel
()
{
AKANTU_DEBUG_IN
();
SolidMechanicsModel
::
initModel
();
registerFEEngineObject
<
MyFEEngineCohesiveType
>
(
"CohesiveFEEngine"
,
mesh
,
Model
::
spatial_dimension
);
/// add cohesive type connectivity
ElementType
type
=
_not_defined
;
for
(
auto
&&
type_ghost
:
ghost_types
)
{
for
(
const
auto
&
tmp_type
:
mesh
.
elementTypes
(
spatial_dimension
,
type_ghost
))
{
const
auto
&
connectivity
=
mesh
.
getConnectivity
(
tmp_type
,
type_ghost
);
if
(
connectivity
.
size
()
==
0
)
continue
;
type
=
tmp_type
;
auto
type_facet
=
Mesh
::
getFacetType
(
type
);
auto
type_cohesive
=
FEEngine
::
getCohesiveElementType
(
type_facet
);
mesh
.
addConnectivityType
(
type_cohesive
,
type_ghost
);
}
}
AKANTU_DEBUG_ASSERT
(
type
!=
_not_defined
,
"No elements in the mesh"
);
getFEEngine
(
"CohesiveFEEngine"
).
initShapeFunctions
(
_not_ghost
);
getFEEngine
(
"CohesiveFEEngine"
).
initShapeFunctions
(
_ghost
);
registerFEEngineObject
<
MyFEEngineFacetType
>
(
"FacetsFEEngine"
,
mesh
.
getMeshFacets
(),
Model
::
spatial_dimension
-
1
);
if
(
is_extrinsic
)
{
getFEEngine
(
"FacetsFEEngine"
).
initShapeFunctions
(
_not_ghost
);
getFEEngine
(
"FacetsFEEngine"
).
initShapeFunctions
(
_ghost
);
}
AKANTU_DEBUG_OUT
();
}
/* -------------------------------------------------------------------------- */
void
SolidMechanicsModelCohesive
::
insertIntrinsicElements
()
{
AKANTU_DEBUG_IN
();
inserter
->
insertIntrinsicElements
();
AKANTU_DEBUG_OUT
();
}
/* -------------------------------------------------------------------------- */
void
SolidMechanicsModelCohesive
::
updateFacetStressSynchronizer
()
{
if
(
facet_stress_synchronizer
!=
nullptr
)
{
const
auto
&
rank_to_element
=
mesh
.
getElementSynchronizer
().
getElementToRank
();
const
auto
&
facet_checks
=
inserter
->
getCheckFacets
();
const
auto
&
element_to_facet
=
mesh
.
getElementToSubelement
();
UInt
rank
=
mesh
.
getCommunicator
().
whoAmI
();
facet_stress_synchronizer
->
updateSchemes
(
[
&
](
auto
&
scheme
,
auto
&
proc
,
auto
&
/*direction*/
)
{
UInt
el
=
0
;
for
(
auto
&&
element
:
scheme
)
{
if
(
not
facet_checks
(
element
))
continue
;
const
auto
&
next_el
=
element_to_facet
(
element
);
UInt
rank_left
=
rank_to_element
(
next_el
[
0
]);
UInt
rank_right
=
rank_to_element
(
next_el
[
1
]);
if
((
rank_left
==
rank
and
rank_right
==
proc
)
or
(
rank_left
==
proc
and
rank_right
==
rank
))
{
scheme
[
el
]
=
element
;
++
el
;
}
}
scheme
.
resize
(
el
);
});
}
}
/* -------------------------------------------------------------------------- */
void
SolidMechanicsModelCohesive
::
initAutomaticInsertion
()
{
AKANTU_DEBUG_IN
();
this
->
inserter
->
limitCheckFacets
();
this
->
updateFacetStressSynchronizer
();
this
->
resizeFacetStress
();
/// compute normals on facets
this
->
computeNormals
();
this
->
initStressInterpolation
();
}
/* -------------------------------------------------------------------------- */
void
SolidMechanicsModelCohesive
::
updateAutomaticInsertion
()
{
AKANTU_DEBUG_IN
();
this
->
inserter
->
limitCheckFacets
();
this
->
updateFacetStressSynchronizer
();
AKANTU_DEBUG_OUT
();
}
/* -------------------------------------------------------------------------- */
void
SolidMechanicsModelCohesive
::
initStressInterpolation
()
{
Mesh
&
mesh_facets
=
inserter
->
getMeshFacets
();
/// compute quadrature points coordinates on facets
Array
<
Real
>
&
position
=
mesh
.
getNodes
();
ElementTypeMapArray
<
Real
>
quad_facets
(
"quad_facets"
,
id
);
quad_facets
.
initialize
(
mesh_facets
,
_nb_component
=
Model
::
spatial_dimension
,
_spatial_dimension
=
Model
::
spatial_dimension
-
1
);
// mesh_facets.initElementTypeMapArray(quad_facets, Model::spatial_dimension,
// Model::spatial_dimension - 1);
getFEEngine
(
"FacetsFEEngine"
)
.
interpolateOnIntegrationPoints
(
position
,
quad_facets
);
/// compute elements quadrature point positions and build
/// element-facet quadrature points data structure
ElementTypeMapArray
<
Real
>
elements_quad_facets
(
"elements_quad_facets"
,
id
);
elements_quad_facets
.
initialize
(
mesh
,
_nb_component
=
Model
::
spatial_dimension
,
_spatial_dimension
=
Model
::
spatial_dimension
);
// mesh.initElementTypeMapArray(elements_quad_facets,
// Model::spatial_dimension,
// Model::spatial_dimension);
for
(
auto
elem_gt
:
ghost_types
)
{
for
(
auto
&
type
:
mesh
.
elementTypes
(
Model
::
spatial_dimension
,
elem_gt
))
{
UInt
nb_element
=
mesh
.
getNbElement
(
type
,
elem_gt
);
if
(
nb_element
==
0
)
continue
;
/// compute elements' quadrature points and list of facet
/// quadrature points positions by element
Array
<
Element
>
&
facet_to_element
=
mesh_facets
.
getSubelementToElement
(
type
,
elem_gt
);
UInt
nb_facet_per_elem
=
facet_to_element
.
getNbComponent
();
Array
<
Real
>
&
el_q_facet
=
elements_quad_facets
(
type
,
elem_gt
);
ElementType
facet_type
=
Mesh
::
getFacetType
(
type
);
UInt
nb_quad_per_facet
=
getFEEngine
(
"FacetsFEEngine"
).
getNbIntegrationPoints
(
facet_type
);
el_q_facet
.
resize
(
nb_element
*
nb_facet_per_elem
*
nb_quad_per_facet
);
for
(
UInt
el
=
0
;
el
<
nb_element
;
++
el
)
{
for
(
UInt
f
=
0
;
f
<
nb_facet_per_elem
;
++
f
)
{
Element
global_facet_elem
=
facet_to_element
(
el
,
f
);
UInt
global_facet
=
global_facet_elem
.
element
;
GhostType
facet_gt
=
global_facet_elem
.
ghost_type
;
const
Array
<
Real
>
&
quad_f
=
quad_facets
(
facet_type
,
facet_gt
);
for
(
UInt
q
=
0
;
q
<
nb_quad_per_facet
;
++
q
)
{
for
(
UInt
s
=
0
;
s
<
Model
::
spatial_dimension
;
++
s
)
{
el_q_facet
(
el
*
nb_facet_per_elem
*
nb_quad_per_facet
+
f
*
nb_quad_per_facet
+
q
,
s
)
=
quad_f
(
global_facet
*
nb_quad_per_facet
+
q
,
s
);
}
}
}
}
}
}
/// loop over non cohesive materials
for
(
auto
&&
material
:
materials
)
{
if
(
dynamic_cast
<
MaterialCohesive
*>
(
material
.
get
()))
continue
;
/// initialize the interpolation function
material
->
initElementalFieldInterpolation
(
elements_quad_facets
);
}
AKANTU_DEBUG_OUT
();
}
/* -------------------------------------------------------------------------- */
void
SolidMechanicsModelCohesive
::
assembleInternalForces
()
{
AKANTU_DEBUG_IN
();
// f_int += f_int_cohe
for
(
auto
&
material
:
this
->
materials
)
{
try
{
auto
&
mat
=
dynamic_cast
<
MaterialCohesive
&>
(
*
material
);
mat
.
computeTraction
(
_not_ghost
);
}
catch
(
std
::
bad_cast
&
bce
)
{
}
}
SolidMechanicsModel
::
assembleInternalForces
();
if
(
isDefaultSolverExplicit
())
{
for
(
auto
&
material
:
materials
)
{
try
{
auto
&
mat
=
dynamic_cast
<
MaterialCohesive
&>
(
*
material
);
mat
.
computeEnergies
();
}
catch
(
std
::
bad_cast
&
bce
)
{
}
}
}
AKANTU_DEBUG_OUT
();
}
/* -------------------------------------------------------------------------- */
void
SolidMechanicsModelCohesive
::
computeNormals
()
{
AKANTU_DEBUG_IN
();
Mesh
&
mesh_facets
=
this
->
inserter
->
getMeshFacets
();
this
->
getFEEngine
(
"FacetsFEEngine"
)
.
computeNormalsOnIntegrationPoints
(
_not_ghost
);
/**
* @todo store tangents while computing normals instead of
* recomputing them as follows:
*/
/* ------------------------------------------------------------------------ */
UInt
tangent_components
=
Model
::
spatial_dimension
*
(
Model
::
spatial_dimension
-
1
);
tangents
.
initialize
(
mesh_facets
,
_nb_component
=
tangent_components
,
_spatial_dimension
=
Model
::
spatial_dimension
-
1
);
// mesh_facets.initElementTypeMapArray(tangents, tangent_components,
// Model::spatial_dimension - 1);
for
(
auto
facet_type
:
mesh_facets
.
elementTypes
(
Model
::
spatial_dimension
-
1
))
{
const
Array
<
Real
>
&
normals
=
this
->
getFEEngine
(
"FacetsFEEngine"
)
.
getNormalsOnIntegrationPoints
(
facet_type
);
Array
<
Real
>
&
tangents
=
this
->
tangents
(
facet_type
);
Math
::
compute_tangents
(
normals
,
tangents
);
}
AKANTU_DEBUG_OUT
();
}
/* -------------------------------------------------------------------------- */
void
SolidMechanicsModelCohesive
::
interpolateStress
()
{
ElementTypeMapArray
<
Real
>
by_elem_result
(
"temporary_stress_by_facets"
,
id
);
for
(
auto
&
material
:
materials
)
{
try
{
auto
&
mat
__attribute__
((
unused
))
=
dynamic_cast
<
MaterialCohesive
&>
(
*
material
);
}
catch
(
std
::
bad_cast
&
)
{
/// interpolate stress on facet quadrature points positions
material
->
interpolateStressOnFacets
(
facet_stress
,
by_elem_result
);
}
}
#if defined(AKANTU_DEBUG_TOOLS)
debug
::
element_manager
.
printData
(
debug
::
_dm_model_cohesive
,
"Interpolated stresses before"
,
facet_stress
);
#endif
this
->
synchronize
(
_gst_smmc_facets_stress
);
#if defined(AKANTU_DEBUG_TOOLS)
debug
::
element_manager
.
printData
(
debug
::
_dm_model_cohesive
,
"Interpolated stresses"
,
facet_stress
);
#endif
}
/* -------------------------------------------------------------------------- */
UInt
SolidMechanicsModelCohesive
::
checkCohesiveStress
()
{
interpolateStress
();
for
(
auto
&
mat
:
materials
)
{
try
{
auto
&
mat_cohesive
=
dynamic_cast
<
MaterialCohesive
&>
(
*
mat
);
/// check which not ghost cohesive elements are to be created
mat_cohesive
.
checkInsertion
();
}
catch
(
std
::
bad_cast
&
)
{
}
}
/// communicate data among processors
this
->
synchronize
(
_gst_smmc_facets
);
/// insert cohesive elements
UInt
nb_new_elements
=
inserter
->
insertElements
();
// if (nb_new_elements > 0) {
// this->reinitializeSolver();
// }
return
nb_new_elements
;
}
/* -------------------------------------------------------------------------- */
void
SolidMechanicsModelCohesive
::
onElementsAdded
(
const
Array
<
Element
>
&
element_list
,
const
NewElementsEvent
&
event
)
{
AKANTU_DEBUG_IN
();
updateCohesiveSynchronizers
();
SolidMechanicsModel
::
onElementsAdded
(
element_list
,
event
);
if
(
cohesive_synchronizer
!=
nullptr
)
cohesive_synchronizer
->
computeAllBufferSizes
(
*
this
);
if
(
is_extrinsic
)
resizeFacetStress
();
/// if (method != _explicit_lumped_mass) {
/// this->initSolver();
/// }
AKANTU_DEBUG_OUT
();
}
/* -------------------------------------------------------------------------- */
void
SolidMechanicsModelCohesive
::
onNodesAdded
(
const
Array
<
UInt
>
&
new_nodes
,
const
NewNodesEvent
&
event
)
{
AKANTU_DEBUG_IN
();
// Array<UInt> nodes_list(nb_new_nodes);
// for (UInt n = 0; n < nb_new_nodes; ++n)
// nodes_list(n) = doubled_nodes(n, 1);
SolidMechanicsModel
::
onNodesAdded
(
new_nodes
,
event
);
UInt
new_node
,
old_node
;
try
{
const
auto
&
cohesive_event
=
dynamic_cast
<
const
CohesiveNewNodesEvent
&>
(
event
);
const
auto
&
old_nodes
=
cohesive_event
.
getOldNodesList
();
auto
copy
=
[
this
,
&
new_node
,
&
old_node
](
auto
&
arr
)
{
for
(
UInt
s
=
0
;
s
<
spatial_dimension
;
++
s
)
{
arr
(
new_node
,
s
)
=
arr
(
old_node
,
s
);
}
};
for
(
auto
&&
pair
:
zip
(
new_nodes
,
old_nodes
))
{
std
::
tie
(
new_node
,
old_node
)
=
pair
;
copy
(
*
displacement
);
copy
(
*
blocked_dofs
);
if
(
velocity
)
copy
(
*
velocity
);
if
(
acceleration
)
copy
(
*
acceleration
);
if
(
current_position
)
copy
(
*
current_position
);
if
(
previous_displacement
)
copy
(
*
previous_displacement
);
}
// if (this->getDOFManager().hasMatrix("M")) {
// this->assembleMass(old_nodes);
// }
// if (this->getDOFManager().hasLumpedMatrix("M")) {
// this->assembleMassLumped(old_nodes);
// }
}
catch
(
std
::
bad_cast
&
)
{
}
AKANTU_DEBUG_OUT
();
}
/* -------------------------------------------------------------------------- */
void
SolidMechanicsModelCohesive
::
onEndSolveStep
(
const
AnalysisMethod
&
)
{
AKANTU_DEBUG_IN
();
/*
* This is required because the Cauchy stress is the stress measure that
* is used to check the insertion of cohesive elements
*/
for
(
auto
&
mat
:
materials
)
{
if
(
mat
->
isFiniteDeformation
())
mat
->
computeAllCauchyStresses
(
_not_ghost
);
}
AKANTU_DEBUG_OUT
();
}
/* -------------------------------------------------------------------------- */
void
SolidMechanicsModelCohesive
::
printself
(
std
::
ostream
&
stream
,
int
indent
)
const
{
std
::
string
space
;
for
(
Int
i
=
0
;
i
<
indent
;
i
++
,
space
+=
AKANTU_INDENT
)
;
stream
<<
space
<<
"SolidMechanicsModelCohesive ["
<<
std
::
endl
;
SolidMechanicsModel
::
printself
(
stream
,
indent
+
1
);
stream
<<
space
<<
"]"
<<
std
::
endl
;
}
/* -------------------------------------------------------------------------- */
void
SolidMechanicsModelCohesive
::
resizeFacetStress
()
{
AKANTU_DEBUG_IN
();
this
->
facet_stress
.
initialize
(
getFEEngine
(
"FacetsFEEngine"
),
_nb_component
=
2
*
spatial_dimension
*
spatial_dimension
,
_spatial_dimension
=
spatial_dimension
-
1
);
// for (auto && ghost_type : ghost_types) {
// for (const auto & type :
// mesh_facets.elementTypes(spatial_dimension - 1, ghost_type)) {
// UInt nb_facet = mesh_facets.getNbElement(type, ghost_type);
// UInt nb_quadrature_points = getFEEngine("FacetsFEEngine")
// .getNbIntegrationPoints(type,
// ghost_type);
// UInt new_size = nb_facet * nb_quadrature_points;
// facet_stress(type, ghost_type).resize(new_size);
// }
// }
AKANTU_DEBUG_OUT
();
}
/* -------------------------------------------------------------------------- */
void
SolidMechanicsModelCohesive
::
addDumpGroupFieldToDumper
(
const
std
::
string
&
dumper_name
,
const
std
::
string
&
field_id
,
const
std
::
string
&
group_name
,
const
ElementKind
&
element_kind
,
bool
padding_flag
)
{
AKANTU_DEBUG_IN
();
UInt
spatial_dimension
=
Model
::
spatial_dimension
;
ElementKind
_element_kind
=
element_kind
;
if
(
dumper_name
==
"cohesive elements"
)
{
_element_kind
=
_ek_cohesive
;
}
else
if
(
dumper_name
==
"facets"
)
{
spatial_dimension
=
Model
::
spatial_dimension
-
1
;
}
SolidMechanicsModel
::
addDumpGroupFieldToDumper
(
dumper_name
,
field_id
,
group_name
,
spatial_dimension
,
_element_kind
,
padding_flag
);
AKANTU_DEBUG_OUT
();
}
/* -------------------------------------------------------------------------- */
void
SolidMechanicsModelCohesive
::
onDump
()
{
this
->
flattenAllRegisteredInternals
(
_ek_cohesive
);
SolidMechanicsModel
::
onDump
();
}
/* -------------------------------------------------------------------------- */
}
// namespace akantu
Event Timeline
Log In to Comment