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 "ComputePlaneFiniteStrain.h"
      11             : #include "UserObject.h"
      12             : 
      13             : registerMooseObject("SolidMechanicsApp", ComputePlaneFiniteStrain);
      14             : 
      15             : InputParameters
      16        1416 : ComputePlaneFiniteStrain::validParams()
      17             : {
      18        1416 :   InputParameters params = Compute2DFiniteStrain::validParams();
      19        1416 :   params.addClassDescription("Compute strain increment and rotation increment for finite strain "
      20             :                              "under 2D planar assumptions.");
      21        2832 :   params.addParam<UserObjectName>("subblock_index_provider",
      22             :                                   "SubblockIndexProvider user object name");
      23        2832 :   params.addCoupledVar("scalar_out_of_plane_strain",
      24             :                        "Scalar variable for generalized plane strain");
      25        2832 :   params.addCoupledVar("out_of_plane_strain", "Nonlinear variable for plane stress condition");
      26             : 
      27        1416 :   return params;
      28           0 : }
      29             : 
      30        1062 : ComputePlaneFiniteStrain::ComputePlaneFiniteStrain(const InputParameters & parameters)
      31             :   : Compute2DFiniteStrain(parameters),
      32        3186 :     _subblock_id_provider(isParamValid("subblock_index_provider")
      33        1062 :                               ? &getUserObject<SubblockIndexProvider>("subblock_index_provider")
      34             :                               : nullptr),
      35        1062 :     _scalar_out_of_plane_strain_coupled(isCoupledScalar("scalar_out_of_plane_strain")),
      36        1062 :     _out_of_plane_strain_coupled(isCoupled("out_of_plane_strain")),
      37        1062 :     _out_of_plane_strain(_out_of_plane_strain_coupled ? coupledValue("out_of_plane_strain")
      38             :                                                       : _zero),
      39        1062 :     _out_of_plane_strain_old(_out_of_plane_strain_coupled ? coupledValueOld("out_of_plane_strain")
      40        1062 :                                                           : _zero)
      41             : {
      42        1062 :   if (_out_of_plane_strain_coupled && _scalar_out_of_plane_strain_coupled)
      43           0 :     mooseError("Must define only one of out_of_plane_strain or scalar_out_of_plane_strain");
      44             : 
      45        1062 :   if (_scalar_out_of_plane_strain_coupled)
      46             :   {
      47         126 :     const auto nscalar_strains = coupledScalarComponents("scalar_out_of_plane_strain");
      48         126 :     _scalar_out_of_plane_strain.resize(nscalar_strains);
      49         126 :     _scalar_out_of_plane_strain_old.resize(nscalar_strains);
      50         252 :     for (unsigned int i = 0; i < nscalar_strains; ++i)
      51             :     {
      52         126 :       _scalar_out_of_plane_strain[i] = &coupledScalarValue("scalar_out_of_plane_strain", i);
      53         126 :       _scalar_out_of_plane_strain_old[i] = &coupledScalarValueOld("scalar_out_of_plane_strain", i);
      54             :     }
      55             :   }
      56        1062 : }
      57             : 
      58             : Real
      59    10653006 : ComputePlaneFiniteStrain::computeOutOfPlaneGradDisp()
      60             : {
      61             :   /**
      62             :    * This is consistent with the approximation of stretch rate tensor
      63             :    * D = log(sqrt(Fhat^T * Fhat)) / dt
      64             :    */
      65    10653006 :   if (_scalar_out_of_plane_strain_coupled)
      66      116356 :     return std::exp((*_scalar_out_of_plane_strain[getCurrentSubblockIndex()])[0]) - 1.0;
      67             :   else
      68    10536650 :     return std::exp(_out_of_plane_strain[_qp]) - 1.0;
      69             : }
      70             : 
      71             : Real
      72    10653006 : ComputePlaneFiniteStrain::computeOutOfPlaneGradDispOld()
      73             : {
      74    10653006 :   if (_scalar_out_of_plane_strain_coupled)
      75      116356 :     return std::exp((*_scalar_out_of_plane_strain_old[getCurrentSubblockIndex()])[0]) - 1.0;
      76             :   else
      77    10536650 :     return std::exp(_out_of_plane_strain_old[_qp]) - 1.0;
      78             : }