ComputeIncrementalBeamStrain Class Reference

#include <ComputeIncrementalBeamStrain.h>

Inheritance diagram for ComputeIncrementalBeamStrain:

Public Member Functions
	ComputeIncrementalBeamStrain (const InputParameters &parameters)
virtual void	computeProperties () override

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
virtual void	initQpStatefulProperties () override
void	computeQpStrain ()
	Computes the displacement and rotation strain increments. More...
void	computeStiffnessMatrix ()
	Computes the stiffness matrices. More...
virtual void	computeRotation ()
	Computes the rotation matrix at time t. For small rotation scenarios, the rotation matrix at time t is same as the intiial rotation matrix. More...

Protected Attributes
const bool	_has_Ix
	Booleans for validity of params. More...
unsigned int	_nrot
	Number of coupled rotational variables. More...
unsigned int	_ndisp
	Number of coupled displacement variables. More...
std::vector< unsigned int >	_rot_num
	Variable numbers corresponding to the rotational variables. More...
std::vector< unsigned int >	_disp_num
	Variable numbers corresponding to the displacement variables. More...
const VariableValue &	_area
	Coupled variable for the beam cross-sectional area. More...
const VariableValue &	_Ay
	Coupled variable for the first moment of area in y direction, i.e., integral of y*dA over the cross-section. More...
const VariableValue &	_Az
	Coupled variable for the first moment of area in z direction, i.e., integral of z*dA over the cross-section. More...
const VariableValue &	_Iy
	Coupled variable for the second moment of area in y direction, i.e., integral of y^2*dA over the cross-section. More...
const VariableValue &	_Iz
	Coupled variable for the second moment of area in z direction, i.e., integral of z^2*dA over the cross-section. More...
const VariableValue &	_Ix
	Coupled variable for the second moment of area in x direction, i.e., integral of (y^2 + z^2)*dA over the cross-section. More...
RankTwoTensor	_original_local_config
	Rotational transformation from global coordinate system to initial beam local configuration. More...
MaterialProperty< Real > &	_original_length
	Initial length of the beam. More...
MaterialProperty< RankTwoTensor > &	_total_rotation
	Rotational transformation from global coordinate system to beam local configuration at time t. More...
MaterialProperty< RealVectorValue > &	_total_disp_strain
	Current total displacement strain integrated over the cross-section in global coordinate system. More...
MaterialProperty< RealVectorValue > &	_total_rot_strain
	Current total rotational strain integrated over the cross-section in global coordinate system. More...
const MaterialProperty< RealVectorValue > &	_total_disp_strain_old
	Old total displacement strain integrated over the cross-section in global coordinate system. More...
const MaterialProperty< RealVectorValue > &	_total_rot_strain_old
	Old total rotational strain integrated over the cross-section in global coordinate system. More...
MaterialProperty< RealVectorValue > &	_mech_disp_strain_increment
	Mechanical displacement strain increment (after removal of eigenstrains) integrated over the cross-section. More...
MaterialProperty< RealVectorValue > &	_mech_rot_strain_increment
	Mechanical rotation strain increment (after removal of eigenstrains) integrated over the cross-section. More...
const MaterialProperty< RealVectorValue > &	_material_stiffness
	Material stiffness vector that relates displacement strain increments to force increments. More...
MaterialProperty< RankTwoTensor > &	_K11
	Stiffness matrix between displacement DOFs of same node or across nodes. More...
MaterialProperty< RankTwoTensor > &	_K21_cross
	Stiffness matrix between displacement DOFs of one node to rotational DOFs of another node. More...
MaterialProperty< RankTwoTensor > &	_K21
	Stiffness matrix between displacement DOFs and rotation DOFs of the same node. More...
MaterialProperty< RankTwoTensor > &	_K22
	Stiffness matrix between rotation DOFs of the same node. More...
MaterialProperty< RankTwoTensor > &	_K22_cross
	Stiffness matrix between rotation DOFs of different nodes. More...
const bool	_large_strain
	Boolean flag to turn on large strain calculation. More...
RealVectorValue	_grad_disp_0_local_t
	Gradient of displacement calculated in the beam local configuration at time t. More...
RealVectorValue	_grad_rot_0_local_t
	Gradient of rotation calculated in the beam local configuration at time t. More...
RealVectorValue	_avg_rot_local_t
	Average rotation calculated in the beam local configuration at time t. More...
std::vector< MaterialPropertyName >	_eigenstrain_names
	Vector of beam eigenstrain names. More...
std::vector< const MaterialProperty< RealVectorValue > * >	_disp_eigenstrain
	Vector of current displacement eigenstrains. More...
std::vector< const MaterialProperty< RealVectorValue > * >	_rot_eigenstrain
	Vector of current rotational eigenstrains. More...
std::vector< const MaterialProperty< RealVectorValue > * >	_disp_eigenstrain_old
	Vector of old displacement eigenstrains. More...
std::vector< const MaterialProperty< RealVectorValue > * >	_rot_eigenstrain_old
	Vector of old rotational eigenstrains. More...
RealVectorValue	_disp0
	Displacement and rotations at the two nodes of the beam in the global coordinate system. More...
RealVectorValue	_disp1
RealVectorValue	_rot0
RealVectorValue	_rot1
NonlinearSystemBase &	_nonlinear_sys
	Reference to the nonlinear system object. More...
std::vector< unsigned int >	_soln_disp_index_0
	Indices of solution vector corresponding to displacement DOFs at the node 0. More...
std::vector< unsigned int >	_soln_disp_index_1
	Indices of solution vector corresponding to displacement DOFs at the node 1. More...
std::vector< unsigned int >	_soln_rot_index_0
	Indices of solution vector corresponding to rotation DOFs at the node 0. More...
std::vector< unsigned int >	_soln_rot_index_1
	Indices of solution vector corresponding to rotation DOFs at the node 1. More...
MaterialProperty< RankTwoTensor > &	_initial_rotation
	Rotational transformation from global coordinate system to initial beam local configuration. More...
MaterialProperty< Real > &	_effective_stiffness
	Psuedo stiffness for critical time step computation. More...
const Function *const	_prefactor_function
	Prefactor function to multiply the elasticity tensor with. More...

Detailed Description

Definition at line 27 of file ComputeIncrementalBeamStrain.h.

Constructor & Destructor Documentation

ComputeIncrementalBeamStrain()

ComputeIncrementalBeamStrain::ComputeIncrementalBeamStrain

(

const InputParameters &

parameters

)

Definition at line 67 of file ComputeIncrementalBeamStrain.C.

  : Material(parameters),
    _has_Ix(isParamValid("Ix")),
    _nrot(coupledComponents("rotations")),
    _ndisp(coupledComponents("displacements")),
    _rot_num(_nrot),
    _disp_num(_ndisp),
    _area(coupledValue("area")),
    _Ay(coupledValue("Ay")),
    _Az(coupledValue("Az")),
    _Iy(coupledValue("Iy")),
    _Iz(coupledValue("Iz")),
    _Ix(_has_Ix ? coupledValue("Ix") : _zero),
    _original_length(declareProperty<Real>("original_length")),
    _total_rotation(declareProperty<RankTwoTensor>("total_rotation")),
    _total_disp_strain(declareProperty<RealVectorValue>("total_disp_strain")),
    _total_rot_strain(declareProperty<RealVectorValue>("total_rot_strain")),
    _total_disp_strain_old(getMaterialPropertyOld<RealVectorValue>("total_disp_strain")),
    _total_rot_strain_old(getMaterialPropertyOld<RealVectorValue>("total_rot_strain")),
    _mech_disp_strain_increment(declareProperty<RealVectorValue>("mech_disp_strain_increment")),
    _mech_rot_strain_increment(declareProperty<RealVectorValue>("mech_rot_strain_increment")),
    _material_stiffness(getMaterialPropertyByName<RealVectorValue>("material_stiffness")),
    _K11(declareProperty<RankTwoTensor>("Jacobian_11")),
    _K21_cross(declareProperty<RankTwoTensor>("Jacobian_12")),
    _K21(declareProperty<RankTwoTensor>("Jacobian_21")),
    _K22(declareProperty<RankTwoTensor>("Jacobian_22")),
    _K22_cross(declareProperty<RankTwoTensor>("Jacobian_22_cross")),
    _large_strain(getParam<bool>("large_strain")),
    _eigenstrain_names(getParam<std::vector<MaterialPropertyName>>("eigenstrain_names")),
    _disp_eigenstrain(_eigenstrain_names.size()),
    _rot_eigenstrain(_eigenstrain_names.size()),
    _disp_eigenstrain_old(_eigenstrain_names.size()),
    _rot_eigenstrain_old(_eigenstrain_names.size()),
    _nonlinear_sys(_fe_problem.getNonlinearSystemBase()),
    _soln_disp_index_0(_ndisp),
    _soln_disp_index_1(_ndisp),
    _soln_rot_index_0(_ndisp),
    _soln_rot_index_1(_ndisp),
    _initial_rotation(declareProperty<RankTwoTensor>("initial_rotation")),
    _effective_stiffness(declareProperty<Real>("effective_stiffness")),
    _prefactor_function(isParamValid("elasticity_prefactor") ? &getFunction("elasticity_prefactor")
                                                             : nullptr)
{
  // Checking for consistency between length of the provided displacements and rotations vector
  if (_ndisp != _nrot)
    mooseError("ComputeIncrementalBeamStrain: The number of variables supplied in 'displacements' "
               "and 'rotations' must match.");
 
  // fetch coupled variables and gradients (as stateful properties if necessary)
  for (unsigned int i = 0; i < _ndisp; ++i)
  {
    MooseVariable * disp_variable = getVar("displacements", i);
    _disp_num[i] = disp_variable->number();
 
    MooseVariable * rot_variable = getVar("rotations", i);
    _rot_num[i] = rot_variable->number();
  }
 
  if (_large_strain && (_Ay[0] > 0.0 || _Ay[1] > 0.0 || _Az[0] > 0.0 || _Az[1] > 0.0))
    mooseError("ComputeIncrementalBeamStrain: Large strain calculation does not currently "
               "support asymmetric beam configurations with non-zero first or third moments of "
               "area.");
 
  for (unsigned int i = 0; i < _eigenstrain_names.size(); ++i)
  {
    _disp_eigenstrain[i] = &getMaterialProperty<RealVectorValue>("disp_" + _eigenstrain_names[i]);
    _rot_eigenstrain[i] = &getMaterialProperty<RealVectorValue>("rot_" + _eigenstrain_names[i]);
    _disp_eigenstrain_old[i] =
        &getMaterialPropertyOld<RealVectorValue>("disp_" + _eigenstrain_names[i]);
    _rot_eigenstrain_old[i] =
        &getMaterialPropertyOld<RealVectorValue>("rot_" + _eigenstrain_names[i]);
  }
}

Member Function Documentation

computeProperties()

void ComputeIncrementalBeamStrain::computeProperties ( )

overridevirtual

Definition at line 184 of file ComputeIncrementalBeamStrain.C.

{
  // fetch the two end nodes for current element
  std::vector<const Node *> node;
  for (unsigned int i = 0; i < 2; ++i)
    node.push_back(_current_elem->node_ptr(i));
 
  // calculate original length of a beam element
  // Nodal positions do not change with time as undisplaced mesh is used by material classes by
  // default
  RealGradient dxyz;
  for (unsigned int i = 0; i < _ndisp; ++i)
    dxyz(i) = (*node[1])(i) - (*node[0])(i);
 
  _original_length[0] = dxyz.norm();
 
  // Fetch the solution for the two end nodes at time t
  const NumericVector<Number> & sol = *_nonlinear_sys.currentSolution();
  const NumericVector<Number> & sol_old = _nonlinear_sys.solutionOld();
 
  for (unsigned int i = 0; i < _ndisp; ++i)
  {
    _soln_disp_index_0[i] = node[0]->dof_number(_nonlinear_sys.number(), _disp_num[i], 0);
    _soln_disp_index_1[i] = node[1]->dof_number(_nonlinear_sys.number(), _disp_num[i], 0);
    _soln_rot_index_0[i] = node[0]->dof_number(_nonlinear_sys.number(), _rot_num[i], 0);
    _soln_rot_index_1[i] = node[1]->dof_number(_nonlinear_sys.number(), _rot_num[i], 0);
 
    _disp0(i) = sol(_soln_disp_index_0[i]) - sol_old(_soln_disp_index_0[i]);
    _disp1(i) = sol(_soln_disp_index_1[i]) - sol_old(_soln_disp_index_1[i]);
    _rot0(i) = sol(_soln_rot_index_0[i]) - sol_old(_soln_rot_index_0[i]);
    _rot1(i) = sol(_soln_rot_index_1[i]) - sol_old(_soln_rot_index_1[i]);
  }
 
  // For small rotation problems, the rotation matrix is essentially the transformation from the
  // global to original beam local configuration and is never updated. This method has to be
  // overriden for scenarios with finite rotation
  computeRotation();
  _initial_rotation[0] = _original_local_config;
 
  for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
    computeQpStrain();
 
  if (_fe_problem.currentlyComputingJacobian())
    computeStiffnessMatrix();
}

computeQpStrain()

void ComputeIncrementalBeamStrain::computeQpStrain ( )

protected

Computes the displacement and rotation strain increments.

Definition at line 231 of file ComputeIncrementalBeamStrain.C.

{
  const Real A_avg = (_area[0] + _area[1]) / 2.0;
  const Real Iz_avg = (_Iz[0] + _Iz[1]) / 2.0;
  Real Ix = _Ix[_qp];
  if (!_has_Ix)
    Ix = _Iy[_qp] + _Iz[_qp];
 
  // Rotate the gradient of displacements and rotations at t+delta t from global coordinate
  // frame to beam local coordinate frame
  const RealVectorValue grad_disp_0(1.0 / _original_length[0] * (_disp1 - _disp0));
  const RealVectorValue grad_rot_0(1.0 / _original_length[0] * (_rot1 - _rot0));
  const RealVectorValue avg_rot(
      0.5 * (_rot0(0) + _rot1(0)), 0.5 * (_rot0(1) + _rot1(1)), 0.5 * (_rot0(2) + _rot1(2)));
 
  _grad_disp_0_local_t = _total_rotation[0] * grad_disp_0;
  _grad_rot_0_local_t = _total_rotation[0] * grad_rot_0;
  _avg_rot_local_t = _total_rotation[0] * avg_rot;
 
  // displacement at any location on beam in local coordinate system at t
  // u_1 = u_n1 - rot_3 * y + rot_2 * z
  // u_2 = u_n2 - rot_1 * z
  // u_3 = u_n3 + rot_1 * y
  // where u_n1, u_n2, u_n3 are displacements at neutral axis
 
  // small strain
  // e_11 = u_1,1 = u_n1, 1 - rot_3, 1 * y + rot_2, 1 * z
  // e_12 = 2 * 0.5 * (u_1,2 + u_2,1) = (- rot_3 + u_n2,1 - rot_1,1 * z)
  // e_13 = 2 * 0.5 * (u_1,3 + u_3,1) = (rot_2 + u_n3,1 + rot_1,1 * y)
 
  // axial and shearing strains at each qp along the length of the beam
  _mech_disp_strain_increment[_qp](0) = _grad_disp_0_local_t(0) * _area[_qp] -
                                        _grad_rot_0_local_t(2) * _Ay[_qp] +
                                        _grad_rot_0_local_t(1) * _Az[_qp];
  _mech_disp_strain_increment[_qp](1) = -_avg_rot_local_t(2) * _area[_qp] +
                                        _grad_disp_0_local_t(1) * _area[_qp] -
                                        _grad_rot_0_local_t(0) * _Az[_qp];
  _mech_disp_strain_increment[_qp](2) = _avg_rot_local_t(1) * _area[_qp] +
                                        _grad_disp_0_local_t(2) * _area[_qp] +
                                        _grad_rot_0_local_t(0) * _Ay[_qp];
 
  // rotational strains at each qp along the length of the beam
  // rot_strain_1 = integral(e_13 * y - e_12 * z) dA
  // rot_strain_2 = integral(e_11 * z) dA
  // rot_strain_3 = integral(e_11 * y) dA
  // Iyz is the product moment of inertia which is zero for most cross-sections so it is assumed to
  // be zero for this analysis
  const Real Iyz = 0;
  _mech_rot_strain_increment[_qp](0) =
      _avg_rot_local_t(1) * _Ay[_qp] + _grad_disp_0_local_t(2) * _Ay[_qp] +
      _grad_rot_0_local_t(0) * Ix + _avg_rot_local_t(2) * _Az[_qp] -
      _grad_disp_0_local_t(1) * _Az[_qp];
  _mech_rot_strain_increment[_qp](1) = _grad_disp_0_local_t(0) * _Az[_qp] -
                                       _grad_rot_0_local_t(2) * Iyz +
                                       _grad_rot_0_local_t(1) * _Iz[_qp];
  _mech_rot_strain_increment[_qp](2) = _grad_disp_0_local_t(0) * _Ay[_qp] -
                                       _grad_rot_0_local_t(2) * _Iy[_qp] +
                                       _grad_rot_0_local_t(1) * Iyz;
 
  if (_large_strain)
  {
    _mech_disp_strain_increment[_qp](0) +=
        0.5 *
        ((Utility::pow<2>(_grad_disp_0_local_t(0)) + Utility::pow<2>(_grad_disp_0_local_t(1)) +
          Utility::pow<2>(_grad_disp_0_local_t(2))) *
             _area[_qp] +
         Utility::pow<2>(_grad_rot_0_local_t(2)) * _Iy[_qp] +
         Utility::pow<2>(_grad_rot_0_local_t(1)) * _Iz[_qp] +
         Utility::pow<2>(_grad_rot_0_local_t(0)) * Ix);
    _mech_disp_strain_increment[_qp](1) += (-_avg_rot_local_t(2) * _grad_disp_0_local_t(0) +
                                            _avg_rot_local_t(0) * _grad_disp_0_local_t(2)) *
                                           _area[_qp];
    _mech_disp_strain_increment[_qp](2) += (_avg_rot_local_t(1) * _grad_disp_0_local_t(0) -
                                            _avg_rot_local_t(0) * _grad_disp_0_local_t(1)) *
                                           _area[_qp];
 
    _mech_rot_strain_increment[_qp](0) += -_avg_rot_local_t(1) * _grad_rot_0_local_t(2) * _Iy[_qp] +
                                          _avg_rot_local_t(2) * _grad_rot_0_local_t(1) * _Iz[_qp];
    _mech_rot_strain_increment[_qp](1) += (_grad_disp_0_local_t(0) * _grad_rot_0_local_t(1) -
                                           _grad_disp_0_local_t(1) * _grad_rot_0_local_t(0)) *
                                          _Iz[_qp];
    _mech_rot_strain_increment[_qp](2) += (_grad_disp_0_local_t(2) * _grad_rot_0_local_t(0) -
                                           _grad_disp_0_local_t(0) * _grad_rot_0_local_t(2)) *
                                          _Iy[_qp];
  }
 
  _total_disp_strain[_qp] = _total_rotation[0].transpose() * _mech_disp_strain_increment[_qp] +
                            _total_disp_strain_old[_qp];
  _total_rot_strain[_qp] = _total_rotation[0].transpose() * _mech_rot_strain_increment[_qp] +
                           _total_disp_strain_old[_qp];
 
  // Convert eigenstrain increment from global to beam local coordinate system and remove eigen
  // strain increment
  for (unsigned int i = 0; i < _eigenstrain_names.size(); ++i)
  {
    _mech_disp_strain_increment[_qp] -=
        _total_rotation[0] * ((*_disp_eigenstrain[i])[_qp] - (*_disp_eigenstrain_old[i])[_qp]) *
        _area[_qp];
    _mech_rot_strain_increment[_qp] -=
        _total_rotation[0] * ((*_rot_eigenstrain[i])[_qp] - (*_rot_eigenstrain_old[i])[_qp]);
  }
 
  Real c1_paper = std::sqrt(_material_stiffness[0](0));
  Real c2_paper = std::sqrt(_material_stiffness[0](1));
 
  Real effec_stiff_1 = std::max(c1_paper, c2_paper);
 
  Real effec_stiff_2 = 2 / (c2_paper * std::sqrt(A_avg / Iz_avg));
 
  _effective_stiffness[_qp] = std::max(effec_stiff_1, _original_length[0] / effec_stiff_2);
 
  if (_prefactor_function)
    _effective_stiffness[_qp] *= std::sqrt(_prefactor_function->value(_t, _q_point[_qp]));
}

Referenced by computeProperties().

computeRotation()

void ComputeIncrementalBeamStrain::computeRotation ( )

protectedvirtual

Computes the rotation matrix at time t. For small rotation scenarios, the rotation matrix at time t is same as the intiial rotation matrix.

Reimplemented in ComputeFiniteBeamStrain.

Definition at line 616 of file ComputeIncrementalBeamStrain.C.

{
  _total_rotation[0] = _original_local_config;
}

Referenced by computeProperties().

computeStiffnessMatrix()

void ComputeIncrementalBeamStrain::computeStiffnessMatrix ( )

protected

Computes the stiffness matrices.

Definition at line 347 of file ComputeIncrementalBeamStrain.C.

{
  const Real youngs_modulus = _material_stiffness[0](0);
  const Real shear_modulus = _material_stiffness[0](1);
 
  const Real A_avg = (_area[0] + _area[1]) / 2.0;
  const Real Iy_avg = (_Iy[0] + _Iy[1]) / 2.0;
  const Real Iz_avg = (_Iz[0] + _Iz[1]) / 2.0;
  Real Ix_avg = (_Ix[0] + _Ix[1]) / 2.0;
  if (!_has_Ix)
    Ix_avg = Iy_avg + Iz_avg;
 
  // K = |K11 K12|
  //     |K21 K22|
 
  // relation between translational displacements at node 0 and translational forces at node 0
  RankTwoTensor K11_local;
  K11_local.zero();
  K11_local(0, 0) = youngs_modulus * A_avg / _original_length[0];
  K11_local(1, 1) = shear_modulus * A_avg / _original_length[0];
  K11_local(2, 2) = shear_modulus * A_avg / _original_length[0];
  _K11[0] = _total_rotation[0].transpose() * K11_local * _total_rotation[0];
 
  // relation between displacements at node 0 and rotational moments at node 0
  RankTwoTensor K21_local;
  K21_local.zero();
  K21_local(2, 1) = shear_modulus * A_avg * 0.5;
  K21_local(1, 2) = -shear_modulus * A_avg * 0.5;
  _K21[0] = _total_rotation[0].transpose() * K21_local * _total_rotation[0];
 
  // relation between rotations at node 0 and rotational moments at node 0
  RankTwoTensor K22_local;
  K22_local.zero();
  K22_local(0, 0) = shear_modulus * Ix_avg / _original_length[0];
  K22_local(1, 1) = youngs_modulus * Iz_avg / _original_length[0] +
                    shear_modulus * A_avg * _original_length[0] / 4.0;
  K22_local(2, 2) = youngs_modulus * Iy_avg / _original_length[0] +
                    shear_modulus * A_avg * _original_length[0] / 4.0;
  _K22[0] = _total_rotation[0].transpose() * K22_local * _total_rotation[0];
 
  // relation between rotations at node 0 and rotational moments at node 1
  RankTwoTensor K22_local_cross = -K22_local;
  K22_local_cross(1, 1) += 2.0 * shear_modulus * A_avg * _original_length[0] / 4.0;
  K22_local_cross(2, 2) += 2.0 * shear_modulus * A_avg * _original_length[0] / 4.0;
  _K22_cross[0] = _total_rotation[0].transpose() * K22_local_cross * _total_rotation[0];
 
  // relation between displacements at node 0 and rotational moments at node 1
  _K21_cross[0] = -_K21[0];
 
  // stiffness matrix for large strain
  if (_large_strain)
  {
    // k1_large is the stiffness matrix obtained from sigma_xx * d(epsilon_xx)
    RankTwoTensor k1_large_11;
    // row 1
    k1_large_11(0, 0) = Utility::pow<2>(_grad_disp_0_local_t(0)) +
                        1.5 * Utility::pow<2>(_grad_rot_0_local_t(2)) * Iy_avg +
                        1.5 * Utility::pow<2>(_grad_rot_0_local_t(1)) * Iz_avg +
                        0.5 * Utility::pow<2>(_grad_disp_0_local_t(1)) +
                        0.5 * Utility::pow<2>(_grad_disp_0_local_t(2)) +
                        0.5 * Utility::pow<2>(_grad_rot_0_local_t(0)) * Ix_avg;
    k1_large_11(1, 0) = 0.5 * _grad_disp_0_local_t(0) * _grad_disp_0_local_t(1) -
                        1.0 / 3.0 * _grad_rot_0_local_t(0) * _grad_rot_0_local_t(1) * Iz_avg;
    k1_large_11(2, 0) = 0.5 * _grad_disp_0_local_t(0) * _grad_disp_0_local_t(2) -
                        1.0 / 3.0 * _grad_rot_0_local_t(0) * _grad_rot_0_local_t(2) * Iy_avg;
 
    // row 2
    k1_large_11(0, 1) = k1_large_11(1, 0);
    k1_large_11(1, 1) = Utility::pow<2>(_grad_disp_0_local_t(1)) +
                        1.5 * Utility::pow<2>(_grad_rot_0_local_t(0)) * Iz_avg +
                        0.5 * Utility::pow<2>(_grad_disp_0_local_t(0)) +
                        0.5 * Utility::pow<2>(_grad_disp_0_local_t(2)) +
                        0.5 * Utility::pow<2>(_grad_rot_0_local_t(2)) * Iy_avg +
                        0.5 * Utility::pow<2>(_grad_rot_0_local_t(1)) * Iz_avg +
                        0.5 * Utility::pow<2>(_grad_rot_0_local_t(0)) * Iy_avg;
    k1_large_11(2, 1) = 0.5 * _grad_disp_0_local_t(1) * _grad_disp_0_local_t(2);
 
    // row 3
    k1_large_11(0, 2) = k1_large_11(2, 0);
    k1_large_11(1, 2) = k1_large_11(2, 1);
    k1_large_11(2, 2) = Utility::pow<2>(_grad_disp_0_local_t(2)) +
                        1.5 * Utility::pow<2>(_grad_rot_0_local_t(0)) * Iy_avg +
                        0.5 * Utility::pow<2>(_grad_disp_0_local_t(0)) +
                        0.5 * Utility::pow<2>(_grad_disp_0_local_t(1)) +
                        0.5 * Utility::pow<2>(_grad_rot_0_local_t(0)) * Iz_avg +
                        0.5 * Utility::pow<2>(_grad_rot_0_local_t(2)) * Iy_avg +
                        0.5 * Utility::pow<2>(_grad_rot_0_local_t(2)) * Iz_avg;
 
    k1_large_11 *= 1.0 / 4.0 / Utility::pow<2>(_original_length[0]);
 
    RankTwoTensor k1_large_21;
    // row 1
    k1_large_21(0, 0) = 0.5 * _grad_disp_0_local_t(0) * _grad_rot_0_local_t(0) * (Ix_avg)-1.0 /
                            3.0 * _grad_disp_0_local_t(1) * _grad_rot_0_local_t(1) * Iz_avg -
                        1.0 / 3.0 * _grad_disp_0_local_t(2) * _grad_rot_0_local_t(2);
    k1_large_21(1, 0) = 1.5 * _grad_disp_0_local_t(0) * _grad_rot_0_local_t(1) * Iz_avg -
                        1.0 / 3.0 * _grad_disp_0_local_t(1) * _grad_rot_0_local_t(0) * Iz_avg;
    k1_large_21(2, 0) = 1.5 * _grad_disp_0_local_t(0) * _grad_rot_0_local_t(2) * Iy_avg -
                        1.0 / 3.0 * _grad_disp_0_local_t(2) * _grad_rot_0_local_t(0) * Iy_avg;
 
    // row 2
    k1_large_21(0, 1) = k1_large_21(1, 0);
    k1_large_21(1, 1) = 0.5 * _grad_disp_0_local_t(1) * _grad_rot_0_local_t(1) * Iz_avg -
                        1.0 / 3.0 * _grad_disp_0_local_t(0) * _grad_rot_0_local_t(0) * Iz_avg;
    k1_large_21(2, 1) = 0.5 * _grad_disp_0_local_t(1) * _grad_rot_0_local_t(2) * Iy_avg;
 
    // row 3
    k1_large_21(0, 2) = k1_large_21(2, 0);
    k1_large_21(1, 2) = k1_large_21(2, 1);
    k1_large_21(2, 2) = 0.5 * _grad_disp_0_local_t(2) * _grad_rot_0_local_t(2) * Iy_avg -
                        1.0 / 3.0 * _grad_disp_0_local_t(0) * _grad_rot_0_local_t(0) * Iy_avg;
    k1_large_21 *= 1.0 / 4.0 / Utility::pow<2>(_original_length[0]);
 
    RankTwoTensor k1_large_22;
    // row 1
    k1_large_22(0, 0) = Utility::pow<2>(_grad_rot_0_local_t(0)) * Utility::pow<2>(Ix_avg) +
                        1.5 * Utility::pow<2>(_grad_disp_0_local_t(1)) * Iz_avg +
                        1.5 * Utility::pow<2>(_grad_disp_0_local_t(2)) * Iy_avg +
                        0.5 * Utility::pow<2>(_grad_disp_0_local_t(0)) * Ix_avg +
                        0.5 * Utility::pow<2>(_grad_disp_0_local_t(2)) * Iz_avg +
                        0.5 * Utility::pow<2>(_grad_disp_0_local_t(1)) * Iy_avg +
                        0.5 * Utility::pow<2>(_grad_rot_0_local_t(2)) * Iy_avg * Ix_avg +
                        0.5 * Utility::pow<2>(_grad_rot_0_local_t(1)) * Iz_avg * Ix_avg;
    k1_large_22(1, 0) = 0.5 * _grad_rot_0_local_t(0) * _grad_rot_0_local_t(1) * Iz_avg * Ix_avg -
                        1.0 / 3.0 * _grad_disp_0_local_t(0) * _grad_disp_0_local_t(1) * Iz_avg;
    k1_large_22(2, 0) = 0.5 * _grad_rot_0_local_t(0) * _grad_rot_0_local_t(2) * Iy_avg * Ix_avg -
                        1.0 / 3.0 * _grad_disp_0_local_t(0) * _grad_disp_0_local_t(2) * Iy_avg;
 
    // row 2
    k1_large_22(0, 1) = k1_large_22(1, 0);
    k1_large_22(1, 1) = Utility::pow<2>(_grad_rot_0_local_t(1)) * Iz_avg * Iz_avg +
                        1.5 * Utility::pow<2>(_grad_disp_0_local_t(0)) * Iz_avg +
                        1.5 * Utility::pow<2>(_grad_rot_0_local_t(2)) * Iy_avg * Iz_avg +
                        0.5 * Utility::pow<2>(_grad_disp_0_local_t(1)) * Iz_avg +
                        0.5 * Utility::pow<2>(_grad_disp_0_local_t(2)) * Iz_avg +
                        0.5 * Utility::pow<2>(_grad_rot_0_local_t(0)) * Iz_avg * Ix_avg;
    k1_large_22(2, 1) = 1.5 * _grad_rot_0_local_t(1) * _grad_rot_0_local_t(2) * Iy_avg * Iz_avg;
 
    // row 3
    k1_large_22(0, 2) = k1_large_22(2, 0);
    k1_large_22(1, 2) = k1_large_22(2, 1);
    k1_large_22(2, 2) = Utility::pow<2>(_grad_rot_0_local_t(2)) * Iy_avg * Iy_avg +
                        1.5 * Utility::pow<2>(_grad_disp_0_local_t(0)) * Iy_avg +
                        1.5 * Utility::pow<2>(_grad_rot_0_local_t(1)) * Iy_avg * Iz_avg +
                        0.5 * Utility::pow<2>(_grad_disp_0_local_t(1)) * Iy_avg +
                        0.5 * Utility::pow<2>(_grad_disp_0_local_t(2)) * Iy_avg +
                        0.5 * Utility::pow<2>(_grad_rot_0_local_t(0)) * Iz_avg * Ix_avg;
 
    k1_large_22 *= 1.0 / 4.0 / Utility::pow<2>(_original_length[0]);
 
    // k2_large and k3_large are contributions from tau_xy * d(gamma_xy) and tau_xz * d(gamma_xz)
    // k2_large for node 1 is negative of that for node 0
    RankTwoTensor k2_large_11;
    // col 1
    k2_large_11(0, 0) =
        0.25 * Utility::pow<2>(_avg_rot_local_t(2)) + 0.25 * Utility::pow<2>(_avg_rot_local_t(1));
    k2_large_11(1, 0) = -1.0 / 6.0 * _avg_rot_local_t(0) * _avg_rot_local_t(1);
    k2_large_11(2, 0) = -1.0 / 6.0 * _avg_rot_local_t(0) * _avg_rot_local_t(2);
 
    // col 2
    k2_large_11(0, 1) = k2_large_11(1, 0);
    k2_large_11(1, 1) = 0.25 * _avg_rot_local_t(0);
 
    // col 3
    k2_large_11(0, 2) = k2_large_11(2, 0);
    k2_large_11(2, 2) = 0.25 * Utility::pow<2>(_avg_rot_local_t(0));
 
    k2_large_11 *= 1.0 / 4.0 / Utility::pow<2>(_original_length[0]);
 
    RankTwoTensor k2_large_22;
    // col1
    k2_large_22(0, 0) = 0.25 * Utility::pow<2>(_avg_rot_local_t(0)) * Ix_avg;
    k2_large_22(1, 0) = 1.0 / 6.0 * _avg_rot_local_t(0) * _avg_rot_local_t(1) * Iz_avg;
    k2_large_22(2, 0) = 1.0 / 6.0 * _avg_rot_local_t(0) * _avg_rot_local_t(2) * Iy_avg;
 
    // col2
    k2_large_22(0, 1) = k2_large_22(1, 0);
    k2_large_22(1, 1) = 0.25 * Utility::pow<2>(_avg_rot_local_t(2)) * Iz_avg +
                        0.25 * Utility::pow<2>(_avg_rot_local_t(1)) * Iz_avg;
 
    // col3
    k2_large_22(0, 2) = k2_large_22(2, 0);
    k2_large_22(2, 2) = 0.25 * Utility::pow<2>(_avg_rot_local_t(2)) * Iy_avg +
                        0.25 * Utility::pow<2>(_avg_rot_local_t(1)) * Iy_avg;
 
    k2_large_22 *= 1.0 / 4.0 / Utility::pow<2>(_original_length[0]);
 
    // k3_large for node 1 is same as that for node 0
    RankTwoTensor k3_large_22;
    // col1
    k3_large_22(0, 0) = 0.25 * Utility::pow<2>(_grad_disp_0_local_t(2)) +
                        0.25 * _grad_rot_0_local_t(0) * Ix_avg +
                        0.25 * Utility::pow<2>(_grad_disp_0_local_t(1));
    k3_large_22(1, 0) = -1.0 / 6.0 * _grad_disp_0_local_t(0) * _grad_disp_0_local_t(1) +
                        1.0 / 6.0 * _grad_rot_0_local_t(0) * _grad_rot_0_local_t(1) * Iz_avg;
    k3_large_22(2, 0) = -1.0 / 6.0 * _grad_disp_0_local_t(0) * _grad_disp_0_local_t(2) +
                        1.0 / 6.0 * _grad_rot_0_local_t(0) * _grad_rot_0_local_t(2) * Iy_avg;
 
    // col2
    k3_large_22(0, 1) = k3_large_22(1, 0);
    k3_large_22(2, 2) = 0.25 * Utility::pow<2>(_grad_disp_0_local_t(0)) +
                        0.25 * _grad_rot_0_local_t(2) * Iy_avg +
                        0.25 * _grad_rot_0_local_t(1) * Iz_avg;
 
    // col3
    k3_large_22(0, 2) = k3_large_22(2, 0);
    k3_large_22(2, 2) = 0.25 * Utility::pow<2>(_grad_disp_0_local_t(0)) +
                        0.25 * _grad_rot_0_local_t(2) * Iy_avg +
                        0.25 * _grad_rot_0_local_t(1) * Iz_avg;
 
    k3_large_22 *= 1.0 / 16.0;
 
    RankTwoTensor k3_large_21;
    // col1
    k3_large_21(0, 0) = -1.0 / 6.0 *
                        (_grad_disp_0_local_t(2) * _avg_rot_local_t(2) +
                         _grad_disp_0_local_t(1) * _avg_rot_local_t(1));
    k3_large_21(1, 0) = 0.25 * _grad_disp_0_local_t(0) * _avg_rot_local_t(1) -
                        1.0 / 6.0 * _grad_disp_0_local_t(1) * _avg_rot_local_t(0);
    k3_large_21(2, 0) = 0.25 * _grad_disp_0_local_t(0) * _avg_rot_local_t(2) -
                        1.0 / 6.0 * _grad_disp_0_local_t(2) * _avg_rot_local_t(0);
 
    // col2
    k3_large_21(0, 1) = 0.25 * _grad_disp_0_local_t(1) * _avg_rot_local_t(0) -
                        1.0 / 6.0 * _grad_disp_0_local_t(0) * _avg_rot_local_t(1);
    k3_large_21(1, 1) = -1.0 / 6.0 * _grad_disp_0_local_t(0) * _avg_rot_local_t(0);
 
    // col3
    k3_large_21(0, 2) = 0.25 * _grad_disp_0_local_t(2) * _avg_rot_local_t(0) -
                        1.0 / 6.0 * _grad_disp_0_local_t(0) * _avg_rot_local_t(2);
    k3_large_21(2, 2) = -1.0 / 6.0 * _grad_disp_0_local_t(0) * _avg_rot_local_t(0);
 
    k3_large_21 *= 1.0 / 8.0 / _original_length[0];
 
    RankTwoTensor k4_large_22;
    // col 1
    k4_large_22(0, 0) = 0.25 * _grad_rot_0_local_t(0) * _avg_rot_local_t(0) * Ix_avg +
                        1.0 / 6.0 * _grad_rot_0_local_t(2) * _avg_rot_local_t(2) * Iy_avg +
                        1.0 / 6.0 * _grad_rot_0_local_t(1) * _avg_rot_local_t(1) * Iz_avg;
    k4_large_22(1, 0) = 1.0 / 6.0 * _grad_rot_0_local_t(1) * _avg_rot_local_t(0) * Iz_avg;
    k4_large_22(2, 0) = 1.0 / 6.0 * _grad_rot_0_local_t(2) * _avg_rot_local_t(0) * Iy_avg;
 
    // col2
    k4_large_22(0, 1) = 1.0 / 6.0 * _grad_rot_0_local_t(0) * _avg_rot_local_t(1) * Iz_avg;
    k4_large_22(1, 1) = 0.25 * _grad_rot_0_local_t(1) * _avg_rot_local_t(1) * Iz_avg +
                        1.0 / 6.0 * _grad_rot_0_local_t(0) * _avg_rot_local_t(0) * Iz_avg;
    k4_large_22(2, 1) = 0.25 * _grad_rot_0_local_t(1) * _avg_rot_local_t(2) * Iz_avg;
 
    // col 3
    k4_large_22(0, 2) = 1.0 / 6.0 * _grad_rot_0_local_t(0) * _avg_rot_local_t(2) * Iy_avg;
    k4_large_22(1, 2) = 0.25 * _grad_rot_0_local_t(2) * _avg_rot_local_t(1) * Iy_avg;
    k4_large_22(2, 2) = 0.25 * _grad_rot_0_local_t(2) * _avg_rot_local_t(2) * Iy_avg +
                        1.0 / 6.0 * _grad_rot_0_local_t(0) * _avg_rot_local_t(0) * Iy_avg;
 
    k3_large_22 += 1.0 / 8.0 / _original_length[0] * (k4_large_22 + k4_large_22.transpose());
 
    // Assembling final matrix
    _K11[0] += _total_rotation[0].transpose() * (k1_large_11 + k2_large_11) * _total_rotation[0];
    _K22[0] += _total_rotation[0].transpose() * (k1_large_22 + k2_large_22 + k3_large_22) *
               _total_rotation[0];
    _K21[0] += _total_rotation[0].transpose() * (k1_large_21 + k3_large_21) * _total_rotation[0];
    _K21_cross[0] +=
        _total_rotation[0].transpose() * (-k1_large_21 + k3_large_21) * _total_rotation[0];
    _K22_cross[0] += _total_rotation[0].transpose() * (-k1_large_22 - k2_large_22 + k3_large_22) *
                     _total_rotation[0];
  }
}

Referenced by computeProperties().

initQpStatefulProperties()

void ComputeIncrementalBeamStrain::initQpStatefulProperties ( )

overrideprotectedvirtual

Definition at line 142 of file ComputeIncrementalBeamStrain.C.

{
  // compute initial orientation of the beam for calculating initial rotation matrix
  const std::vector<RealGradient> * orientation =
      &_subproblem.assembly(_tid).getFE(FEType(), 1)->get_dxyzdxi();
  RealGradient x_orientation = (*orientation)[0];
  x_orientation /= x_orientation.norm();
 
  RealGradient y_orientation = getParam<RealGradient>("y_orientation");
  y_orientation /= y_orientation.norm();
  Real sum = x_orientation(0) * y_orientation(0) + x_orientation(1) * y_orientation(1) +
             x_orientation(2) * y_orientation(2);
 
  if (std::abs(sum) > 1e-4)
    mooseError("ComputeIncrementalBeamStrain: y_orientation should be perpendicular to "
               "the axis of the beam.");
 
  // Calculate z orientation as a cross product of the x and y orientations
  RealGradient z_orientation;
  z_orientation(0) = (x_orientation(1) * y_orientation(2) - x_orientation(2) * y_orientation(1));
  z_orientation(1) = (x_orientation(2) * y_orientation(0) - x_orientation(0) * y_orientation(2));
  z_orientation(2) = (x_orientation(0) * y_orientation(1) - x_orientation(1) * y_orientation(0));
 
  // Rotation matrix from global to original beam local configuration
  _original_local_config(0, 0) = x_orientation(0);
  _original_local_config(0, 1) = x_orientation(1);
  _original_local_config(0, 2) = x_orientation(2);
  _original_local_config(1, 0) = y_orientation(0);
  _original_local_config(1, 1) = y_orientation(1);
  _original_local_config(1, 2) = y_orientation(2);
  _original_local_config(2, 0) = z_orientation(0);
  _original_local_config(2, 1) = z_orientation(1);
  _original_local_config(2, 2) = z_orientation(2);
 
  _total_rotation[_qp] = _original_local_config;
 
  RealVectorValue temp;
  _total_disp_strain[_qp] = temp;
  _total_rot_strain[_qp] = temp;
}

validParams()

InputParameters ComputeIncrementalBeamStrain::validParams ( )

static

Definition at line 25 of file ComputeIncrementalBeamStrain.C.

{
  InputParameters params = Material::validParams();
  params.addClassDescription("Compute a infinitesimal/large strain increment for the beam.");
  params.addRequiredCoupledVar(
      "rotations", "The rotations appropriate for the simulation geometry and coordinate system");
  params.addRequiredCoupledVar(
      "displacements",
      "The displacements appropriate for the simulation geometry and coordinate system");
  params.addRequiredParam<RealGradient>("y_orientation",
                                        "Orientation of the y direction along "
                                        "with Iyy is provided. This should be "
                                        "perpendicular to the axis of the beam.");
  params.addRequiredCoupledVar(
      "area",
      "Cross-section area of the beam. Can be supplied as either a number or a variable name.");
  params.addCoupledVar("Ay",
                       0.0,
                       "First moment of area of the beam about y axis. Can be supplied "
                       "as either a number or a variable name.");
  params.addCoupledVar("Az",
                       0.0,
                       "First moment of area of the beam about z axis. Can be supplied "
                       "as either a number or a variable name.");
  params.addCoupledVar("Ix",
                       "Second moment of area of the beam about x axis. Can be "
                       "supplied as either a number or a variable name. Defaults to Iy+Iz.");
  params.addRequiredCoupledVar("Iy",
                               "Second moment of area of the beam about y axis. Can be "
                               "supplied as either a number or a variable name.");
  params.addRequiredCoupledVar("Iz",
                               "Second moment of area of the beam about z axis. Can be "
                               "supplied as either a number or a variable name.");
  params.addParam<bool>("large_strain", false, "Set to true if large strain are to be calculated.");
  params.addParam<std::vector<MaterialPropertyName>>(
      "eigenstrain_names", "List of beam eigenstrains to be applied in this strain calculation.");
  params.addParam<FunctionName>(
      "elasticity_prefactor",
      "Optional function to use as a scalar prefactor on the elasticity vector for the beam.");
  return params;
}

Referenced by ComputeFiniteBeamStrain::validParams().

Member Data Documentation

_area

const VariableValue& ComputeIncrementalBeamStrain::_area

protected

Coupled variable for the beam cross-sectional area.

Definition at line 64 of file ComputeIncrementalBeamStrain.h.

Referenced by computeQpStrain(), and computeStiffnessMatrix().

_avg_rot_local_t

RealVectorValue ComputeIncrementalBeamStrain::_avg_rot_local_t

protected

Average rotation calculated in the beam local configuration at time t.

Definition at line 136 of file ComputeIncrementalBeamStrain.h.

Referenced by computeQpStrain(), and computeStiffnessMatrix().

_Ay

const VariableValue& ComputeIncrementalBeamStrain::_Ay

protected

Coupled variable for the first moment of area in y direction, i.e., integral of y*dA over the cross-section.

Definition at line 67 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain(), and computeQpStrain().

_Az

const VariableValue& ComputeIncrementalBeamStrain::_Az

protected

Coupled variable for the first moment of area in z direction, i.e., integral of z*dA over the cross-section.

Definition at line 70 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain(), and computeQpStrain().

_disp0

RealVectorValue ComputeIncrementalBeamStrain::_disp0

protected

Displacement and rotations at the two nodes of the beam in the global coordinate system.

Definition at line 154 of file ComputeIncrementalBeamStrain.h.

Referenced by computeProperties(), computeQpStrain(), and ComputeFiniteBeamStrain::computeRotation().

_disp1

RealVectorValue ComputeIncrementalBeamStrain::_disp1

protected

Definition at line 154 of file ComputeIncrementalBeamStrain.h.

Referenced by computeProperties(), computeQpStrain(), and ComputeFiniteBeamStrain::computeRotation().

_disp_eigenstrain

std::vector<const MaterialProperty<RealVectorValue> *> ComputeIncrementalBeamStrain::_disp_eigenstrain

protected

Vector of current displacement eigenstrains.

Definition at line 142 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain(), and computeQpStrain().

_disp_eigenstrain_old

std::vector<const MaterialProperty<RealVectorValue> *> ComputeIncrementalBeamStrain::_disp_eigenstrain_old

protected

Vector of old displacement eigenstrains.

Definition at line 148 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain(), and computeQpStrain().

_disp_num

std::vector<unsigned int> ComputeIncrementalBeamStrain::_disp_num

protected

Variable numbers corresponding to the displacement variables.

Definition at line 61 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain(), and computeProperties().

_effective_stiffness

MaterialProperty<Real>& ComputeIncrementalBeamStrain::_effective_stiffness

protected

Psuedo stiffness for critical time step computation.

Definition at line 174 of file ComputeIncrementalBeamStrain.h.

Referenced by computeQpStrain().

_eigenstrain_names

std::vector<MaterialPropertyName> ComputeIncrementalBeamStrain::_eigenstrain_names

protected

Vector of beam eigenstrain names.

Definition at line 139 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain(), and computeQpStrain().

_grad_disp_0_local_t

RealVectorValue ComputeIncrementalBeamStrain::_grad_disp_0_local_t

protected

Gradient of displacement calculated in the beam local configuration at time t.

Definition at line 130 of file ComputeIncrementalBeamStrain.h.

Referenced by computeQpStrain(), and computeStiffnessMatrix().

_grad_rot_0_local_t

RealVectorValue ComputeIncrementalBeamStrain::_grad_rot_0_local_t

protected

Gradient of rotation calculated in the beam local configuration at time t.

Definition at line 133 of file ComputeIncrementalBeamStrain.h.

Referenced by computeQpStrain(), and computeStiffnessMatrix().

_has_Ix

const bool ComputeIncrementalBeamStrain::_has_Ix

protected

Booleans for validity of params.

Definition at line 49 of file ComputeIncrementalBeamStrain.h.

Referenced by computeQpStrain(), and computeStiffnessMatrix().

_initial_rotation

MaterialProperty<RankTwoTensor>& ComputeIncrementalBeamStrain::_initial_rotation

protected

Rotational transformation from global coordinate system to initial beam local configuration.

Definition at line 171 of file ComputeIncrementalBeamStrain.h.

Referenced by computeProperties().

_Ix

const VariableValue& ComputeIncrementalBeamStrain::_Ix

protected

Coupled variable for the second moment of area in x direction, i.e., integral of (y^2 + z^2)*dA over the cross-section.

Definition at line 79 of file ComputeIncrementalBeamStrain.h.

Referenced by computeQpStrain(), and computeStiffnessMatrix().

_Iy

const VariableValue& ComputeIncrementalBeamStrain::_Iy

protected

Coupled variable for the second moment of area in y direction, i.e., integral of y^2*dA over the cross-section.

Definition at line 73 of file ComputeIncrementalBeamStrain.h.

Referenced by computeQpStrain(), and computeStiffnessMatrix().

_Iz

const VariableValue& ComputeIncrementalBeamStrain::_Iz

protected

Coupled variable for the second moment of area in z direction, i.e., integral of z^2*dA over the cross-section.

Definition at line 76 of file ComputeIncrementalBeamStrain.h.

Referenced by computeQpStrain(), and computeStiffnessMatrix().

_K11

MaterialProperty<RankTwoTensor>& ComputeIncrementalBeamStrain::_K11

protected

Stiffness matrix between displacement DOFs of same node or across nodes.

Definition at line 112 of file ComputeIncrementalBeamStrain.h.

Referenced by computeStiffnessMatrix().

_K21

MaterialProperty<RankTwoTensor>& ComputeIncrementalBeamStrain::_K21

protected

Stiffness matrix between displacement DOFs and rotation DOFs of the same node.

Definition at line 118 of file ComputeIncrementalBeamStrain.h.

Referenced by computeStiffnessMatrix().

_K21_cross

MaterialProperty<RankTwoTensor>& ComputeIncrementalBeamStrain::_K21_cross

protected

Stiffness matrix between displacement DOFs of one node to rotational DOFs of another node.

Definition at line 115 of file ComputeIncrementalBeamStrain.h.

Referenced by computeStiffnessMatrix().

_K22

MaterialProperty<RankTwoTensor>& ComputeIncrementalBeamStrain::_K22

protected

Stiffness matrix between rotation DOFs of the same node.

Definition at line 121 of file ComputeIncrementalBeamStrain.h.

Referenced by computeStiffnessMatrix().

_K22_cross

MaterialProperty<RankTwoTensor>& ComputeIncrementalBeamStrain::_K22_cross

protected

Stiffness matrix between rotation DOFs of different nodes.

Definition at line 124 of file ComputeIncrementalBeamStrain.h.

Referenced by computeStiffnessMatrix().

_large_strain

const bool ComputeIncrementalBeamStrain::_large_strain

protected

Boolean flag to turn on large strain calculation.

Definition at line 127 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain(), computeQpStrain(), and computeStiffnessMatrix().

_material_stiffness

const MaterialProperty<RealVectorValue>& ComputeIncrementalBeamStrain::_material_stiffness

protected

Material stiffness vector that relates displacement strain increments to force increments.

Definition at line 109 of file ComputeIncrementalBeamStrain.h.

Referenced by computeQpStrain(), and computeStiffnessMatrix().

_mech_disp_strain_increment

MaterialProperty<RealVectorValue>& ComputeIncrementalBeamStrain::_mech_disp_strain_increment

protected

Mechanical displacement strain increment (after removal of eigenstrains) integrated over the cross-section.

Definition at line 103 of file ComputeIncrementalBeamStrain.h.

Referenced by computeQpStrain().

_mech_rot_strain_increment

MaterialProperty<RealVectorValue>& ComputeIncrementalBeamStrain::_mech_rot_strain_increment

protected

Mechanical rotation strain increment (after removal of eigenstrains) integrated over the cross-section.

Definition at line 106 of file ComputeIncrementalBeamStrain.h.

Referenced by computeQpStrain().

_ndisp

unsigned int ComputeIncrementalBeamStrain::_ndisp

protected

Number of coupled displacement variables.

Definition at line 55 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain(), and computeProperties().

_nonlinear_sys

NonlinearSystemBase& ComputeIncrementalBeamStrain::_nonlinear_sys

protected

Reference to the nonlinear system object.

Definition at line 157 of file ComputeIncrementalBeamStrain.h.

Referenced by computeProperties().

_nrot

unsigned int ComputeIncrementalBeamStrain::_nrot

protected

Number of coupled rotational variables.

Definition at line 52 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain().

_original_length

MaterialProperty<Real>& ComputeIncrementalBeamStrain::_original_length

protected

Initial length of the beam.

Definition at line 85 of file ComputeIncrementalBeamStrain.h.

Referenced by computeProperties(), computeQpStrain(), ComputeFiniteBeamStrain::computeRotation(), and computeStiffnessMatrix().

_original_local_config

RankTwoTensor ComputeIncrementalBeamStrain::_original_local_config

protected

Rotational transformation from global coordinate system to initial beam local configuration.

Definition at line 82 of file ComputeIncrementalBeamStrain.h.

Referenced by computeProperties(), computeRotation(), and initQpStatefulProperties().

_prefactor_function

const Function* const ComputeIncrementalBeamStrain::_prefactor_function

protected

Prefactor function to multiply the elasticity tensor with.

Definition at line 177 of file ComputeIncrementalBeamStrain.h.

Referenced by computeQpStrain().

_rot0

RealVectorValue ComputeIncrementalBeamStrain::_rot0

protected

Definition at line 154 of file ComputeIncrementalBeamStrain.h.

Referenced by computeProperties(), computeQpStrain(), and ComputeFiniteBeamStrain::computeRotation().

_rot1

RealVectorValue ComputeIncrementalBeamStrain::_rot1

protected

Definition at line 154 of file ComputeIncrementalBeamStrain.h.

Referenced by computeProperties(), computeQpStrain(), and ComputeFiniteBeamStrain::computeRotation().

_rot_eigenstrain

std::vector<const MaterialProperty<RealVectorValue> *> ComputeIncrementalBeamStrain::_rot_eigenstrain

protected

Vector of current rotational eigenstrains.

Definition at line 145 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain(), and computeQpStrain().

_rot_eigenstrain_old

std::vector<const MaterialProperty<RealVectorValue> *> ComputeIncrementalBeamStrain::_rot_eigenstrain_old

protected

Vector of old rotational eigenstrains.

Definition at line 151 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain(), and computeQpStrain().

_rot_num

std::vector<unsigned int> ComputeIncrementalBeamStrain::_rot_num

protected

Variable numbers corresponding to the rotational variables.

Definition at line 58 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain(), and computeProperties().

_soln_disp_index_0

std::vector<unsigned int> ComputeIncrementalBeamStrain::_soln_disp_index_0

protected

Indices of solution vector corresponding to displacement DOFs at the node 0.

Definition at line 159 of file ComputeIncrementalBeamStrain.h.

Referenced by computeProperties().

_soln_disp_index_1

std::vector<unsigned int> ComputeIncrementalBeamStrain::_soln_disp_index_1

protected

Indices of solution vector corresponding to displacement DOFs at the node 1.

Definition at line 162 of file ComputeIncrementalBeamStrain.h.

Referenced by computeProperties().

_soln_rot_index_0

std::vector<unsigned int> ComputeIncrementalBeamStrain::_soln_rot_index_0

protected

Indices of solution vector corresponding to rotation DOFs at the node 0.

Definition at line 165 of file ComputeIncrementalBeamStrain.h.

Referenced by computeProperties().

_soln_rot_index_1

std::vector<unsigned int> ComputeIncrementalBeamStrain::_soln_rot_index_1

protected

Indices of solution vector corresponding to rotation DOFs at the node 1.

Definition at line 168 of file ComputeIncrementalBeamStrain.h.

Referenced by computeProperties().

_total_disp_strain

MaterialProperty<RealVectorValue>& ComputeIncrementalBeamStrain::_total_disp_strain

protected

Current total displacement strain integrated over the cross-section in global coordinate system.

Definition at line 91 of file ComputeIncrementalBeamStrain.h.

Referenced by computeQpStrain(), and initQpStatefulProperties().

_total_disp_strain_old

const MaterialProperty<RealVectorValue>& ComputeIncrementalBeamStrain::_total_disp_strain_old

protected

Old total displacement strain integrated over the cross-section in global coordinate system.

Definition at line 97 of file ComputeIncrementalBeamStrain.h.

Referenced by computeQpStrain().

_total_rot_strain

MaterialProperty<RealVectorValue>& ComputeIncrementalBeamStrain::_total_rot_strain

protected

Current total rotational strain integrated over the cross-section in global coordinate system.

Definition at line 94 of file ComputeIncrementalBeamStrain.h.

Referenced by computeQpStrain(), and initQpStatefulProperties().

_total_rot_strain_old

const MaterialProperty<RealVectorValue>& ComputeIncrementalBeamStrain::_total_rot_strain_old

protected

Old total rotational strain integrated over the cross-section in global coordinate system.

Definition at line 100 of file ComputeIncrementalBeamStrain.h.

_total_rotation

MaterialProperty<RankTwoTensor>& ComputeIncrementalBeamStrain::_total_rotation

protected

Rotational transformation from global coordinate system to beam local configuration at time t.

Definition at line 88 of file ComputeIncrementalBeamStrain.h.

Referenced by computeQpStrain(), ComputeFiniteBeamStrain::computeRotation(), computeRotation(), computeStiffnessMatrix(), and initQpStatefulProperties().

---

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

• tensor_mechanics/include/materials/ComputeIncrementalBeamStrain.h
• tensor_mechanics/src/materials/ComputeIncrementalBeamStrain.C