#include <ADComputeIsotropicElasticityTensorShell.h>

Inheritance diagram for ADComputeIsotropicElasticityTensorShell< compute_stage >:

Public Member Functions
	ADComputeIsotropicElasticityTensorShell (const InputParameters &parameters)

Protected Member Functions
virtual void	computeQpProperties () override

Protected Attributes
Real	_poissons_ratio
	Lame' paramters. More...

Real	_shear_modulus

Real	_youngs_modulus

RankFourTensor	_Cijkl
	Individual elasticity tensor. More...

std::vector< Point >	_t_points
	Quadrature points along thickness. More...

std::vector< ADMaterialProperty(RankFourTensor) * >	_elasticity_tensor
	Material property elasticity tensor. More...

std::vector< const ADMaterialProperty(RankTwoTensor) * >	_ge
	Material property for ge matrix. More...

	usingMaterialMembers

Detailed Description

template<ComputeStage compute_stage>
class ADComputeIsotropicElasticityTensorShell< compute_stage >

Definition at line 18 of file ADComputeIsotropicElasticityTensorShell.h.

Constructor & Destructor Documentation

◆ ADComputeIsotropicElasticityTensorShell()

template<ComputeStage compute_stage>

ADComputeIsotropicElasticityTensorShell< compute_stage >::ADComputeIsotropicElasticityTensorShell ( const InputParameters & parameters )

Definition at line 36 of file ADComputeIsotropicElasticityTensorShell.C.

   : ADMaterial<compute_stage>(parameters),
     _poissons_ratio(getParam<Real>("poissons_ratio")),
     _youngs_modulus(getParam<Real>("youngs_modulus"))
 {
   _Cijkl.fillSymmetricIsotropicEandNu(_youngs_modulus, _poissons_ratio);
  
   // correction for plane stress
   _Cijkl(0, 0, 0, 0) = _youngs_modulus / (1.0 - _poissons_ratio * _poissons_ratio);
   _Cijkl(1, 1, 1, 1) = _Cijkl(0, 0, 0, 0);
   _Cijkl(0, 0, 1, 1) = _Cijkl(0, 0, 0, 0) * _poissons_ratio;
   _Cijkl(1, 1, 0, 0) = _Cijkl(0, 0, 1, 1);
   _Cijkl(0, 0, 2, 2) = 0.0;
   _Cijkl(1, 1, 2, 2) = 0.0;
   _Cijkl(2, 2, 2, 2) = 0.0;
   _Cijkl(2, 2, 0, 0) = 0.0;
   _Cijkl(2, 2, 1, 1) = 0.0;
  
   // get number of quadrature points along thickness based on order
   std::unique_ptr<QGauss> t_qrule = libmesh_make_unique<QGauss>(
       1, Utility::string_to_enum<Order>(getParam<std::string>("through_thickness_order")));
   _t_points = t_qrule->get_points();
   _elasticity_tensor.resize(_t_points.size());
   _ge.resize(_t_points.size());
   for (unsigned int t = 0; t < _t_points.size(); ++t)
   {
     _elasticity_tensor[t] =
         &declareADProperty<RankFourTensor>("elasticity_tensor_t_points_" + std::to_string(t));
     _ge[t] = &getADMaterialProperty<RankTwoTensor>("ge_t_points_" + std::to_string(t));
   }
 }

Member Function Documentation

◆ computeQpProperties()

template<ComputeStage compute_stage>

void ADComputeIsotropicElasticityTensorShell< compute_stage >::computeQpProperties ( )

overrideprotectedvirtual

Definition at line 71 of file ADComputeIsotropicElasticityTensorShell.C.

 {
   for (unsigned int t = 0; t < _t_points.size(); ++t)
   {
     (*_elasticity_tensor[t])[_qp].zero();
     // compute contravariant elasticity tensor
     for (unsigned int i = 0; i < 3; ++i)
       for (unsigned int j = 0; j < 3; ++j)
         for (unsigned int k = 0; k < 3; ++k)
           for (unsigned int l = 0; l < 3; ++l)
             for (unsigned int m = 0; m < 3; ++m)
               for (unsigned int n = 0; n < 3; ++n)
                 for (unsigned int o = 0; o < 3; ++o)
                   for (unsigned int p = 0; p < 3; ++p)
                     (*_elasticity_tensor[t])[_qp](i, j, k, l) +=
                         (*_ge[t])[_qp](i, m) * (*_ge[t])[_qp](j, n) * (*_ge[t])[_qp](k, o) *
                         (*_ge[t])[_qp](l, p) * _Cijkl(m, n, o, p);
   }
 }