CrystalPlasticitySlipRateGSS Class Reference

Phenomenological constitutive model slip rate userobject class. More...

#include <CrystalPlasticitySlipRateGSS.h>

Inheritance diagram for CrystalPlasticitySlipRateGSS:

Public Member Functions
	CrystalPlasticitySlipRateGSS (const InputParameters &parameters)
virtual bool	calcSlipRate (unsigned int qp, Real dt, std::vector< Real > &val) const
virtual bool	calcSlipRateDerivative (unsigned int qp, Real, std::vector< Real > &val) const
virtual void	calcFlowDirection (unsigned int qp, std::vector< RankTwoTensor > &flow_direction) const
virtual bool	calcSlipRate (unsigned qp, Real dt, std::vector< Real > &val) const =0
virtual bool	calcSlipRateDerivative (unsigned qp, Real dt, std::vector< Real > &val) const =0
virtual unsigned int	variableSize () const
	Returns the size of variable. More...

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
virtual void	readFileFlowRateParams ()
virtual void	getFlowRateParams ()
virtual void	getSlipSystems ()

Protected Attributes
const MaterialProperty< std::vector< Real > > &	_mat_prop_state_var
const MaterialProperty< RankTwoTensor > &	_pk2
DenseVector< Real >	_a0
DenseVector< Real >	_xm
const MaterialProperty< std::vector< RankTwoTensor > > &	_flow_direction
unsigned int	_num_slip_sys_props
	Number of slip system specific properties provided in the file containing slip system normals and directions. More...
std::vector< Real >	_flowprops
std::string	_slip_sys_file_name
	File should contain slip plane normal and direction. More...
std::string	_slip_sys_flow_prop_file_name
	File should contain values of the flow rate equation parameters. More...
unsigned int	_num_slip_sys_flowrate_props
	Number of slip system flow rate parameters. More...
Real	_slip_incr_tol
	Slip increment tolerance. More...
DenseVector< Real >	_mo
DenseVector< Real >	_no
const MaterialProperty< RankTwoTensor > &	_crysrot
	Crystal rotation. More...
unsigned int	_variable_size

Detailed Description

Phenomenological constitutive model slip rate userobject class.

Definition at line 23 of file CrystalPlasticitySlipRateGSS.h.

Constructor & Destructor Documentation

CrystalPlasticitySlipRateGSS()

CrystalPlasticitySlipRateGSS::CrystalPlasticitySlipRateGSS

(

const InputParameters &

parameters

)

Definition at line 31 of file CrystalPlasticitySlipRateGSS.C.

  : CrystalPlasticitySlipRate(parameters),
    _mat_prop_state_var(
        getMaterialProperty<std::vector<Real>>(parameters.get<std::string>("uo_state_var_name"))),
    _pk2(getMaterialPropertyByName<RankTwoTensor>("pk2")),
    _a0(_variable_size),
    _xm(_variable_size),
    _flow_direction(getMaterialProperty<std::vector<RankTwoTensor>>(_name + "_flow_direction"))
{
  if (_slip_sys_flow_prop_file_name.length() != 0)
    readFileFlowRateParams();
  else
    getFlowRateParams();
}

Member Function Documentation

calcFlowDirection()

void CrystalPlasticitySlipRateGSS::calcFlowDirection ( unsigned int  qp, std::vector< RankTwoTensor > &  flow_direction  ) const

virtual

Implements CrystalPlasticitySlipRate.

Definition at line 133 of file CrystalPlasticitySlipRateGSS.C.

{
  DenseVector<Real> mo(LIBMESH_DIM * _variable_size), no(LIBMESH_DIM * _variable_size);
 
  // Update slip direction and normal with crystal orientation
  for (unsigned int i = 0; i < _variable_size; ++i)
  {
    for (unsigned int j = 0; j < LIBMESH_DIM; ++j)
    {
      mo(i * LIBMESH_DIM + j) = 0.0;
      for (unsigned int k = 0; k < LIBMESH_DIM; ++k)
        mo(i * LIBMESH_DIM + j) =
            mo(i * LIBMESH_DIM + j) + _crysrot[qp](j, k) * _mo(i * LIBMESH_DIM + k);
    }
 
    for (unsigned int j = 0; j < LIBMESH_DIM; ++j)
    {
      no(i * LIBMESH_DIM + j) = 0.0;
      for (unsigned int k = 0; k < LIBMESH_DIM; ++k)
        no(i * LIBMESH_DIM + j) =
            no(i * LIBMESH_DIM + j) + _crysrot[qp](j, k) * _no(i * LIBMESH_DIM + k);
    }
  }
 
  // Calculate Schmid tensor and resolved shear stresses
  for (unsigned int i = 0; i < _variable_size; ++i)
    for (unsigned int j = 0; j < LIBMESH_DIM; ++j)
      for (unsigned int k = 0; k < LIBMESH_DIM; ++k)
        flow_direction[i](j, k) = mo(i * LIBMESH_DIM + j) * no(i * LIBMESH_DIM + k);
}

calcSlipRate() [1/2]

bool CrystalPlasticitySlipRateGSS::calcSlipRate ( unsigned int  qp, Real  dt, std::vector< Real > &  val  ) const

virtual

Definition at line 166 of file CrystalPlasticitySlipRateGSS.C.

{
  DenseVector<Real> tau(_variable_size);
 
  for (unsigned int i = 0; i < _variable_size; ++i)
    tau(i) = _pk2[qp].doubleContraction(_flow_direction[qp][i]);
 
  for (unsigned int i = 0; i < _variable_size; ++i)
  {
    val[i] = _a0(i) * std::pow(std::abs(tau(i) / _mat_prop_state_var[qp][i]), 1.0 / _xm(i)) *
             std::copysign(1.0, tau(i));
    if (std::abs(val[i] * dt) > _slip_incr_tol)
    {
#ifdef DEBUG
      mooseWarning("Maximum allowable slip increment exceeded ", std::abs(val[i]) * dt);
#endif
      return false;
    }
  }
 
  return true;
}

calcSlipRate() [2/2]

virtual bool CrystalPlasticitySlipRate::calcSlipRate ( unsigned  qp, Real  dt, std::vector< Real > &  val  ) const

pure virtualinherited

calcSlipRateDerivative() [1/2]

bool CrystalPlasticitySlipRateGSS::calcSlipRateDerivative ( unsigned int  qp, Real  , std::vector< Real > &  val  ) const

virtual

Definition at line 190 of file CrystalPlasticitySlipRateGSS.C.

{
  DenseVector<Real> tau(_variable_size);
 
  for (unsigned int i = 0; i < _variable_size; ++i)
    tau(i) = _pk2[qp].doubleContraction(_flow_direction[qp][i]);
 
  for (unsigned int i = 0; i < _variable_size; ++i)
    val[i] = _a0(i) / _xm(i) *
             std::pow(std::abs(tau(i) / _mat_prop_state_var[qp][i]), 1.0 / _xm(i) - 1.0) /
             _mat_prop_state_var[qp][i];
 
  return true;
}

calcSlipRateDerivative() [2/2]

virtual bool CrystalPlasticitySlipRate::calcSlipRateDerivative ( unsigned  qp, Real  dt, std::vector< Real > &  val  ) const

pure virtualinherited

getFlowRateParams()

void CrystalPlasticitySlipRateGSS::getFlowRateParams ( )

protectedvirtual

Reimplemented from CrystalPlasticitySlipRate.

Definition at line 72 of file CrystalPlasticitySlipRateGSS.C.

{
  if (_flowprops.size() <= 0)
    mooseError("CrystalPlasticitySlipRateGSS: Error in reading flow rate  properties: Specify "
               "input in .i file or a slip_sys_flow_prop_file_name");
 
  _a0.resize(_variable_size);
  _xm.resize(_variable_size);
 
  unsigned int num_data_grp = 2 + _num_slip_sys_flowrate_props; // Number of data per group e.g.
                                                                // start_slip_sys, end_slip_sys,
                                                                // value1, value2, ..
 
  for (unsigned int i = 0; i < _flowprops.size() / num_data_grp; ++i)
  {
    Real vs, ve;
    unsigned int is, ie;
 
    vs = _flowprops[i * num_data_grp];
    ve = _flowprops[i * num_data_grp + 1];
 
    if (vs <= 0 || ve <= 0)
      mooseError("CrystalPlasticitySlipRateGSS: Indices in flow rate parameter read must be "
                 "positive integers: is = ",
                 vs,
                 " ie = ",
                 ve);
 
    if (vs != std::floor(vs) || ve != std::floor(ve))
      mooseError("CrystalPlasticitySlipRateGSS: Error in reading flow props: Values specifying "
                 "start and end number of slip system groups should be integer");
 
    is = static_cast<unsigned int>(vs);
    ie = static_cast<unsigned int>(ve);
 
    if (is > ie)
      mooseError("CrystalPlasticitySlipRateGSS: Start index is = ",
                 is,
                 " should be greater than end index ie = ",
                 ie,
                 " in flow rate parameter read");
 
    for (unsigned int j = is; j <= ie; ++j)
    {
      _a0(j - 1) = _flowprops[i * num_data_grp + 2];
      _xm(j - 1) = _flowprops[i * num_data_grp + 3];
    }
  }
 
  for (unsigned int i = 0; i < _variable_size; ++i)
  {
    if (!(_a0(i) > 0.0 && _xm(i) > 0.0))
    {
      mooseWarning(
          "CrystalPlasticitySlipRateGSS: Non-positive flow rate parameters ", _a0(i), ",", _xm(i));
      break;
    }
  }
}

Referenced by CrystalPlasticitySlipRateGSS().

getSlipSystems()

void CrystalPlasticitySlipRate::getSlipSystems ( )

protectedvirtualinherited

Definition at line 68 of file CrystalPlasticitySlipRate.C.

{
  Real vec[LIBMESH_DIM];
  std::ifstream fileslipsys;
 
  MooseUtils::checkFileReadable(_slip_sys_file_name);
 
  fileslipsys.open(_slip_sys_file_name.c_str());
 
  for (unsigned int i = 0; i < _variable_size; ++i)
  {
    // Read the slip normal
    for (unsigned int j = 0; j < LIBMESH_DIM; ++j)
      if (!(fileslipsys >> vec[j]))
        mooseError(
            "CrystalPlasticitySlipRate Error: Premature end of file reading slip system file \n");
 
    // Normalize the vectors
    Real mag;
    mag = Utility::pow<2>(vec[0]) + Utility::pow<2>(vec[1]) + Utility::pow<2>(vec[2]);
    mag = std::sqrt(mag);
 
    for (unsigned j = 0; j < LIBMESH_DIM; ++j)
      _no(i * LIBMESH_DIM + j) = vec[j] / mag;
 
    // Read the slip direction
    for (unsigned int j = 0; j < LIBMESH_DIM; ++j)
      if (!(fileslipsys >> vec[j]))
        mooseError(
            "CrystalPlasticitySlipRate Error: Premature end of file reading slip system file \n");
 
    // Normalize the vectors
    mag = Utility::pow<2>(vec[0]) + Utility::pow<2>(vec[1]) + Utility::pow<2>(vec[2]);
    mag = std::sqrt(mag);
 
    for (unsigned int j = 0; j < LIBMESH_DIM; ++j)
      _mo(i * LIBMESH_DIM + j) = vec[j] / mag;
 
    mag = 0.0;
    for (unsigned int j = 0; j < LIBMESH_DIM; ++j)
      mag += _mo(i * LIBMESH_DIM + j) * _no(i * LIBMESH_DIM + j);
 
    if (std::abs(mag) > 1e-8)
      mooseError("CrystalPlasticitySlipRate Error: Slip direction and normal not orthonormal, "
                 "System number = ",
                 i,
                 "\n");
  }
 
  fileslipsys.close();
}

Referenced by CrystalPlasticitySlipRate::CrystalPlasticitySlipRate().

readFileFlowRateParams()

void CrystalPlasticitySlipRateGSS::readFileFlowRateParams ( )

protectedvirtual

Reimplemented from CrystalPlasticitySlipRate.

Definition at line 47 of file CrystalPlasticitySlipRateGSS.C.

{
  MooseUtils::checkFileReadable(_slip_sys_flow_prop_file_name);
 
  std::ifstream file;
  file.open(_slip_sys_flow_prop_file_name.c_str());
 
  std::vector<Real> vec;
  vec.resize(_num_slip_sys_flowrate_props);
 
  for (unsigned int i = 0; i < _variable_size; ++i)
  {
    for (unsigned int j = 0; j < _num_slip_sys_flowrate_props; ++j)
      if (!(file >> vec[j]))
        mooseError(
            "Error CrystalPlasticitySlipRateGSS: Premature end of slip_sys_flow_rate_param file");
 
    _a0(i) = vec[0];
    _xm(i) = vec[1];
  }
 
  file.close();
}

Referenced by CrystalPlasticitySlipRateGSS().

validParams()

InputParameters CrystalPlasticitySlipRateGSS::validParams ( )

static

Definition at line 19 of file CrystalPlasticitySlipRateGSS.C.

{
  InputParameters params = CrystalPlasticitySlipRate::validParams();
  params.addParam<std::string>("uo_state_var_name",
                               "Name of state variable property: Same as "
                               "state variable user object specified in input "
                               "file.");
  params.addClassDescription("Phenomenological constitutive model slip rate class.  Override the "
                             "virtual functions in your class");
  return params;
}

variableSize()

unsigned int CrystalPlasticityUOBase::variableSize ( ) const

virtualinherited

Returns the size of variable.

Definition at line 30 of file CrystalPlasticityUOBase.C.

{
  return _variable_size;
}

Member Data Documentation

_a0

DenseVector<Real> CrystalPlasticitySlipRateGSS::_a0

protected

Definition at line 43 of file CrystalPlasticitySlipRateGSS.h.

Referenced by calcSlipRate(), calcSlipRateDerivative(), getFlowRateParams(), and readFileFlowRateParams().

_crysrot

const MaterialProperty<RankTwoTensor>& CrystalPlasticitySlipRate::_crysrot

protectedinherited

Crystal rotation.

Definition at line 70 of file CrystalPlasticitySlipRate.h.

Referenced by calcFlowDirection().

_flow_direction

const MaterialProperty<std::vector<RankTwoTensor> >& CrystalPlasticitySlipRateGSS::_flow_direction

protected

Definition at line 46 of file CrystalPlasticitySlipRateGSS.h.

Referenced by calcSlipRate(), and calcSlipRateDerivative().

_flowprops

std::vector<Real> CrystalPlasticitySlipRate::_flowprops

protectedinherited

Definition at line 47 of file CrystalPlasticitySlipRate.h.

Referenced by getFlowRateParams().

_mat_prop_state_var

const MaterialProperty<std::vector<Real> >& CrystalPlasticitySlipRateGSS::_mat_prop_state_var

protected

Definition at line 39 of file CrystalPlasticitySlipRateGSS.h.

Referenced by calcSlipRate(), and calcSlipRateDerivative().

_mo

DenseVector<Real> CrystalPlasticitySlipRate::_mo

protectedinherited

Definition at line 66 of file CrystalPlasticitySlipRate.h.

Referenced by calcFlowDirection(), and CrystalPlasticitySlipRate::getSlipSystems().

_no

DenseVector<Real> CrystalPlasticitySlipRate::_no

protectedinherited

Definition at line 67 of file CrystalPlasticitySlipRate.h.

Referenced by calcFlowDirection(), and CrystalPlasticitySlipRate::getSlipSystems().

_num_slip_sys_flowrate_props

unsigned int CrystalPlasticitySlipRate::_num_slip_sys_flowrate_props

protectedinherited

Number of slip system flow rate parameters.

Definition at line 61 of file CrystalPlasticitySlipRate.h.

Referenced by getFlowRateParams(), and readFileFlowRateParams().

_num_slip_sys_props

unsigned int CrystalPlasticitySlipRate::_num_slip_sys_props

protectedinherited

Number of slip system specific properties provided in the file containing slip system normals and directions.

Definition at line 45 of file CrystalPlasticitySlipRate.h.

_pk2

const MaterialProperty<RankTwoTensor>& CrystalPlasticitySlipRateGSS::_pk2

protected

Definition at line 41 of file CrystalPlasticitySlipRateGSS.h.

Referenced by calcSlipRate(), and calcSlipRateDerivative().

_slip_incr_tol

Real CrystalPlasticitySlipRate::_slip_incr_tol

protectedinherited

Slip increment tolerance.

Definition at line 64 of file CrystalPlasticitySlipRate.h.

Referenced by calcSlipRate().

_slip_sys_file_name

std::string CrystalPlasticitySlipRate::_slip_sys_file_name

protectedinherited

File should contain slip plane normal and direction.

Definition at line 50 of file CrystalPlasticitySlipRate.h.

Referenced by CrystalPlasticitySlipRate::getSlipSystems().

_slip_sys_flow_prop_file_name

std::string CrystalPlasticitySlipRate::_slip_sys_flow_prop_file_name

protectedinherited

File should contain values of the flow rate equation parameters.

Values for every slip system must be provided. Should have the same order of slip systens as in slip_sys_file. The option of reading all the properties from .i is still present.

Definition at line 58 of file CrystalPlasticitySlipRate.h.

Referenced by CrystalPlasticitySlipRateGSS(), and readFileFlowRateParams().

_variable_size

unsigned int CrystalPlasticityUOBase::_variable_size

protectedinherited

Definition at line 33 of file CrystalPlasticityUOBase.h.

Referenced by calcFlowDirection(), calcSlipRate(), calcSlipRateDerivative(), CrystalPlasticitySlipResistanceGSS::calcSlipResistance(), CrystalPlasticityStateVarRateComponentGSS::calcStateVariableEvolutionRateComponent(), CrystalPlasticityStateVarRateComponentVoce::calcStateVariableEvolutionRateComponent(), getFlowRateParams(), CrystalPlasticitySlipRate::getSlipSystems(), CrystalPlasticityStateVariable::initSlipSysProps(), CrystalPlasticityStateVarRateComponentVoce::initSlipSystemGroupID(), CrystalPlasticityStateVarRateComponentVoce::initSlipSystemPlaneID(), readFileFlowRateParams(), CrystalPlasticityStateVariable::readInitialValueFromFile(), CrystalPlasticityStateVariable::updateStateVariable(), and CrystalPlasticityUOBase::variableSize().

_xm

DenseVector<Real> CrystalPlasticitySlipRateGSS::_xm

protected

Definition at line 44 of file CrystalPlasticitySlipRateGSS.h.

Referenced by calcSlipRate(), calcSlipRateDerivative(), getFlowRateParams(), and readFileFlowRateParams().

---

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

• tensor_mechanics/include/userobjects/CrystalPlasticitySlipRateGSS.h
• tensor_mechanics/src/userobjects/CrystalPlasticitySlipRateGSS.C