ComputeFiniteBeamStrain Class Reference

#include <ComputeFiniteBeamStrain.h>

Inheritance diagram for ComputeFiniteBeamStrain:

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

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
void	computeRotation () override
	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...
virtual void	initQpStatefulProperties () override
void	computeQpStrain ()
	Computes the displacement and rotation strain increments. More...
void	computeStiffnessMatrix ()
	Computes the stiffness matrices. More...

Protected Attributes
const MaterialProperty< RankTwoTensor > &	_total_rotation_old
	Rotational transformation from the global to beam local coordinate system at time t. More...
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 25 of file ComputeFiniteBeamStrain.h.

Constructor & Destructor Documentation

ComputeFiniteBeamStrain()

ComputeFiniteBeamStrain::ComputeFiniteBeamStrain

(

const InputParameters &

parameters

)

Definition at line 32 of file ComputeFiniteBeamStrain.C.

  : ComputeIncrementalBeamStrain(parameters),
    _total_rotation_old(getMaterialPropertyOld<RankTwoTensor>("total_rotation"))
{
}

Member Function Documentation

computeProperties()

void ComputeIncrementalBeamStrain::computeProperties ( )

overridevirtualinherited

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

protectedinherited

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 ComputeIncrementalBeamStrain::computeProperties().

computeRotation()

void ComputeFiniteBeamStrain::computeRotation ( )

overrideprotectedvirtual

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 from ComputeIncrementalBeamStrain.

Definition at line 39 of file ComputeFiniteBeamStrain.C.

{
  // First - convert the differential beam displacement to the beam configuration at previous time
  // step
  const RealVectorValue delta_disp_local(_total_rotation_old[0] * (_disp1 - _disp0));
 
  // Second - calculate the incremental rotation matrix using Euler angles
  const Real intermediate_length_1 =
      std::sqrt(Utility::pow<2>(_original_length[0] + delta_disp_local(0)) +
                Utility::pow<2>(delta_disp_local(2)));
  const Real cos_alpha = (_original_length[0] + delta_disp_local(0)) / intermediate_length_1;
  const Real sin_alpha = std::sqrt(1.0 - Utility::pow<2>(cos_alpha));
 
  const Real intermediate_length_2 =
      std::sqrt(Utility::pow<2>(intermediate_length_1) + Utility::pow<2>(delta_disp_local(1)));
  const Real sin_beta = delta_disp_local(1) / intermediate_length_2;
  const Real cos_beta = std::sqrt(1.0 - Utility::pow<2>(sin_beta));
 
  const RealVectorValue rotation_d_1(cos_alpha * cos_beta, sin_beta, sin_alpha * cos_beta);
  const RealVectorValue rotation_d_2(-cos_alpha * sin_beta, cos_beta, -sin_alpha * sin_beta);
  const RealVectorValue rotation_d_3(-sin_alpha, 0.0, cos_alpha);
 
  const RankTwoTensor rotation_d(rotation_d_1, rotation_d_2, rotation_d_3);
 
  // Convert average rotational displacement to the beam configuration at previous time step
  const RealVectorValue avg_rot_local(_total_rotation_old[0] * (_rot0 + _rot1));
 
  const Real gamma_increment =
      0.5 * (rotation_d_1(0) * avg_rot_local(0) + rotation_d_1(1) * avg_rot_local(1) +
             rotation_d_1(2) * avg_rot_local(2));
 
  RankTwoTensor rotation_a;
  rotation_a(0, 0) = 1.0;
  rotation_a(1, 1) = std::cos(gamma_increment);
  rotation_a(1, 2) = std::sin(gamma_increment);
  rotation_a(2, 1) = -rotation_a(1, 2);
  rotation_a(2, 2) = rotation_a(1, 1);
 
  // Total rotation matrix from global configuration to beam local config at current time
  _total_rotation[0] = rotation_a * rotation_d * _total_rotation_old[0];
}

computeStiffnessMatrix()

void ComputeIncrementalBeamStrain::computeStiffnessMatrix ( )

protectedinherited

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 ComputeIncrementalBeamStrain::computeProperties().

initQpStatefulProperties()

void ComputeIncrementalBeamStrain::initQpStatefulProperties ( )

overrideprotectedvirtualinherited

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 ComputeFiniteBeamStrain::validParams ( )

static

Definition at line 23 of file ComputeFiniteBeamStrain.C.

{
  InputParameters params = ComputeIncrementalBeamStrain::validParams();
  params.addClassDescription("Compute a rotation increment for finite rotations of the beam and "
                             "computes the small/large strain increments in the current rotated "
                             "configuration of the beam.");
  return params;
}

Member Data Documentation

_area

const VariableValue& ComputeIncrementalBeamStrain::_area

protectedinherited

Coupled variable for the beam cross-sectional area.

Definition at line 64 of file ComputeIncrementalBeamStrain.h.

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

_avg_rot_local_t

RealVectorValue ComputeIncrementalBeamStrain::_avg_rot_local_t

protectedinherited

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

Definition at line 136 of file ComputeIncrementalBeamStrain.h.

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

_Ay

const VariableValue& ComputeIncrementalBeamStrain::_Ay

protectedinherited

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::ComputeIncrementalBeamStrain(), and ComputeIncrementalBeamStrain::computeQpStrain().

_Az

const VariableValue& ComputeIncrementalBeamStrain::_Az

protectedinherited

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::ComputeIncrementalBeamStrain(), and ComputeIncrementalBeamStrain::computeQpStrain().

_disp0

RealVectorValue ComputeIncrementalBeamStrain::_disp0

protectedinherited

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 ComputeIncrementalBeamStrain::computeProperties(), ComputeIncrementalBeamStrain::computeQpStrain(), and computeRotation().

_disp1

RealVectorValue ComputeIncrementalBeamStrain::_disp1

protectedinherited

Definition at line 154 of file ComputeIncrementalBeamStrain.h.

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

_disp_eigenstrain

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

protectedinherited

Vector of current displacement eigenstrains.

Definition at line 142 of file ComputeIncrementalBeamStrain.h.

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

_disp_eigenstrain_old

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

protectedinherited

Vector of old displacement eigenstrains.

Definition at line 148 of file ComputeIncrementalBeamStrain.h.

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

_disp_num

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

protectedinherited

Variable numbers corresponding to the displacement variables.

Definition at line 61 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain::ComputeIncrementalBeamStrain(), and ComputeIncrementalBeamStrain::computeProperties().

_effective_stiffness

MaterialProperty<Real>& ComputeIncrementalBeamStrain::_effective_stiffness

protectedinherited

Psuedo stiffness for critical time step computation.

Definition at line 174 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain::computeQpStrain().

_eigenstrain_names

std::vector<MaterialPropertyName> ComputeIncrementalBeamStrain::_eigenstrain_names

protectedinherited

Vector of beam eigenstrain names.

Definition at line 139 of file ComputeIncrementalBeamStrain.h.

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

_grad_disp_0_local_t

RealVectorValue ComputeIncrementalBeamStrain::_grad_disp_0_local_t

protectedinherited

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

Definition at line 130 of file ComputeIncrementalBeamStrain.h.

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

_grad_rot_0_local_t

RealVectorValue ComputeIncrementalBeamStrain::_grad_rot_0_local_t

protectedinherited

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

Definition at line 133 of file ComputeIncrementalBeamStrain.h.

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

_has_Ix

const bool ComputeIncrementalBeamStrain::_has_Ix

protectedinherited

Booleans for validity of params.

Definition at line 49 of file ComputeIncrementalBeamStrain.h.

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

_initial_rotation

MaterialProperty<RankTwoTensor>& ComputeIncrementalBeamStrain::_initial_rotation

protectedinherited

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

Definition at line 171 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain::computeProperties().

_Ix

const VariableValue& ComputeIncrementalBeamStrain::_Ix

protectedinherited

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 ComputeIncrementalBeamStrain::computeQpStrain(), and ComputeIncrementalBeamStrain::computeStiffnessMatrix().

_Iy

const VariableValue& ComputeIncrementalBeamStrain::_Iy

protectedinherited

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 ComputeIncrementalBeamStrain::computeQpStrain(), and ComputeIncrementalBeamStrain::computeStiffnessMatrix().

_Iz

const VariableValue& ComputeIncrementalBeamStrain::_Iz

protectedinherited

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 ComputeIncrementalBeamStrain::computeQpStrain(), and ComputeIncrementalBeamStrain::computeStiffnessMatrix().

_K11

MaterialProperty<RankTwoTensor>& ComputeIncrementalBeamStrain::_K11

protectedinherited

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

Definition at line 112 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain::computeStiffnessMatrix().

_K21

MaterialProperty<RankTwoTensor>& ComputeIncrementalBeamStrain::_K21

protectedinherited

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

Definition at line 118 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain::computeStiffnessMatrix().

_K21_cross

MaterialProperty<RankTwoTensor>& ComputeIncrementalBeamStrain::_K21_cross

protectedinherited

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

Definition at line 115 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain::computeStiffnessMatrix().

_K22

MaterialProperty<RankTwoTensor>& ComputeIncrementalBeamStrain::_K22

protectedinherited

Stiffness matrix between rotation DOFs of the same node.

Definition at line 121 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain::computeStiffnessMatrix().

_K22_cross

MaterialProperty<RankTwoTensor>& ComputeIncrementalBeamStrain::_K22_cross

protectedinherited

Stiffness matrix between rotation DOFs of different nodes.

Definition at line 124 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain::computeStiffnessMatrix().

_large_strain

const bool ComputeIncrementalBeamStrain::_large_strain

protectedinherited

Boolean flag to turn on large strain calculation.

Definition at line 127 of file ComputeIncrementalBeamStrain.h.

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

_material_stiffness

const MaterialProperty<RealVectorValue>& ComputeIncrementalBeamStrain::_material_stiffness

protectedinherited

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

Definition at line 109 of file ComputeIncrementalBeamStrain.h.

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

_mech_disp_strain_increment

MaterialProperty<RealVectorValue>& ComputeIncrementalBeamStrain::_mech_disp_strain_increment

protectedinherited

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

Definition at line 103 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain::computeQpStrain().

_mech_rot_strain_increment

MaterialProperty<RealVectorValue>& ComputeIncrementalBeamStrain::_mech_rot_strain_increment

protectedinherited

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

Definition at line 106 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain::computeQpStrain().

_ndisp

unsigned int ComputeIncrementalBeamStrain::_ndisp

protectedinherited

Number of coupled displacement variables.

Definition at line 55 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain::ComputeIncrementalBeamStrain(), and ComputeIncrementalBeamStrain::computeProperties().

_nonlinear_sys

NonlinearSystemBase& ComputeIncrementalBeamStrain::_nonlinear_sys

protectedinherited

Reference to the nonlinear system object.

Definition at line 157 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain::computeProperties().

_nrot

unsigned int ComputeIncrementalBeamStrain::_nrot

protectedinherited

Number of coupled rotational variables.

Definition at line 52 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain::ComputeIncrementalBeamStrain().

_original_length

MaterialProperty<Real>& ComputeIncrementalBeamStrain::_original_length

protectedinherited

Initial length of the beam.

Definition at line 85 of file ComputeIncrementalBeamStrain.h.

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

_original_local_config

RankTwoTensor ComputeIncrementalBeamStrain::_original_local_config

protectedinherited

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

Definition at line 82 of file ComputeIncrementalBeamStrain.h.

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

_prefactor_function

const Function* const ComputeIncrementalBeamStrain::_prefactor_function

protectedinherited

Prefactor function to multiply the elasticity tensor with.

Definition at line 177 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain::computeQpStrain().

_rot0

RealVectorValue ComputeIncrementalBeamStrain::_rot0

protectedinherited

Definition at line 154 of file ComputeIncrementalBeamStrain.h.

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

_rot1

RealVectorValue ComputeIncrementalBeamStrain::_rot1

protectedinherited

Definition at line 154 of file ComputeIncrementalBeamStrain.h.

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

_rot_eigenstrain

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

protectedinherited

Vector of current rotational eigenstrains.

Definition at line 145 of file ComputeIncrementalBeamStrain.h.

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

_rot_eigenstrain_old

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

protectedinherited

Vector of old rotational eigenstrains.

Definition at line 151 of file ComputeIncrementalBeamStrain.h.

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

_rot_num

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

protectedinherited

Variable numbers corresponding to the rotational variables.

Definition at line 58 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain::ComputeIncrementalBeamStrain(), and ComputeIncrementalBeamStrain::computeProperties().

_soln_disp_index_0

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

protectedinherited

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

Definition at line 159 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain::computeProperties().

_soln_disp_index_1

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

protectedinherited

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

Definition at line 162 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain::computeProperties().

_soln_rot_index_0

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

protectedinherited

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

Definition at line 165 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain::computeProperties().

_soln_rot_index_1

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

protectedinherited

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

Definition at line 168 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain::computeProperties().

_total_disp_strain

MaterialProperty<RealVectorValue>& ComputeIncrementalBeamStrain::_total_disp_strain

protectedinherited

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

Definition at line 91 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain::computeQpStrain(), and ComputeIncrementalBeamStrain::initQpStatefulProperties().

_total_disp_strain_old

const MaterialProperty<RealVectorValue>& ComputeIncrementalBeamStrain::_total_disp_strain_old

protectedinherited

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

Definition at line 97 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain::computeQpStrain().

_total_rot_strain

MaterialProperty<RealVectorValue>& ComputeIncrementalBeamStrain::_total_rot_strain

protectedinherited

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

Definition at line 94 of file ComputeIncrementalBeamStrain.h.

Referenced by ComputeIncrementalBeamStrain::computeQpStrain(), and ComputeIncrementalBeamStrain::initQpStatefulProperties().

_total_rot_strain_old

const MaterialProperty<RealVectorValue>& ComputeIncrementalBeamStrain::_total_rot_strain_old

protectedinherited

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

protectedinherited

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

Definition at line 88 of file ComputeIncrementalBeamStrain.h.

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

_total_rotation_old

const MaterialProperty<RankTwoTensor>& ComputeFiniteBeamStrain::_total_rotation_old

protected

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

Definition at line 37 of file ComputeFiniteBeamStrain.h.

Referenced by computeRotation().

---

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

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