HyperElasticIsotropic

Hyperelastic isotropic constitutive model.

HyperElasticIsotropic implements a small set of hyperelastic isotropic relations. Given deformation gradient $var element = document.getElementById("moose-equation-28786d13-cf77-4f89-8743-57ab71dfedda");katex.render("\\mathbf{F}", element, {displayMode:false,throwOnError:false});$ , shear modulus $var element = document.getElementById("moose-equation-47cd9fe9-abbd-41d8-9a06-8f81571bd0d0");katex.render("\\mu", element, {displayMode:false,throwOnError:false});$ , and bulk modulus $var element = document.getElementById("moose-equation-2bf8d16e-ed87-4797-bef3-ed957c6a9055");katex.render("K", element, {displayMode:false,throwOnError:false});$ , it computes stress and an optional consistent tangent operator.

Overview

Inputs: - "F": deformation gradient tensor $var element = document.getElementById("moose-equation-99c2cdb1-37e3-4bec-8def-f9d71ceb7eb2");katex.render("\\mathbf{F}", element, {displayMode:false,throwOnError:false});$ (rank-2). - "mu": shear modulus $var element = document.getElementById("moose-equation-396de6f7-1577-427f-9b20-1bc66edd8698");katex.render("\\mu", element, {displayMode:false,throwOnError:false});$ . - "K": bulk modulus $var element = document.getElementById("moose-equation-f1027762-7620-494f-a55d-5ca4f2023a60");katex.render("K", element, {displayMode:false,throwOnError:false});$ .

Outputs: - Primary output: first Piola-Kirchhoff-like stress $var element = document.getElementById("moose-equation-983ef4ec-2608-4a67-b8e4-1629c6f6f2ed");katex.render("\\mathbf{P}", element, {displayMode:false,throwOnError:false});$ computed from $var element = document.getElementById("moose-equation-07f08602-42ab-4aee-a90b-9892eea9ce80");katex.render("\\mathbf{F}", element, {displayMode:false,throwOnError:false});$ and isotropic 4th-order moduli. - Optional: "tangent_operator" to receive the 4th-order stiffness.

Example Input File Syntax

[TensorComputes<<<{"href": "../../syntax/TensorComputes/index.html"}>>>]
  [Solve<<<{"href": "../../syntax/TensorComputes/Solve/index.html"}>>>]
    [hyper_elasticity]
      type = HyperElasticIsotropic<<<{"description": "Hyperelastic isotropic constitutive model.", "href": "HyperElasticIsotropic.html"}>>>
      buffer<<<{"description": "The buffer this compute is writing to"}>>> = stress
      F<<<{"description": "Deformation gradient tensor"}>>> = Fnew
      K<<<{"description": "Deformation gradient tensor"}>>> = K
      mu<<<{"description": "Deformation gradient tensor"}>>> = mu
    []
  []
[]

Input Parameters

FDeformation gradient tensor
C++ Type:std::string
Controllable:No
Description:Deformation gradient tensor
KDeformation gradient tensor
C++ Type:std::string
Controllable:No
Description:Deformation gradient tensor
bufferThe buffer this compute is writing to
C++ Type:std::string
Controllable:No
Description:The buffer this compute is writing to
muDeformation gradient tensor
C++ Type:std::string
Controllable:No
Description:Deformation gradient tensor

Required Parameters

tangent_operatordstressdstrainStiffness tensor
Default:dstressdstrain
C++ Type:std::string
Controllable:No
Description:Stiffness tensor

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:No
Description:Set the enabled status of the MooseObject.

Advanced Parameters

(test/tests/mechanics/mech.i)

[Domain]
  dim = 2
  nx = 32
  ny = 32
  xmax = ${fparse 2*pi}
  ymax = ${fparse 2*pi}
  zmax = ${fparse 2*pi}
  mesh_mode = DUMMY
  parallel_mode = FFT_AUTO
[]

[TensorComputes]
  [Initialize]
    [phase]
      type = ParsedCompute
      expression = '(cos(x)/2+0.5)^1*(cos(y)/2+0.5)^1*(cos(z)/2+0.5)^1'
      extra_symbols = true
      buffer = phase
    []
    [K]
      type = ParsedCompute
      buffer = K
      expression = '(1-phase)*Ka + phase*Kb'
      inputs = phase
      constant_names = 'Ka Kb'
      constant_expressions = '1 10'
    []
    [mu]
      type = ParsedCompute
      buffer = mu
      expression = '(1-phase)*mua + phase*mub'
      inputs = phase
      constant_names = 'mua mub'
      constant_expressions = '0.5 5'
    []
    [Finit]
      type = RankTwoIdentity
      buffer = F
    []
  []

  [Solve]
    [hyper_elasticity]
      type = HyperElasticIsotropic
      buffer = stress
      F = Fnew
      K = K
      mu = mu
    []

    [root]
      [applied_strain]
        type = MacroscopicShearTensor
        buffer = applied_strain
      []
      [mech]
        type = FFTMechanics
        buffer = Fnew
        F = F
        K = K
        mu = mu
        l_max_its = 40
        l_tol = 1e-5
        nl_rel_tol = 2e-4
        nl_abs_tol = 2e-3
        constitutive_model = hyper_elasticity
        stress = stress
        applied_macroscopic_strain = applied_strain
      []
    []
  []

  [Postprocess]
    [displacements]
      type = ComputeDisplacements
      buffer = disp
      F = F
    []
    [vonmises]
      type = ComputeVonMisesStress
      buffer = sV
    []
  []
[]

[TensorSolver]
  # no variables are integrated by this solver (FFTMechanics performs a steady state mechanics solve)
  type = ForwardEulerSolver
  root_compute = root
  # deformation tensor is just forwarded Fnew -> F
  forward_buffer = F
  forward_buffer_new = Fnew
  substeps = 3
[]

[TensorOutputs]
  active = 'parallel' # serial

  [serial]
    type = XDMFTensorOutput
    file_base = mech_serial
    buffer = 'disp sV F phase'
    output_mode = 'OVERSIZED_NODAL CELL CELL NODE'
    enable_hdf5 = true
    execute_on = 'TIMESTEP_END'
  []
  [parallel]
    # NODE and OVERSIZED_NODAL are not yet available in parallel output
    type = XDMFTensorOutput
    file_base = mech_parallel
    buffer = 'sV F phase'
    output_mode = 'CELL CELL CELL'
    enable_hdf5 = true
    execute_on = 'TIMESTEP_END'
  []
[]


[Executioner]
  type = Transient
  num_steps = 3
  dt = 0.02
[]

(examples/degeus_mechanics/mech.i)

[Domain]
  dim = 3
  nx = 32
  ny = 32
  nz = 32
  xmax = ${fparse 2*pi}
  ymax = ${fparse 2*pi}
  zmax = ${fparse 2*pi}
  mesh_mode = DUMMY
[]

[TensorComputes]
  [Initialize]
    [Finit]
      type = RankTwoIdentity
      buffer = F
    []

    [phase]
      type = PhaseMechanicsTest
      buffer = phase
    []
    [K]
      type = ParsedCompute
      buffer = K
      expression = '(1-phase)*Ka + phase*Kb'
      inputs = phase
      constant_names = 'Ka Kb'
      constant_expressions = '0.833 8.33'
    []
    [mu]
      type = ParsedCompute
      buffer = mu
      expression = '(1-phase)*mua + phase*mub'
      inputs = phase
      constant_names = 'mua mub'
      constant_expressions = '0.386 3.86'
    []
  []

  [Solve]
    [hyper_elasticity]
      type = HyperElasticIsotropic
      buffer = stress
      F = Fnew
      K = K
      mu = mu
    []

    [root]
      [applied_strain]
        type = MacroscopicShearTensor
        buffer = applied_strain
      []
      [mech]
        type = FFTMechanics
        buffer = Fnew
        F = F
        K = K
        mu = mu
        l_tol = 1e-2
        nl_rel_tol = 2e-2
        nl_abs_tol = 2e-2
        constitutive_model = hyper_elasticity
        stress = stress
        applied_macroscopic_strain = applied_strain
      []
    []
  []

  [Postprocess]
    [displacements]
      type = ComputeDisplacements
      buffer = disp
      F = F
    []
    [vonmises]
      type = ComputeVonMisesStress
      buffer = sV
    []
  []
[]

[TensorSolver]
  # no variables are integrated by this solver (FFTMechanics performs a steady state mechanics solve)
  type = ForwardEulerSolver
  root_compute = root
  # deformation tensor is just forwarded Fnew -> F
  forward_buffer = F
  forward_buffer_new = Fnew
  substeps = 10
[]

[TensorOutputs]
  [deformation_tensor]
    type = XDMFTensorOutput
    buffer = 'disp sV F'
    output_mode = 'OVERSIZED_NODAL CELL CELL'
    enable_hdf5 = true
  []
[]

[Executioner]
  type = Transient
  num_steps = 100
  dt = 0.01
[]

(test/tests/mechanics/mech.i)

[Domain]
  dim = 2
  nx = 32
  ny = 32
  xmax = ${fparse 2*pi}
  ymax = ${fparse 2*pi}
  zmax = ${fparse 2*pi}
  mesh_mode = DUMMY
  parallel_mode = FFT_AUTO
[]

[TensorComputes]
  [Initialize]
    [phase]
      type = ParsedCompute
      expression = '(cos(x)/2+0.5)^1*(cos(y)/2+0.5)^1*(cos(z)/2+0.5)^1'
      extra_symbols = true
      buffer = phase
    []
    [K]
      type = ParsedCompute
      buffer = K
      expression = '(1-phase)*Ka + phase*Kb'
      inputs = phase
      constant_names = 'Ka Kb'
      constant_expressions = '1 10'
    []
    [mu]
      type = ParsedCompute
      buffer = mu
      expression = '(1-phase)*mua + phase*mub'
      inputs = phase
      constant_names = 'mua mub'
      constant_expressions = '0.5 5'
    []
    [Finit]
      type = RankTwoIdentity
      buffer = F
    []
  []

  [Solve]
    [hyper_elasticity]
      type = HyperElasticIsotropic
      buffer = stress
      F = Fnew
      K = K
      mu = mu
    []

    [root]
      [applied_strain]
        type = MacroscopicShearTensor
        buffer = applied_strain
      []
      [mech]
        type = FFTMechanics
        buffer = Fnew
        F = F
        K = K
        mu = mu
        l_max_its = 40
        l_tol = 1e-5
        nl_rel_tol = 2e-4
        nl_abs_tol = 2e-3
        constitutive_model = hyper_elasticity
        stress = stress
        applied_macroscopic_strain = applied_strain
      []
    []
  []

  [Postprocess]
    [displacements]
      type = ComputeDisplacements
      buffer = disp
      F = F
    []
    [vonmises]
      type = ComputeVonMisesStress
      buffer = sV
    []
  []
[]

[TensorSolver]
  # no variables are integrated by this solver (FFTMechanics performs a steady state mechanics solve)
  type = ForwardEulerSolver
  root_compute = root
  # deformation tensor is just forwarded Fnew -> F
  forward_buffer = F
  forward_buffer_new = Fnew
  substeps = 3
[]

[TensorOutputs]
  active = 'parallel' # serial

  [serial]
    type = XDMFTensorOutput
    file_base = mech_serial
    buffer = 'disp sV F phase'
    output_mode = 'OVERSIZED_NODAL CELL CELL NODE'
    enable_hdf5 = true
    execute_on = 'TIMESTEP_END'
  []
  [parallel]
    # NODE and OVERSIZED_NODAL are not yet available in parallel output
    type = XDMFTensorOutput
    file_base = mech_parallel
    buffer = 'sV F phase'
    output_mode = 'CELL CELL CELL'
    enable_hdf5 = true
    execute_on = 'TIMESTEP_END'
  []
[]


[Executioner]
  type = Transient
  num_steps = 3
  dt = 0.02
[]

(test/tests/mechanics/mech3d.i)

[Domain]
  dim = 3
  nx = 16
  ny = 16
  nz = 16
  xmax = ${fparse 2*pi}
  ymax = ${fparse 2*pi}
  zmax = ${fparse 2*pi}
  mesh_mode = DUMMY
  parallel_mode = FFT_AUTO
[]

[TensorComputes]
  [Initialize]
    [phase]
      type = ParsedCompute
      expression = '(cos(x)/2+0.5)^1*(cos(y)/2+0.5)^1*(cos(z)/2+0.5)^1'
      extra_symbols = true
      buffer = phase
    []
    [K]
      type = ParsedCompute
      buffer = K
      expression = '(1-phase)*Ka + phase*Kb'
      inputs = phase
      constant_names = 'Ka Kb'
      constant_expressions = '1 10'
    []
    [mu]
      type = ParsedCompute
      buffer = mu
      expression = '(1-phase)*mua + phase*mub'
      inputs = phase
      constant_names = 'mua mub'
      constant_expressions = '0.5 5'
    []
    [Finit]
      type = RankTwoIdentity
      buffer = F
    []
  []

  [Solve]
    [hyper_elasticity]
      type = HyperElasticIsotropic
      buffer = stress
      F = Fnew
      K = K
      mu = mu
    []

    [root]
      [applied_strain]
        type = MacroscopicShearTensor
        buffer = applied_strain
      []
      [mech]
        type = FFTMechanics
        buffer = Fnew
        F = F
        K = K
        mu = mu
        l_tol = 1e-2
        nl_rel_tol = 2e-2
        nl_abs_tol = 2e-2
        constitutive_model = hyper_elasticity
        stress = stress
        applied_macroscopic_strain = applied_strain
      []
    []
  []

  [Postprocess]
    [displacements]
      type = ComputeDisplacements
      buffer = disp
      F = F
    []
    [vonmises]
      type = ComputeVonMisesStress
      buffer = sV
    []
  []
[]

[TensorSolver]
  # no variables are integrated by this solver (FFTMechanics performs a steady state mechanics solve)
  type = ForwardEulerSolver
  root_compute = root
  # deformation tensor is just forwarded Fnew -> F
  forward_buffer = F
  forward_buffer_new = Fnew
  substeps = 10
[]

[TensorOutputs]
  active = 'parallel' # serial

  [serial]
    type = XDMFTensorOutput
    file_base = mech3d_serial
    buffer = 'disp sV F phase'
    output_mode = 'OVERSIZED_NODAL CELL CELL NODE'
    enable_hdf5 = true
    execute_on = 'TIMESTEP_END'
  []
  [parallel]
    # NODE and OVERSIZED_NODAL are not yet available in parallel output
    type = XDMFTensorOutput
    file_base = mech3d_parallel
    buffer = 'sV F phase'
    output_mode = 'CELL CELL CELL'
    enable_hdf5 = true
    execute_on = 'TIMESTEP_END'
  []
[]

[Executioner]
  type = Transient
  num_steps = 3
  dt = 0.01
[]

Overview
Example Input File Syntax
Input Parameters
Input Files