Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F79721282
material_non_local.hh
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
Tue, Aug 27, 22:15
Size
6 KB
Mime Type
text/x-c++
Expires
Thu, Aug 29, 22:15 (1 d, 22 h)
Engine
blob
Format
Raw Data
Handle
20134019
Attached To
rAKA akantu
material_non_local.hh
View Options
/**
* @file material_non_local.hh
*
* @author Nicolas Richart <nicolas.richart@epfl.ch>
*
* @date creation: Wed Aug 31 2011
* @date last modification: Thu Jun 05 2014
*
* @brief Material class that handle the non locality of a law for example damage.
*
* @section LICENSE
*
* Copyright (©) 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 "material.hh"
#include "aka_grid_dynamic.hh"
#include "fe_engine.hh"
#include "weight_function.hh"
namespace
akantu
{
class
GridSynchronizer
;
}
/* -------------------------------------------------------------------------- */
#ifndef __AKANTU_MATERIAL_NON_LOCAL_HH__
#define __AKANTU_MATERIAL_NON_LOCAL_HH__
__BEGIN_AKANTU__
/* -------------------------------------------------------------------------- */
template
<
UInt
dim
,
template
<
UInt
>
class
WeightFunction
=
BaseWeightFunction
>
class
MaterialNonLocal
:
public
virtual
Material
{
/* ------------------------------------------------------------------------ */
/* Constructors/Destructors */
/* ------------------------------------------------------------------------ */
public
:
MaterialNonLocal
(
SolidMechanicsModel
&
model
,
const
ID
&
id
=
""
);
virtual
~
MaterialNonLocal
();
typedef
std
::
vector
<
std
::
pair
<
QuadraturePoint
,
QuadraturePoint
>
>
PairList
;
/* ------------------------------------------------------------------------ */
/* Methods */
/* ------------------------------------------------------------------------ */
public
:
/// initialize the material computed parameter
virtual
void
initMaterial
();
virtual
void
updateResidual
(
GhostType
ghost_type
);
virtual
void
computeAllNonLocalStresses
(
GhostType
ghost_type
=
_not_ghost
);
void
savePairs
(
const
std
::
string
&
filename
)
const
;
void
neighbourhoodStatistics
(
const
std
::
string
&
filename
)
const
;
protected
:
void
updatePairList
(
const
ElementTypeMapArray
<
Real
>
&
quadrature_points_coordinates
);
void
computeWeights
(
const
ElementTypeMapArray
<
Real
>
&
quadrature_points_coordinates
);
void
createCellList
(
ElementTypeMapArray
<
Real
>
&
quadrature_points_coordinates
);
void
cleanupExtraGhostElement
(
const
ElementTypeMap
<
UInt
>
&
nb_ghost_protected
);
void
fillCellList
(
const
ElementTypeMapArray
<
Real
>
&
quadrature_points_coordinates
,
const
GhostType
&
ghost_type
);
/// constitutive law
virtual
void
computeNonLocalStresses
(
GhostType
ghost_type
=
_not_ghost
)
=
0
;
template
<
typename
T
>
void
weightedAvergageOnNeighbours
(
const
InternalField
<
T
>
&
to_accumulate
,
InternalField
<
T
>
&
accumulated
,
UInt
nb_degree_of_freedom
,
GhostType
ghost_type2
=
_not_ghost
)
const
;
virtual
inline
UInt
getNbDataForElements
(
const
Array
<
Element
>
&
elements
,
SynchronizationTag
tag
)
const
;
virtual
inline
void
packElementData
(
CommunicationBuffer
&
buffer
,
const
Array
<
Element
>
&
elements
,
SynchronizationTag
tag
)
const
;
virtual
inline
void
unpackElementData
(
CommunicationBuffer
&
buffer
,
const
Array
<
Element
>
&
elements
,
SynchronizationTag
tag
);
virtual
inline
void
onElementsAdded
(
const
Array
<
Element
>
&
element_list
);
virtual
inline
void
onElementsRemoved
(
const
Array
<
Element
>
&
element_list
,
const
ElementTypeMapArray
<
UInt
>
&
new_numbering
,
const
RemovedElementsEvent
&
event
);
/* ------------------------------------------------------------------------ */
/* Accessors */
/* ------------------------------------------------------------------------ */
public
:
void
registerNonLocalVariable
(
InternalField
<
Real
>
&
local
,
InternalField
<
Real
>
&
non_local
,
UInt
nb_degree_of_freedom
)
{
ID
id
=
local
.
getID
();
NonLocalVariable
&
non_local_variable
=
non_local_variables
[
id
];
non_local_variable
.
local
=
&
local
;
non_local_variable
.
non_local
=
&
non_local
;
non_local_variable
.
nb_component
=
nb_degree_of_freedom
;
}
AKANTU_GET_MACRO
(
PairList
,
pair_list
,
const
PairList
&
)
Real
getRadius
()
const
{
return
weight_func
->
getRadius
();
}
AKANTU_GET_MACRO
(
CellList
,
*
spatial_grid
,
const
SpatialGrid
<
QuadraturePoint
>
&
)
/* ------------------------------------------------------------------------ */
/* Class Members */
/* ------------------------------------------------------------------------ */
protected
:
/// the weight function used
WeightFunction
<
dim
>
*
weight_func
;
private
:
/**
* the pairs of quadrature points
* 0: not ghost to not ghost
* 1: not ghost to ghost
*/
PairList
pair_list
[
2
];
/// the weights associated to the pairs
Array
<
Real
>
*
pair_weight
[
2
];
/// the regular grid to construct/update the pair lists
SpatialGrid
<
QuadraturePoint
>
*
spatial_grid
;
/// the types of the existing pairs
typedef
std
::
set
<
std
::
pair
<
ElementType
,
ElementType
>
>
pair_type
;
pair_type
existing_pairs
[
2
];
/// count the number of calls of computeStress
UInt
compute_stress_calls
;
struct
NonLocalVariable
{
InternalField
<
Real
>
*
local
;
InternalField
<
Real
>
*
non_local
;
UInt
nb_component
;
};
std
::
map
<
ID
,
NonLocalVariable
>
non_local_variables
;
bool
is_creating_grid
;
GridSynchronizer
*
grid_synchronizer
;
};
/* -------------------------------------------------------------------------- */
/* inline functions */
/* -------------------------------------------------------------------------- */
#include "material_non_local_inline_impl.cc"
__END_AKANTU__
#endif
/* __AKANTU_MATERIAL_NON_LOCAL_HH__ */
Event Timeline
Log In to Comment