#include <InertialForce.h>
Definition at line 22 of file InertialForce.h.
◆ InertialForce()
InertialForce::InertialForce |
( |
const InputParameters & |
parameters | ) |
|
Definition at line 44 of file InertialForce.C.
45 : TimeKernel(parameters),
46 _density(getMaterialProperty<Real>(
"density")),
47 _beta(isParamValid(
"beta") ? getParam<Real>(
"beta") : 0.1),
48 _gamma(isParamValid(
"gamma") ? getParam<Real>(
"gamma") : 0.1),
49 _eta(getMaterialProperty<Real>(
"eta")),
50 _alpha(getParam<Real>(
"alpha"))
52 if (isParamValid(
"beta") && isParamValid(
"gamma") && isParamValid(
"velocity") &&
53 isParamValid(
"acceleration"))
55 _vel_old = &coupledValueOld(
"velocity");
59 else if (!isParamValid(
"beta") && !isParamValid(
"gamma") && !isParamValid(
"velocity") &&
60 !isParamValid(
"acceleration"))
69 mooseError(
"InertialForce: Either all or none of `beta`, `gamma`, `velocity`and `acceleration` " 70 "should be provided as input.");
const VariableValue * _u_dot
const VariableValue * _du_dotdot_du
const VariableValue * _accel_old
const VariableValue * _u_old
const MaterialProperty< Real > & _density
const VariableValue * _du_dot_du
const VariableValue * _u_dot_old
const VariableValue * _u_dotdot
const MaterialProperty< Real > & _eta
const VariableValue * _vel_old
◆ computeQpJacobian()
Real InertialForce::computeQpJacobian |
( |
| ) |
|
|
protectedvirtual |
Definition at line 94 of file InertialForce.C.
98 else if (isParamValid(
"beta"))
99 return _test[_i][_qp] *
_density[_qp] / (
_beta * _dt * _dt) * _phi[_j][_qp] +
103 return _test[_i][_qp] *
_density[_qp] * (*_du_dotdot_du)[_qp] * _phi[_j][_qp] +
const MaterialProperty< Real > & _density
const VariableValue * _du_dot_du
const MaterialProperty< Real > & _eta
◆ computeQpResidual()
Real InertialForce::computeQpResidual |
( |
| ) |
|
|
protectedvirtual |
Definition at line 74 of file InertialForce.C.
78 else if (isParamValid(
"beta"))
80 Real accel = 1. /
_beta *
81 (((_u[_qp] - (*_u_old)[_qp]) / (_dt * _dt)) - (*
_vel_old)[_qp] / _dt -
82 (*_accel_old)[_qp] * (0.5 -
_beta));
83 Real vel = (*_vel_old)[_qp] + (_dt * (1 -
_gamma)) * (*_accel_old)[_qp] +
_gamma * _dt * accel;
84 return _test[_i][_qp] *
_density[_qp] *
88 return _test[_i][_qp] *
_density[_qp] *
89 ((*_u_dotdot)[_qp] + (*_u_dot)[_qp] *
_eta[_qp] * (1.0 +
_alpha) -
const MaterialProperty< Real > & _density
const VariableValue * _u_dot_old
const MaterialProperty< Real > & _eta
const VariableValue * _vel_old
◆ _accel_old
const VariableValue* InertialForce::_accel_old |
|
private |
◆ _alpha
const Real InertialForce::_alpha |
|
private |
◆ _beta
const Real InertialForce::_beta |
|
private |
◆ _density
const MaterialProperty<Real>& InertialForce::_density |
|
private |
◆ _du_dot_du
const VariableValue* InertialForce::_du_dot_du |
|
private |
◆ _du_dotdot_du
const VariableValue* InertialForce::_du_dotdot_du |
|
private |
◆ _eta
const MaterialProperty<Real>& InertialForce::_eta |
|
private |
◆ _gamma
const Real InertialForce::_gamma |
|
private |
◆ _u_dot
const VariableValue* InertialForce::_u_dot |
|
private |
◆ _u_dot_old
const VariableValue* InertialForce::_u_dot_old |
|
private |
◆ _u_dotdot
const VariableValue* InertialForce::_u_dotdot |
|
private |
◆ _u_old
const VariableValue* InertialForce::_u_old |
|
private |
◆ _vel_old
const VariableValue* InertialForce::_vel_old |
|
private |
The documentation for this class was generated from the following files: