20 "Residual due to stress related Rayleigh damping and HHT time integration terms ");
21 params.
addParam<MaterialPropertyName>(
"zeta",
23 "Name of material property or a constant real " 24 "number defining the zeta parameter for the " 26 params.
addParam<
Real>(
"alpha", 0,
"alpha parameter for HHT time integration");
27 params.
addParam<
bool>(
"static_initialization",
29 "Set to true to get the system to " 30 "equilibrium under gravity by running a " 31 "quasi-static analysis (by solving Ku = F) " 32 "in the first time step");
38 _stress_older(getMaterialPropertyOlderByName<
RankTwoTensor>(_base_name +
"stress")),
39 _stress_old(getMaterialPropertyOldByName<
RankTwoTensor>(_base_name +
"stress")),
40 _zeta(getMaterialProperty<
Real>(
"zeta")),
41 _alpha(getParam<
Real>(
"alpha")),
42 _static_initialization(getParam<bool>(
"static_initialization"))
71 residual +=
_stress[_qp].trace() / 3.0 *
110 for (
unsigned int i = 0; i <
_ndisp; ++i)
static InputParameters validParams()
virtual Real computeQpJacobian()
bool _volumetric_locking_correction
Flag for volumetric locking correction.
unsigned int _ndisp
Coupled displacement variables.
std::vector< std::vector< Real > > _avg_grad_test
Gradient of test function averaged over the element. Used in volumetric locking correction calculatio...
registerMooseObject("SolidMechanicsApp", DynamicStressDivergenceTensors)
const unsigned int _component
An integer corresponding to the direction this kernel acts in.
const MaterialProperty< Real > & _zeta
const bool _static_initialization
const MaterialProperty< RankTwoTensor > & _stress_old
virtual Real computeQpJacobian() override
std::vector< unsigned int > _disp_var
Displacement variables IDs.
virtual Real computeQpResidual()
DynamicStressDivergenceTensors(const InputParameters ¶meters)
virtual Real computeQpOffDiagJacobian(unsigned int jvar) override
virtual Real computeQpOffDiagJacobian(unsigned int jvar)
DynamicStressDivergenceTensors derives from StressDivergenceTensors and adds stress related Rayleigh ...
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
const MaterialProperty< RankTwoTensor > & _stress_older
{@ The old and older states of the stress tensor that the divergence operator operates on ...
static InputParameters validParams()
const MaterialProperty< RankTwoTensor > & _stress
The stress tensor that the divergence operator operates on.
StressDivergenceTensors mostly copies from StressDivergence.