Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F64927792
ntn_friclaw_linear_slip_weakening_no_healing_tmpl.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
Thu, May 30, 12:19
Size
2 KB
Mime Type
text/x-c++
Expires
Sat, Jun 1, 12:19 (2 d)
Engine
blob
Format
Raw Data
Handle
17979313
Attached To
rAKA akantu
ntn_friclaw_linear_slip_weakening_no_healing_tmpl.hh
View Options
/**
* @file ntn_friclaw_linear_slip_weakening_no_healing_tmpl.hh
*
* @author David Simon Kammer <david.kammer@epfl.ch>
*
*
* @brief implementation of linear slip weakening
*
* @section LICENSE
*
* Copyright (©) 2010-2012, 2014 EPFL (Ecole Polytechnique Fédérale de Lausanne)
* Laboratory (LSMS - Laboratoire de Simulation en Mécanique des Solides)
*
*/
/* -------------------------------------------------------------------------- */
__BEGIN_AKANTU__
/* -------------------------------------------------------------------------- */
template
<
class
Regularisation
>
NTNFricLawLinearSlipWeakeningNoHealing
<
Regularisation
>::
NTNFricLawLinearSlipWeakeningNoHealing
(
NTNBaseContact
*
contact
,
const
FrictionID
&
id
,
const
MemoryID
&
memory_id
)
:
NTNFricLawLinearSlipWeakening
<
Regularisation
>
(
contact
,
id
,
memory_id
)
{
AKANTU_DEBUG_IN
();
AKANTU_DEBUG_OUT
();
}
/* -------------------------------------------------------------------------- */
template
<
class
Regularisation
>
void
NTNFricLawLinearSlipWeakeningNoHealing
<
Regularisation
>::
computeFrictionCoefficient
()
{
AKANTU_DEBUG_IN
();
// get arrays
const
SynchronizedArray
<
Real
>
&
slip
=
this
->
internalGetCumulativeSlip
();
UInt
nb_contact_nodes
=
this
->
contact
->
getNbContactNodes
();
for
(
UInt
n
=
0
;
n
<
nb_contact_nodes
;
++
n
)
{
if
(
slip
(
n
)
>=
this
->
d_c
(
n
))
{
this
->
mu
(
n
)
=
this
->
mu_k
(
n
);
}
else
{
// mu = mu_k + (1 - slip / Dc) * (mu_s - mu_k)
this
->
mu
(
n
)
=
this
->
mu_k
(
n
)
+
(
1
-
(
slip
(
n
)
/
this
->
d_c
(
n
)))
*
(
this
->
mu_s
(
n
)
-
this
->
mu_k
(
n
));
}
}
AKANTU_DEBUG_OUT
();
}
/* -------------------------------------------------------------------------- */
template
<
class
Regularisation
>
void
NTNFricLawLinearSlipWeakeningNoHealing
<
Regularisation
>::
printself
(
std
::
ostream
&
stream
,
int
indent
)
const
{
AKANTU_DEBUG_IN
();
std
::
string
space
;
for
(
Int
i
=
0
;
i
<
indent
;
i
++
,
space
+=
AKANTU_INDENT
);
stream
<<
space
<<
"NTNFricLawLinearSlipWeakeningNoHealing ["
<<
std
::
endl
;
NTNFricLawLinearSlipWeakening
<
Regularisation
>::
printself
(
stream
,
++
indent
);
stream
<<
space
<<
"]"
<<
std
::
endl
;
AKANTU_DEBUG_OUT
();
}
__END_AKANTU__
Event Timeline
Log In to Comment