Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F84788049
material_brittle_inline_impl.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
Tue, Sep 24, 21:36
Size
3 KB
Mime Type
text/x-c++
Expires
Thu, Sep 26, 21:36 (2 d)
Engine
blob
Format
Raw Data
Handle
21077782
Attached To
rAKA akantu
material_brittle_inline_impl.cc
View Options
/**
* @file material_brittle_inline_impl.cc
*
* @author Aranda Ruiz Josue <josue.arandaruiz@epfl.ch>
* @author Daniel Pino Muñoz <daniel.pinomunoz@epfl.ch>
*
*
* @brief Implementation of the inline functions of the material brittle
*
* @section LICENSE
*
* Copyright (©) 2010-2012, 2014 EPFL (Ecole Polytechnique Fédérale de Lausanne)
* Laboratory (LSMS - Laboratoire de Simulation en Mécanique des Solides)
*
*/
/* -------------------------------------------------------------------------- */
template
<
UInt
spatial_dimension
>
inline
void
MaterialBrittle
<
spatial_dimension
>::
computeStressOnQuad
(
Matrix
<
Real
>
&
grad_u
,
Matrix
<
Real
>
&
grad_v
,
Matrix
<
Real
>
&
sigma
,
Real
&
dam
,
Real
&
sigma_equivalent
,
Real
&
fracture_stress
)
{
MaterialElastic
<
spatial_dimension
>::
computeStressOnQuad
(
grad_u
,
sigma
);
Real
equiv_strain_rate
=
0.
;
Real
volume_change_rate
=
grad_v
.
trace
();
if
(
spatial_dimension
==
2
)
{
equiv_strain_rate
+=
2.
/
3.
*
pow
(
volume_change_rate
/
3.
,
2.
);
}
for
(
UInt
i
=
0
;
i
<
spatial_dimension
;
++
i
)
for
(
UInt
j
=
0
;
j
<
spatial_dimension
;
++
j
)
equiv_strain_rate
+=
2.
/
3.
*
pow
(
0.5
*
(
grad_v
(
i
,
j
)
+
grad_v
(
j
,
i
))
-
(
i
==
j
)
*
volume_change_rate
/
3.
,
2.
);
equiv_strain_rate
=
sqrt
(
equiv_strain_rate
);
fracture_stress
=
S_0
;
if
(
equiv_strain_rate
>
E_0
)
fracture_stress
=
A
;
Vector
<
Real
>
principal_stress
(
spatial_dimension
);
sigma
.
eig
(
principal_stress
);
sigma_equivalent
=
principal_stress
(
0
);
for
(
UInt
i
=
1
;
i
<
spatial_dimension
;
++
i
)
sigma_equivalent
=
std
::
max
(
sigma_equivalent
,
principal_stress
(
i
));
if
(
!
this
->
is_non_local
)
{
computeDamageAndStressOnQuad
(
sigma
,
dam
,
sigma_equivalent
,
fracture_stress
);
}
}
/* -------------------------------------------------------------------------- */
template
<
UInt
spatial_dimension
>
inline
void
MaterialBrittle
<
spatial_dimension
>::
computeDamageAndStressOnQuad
(
Matrix
<
Real
>
&
sigma
,
Real
&
dam
,
Real
&
sigma_c
,
Real
&
fracture_stress
)
{
if
(
sigma_c
>
fracture_stress
)
dam
=
1.
;
dam
=
std
::
min
(
dam
,
1.
);
sigma
*=
1
-
dam
;
}
/* -------------------------------------------------------------------------- */
template
<
UInt
spatial_dimension
>
inline
UInt
MaterialBrittle
<
spatial_dimension
>::
getNbData
(
const
Array
<
Element
>
&
elements
,
const
SynchronizationTag
&
tag
)
const
{
AKANTU_DEBUG_IN
();
UInt
size
=
MaterialDamage
<
spatial_dimension
>::
getNbData
(
elements
,
tag
);
AKANTU_DEBUG_OUT
();
return
size
;
}
/* -------------------------------------------------------------------------- */
template
<
UInt
spatial_dimension
>
inline
void
MaterialBrittle
<
spatial_dimension
>::
packData
(
CommunicationBuffer
&
buffer
,
const
Array
<
Element
>
&
elements
,
const
SynchronizationTag
&
tag
)
const
{
AKANTU_DEBUG_IN
();
MaterialDamage
<
spatial_dimension
>::
packData
(
buffer
,
elements
,
tag
);
AKANTU_DEBUG_OUT
();
}
/* -------------------------------------------------------------------------- */
template
<
UInt
spatial_dimension
>
inline
void
MaterialBrittle
<
spatial_dimension
>::
unpackData
(
CommunicationBuffer
&
buffer
,
const
Array
<
Element
>
&
elements
,
const
SynchronizationTag
&
tag
)
{
AKANTU_DEBUG_IN
();
MaterialDamage
<
spatial_dimension
>::
unpackData
(
buffer
,
elements
,
tag
);
AKANTU_DEBUG_OUT
();
}
Event Timeline
Log In to Comment