LynxMass.h 2.21 KB
Newer Older
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56
/******************************************************************************/
/*                       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/>     */
/******************************************************************************/

#ifndef LYNXMASS_H
#define LYNXMASS_H

#include "Kernel.h"
#include "DerivativeMaterialInterface.h"

class LynxMass;

template <>
InputParameters validParams<LynxMass>();

class LynxMass : public DerivativeMaterialInterface<Kernel>
{
public:
  LynxMass(const InputParameters & parameters);
  virtual ~LynxMass() {}
  enum PenaltyType
  {
    LINEAR,
    LAPLACE
  };

protected:
  virtual Real computeQpResidual();
  virtual Real computeQpJacobian();
  virtual Real computeQpOffDiagJacobian(unsigned jvar);

  PenaltyType _penalty_type;

  unsigned _ndisp;
  std::vector<unsigned> _disp_var;

  const MaterialProperty<Real> & _penalty;
  const MaterialProperty<RankTwoTensor> & _strain_increment;
};

57
#endif // LYNXMASS_H