Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F80770265
material_hyper_elastic1.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
Mon, Sep 2, 12:34
Size
2 KB
Mime Type
text/x-c
Expires
Wed, Sep 4, 12:34 (2 d)
Engine
blob
Format
Raw Data
Handle
20437621
Attached To
rMUSPECTRE µSpectre
material_hyper_elastic1.cc
View Options
/**
* file material_hyper_elastic1.cc
*
* @author Till Junge <till.junge@epfl.ch>
*
* @date 14 Nov 2017
*
* @brief Implementation for materialhyperelastic1
*
* @section LICENCE
*
* Copyright (C) 2017 Till Junge
*
* µSpectre is free software; you can redistribute it and/or
* modify it under the terms of the GNU General Public License as
* published by the Free Software Foundation, either version 3, or (at
* your option) any later version.
*
* µSpectre 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
* General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with GNU Emacs; see the file COPYING. If not, write to the
* Free Software Foundation, Inc., 59 Temple Place - Suite 330,
* Boston, MA 02111-1307, USA.
*/
#include "materials/material_hyper_elastic1.hh"
#include "common/tensor_algebra.hh"
#include <tuple>
namespace
muSpectre
{
/* ---------------------------------------------------------------------- */
template
<
Dim_t
DimS
,
Dim_t
DimM
>
MaterialHyperElastic1
<
DimS
,
DimM
>::
MaterialHyperElastic1
(
std
::
string
name
,
Real
young
,
Real
poisson
)
:
Parent
(
name
),
young
{
young
},
poisson
{
poisson
},
lambda
{
young
*
poisson
/
((
1
+
poisson
)
*
(
1
-
2
*
poisson
))},
mu
{
young
/
(
2
*
(
1
+
poisson
))},
C
{
lambda
*
Tensors
::
outer
<
DimM
>
(
Tensors
::
I2
<
DimM
>
(),
Tensors
::
I2
<
DimM
>
())
+
2
*
mu
*
Tensors
::
I4S
<
DimM
>
()}
{}
/* ---------------------------------------------------------------------- */
template
<
Dim_t
DimS
,
Dim_t
DimM
>
template
<
class
s_t
>
decltype
(
auto
)
MaterialHyperElastic1
<
DimS
,
DimM
>::
evaluate_stress
(
s_t
&&
E
)
{
return
E
.
trace
()
*
lambda
*
Strain_t
::
Identity
()
+
2
*
mu
*
E
;
}
/* ---------------------------------------------------------------------- */
template
<
Dim_t
DimS
,
Dim_t
DimM
>
template
<
class
s_t
>
decltype
(
auto
)
MaterialHyperElastic1
<
DimS
,
DimM
>::
evaluate_stress_tangent
(
s_t
&&
E
)
{
return
std
::
forward_as_tuple
(
this
->
evaluate_stress
(
std
::
move
(
E
)),
this
->
C
);
}
template
class
MaterialHyperElastic1
<
2
,
2
>
;
template
class
MaterialHyperElastic1
<
2
,
3
>
;
template
class
MaterialHyperElastic1
<
3
,
3
>
;
}
// muSpectre
Event Timeline
Log In to Comment