ComputeIncrementalBeamStrain Class Reference

#include <ComputeIncrementalBeamStrain.h>

Inheritance diagram for ComputeIncrementalBeamStrain:

Public Member Functions
	ComputeIncrementalBeamStrain (const InputParameters &parameters)

Protected Member Functions
virtual void	computeProperties () override
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
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...

Detailed Description

Definition at line 23 of file ComputeIncrementalBeamStrain.h.

Constructor & Destructor Documentation

ComputeIncrementalBeamStrain()

ComputeIncrementalBeamStrain::ComputeIncrementalBeamStrain

(

const InputParameters &

parameters

)

Definition at line 60 of file ComputeIncrementalBeamStrain.C.

  : Material(parameters),
    _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(isParamValid("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"))
{
  // 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 ( )

overrideprotectedvirtual

Definition at line 173 of file ComputeIncrementalBeamStrain.C.

{
  // fetch the two end nodes for current element
  std::vector<Node *> node;
  for (unsigned int i = 0; i < 2; ++i)
    node.push_back(_current_elem->get_node(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 220 of file ComputeIncrementalBeamStrain.C.

Referenced by computeProperties().

{
  Real Ix = _Ix[_qp];
  if (!isParamValid("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]);
  }
}

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 591 of file ComputeIncrementalBeamStrain.C.

Referenced by computeProperties().

{
  _total_rotation[0] = _original_local_config;
}

computeStiffnessMatrix()

void ComputeIncrementalBeamStrain::computeStiffnessMatrix ( )

protected

Computes the stiffness matrices.

Definition at line 322 of file ComputeIncrementalBeamStrain.C.

Referenced by computeProperties().

{
  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 (!isParamValid("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];
  }
}

initQpStatefulProperties()

void ComputeIncrementalBeamStrain::initQpStatefulProperties ( )

overrideprotectedvirtual

Definition at line 131 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;
}

Member Data Documentation

_area

const VariableValue& ComputeIncrementalBeamStrain::_area

protected

Coupled variable for the beam cross-sectional area.

Definition at line 54 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 126 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 57 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 60 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 144 of file ComputeIncrementalBeamStrain.h.

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

_disp1

RealVectorValue ComputeIncrementalBeamStrain::_disp1

protected

Definition at line 144 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 132 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 138 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 51 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain(), and computeProperties().

_eigenstrain_names

std::vector<MaterialPropertyName> ComputeIncrementalBeamStrain::_eigenstrain_names

protected

Vector of beam eigenstrain names.

Definition at line 129 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 120 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 123 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 161 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 69 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 63 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 66 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 102 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 108 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 105 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 111 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 114 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 117 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 99 of file ComputeIncrementalBeamStrain.h.

Referenced by 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 93 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 96 of file ComputeIncrementalBeamStrain.h.

Referenced by computeQpStrain().

_ndisp

unsigned int ComputeIncrementalBeamStrain::_ndisp

protected

Number of coupled displacement variables.

Definition at line 45 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 147 of file ComputeIncrementalBeamStrain.h.

Referenced by computeProperties().

_nrot

unsigned int ComputeIncrementalBeamStrain::_nrot

protected

Number of coupled rotational variables.

Definition at line 42 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain().

_original_length

MaterialProperty<Real>& ComputeIncrementalBeamStrain::_original_length

protected

Initial length of the beam.

Definition at line 75 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 72 of file ComputeIncrementalBeamStrain.h.

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

_rot0

RealVectorValue ComputeIncrementalBeamStrain::_rot0

protected

Definition at line 144 of file ComputeIncrementalBeamStrain.h.

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

_rot1

RealVectorValue ComputeIncrementalBeamStrain::_rot1

protected

Definition at line 144 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 135 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 141 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 48 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 149 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 152 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 155 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 158 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 81 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 87 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 84 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 90 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 78 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