Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F85093155
isotropic_hardening.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, Sep 26, 17:51
Size
3 KB
Mime Type
text/x-c++
Expires
Sat, Sep 28, 17:51 (2 d)
Engine
blob
Format
Raw Data
Handle
21127930
Attached To
rTAMAAS tamaas
isotropic_hardening.hh
View Options
/*
* SPDX-License-Indentifier: AGPL-3.0-or-later
*
* Copyright (©) 2016-2022 EPFL (École Polytechnique Fédérale de Lausanne),
* Laboratory (LSMS - Laboratoire de Simulation en Mécanique des Solides)
*
* This program is free software: you can redistribute it and/or modify
* it under the terms of the GNU Affero General Public License as published
* by the Free Software Foundation, either version 3 of the License, or
* (at your option) any later version.
*
* This program 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 Affero General Public License for more details.
*
* You should have received a copy of the GNU Affero General Public License
* along with this program. If not, see <https://www.gnu.org/licenses/>.
*
*/
/* -------------------------------------------------------------------------- */
#ifndef ISOTROPIC_HARDENING_HH
#define ISOTROPIC_HARDENING_HH
/* -------------------------------------------------------------------------- */
#include "grid.hh"
#include "influence.hh"
#include "internal.hh"
#include "material.hh"
#include "model.hh"
#include "model_type.hh"
#include "static_types.hh"
/* -------------------------------------------------------------------------- */
namespace
tamaas
{
/* -------------------------------------------------------------------------- */
class
IsotropicHardening
:
Material
{
using
parent
=
Material
;
static
constexpr
auto
type
=
model_type
::
volume_2d
;
using
trait
=
model_type_traits
<
type
>
;
static
constexpr
UInt
dim
=
trait
::
dimension
;
using
Field
=
typename
parent
::
Field
;
using
Mat
=
SymMatrixProxy
<
Real
,
dim
>
;
using
CMat
=
SymMatrixProxy
<
const
Real
,
dim
>
;
public
:
/// Constructor
IsotropicHardening
(
Model
*
model
,
Real
sigma_0
,
Real
h
);
/// Compute plastic strain increment with radial return algorithm
void
computeInelasticDeformationIncrement
(
Field
&
increment
,
const
Field
&
strain
,
const
Field
&
strain_increment
)
override
;
/// Compute stress
void
computeStress
(
Field
&
stress
,
const
Field
&
strain
,
const
Field
&
strain_increment
)
override
;
/// Compute stress due to plastic strain increment
void
computeEigenStress
(
Field
&
stress
,
const
Field
&
strain
,
const
Field
&
strain_increment
)
override
;
/// Update internal variables
void
update
()
override
;
void
applyTangentIncrement
(
Field
&
output
,
const
Field
&
input
,
const
Field
&
strain
,
const
Field
&
strain_increment
)
const
;
/// Linear hardening function
static
__device__
__host__
Real
hardening
(
Real
p
,
Real
h
,
Real
sigma_0
)
{
return
sigma_0
+
h
*
p
;
}
Real
getHardeningModulus
()
const
{
return
h
;
}
Real
getYieldStress
()
const
{
return
sigma_0
;
}
const
GridBase
<
Real
>&
getPlasticStrain
()
const
{
return
*
plastic_strain
;
}
GridBase
<
Real
>&
getPlasticStrain
()
{
return
*
plastic_strain
;
}
void
setHardeningModulus
(
Real
h_
)
{
if
(
h_
<
0
)
TAMAAS_EXCEPTION
(
"Hardening modulus should be positive"
);
h
=
h_
;
}
void
setYieldStress
(
Real
sigma_0_
)
{
if
(
sigma_0_
<
0
)
TAMAAS_EXCEPTION
(
"Yield stress should be positive"
);
sigma_0
=
sigma_0_
;
}
protected
:
Real
sigma_0
;
/// < initial yield stress
Real
h
;
/// < hardening modulus
Internal
<
Real
,
dim
>
plastic_strain
,
cumulated_plastic_strain
;
};
}
// namespace tamaas
/* -------------------------------------------------------------------------- */
#endif
// ISOTROPIC_HARDENING_HH
Event Timeline
Log In to Comment