PyCElasticityTensor

Computes PyC elasticity tensor

Description

The PyC material is assumed to be transverse isotropic. The Poisson's ratio and elastic modulus of PyC are given by Miller et al. (2018) and Ho (1993). The radial and tangential components, and , of the elastic modulus (GPa) are computed as: where (kg/m) is the density, (dimensionless) is the as-fabricated anisotropy, (= 30 Angstroms ) is the crystallite diameter, (10 n/m, >0.18 MeV) is the fast neutron fluence, and (C) is the temperature. The default value of is set to be 30 Angstroms such that the factor of that term takes a value of 1. The Poisson's ratios and are assumed to have a constant value of 0.33 and the Poisson's ratio is calculated to be .

The correlation is given for fast fluence with neutron energy threshold >0.18 MeV. The model performs the fast fluence conversion from >0.10 MeV to >0.18 MeV using the flux_conversion_factor parameter. The shear modulus for the planes perpendicular to the PyC layer are assumed to be equal to the . The average Young's modulus is computed as and the average Poisson's ratio is computed as .

Instead of the default anisotropic model, the user can also choose the isotropic model. For this isotropic model, the Poisson's ratio is set to a constant value of 0.33 and the elastic modulus, (GPa) is computed as . The user should set anisotropy = false to use the isotropic model.

These models are assumed valid at all temperatures, density: 1.8 - 2.0 g/cm, fast fluence: 0 - 3.9610 n/m (>0.18 MeV), and all Bacon Anisotropy Factor (BAF) values. The function is evaluated for all density values below 2.25 g/cm (theoretical density), including outside validity range. It is calculated with density value of 2.25 g/cm for densities exceeding the theoretical density, but results should be considered invalid. The functionality is bounded by the fast fluence value of 3.9610 n/m.

Example Input Syntax

[Materials<<<{"href": "../../../syntax/Materials/index.html"}>>>]
  [elasticity_tensor]
    type = PyCElasticityTensor<<<{"description": "Computes PyC elasticity tensor", "href": "PyCElasticityTensor.html"}>>>
    block<<<{"description": "The list of blocks (ids or names) that this object will be applied"}>>> = PyC_quad
    initial_BAF<<<{"description": "Initial Bacon Anisotropy Factor (BAF) must be greater than or equal to 1.0."}>>> = 1.05
    temperature<<<{"description": "Coupled temperature (K)"}>>> = temperature
  []
[]
(test/tests/triso/pyc_elasticity_tensor/anisotropic_exact.i)

Input Parameters

  • flux_conversion_factorConvert fast neutron flux E>0.10 to E>0.18 MeV

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Convert fast neutron flux E>0.10 to E>0.18 MeV

  • temperatureCoupled temperature (K)

    C++ Type:std::vector<VariableName>

    Unit:(no unit assumed)

    Controllable:No

    Description:Coupled temperature (K)

Required Parameters

  • anisotropyTruetrue: use anisotropic PyC elasticity model; false: use isotropic PyC elasticity model.

    Default:True

    C++ Type:bool

    Controllable:No

    Description:true: use anisotropic PyC elasticity model; false: use isotropic PyC elasticity model.

  • base_nameOptional parameter that allows the user to define multiple mechanics material systems on the same block, i.e. for multiple phases

    C++ Type:std::string

    Controllable:No

    Description:Optional parameter that allows the user to define multiple mechanics material systems on the same block, i.e. for multiple phases

  • blockThe list of blocks (ids or names) that this object will be applied

    C++ Type:std::vector<SubdomainName>

    Controllable:No

    Description:The list of blocks (ids or names) that this object will be applied

  • boundaryThe list of boundaries (ids or names) from the mesh where this object applies

    C++ Type:std::vector<BoundaryName>

    Controllable:No

    Description:The list of boundaries (ids or names) from the mesh where this object applies

  • computeTrueWhen false, MOOSE will not call compute methods on this material. The user must call computeProperties() after retrieving the MaterialBase via MaterialBasePropertyInterface::getMaterialBase(). Non-computed MaterialBases are not sorted for dependencies.

    Default:True

    C++ Type:bool

    Controllable:No

    Description:When false, MOOSE will not call compute methods on this material. The user must call computeProperties() after retrieving the MaterialBase via MaterialBasePropertyInterface::getMaterialBase(). Non-computed MaterialBases are not sorted for dependencies.

  • constant_onNONEWhen ELEMENT, MOOSE will only call computeQpProperties() for the 0th quadrature point, and then copy that value to the other qps.When SUBDOMAIN, MOOSE will only call computeQpProperties() for the 0th quadrature point, and then copy that value to the other qps. Evaluations on element qps will be skipped

    Default:NONE

    C++ Type:MooseEnum

    Options:NONE, ELEMENT, SUBDOMAIN

    Controllable:No

    Description:When ELEMENT, MOOSE will only call computeQpProperties() for the 0th quadrature point, and then copy that value to the other qps.When SUBDOMAIN, MOOSE will only call computeQpProperties() for the 0th quadrature point, and then copy that value to the other qps. Evaluations on element qps will be skipped

  • crystallite_diameter3e-09Crystallite diameter (m)

    Default:3e-09

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Crystallite diameter (m)

  • declare_suffixAn optional suffix parameter that can be appended to any declared properties. The suffix will be prepended with a '_' character.

    C++ Type:MaterialPropertyName

    Unit:(no unit assumed)

    Controllable:No

    Description:An optional suffix parameter that can be appended to any declared properties. The suffix will be prepended with a '_' character.

  • elastic_constants_scale_factor1Scale factor to be applied to the elastic constants for anisotropic model

    Default:1

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Scale factor to be applied to the elastic constants for anisotropic model

  • elasticity_tensor_prefactorOptional function to use as a scalar prefactor on the elasticity tensor.

    C++ Type:FunctionName

    Unit:(no unit assumed)

    Controllable:No

    Description:Optional function to use as a scalar prefactor on the elasticity tensor.

  • fast_neutron_fluencefast_neutron_fluenceCoupled fast (E>0.10 MeV) neutron fluence

    Default:fast_neutron_fluence

    C++ Type:MaterialPropertyName

    Unit:(no unit assumed)

    Controllable:No

    Description:Coupled fast (E>0.10 MeV) neutron fluence

  • initial_BAF1Initial Bacon Anisotropy Factor (BAF) must be greater than or equal to 1.0.

    Default:1

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Initial Bacon Anisotropy Factor (BAF) must be greater than or equal to 1.0.

  • poissons_ratio0.33Poisson's ratio

    Default:0.33

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Poisson's ratio

Optional Parameters

  • control_tagsAdds user-defined labels for accessing object parameters via control logic.

    C++ Type:std::vector<std::string>

    Controllable:No

    Description:Adds user-defined labels for accessing object parameters via control logic.

  • enableTrueSet the enabled status of the MooseObject.

    Default:True

    C++ Type:bool

    Controllable:Yes

    Description:Set the enabled status of the MooseObject.

  • implicitTrueDetermines whether this object is calculated using an implicit or explicit form

    Default:True

    C++ Type:bool

    Controllable:No

    Description:Determines whether this object is calculated using an implicit or explicit form

  • seed0The seed for the master random number generator

    Default:0

    C++ Type:unsigned int

    Controllable:No

    Description:The seed for the master random number generator

  • use_displaced_meshFalseWhether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.

    Default:False

    C++ Type:bool

    Controllable:No

    Description:Whether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.

Advanced Parameters

  • output_propertiesList of material properties, from this material, to output (outputs must also be defined to an output type)

    C++ Type:std::vector<std::string>

    Controllable:No

    Description:List of material properties, from this material, to output (outputs must also be defined to an output type)

  • outputsnone Vector of output names where you would like to restrict the output of variables(s) associated with this object

    Default:none

    C++ Type:std::vector<OutputName>

    Controllable:No

    Description:Vector of output names where you would like to restrict the output of variables(s) associated with this object

Outputs Parameters

  • poissons_ratio_scale_factor1Scale factor to be applied to the Poisson's ratio for isotropic model

    Default:1

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Scale factor to be applied to the Poisson's ratio for isotropic model

  • youngs_modulus_scale_factor1Scale factor to be applied to the Young's modulus for isotropic model

    Default:1

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Scale factor to be applied to the Young's modulus for isotropic model

Advanced: Scaling Factors Parameters

  • prop_getter_suffixAn optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.

    C++ Type:MaterialPropertyName

    Unit:(no unit assumed)

    Controllable:No

    Description:An optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.

  • use_interpolated_stateFalseFor the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.

    Default:False

    C++ Type:bool

    Controllable:No

    Description:For the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.

Material Property Retrieval Parameters

Input Files

References

  1. F. Ho. NP-MHTGR: Material Models of Pyrocarbon and Pyrolytic Silicon Carbide. Report CEGA-002820 Rev. 1, CEGA Corporation, July 1993.[BibTeX]
  2. G.K. Miller, D.A. Petti, J.T. Maki, D.L. Knudson, and W.F. Skerjanc. PARFUME Theory and Model Basis Report. Report INL/EXT-08-14497 (Rev.1), Idaho National Laboratory, September 2018.[BibTeX]