SymmAnisotropicElasticityTensor Class Reference

#include <SymmAnisotropicElasticityTensor.h>

Inheritance diagram for SymmAnisotropicElasticityTensor:

Public Member Functions
	SymmAnisotropicElasticityTensor ()
	SymmAnisotropicElasticityTensor (std::vector< Real > &init_list, bool all_21)
	SymmAnisotropicElasticityTensor (const SymmAnisotropicElasticityTensor &a)
SymmAnisotropicElasticityTensor &	operator= (const SymmAnisotropicElasticityTensor &a)=default
virtual	~SymmAnisotropicElasticityTensor ()=default
void	setFirstEulerAngle (const Real a1)
	Set the first euler angle. More...
void	setSecondEulerAngle (const Real a2)
	Set the second euler angle. More...
void	setThirdEulerAngle (const Real a3)
	Set the third euler angle. More...
void	setMaterialConstantc11 (const Real c11)
	Set the material constant c11; assumes cubic material. More...
void	setMaterialConstantc12 (const Real c12)
	Set the material constant c22; assumes cubic material. More...
void	setMaterialConstantc44 (const Real c44)
	Set the material constant c44; assumes cubic material. More...
virtual void	rotate (const Real a1, const Real a2, const Real a3)
	Perform rotation around three axes. More...
Real	firstEulerAngle ()
Real	secondEulerAngle ()
Real	thirdEulerAngle ()
void	show_dt_matrix ()
void	show_r_matrix ()
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...
virtual void	multiply (const SymmTensor &x, SymmTensor &b) const
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 void	adjustForCracking (const RealVectorValue &crack_flags)
virtual void	adjustForCrackingWithShearRetention (const RealVectorValue &crack_flags)
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
void	form_r_matrix ()
void	initialize_material_dt_matrix ()
void	initialize_anisotropic_material_dt_matrix ()
void	form_rotational_q_matrix ()
void	form_transformation_t_matrix ()
void	form_transformed_material_dmat_matrix ()
void	form_transformed_material_dt_matrix ()
void	form_rotated_material_qdmat_matrix ()
virtual void	calculateEntries (unsigned int qp)
	Fill in the matrix. More...

Protected Attributes
DenseMatrix< Real >	_dmat
DenseMatrix< Real >	_qdmat
DenseMatrix< Real >	_dt
DenseMatrix< Real >	_qdt
DenseMatrix< Real >	_r
DenseMatrix< Real >	_q
DenseMatrix< Real >	_qt
std::vector< Real >	_euler_angle
DenseMatrix< Real >	_trans_d6_to_d9
DenseMatrix< Real >	_trans_d9_to_d6
Real	_c11
Real	_c12
Real	_c44
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

Definition at line 14 of file SymmAnisotropicElasticityTensor.h.

Constructor & Destructor Documentation

SymmAnisotropicElasticityTensor() [1/3]

SymmAnisotropicElasticityTensor::SymmAnisotropicElasticityTensor

(

)

Definition at line 14 of file SymmAnisotropicElasticityTensor.C.

  : SymmElasticityTensor(false),
    _dmat(9, 9),
    _qdmat(9, 9),
    _dt(6, 6),
    _qdt(6, 6),
    _r(3, 3),
    _q(9, 9),
    _qt(9, 9),
    _euler_angle(3),
    _trans_d6_to_d9(9, 6),
    //_transpose_trans_d6_to_d9(6,9),
    _trans_d9_to_d6(6, 9),
    //_transpose_trans_d9_to_d6(9,6),
    _c11(0),
    _c12(0),
    _c44(0)
{
  // form_transformation_t_matrix();
}

SymmAnisotropicElasticityTensor() [2/3]

SymmAnisotropicElasticityTensor::SymmAnisotropicElasticityTensor	(	std::vector< Real > &	init_list,
		bool	all_21
	)

Definition at line 35 of file SymmAnisotropicElasticityTensor.C.

  : SymmElasticityTensor(false),
    _dmat(9, 9),
    _qdmat(9, 9),
    _dt(6, 6),
    _qdt(6, 6),
    _r(3, 3),
    _q(9, 9),
    _qt(9, 9),
    _euler_angle(3),
    _trans_d6_to_d9(9, 6),
    //_transpose_trans_d6_to_d9(6,9),
    _trans_d9_to_d6(6, 9),
    //_transpose_trans_d9_to_d6(9,6),
    _c11(0),
    _c12(0),
    _c44(0)
{
  // test the input vector length to make sure it's correct
  if ((all_21 == true && init_list.size() != 21) || (all_21 == false && init_list.size() != 9))
    mooseError("Please correct the number of entries in the stiffness input.");
 
  if (all_21 == true)
  {
    for (int i = 0; i < 21; i++)
      _val[i] = init_list[i];
  }
  else
  {
    _val[0] = init_list[0];  // C1111
    _val[1] = init_list[1];  // C1122
    _val[2] = init_list[2];  // C1133
    _val[6] = init_list[3];  // C2222
    _val[7] = init_list[4];  // C2233
    _val[11] = init_list[5]; // C3333
    _val[15] = init_list[6]; // C2323
    _val[18] = init_list[7]; // C1313
    _val[20] = init_list[8]; // C1212
  }
 
  // form_transformation_t_matrix();
}

SymmAnisotropicElasticityTensor() [3/3]

SymmAnisotropicElasticityTensor::SymmAnisotropicElasticityTensor

(

const SymmAnisotropicElasticityTensor &

a

)

Definition at line 79 of file SymmAnisotropicElasticityTensor.C.

  : SymmElasticityTensor(false)
{
  *this = a;
}

~SymmAnisotropicElasticityTensor()

virtual SymmAnisotropicElasticityTensor::~SymmAnisotropicElasticityTensor ( )

virtualdefault

Member Function Documentation

adjustForCracking()

void SymmElasticityTensor::adjustForCracking ( const RealVectorValue &  crack_flags )

virtualinherited

Reimplemented in SymmIsotropicElasticityTensor.

Definition at line 372 of file SymmElasticityTensor.C.

{
  mooseError("adjustForCracking method not defined");
}

Referenced by SolidModel::crackingStrainDirections().

adjustForCrackingWithShearRetention()

void SymmElasticityTensor::adjustForCrackingWithShearRetention ( const RealVectorValue &  crack_flags )

virtualinherited

Reimplemented in SymmIsotropicElasticityTensor.

Definition at line 378 of file SymmElasticityTensor.C.

{
  mooseError("adjustForCrackingWithShearRetention method not defined");
}

Referenced by SolidModel::crackingStrainDirections().

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 SymmAnisotropicElasticityTensor::calculateEntries ( unsigned int  qp )

protectedvirtual

Fill in the matrix.

Reimplemented from SymmElasticityTensor.

Definition at line 365 of file SymmAnisotropicElasticityTensor.C.

{
 
  // The following four lines of code force the calculateEntries function to be useful
  // only for CUBIC ANISOTROPIC materials.
  zero();
  setMaterialConstantc11(_c11);
  setMaterialConstantc12(_c12);
  setMaterialConstantc44(_c44);
 
  form_r_matrix();
  // initialize_material_anisotropic_dt_matrix();
  form_rotational_q_matrix();
  // form_transformation_t_matrix();
  form_transformed_material_dmat_matrix();
  form_rotated_material_qdmat_matrix();
  form_transformed_material_dt_matrix();
 
  unsigned count(0);
 
  for (int j(0); j < 6; ++j)
  {
    for (int i(j); i < 6; ++i)
    {
      _val[count++] = _dt(i, j);
    }
  }
}

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 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 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 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 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().

firstEulerAngle()

Real SymmAnisotropicElasticityTensor::firstEulerAngle

(

)

Returns

the first Euler angle

Definition at line 105 of file SymmAnisotropicElasticityTensor.C.

{
  return _euler_angle[0];
}

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().

form_r_matrix()

void SymmAnisotropicElasticityTensor::form_r_matrix ( )

protected

Definition at line 171 of file SymmAnisotropicElasticityTensor.C.

{
  Real phi1 = _euler_angle[0] * (libMesh::pi / 180.0);
  Real phi = _euler_angle[1] * (libMesh::pi / 180.0);
  Real phi2 = _euler_angle[2] * (libMesh::pi / 180.0);
 
  Real cp1 = std::cos(phi1);
  Real cp2 = std::cos(phi2);
  Real cp = std::cos(phi);
 
  Real sp1 = std::sin(phi1);
  Real sp2 = std::sin(phi2);
  Real sp = std::sin(phi);
 
  _r(0, 0) = cp1 * cp2 - sp1 * sp2 * cp;
  _r(0, 1) = sp1 * cp2 + cp1 * sp2 * cp;
  _r(0, 2) = sp2 * sp;
  _r(1, 0) = -cp1 * sp2 - sp1 * cp2 * cp;
  _r(1, 1) = -sp1 * sp2 + cp1 * cp2 * cp;
  _r(1, 2) = cp2 * sp;
  _r(2, 0) = sp1 * sp;
  _r(2, 1) = -cp1 * sp;
  _r(2, 2) = cp;
}

Referenced by calculateEntries(), and rotate().

form_rotated_material_qdmat_matrix()

void SymmAnisotropicElasticityTensor::form_rotated_material_qdmat_matrix ( )

protected

Definition at line 349 of file SymmAnisotropicElasticityTensor.C.

{
  // The function makes use of Q matrix to rotate
  // Dmat[9][9] to QDmat[9][9]
  // QDmat = QT * Dmat * Q
 
  DenseMatrix<Real> outputMatrix(9, 9);
 
  _q.get_transpose(outputMatrix);
  outputMatrix.right_multiply(_dmat);
 
  _qdmat = outputMatrix;
  _qdmat.right_multiply(_q);
}

Referenced by calculateEntries(), and rotate().

form_rotational_q_matrix()

void SymmAnisotropicElasticityTensor::form_rotational_q_matrix ( )

protected

Definition at line 225 of file SymmAnisotropicElasticityTensor.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)
          _q(((i * 3) + k), ((j * 3) + l)) = _r(i, j) * _r(k, l);
 
  /*for (int p = 0; p < 9; ++p)
    for (int q = 0; q < 9; ++q)
    _qt(q,p) = _q(p,q);*/
}

Referenced by calculateEntries(), and rotate().

form_transformation_t_matrix()

void SymmAnisotropicElasticityTensor::form_transformation_t_matrix ( )

protected

Definition at line 240 of file SymmAnisotropicElasticityTensor.C.

{
  // Forms to two kinds of transformation matrix
  // TransD6toD9 transfroms Dt[6][6] to Dmat[9][9]
  // TransD9toD6 transforms Dmat[9][9] to Dt[6][6]
 
  //  Real sqrt2 = std::sqrt(2.0);
  //  Real a = 1/sqrt2;
  Real a = 1.0;
 
  _trans_d6_to_d9(0, 0) = _trans_d6_to_d9(4, 1) = _trans_d6_to_d9(8, 2) = 1.0;
  _trans_d6_to_d9(1, 3) = _trans_d6_to_d9(3, 3) = a;
  _trans_d6_to_d9(5, 4) = _trans_d6_to_d9(7, 4) = a;
  _trans_d6_to_d9(2, 5) = _trans_d6_to_d9(6, 5) = a;
 
  /*for (int i = 0; i < 9; ++i)
    {
    for (int j = 0; j < 6; ++j)
    {
    _transpose_trans_d6_to_d9(j,i) = _trans_d6_to_d9(i,j);
    }
    }*/
 
  _trans_d9_to_d6(0, 0) = _trans_d9_to_d6(1, 4) = _trans_d9_to_d6(2, 8) = 1.0;
  //  _trans_d9_to_d6(3,3) = _trans_d9_to_d6(4,7) =  _trans_d9_to_d6(5,6) = sqrt2;
  _trans_d9_to_d6(3, 3) = _trans_d9_to_d6(4, 7) = _trans_d9_to_d6(5, 6) = 1.0;
 
  /*for (int i = 0; i < 6; ++i)
    {
    for (int j = 0; j < 9; ++j)
    {
    _transpose_trans_d9_to_d6(j,i) =  _trans_d9_to_d6(i,j);
    }
    }*/
}

form_transformed_material_dmat_matrix()

void SymmAnisotropicElasticityTensor::form_transformed_material_dmat_matrix ( )

protected

Definition at line 277 of file SymmAnisotropicElasticityTensor.C.

{
 
  // THIS TRANSFORMATION IS VALID ONLY WHEN THE INCOMING MATRIX HAS NOT BEEN ROTATED
 
  // The function makes use of TransD6toD9 matrix to transfrom
  // Dt[6][6] to Dmat[9][9]
  // Dmat = T * Dt * TT
 
  // DenseMatrix<Real> outputMatrix(9,6);
 
  // outputMatrix = _trans_d6_to_d9;
  // outputMatrix.right_multiply(_dt);
 
  //_dmat = outputMatrix;
  //_dmat.right_multiply_transpose(_trans_d6_to_d9);
 
  // Use the plug-and-chug functions given in SymmElasticityTensor
  // to take the existing _val[] and put them into the 9x9 _dmat matrix.
  SymmElasticityTensor temp_dt;
  copyValues(temp_dt);
 
  ColumnMajorMatrix temp_dmat = temp_dt.columnMajorMatrix9x9();
 
  for (unsigned j(0); j < 9; ++j)
  {
    for (unsigned i(0); i < 9; ++i)
    {
      _dmat(i, j) = temp_dmat(i, j);
    }
  }
}

Referenced by calculateEntries(), and rotate().

form_transformed_material_dt_matrix()

void SymmAnisotropicElasticityTensor::form_transformed_material_dt_matrix ( )

protected

Definition at line 311 of file SymmAnisotropicElasticityTensor.C.

{
  // The function makes use of TransD6toD9 matrix to transfrom
  // QDmat[9][9] to Dt[6][6]
  // Dt = TT * QDmat * T
 
  //   DenseMatrix<Real> outputMatrix(6,9);
 
  //   outputMatrix = _trans_d9_to_d6;
  //   outputMatrix.right_multiply(_qdmat);
  //   _dt = outputMatrix;
  //   _dt.right_multiply(_transpose_trans_d9_to_d6);
 
  // The transformation below is general and should work whether or not the
  //   incoming tensor has been rotated.
 
  ColumnMajorMatrix tmp(9, 9);
  for (unsigned j(0); j < 9; ++j)
  {
    for (unsigned i(0); i < 9; ++i)
    {
      tmp(i, j) = _qdmat(i, j);
    }
  }
 
  SymmElasticityTensor fred;
  fred.convertFrom9x9(tmp);
  ColumnMajorMatrix wilma = fred.columnMajorMatrix6x6();
  for (unsigned j(0); j < 6; ++j)
  {
    for (unsigned i(0); i < 6; ++i)
    {
      _dt(i, j) = wilma(i, j);
    }
  }
}

Referenced by calculateEntries(), and rotate().

initialize_anisotropic_material_dt_matrix()

void SymmAnisotropicElasticityTensor::initialize_anisotropic_material_dt_matrix ( )

protected

Definition at line 210 of file SymmAnisotropicElasticityTensor.C.

{
  // This function initializes the 6 x 6 material Dt matrix for an anisotropic material
 
  int k = 0;
  for (int i = 0; i < 6; i++)
  {
    for (int j = i; j < 6; j++)
    {
      _dt(i, j) = _dt(j, i) = _val[k++];
    }
  }
}

initialize_material_dt_matrix()

void SymmAnisotropicElasticityTensor::initialize_material_dt_matrix ( )

protected

Definition at line 198 of file SymmAnisotropicElasticityTensor.C.

{
  // This function initializes the 6 x 6 material Dt matrix for a cubic material
 
  _dt(0, 0) = _dt(1, 1) = _dt(2, 2) = _c11;
  _dt(0, 1) = _dt(0, 2) = _dt(1, 0) = _dt(2, 0) = _dt(1, 2) = _dt(2, 1) = _c12;
  // beware the factor of two here
  _dt(3, 3) = _dt(4, 4) = _dt(5, 5) = 2 * _c44;
}

multiply()

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

virtualinherited

Reimplemented in SymmIsotropicElasticityTensor.

Definition at line 56 of file SymmElasticityTensor.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 + 2 * (_val[3] * xy + _val[4] * yz + _val[5] * zx);
  b.yy() = _val[1] * xx + _val[6] * yy + _val[7] * zz +
           2 * (_val[8] * xy + _val[9] * yz + _val[10] * zx);
  b.zz() = _val[2] * xx + _val[7] * yy + _val[11] * zz +
           2 * (_val[12] * xy + _val[13] * yz + _val[14] * zx);
  b.xy() = _val[3] * xx + _val[8] * yy + _val[12] * zz +
           2 * (_val[15] * xy + _val[16] * yz + _val[17] * zx);
  b.yz() = _val[4] * xx + _val[9] * yy + _val[13] * zz +
           2 * (_val[16] * xy + _val[18] * yz + _val[19] * zx);
  b.zx() = _val[5] * xx + _val[10] * yy + _val[14] * zz +
           2 * (_val[17] * xy + _val[19] * yz + _val[20] * zx);
}

Referenced by SymmElasticityTensor::operator*().

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;
    }
  }

operator=()

SymmAnisotropicElasticityTensor& SymmAnisotropicElasticityTensor::operator= ( const SymmAnisotropicElasticityTensor &  a )

default

rotate()

void SymmAnisotropicElasticityTensor::rotate ( const Real  a1, const Real  a2, const Real  a3  )

virtual

Perform rotation around three axes.

Definition at line 144 of file SymmAnisotropicElasticityTensor.C.

{
  setFirstEulerAngle(a1);
  setSecondEulerAngle(a2);
  setThirdEulerAngle(a3);
 
  // pulled from calculateEntries to sub in the initialize_anisotropic_material_dt_matrix() call
  //  calculateEntries(0);
 
  form_r_matrix();
  form_rotational_q_matrix();
  form_transformed_material_dmat_matrix();
  form_rotated_material_qdmat_matrix();
  form_transformed_material_dt_matrix();
 
  unsigned count(0);
 
  for (int j(0); j < 6; ++j)
  {
    for (int i(j); i < 6; ++i)
    {
      _val[count++] = _dt(i, j);
    }
  }
}

Referenced by LinearGeneralAnisotropicMaterial::LinearGeneralAnisotropicMaterial().

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().

secondEulerAngle()

Real SymmAnisotropicElasticityTensor::secondEulerAngle

(

)

Returns

the first Euler angle

Definition at line 111 of file SymmAnisotropicElasticityTensor.C.

{
  return _euler_angle[1];
}

setFirstEulerAngle()

void SymmAnisotropicElasticityTensor::setFirstEulerAngle

(

const Real

a1

)

Set the first euler angle.

Definition at line 87 of file SymmAnisotropicElasticityTensor.C.

{
  _euler_angle[0] = a1;
}

Referenced by LinearAnisotropicMaterial::LinearAnisotropicMaterial(), and rotate().

setMaterialConstantc11()

void SymmAnisotropicElasticityTensor::setMaterialConstantc11

(

const Real

c11

)

Set the material constant c11; assumes cubic material.

Definition at line 123 of file SymmAnisotropicElasticityTensor.C.

{
  _c11 = c11;
  _val[0] = _val[6] = _val[11] = _c11;
}

Referenced by calculateEntries(), and LinearAnisotropicMaterial::LinearAnisotropicMaterial().

setMaterialConstantc12()

void SymmAnisotropicElasticityTensor::setMaterialConstantc12

(

const Real

c12

)

Set the material constant c22; assumes cubic material.

Definition at line 130 of file SymmAnisotropicElasticityTensor.C.

{
  _c12 = c12;
  _val[1] = _val[2] = _val[7] = _c12;
}

Referenced by calculateEntries(), and LinearAnisotropicMaterial::LinearAnisotropicMaterial().

setMaterialConstantc44()

void SymmAnisotropicElasticityTensor::setMaterialConstantc44

(

const Real

c44

)

Set the material constant c44; assumes cubic material.

Definition at line 137 of file SymmAnisotropicElasticityTensor.C.

{
  _c44 = c44;
  _val[15] = _val[18] = _val[20] = _c44;
}

Referenced by calculateEntries(), and LinearAnisotropicMaterial::LinearAnisotropicMaterial().

setSecondEulerAngle()

void SymmAnisotropicElasticityTensor::setSecondEulerAngle

(

const Real

a2

)

Set the second euler angle.

Definition at line 93 of file SymmAnisotropicElasticityTensor.C.

{
  _euler_angle[1] = a2;
}

Referenced by LinearAnisotropicMaterial::LinearAnisotropicMaterial(), and rotate().

setThirdEulerAngle()

void SymmAnisotropicElasticityTensor::setThirdEulerAngle

(

const Real

a3

)

Set the third euler angle.

Definition at line 99 of file SymmAnisotropicElasticityTensor.C.

{
  _euler_angle[2] = a3;
}

Referenced by LinearAnisotropicMaterial::LinearAnisotropicMaterial(), and rotate().

show_dt_matrix()

void SymmAnisotropicElasticityTensor::show_dt_matrix

(

)

Definition at line 395 of file SymmAnisotropicElasticityTensor.C.

{
  printf("\nSymmAnisotropicElasticityTensor::show_dt_matrix()\n");
 
  for (int j = 0; j < 6; ++j)
  {
    printf("  ");
    for (int i = 0; i < 6; ++i)
    {
      printf("%12.4f  ", _dt(i, j));
    }
    printf("\n");
  }
}

show_r_matrix()

void SymmAnisotropicElasticityTensor::show_r_matrix

(

)

Definition at line 411 of file SymmAnisotropicElasticityTensor.C.

{
  printf("\nSymmAnisotropicElasticityTensor::show_r_matrix()  Euler angles are (%f, %f, %f)\n",
         _euler_angle[0],
         _euler_angle[1],
         _euler_angle[2]);
 
  for (int j = 0; j < 3; ++j)
  {
    printf("  ");
    for (int i = 0; i < 3; ++i)
    {
      printf("%8.4f  ", _r(i, j));
    }
    printf("\n");
  }
}

stiffness()

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];
}

thirdEulerAngle()

Real SymmAnisotropicElasticityTensor::thirdEulerAngle

(

)

Returns

the first Euler angle

Definition at line 117 of file SymmAnisotropicElasticityTensor.C.

{
  return _euler_angle[2];
}

valueAtIndex()

Real SymmElasticityTensor::valueAtIndex ( int  i ) const

inherited

Definition at line 409 of file SymmElasticityTensor.C.

{
  return _val[i];
}

Referenced by RateDepSmearCrackModel::computeStress().

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 calculateEntries().

Member Data Documentation

_c11

Real SymmAnisotropicElasticityTensor::_c11

protected

Definition at line 104 of file SymmAnisotropicElasticityTensor.h.

Referenced by calculateEntries(), initialize_material_dt_matrix(), and setMaterialConstantc11().

_c12

Real SymmAnisotropicElasticityTensor::_c12

protected

Definition at line 104 of file SymmAnisotropicElasticityTensor.h.

Referenced by calculateEntries(), initialize_material_dt_matrix(), and setMaterialConstantc12().

_c44

Real SymmAnisotropicElasticityTensor::_c44

protected

Definition at line 104 of file SymmAnisotropicElasticityTensor.h.

Referenced by calculateEntries(), initialize_material_dt_matrix(), and setMaterialConstantc44().

_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().

_dmat

DenseMatrix<Real> SymmAnisotropicElasticityTensor::_dmat

protected

Definition at line 89 of file SymmAnisotropicElasticityTensor.h.

Referenced by form_rotated_material_qdmat_matrix(), and form_transformed_material_dmat_matrix().

_dt

DenseMatrix<Real> SymmAnisotropicElasticityTensor::_dt

protected

Definition at line 91 of file SymmAnisotropicElasticityTensor.h.

Referenced by calculateEntries(), form_transformed_material_dt_matrix(), initialize_anisotropic_material_dt_matrix(), initialize_material_dt_matrix(), rotate(), and show_dt_matrix().

_euler_angle

std::vector<Real> SymmAnisotropicElasticityTensor::_euler_angle

protected

Definition at line 96 of file SymmAnisotropicElasticityTensor.h.

Referenced by firstEulerAngle(), form_r_matrix(), secondEulerAngle(), setFirstEulerAngle(), setSecondEulerAngle(), setThirdEulerAngle(), show_r_matrix(), and thirdEulerAngle().

_q

DenseMatrix<Real> SymmAnisotropicElasticityTensor::_q

protected

Definition at line 94 of file SymmAnisotropicElasticityTensor.h.

Referenced by form_rotated_material_qdmat_matrix(), and form_rotational_q_matrix().

_qdmat

DenseMatrix<Real> SymmAnisotropicElasticityTensor::_qdmat

protected

Definition at line 90 of file SymmAnisotropicElasticityTensor.h.

Referenced by form_rotated_material_qdmat_matrix(), and form_transformed_material_dt_matrix().

_qdt

DenseMatrix<Real> SymmAnisotropicElasticityTensor::_qdt

protected

Definition at line 92 of file SymmAnisotropicElasticityTensor.h.

_qt

DenseMatrix<Real> SymmAnisotropicElasticityTensor::_qt

protected

Definition at line 95 of file SymmAnisotropicElasticityTensor.h.

_r

DenseMatrix<Real> SymmAnisotropicElasticityTensor::_r

protected

Definition at line 93 of file SymmAnisotropicElasticityTensor.h.

Referenced by form_r_matrix(), form_rotational_q_matrix(), and show_r_matrix().

_trans_d6_to_d9

DenseMatrix<Real> SymmAnisotropicElasticityTensor::_trans_d6_to_d9

protected

Definition at line 98 of file SymmAnisotropicElasticityTensor.h.

Referenced by form_transformation_t_matrix().

_trans_d9_to_d6

DenseMatrix<Real> SymmAnisotropicElasticityTensor::_trans_d9_to_d6

protected

Definition at line 101 of file SymmAnisotropicElasticityTensor.h.

Referenced by form_transformation_t_matrix().

_val

Real SymmElasticityTensor::_val[21]

protectedinherited

Definition at line 188 of file SymmElasticityTensor.h.

Referenced by SymmIsotropicElasticityTensor::adjustForCracking(), SymmIsotropicElasticityTensor::adjustForCrackingWithShearRetention(), calculateEntries(), SymmElasticityTensor::columnMajorMatrix6x6(), SymmElasticityTensor::columnMajorMatrix9x9(), SymmElasticityTensor::convertFrom6x6(), SymmElasticityTensor::convertFrom9x9(), SymmElasticityTensor::copyValues(), dataLoad(), dataStore(), SymmElasticityTensor::fillFromInputVector(), initialize_anisotropic_material_dt_matrix(), SymmElasticityTensor::multiply(), SymmIsotropicElasticityTensor::multiply(), SymmElasticityTensor::operator*=(), SymmElasticityTensor::operator+=(), SymmElasticityTensor::operator-=(), SymmElasticityTensor::operator/=(), operator<<(), rotate(), SymmIsotropicElasticityTensor::setEntries(), setMaterialConstantc11(), setMaterialConstantc12(), setMaterialConstantc44(), SymmIsotropicElasticityTensor::stiffness(), SymmElasticityTensor::stiffness(), SymmElasticityTensor::sum_3x1(), SymmElasticityTensor::sum_3x3(), 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/SymmAnisotropicElasticityTensor.h
• solid_mechanics/src/materials/SymmAnisotropicElasticityTensor.C