Al6061PlasticityStressUpdate

Computes the stress as a function of temperature and plastic strain from experimentally-based hardening functions. Note: This material model must be run in conjunction with ComputeMultipleInelasticStress.

Description

The Al6061PlasticityStressUpdate model computes the plastic response of Al6061 as a function of fast neutron fluence and temperature. The model is based on recommendations in Rabin et al. (2020).

Model Details

Pages 456-458 of Rabin et al. (2020) give stress vs. strain tables at a set of fast neutron fluence and temperature levels. For each table, it is assumed that Young's modulus is the first non-zero stress entry divided by the first non-zero strain entry. This data is reflected in Al6061ElasticityTensor with the USHPRR option. Likewise, it is assumed for each table that the yield stress is the first non-zero stress entry.

With that in mind, the tables were manipulated to subtract the elastic strain from each strain entry. The elastic strain is simply the stress value divided by Young's modulus. The result is a set of data tables that give stress as a function of plastic strain, with the first entry in each table being the yield stress and a plastic strain of zero.

The model interpolates between fast neutron fluence and temperature values to compute a stress level for a plastic strain. The plastic strain used is determined based on the radial return approach. See IsotropicPlasticityStressUpdate for details.

warning:Data consistency

Data for Young's modulus, yield stress, and plasticity for Al6061 is separately available in Rabin et al. (2020). However, the data is not consistent in that the Young's modulus and yield stress data do not agree with the equivalent data incorporated in the plasticity data. For consistent data, it is highly recommended to use the USHPRR option in the Al6061ElasticityTensor. The inconsistency in data is the reason this model does not use the independent yield stress data in Rabin et al. (2020).

Example Input Syntax

[Materials<<<{"href": "../../../syntax/Materials/index.html"}>>>]
  [plasticity]
    type = Al6061PlasticityStressUpdate<<<{"description": "Computes the stress as a function of temperature and plastic strain from experimentally-based hardening functions. Note: This material model must be run in conjunction with ComputeMultipleInelasticStress.", "href": "Al6061PlasticityStressUpdate.html"}>>>
    property<<<{"description": "List of real-valued material properties to be used in interpolation."}>>> = fluence
    temperature<<<{"description": "Coupled Temperature"}>>> = temperature
    relative_tolerance<<<{"description": "Relative convergence tolerance for Newton iteration"}>>> = 1e-10
  []
[]

Input Parameters

data0 4.6e+25 9.6e+26 1.2e+27 294.15 323.15 423.15 323.15 423.15 0 0.0007766 0.002864 0.00454 0.006872 0.00876 0.01058 0.01306 0.01543 0.02038 0.0238 0.02662 0.03034 0.03441 0.04184 0.04595 0.04976 0.06164 0.0742 0.08081 0.08917 0.09662 0.1043 0.1125 0.1199 0.1305 0.139 0.1459 0.1541 0.162 0.1707 0.1788 0.1899 0.1998 0.207 0.2143 0.2286 0.2387 0.2491 0.2588 0.2698 0.2784 0.2859 0.2966 0.3048 0.3124 0.3177 0.3225 0.3275 0.3319 0.3347 0 0.0002959 0.001094 0.002442 0.003791 0.005689 0.008238 0.01089 0.01369 0.01683 0.02073 0.02363 0.02808 0.03113 0.03428 0.03728 0.04073 0.04423 0.04733 0.05117 0.05497 0.05812 0.06097 0.06397 0.06702 0.07007 0.07327 0.07627 0.08027 0.08337 0.08692 0.08977 0.09257 0.09657 0.09957 0.1026 0.1066 0.1104 0.1135 0.1173 0.1209 0.1243 0.1277 0.1313 0.1354 0.1392 0.1434 0.1474 0.1503 0.1527 0.1546 0.1555 0 0.0002857 0.0007 0.002214 0.004143 0.006229 0.008357 0.01084 0.01397 0.01704 0.02067 0.02351 0.0267 0.02989 0.03323 0.03601 0.03951 0.04271 0.04586 0.04986 0.0527 0.0563 0.0594 0.0627 0.066 0.07 0.073 0.07656 0.08021 0.08331 0.08631 0.09021 0.09371 0.09667 0.09967 0.1026 0.1065 0.1102 0.1137 0.1172 0.1212 0.1242 0.1279 0.1306 0.1346 0.1379 0.1412 0.1445 0.1481 0.1517 0.1548 0.1578 0.1609 0.1641 0.1674 0.1706 0.1723 0.1742 0.1757 0.1767 0.1782 0 0.0001287 0.0006575 0.001791 0.003029 0.00412 0.006953 0.009887 0.01277 0.01591 0.01955 0.02234 0.02523 0.02803 0.03153 0.03452 0.03802 0.04151 0.04511 0.04886 0.05236 0.05541 0.05911 0.06186 0.06577 0.06923 0.07319 0.07691 0.08062 0.08338 0.08705 0.09072 0.09424 0.09852 0.1032 0.1069 0.1101 0.1128 0 0.0003685 0.0009901 0.002318 0.004159 0.006752 0.009546 0.01274 0.01653 0.01953 0.02338 0.02648 0.02978 0.03318 0.03628 0.03968 0.04288 0.04648 0.04973 0.05319 0.05624 0.05994 0.06355 0.0671 0.07055 0.07436 0.07762 0.08159 0.0849 0.08827 0.09329 0.09717 0.1004 0.1034 0.1071 0.1094 0.111 4.7e+07 5.2e+07 5.9e+07 6.4e+07 7e+07 7.7e+07 8.1e+07 8.6e+07 9.1e+07 9.9e+07 1.03e+08 1.07e+08 1.11e+08 1.14e+08 1.2e+08 1.22e+08 1.24e+08 1.29e+08 1.32e+08 1.34e+08 1.35e+08 1.36e+08 1.37e+08 1.38e+08 1.39e+08 1.4e+08 1.4e+08 1.4e+08 1.41e+08 1.41e+08 1.41e+08 1.41e+08 1.41e+08 1.41e+08 1.41e+08 1.41e+08 1.4e+08 1.4e+08 1.39e+08 1.38e+08 1.37e+08 1.35e+08 1.34e+08 1.32e+08 1.3e+08 1.28e+08 1.26e+08 1.24e+08 1.22e+08 1.18e+08 1.13e+08 1.71e+08 1.85e+08 1.93e+08 1.98e+08 2.03e+08 2.07e+08 2.12e+08 2.17e+08 2.21e+08 2.26e+08 2.3e+08 2.34e+08 2.39e+08 2.42e+08 2.45e+08 2.47e+08 2.5e+08 2.52e+08 2.54e+08 2.57e+08 2.59e+08 2.6e+08 2.61e+08 2.63e+08 2.64e+08 2.65e+08 2.65e+08 2.65e+08 2.65e+08 2.65e+08 2.66e+08 2.67e+08 2.67e+08 2.67e+08 2.67e+08 2.66e+08 2.66e+08 2.65e+08 2.63e+08 2.61e+08 2.57e+08 2.53e+08 2.49e+08 2.42e+08 2.34e+08 2.26e+08 2.17e+08 2.06e+08 1.97e+08 1.89e+08 1.82e+08 1.75e+08 1.33e+08 1.42e+08 1.47e+08 1.52e+08 1.55e+08 1.57e+08 1.6e+08 1.62e+08 1.65e+08 1.69e+08 1.72e+08 1.73e+08 1.75e+08 1.77e+08 1.78e+08 1.8e+08 1.8e+08 1.8e+08 1.81e+08 1.81e+08 1.82e+08 1.82e+08 1.82e+08 1.82e+08 1.82e+08 1.82e+08 1.82e+08 1.81e+08 1.8e+08 1.8e+08 1.8e+08 1.8e+08 1.8e+08 1.79e+08 1.79e+08 1.79e+08 1.78e+08 1.77e+08 1.75e+08 1.74e+08 1.73e+08 1.72e+08 1.7e+08 1.68e+08 1.65e+08 1.62e+08 1.59e+08 1.55e+08 1.5e+08 1.44e+08 1.39e+08 1.35e+08 1.29e+08 1.24e+08 1.18e+08 1.12e+08 1.08e+08 1.01e+08 9.2e+07 8.3e+07 6.4e+07 2.75e+08 2.84e+08 2.93e+08 3e+08 3.05e+08 3.09e+08 3.16e+08 3.23e+08 3.29e+08 3.34e+08 3.37e+08 3.41e+08 3.45e+08 3.47e+08 3.49e+08 3.51e+08 3.53e+08 3.55e+08 3.55e+08 3.56e+08 3.56e+08 3.55e+08 3.55e+08 3.54e+08 3.5e+08 3.47e+08 3.42e+08 3.36e+08 3.3e+08 3.25e+08 3.18e+08 3.09e+08 2.99e+08 2.86e+08 2.71e+08 2.57e+08 2.44e+08 2.32e+08 2.22e+08 2.32e+08 2.41e+08 2.48e+08 2.51e+08 2.53e+08 2.55e+08 2.57e+08 2.59e+08 2.59e+08 2.6e+08 2.6e+08 2.6e+08 2.6e+08 2.6e+08 2.6e+08 2.6e+08 2.6e+08 2.59e+08 2.58e+08 2.57e+08 2.57e+08 2.55e+08 2.54e+08 2.53e+08 2.49e+08 2.46e+08 2.41e+08 2.37e+08 2.32e+08 2.24e+08 2.15e+08 2.06e+08 1.98e+08 1.85e+08 1.74e+08 1.63e+08 Vector of vectors of tabulated data. The first vector is the outermost set of values for interpolation. If there are only two vectors, the second vector is the tabulated output values. Otherwise, the number of the next set of vectors equals the number of entries in the first vector. This pattern is repeated, with a number of vectors associated with each entry in a preceding vector, until the tabulated output values are given. In every case, the number of tabulated output vectors must equal the number of vectors in the preceding set of values.
Default:0 4.6e+25 9.6e+26 1.2e+27 294.15 323.15 423.15 323.15 423.15 0 0.0007766 0.002864 0.00454 0.006872 0.00876 0.01058 0.01306 0.01543 0.02038 0.0238 0.02662 0.03034 0.03441 0.04184 0.04595 0.04976 0.06164 0.0742 0.08081 0.08917 0.09662 0.1043 0.1125 0.1199 0.1305 0.139 0.1459 0.1541 0.162 0.1707 0.1788 0.1899 0.1998 0.207 0.2143 0.2286 0.2387 0.2491 0.2588 0.2698 0.2784 0.2859 0.2966 0.3048 0.3124 0.3177 0.3225 0.3275 0.3319 0.3347 0 0.0002959 0.001094 0.002442 0.003791 0.005689 0.008238 0.01089 0.01369 0.01683 0.02073 0.02363 0.02808 0.03113 0.03428 0.03728 0.04073 0.04423 0.04733 0.05117 0.05497 0.05812 0.06097 0.06397 0.06702 0.07007 0.07327 0.07627 0.08027 0.08337 0.08692 0.08977 0.09257 0.09657 0.09957 0.1026 0.1066 0.1104 0.1135 0.1173 0.1209 0.1243 0.1277 0.1313 0.1354 0.1392 0.1434 0.1474 0.1503 0.1527 0.1546 0.1555 0 0.0002857 0.0007 0.002214 0.004143 0.006229 0.008357 0.01084 0.01397 0.01704 0.02067 0.02351 0.0267 0.02989 0.03323 0.03601 0.03951 0.04271 0.04586 0.04986 0.0527 0.0563 0.0594 0.0627 0.066 0.07 0.073 0.07656 0.08021 0.08331 0.08631 0.09021 0.09371 0.09667 0.09967 0.1026 0.1065 0.1102 0.1137 0.1172 0.1212 0.1242 0.1279 0.1306 0.1346 0.1379 0.1412 0.1445 0.1481 0.1517 0.1548 0.1578 0.1609 0.1641 0.1674 0.1706 0.1723 0.1742 0.1757 0.1767 0.1782 0 0.0001287 0.0006575 0.001791 0.003029 0.00412 0.006953 0.009887 0.01277 0.01591 0.01955 0.02234 0.02523 0.02803 0.03153 0.03452 0.03802 0.04151 0.04511 0.04886 0.05236 0.05541 0.05911 0.06186 0.06577 0.06923 0.07319 0.07691 0.08062 0.08338 0.08705 0.09072 0.09424 0.09852 0.1032 0.1069 0.1101 0.1128 0 0.0003685 0.0009901 0.002318 0.004159 0.006752 0.009546 0.01274 0.01653 0.01953 0.02338 0.02648 0.02978 0.03318 0.03628 0.03968 0.04288 0.04648 0.04973 0.05319 0.05624 0.05994 0.06355 0.0671 0.07055 0.07436 0.07762 0.08159 0.0849 0.08827 0.09329 0.09717 0.1004 0.1034 0.1071 0.1094 0.111 4.7e+07 5.2e+07 5.9e+07 6.4e+07 7e+07 7.7e+07 8.1e+07 8.6e+07 9.1e+07 9.9e+07 1.03e+08 1.07e+08 1.11e+08 1.14e+08 1.2e+08 1.22e+08 1.24e+08 1.29e+08 1.32e+08 1.34e+08 1.35e+08 1.36e+08 1.37e+08 1.38e+08 1.39e+08 1.4e+08 1.4e+08 1.4e+08 1.41e+08 1.41e+08 1.41e+08 1.41e+08 1.41e+08 1.41e+08 1.41e+08 1.41e+08 1.4e+08 1.4e+08 1.39e+08 1.38e+08 1.37e+08 1.35e+08 1.34e+08 1.32e+08 1.3e+08 1.28e+08 1.26e+08 1.24e+08 1.22e+08 1.18e+08 1.13e+08 1.71e+08 1.85e+08 1.93e+08 1.98e+08 2.03e+08 2.07e+08 2.12e+08 2.17e+08 2.21e+08 2.26e+08 2.3e+08 2.34e+08 2.39e+08 2.42e+08 2.45e+08 2.47e+08 2.5e+08 2.52e+08 2.54e+08 2.57e+08 2.59e+08 2.6e+08 2.61e+08 2.63e+08 2.64e+08 2.65e+08 2.65e+08 2.65e+08 2.65e+08 2.65e+08 2.66e+08 2.67e+08 2.67e+08 2.67e+08 2.67e+08 2.66e+08 2.66e+08 2.65e+08 2.63e+08 2.61e+08 2.57e+08 2.53e+08 2.49e+08 2.42e+08 2.34e+08 2.26e+08 2.17e+08 2.06e+08 1.97e+08 1.89e+08 1.82e+08 1.75e+08 1.33e+08 1.42e+08 1.47e+08 1.52e+08 1.55e+08 1.57e+08 1.6e+08 1.62e+08 1.65e+08 1.69e+08 1.72e+08 1.73e+08 1.75e+08 1.77e+08 1.78e+08 1.8e+08 1.8e+08 1.8e+08 1.81e+08 1.81e+08 1.82e+08 1.82e+08 1.82e+08 1.82e+08 1.82e+08 1.82e+08 1.82e+08 1.81e+08 1.8e+08 1.8e+08 1.8e+08 1.8e+08 1.8e+08 1.79e+08 1.79e+08 1.79e+08 1.78e+08 1.77e+08 1.75e+08 1.74e+08 1.73e+08 1.72e+08 1.7e+08 1.68e+08 1.65e+08 1.62e+08 1.59e+08 1.55e+08 1.5e+08 1.44e+08 1.39e+08 1.35e+08 1.29e+08 1.24e+08 1.18e+08 1.12e+08 1.08e+08 1.01e+08 9.2e+07 8.3e+07 6.4e+07 2.75e+08 2.84e+08 2.93e+08 3e+08 3.05e+08 3.09e+08 3.16e+08 3.23e+08 3.29e+08 3.34e+08 3.37e+08 3.41e+08 3.45e+08 3.47e+08 3.49e+08 3.51e+08 3.53e+08 3.55e+08 3.55e+08 3.56e+08 3.56e+08 3.55e+08 3.55e+08 3.54e+08 3.5e+08 3.47e+08 3.42e+08 3.36e+08 3.3e+08 3.25e+08 3.18e+08 3.09e+08 2.99e+08 2.86e+08 2.71e+08 2.57e+08 2.44e+08 2.32e+08 2.22e+08 2.32e+08 2.41e+08 2.48e+08 2.51e+08 2.53e+08 2.55e+08 2.57e+08 2.59e+08 2.59e+08 2.6e+08 2.6e+08 2.6e+08 2.6e+08 2.6e+08 2.6e+08 2.6e+08 2.6e+08 2.59e+08 2.58e+08 2.57e+08 2.57e+08 2.55e+08 2.54e+08 2.53e+08 2.49e+08 2.46e+08 2.41e+08 2.37e+08 2.32e+08 2.24e+08 2.15e+08 2.06e+08 1.98e+08 1.85e+08 1.74e+08 1.63e+08
C++ Type:std::vector<std::vector<double>>
Unit:(no unit assumed)
Controllable:No
Description:Vector of vectors of tabulated data. The first vector is the outermost set of values for interpolation. If there are only two vectors, the second vector is the tabulated output values. Otherwise, the number of the next set of vectors equals the number of entries in the first vector. This pattern is repeated, with a number of vectors associated with each entry in a preceding vector, until the tabulated output values are given. In every case, the number of tabulated output vectors must equal the number of vectors in the preceding set of values.
propertyList of real-valued material properties to be used in interpolation.
C++ Type:std::vector<MaterialPropertyName>
Unit:(no unit assumed)
Controllable:No
Description:List of real-valued material properties to be used in interpolation.

Required Parameters

absolute_tolerance1e-11Absolute convergence tolerance for Newton iteration
Default:1e-11
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Absolute convergence tolerance for Newton iteration
acceptable_multiplier10Factor applied to relative and absolute tolerance for acceptable convergence if iterations are no longer making progress
Default:10
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Factor applied to relative and absolute tolerance for acceptable convergence if iterations are no longer making progress
adaptive_substeppingFalseUse adaptive substepping, where the number of substeps is successively doubled until the return mapping model successfully converges or the maximum number of substeps is reached.
Default:False
C++ Type:bool
Controllable:No
Description:Use adaptive substepping, where the number of substeps is successively doubled until the return mapping model successfully converges or the maximum number of substeps is reached.
automatic_differentiation_return_mappingFalseWhether to use automatic differentiation to compute the derivative.
Default:False
C++ Type:bool
Controllable:No
Description:Whether to use automatic differentiation to compute the derivative.
base_nameOptional parameter that defines a prefix for all material properties related to this stress update model. This allows for multiple models of the same type to be used without naming conflicts.
C++ Type:std::string
Controllable:No
Description:Optional parameter that defines a prefix for all material properties related to this stress update model. This allows for multiple models of the same type to be used without naming conflicts.
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
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
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.
max_inelastic_increment0.0001The maximum inelastic strain increment allowed in a time step
Default:0.0001
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:The maximum inelastic strain increment allowed in a time step
maximum_number_substeps25The maximum number of substeps allowed before cutting the time step.
Default:25
C++ Type:unsigned int
Controllable:No
Description:The maximum number of substeps allowed before cutting the time step.
relative_tolerance1e-08Relative convergence tolerance for Newton iteration
Default:1e-08
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Relative convergence tolerance for Newton iteration
temperatureCoupled Temperature
C++ Type:std::vector<VariableName>
Unit:(no unit assumed)
Controllable:No
Description:Coupled Temperature
use_substep_integration_errorFalseIf true, it establishes a substep size that will yield, at most,the creep numerical integration error given by substep_strain_tolerance.
Default:False
C++ Type:bool
Controllable:No
Description:If true, it establishes a substep size that will yield, at most,the creep numerical integration error given by substep_strain_tolerance.
use_substeppingNONEWhether and how to use substepping
Default:NONE
C++ Type:MooseEnum
Options:NONE, ERROR_BASED, INCREMENT_BASED
Controllable:No
Description:Whether and how to use substepping

Optional Parameters

apply_strainTrueFlag to apply strain. Used for testing.
Default:True
C++ Type:bool
Controllable:No
Description:Flag to apply strain. Used for testing.
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.
effective_inelastic_strain_nameeffective_plastic_strainName of the material property that stores the effective inelastic strain
Default:effective_plastic_strain
C++ Type:std::string
Controllable:No
Description:Name of the material property that stores the effective inelastic strain
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
substep_strain_tolerance0.1Maximum ratio of the initial elastic strain increment at start of the return mapping solve to the maximum inelastic strain allowable in a single substep. Reduce this value to increase the number of substeps
Default:0.1
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Maximum ratio of the initial elastic strain increment at start of the return mapping solve to the maximum inelastic strain allowable in a single substep. Reduce this value to increase the number of substeps
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

internal_solve_full_iteration_historyFalseSet true to output full internal Newton iteration history at times determined by `internal_solve_output_on`. If false, only a summary is output.
Default:False
C++ Type:bool
Controllable:No
Description:Set true to output full internal Newton iteration history at times determined by `internal_solve_output_on`. If false, only a summary is output.
internal_solve_output_onon_errorWhen to output internal Newton solve information
Default:on_error
C++ Type:MooseEnum
Options:never, on_error, always
Controllable:No
Description:When to output internal Newton solve information

Debug 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

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

(test/tests/solid_mechanics/Al6061_plasticity/plasticity.i)
(test/tests/solid_mechanics/Al6061_plasticity/interpolated.i)

References

B. Rabin, M. Meyer, J. Cole, I. Glagolenko, W. Jones, J-F. Jue, D. Keiser Jr, C. Miller, G. Moore, H. Ozaltun, F. Rice, A. Robinson, J. Smith, D. Wachs, W. Williams, N. Woolstenhulme, G. Hofman, and Y. Kim. Preliminary report on U-Mo monolithic fuel for research reactors. Technical Report INL EXT-17-40975, Idaho National Laboratory, 4 2020.

@techreport{EXT-17-40975R4,
    author = "Rabin, B. and Meyer, M. and Cole, J. and Glagolenko, I. and Jones, W. and Jue, J-F. and Jr, D. Keiser and Miller, C. and Moore, G. and Ozaltun, H. and Rice, F. and Robinson, A. and Smith, J. and Wachs, D. and Williams, W. and Woolstenhulme, N. and Hofman, G. and Kim, Y.",
    title = "Preliminary Report on {U-Mo} Monolithic Fuel for Research Reactors",
    institution = "Idaho National Laboratory",
    number = "INL EXT-17-40975",
    month = "4",
    year = "2020"
}

(test/tests/solid_mechanics/Al6061_plasticity/plasticity.i)

#
# This test checks plasticity based on yield stress and hardening data
# recommended by the USHPRR program.
#
# The mesh involves separate 3D elements, one for each hardening curve
# to be tested.
#
[GlobalParams]
  displacements = 'disp_x disp_y disp_z'
[]

[Mesh]
  use_displaced_mesh = false
  [b0]
    type = GeneratedMeshGenerator
    dim = 3
    nx = 1
    ny = 1
    nz = 1
    xmin = 0
    xmax = 1
    ymin = 0
    ymax = 0.25
    zmin = 0
    zmax = 0.25
  []
  [all]
    type = CombinerGenerator
    inputs = b0
    positions = '0 0 0
                 0 1 0
                 0 2 0
                 0 3 0
                 0 4 0'
  []
  [final]
    type = SubdomainPerElementGenerator
    input = all
    subdomain_ids = '0 1 2 3 4'
  []
[]

[Functions]
  [fluence]
    type = ParsedFunction
    expression = 'if(y<0.5,0,
             if(y<2.5,4.6e25,
             if(y<3.5,9.6e26,1.2e27)))'
  []
  [temperature]
    type = ParsedFunction
    expression = 'if(y<0.5,294.15,
             if(y<1.5,323.15,
             if(y<2.5,423.15,
             if(y<3.5,323.15,423.15))))'
  []
  [disp_x]
    type = ParsedFunction
    expression = 't/300'
  []
[]

[AuxVariables]
  [temperature]
  []
[]

[ICs]
  [temperature]
    type = FunctionIC
    variable = temperature
    function = temperature
  []
[]

[AuxKernels]
  [temperature]
    type = FunctionAux
    function = temperature
    variable = temperature
  []
[]

[Physics/SolidMechanics/QuasiStatic]
  [all]
    add_variables = true
    strain = SMALL
    incremental = true
    generate_output = 'stress_xx strain_xx plastic_strain_xx strain_yy strain_zz stress_yy stress_zz'
  []
[]

[BCs]
  [fixed_x]
    type = DirichletBC
    value = 0
    boundary = left
    variable = disp_x
  []
  [fixed_y]
    type = DirichletBC
    value = 0
    boundary = bottom
    variable = disp_y
  []
  [fixed_z]
    type = DirichletBC
    value = 0
    boundary = back
    variable = disp_z
  []
  [pull]
    type = FunctionDirichletBC
    function = disp_x
    boundary = right
    variable = disp_x
  []
[]

[Materials]
  [fluence]
    type = GenericFunctionMaterial
    prop_names = fluence
    prop_values = fluence
  []
  [elasticity_tensor]
    type = Al6061ElasticityTensor
    fast_neutron_fluence = fluence
    temperature = temperature
  []
  [stress]
    type = ComputeMultipleInelasticStress
    inelastic_models = plasticity
  []
  [plasticity]
    type = Al6061PlasticityStressUpdate
    property = fluence
    temperature = temperature
    relative_tolerance = 1e-10
  []
[]

[Executioner]
  type = Transient
  solve_type = Newton
  line_search = 'none'
  nl_rel_tol = 1e-12
  end_time = 20
  dt = 1
[]

[Postprocessors]
  [xx_strain]
    type = ElementAverageValue
    variable = strain_xx
    execute_on = 'initial timestep_end'
  []
  [yy_strain]
    type = ElementAverageValue
    variable = strain_yy
    execute_on = 'initial timestep_end'
  []
  [zz_strain]
    type = ElementAverageValue
    variable = strain_zz
    execute_on = 'initial timestep_end'
  []
  [stress_0]
    type = ElementalVariableValue
    elementid = 0
    variable = stress_xx
  []
  [stress_1]
    type = ElementalVariableValue
    elementid = 1
    variable = stress_xx
  []
  [stress_2]
    type = ElementalVariableValue
    elementid = 2
    variable = stress_xx
  []
  [stress_3]
    type = ElementalVariableValue
    elementid = 3
    variable = stress_xx
  []
  [stress_4]
    type = ElementalVariableValue
    elementid = 4
    variable = stress_xx
  []
[]

[Outputs]
  csv = true
[]

(test/tests/solid_mechanics/Al6061_plasticity/plasticity.i)

#
# This test checks plasticity based on yield stress and hardening data
# recommended by the USHPRR program.
#
# The mesh involves separate 3D elements, one for each hardening curve
# to be tested.
#
[GlobalParams]
  displacements = 'disp_x disp_y disp_z'
[]

[Mesh]
  use_displaced_mesh = false
  [b0]
    type = GeneratedMeshGenerator
    dim = 3
    nx = 1
    ny = 1
    nz = 1
    xmin = 0
    xmax = 1
    ymin = 0
    ymax = 0.25
    zmin = 0
    zmax = 0.25
  []
  [all]
    type = CombinerGenerator
    inputs = b0
    positions = '0 0 0
                 0 1 0
                 0 2 0
                 0 3 0
                 0 4 0'
  []
  [final]
    type = SubdomainPerElementGenerator
    input = all
    subdomain_ids = '0 1 2 3 4'
  []
[]

[Functions]
  [fluence]
    type = ParsedFunction
    expression = 'if(y<0.5,0,
             if(y<2.5,4.6e25,
             if(y<3.5,9.6e26,1.2e27)))'
  []
  [temperature]
    type = ParsedFunction
    expression = 'if(y<0.5,294.15,
             if(y<1.5,323.15,
             if(y<2.5,423.15,
             if(y<3.5,323.15,423.15))))'
  []
  [disp_x]
    type = ParsedFunction
    expression = 't/300'
  []
[]

[AuxVariables]
  [temperature]
  []
[]

[ICs]
  [temperature]
    type = FunctionIC
    variable = temperature
    function = temperature
  []
[]

[AuxKernels]
  [temperature]
    type = FunctionAux
    function = temperature
    variable = temperature
  []
[]

[Physics/SolidMechanics/QuasiStatic]
  [all]
    add_variables = true
    strain = SMALL
    incremental = true
    generate_output = 'stress_xx strain_xx plastic_strain_xx strain_yy strain_zz stress_yy stress_zz'
  []
[]

[BCs]
  [fixed_x]
    type = DirichletBC
    value = 0
    boundary = left
    variable = disp_x
  []
  [fixed_y]
    type = DirichletBC
    value = 0
    boundary = bottom
    variable = disp_y
  []
  [fixed_z]
    type = DirichletBC
    value = 0
    boundary = back
    variable = disp_z
  []
  [pull]
    type = FunctionDirichletBC
    function = disp_x
    boundary = right
    variable = disp_x
  []
[]

[Materials]
  [fluence]
    type = GenericFunctionMaterial
    prop_names = fluence
    prop_values = fluence
  []
  [elasticity_tensor]
    type = Al6061ElasticityTensor
    fast_neutron_fluence = fluence
    temperature = temperature
  []
  [stress]
    type = ComputeMultipleInelasticStress
    inelastic_models = plasticity
  []
  [plasticity]
    type = Al6061PlasticityStressUpdate
    property = fluence
    temperature = temperature
    relative_tolerance = 1e-10
  []
[]

[Executioner]
  type = Transient
  solve_type = Newton
  line_search = 'none'
  nl_rel_tol = 1e-12
  end_time = 20
  dt = 1
[]

[Postprocessors]
  [xx_strain]
    type = ElementAverageValue
    variable = strain_xx
    execute_on = 'initial timestep_end'
  []
  [yy_strain]
    type = ElementAverageValue
    variable = strain_yy
    execute_on = 'initial timestep_end'
  []
  [zz_strain]
    type = ElementAverageValue
    variable = strain_zz
    execute_on = 'initial timestep_end'
  []
  [stress_0]
    type = ElementalVariableValue
    elementid = 0
    variable = stress_xx
  []
  [stress_1]
    type = ElementalVariableValue
    elementid = 1
    variable = stress_xx
  []
  [stress_2]
    type = ElementalVariableValue
    elementid = 2
    variable = stress_xx
  []
  [stress_3]
    type = ElementalVariableValue
    elementid = 3
    variable = stress_xx
  []
  [stress_4]
    type = ElementalVariableValue
    elementid = 4
    variable = stress_xx
  []
[]

[Outputs]
  csv = true
[]

(test/tests/solid_mechanics/Al6061_plasticity/interpolated.i)

#
# This test checks plasticity based on yield stress and hardening data
# recommended by the USHPRR program.
#
# The test pulls on an element while increasing temperature and fluence.
# This causes the yield stress and hardening curves to be interpolated
# between datasets.
#
# The usefulness of this test is in showing that the interpolation
# works and that the resulting response is robust.  The computed
# stresses have no inherent usefulness.
#
[GlobalParams]
  displacements = 'disp_x disp_y disp_z'
[]

[Mesh]
  use_displaced_mesh = false
  [b0]
    type = GeneratedMeshGenerator
    dim = 3
    nx = 1
    ny = 1
    nz = 1
    xmin = 0
    xmax = 1
    ymin = 0
    ymax = 0.25
    zmin = 0
    zmax = 0.25
  []
[]

[Functions]
  [fluence]
    type = ParsedFunction
    expression = '1.2e27*t/105'
  []
  [temperature]
    type = ParsedFunction
    expression = '294.15 + (423.15-294.15)*t/105'
  []
  [disp_x]
    type = ParsedFunction
    expression = 't/300'
  []
[]

[AuxVariables]
  [temperature]
  []
[]

[ICs]
  [temperature]
    type = FunctionIC
    variable = temperature
    function = temperature
  []
[]

[AuxKernels]
  [temperature]
    type = FunctionAux
    function = temperature
    variable = temperature
  []
[]

[Physics/SolidMechanics/QuasiStatic]
  [all]
    add_variables = true
    strain = SMALL
    incremental = true
    generate_output = 'stress_xx strain_xx plastic_strain_xx strain_yy strain_zz stress_yy stress_zz'
  []
[]

[BCs]
  [fixed_x]
    type = DirichletBC
    value = 0
    boundary = left
    variable = disp_x
  []
  [fixed_y]
    type = DirichletBC
    value = 0
    boundary = bottom
    variable = disp_y
  []
  [fixed_z]
    type = DirichletBC
    value = 0
    boundary = back
    variable = disp_z
  []
  [pull]
    type = FunctionDirichletBC
    function = disp_x
    boundary = right
    variable = disp_x
  []
[]

[Materials]
  [fluence]
    type = GenericFunctionMaterial
    prop_names = fluence
    prop_values = fluence
  []
  [elasticity_tensor]
    type = Al6061ElasticityTensor
    fast_neutron_fluence = fluence
    temperature = temperature
  []
  [stress]
    type = ComputeMultipleInelasticStress
    inelastic_models = plasticity
  []
  [plasticity]
    type = Al6061PlasticityStressUpdate
    property = fluence
    temperature = temperature
    relative_tolerance = 1e-10
  []
[]

[Executioner]
  type = Transient
  solve_type = Newton
  line_search = 'none'
  nl_rel_tol = 1e-12
  end_time = 105
  dt = 1
[]

[Postprocessors]
  [xx_strain]
    type = ElementAverageValue
    variable = strain_xx
    execute_on = 'initial timestep_end'
  []
  [yy_strain]
    type = ElementAverageValue
    variable = strain_yy
    execute_on = 'initial timestep_end'
  []
  [zz_strain]
    type = ElementAverageValue
    variable = strain_zz
    execute_on = 'initial timestep_end'
  []
  [stress_xx]
    type = ElementalVariableValue
    elementid = 0
    variable = stress_xx
  []
[]

[Outputs]
  csv = true
[]

Description
Example Input Syntax
Input Parameters
Input Files
References