www.mooseframework.org
Public Member Functions | Static Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
StressDivergenceRSphericalTensors Class Reference

StressDivergenceRSphericalTensors is a modification of StressDivergenceTensors for 1D spherically symmetric problems. More...

#include <StressDivergenceRSphericalTensors.h>

Inheritance diagram for StressDivergenceRSphericalTensors:
[legend]

Public Member Functions

 StressDivergenceRSphericalTensors (const InputParameters &parameters)
 
virtual void computeJacobian () override
 
virtual void computeOffDiagJacobian (MooseVariableFEBase &jvar) override
 

Static Public Member Functions

static InputParameters validParams ()
 

Protected Member Functions

virtual void initialSetup () override
 
virtual Real computeQpResidual () override
 
virtual Real computeQpJacobian () override
 
virtual Real computeQpOffDiagJacobian (unsigned int jvar) override
 
Real calculateJacobian (unsigned int ivar, unsigned int jvar)
 
virtual void computeResidual () override
 
virtual void computeFiniteDeformJacobian ()
 
virtual void computeAverageGradientTest ()
 
virtual void computeAverageGradientPhi ()
 

Protected Attributes

const std::string _base_name
 
bool _use_finite_deform_jacobian
 
const MaterialProperty< RankTwoTensor > & _stress
 
const MaterialProperty< RankFourTensor > & _Jacobian_mult
 
std::vector< RankFourTensor_finite_deform_Jacobian_mult
 
const MaterialProperty< RankTwoTensor > * _deformation_gradient
 
const MaterialProperty< RankTwoTensor > * _deformation_gradient_old
 
const MaterialProperty< RankTwoTensor > * _rotation_increment
 
const unsigned int _component
 
unsigned int _ndisp
 Coupled displacement variables. More...
 
std::vector< unsigned int > _disp_var
 
const bool _temp_coupled
 
const unsigned int _temp_var
 
const MaterialProperty< RankTwoTensor > *const _deigenstrain_dT
 d(strain)/d(temperature), if computed by ComputeThermalExpansionEigenstrain More...
 
const bool _out_of_plane_strain_coupled
 
const VariableValue * _out_of_plane_strain
 
const unsigned int _out_of_plane_strain_var
 
const unsigned int _out_of_plane_direction
 
const bool _use_displaced_mesh
 Whether this object is acting on the displaced mesh. More...
 
std::vector< std::vector< Real > > _avg_grad_test
 Gradient of test function averaged over the element. Used in volumetric locking correction calculation. More...
 
std::vector< std::vector< Real > > _avg_grad_phi
 Gradient of phi function averaged over the element. Used in volumetric locking correction calculation. More...
 
bool _volumetric_locking_correction
 Flag for volumetric locking correction. More...
 
Assembly & _assembly_undisplaced
 undisplaced problem More...
 
MooseVariable & _var_undisplaced
 Reference to this Kernel's undisplaced MooseVariable object. More...
 
const VariablePhiGradient & _grad_phi_undisplaced
 Shape and test functions on the undisplaced mesh. More...
 
const VariableTestGradient & _grad_test_undisplaced
 

Detailed Description

StressDivergenceRSphericalTensors is a modification of StressDivergenceTensors for 1D spherically symmetric problems.

The main modifications from the original StressDivergenceTensors code are requirements from the dependence of stress in the polar and azimuthal stresses on displacement and position in the radial direction. This kernel is for symmetrical loading only. If solving an anisotropic material problem, recall that the orientation of the basis vectors (\hat{e}_r) change with position, so the components of the elasticity tensor are functions of position. Reference: Bower, A.F. Applied Mechanics of Solids (2012). Chapter 4. Available online at solidmechanics.org Within this kernel, '_disp_x' refers to displacement in the radial direction. The COORD_TYPE in the Problem block must be set to RSpherical.

Definition at line 34 of file StressDivergenceRSphericalTensors.h.

Constructor & Destructor Documentation

◆ StressDivergenceRSphericalTensors()

StressDivergenceRSphericalTensors::StressDivergenceRSphericalTensors ( const InputParameters &  parameters)

Definition at line 30 of file StressDivergenceRSphericalTensors.C.

32  : StressDivergenceTensors(parameters)
33 {
34  if (_component != 0)
35  mooseError("Invalid component for this 1D RSpherical problem.");
36 }

Member Function Documentation

◆ calculateJacobian()

Real StressDivergenceRSphericalTensors::calculateJacobian ( unsigned int  ivar,
unsigned int  jvar 
)
protected

Definition at line 74 of file StressDivergenceRSphericalTensors.C.

75 {
76  RealGradient test_r, phi_r;
77 
78  mooseAssert(ivar == 0 && jvar == 0,
79  "Invalid component in Jacobian Calculation"); // Only nonzero case for a 1D simulation
80 
81  if (ivar == 0) // Case grad_test for r, requires contributions from stress_{rr}, stress_{\theta
82  // \theta}, and stress_{\phi \phi}
83  {
84  test_r(0) = _grad_test[_i][_qp](0);
85  test_r(1) = _test[_i][_qp] / _q_point[_qp](0);
86  test_r(2) = _test[_i][_qp] / _q_point[_qp](0);
87  }
88 
89  if (jvar == 0)
90  {
91  phi_r(0) = _grad_phi[_j][_qp](0);
92  phi_r(1) = _phi[_j][_qp] / _q_point[_qp](0);
93  phi_r(2) = _phi[_j][_qp] / _q_point[_qp](0);
94  }
95 
96  return ElasticityTensorTools::elasticJacobian(_Jacobian_mult[_qp], ivar, jvar, test_r, phi_r);
97 }

Referenced by computeQpJacobian(), and computeQpOffDiagJacobian().

◆ computeAverageGradientPhi()

void StressDivergenceTensors::computeAverageGradientPhi ( )
protectedvirtualinherited

Reimplemented in StressDivergenceRZTensors.

Definition at line 423 of file StressDivergenceTensors.C.

424 {
425  // Calculate volume average derivatives for phi
426  _avg_grad_phi.resize(_phi.size());
427  for (_i = 0; _i < _phi.size(); ++_i)
428  {
429  _avg_grad_phi[_i].resize(3);
430  for (unsigned int component = 0; component < _mesh.dimension(); ++component)
431  {
432  _avg_grad_phi[_i][component] = 0.0;
433  for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
434  _avg_grad_phi[_i][component] += _grad_phi[_i][_qp](component) * _JxW[_qp] * _coord[_qp];
435 
436  _avg_grad_phi[_i][component] /= _current_elem_volume;
437  }
438  }
439 }

Referenced by StressDivergenceTensors::computeJacobian(), and StressDivergenceTensors::computeOffDiagJacobian().

◆ computeAverageGradientTest()

void StressDivergenceTensors::computeAverageGradientTest ( )
protectedvirtualinherited

Reimplemented in StressDivergenceRZTensors.

Definition at line 407 of file StressDivergenceTensors.C.

408 {
409  // Calculate volume averaged value of shape function derivative
410  _avg_grad_test.resize(_test.size());
411  for (_i = 0; _i < _test.size(); ++_i)
412  {
413  _avg_grad_test[_i].resize(3);
414  _avg_grad_test[_i][_component] = 0.0;
415  for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
416  _avg_grad_test[_i][_component] += _grad_test[_i][_qp](_component) * _JxW[_qp] * _coord[_qp];
417 
418  _avg_grad_test[_i][_component] /= _current_elem_volume;
419  }
420 }

Referenced by StressDivergenceTensors::computeJacobian(), StressDivergenceTensors::computeOffDiagJacobian(), and StressDivergenceTensors::computeResidual().

◆ computeFiniteDeformJacobian()

void StressDivergenceTensors::computeFiniteDeformJacobian ( )
protectedvirtualinherited

Definition at line 350 of file StressDivergenceTensors.C.

351 {
352  const RankTwoTensor I(RankTwoTensor::initIdentity);
353  const RankFourTensor II_ijkl = I.mixedProductIkJl(I);
354 
355  // Bring back to unrotated config
356  const RankTwoTensor unrotated_stress =
357  (*_rotation_increment)[_qp].transpose() * _stress[_qp] * (*_rotation_increment)[_qp];
358 
359  // Incremental deformation gradient Fhat
360  const RankTwoTensor Fhat =
361  (*_deformation_gradient)[_qp] * (*_deformation_gradient_old)[_qp].inverse();
362  const RankTwoTensor Fhatinv = Fhat.inverse();
363 
364  const RankTwoTensor rot_times_stress = (*_rotation_increment)[_qp] * unrotated_stress;
365  const RankFourTensor dstress_drot =
366  I.mixedProductIkJl(rot_times_stress) + I.mixedProductJkIl(rot_times_stress);
367  const RankFourTensor rot_rank_four =
368  (*_rotation_increment)[_qp].mixedProductIkJl((*_rotation_increment)[_qp]);
369  const RankFourTensor drot_dUhatinv = Fhat.mixedProductIkJl(I);
370 
371  const RankTwoTensor A = I - Fhatinv;
372 
373  // Ctilde = Chat^-1 - I
374  const RankTwoTensor Ctilde = A * A.transpose() - A - A.transpose();
375  const RankFourTensor dCtilde_dFhatinv =
376  -I.mixedProductIkJl(A) - I.mixedProductJkIl(A) + II_ijkl + I.mixedProductJkIl(I);
377 
378  // Second order approximation of Uhat - consistent with strain increment definition
379  // const RankTwoTensor Uhat = I - 0.5 * Ctilde - 3.0/8.0 * Ctilde * Ctilde;
380 
381  RankFourTensor dUhatinv_dCtilde =
382  0.5 * II_ijkl - 1.0 / 8.0 * (I.mixedProductIkJl(Ctilde) + Ctilde.mixedProductIkJl(I));
383  RankFourTensor drot_dFhatinv = drot_dUhatinv * dUhatinv_dCtilde * dCtilde_dFhatinv;
384 
385  drot_dFhatinv -= Fhat.mixedProductIkJl((*_rotation_increment)[_qp].transpose());
386  _finite_deform_Jacobian_mult[_qp] = dstress_drot * drot_dFhatinv;
387 
388  const RankFourTensor dstrain_increment_dCtilde =
389  -0.5 * II_ijkl + 0.25 * (I.mixedProductIkJl(Ctilde) + Ctilde.mixedProductIkJl(I));
391  rot_rank_four * _Jacobian_mult[_qp] * dstrain_increment_dCtilde * dCtilde_dFhatinv;
392  _finite_deform_Jacobian_mult[_qp] += Fhat.mixedProductJkIl(_stress[_qp]);
393 
394  const RankFourTensor dFhat_dFhatinv = -Fhat.mixedProductIkJl(Fhat.transpose());
395  const RankTwoTensor dJ_dFhatinv = dFhat_dFhatinv.innerProductTranspose(Fhat.ddet());
396 
397  // Component from Jacobian derivative
398  _finite_deform_Jacobian_mult[_qp] += _stress[_qp].outerProduct(dJ_dFhatinv);
399 
400  // Derivative of Fhatinv w.r.t. undisplaced coordinates
401  const RankTwoTensor Finv = (*_deformation_gradient)[_qp].inverse();
402  const RankFourTensor dFhatinv_dGradu = -Fhatinv.mixedProductIkJl(Finv.transpose());
403  _finite_deform_Jacobian_mult[_qp] = _finite_deform_Jacobian_mult[_qp] * dFhatinv_dGradu;
404 }

Referenced by StressDivergenceTensors::computeJacobian(), and StressDivergenceTensors::computeOffDiagJacobian().

◆ computeJacobian()

void StressDivergenceTensors::computeJacobian ( )
overridevirtualinherited

Reimplemented from ALEKernel.

Definition at line 168 of file StressDivergenceTensors.C.

169 {
171  {
174  }
175 
177  {
178  _finite_deform_Jacobian_mult.resize(_qrule->n_points());
179 
180  for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
182 
184  }
185  else
186  Kernel::computeJacobian();
187 }

◆ computeOffDiagJacobian()

void StressDivergenceTensors::computeOffDiagJacobian ( MooseVariableFEBase &  jvar)
overridevirtualinherited

Reimplemented from ALEKernel.

Definition at line 190 of file StressDivergenceTensors.C.

191 {
193  {
196  }
197 
199  {
200  _finite_deform_Jacobian_mult.resize(_qrule->n_points());
201 
202  for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
204 
206  }
207  else
208  Kernel::computeOffDiagJacobian(jvar);
209 }

◆ computeQpJacobian()

Real StressDivergenceRSphericalTensors::computeQpJacobian ( )
overrideprotectedvirtual

Reimplemented from StressDivergenceTensors.

Definition at line 55 of file StressDivergenceRSphericalTensors.C.

56 {
58 }

◆ computeQpOffDiagJacobian()

Real StressDivergenceRSphericalTensors::computeQpOffDiagJacobian ( unsigned int  jvar)
overrideprotectedvirtual

Reimplemented from StressDivergenceTensors.

Definition at line 61 of file StressDivergenceRSphericalTensors.C.

62 {
63  for (unsigned int i = 0; i < _ndisp; ++i)
64  if (jvar == _disp_var[i])
65  return calculateJacobian(_component, i);
66 
67  if (_temp_coupled && jvar == _temp_var)
68  return 0.0;
69 
70  return 0.0;
71 }

◆ computeQpResidual()

Real StressDivergenceRSphericalTensors::computeQpResidual ( )
overrideprotectedvirtual

Reimplemented from StressDivergenceTensors.

Definition at line 47 of file StressDivergenceRSphericalTensors.C.

48 {
49  return _grad_test[_i][_qp](0) * _stress[_qp](0, 0) + // stress_{rr} part 1
50  +(_test[_i][_qp] / _q_point[_qp](0)) * _stress[_qp](1, 1) + // stress_{\theta \theta}
51  +(_test[_i][_qp] / _q_point[_qp](0)) * _stress[_qp](2, 2); // stress_{\phi \phi}
52 }

◆ computeResidual()

void StressDivergenceTensors::computeResidual ( )
overrideprotectedvirtualinherited

Definition at line 126 of file StressDivergenceTensors.C.

127 {
128  prepareVectorTag(_assembly, _var.number());
129 
132 
133  precalculateResidual();
134  for (_i = 0; _i < _test.size(); ++_i)
135  for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
136  _local_re(_i) += _JxW[_qp] * _coord[_qp] * computeQpResidual();
137 
138  accumulateTaggedLocalResidual();
139 
140  if (_has_save_in)
141  {
142  Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
143  for (const auto & var : _save_in)
144  var->sys().solution().add_vector(_local_re, var->dofIndices());
145  }
146 }

◆ initialSetup()

void StressDivergenceRSphericalTensors::initialSetup ( )
overrideprotectedvirtual

Reimplemented from StressDivergenceTensors.

Definition at line 39 of file StressDivergenceRSphericalTensors.C.

40 {
41  if (getBlockCoordSystem() != Moose::COORD_RSPHERICAL)
42  mooseError("The coordinate system in the Problem block must be set to RSPHERICAL for 1D "
43  "spherically symmetric geometries.");
44 }

◆ validParams()

InputParameters StressDivergenceRSphericalTensors::validParams ( )
static

Definition at line 20 of file StressDivergenceRSphericalTensors.C.

21 {
22  InputParameters params = StressDivergenceTensors::validParams();
23  params.addClassDescription(
24  "Calculate stress divergence for a spherically symmetric 1D problem in polar coordinates.");
25  params.set<unsigned int>("component") = 0;
26  params.set<bool>("use_displaced_mesh") = true;
27  return params;
28 }

Member Data Documentation

◆ _assembly_undisplaced

Assembly& ALEKernel::_assembly_undisplaced
protectedinherited

undisplaced problem

Definition at line 34 of file ALEKernel.h.

◆ _avg_grad_phi

std::vector<std::vector<Real> > StressDivergenceTensors::_avg_grad_phi
protectedinherited

◆ _avg_grad_test

std::vector<std::vector<Real> > StressDivergenceTensors::_avg_grad_test
protectedinherited

◆ _base_name

const std::string StressDivergenceTensors::_base_name
protectedinherited

◆ _component

const unsigned int StressDivergenceTensors::_component
protectedinherited

◆ _deformation_gradient

const MaterialProperty<RankTwoTensor>* StressDivergenceTensors::_deformation_gradient
protectedinherited

◆ _deformation_gradient_old

const MaterialProperty<RankTwoTensor>* StressDivergenceTensors::_deformation_gradient_old
protectedinherited

◆ _deigenstrain_dT

const MaterialProperty<RankTwoTensor>* const StressDivergenceTensors::_deigenstrain_dT
protectedinherited

d(strain)/d(temperature), if computed by ComputeThermalExpansionEigenstrain

Definition at line 72 of file StressDivergenceTensors.h.

◆ _disp_var

std::vector<unsigned int> StressDivergenceTensors::_disp_var
protectedinherited

◆ _finite_deform_Jacobian_mult

std::vector<RankFourTensor> StressDivergenceTensors::_finite_deform_Jacobian_mult
protectedinherited

◆ _grad_phi_undisplaced

const VariablePhiGradient& ALEKernel::_grad_phi_undisplaced
protectedinherited

Shape and test functions on the undisplaced mesh.

Definition at line 40 of file ALEKernel.h.

Referenced by StressDivergenceTensors::computeQpJacobian(), and StressDivergenceTensors::computeQpOffDiagJacobian().

◆ _grad_test_undisplaced

const VariableTestGradient& ALEKernel::_grad_test_undisplaced
protectedinherited

Definition at line 41 of file ALEKernel.h.

◆ _Jacobian_mult

const MaterialProperty<RankFourTensor>& StressDivergenceTensors::_Jacobian_mult
protectedinherited

◆ _ndisp

unsigned int StressDivergenceTensors::_ndisp
protectedinherited

◆ _out_of_plane_direction

const unsigned int StressDivergenceTensors::_out_of_plane_direction
protectedinherited

◆ _out_of_plane_strain

const VariableValue* StressDivergenceTensors::_out_of_plane_strain
protectedinherited

◆ _out_of_plane_strain_coupled

const bool StressDivergenceTensors::_out_of_plane_strain_coupled
protectedinherited

◆ _out_of_plane_strain_var

const unsigned int StressDivergenceTensors::_out_of_plane_strain_var
protectedinherited

◆ _rotation_increment

const MaterialProperty<RankTwoTensor>* StressDivergenceTensors::_rotation_increment
protectedinherited

◆ _stress

const MaterialProperty<RankTwoTensor>& StressDivergenceTensors::_stress
protectedinherited

◆ _temp_coupled

const bool StressDivergenceTensors::_temp_coupled
protectedinherited

◆ _temp_var

const unsigned int StressDivergenceTensors::_temp_var
protectedinherited

◆ _use_displaced_mesh

const bool StressDivergenceTensors::_use_displaced_mesh
protectedinherited

Whether this object is acting on the displaced mesh.

Definition at line 80 of file StressDivergenceTensors.h.

Referenced by StressDivergenceTensors::computeQpJacobian(), and StressDivergenceTensors::computeQpResidual().

◆ _use_finite_deform_jacobian

bool StressDivergenceTensors::_use_finite_deform_jacobian
protectedinherited

◆ _var_undisplaced

MooseVariable& ALEKernel::_var_undisplaced
protectedinherited

Reference to this Kernel's undisplaced MooseVariable object.

Definition at line 37 of file ALEKernel.h.

◆ _volumetric_locking_correction

bool StressDivergenceTensors::_volumetric_locking_correction
protectedinherited

The documentation for this class was generated from the following files:
ElasticityTensorTools::elasticJacobian
Real elasticJacobian(const RankFourTensor &r4t, unsigned int i, unsigned int k, const RealGradient &grad_test, const RealGradient &grad_phi)
This is used for the standard kernel stress_ij*d(test)/dx_j, when varied wrt u_k Jacobian entry: d(st...
Definition: ElasticityTensorTools.C:21
ALEKernel::computeOffDiagJacobian
virtual void computeOffDiagJacobian(MooseVariableFEBase &jvar) override
Definition: ALEKernel.C:43
StressDivergenceTensors::StressDivergenceTensors
StressDivergenceTensors(const InputParameters &parameters)
Definition: StressDivergenceTensors.C:62
StressDivergenceTensors::_component
const unsigned int _component
Definition: StressDivergenceTensors.h:62
StressDivergenceTensors::computeAverageGradientPhi
virtual void computeAverageGradientPhi()
Definition: StressDivergenceTensors.C:423
StressDivergenceTensors::_disp_var
std::vector< unsigned int > _disp_var
Definition: StressDivergenceTensors.h:66
libMesh::RealGradient
VectorValue< Real > RealGradient
Definition: GrainForceAndTorqueInterface.h:17
StressDivergenceTensors::_rotation_increment
const MaterialProperty< RankTwoTensor > * _rotation_increment
Definition: StressDivergenceTensors.h:59
StressDivergenceTensors::_temp_coupled
const bool _temp_coupled
Definition: StressDivergenceTensors.h:68
StressDivergenceTensors::_avg_grad_phi
std::vector< std::vector< Real > > _avg_grad_phi
Gradient of phi function averaged over the element. Used in volumetric locking correction calculation...
Definition: StressDivergenceTensors.h:86
ALEKernel::computeJacobian
virtual void computeJacobian() override
Definition: ALEKernel.C:36
StressDivergenceTensors::computeFiniteDeformJacobian
virtual void computeFiniteDeformJacobian()
Definition: StressDivergenceTensors.C:350
StressDivergenceTensors::_volumetric_locking_correction
bool _volumetric_locking_correction
Flag for volumetric locking correction.
Definition: StressDivergenceTensors.h:89
StressDivergenceTensors::computeQpResidual
virtual Real computeQpResidual() override
Definition: StressDivergenceTensors.C:149
StressDivergenceTensors::validParams
static InputParameters validParams()
Definition: StressDivergenceTensors.C:26
MaterialTensorCalculatorTools::component
Real component(const SymmTensor &symm_tensor, unsigned int index)
Definition: MaterialTensorCalculatorTools.C:16
StressDivergenceTensors::_Jacobian_mult
const MaterialProperty< RankFourTensor > & _Jacobian_mult
Definition: StressDivergenceTensors.h:54
StressDivergenceTensors::_finite_deform_Jacobian_mult
std::vector< RankFourTensor > _finite_deform_Jacobian_mult
Definition: StressDivergenceTensors.h:56
StressDivergenceTensors::_stress
const MaterialProperty< RankTwoTensor > & _stress
Definition: StressDivergenceTensors.h:53
StressDivergenceRSphericalTensors::calculateJacobian
Real calculateJacobian(unsigned int ivar, unsigned int jvar)
Definition: StressDivergenceRSphericalTensors.C:74
RankFourTensorTempl
Definition: ACGrGrElasticDrivingForce.h:20
StressDivergenceTensors::_avg_grad_test
std::vector< std::vector< Real > > _avg_grad_test
Gradient of test function averaged over the element. Used in volumetric locking correction calculatio...
Definition: StressDivergenceTensors.h:83
StressDivergenceTensors::computeAverageGradientTest
virtual void computeAverageGradientTest()
Definition: StressDivergenceTensors.C:407
StressDivergenceTensors::_ndisp
unsigned int _ndisp
Coupled displacement variables.
Definition: StressDivergenceTensors.h:65
StressDivergenceTensors::_use_finite_deform_jacobian
bool _use_finite_deform_jacobian
Definition: StressDivergenceTensors.h:51
RankTwoTensorTempl< Real >
StressDivergenceTensors::_temp_var
const unsigned int _temp_var
Definition: StressDivergenceTensors.h:69