ComputeReducedOrderEigenstrain.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 "ComputeEigenstrainBase.h"
#include "RankTwoTensor.h"
 
// libmesh includes
#include "libmesh/dense_matrix.h"
#include "libmesh/dense_vector.h"
 
class Assembly;
class InputParameters;
class MooseObject;
class SubProblem;
 
class ComputeReducedOrderEigenstrain;
 
template <>
InputParameters validParams<ComputeReducedOrderEigenstrain>();
 
class ComputeReducedOrderEigenstrain : public ComputeEigenstrainBase
{
public:
  static InputParameters validParams();
 
  ComputeReducedOrderEigenstrain(const InputParameters & parameters);
 
  void initQpStatefulProperties();
  void computeProperties();
  void computeQpEigenstrain();
 
private:
  void applyEigenstrain(MaterialProperty<RankTwoTensor> & strain);
 
  void sumEigenstrain();
 
  void prepareEigenstrain();
 
  std::vector<MaterialPropertyName> _input_eigenstrain_names;
  std::vector<const MaterialProperty<RankTwoTensor> *> _eigenstrains;
 
  SubProblem & _subproblem;
  const unsigned _ncols;
  const bool _second_order;
  std::vector<RankTwoTensor> _eigsum;
  DenseMatrix<Real> _A;
  std::vector<DenseVector<Real>> _b;
  DenseMatrix<Real> _AT;
  DenseVector<Real> _ATb;
  std::vector<DenseVector<Real>> _x;
  std::vector<Real> _vals;
  RankTwoTensor _adjusted_eigenstrain;
};
 
template <>
InputParameters validParams<ComputeReducedOrderEigenstrain>();