PoroMechanicsAction.C

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
 
#include "PoroMechanicsAction.h"
 
#include "Factory.h"
#include "FEProblem.h"
#include "Parser.h"
#include "Conversion.h"
 
registerMooseAction("SolidMechanicsApp", PoroMechanicsAction, "setup_mesh_complete");
 
registerMooseAction("SolidMechanicsApp", PoroMechanicsAction, "validate_coordinate_systems");
 
registerMooseAction("SolidMechanicsApp", PoroMechanicsAction, "add_kernel");
 
InputParameters
PoroMechanicsAction::validParams()
{
  InputParameters params = QuasiStaticSolidMechanicsPhysics::validParams();
  params.addRequiredParam<VariableName>("porepressure", "The pore pressure variable");
  params.addClassDescription("Adds the poro-mechanics coupling term");
  return params;
}
 
PoroMechanicsAction::PoroMechanicsAction(const InputParameters & params)
  : QuasiStaticSolidMechanicsPhysics(params)
{
  if (_use_ad)
    paramError("use_ad", "AD not setup for use with PoroMechanicsAction");
}
 
void
PoroMechanicsAction::act()
{
  QuasiStaticSolidMechanicsPhysics::act();
 
  if (_current_task == "add_kernel")
  {
    // Prepare displacements and set value for dim
    std::vector<VariableName> displacements = getParam<std::vector<VariableName>>("displacements");
    unsigned int dim = displacements.size();
 
    // all the kernels added below have pore pressure as a coupled variable
    // add this to the kernel's params
    std::string type("PoroMechanicsCoupling");
    InputParameters params = _factory.getValidParams(type);
    VariableName pp_var(getParam<VariableName>("porepressure"));
    params.addCoupledVar("porepressure", "");
    params.set<std::vector<VariableName>>("porepressure") = {pp_var};
 
    // now add the kernels
    for (unsigned int i = 0; i < dim; ++i)
    {
      std::string kernel_name = "PoroMechanics" + Moose::stringify(i);
 
      params.set<unsigned int>("component") = i;
      params.set<NonlinearVariableName>("variable") = displacements[i];
 
      _problem->addKernel(type, kernel_name, params);
    }
  }
}