LinearIsotropic material for use in simple applications that don't need material properties. More...

#include <LinearAnisotropicMaterial.h>

Inheritance diagram for LinearAnisotropicMaterial:

Public Member Functions
	LinearAnisotropicMaterial (const InputParameters &parameters)

Protected Member Functions
virtual void	computeProperties ()

virtual void	computeStress (const SymmTensor &strain, SymmTensor &stress)

virtual void	computeStrain (const SymmTensor &total_strain, SymmTensor &elastic_strain)
	Will always be passed to full symmetric strain tensor. More...

virtual Real	computeAlpha ()

virtual void	initQpStatefulProperties () override

template<typename T >
MaterialProperty< T > &	createProperty (const std::string &prop_name)

template<typename T >
const MaterialProperty< T > &	getPropertyOld (const std::string &prop_name)

Protected Attributes
Real	_material_constant_c11

Real	_material_constant_c12

Real	_material_constant_c44

Real	_euler_angle_1

Real	_euler_angle_2

Real	_euler_angle_3

Real	_youngs_modulus

Real	_poissons_ratio

Real	_t_ref

Real	_alpha

SymmElasticityTensor *	_local_elasticity_tensor

Real	_pi

Real	_tol

const std::string	_appended_property_name

const VariableGradient &	_grad_disp_x

const VariableGradient &	_grad_disp_y

const VariableGradient &	_grad_disp_z

bool	_has_temp

const VariableValue &	_temp

bool	_has_c

const VariableValue &	_c

std::vector< VolumetricModel * >	_volumetric_models

MaterialProperty< SymmTensor > &	_stress

MaterialProperty< SymmElasticityTensor > &	_elasticity_tensor

MaterialProperty< SymmElasticityTensor > &	_Jacobian_mult

SymmTensor	_d_strain_dT

MaterialProperty< SymmTensor > &	_d_stress_dT

MaterialProperty< SymmTensor > &	_elastic_strain

Detailed Description

LinearIsotropic material for use in simple applications that don't need material properties.

Definition at line 24 of file LinearAnisotropicMaterial.h.

Constructor & Destructor Documentation

◆ LinearAnisotropicMaterial()

LinearAnisotropicMaterial::LinearAnisotropicMaterial ( const InputParameters & parameters )

Definition at line 30 of file LinearAnisotropicMaterial.C.

   : LinearIsotropicMaterial(parameters),
     _material_constant_c11(getParam<Real>("material_constant_c11")),
     _material_constant_c12(getParam<Real>("material_constant_c12")),
     _material_constant_c44(getParam<Real>("material_constant_c44")),
     _euler_angle_1(getParam<Real>("euler_angle_1")),
     _euler_angle_2(getParam<Real>("euler_angle_2")),
     _euler_angle_3(getParam<Real>("euler_angle_3"))
 {
   SymmAnisotropicElasticityTensor * aniso_elasticity_tensor = new SymmAnisotropicElasticityTensor;
   aniso_elasticity_tensor->setMaterialConstantc11(_material_constant_c11);
   aniso_elasticity_tensor->setMaterialConstantc12(_material_constant_c12);
   aniso_elasticity_tensor->setMaterialConstantc44(_material_constant_c44);
  
   /* AMJ: I believe I have identified a logic bug with using inherited classes (e.g.,
      LinearAnisotropicMaterial and non-zero Euler angles in conjunction with/inheriting from
      LinearIsotropicMaterial.  When using Euler angles = 0.0, no problem occurs.  However,
      using Euler angles != zero causes the _elasticity_tensor[_qp]=*_local_elasticity_tensor
      to rotate with every single quadrature point, every time the Material class is
      computed.  This is due to the _local_elasticity_tensor->calculate(_qp) call.  Because
      we are dereferencing the _local_elasticity_tensor, we are actually changing the
      original elasticity tensor that is supplied (from the input or material class
      construction, etc).  I've attempted to fix this by moving the relevant information
      into a local copy_local_tensor variable, but (perhaps because I don't know enough
      C++) I cannot get this to work in LinearIsotropicMaterial for an arbitrary derived class that
      may use different derived kinds of Symm***ElasticityTensors.  So... good luck, Chuck.
      I might not be able to fix the problem, but I think I've identified it correctly. */
  
   aniso_elasticity_tensor->setFirstEulerAngle(_euler_angle_1);
   aniso_elasticity_tensor->setSecondEulerAngle(_euler_angle_2);
   aniso_elasticity_tensor->setThirdEulerAngle(_euler_angle_3);
  
   delete _local_elasticity_tensor;
   _local_elasticity_tensor = aniso_elasticity_tensor;
 }

Member Function Documentation

◆ computeAlpha()

Real LinearIsotropicMaterial::computeAlpha ( )

protectedvirtualinherited

Definition at line 143 of file LinearIsotropicMaterial.C.

 {
   return _alpha;
 }

Referenced by LinearIsotropicMaterial::computeProperties().

◆ computeProperties()

void LinearIsotropicMaterial::computeProperties ( )

protectedvirtualinherited

Definition at line 54 of file LinearIsotropicMaterial.C.

 {
   for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
   {
     Real alpha = computeAlpha();
  
     _local_elasticity_tensor->calculate(_qp);
  
     _elasticity_tensor[_qp] = *_local_elasticity_tensor;
  
     SymmTensor strn(_grad_disp_x[_qp](0),
                     _grad_disp_y[_qp](1),
                     _grad_disp_z[_qp](2),
                     0.5 * (_grad_disp_x[_qp](1) + _grad_disp_y[_qp](0)),
                     0.5 * (_grad_disp_y[_qp](2) + _grad_disp_z[_qp](1)),
                     0.5 * (_grad_disp_z[_qp](0) + _grad_disp_x[_qp](2)));
  
     // Add in Isotropic Thermal Strain
     if (_has_temp)
     {
       Real isotropic_strain = alpha * (_temp[_qp] - _t_ref);
  
       strn.addDiag(-isotropic_strain);
  
       _d_strain_dT.zero();
       _d_strain_dT.addDiag(-alpha);
     }
  
     SymmTensor strain(strn);
  
     computeStress(strain, _stress[_qp]);
   }
 }

◆ computeStrain()

void LinearIsotropicMaterial::computeStrain	(	const SymmTensor &	total_strain,
		SymmTensor &	elastic_strain
	)

protectedvirtualinherited

Will always be passed to full symmetric strain tensor.

What should come out is a modified strain tensor.

Definition at line 117 of file LinearIsotropicMaterial.C.

 {
  
   // Multiplier that zeros out stiffness
   // Real h = (1.0 - _c[_qp]*_c[_qp]);
   /*Real vl = -0.05 + _c[_qp]*1.1;
   if (vl < 0.0)
     vl = 0.0;
   if (vl > 1.0)
   vl = 1.0;*/
  
   Real h = (1.0 + std::cos(_pi * _c[_qp])) / 2.0;
   if (h < _tol)
     h = _tol;
  
   elastic_strain = total_strain;
   // Jacobian multiplier of the stress
   _Jacobian_mult[_qp] = *_local_elasticity_tensor * h;
  
   SymmTensor d_stress_dT(*_local_elasticity_tensor * _d_strain_dT);
   //   d_stress_dT *= _dt;
   d_stress_dT *= h;
   _d_stress_dT[_qp] = d_stress_dT;
 }

Referenced by LinearIsotropicMaterial::computeStress().

◆ computeStress()

void LinearIsotropicMaterial::computeStress	(	const SymmTensor &	strain,
		SymmTensor &	stress
	)

protectedvirtualinherited

Definition at line 89 of file LinearIsotropicMaterial.C.

 {
   // Add in any extra strain components
   SymmTensor elastic_strain;
  
   computeStrain(strain, elastic_strain);
  
   // Multiplier that zeros out stiffness
   // Real h = (1.0 - _c[_qp]*_c[_qp]);
   /*Real vl = -0.05 + _c[_qp]*1.1;
   if (vl < 0.0)
     vl = 0.0;
   if (vl > 1.0)
   vl = 1.0;*/
  
   Real h = (1.0 + std::cos(_pi * _c[_qp])) / 2.0;
   if (h < _tol)
     h = _tol;
  
   // Save that off as the elastic strain
   _elastic_strain[_qp] = elastic_strain;
  
   // Create column vector
   // C * e
   stress = (*_local_elasticity_tensor) * h * elastic_strain;
 }

Referenced by LinearIsotropicMaterial::computeProperties().

◆ createProperty()

template<typename T >

MaterialProperty<T>& SolidMechanicsMaterial::createProperty ( const std::string & prop_name )

inlineprotectedinherited

Definition at line 58 of file SolidMechanicsMaterial.h.

   {
     std::string name(prop_name + _appended_property_name);
     return declareProperty<T>(name);
   }

◆ getPropertyOld()

template<typename T >

const MaterialProperty<T>& SolidMechanicsMaterial::getPropertyOld ( const std::string & prop_name )

inlineprotectedinherited

Definition at line 65 of file SolidMechanicsMaterial.h.

   {
     std::string name(prop_name + _appended_property_name);
     return getMaterialPropertyOld<T>(name);
   }

◆ initQpStatefulProperties()

void SolidMechanicsMaterial::initQpStatefulProperties ( )

overrideprotectedvirtualinherited

Definition at line 32 of file SolidMechanicsMaterial.C.

33 {

34 }

Member Data Documentation

◆ _alpha

Real LinearIsotropicMaterial::_alpha

protectedinherited

Definition at line 48 of file LinearIsotropicMaterial.h.

Referenced by LinearIsotropicMaterial::computeAlpha().

◆ _appended_property_name

const std::string SolidMechanicsMaterial::_appended_property_name

protectedinherited

Definition at line 32 of file SolidMechanicsMaterial.h.

Referenced by SolidMechanicsMaterial::createProperty(), and SolidMechanicsMaterial::getPropertyOld().

◆ _c

const VariableValue& SolidMechanicsMaterial::_c

protectedinherited

Definition at line 41 of file SolidMechanicsMaterial.h.

Referenced by LinearIsotropicMaterial::computeStrain(), and LinearIsotropicMaterial::computeStress().

◆ _d_strain_dT

SymmTensor SolidMechanicsMaterial::_d_strain_dT

protectedinherited

Definition at line 50 of file SolidMechanicsMaterial.h.

Referenced by LinearIsotropicMaterial::computeProperties(), and LinearIsotropicMaterial::computeStrain().

◆ _d_stress_dT

MaterialProperty<SymmTensor>& SolidMechanicsMaterial::_d_stress_dT

protectedinherited

Definition at line 53 of file SolidMechanicsMaterial.h.

Referenced by LinearIsotropicMaterial::computeStrain().

◆ _elastic_strain

MaterialProperty<SymmTensor>& SolidMechanicsMaterial::_elastic_strain

protectedinherited

Definition at line 55 of file SolidMechanicsMaterial.h.

Referenced by LinearGeneralAnisotropicMaterial::computeQpStrain(), LinearGeneralAnisotropicMaterial::computeQpStress(), and LinearIsotropicMaterial::computeStress().

◆ _elasticity_tensor

MaterialProperty<SymmElasticityTensor>& SolidMechanicsMaterial::_elasticity_tensor

protectedinherited

Definition at line 46 of file SolidMechanicsMaterial.h.

Referenced by LinearIsotropicMaterial::computeProperties(), LinearGeneralAnisotropicMaterial::computeQpElasticityTensor(), and LinearGeneralAnisotropicMaterial::computeQpStress().

◆ _euler_angle_1

Real LinearAnisotropicMaterial::_euler_angle_1

protected

Definition at line 33 of file LinearAnisotropicMaterial.h.

Referenced by LinearAnisotropicMaterial().

◆ _euler_angle_2

Real LinearAnisotropicMaterial::_euler_angle_2

protected

Definition at line 34 of file LinearAnisotropicMaterial.h.

Referenced by LinearAnisotropicMaterial().

◆ _euler_angle_3

Real LinearAnisotropicMaterial::_euler_angle_3

protected

Definition at line 35 of file LinearAnisotropicMaterial.h.

Referenced by LinearAnisotropicMaterial().

◆ _grad_disp_x

const VariableGradient& SolidMechanicsMaterial::_grad_disp_x

protectedinherited

Definition at line 33 of file SolidMechanicsMaterial.h.

Referenced by LinearIsotropicMaterial::computeProperties(), and LinearGeneralAnisotropicMaterial::computeQpStrain().

◆ _grad_disp_y

const VariableGradient& SolidMechanicsMaterial::_grad_disp_y

protectedinherited

Definition at line 34 of file SolidMechanicsMaterial.h.

Referenced by LinearIsotropicMaterial::computeProperties(), and LinearGeneralAnisotropicMaterial::computeQpStrain().

◆ _grad_disp_z

const VariableGradient& SolidMechanicsMaterial::_grad_disp_z

protectedinherited

Definition at line 35 of file SolidMechanicsMaterial.h.

Referenced by LinearIsotropicMaterial::computeProperties(), and LinearGeneralAnisotropicMaterial::computeQpStrain().

◆ _has_c

bool SolidMechanicsMaterial::_has_c

protectedinherited

Definition at line 40 of file SolidMechanicsMaterial.h.

◆ _has_temp

bool SolidMechanicsMaterial::_has_temp

protectedinherited

Definition at line 37 of file SolidMechanicsMaterial.h.

Referenced by LinearIsotropicMaterial::computeProperties().

◆ _Jacobian_mult

MaterialProperty<SymmElasticityTensor>& SolidMechanicsMaterial::_Jacobian_mult

protectedinherited

Definition at line 47 of file SolidMechanicsMaterial.h.

Referenced by LinearGeneralAnisotropicMaterial::computeQpElasticityTensor(), and LinearIsotropicMaterial::computeStrain().

◆ _local_elasticity_tensor

SymmElasticityTensor* LinearIsotropicMaterial::_local_elasticity_tensor

protectedinherited

Definition at line 50 of file LinearIsotropicMaterial.h.

Referenced by LinearIsotropicMaterial::computeProperties(), LinearIsotropicMaterial::computeStrain(), LinearAnisotropicMaterial(), LinearIsotropicMaterial::LinearIsotropicMaterial(), and LinearIsotropicMaterial::~LinearIsotropicMaterial().

◆ _material_constant_c11

Real LinearAnisotropicMaterial::_material_constant_c11

protected

Definition at line 30 of file LinearAnisotropicMaterial.h.

Referenced by LinearAnisotropicMaterial().

◆ _material_constant_c12

Real LinearAnisotropicMaterial::_material_constant_c12

protected

Definition at line 31 of file LinearAnisotropicMaterial.h.

Referenced by LinearAnisotropicMaterial().

◆ _material_constant_c44

Real LinearAnisotropicMaterial::_material_constant_c44

protected

Definition at line 32 of file LinearAnisotropicMaterial.h.

Referenced by LinearAnisotropicMaterial().

◆ _pi

Real LinearIsotropicMaterial::_pi

protectedinherited

Definition at line 52 of file LinearIsotropicMaterial.h.

Referenced by LinearIsotropicMaterial::computeStrain(), LinearIsotropicMaterial::computeStress(), and LinearIsotropicMaterial::LinearIsotropicMaterial().

◆ _poissons_ratio

Real LinearIsotropicMaterial::_poissons_ratio

protectedinherited

Definition at line 45 of file LinearIsotropicMaterial.h.

Referenced by LinearIsotropicMaterial::LinearIsotropicMaterial().

◆ _stress

MaterialProperty<SymmTensor>& SolidMechanicsMaterial::_stress

protectedinherited

Definition at line 45 of file SolidMechanicsMaterial.h.

Referenced by LinearIsotropicMaterial::computeProperties(), and LinearGeneralAnisotropicMaterial::computeQpStress().

◆ _t_ref

Real LinearIsotropicMaterial::_t_ref

protectedinherited

Definition at line 47 of file LinearIsotropicMaterial.h.

Referenced by LinearIsotropicMaterial::computeProperties().

◆ _temp

const VariableValue& SolidMechanicsMaterial::_temp

protectedinherited

Definition at line 38 of file SolidMechanicsMaterial.h.

Referenced by LinearIsotropicMaterial::computeProperties().

◆ _tol

Real LinearIsotropicMaterial::_tol

protectedinherited

Definition at line 53 of file LinearIsotropicMaterial.h.

Referenced by LinearIsotropicMaterial::computeStrain(), LinearIsotropicMaterial::computeStress(), and LinearIsotropicMaterial::LinearIsotropicMaterial().

◆ _volumetric_models

std::vector<VolumetricModel *> SolidMechanicsMaterial::_volumetric_models

protectedinherited

Definition at line 43 of file SolidMechanicsMaterial.h.

◆ _youngs_modulus

Real LinearIsotropicMaterial::_youngs_modulus

protectedinherited

Definition at line 44 of file LinearIsotropicMaterial.h.

Referenced by LinearIsotropicMaterial::LinearIsotropicMaterial().

The documentation for this class was generated from the following files:

solid_mechanics/include/materials/total_strain/LinearAnisotropicMaterial.h
solid_mechanics/src/materials/total_strain/LinearAnisotropicMaterial.C

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ LinearAnisotropicMaterial()

Member Function Documentation

◆ computeAlpha()

◆ computeProperties()

◆ computeStrain()

◆ computeStress()

◆ createProperty()

◆ getPropertyOld()

◆ initQpStatefulProperties()

Member Data Documentation

◆ _alpha

◆ _appended_property_name

◆ _c

◆ _d_strain_dT

◆ _d_stress_dT

◆ _elastic_strain

◆ _elasticity_tensor

◆ _euler_angle_1

◆ _euler_angle_2

◆ _euler_angle_3

◆ _grad_disp_x

◆ _grad_disp_y

◆ _grad_disp_z

◆ _has_c

◆ _has_temp

◆ _Jacobian_mult

◆ _local_elasticity_tensor

◆ _material_constant_c11

◆ _material_constant_c12

◆ _material_constant_c44

◆ _pi

◆ _poissons_ratio

◆ _stress

◆ _t_ref

◆ _temp

◆ _tol

◆ _volumetric_models

◆ _youngs_modulus