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 "ADComputeRSphericalIncrementalStrain.h"
      11             : #include "Assembly.h"
      12             : #include "FEProblem.h"
      13             : #include "MooseMesh.h"
      14             : 
      15             : #include "libmesh/quadrature.h"
      16             : 
      17             : registerMooseObject("TensorMechanicsApp", ADComputeRSphericalIncrementalStrain);
      18             : 
      19             : InputParameters
      20          48 : ADComputeRSphericalIncrementalStrain::validParams()
      21             : {
      22          48 :   InputParameters params = ADComputeIncrementalSmallStrain::validParams();
      23          48 :   params.addClassDescription(
      24             :       "Compute a strain increment for incremental strains in 1D spherical symmetry problems.");
      25          48 :   return params;
      26           0 : }
      27             : 
      28          36 : ADComputeRSphericalIncrementalStrain::ADComputeRSphericalIncrementalStrain(
      29          36 :     const InputParameters & parameters)
      30          36 :   : ADComputeIncrementalSmallStrain(parameters), _disp_old_0(coupledValueOld("displacements", 0))
      31             : {
      32          36 : }
      33             : 
      34             : void
      35          36 : ADComputeRSphericalIncrementalStrain::initialSetup()
      36             : {
      37          36 :   ADComputeIncrementalStrainBase::initialSetup();
      38             : 
      39          36 :   const auto & subdomainIDs = _mesh.meshSubdomains();
      40          72 :   for (auto subdomainID : subdomainIDs)
      41          36 :     if (_fe_problem.getCoordSystem(subdomainID) != Moose::COORD_RSPHERICAL)
      42           0 :       mooseError("The coordinate system must be set to RSPHERICAL for 1D R spherical simulations.");
      43          36 : }
      44             : 
      45             : void
      46         726 : ADComputeRSphericalIncrementalStrain::computeTotalStrainIncrement(
      47             :     ADRankTwoTensor & total_strain_increment)
      48             : {
      49             :   // Deformation gradient calculation in cylindrical coordinates
      50         726 :   ADRankTwoTensor A;  // Deformation gradient
      51         726 :   RankTwoTensor Fbar; // Old Deformation gradient
      52             : 
      53             :   // Step through calculating the current and old deformation gradients
      54             :   // Only diagonal components are nonzero because this is a 1D material
      55             :   // Note: x_disp is the radial displacement
      56         726 :   A(0, 0) = (*_grad_disp[0])[_qp](0);
      57         726 :   Fbar(0, 0) = (*_grad_disp_old[0])[_qp](0);
      58             : 
      59             :   // The polar and azimuthal strains are functions of radial displacement
      60         726 :   if (!MooseUtils::relativeFuzzyEqual(_q_point[_qp](0), 0.0))
      61             :   {
      62        1452 :     A(1, 1) = (*_disp[0])[_qp] / _q_point[_qp](0);
      63         726 :     Fbar(1, 1) = _disp_old_0[_qp] / _q_point[_qp](0);
      64             :   }
      65             : 
      66             :   // The polar and azimuthal strains are equivalent in this 1D problem
      67         726 :   A(2, 2) = A(1, 1);
      68         726 :   Fbar(2, 2) = Fbar(1, 1);
      69             : 
      70             :   // very nearly A = gradU - gradUold, adapted to cylindrical coords
      71         726 :   A -= Fbar;
      72             : 
      73        1452 :   total_strain_increment = 0.5 * (A + A.transpose());
      74         726 : }