MFEMLinearElasticityKernel.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

#ifdef MFEM_ENABLED

#pragma once
#include "MFEMKernel.h"

// clang-format off
/*
 * \f[
 * (c_{ikjl} \nabla u_j, \nabla v_i),
 * c_{ikjl} = \lamba \delta_{ik} \delta_{jl} + \mu (\delta_{ij} \delta_{kl} + \delta_{il} \delta_{jk}),
 * \lambda = (E\nu)/((1-2\nu)(1+\nu)),
 * \mu = E/(2(1+\nu)),
 * E is Young's modulus,
 * \nu is Poisson's ratio
 * \f]
*/
// clang-format on
class MFEMLinearElasticityKernel : public MFEMKernel
{
public:
  static InputParameters validParams();

  MFEMLinearElasticityKernel(const InputParameters & parameters);

  virtual mfem::BilinearFormIntegrator * createBFIntegrator() override;

protected:
  const MFEMScalarCoefficientName & _lambda_name;
  const MFEMScalarCoefficientName & _mu_name;
  mfem::Coefficient & _lambda;
  mfem::Coefficient & _mu;
};

#endif