SymmIsotropicElasticityTensor Class Reference

Defines an Isotropic Elasticity Tensor. More...

#include <SymmIsotropicElasticityTensor.h>

Inheritance diagram for SymmIsotropicElasticityTensor:

Public Member Functions
	SymmIsotropicElasticityTensor (const bool constant=true)
virtual	~SymmIsotropicElasticityTensor ()
void	unsetConstants ()
void	setLambda (const Real lambda)
	Set the first Lame Coefficient. More...
void	setMu (const Real mu)
	Set the second Lame Coefficient. More...
void	setYoungsModulus (const Real E)
	Set the Young's Modulus. More...
void	setPoissonsRatio (const Real nu)
	Set Poissons Ratio. More...
void	setBulkModulus (const Real k)
	Set the Bulk Modulus. More...
void	setShearModulus (const Real k)
	Set the shear modulus... More...
Real	mu () const
	Return Mu. More...
Real	shearModulus () const
	Return the shear modulus... More...
Real	youngsModulus () const
	Return the youngs modulus. More...
virtual Real	stiffness (const unsigned i, const unsigned j, const RealGradient &test, const RealGradient &phi) const
virtual void	multiply (const SymmTensor &x, SymmTensor &b) const
virtual void	adjustForCracking (const RealVectorValue &crack_flags)
virtual void	adjustForCrackingWithShearRetention (const RealVectorValue &crack_flags)
void	constant (bool c)
void	copyValues (SymmElasticityTensor &rhs) const
void	calculate (unsigned int qp)
	Public function that will be called whenever the values for this matrix need to be filled in. More...
SymmTensor	operator* (const SymmTensor &x) const
SymmElasticityTensor	operator* (Real x) const
virtual Real	stiffness (const unsigned int i, const unsigned int j, const RealGradient &test, const RealGradient &phi) const
SymmElasticityTensor	operator+ (const SymmElasticityTensor &rhs) const
void	operator+= (const SymmElasticityTensor &rhs)
void	operator-= (const SymmElasticityTensor &rhs)
void	operator*= (Real rhs)
void	operator/= (Real rhs)
void	zero ()
void	convertFrom9x9 (const ColumnMajorMatrix &cmm)
void	convertFrom6x6 (const ColumnMajorMatrix &cmm)
ColumnMajorMatrix	columnMajorMatrix9x9 () const
ColumnMajorMatrix	columnMajorMatrix6x6 () const
void	form9x9Rotation (const ColumnMajorMatrix &R_3x3, ColumnMajorMatrix &R_9x9) const
void	rotateFromGlobalToLocal (const ColumnMajorMatrix &R)
void	rotateFromLocalToGlobal (const ColumnMajorMatrix &R)
virtual SymmElasticityTensor	calculateDerivative (unsigned int qp, unsigned int i)
void	fillFromInputVector (std::vector< Real > input, bool all)
Real	sum_3x3 () const
RealGradient	sum_3x1 () const
Real	valueAtIndex (int i) const

Protected Member Functions
virtual void	calculateEntries (unsigned int qp)
	Fill in the matrix. More...
void	setEntries (Real C11, Real C12, Real C44)
void	calculateLameCoefficients ()
	Calculates lambda and mu based on what has been set. More...

Protected Attributes
bool	_lambda_set
bool	_mu_set
bool	_E_set
bool	_nu_set
bool	_k_set
Real	_lambda
Real	_mu
Real	_E
Real	_nu
Real	_k
bool	_constant
	Whether or not the matrix is constant for all of time and space. More...
bool	_values_computed
	Whether or not the values have been computed once. More...
Real	_val [21]

Detailed Description

Defines an Isotropic Elasticity Tensor.

The input is any two of the following:

Youngs Modulus Poissons Ration First and Second Lame Coefficients (lambda and mu) Bulk Modulus

Internally this class uses the Lame Coefficients. Everything is is transformed to these coefficients.

Note that by default this tensor is constant... meaning that it never changes after the first time it is computed.

If you want to modify this behavior you can pass in false to the constructor.

Definition at line 33 of file SymmIsotropicElasticityTensor.h.

Constructor & Destructor Documentation

SymmIsotropicElasticityTensor()

SymmIsotropicElasticityTensor::SymmIsotropicElasticityTensor

(

const bool

constant = true

)

Definition at line 12 of file SymmIsotropicElasticityTensor.C.

  : SymmElasticityTensor(constant),
    _lambda_set(false),
    _mu_set(false),
    _E_set(false),
    _nu_set(false),
    _k_set(false),
    _lambda(0),
    _mu(0),
    _E(0),
    _nu(0),
    _k(0)
{
}

~SymmIsotropicElasticityTensor()

virtual SymmIsotropicElasticityTensor::~SymmIsotropicElasticityTensor ( )

inlinevirtual

Definition at line 38 of file SymmIsotropicElasticityTensor.h.

{}

Member Function Documentation

adjustForCracking()

void SymmIsotropicElasticityTensor::adjustForCracking ( const RealVectorValue &  crack_flags )

virtual

Reimplemented from SymmElasticityTensor.

Definition at line 251 of file SymmIsotropicElasticityTensor.C.

{
  const RealVectorValue & c(crack_flags);
  const Real c0(c(0));
  const Real c0_coupled(c0 < 1 ? 0 : 1);
  const Real c1(c(1));
  const Real c1_coupled(c1 < 1 ? 0 : 1);
  const Real c2(c(2));
  const Real c2_coupled(c2 < 1 ? 0 : 1);
 
  const Real c01(c0_coupled * c1_coupled);
  const Real c02(c0_coupled * c2_coupled);
  const Real c12(c1_coupled * c2_coupled);
  const Real c012(c0_coupled * c12);
 
  const Real ym = _mu * (3 * _lambda + 2 * _mu) / (_lambda + _mu);
 
  // Assume Poisson's ratio goes to zero for the cracked direction.
 
  _val[0] = (c0 < 1 ? c0 * ym : _val[0]);
  _val[1] *= c01;
  _val[2] *= c02;
  _val[3] *= c01;
  _val[4] *= c012;
  _val[5] *= c02;
 
  _val[6] = (c1 < 1 ? c1 * ym : _val[6]);
  _val[7] *= c12;
  _val[8] *= c01;
  _val[9] *= c12;
  _val[10] *= c012;
 
  _val[11] = (c2 < 1 ? c2 * ym : _val[11]);
  _val[12] *= c012;
  _val[13] *= c12;
  _val[14] *= c02;
}

Referenced by adjustForCrackingWithShearRetention().

adjustForCrackingWithShearRetention()

void SymmIsotropicElasticityTensor::adjustForCrackingWithShearRetention ( const RealVectorValue &  crack_flags )

virtual

Reimplemented from SymmElasticityTensor.

Definition at line 290 of file SymmIsotropicElasticityTensor.C.

{
  const RealVectorValue & c = crack_flags;
  const Real c0 = c(0);
  const Real c0_coupled = (c0 < 1 ? 0 : 1);
  const Real c1 = c(1);
  const Real c1_coupled = (c1 < 1 ? 0 : 1);
  const Real c2 = c(2);
  const Real c2_coupled = (c2 < 1 ? 0 : 1);
  const Real c01 = c0_coupled * c1_coupled;
  const Real c02 = c0_coupled * c2_coupled;
  const Real c12 = c1_coupled * c2_coupled;
  const Real c012 = c0_coupled * c12;
  adjustForCracking(crack_flags);
  _val[15] *= c01;
  _val[16] *= c012;
  _val[17] *= c012;
  _val[18] *= c12;
  _val[19] *= c012;
  _val[20] *= c02;
}

calculate()

void SymmElasticityTensor::calculate ( unsigned int  qp )

inherited

Public function that will be called whenever the values for this matrix need to be filled in.

Definition at line 41 of file SymmElasticityTensor.C.

{
  if (!_constant || !_values_computed)
  {
    calculateEntries(qp);
    _values_computed = true;
  }
}

Referenced by SolidModel::computeElasticityTensor(), LinearIsotropicMaterial::computeProperties(), and SolidModel::createElasticityTensor().

calculateDerivative()

SymmElasticityTensor SymmElasticityTensor::calculateDerivative ( unsigned int  qp, unsigned int  i  )

virtualinherited

Definition at line 328 of file SymmElasticityTensor.C.

{
  return SymmElasticityTensor();
}

calculateEntries()

void SymmIsotropicElasticityTensor::calculateEntries ( unsigned int  qp )

protectedvirtual

Fill in the matrix.

Reimplemented from SymmElasticityTensor.

Definition at line 131 of file SymmIsotropicElasticityTensor.C.

{
  calculateLameCoefficients();
 
  const Real C12(_lambda);
  const Real C44(_mu);
  const Real C11(2 * C44 + C12);
 
  setEntries(C11, C12, C44);
}

calculateLameCoefficients()

void SymmIsotropicElasticityTensor::calculateLameCoefficients ( )

protected

Calculates lambda and mu based on what has been set.

These are based on Michael Tonks's's notes

Definition at line 94 of file SymmIsotropicElasticityTensor.C.

{
  if (_lambda_set && _mu_set) // First and second Lame
    return;
  else if (_lambda_set && _nu_set)
    _mu = (_lambda * (1.0 - 2.0 * _nu)) / (2.0 * _nu);
  else if (_lambda_set && _k_set)
    _mu = (3.0 * (_k - _lambda)) / 2.0;
  else if (_lambda_set && _E_set)
    _mu = ((_E - 3.0 * _lambda) / 4.0) +
          (std::sqrt((_E - 3.0 * _lambda) * (_E - 3.0 * _lambda) + 8.0 * _lambda * _E) / 4.0);
  else if (_mu_set && _nu_set)
    _lambda = (2.0 * _mu * _nu) / (1.0 - 2.0 * _nu);
  else if (_mu_set && _k_set)
    _lambda = (3.0 * _k - 2.0 * _mu) / 3.0;
  else if (_mu_set && _E_set)
    _lambda = ((2.0 * _mu - _E) * _mu) / (_E - 3.0 * _mu);
  else if (_nu_set && _k_set)
  {
    _lambda = (3.0 * _k * _nu) / (1.0 + _nu);
    _mu = (3.0 * _k * (1.0 - 2.0 * _nu)) / (2.0 * (1.0 + _nu));
  }
  else if (_E_set && _nu_set) // Young's Modulus and Poisson's Ratio
  {
    _lambda = (_nu * _E) / ((1.0 + _nu) * (1 - 2.0 * _nu));
    _mu = _E / (2.0 * (1.0 + _nu));
  }
  else if (_E_set && _k_set)
  {
    _lambda = (3.0 * _k * (3.0 * _k - _E)) / (9.0 * _k - _E);
    _mu = (3.0 * _E * _k) / (9.0 * _k - _E);
  }
  _lambda_set = true;
  _mu_set = true;
}

Referenced by calculateEntries().

columnMajorMatrix6x6()

ColumnMajorMatrix SymmElasticityTensor::columnMajorMatrix6x6 ( ) const

inherited

Definition at line 230 of file SymmElasticityTensor.C.

{
  ColumnMajorMatrix cmm(6, 6);
  unsigned count(0);
  for (unsigned i(0); i < 6; ++i)
  {
    for (unsigned j(i); j < 6; ++j)
    {
      cmm(i, j) = cmm(j, i) = _val[count++];
    }
  }
  return cmm;
}

Referenced by SymmAnisotropicElasticityTensor::form_transformed_material_dt_matrix().

columnMajorMatrix9x9()

ColumnMajorMatrix SymmElasticityTensor::columnMajorMatrix9x9 ( ) const

inherited

Definition at line 245 of file SymmElasticityTensor.C.

{
  ColumnMajorMatrix cmm(9, 9);
  cmm(0, 0) = _val[0];
  cmm(0, 1) = cmm(1, 0) = _val[3];
  cmm(0, 2) = cmm(2, 0) = _val[5];
  cmm(0, 3) = cmm(3, 0) = _val[3];
  cmm(0, 4) = cmm(4, 0) = _val[1];
  cmm(0, 5) = cmm(5, 0) = _val[4];
  cmm(0, 6) = cmm(6, 0) = _val[5];
  cmm(0, 7) = cmm(7, 0) = _val[4];
  cmm(0, 8) = cmm(8, 0) = _val[2];
 
  cmm(1, 1) = _val[15];
  cmm(1, 2) = cmm(2, 1) = _val[17];
  cmm(1, 3) = cmm(3, 1) = _val[15];
  cmm(1, 4) = cmm(4, 1) = _val[8];
  cmm(1, 5) = cmm(5, 1) = _val[16];
  cmm(1, 6) = cmm(6, 1) = _val[17];
  cmm(1, 7) = cmm(7, 1) = _val[16];
  cmm(1, 8) = cmm(8, 1) = _val[12];
 
  cmm(2, 2) = _val[20];
  cmm(2, 3) = cmm(3, 2) = _val[17];
  cmm(2, 4) = cmm(4, 2) = _val[10];
  cmm(2, 5) = cmm(5, 2) = _val[19];
  cmm(2, 6) = cmm(6, 2) = _val[20];
  cmm(2, 7) = cmm(7, 2) = _val[19];
  cmm(2, 8) = cmm(8, 2) = _val[14];
 
  cmm(3, 3) = _val[15];
  cmm(3, 4) = cmm(4, 3) = _val[8];
  cmm(3, 5) = cmm(5, 3) = _val[16];
  cmm(3, 6) = cmm(6, 3) = _val[17];
  cmm(3, 7) = cmm(7, 3) = _val[16];
  cmm(3, 8) = cmm(8, 3) = _val[12];
 
  cmm(4, 4) = _val[6];
  cmm(4, 5) = cmm(5, 4) = _val[9];
  cmm(4, 6) = cmm(6, 4) = _val[10];
  cmm(4, 7) = cmm(7, 4) = _val[9];
  cmm(4, 8) = cmm(8, 4) = _val[7];
 
  cmm(5, 5) = _val[18];
  cmm(5, 6) = cmm(6, 5) = _val[19];
  cmm(5, 7) = cmm(7, 5) = _val[18];
  cmm(5, 8) = cmm(8, 5) = _val[13];
 
  cmm(6, 6) = _val[20];
  cmm(6, 7) = cmm(7, 6) = _val[19];
  cmm(6, 8) = cmm(8, 6) = _val[14];
 
  cmm(7, 7) = _val[18];
  cmm(7, 8) = cmm(8, 7) = _val[13];
 
  cmm(8, 8) = _val[11];
 
  return cmm;
}

Referenced by SymmAnisotropicElasticityTensor::form_transformed_material_dmat_matrix(), SymmElasticityTensor::rotateFromGlobalToLocal(), and SymmElasticityTensor::rotateFromLocalToGlobal().

constant()

void SymmElasticityTensor::constant ( bool  c )

inlineinherited

Definition at line 66 of file SymmElasticityTensor.h.

{ _constant = c; }

convertFrom6x6()

void SymmElasticityTensor::convertFrom6x6 ( const ColumnMajorMatrix &  cmm )

inherited

Definition at line 194 of file SymmElasticityTensor.C.

{
  if (input.numEntries() != 36)
  {
    mooseError("Cannot convert from ColumnMajorMatrix (wrong size)");
  }
 
  _val[0] = input(0, 0);
  _val[1] = input(0, 1);
  _val[2] = input(0, 2);
  _val[3] = input(0, 3);
  _val[4] = input(0, 4);
  _val[5] = input(0, 5);
 
  _val[6] = input(1, 1);
  _val[7] = input(1, 2);
  _val[8] = input(1, 3);
  _val[9] = input(1, 4);
  _val[10] = input(1, 5);
 
  _val[11] = input(2, 2);
  _val[12] = input(2, 3);
  _val[13] = input(2, 4);
  _val[14] = input(2, 5);
 
  _val[15] = input(3, 3);
  _val[16] = input(3, 4);
  _val[17] = input(3, 5);
 
  _val[18] = input(4, 4);
  _val[19] = input(4, 5);
 
  _val[20] = input(5, 5);
}

convertFrom9x9()

void SymmElasticityTensor::convertFrom9x9 ( const ColumnMajorMatrix &  cmm )

inherited

Definition at line 158 of file SymmElasticityTensor.C.

{
  if (input.numEntries() != 81)
  {
    mooseError("Cannot convert from ColumnMajorMatrix (wrong size)");
  }
 
  _val[0] = input(0, 0);
  _val[1] = input(0, 4);
  _val[2] = input(0, 8);
  _val[3] = input(0, 1);
  _val[4] = input(0, 5);
  _val[5] = input(0, 2);
 
  _val[6] = input(4, 4);
  _val[7] = input(4, 8);
  _val[8] = input(4, 3);
  _val[9] = input(4, 5);
  _val[10] = input(4, 6);
 
  _val[11] = input(8, 8);
  _val[12] = input(8, 3);
  _val[13] = input(8, 5);
  _val[14] = input(8, 6);
 
  _val[15] = input(1, 1);
  _val[16] = input(1, 5);
  _val[17] = input(1, 2);
 
  _val[18] = input(5, 5);
  _val[19] = input(5, 6);
 
  _val[20] = input(2, 2);
}

Referenced by SymmAnisotropicElasticityTensor::form_transformed_material_dt_matrix(), SymmElasticityTensor::rotateFromGlobalToLocal(), and SymmElasticityTensor::rotateFromLocalToGlobal().

copyValues()

void SymmElasticityTensor::copyValues ( SymmElasticityTensor &  rhs ) const

inlineinherited

Definition at line 70 of file SymmElasticityTensor.h.

  {
    for (unsigned i(0); i < 21; ++i)
    {
      rhs._val[i] = _val[i];
    }
  }

Referenced by SymmAnisotropicElasticityTensor::form_transformed_material_dmat_matrix().

fillFromInputVector()

void SymmElasticityTensor::fillFromInputVector ( std::vector< Real >  input, bool  all  )

inherited

Definition at line 384 of file SymmElasticityTensor.C.

{
  if ((all == true && input.size() != 21) || (all == false && input.size() != 9))
    mooseError("Please check the number of entries in the stiffness input vector.");
 
  if (all == true)
  {
    for (int i = 0; i < 21; i++)
      _val[i] = input[i];
  }
  else
  {
    _val[0] = input[0];  // C1111
    _val[1] = input[1];  // C1122
    _val[2] = input[2];  // C1133
    _val[6] = input[3];  // C2222
    _val[7] = input[4];  // C2233
    _val[11] = input[5]; // C3333
    _val[15] = input[6]; // C2323
    _val[18] = input[7]; // C1313
    _val[20] = input[8]; // C1212
  }
}

Referenced by LinearGeneralAnisotropicMaterial::LinearGeneralAnisotropicMaterial(), and MacroElastic::updateElasticityTensor().

form9x9Rotation()

void SymmElasticityTensor::form9x9Rotation ( const ColumnMajorMatrix &  R_3x3, ColumnMajorMatrix &  R_9x9  ) const

inherited

Definition at line 341 of file SymmElasticityTensor.C.

{
  for (int i = 0; i < 3; ++i)
  {
    for (int j = 0; j < 3; ++j)
    {
      for (int k = 0; k < 3; ++k)
      {
        for (int l = 0; l < 3; ++l)
        {
          R_9x9(((i * 3) + k), ((j * 3) + l)) = R_3x3(i, j) * R_3x3(k, l);
        }
      }
    }
  }
}

Referenced by SolidModel::crackingStrainDirections().

mu()

Real SymmIsotropicElasticityTensor::mu

(

)

const

Return Mu.

Definition at line 84 of file SymmIsotropicElasticityTensor.C.

{
  if (!_mu_set)
  {
    mooseError("mu not set");
  }
  return _mu;
}

Referenced by setMu(), and shearModulus().

multiply()

void SymmIsotropicElasticityTensor::multiply ( const SymmTensor &  x, SymmTensor &  b  ) const

virtual

Reimplemented from SymmElasticityTensor.

Definition at line 224 of file SymmIsotropicElasticityTensor.C.

{
  const Real xx = x.xx();
  const Real yy = x.yy();
  const Real zz = x.zz();
  const Real xy = x.xy();
  const Real yz = x.yz();
  const Real zx = x.zx();
 
  b.xx() = _val[0] * xx + _val[1] * yy + _val[2] * zz;
  b.yy() = _val[1] * xx + _val[6] * yy + _val[7] * zz;
  b.zz() = _val[2] * xx + _val[7] * yy + _val[11] * zz;
  b.xy() = 2 * _val[15] * xy;
  b.yz() = 2 * _val[18] * yz;
  b.zx() = 2 * _val[20] * zx;
 
  b.xx() += 2 * (_val[3] * xy + _val[4] * yz + _val[5] * zx);
  b.yy() += 2 * (_val[8] * xy + _val[9] * yz + _val[10] * zx);
  b.zz() += 2 * (_val[12] * xy + _val[13] * yz + _val[14] * zx);
  b.xy() += _val[3] * xx + _val[8] * yy + _val[12] * zz;
  b.yz() += _val[4] * xx + _val[9] * yy + _val[13] * zz;
  b.zx() += _val[5] * xx + _val[10] * yy + _val[14] * zz;
  b.yz() += 2 * _val[16] * xy;
  b.zx() += 2 * _val[17] * xy + 2 * _val[19] * yz;
}

operator*() [1/2]

SymmTensor SymmElasticityTensor::operator* ( const SymmTensor &  x ) const

inherited

Definition at line 79 of file SymmElasticityTensor.C.

{
  SymmTensor b;
  multiply(x, b);
  return b;
}

operator*() [2/2]

SymmElasticityTensor SymmElasticityTensor::operator* ( Real  x ) const

inherited

Definition at line 333 of file SymmElasticityTensor.C.

{
  SymmElasticityTensor fred(*this);
  fred *= x;
  return fred;
}

operator*=()

void SymmElasticityTensor::operator*= ( Real  rhs )

inlineinherited

Definition at line 118 of file SymmElasticityTensor.h.

  {
    for (unsigned i(0); i < 21; ++i)
    {
      _val[i] *= rhs;
    }
  }

operator+()

SymmElasticityTensor SymmElasticityTensor::operator+ ( const SymmElasticityTensor &  rhs ) const

inlineinherited

Definition at line 93 of file SymmElasticityTensor.h.

  {
    SymmElasticityTensor t = *this;
 
    t += rhs;
 
    return t;
  }

operator+=()

void SymmElasticityTensor::operator+= ( const SymmElasticityTensor &  rhs )

inlineinherited

Definition at line 102 of file SymmElasticityTensor.h.

  {
    for (unsigned i(0); i < 21; ++i)
    {
      _val[i] += rhs._val[i];
    }
  }

operator-=()

void SymmElasticityTensor::operator-= ( const SymmElasticityTensor &  rhs )

inlineinherited

Definition at line 110 of file SymmElasticityTensor.h.

  {
    for (unsigned i(0); i < 21; ++i)
    {
      _val[i] -= rhs._val[i];
    }
  }

operator/=()

void SymmElasticityTensor::operator/= ( Real  rhs )

inlineinherited

Definition at line 126 of file SymmElasticityTensor.h.

  {
    for (unsigned i(0); i < 21; ++i)
    {
      _val[i] /= rhs;
    }
  }

rotateFromGlobalToLocal()

void SymmElasticityTensor::rotateFromGlobalToLocal ( const ColumnMajorMatrix &  R )

inherited

Definition at line 366 of file SymmElasticityTensor.C.

{
  convertFrom9x9(R.transpose() * (columnMajorMatrix9x9() * R));
}

rotateFromLocalToGlobal()

void SymmElasticityTensor::rotateFromLocalToGlobal ( const ColumnMajorMatrix &  R )

inherited

Definition at line 360 of file SymmElasticityTensor.C.

{
  convertFrom9x9((R * columnMajorMatrix9x9()) * R.transpose());
}

Referenced by SolidModel::crackingStrainDirections().

setBulkModulus()

void SymmIsotropicElasticityTensor::setBulkModulus

(

const Real

k

)

Set the Bulk Modulus.

Definition at line 56 of file SymmIsotropicElasticityTensor.C.

{
  _k = k;
  _k_set = true;
}

setEntries()

void SymmIsotropicElasticityTensor::setEntries ( Real  C11, Real  C12, Real  C44  )

protected

Definition at line 143 of file SymmIsotropicElasticityTensor.C.

{
  _val[0] = _val[6] = _val[11] = C11;
  _val[1] = _val[2] = _val[7] = C12;
  _val[15] = _val[18] = _val[20] = C44;
  _val[3] = _val[4] = _val[5] = 0;
  _val[8] = _val[9] = _val[10] = 0;
  _val[12] = _val[13] = _val[14] = 0;
  _val[16] = _val[17] = 0;
  _val[19] = 0;
}

Referenced by calculateEntries().

setLambda()

void SymmIsotropicElasticityTensor::setLambda

(

const Real

lambda

)

Set the first Lame Coefficient.

Definition at line 28 of file SymmIsotropicElasticityTensor.C.

{
  _lambda = lambda;
  _lambda_set = true;
}

setMu()

void SymmIsotropicElasticityTensor::setMu

(

const Real

mu

)

Set the second Lame Coefficient.

Definition at line 35 of file SymmIsotropicElasticityTensor.C.

{
  _mu = mu;
  _mu_set = true;
}

Referenced by setShearModulus().

setPoissonsRatio()

void SymmIsotropicElasticityTensor::setPoissonsRatio

(

const Real

nu

)

Set Poissons Ratio.

Definition at line 49 of file SymmIsotropicElasticityTensor.C.

{
  _nu = nu;
  _nu_set = true;
}

Referenced by SolidModel::createElasticityTensor(), LinearIsotropicMaterial::LinearIsotropicMaterial(), and SolidModel::updateElasticityTensor().

setShearModulus()

void SymmIsotropicElasticityTensor::setShearModulus

(

const Real

k

)

Set the shear modulus...

same thing as Mu

Definition at line 63 of file SymmIsotropicElasticityTensor.C.

{
  setMu(k);
}

setYoungsModulus()

void SymmIsotropicElasticityTensor::setYoungsModulus

(

const Real

E

)

Set the Young's Modulus.

Definition at line 42 of file SymmIsotropicElasticityTensor.C.

{
  _E = E;
  _E_set = true;
}

Referenced by SolidModel::createElasticityTensor(), LinearIsotropicMaterial::LinearIsotropicMaterial(), and SolidModel::updateElasticityTensor().

shearModulus()

Real SymmIsotropicElasticityTensor::shearModulus

(

)

const

Return the shear modulus...

same thing as Mu

Definition at line 69 of file SymmIsotropicElasticityTensor.C.

{
  return mu();
}

Referenced by ReturnMappingModel::computeStress(), PowerLawCreepModel::computeStressInitialize(), IsotropicTempDepHardening::computeStressInitialize(), IsotropicPlasticity::computeStressInitialize(), CLSHPlasticModel::computeStressInitialize(), and IsotropicPowerLawHardening::computeStressInitialize().

stiffness() [1/2]

Real SymmIsotropicElasticityTensor::stiffness ( const unsigned  i, const unsigned  j, const RealGradient &  test, const RealGradient &  phi  ) const

virtual

Definition at line 156 of file SymmIsotropicElasticityTensor.C.

{
  RealGradient b;
  if (0 == i && 0 == j)
  {
    b(0) = _val[0] * phi(0);
    b(1) = _val[15] * phi(1);
    b(2) = _val[20] * phi(2);
  }
  else if (1 == i && 1 == j)
  {
    b(0) = _val[15] * phi(0);
    b(1) = _val[6] * phi(1);
    b(2) = _val[18] * phi(2);
  }
  else if (2 == i && 2 == j)
  {
    b(0) = _val[20] * phi(0);
    b(1) = _val[18] * phi(1);
    b(2) = _val[11] * phi(2);
  }
  else if (0 == i && 1 == j)
  {
    b(0) = _val[1] * phi(1);
    b(1) = _val[15] * phi(0);
    b(2) = 0;
  }
  else if (1 == i && 0 == j)
  {
    b(0) = _val[15] * phi(1);
    b(1) = _val[1] * phi(0);
    b(2) = 0;
  }
  else if (1 == i && 2 == j)
  {
    b(0) = 0;
    b(1) = _val[7] * phi(2);
    b(2) = _val[18] * phi(1);
  }
  else if (2 == i && 1 == j)
  {
    b(0) = 0;
    b(1) = _val[18] * phi(2);
    b(2) = _val[7] * phi(1);
  }
  else if (0 == i && 2 == j)
  {
    b(0) = _val[2] * phi(2);
    b(1) = 0;
    b(2) = _val[20] * phi(0);
  }
  else if (2 == i && 0 == j)
  {
    b(0) = _val[20] * phi(2);
    b(1) = 0;
    b(2) = _val[2] * phi(0);
  }
  else
  {
    mooseError("Wrong index in stiffness calculation");
  }
  return test * b;
}

stiffness() [2/2]

Real SymmElasticityTensor::stiffness ( const unsigned int  i, const unsigned int  j, const RealGradient &  test, const RealGradient &  phi  ) const

virtualinherited

Definition at line 87 of file SymmElasticityTensor.C.

{
  RealGradient b;
  if (0 == i && 0 == j)
  {
    b(0) = _val[0] * phi(0) + _val[3] * phi(1) + _val[5] * phi(2);
    b(1) = _val[3] * phi(0) + _val[15] * phi(1) + _val[17] * phi(2);
    b(2) = _val[5] * phi(0) + _val[17] * phi(1) + _val[20] * phi(2);
  }
  else if (1 == i && 1 == j)
  {
    b(0) = _val[15] * phi(0) + _val[8] * phi(1) + _val[16] * phi(2);
    b(1) = _val[8] * phi(0) + _val[6] * phi(1) + _val[9] * phi(2);
    b(2) = _val[16] * phi(0) + _val[9] * phi(1) + _val[18] * phi(2);
  }
  else if (2 == i && 2 == j)
  {
    b(0) = _val[20] * phi(0) + _val[19] * phi(1) + _val[14] * phi(2);
    b(1) = _val[19] * phi(0) + _val[18] * phi(1) + _val[13] * phi(2);
    b(2) = _val[14] * phi(0) + _val[13] * phi(1) + _val[11] * phi(2);
  }
  else if (0 == i && 1 == j)
  {
    b(0) = _val[3] * phi(0) + _val[1] * phi(1) + _val[4] * phi(2);
    b(1) = _val[15] * phi(0) + _val[8] * phi(1) + _val[16] * phi(2);
    b(2) = _val[17] * phi(0) + _val[10] * phi(1) + _val[19] * phi(2);
  }
  else if (1 == i && 0 == j)
  {
    b(0) = _val[3] * phi(0) + _val[15] * phi(1) + _val[17] * phi(2);
    b(1) = _val[1] * phi(0) + _val[8] * phi(1) + _val[10] * phi(2);
    b(2) = _val[4] * phi(0) + _val[16] * phi(1) + _val[19] * phi(2);
  }
  else if (1 == i && 2 == j)
  {
    b(0) = _val[17] * phi(0) + _val[16] * phi(1) + _val[12] * phi(2);
    b(1) = _val[10] * phi(0) + _val[9] * phi(1) + _val[7] * phi(2);
    b(2) = _val[19] * phi(0) + _val[18] * phi(1) + _val[13] * phi(2);
  }
  else if (2 == i && 1 == j)
  {
    b(0) = _val[17] * phi(0) + _val[10] * phi(1) + _val[19] * phi(2);
    b(1) = _val[16] * phi(0) + _val[9] * phi(1) + _val[18] * phi(2);
    b(2) = _val[12] * phi(0) + _val[7] * phi(1) + _val[13] * phi(2);
  }
  else if (0 == i && 2 == j)
  {
    b(0) = _val[5] * phi(0) + _val[4] * phi(1) + _val[2] * phi(2);
    b(1) = _val[17] * phi(0) + _val[16] * phi(1) + _val[12] * phi(2);
    b(2) = _val[20] * phi(0) + _val[19] * phi(1) + _val[14] * phi(2);
  }
  else if (2 == i && 0 == j)
  {
    b(0) = _val[5] * phi(0) + _val[17] * phi(1) + _val[20] * phi(2);
    b(1) = _val[4] * phi(0) + _val[16] * phi(1) + _val[19] * phi(2);
    b(2) = _val[2] * phi(0) + _val[12] * phi(1) + _val[14] * phi(2);
  }
  else
  {
    std::stringstream s;
    s << "Wrong index in stiffness calculation: ";
    s << i << " " << j;
    mooseError(s.str());
  }
  return test * b;
}

sum_3x1()

RealGradient SymmElasticityTensor::sum_3x1 ( ) const

inherited

Definition at line 422 of file SymmElasticityTensor.C.

{
  // used for volumetric locking correction
  RealGradient a(3);
  a(0) = _val[0] + _val[1] + _val[2];  // C00 + C01 + C02
  a(1) = _val[1] + _val[6] + _val[7];  // C10 + C11 + C12
  a(2) = _val[2] + _val[7] + _val[11]; // C20 + C21 + C22
  return a;
}

sum_3x3()

Real SymmElasticityTensor::sum_3x3 ( ) const

inherited

Definition at line 415 of file SymmElasticityTensor.C.

{
  // summation of Cij for i and j ranging from 0 to 2 - used in the volumetric locking correction
  return _val[0] + 2 * (_val[1] + _val[2] + _val[7]) + _val[6] + _val[11];
}

unsetConstants()

void SymmIsotropicElasticityTensor::unsetConstants ( )

inline

Definition at line 40 of file SymmIsotropicElasticityTensor.h.

{ _lambda_set = _mu_set = _E_set = _nu_set = _k_set = false; }

Referenced by SolidModel::updateElasticityTensor().

valueAtIndex()

Real SymmElasticityTensor::valueAtIndex ( int  i ) const

inherited

Definition at line 409 of file SymmElasticityTensor.C.

{
  return _val[i];
}

Referenced by RateDepSmearCrackModel::computeStress().

youngsModulus()

Real SymmIsotropicElasticityTensor::youngsModulus

(

)

const

Return the youngs modulus.

Definition at line 75 of file SymmIsotropicElasticityTensor.C.

{
  if (!_E_set)
    mooseError("Youngs modulus not set");
 
  return _E;
}

Referenced by IsotropicPowerLawHardening::computeStressInitialize().

zero()

void SymmElasticityTensor::zero ( )

inlineinherited

Definition at line 134 of file SymmElasticityTensor.h.

  {
    for (unsigned i(0); i < 21; ++i)
    {
      _val[i] = 0;
    }
  }

Referenced by SymmAnisotropicElasticityTensor::calculateEntries().

Member Data Documentation

_constant

bool SymmElasticityTensor::_constant

protectedinherited

Whether or not the matrix is constant for all of time and space.

Definition at line 173 of file SymmElasticityTensor.h.

Referenced by SymmElasticityTensor::calculate(), SymmElasticityTensor::constant(), dataLoad(), and dataStore().

_E

Real SymmIsotropicElasticityTensor::_E

protected

Definition at line 100 of file SymmIsotropicElasticityTensor.h.

Referenced by calculateLameCoefficients(), setYoungsModulus(), and youngsModulus().

_E_set

bool SymmIsotropicElasticityTensor::_E_set

protected

Definition at line 98 of file SymmIsotropicElasticityTensor.h.

Referenced by calculateLameCoefficients(), setYoungsModulus(), unsetConstants(), and youngsModulus().

_k

Real SymmIsotropicElasticityTensor::_k

protected

Definition at line 100 of file SymmIsotropicElasticityTensor.h.

Referenced by calculateLameCoefficients(), and setBulkModulus().

_k_set

bool SymmIsotropicElasticityTensor::_k_set

protected

Definition at line 98 of file SymmIsotropicElasticityTensor.h.

Referenced by calculateLameCoefficients(), setBulkModulus(), and unsetConstants().

_lambda

Real SymmIsotropicElasticityTensor::_lambda

protected

Definition at line 100 of file SymmIsotropicElasticityTensor.h.

Referenced by adjustForCracking(), calculateEntries(), calculateLameCoefficients(), and setLambda().

_lambda_set

bool SymmIsotropicElasticityTensor::_lambda_set

protected

Definition at line 98 of file SymmIsotropicElasticityTensor.h.

Referenced by calculateLameCoefficients(), setLambda(), and unsetConstants().

_mu

Real SymmIsotropicElasticityTensor::_mu

protected

Definition at line 100 of file SymmIsotropicElasticityTensor.h.

Referenced by adjustForCracking(), calculateEntries(), calculateLameCoefficients(), mu(), and setMu().

_mu_set

bool SymmIsotropicElasticityTensor::_mu_set

protected

Definition at line 98 of file SymmIsotropicElasticityTensor.h.

Referenced by calculateLameCoefficients(), mu(), setMu(), and unsetConstants().

_nu

Real SymmIsotropicElasticityTensor::_nu

protected

Definition at line 100 of file SymmIsotropicElasticityTensor.h.

Referenced by calculateLameCoefficients(), and setPoissonsRatio().

_nu_set

bool SymmIsotropicElasticityTensor::_nu_set

protected

Definition at line 98 of file SymmIsotropicElasticityTensor.h.

Referenced by calculateLameCoefficients(), setPoissonsRatio(), and unsetConstants().

_val

Real SymmElasticityTensor::_val[21]

protectedinherited

Definition at line 188 of file SymmElasticityTensor.h.

Referenced by adjustForCracking(), adjustForCrackingWithShearRetention(), SymmAnisotropicElasticityTensor::calculateEntries(), SymmElasticityTensor::columnMajorMatrix6x6(), SymmElasticityTensor::columnMajorMatrix9x9(), SymmElasticityTensor::convertFrom6x6(), SymmElasticityTensor::convertFrom9x9(), SymmElasticityTensor::copyValues(), dataLoad(), dataStore(), SymmElasticityTensor::fillFromInputVector(), SymmAnisotropicElasticityTensor::initialize_anisotropic_material_dt_matrix(), SymmElasticityTensor::multiply(), multiply(), SymmElasticityTensor::operator*=(), SymmElasticityTensor::operator+=(), SymmElasticityTensor::operator-=(), SymmElasticityTensor::operator/=(), operator<<(), SymmAnisotropicElasticityTensor::rotate(), setEntries(), SymmAnisotropicElasticityTensor::setMaterialConstantc11(), SymmAnisotropicElasticityTensor::setMaterialConstantc12(), SymmAnisotropicElasticityTensor::setMaterialConstantc44(), stiffness(), SymmElasticityTensor::stiffness(), SymmElasticityTensor::sum_3x1(), SymmElasticityTensor::sum_3x3(), SymmAnisotropicElasticityTensor::SymmAnisotropicElasticityTensor(), SymmElasticityTensor::SymmElasticityTensor(), SymmElasticityTensor::valueAtIndex(), and SymmElasticityTensor::zero().

_values_computed

bool SymmElasticityTensor::_values_computed

protectedinherited

Whether or not the values have been computed once.

Definition at line 178 of file SymmElasticityTensor.h.

Referenced by SymmElasticityTensor::calculate(), dataLoad(), and dataStore().

---

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

• solid_mechanics/include/materials/SymmIsotropicElasticityTensor.h
• solid_mechanics/src/materials/SymmIsotropicElasticityTensor.C