AD Exponential Softening

Softening model with an exponential softening response upon cracking. This class is intended to be used with ADComputeSmearedCrackingStress and relies on automatic differentiation.

Description

The material ADExponentialSoftening computes the reduced stress and stiffness in the direction of a crack according to a exponential function. The computed cracked stiffness ratio softens the tensile response of the material once the principle stress exceeds the cracking stress threshold of the material. It uses automatic differentiation for the computation of its Jacobian.

For more details, refer to the non-AD version ExponentialSoftening.

Example Input File

[./exponential_softening]
  type = ADExponentialSoftening
[../]

ADExponentialSoftening must be run in conjunction with the ad fixed smeared cracking material model as shown below:

[./cracking_stress]
  type = ADComputeSmearedCrackingStress
  shear_retention_factor = 0.1
  cracking_stress = 3.e9
  cracked_elasticity_type = FULL
  softening_models = exponential_softening
[../]

Input Parameters

alpha-1Initial slope of the exponential softening curve at crack initiation. If not specified, it is equal to the negative of the Young's modulus.
Default:-1
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Initial slope of the exponential softening curve at crack initiation. If not specified, it is equal to the negative of the Young's modulus.
beta1Multiplier applied to alpha to control the exponential softening behavior.
Default:1
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Multiplier applied to alpha to control the exponential softening behavior.
blockThe list of blocks (ids or names) that this object will be applied
C++ Type:std::vector<SubdomainName>
Unit:(no unit assumed)
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>
Unit:(no unit assumed)
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
Unit:(no unit assumed)
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.
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.
residual_stress0The fraction of the cracking stress allowed to be maintained following a crack.
Default:0
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:The fraction of the cracking stress allowed to be maintained following a crack.
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
Unit:(no unit assumed)
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.

Optional Parameters

control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector<std::string>
Unit:(no unit assumed)
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
Unit:(no unit assumed)
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
Unit:(no unit assumed)
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
Unit:(no unit assumed)
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
Unit:(no unit assumed)
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>
Unit:(no unit assumed)
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>
Unit:(no unit assumed)
Controllable:No
Description:Vector of output names where you would like to restrict the output of variables(s) associated with this object

Outputs Parameters

References

No citations exist within this document.

(moose/modules/solid_mechanics/test/tests/ad_smeared_cracking/cracking_rotation.i)

# This test is to ensure that the smeared cracking model correctly handles finite
# rotation of cracked elements.

# This consists of a single element that is first  subjected to tensile loading
# in the y-direction via a prescribed displacement. This loading is sufficiently
# high to crack the material in that direction, but not completely unload. The
# prescribed displacement is then reversed so that the element is returned to its
# original configuration.

# In the next phase of the analysis, this element is then rotated 90 degrees by
# prescribing the displacement of the bottom of the element. The prescribed
# displacement BC used to crack the element in the first phase is deactivated.

# Once the element is fully rotated, a new BC is activated on what was originally
# the top surface (but is now the surface on the right hand side) to pull in
# the x-direction.

# If everything is working correctly, the model should re-load on the original
# crack (which should be rotated along with the elemnent) up to the peak stress
# in the first phase of the analysis, and then continue the unloading process
# as the crack strains continue to increase. Throughout this analysis, there should
# only be a single crack, as manifested in the crack_flags variables.

[Mesh]
  file = cracking_test.e
[]

[GlobalParams]
  displacements = 'disp_x disp_y disp_z'
[]

[Physics/SolidMechanics/QuasiStatic]
  [./all]
    strain = FINITE
    add_variables = true
    generate_output = 'stress_xx stress_yy stress_zz stress_xy stress_yz stress_zx'
    use_automatic_differentiation = true
  [../]
[]

[AuxVariables]
  [./crack_flags1]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./crack_flags2]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./crack_flags3]
    order = CONSTANT
    family = MONOMIAL
  [../]
[]

[AuxKernels]
  [./crack_flags1]
    type = ADMaterialRealVectorValueAux
    property = crack_flags
    variable = crack_flags1
    component = 0
  [../]
  [./crack_flags2]
    type = ADMaterialRealVectorValueAux
    property = crack_flags
    variable = crack_flags2
    component = 1
  [../]
  [./crack_flags3]
    type = ADMaterialRealVectorValueAux
    property = crack_flags
    variable = crack_flags3
    component = 2
  [../]
[]

[BCs]
  [./x_pin]
    type = ADDirichletBC
    variable = disp_x
    boundary = '15 16'
    value = 0.0
  [../]
  [./y_pin]
    type = ADDirichletBC
    variable = disp_y
    boundary = '15 16'
    value = 0.0
  [../]
  [./z_all]
    type = ADDirichletBC
    variable = disp_z
    boundary = '11 12 13 14 15 16 17 18'
    value = 0.0
  [../]
  [./x_lb]
    type = ADFunctionDirichletBC
    variable = disp_x
    boundary = '11 12'
    function = 'if(t<10,0,if(t>=100,1,1-cos((t-10)*pi/180)))'
  [../]
  [./y_lb]
    type = ADFunctionDirichletBC
    variable = disp_y
    boundary = '11 12'
    function = 'if(t<10,0,if(t>=100,1,sin((t-10)*pi/180)))'
  [../]
  [./x_lt]
    type = ADFunctionDirichletBC
    variable = disp_x
    boundary = '13 14'
    function = '2+(t-100)*0.01'
  [../]
  [./x_rt]
    type = ADFunctionDirichletBC
    variable = disp_x
    boundary = '17 18'
    function = '1+(t-100)*0.01'
  [../]
  [./top_pull]
    type = ADFunctionDirichletBC
    variable = disp_y
    boundary = '13 14 17 18'
    function = 'if(t<5,t*0.01,0.05-(t-5)*0.01)'
  [../]
[]

[Materials]
  [./elasticity_tensor]
    type = ADComputeIsotropicElasticityTensor
    youngs_modulus = 100.e9
    poissons_ratio = 0.
  [../]
  [./cracking_stress]
    type = ADComputeSmearedCrackingStress
    shear_retention_factor = 0.1
    cracking_stress = 3.e9
    cracked_elasticity_type = FULL
    softening_models = exponential_softening
  [../]
  [./exponential_softening]
    type = ADExponentialSoftening
  [../]
[]

[Executioner]
  type = Transient
  solve_type = Newton
  petsc_options_iname = '-ksp_gmres_restart -pc_type'
  petsc_options_value = '101                lu'
  line_search = 'none'

  l_max_its = 100
  l_tol = 1e-5

  nl_max_its = 100
  nl_abs_tol = 1e-5
  nl_rel_tol = 1e-12

  start_time = 0
  end_time = 110
  dt = 1
[]

[Controls]
  [./p1]
    type = TimePeriod
    start_time = 0.0
    end_time = 10.0
    disable_objects = 'BCs/x_lt BCs/x_rt'
    enable_objects = 'BCs/top_pull'
    reverse_on_false = false
    execute_on = 'initial timestep_begin'
  [../]
  [./p2]
    type = TimePeriod
    start_time = 10.0
    end_time = 101.0
    disable_objects = 'BCs/x_lt BCs/x_rt BCs/top_pull'
    reverse_on_false = false
    execute_on = 'initial timestep_begin'
  [../]
  [./p3]
    type = TimePeriod
    start_time = 101.0
    end_time = 110.0
    enable_objects = 'BCs/x_lt BCs/x_rt'
    disable_objects = 'BCs/top_pull'
    reverse_on_false = false
    execute_on = 'initial timestep_begin'
  [../]
[]

[Outputs]
  exodus = true
[]

(moose/modules/solid_mechanics/test/tests/ad_smeared_cracking/cracking_rotation.i)

# This test is to ensure that the smeared cracking model correctly handles finite
# rotation of cracked elements.

# This consists of a single element that is first  subjected to tensile loading
# in the y-direction via a prescribed displacement. This loading is sufficiently
# high to crack the material in that direction, but not completely unload. The
# prescribed displacement is then reversed so that the element is returned to its
# original configuration.

# In the next phase of the analysis, this element is then rotated 90 degrees by
# prescribing the displacement of the bottom of the element. The prescribed
# displacement BC used to crack the element in the first phase is deactivated.

# Once the element is fully rotated, a new BC is activated on what was originally
# the top surface (but is now the surface on the right hand side) to pull in
# the x-direction.

# If everything is working correctly, the model should re-load on the original
# crack (which should be rotated along with the elemnent) up to the peak stress
# in the first phase of the analysis, and then continue the unloading process
# as the crack strains continue to increase. Throughout this analysis, there should
# only be a single crack, as manifested in the crack_flags variables.

[Mesh]
  file = cracking_test.e
[]

[GlobalParams]
  displacements = 'disp_x disp_y disp_z'
[]

[Physics/SolidMechanics/QuasiStatic]
  [./all]
    strain = FINITE
    add_variables = true
    generate_output = 'stress_xx stress_yy stress_zz stress_xy stress_yz stress_zx'
    use_automatic_differentiation = true
  [../]
[]

[AuxVariables]
  [./crack_flags1]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./crack_flags2]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./crack_flags3]
    order = CONSTANT
    family = MONOMIAL
  [../]
[]

[AuxKernels]
  [./crack_flags1]
    type = ADMaterialRealVectorValueAux
    property = crack_flags
    variable = crack_flags1
    component = 0
  [../]
  [./crack_flags2]
    type = ADMaterialRealVectorValueAux
    property = crack_flags
    variable = crack_flags2
    component = 1
  [../]
  [./crack_flags3]
    type = ADMaterialRealVectorValueAux
    property = crack_flags
    variable = crack_flags3
    component = 2
  [../]
[]

[BCs]
  [./x_pin]
    type = ADDirichletBC
    variable = disp_x
    boundary = '15 16'
    value = 0.0
  [../]
  [./y_pin]
    type = ADDirichletBC
    variable = disp_y
    boundary = '15 16'
    value = 0.0
  [../]
  [./z_all]
    type = ADDirichletBC
    variable = disp_z
    boundary = '11 12 13 14 15 16 17 18'
    value = 0.0
  [../]
  [./x_lb]
    type = ADFunctionDirichletBC
    variable = disp_x
    boundary = '11 12'
    function = 'if(t<10,0,if(t>=100,1,1-cos((t-10)*pi/180)))'
  [../]
  [./y_lb]
    type = ADFunctionDirichletBC
    variable = disp_y
    boundary = '11 12'
    function = 'if(t<10,0,if(t>=100,1,sin((t-10)*pi/180)))'
  [../]
  [./x_lt]
    type = ADFunctionDirichletBC
    variable = disp_x
    boundary = '13 14'
    function = '2+(t-100)*0.01'
  [../]
  [./x_rt]
    type = ADFunctionDirichletBC
    variable = disp_x
    boundary = '17 18'
    function = '1+(t-100)*0.01'
  [../]
  [./top_pull]
    type = ADFunctionDirichletBC
    variable = disp_y
    boundary = '13 14 17 18'
    function = 'if(t<5,t*0.01,0.05-(t-5)*0.01)'
  [../]
[]

[Materials]
  [./elasticity_tensor]
    type = ADComputeIsotropicElasticityTensor
    youngs_modulus = 100.e9
    poissons_ratio = 0.
  [../]
  [./cracking_stress]
    type = ADComputeSmearedCrackingStress
    shear_retention_factor = 0.1
    cracking_stress = 3.e9
    cracked_elasticity_type = FULL
    softening_models = exponential_softening
  [../]
  [./exponential_softening]
    type = ADExponentialSoftening
  [../]
[]

[Executioner]
  type = Transient
  solve_type = Newton
  petsc_options_iname = '-ksp_gmres_restart -pc_type'
  petsc_options_value = '101                lu'
  line_search = 'none'

  l_max_its = 100
  l_tol = 1e-5

  nl_max_its = 100
  nl_abs_tol = 1e-5
  nl_rel_tol = 1e-12

  start_time = 0
  end_time = 110
  dt = 1
[]

[Controls]
  [./p1]
    type = TimePeriod
    start_time = 0.0
    end_time = 10.0
    disable_objects = 'BCs/x_lt BCs/x_rt'
    enable_objects = 'BCs/top_pull'
    reverse_on_false = false
    execute_on = 'initial timestep_begin'
  [../]
  [./p2]
    type = TimePeriod
    start_time = 10.0
    end_time = 101.0
    disable_objects = 'BCs/x_lt BCs/x_rt BCs/top_pull'
    reverse_on_false = false
    execute_on = 'initial timestep_begin'
  [../]
  [./p3]
    type = TimePeriod
    start_time = 101.0
    end_time = 110.0
    enable_objects = 'BCs/x_lt BCs/x_rt'
    disable_objects = 'BCs/top_pull'
    reverse_on_false = false
    execute_on = 'initial timestep_begin'
  [../]
[]

[Outputs]
  exodus = true
[]