Nonlinear.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html
 
#pragma once
 
#include "Element.h"
 
// Forward declarations
class MaterialModel;
class VolumetricModel;
 
namespace SolidMechanics
{
 
class Nonlinear : public Element
{
public:
  Nonlinear(SolidModel & solid_model, const std::string & name, const InputParameters & parameters);
 
  virtual ~Nonlinear();
 
  const ColumnMajorMatrix & incrementalRotation() const { return _incremental_rotation; }
 
  const std::vector<ColumnMajorMatrix> & Fhat() const { return _Fhat; }
 
protected:
  enum DecompMethod
  {
    RashidApprox = 0,
    Eigen = 1
  };
 
  DecompMethod _decomp_method;
 
  ColumnMajorMatrix _incremental_rotation;
  ColumnMajorMatrix _Uhat;
 
  std::vector<ColumnMajorMatrix> _Fhat;
  std::vector<ColumnMajorMatrix> _Fbar;
  ColumnMajorMatrix _F;
 
  virtual void init();
 
  virtual void computeStrain(const unsigned qp,
                             const SymmTensor & total_strain_old,
                             SymmTensor & total_strain_new,
                             SymmTensor & strain_increment);
 
  virtual Real volumeRatioOld(unsigned /*qp*/) const { mooseError("volumeRatioOld not defined"); }
 
  virtual void finalizeStress(std::vector<SymmTensor *> & t);
 
  virtual void computeIncrementalDeformationGradient(std::vector<ColumnMajorMatrix> & Fhat) = 0;
  void computeStrainIncrement(const ColumnMajorMatrix & Fhat, SymmTensor & strain_increment);
  void computePolarDecomposition(const ColumnMajorMatrix & Fhat);
 
  void computeStrainAndRotationIncrement(const ColumnMajorMatrix & Fhat,
                                         SymmTensor & strain_increment);
};
 
} // namespace solid_mechanics