Pressure Action System

Set up Pressure boundary conditions

Description

The Pressure Action, given in the input file as simply Pressure, is designed to simplify the input file when several variables have the same pressure boundary condition magnitude applied in the normal component. Hydrostatic stress is a good example of a use case for the Pressure Action.

Constructed MooseObjects

The Pressure Action is used to construct the Pressure boundary condition in all directions.

Table 1: Correspondence Among Action Functionality and MooseObjects

Functionality	Replaced Classes	Associated Parameters
A pressure traction force	Pressure BC	`displacements` : a string of the displacement variables to which the Pressure BC should be applied

The Pressure Action only applies the pressure traction in the same component direction as the displacements variables are listed. That is, for the argument displacements = 'disp_x disp_y disp_z', the Pressure Action will create three separate Pressure BCs:

For the variable disp_x, the parameter setting component = 0 is used
For the variable disp_y, the parameter setting component = 1 is used
For the variable disp_z, the parameter setting component = 2 is used

note:Displacement Variable-Component Relationship is Relative

Note that the location of each of the variables in the displacements string determines the value of the corresponding component.

As in the Pressure boundary condition, the magnitude of the Pressure boundary condition can be specified as either a scalar (use the input parameter factor) or a function parameter.

The use of the automatic differentiation boundary condition ADPressure can optionally be selected via the use_automatic_differentiation input parameter.

Example Input Syntax

[./Pressure]
  [./Side1]
    boundary = 1
    function = rampConstant
    displacements = 'disp_x disp_y disp_z'
  [../]
  [./Side2]
    boundary = 2
    function = zeroRamp
    displacements = 'disp_x disp_y disp_z'
    factor = 2.0
  [../]
  [./Side3]
    boundary = 3
    function = rampUnramp
    displacements = 'disp_x disp_y disp_z'
  [../]
[../]

/opt/civet/build_0/moose/modules/tensor_mechanics/test/tests/pressure/pressure_test.i

#
# Pressure Test
#
# This test is designed to compute pressure loads on three faces of a unit cube.
#
# The mesh is composed of one block with a single element.  Symmetry bcs are
# applied to the faces opposite the pressures.  Poisson's ratio is zero,
# which makes it trivial to check displacements.
#

[Mesh]
  type = FileMesh
  file = pressure_test.e
  displacements = 'disp_x disp_y disp_z'
[]

[Functions]
  [./rampConstant]
    type = PiecewiseLinear
    x = '0. 1. 2.'
    y = '0. 1. 1.'
    scale_factor = 1.0
  [../]
  [./zeroRamp]
    type = PiecewiseLinear
    x = '0. 1. 2.'
    y = '0. 0. 1.'
    scale_factor = 1.0
  [../]
  [./rampUnramp]
    type = PiecewiseLinear
    x = '0. 1. 2.'
    y = '0. 1. 0.'
    scale_factor = 10.0
  [../]
[]

[Variables]
  [./disp_x]
    order = FIRST
    family = LAGRANGE
  [../]
  [./disp_y]
    order = FIRST
    family = LAGRANGE
  [../]
  [./disp_z]
    order = FIRST
    family = LAGRANGE
  [../]
[]

[Kernels]
  [./TensorMechanics]
    displacements = 'disp_x disp_y disp_z'
  [../]
[]

[BCs]
  [./no_x]
    type = DirichletBC
    variable = disp_x
    boundary = 4
    value = 0.0
  [../]
  [./no_y]
    type = DirichletBC
    variable = disp_y
    boundary = 5
    value = 0.0
  [../]
  [./no_z]
    type = DirichletBC
    variable = disp_z
    boundary = 6
    value = 0.0
  [../]
  [./Pressure]
    [./Side1]
      boundary = 1
      function = rampConstant
      displacements = 'disp_x disp_y disp_z'
    [../]
    [./Side2]
      boundary = 2
      function = zeroRamp
      displacements = 'disp_x disp_y disp_z'
      factor = 2.0
    [../]
    [./Side3]
      boundary = 3
      function = rampUnramp
      displacements = 'disp_x disp_y disp_z'
    [../]
  [../]
[]

[Materials]
  [./Elasticity_tensor]
    type = ComputeElasticityTensor
    block = 1
    fill_method = symmetric_isotropic
    C_ijkl = '0 0.5e6'
  [../]
  [./strain]
    type = ComputeSmallStrain
    displacements = 'disp_x disp_y disp_z'
    block = 1
  [../]
  [./stress]
    type = ComputeLinearElasticStress
    block = 1
  [../]
[]

[Executioner]
  type = Transient
  solve_type = PJFNK
  nl_abs_tol = 1e-10
  l_max_its = 20
  start_time = 0.0
  dt = 1.0
  num_steps = 2
  end_time = 2.0
[]

[Outputs]
  [./out]
    type = Exodus
    elemental_as_nodal = true
  [../]
[]

Input Parameters

boundaryThe list of boundary IDs from the mesh where the pressure will be applied
C++ Type:std::vector
Options:
Description:The list of boundary IDs from the mesh where the pressure will be applied

Required Parameters

active__all__ If specified only the blocks named will be visited and made active
Default:__all__
C++ Type:std::vector
Options:
Description:If specified only the blocks named will be visited and made active
alpha0alpha parameter for HHT time integration
Default:0
C++ Type:double
Options:
Description:alpha parameter for HHT time integration
disp_xThe x displacement
C++ Type:VariableName
Options:
Description:The x displacement
disp_yThe y displacement
C++ Type:VariableName
Options:
Description:The y displacement
disp_zThe z displacement
C++ Type:VariableName
Options:
Description:The z displacement
displacementsThe displacements appropriate for the simulation geometry and coordinate system
C++ Type:std::vector
Options:
Description:The displacements appropriate for the simulation geometry and coordinate system
factor1The factor to use in computing the pressure
Default:1
C++ Type:double
Options:
Description:The factor to use in computing the pressure
functionThe function that describes the pressure
C++ Type:FunctionName
Options:
Description:The function that describes the pressure
inactiveIf specified blocks matching these identifiers will be skipped.
C++ Type:std::vector
Options:
Description:If specified blocks matching these identifiers will be skipped.
save_in_disp_xThe save_in variables for x displacement
C++ Type:std::vector
Options:
Description:The save_in variables for x displacement
save_in_disp_yThe save_in variables for y displacement
C++ Type:std::vector
Options:
Description:The save_in variables for y displacement
save_in_disp_zThe save_in variables for z displacement
C++ Type:std::vector
Options:
Description:The save_in variables for z displacement
use_automatic_differentiationFalseFlag to use automatic differentiation (AD) objects when possible
Default:False
C++ Type:bool
Options:
Description:Flag to use automatic differentiation (AD) objects when possible

Optional Parameters

Associated Actions

Available Actions

Tensor Mechanics App
PressureActionSet up Pressure boundary conditions

Description
Constructed MooseObjects
Example Input Syntax
Input Parameters
Associated Actions
Available Actions

Application Usage

Physics and Syntax

Application Development

Framework Development

MOOSEDocs

Infrastructure

Questions

Information and Tools