ComputeIncrementalBeamStrain.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#pragma once

#include "Material.h"
#include "RankTwoTensor.h"

// Forward Declarations
class Function;

class ComputeIncrementalBeamStrain : public Material
{
public:
  static InputParameters validParams();

  ComputeIncrementalBeamStrain(const InputParameters & parameters);

  virtual void computeProperties() override;

protected:
  virtual void initQpStatefulProperties() override;

  void computeQpStrain();

  void computeStiffnessMatrix();

  virtual void computeRotation();

  const bool _has_Ix;

  unsigned int _nrot;

  unsigned int _ndisp;

  std::vector<unsigned int> _rot_num;

  std::vector<unsigned int> _disp_num;

  const VariableValue & _area;

  const VariableValue & _Ay;

  const VariableValue & _Az;

  const VariableValue & _Iy;

  const VariableValue & _Iz;

  const VariableValue & _Ix;

  RankTwoTensor & _original_local_config;

  MaterialProperty<Real> & _original_length;

  MaterialProperty<RankTwoTensor> & _total_rotation;

  MaterialProperty<RealVectorValue> & _total_disp_strain;

  MaterialProperty<RealVectorValue> & _total_rot_strain;

  const MaterialProperty<RealVectorValue> & _total_disp_strain_old;

  const MaterialProperty<RealVectorValue> & _total_rot_strain_old;

  MaterialProperty<RealVectorValue> & _mech_disp_strain_increment;

  MaterialProperty<RealVectorValue> & _mech_rot_strain_increment;

  const MaterialProperty<RealVectorValue> & _material_stiffness;

  MaterialProperty<RankTwoTensor> & _K11;

  MaterialProperty<RankTwoTensor> & _K21_cross;

  MaterialProperty<RankTwoTensor> & _K21;

  MaterialProperty<RankTwoTensor> & _K22;

  MaterialProperty<RankTwoTensor> & _K22_cross;

  const bool _large_strain;

  RealVectorValue _grad_disp_0_local_t;

  RealVectorValue _grad_rot_0_local_t;

  RealVectorValue _avg_rot_local_t;

  std::vector<MaterialPropertyName> _eigenstrain_names;

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

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

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

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

  RealVectorValue _disp0, _disp1, _rot0, _rot1;

  NonlinearSystemBase & _nonlinear_sys;
  std::vector<unsigned int> _soln_disp_index_0;

  std::vector<unsigned int> _soln_disp_index_1;

  std::vector<unsigned int> _soln_rot_index_0;

  std::vector<unsigned int> _soln_rot_index_1;

  MaterialProperty<RankTwoTensor> & _initial_rotation;

  MaterialProperty<Real> & _effective_stiffness;

  const Function * const _prefactor_function;
};