Line data    Source code

       1             : //* This file is part of the MOOSE framework
       2             : //* https://mooseframework.inl.gov
       3             : //*
       4             : //* All rights reserved, see COPYRIGHT for full restrictions
       5             : //* https://github.com/idaholab/moose/blob/master/COPYRIGHT
       6             : //*
       7             : //* Licensed under LGPL 2.1, please see LICENSE for details
       8             : //* https://www.gnu.org/licenses/lgpl-2.1.html
       9             : 
      10             : #include "ComputeCosseratElasticityTensor.h"
      11             : 
      12             : registerMooseObject("SolidMechanicsApp", ComputeCosseratElasticityTensor);
      13             : 
      14             : InputParameters
      15         272 : ComputeCosseratElasticityTensor::validParams()
      16             : {
      17         272 :   InputParameters params = ComputeElasticityTensorBase::validParams();
      18         272 :   params.addClassDescription("Compute Cosserat elasticity and flexural bending rigidity tensors");
      19         544 :   params.addRequiredParam<std::vector<Real>>("E_ijkl", "Elastic stiffness tensor for material");
      20         544 :   params.addParam<MooseEnum>(
      21         544 :       "fill_method", RankFourTensor::fillMethodEnum() = "symmetric9", "The fill method");
      22         544 :   params.addRequiredParam<std::vector<Real>>("B_ijkl", "Flexural bending rigidity tensor.");
      23         544 :   params.addParam<MooseEnum>("fill_method_bending",
      24         544 :                              RankFourTensor::fillMethodEnum() = "antisymmetric_isotropic",
      25             :                              "The fill method for the 'bending' tensor.");
      26         272 :   return params;
      27           0 : }
      28             : 
      29         204 : ComputeCosseratElasticityTensor::ComputeCosseratElasticityTensor(const InputParameters & parameters)
      30             :   : ComputeElasticityTensorBase(parameters),
      31         612 :     _Eijkl(getParam<std::vector<Real>>("E_ijkl"),
      32         408 :            (RankFourTensor::FillMethod)(int)getParam<MooseEnum>("fill_method")),
      33         612 :     _Bijkl(getParam<std::vector<Real>>("B_ijkl"),
      34         204 :            (RankFourTensor::FillMethod)(int)getParam<MooseEnum>("fill_method_bending")),
      35         204 :     _elastic_flexural_rigidity_tensor(
      36         408 :         declareProperty<RankFourTensor>("elastic_flexural_rigidity_tensor"))
      37             : {
      38         408 :   if (!isParamValid("elasticity_tensor_prefactor"))
      39         408 :     issueGuarantee(_elasticity_tensor_name, Guarantee::CONSTANT_IN_TIME);
      40         204 : }
      41             : 
      42             : void
      43       54240 : ComputeCosseratElasticityTensor::computeQpElasticityTensor()
      44             : {
      45       54240 :   _elasticity_tensor[_qp] = _Eijkl;
      46       54240 :   _elastic_flexural_rigidity_tensor[_qp] = _Bijkl;
      47       54240 : }