ComputeReducedOrderEigenstrain.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* 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 "DenseMatrix.h"
#include "ComputeEigenstrainBase.h"
#include "RankTwoTensor.h"

// libmesh includes
#include "libmesh/dense_vector.h"

class Assembly;
class InputParameters;
class MooseObject;
class SubProblem;

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;
};