Line data    Source code

       1             : //* This file is part of the MOOSE framework
       2             : //* https://www.mooseframework.org
       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        8860 : ComputeIncrementalStrainBase::validParams()
      15             : {
      16        8860 :   InputParameters params = ComputeStrainBase::validParams();
      17        8860 :   return params;
      18             : }
      19             : 
      20        6642 : ComputeIncrementalStrainBase::ComputeIncrementalStrainBase(const InputParameters & parameters)
      21             :   : ComputeStrainBase(parameters),
      22       13284 :     _grad_disp_old(3),
      23        6642 :     _strain_rate(declareProperty<RankTwoTensor>(_base_name + "strain_rate")),
      24        6642 :     _strain_increment(declareProperty<RankTwoTensor>(_base_name + "strain_increment")),
      25        6642 :     _rotation_increment(declareProperty<RankTwoTensor>(_base_name + "rotation_increment")),
      26        6642 :     _deformation_gradient(declareProperty<RankTwoTensor>(_base_name + "deformation_gradient")),
      27       13284 :     _mechanical_strain_old(getMaterialPropertyOld<RankTwoTensor>(_base_name + "mechanical_strain")),
      28       13284 :     _total_strain_old(getMaterialPropertyOld<RankTwoTensor>(_base_name + "total_strain")),
      29       13284 :     _eigenstrains_old(_eigenstrain_names.size())
      30             : {
      31        8031 :   for (unsigned int i = 0; i < _eigenstrains_old.size(); ++i)
      32        1389 :     _eigenstrains_old[i] = &getMaterialPropertyOld<RankTwoTensor>(_eigenstrain_names[i]);
      33        6642 : }
      34             : 
      35             : void
      36        5817 : ComputeIncrementalStrainBase::initialSetup()
      37             : {
      38        5817 :   ComputeStrainBase::initialSetup();
      39       23268 :   for (unsigned int i = 0; i < 3; ++i)
      40             :   {
      41       17451 :     if (_fe_problem.isTransient() && i < _ndisp)
      42       16506 :       _grad_disp_old[i] = &coupledGradientOld("displacements", i);
      43             :     else
      44         945 :       _grad_disp_old[i] = &_grad_zero;
      45             :   }
      46        5817 : }
      47             : 
      48             : void
      49      523208 : ComputeIncrementalStrainBase::initQpStatefulProperties()
      50             : {
      51      523208 :   _mechanical_strain[_qp].zero();
      52      523208 :   _total_strain[_qp].zero();
      53      523208 :   _deformation_gradient[_qp].setToIdentity();
      54      523208 : }
      55             : 
      56             : void
      57    24816513 : ComputeIncrementalStrainBase::subtractEigenstrainIncrementFromStrain(RankTwoTensor & strain)
      58             : {
      59    26501358 :   for (unsigned int i = 0; i < _eigenstrains.size(); ++i)
      60             :   {
      61     1684845 :     strain -= (*_eigenstrains[i])[_qp];
      62     1684845 :     strain += (*_eigenstrains_old[i])[_qp];
      63             :   }
      64    24816513 : }