LynxStrainRatioAux.C 2.17 KB
Newer Older
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
/******************************************************************************/
/*                       LYNX, a MOOSE-based application                      */
/*                                                                            */
/*          Copyright (C) 2017 by Antoine B. Jacquey and Mauro Cacace         */
/*             GFZ Potsdam, German Research Centre for Geosciences            */
/*                                                                            */
/*    This program is free software: you can redistribute it and/or modify    */
/*    it under the terms of the GNU General Public License as published by    */
/*      the Free Software Foundation, either version 3 of the License, or     */
/*                     (at your option) any later version.                    */
/*                                                                            */
/*       This program is distributed in the hope that it will be useful,      */
/*       but WITHOUT ANY WARRANTY; without even the implied warranty of       */
/*        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the       */
/*                GNU General Public License for more details.                */
/*                                                                            */
/*      You should have received a copy of the GNU General Public License     */
/*    along with this program. If not, see <http://www.gnu.org/licenses/>     */
/******************************************************************************/
20

21
#include "LynxStrainRatioAux.h"
22

23
registerMooseObject("LynxApp", LynxStrainRatioAux);
24

25 26 27 28
template <>
InputParameters
validParams<LynxStrainRatioAux>()
{
29
  InputParameters params = validParams<LynxElasticStrainAuxBase>();
30 31 32 33
  params.addClassDescription(
      "Access the strain ratio (volumetric strain on norm of the elastic strain).");
  return params;
}
34

35
LynxStrainRatioAux::LynxStrainRatioAux(const InputParameters & parameters)
36
  : LynxElasticStrainAuxBase(parameters)
37 38
{
}
39

40 41 42 43
Real
LynxStrainRatioAux::computeValue()
{
  Real strain_norm = _elastic_strain[_qp].L2norm();
44

45 46 47 48 49
  if (strain_norm != 0.0)
    return _elastic_strain[_qp].trace() / strain_norm;
  else
    return -std::sqrt(3.0);
}