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       12156 : ComputeIncrementalStrainBase::validParams()
      15             : {
      16       12156 :   InputParameters params = ComputeStrainBase::validParams();
      17       12156 :   return params;
      18             : }
      19             : 
      20        9114 : ComputeIncrementalStrainBase::ComputeIncrementalStrainBase(const InputParameters & parameters)
      21             :   : ComputeStrainBase(parameters),
      22       18228 :     _grad_disp_old(3),
      23        9114 :     _strain_rate(declareProperty<RankTwoTensor>(_base_name + "strain_rate")),
      24        9114 :     _strain_increment(declareProperty<RankTwoTensor>(_base_name + "strain_increment")),
      25        9114 :     _rotation_increment(declareProperty<RankTwoTensor>(_base_name + "rotation_increment")),
      26        9114 :     _deformation_gradient(declareProperty<RankTwoTensor>(_base_name + "deformation_gradient")),
      27       18228 :     _mechanical_strain_old(getMaterialPropertyOld<RankTwoTensor>(_base_name + "mechanical_strain")),
      28       18228 :     _total_strain_old(getMaterialPropertyOld<RankTwoTensor>(_base_name + "total_strain")),
      29       18228 :     _eigenstrains_old(_eigenstrain_names.size())
      30             : {
      31       10893 :   for (unsigned int i = 0; i < _eigenstrains_old.size(); ++i)
      32        1779 :     _eigenstrains_old[i] = &getMaterialPropertyOld<RankTwoTensor>(_eigenstrain_names[i]);
      33        9114 : }
      34             : 
      35             : void
      36        8007 : ComputeIncrementalStrainBase::initialSetup()
      37             : {
      38        8007 :   ComputeStrainBase::initialSetup();
      39       32028 :   for (unsigned int i = 0; i < 3; ++i)
      40             :   {
      41       24021 :     if (_fe_problem.isTransient() && i < _ndisp)
      42       22752 :       _grad_disp_old[i] = &coupledGradientOld("displacements", i);
      43             :     else
      44        1269 :       _grad_disp_old[i] = &_grad_zero;
      45             :   }
      46        8007 : }
      47             : 
      48             : void
      49      705136 : ComputeIncrementalStrainBase::initQpStatefulProperties()
      50             : {
      51      705136 :   _mechanical_strain[_qp].zero();
      52      705136 :   _total_strain[_qp].zero();
      53      705136 :   _deformation_gradient[_qp].setToIdentity();
      54      705136 :   _rotation_increment[_qp].setToIdentity();
      55      705136 : }
      56             : 
      57             : void
      58    36469945 : ComputeIncrementalStrainBase::subtractEigenstrainIncrementFromStrain(RankTwoTensor & strain)
      59             : {
      60    39779802 :   for (unsigned int i = 0; i < _eigenstrains.size(); ++i)
      61             :   {
      62     3309857 :     strain -= (*_eigenstrains[i])[_qp];
      63     3309857 :     strain += (*_eigenstrains_old[i])[_qp];
      64             :   }
      65    36469945 : }