GRINS-0.8.0
stress_strain_law.h
Go to the documentation of this file.
1 //-----------------------------------------------------------------------bl-
2 //--------------------------------------------------------------------------
3 //
4 // GRINS - General Reacting Incompressible Navier-Stokes
5 //
6 // Copyright (C) 2014-2017 Paul T. Bauman, Roy H. Stogner
7 // Copyright (C) 2010-2013 The PECOS Development Team
8 //
9 // This library is free software; you can redistribute it and/or
10 // modify it under the terms of the Version 2.1 GNU Lesser General
11 // Public License as published by the Free Software Foundation.
12 //
13 // This library is distributed in the hope that it will be useful,
14 // but WITHOUT ANY WARRANTY; without even the implied warranty of
15 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
16 // Lesser General Public License for more details.
17 //
18 // You should have received a copy of the GNU Lesser General Public
19 // License along with this library; if not, write to the Free Software
20 // Foundation, Inc. 51 Franklin Street, Fifth Floor,
21 // Boston, MA 02110-1301 USA
22 //
23 //-----------------------------------------------------------------------el-
24 
25 #ifndef GRINS_STRESS_STRAIN_LAW_H
26 #define GRINS_STRESS_STRAIN_LAW_H
27 
28 // libMesh
29 #include "libmesh/libmesh_common.h"
30 
31 // Forward declarations
32 namespace libMesh
33 {
34  template <typename T>
35  class TensorValue;
36 }
37 
38 namespace GRINS
39 {
40  class ElasticityTensor;
41 
42  template <typename Law>
44  {
45  public:
46 
48  virtual ~StressStrainLaw(){};
49 
50  void compute_stress( unsigned int dim,
56 
57  void compute_stress_and_elasticity( unsigned int dim,
63  ElasticityTensor& C );
64 
66  libMesh::Real compute_33_stress( const libMesh::TensorValue<libMesh::Real>& g_contra,
70  };
71 
72  template <typename Law>
73  inline
74  void StressStrainLaw<Law>::compute_stress( unsigned int dim,
80  {
81  static_cast<Law*>(this)->compute_stress_imp(dim,g_contra,g_cov,G_contra,G_cov,stress);
82  return;
83  }
84 
85  template <typename Law>
86  inline
93  ElasticityTensor& C )
94  {
95  static_cast<Law*>(this)->compute_stress_and_elasticity_imp(dim,g_contra,g_cov,G_contra,G_cov,stress,C);
96  return;
97  }
98 
99  template <typename Law>
100  inline
103  const libMesh::TensorValue<libMesh::Real>& G_contra,
105  {
106  return static_cast<Law*>(this)->compute_33_stress_imp( g_contra, g_cov, G_contra, G_cov );
107  }
108 
109 } // end namespace GRINS
110 
111 #endif //GRINS_STRESS_STRAIN_LAW_H
void compute_stress(unsigned int dim, const libMesh::TensorValue< libMesh::Real > &g_contra, const libMesh::TensorValue< libMesh::Real > &g_cov, const libMesh::TensorValue< libMesh::Real > &G_contra, const libMesh::TensorValue< libMesh::Real > &G_cov, libMesh::TensorValue< libMesh::Real > &stress)
void compute_stress_and_elasticity(unsigned int dim, const libMesh::TensorValue< libMesh::Real > &g_contra, const libMesh::TensorValue< libMesh::Real > &g_cov, const libMesh::TensorValue< libMesh::Real > &G_contra, const libMesh::TensorValue< libMesh::Real > &G_cov, libMesh::TensorValue< libMesh::Real > &stress, ElasticityTensor &C)
GRINS namespace.
libMesh::Real compute_33_stress(const libMesh::TensorValue< libMesh::Real > &g_contra, const libMesh::TensorValue< libMesh::Real > &g_cov, const libMesh::TensorValue< libMesh::Real > &G_contra, const libMesh::TensorValue< libMesh::Real > &G_cov)
This is primarily a helper function for the plane stress cases.

Generated on Tue Dec 19 2017 12:47:28 for GRINS-0.8.0 by  doxygen 1.8.9.1