Line data    Source code

       1             : //* This file is part of the MOOSE framework
       2             : //* https://mooseframework.inl.gov
       3             : //*
       4             : //* All rights reserved, see COPYRIGHT for full restrictions
       5             : //* https://github.com/idaholab/moose/blob/master/COPYRIGHT
       6             : //*
       7             : //* Licensed under LGPL 2.1, please see LICENSE for details
       8             : //* https://www.gnu.org/licenses/lgpl-2.1.html
       9             : 
      10             : #include "ComputeIncrementalStrainBase.h"
      11             : #include "MooseMesh.h"
      12             : 
      13             : InputParameters
      14       18736 : ComputeIncrementalStrainBase::validParams()
      15             : {
      16       18736 :   InputParameters params = ComputeStrainBase::validParams();
      17       18736 :   return params;
      18             : }
      19             : 
      20       14046 : ComputeIncrementalStrainBase::ComputeIncrementalStrainBase(const InputParameters & parameters)
      21             :   : ComputeStrainBase(parameters),
      22       28092 :     _grad_disp_old(3),
      23       14046 :     _strain_rate(declareProperty<RankTwoTensor>(_base_name + "strain_rate")),
      24       14046 :     _strain_increment(declareProperty<RankTwoTensor>(_base_name + "strain_increment")),
      25       14046 :     _rotation_increment(declareProperty<RankTwoTensor>(_base_name + "rotation_increment")),
      26       14046 :     _deformation_gradient(declareProperty<RankTwoTensor>(_base_name + "deformation_gradient")),
      27       28092 :     _mechanical_strain_old(getMaterialPropertyOld<RankTwoTensor>(_base_name + "mechanical_strain")),
      28       28092 :     _total_strain_old(getMaterialPropertyOld<RankTwoTensor>(_base_name + "total_strain")),
      29       28092 :     _eigenstrains_old(_eigenstrain_names.size())
      30             : {
      31       16860 :   for (unsigned int i = 0; i < _eigenstrains_old.size(); ++i)
      32        2814 :     _eigenstrains_old[i] = &getMaterialPropertyOld<RankTwoTensor>(_eigenstrain_names[i]);
      33       14046 : }
      34             : 
      35             : void
      36       12318 : ComputeIncrementalStrainBase::initialSetup()
      37             : {
      38       12318 :   ComputeStrainBase::initialSetup();
      39       49272 :   for (unsigned int i = 0; i < 3; ++i)
      40             :   {
      41       36954 :     if (_fe_problem.isTransient() && i < _ndisp)
      42       34974 :       _grad_disp_old[i] = &coupledGradientOld("displacements", i);
      43             :     else
      44        1980 :       _grad_disp_old[i] = &_grad_zero;
      45             :   }
      46       12318 : }
      47             : 
      48             : void
      49     1089816 : ComputeIncrementalStrainBase::initQpStatefulProperties()
      50             : {
      51     1089816 :   _mechanical_strain[_qp].zero();
      52     1089816 :   _total_strain[_qp].zero();
      53     1089816 :   _deformation_gradient[_qp].setToIdentity();
      54     1089816 :   _rotation_increment[_qp].setToIdentity();
      55     1089816 : }
      56             : 
      57             : void
      58    51259350 : ComputeIncrementalStrainBase::subtractEigenstrainIncrementFromStrain(RankTwoTensor & strain)
      59             : {
      60    54787240 :   for (unsigned int i = 0; i < _eigenstrains.size(); ++i)
      61             :   {
      62     3527890 :     strain -= (*_eigenstrains[i])[_qp];
      63     3527890 :     strain += (*_eigenstrains_old[i])[_qp];
      64             :   }
      65    51259350 : }