ComputeFVDamperElasticity

Compute the deformations, forces and the stiffness matrix corresponding to a two-node nonlinear fluid viscous damper element.

Description

This material class computes the deformations, forces and stiffness matrix of a nonlinear fluid viscous damper element. Further information on the theory of the FVD element can be found at FV damper theory. For detailed information on the usage of the element with examples, refer to FV damper user.

Input Parameters

alphaVelocity exponent of the damper.
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Velocity exponent of the damper.
betaBeta parameter of Newmark algorithm.
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Beta parameter of Newmark algorithm.
cdDamping co-efficient.
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Damping co-efficient.
displacementsThe displacement variables appropriate for the simulation of geometry and coordinate system.
C++ Type:std::vector<VariableName>
Unit:(no unit assumed)
Controllable:No
Description:The displacement variables appropriate for the simulation of geometry and coordinate system.
gammaGamma parameter of Newmark algorithm.
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Gamma parameter of Newmark algorithm.
kAxial stiffness of the damper assembly.
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Axial stiffness of the damper assembly.
y_orientationOrientation of the y direction along with Ky is provided. This should be perpendicular to the axis of the Damper.
C++ Type:libMesh::VectorValue<double>
Unit:(no unit assumed)
Controllable:No
Description:Orientation of the y direction along with Ky is provided. This should be perpendicular to the axis of the Damper.

Required Parameters

abs_tol1e-06Absolute tolerance for error in adaptive algorithm
Default:1e-06
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Absolute tolerance for error in adaptive algorithm
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_onELEMENTWhen 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:ELEMENT
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.
rel_tol1e-06Relative tolerance for error in adaptive algorithm
Default:1e-06
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Relative tolerance for error in adaptive algorithm

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

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/materials/fv_damper/fv_damper_seismic_frame.i)
(test/tests/materials/fv_damper/fv_damper_axial.i)

References

No citations exist within this document.

(test/tests/materials/fv_damper/fv_damper_seismic_frame.i)

# Test for a protal frame installed with a diagonal Fluid Viscous Damper element
# The mesh for the geometry is imported from exodus file "fv_damper_frame.e"
# Beam-column elements are modeled as Euler-Bernouli beam elements

# Dimesnions
# columns
# E = 20e9, G = 100e9, A = 9e-2, Iy = Iz = 6.75e-4

# Beams
# E = 20e9, G = 100e9, A = 9e-2, Iy = Iz = 6.75e-1 (I_beam = 1000*I_Column to simulate rigid beam)

[Mesh]
  type = FileMesh
  file = 'fv_damper_frame.e'
[]

[Variables]
  [./disp_x]
    order = FIRST
    family = LAGRANGE
  [../]
  [./disp_y]
    order = FIRST
    family = LAGRANGE
  [../]
  [./disp_z]
    order = FIRST
    family = LAGRANGE
  [../]
  [./rot_x]
    order = FIRST
    family = LAGRANGE
    block = '1 2'
  [../]
  [./rot_y]
    order = FIRST
    family = LAGRANGE
    block = '1 2'
  [../]
  [./rot_z]
    order = FIRST
    family = LAGRANGE
    block = '1 2'
  [../]
[]

[AuxVariables]
  [./vel_x]
    order = FIRST
    family = LAGRANGE
  [../]
  [./vel_y]
    order = FIRST
    family = LAGRANGE
  [../]
  [./vel_z]
    order = FIRST
    family = LAGRANGE
  [../]
  [./accel_x]
    order = FIRST
    family = LAGRANGE
  [../]
  [./accel_y]
    order = FIRST
    family = LAGRANGE
  [../]
  [./accel_z]
    order = FIRST
    family = LAGRANGE
  [../]
  [./reaction_x]
    block = '1 2 3'
  [../]
  [./reaction_y]
    block = '1 2 3'
  [../]
  [./reaction_z]
    block = '1 2 3'
  [../]
  [./reaction_xx]
    block = '1 2 3'
  [../]
  [./reaction_yy]
    block = '1 2 3'
  [../]
  [./reaction_zz]
    block = '1 2 3'
  [../]
[]

[AuxKernels]
  [./accel_x]
    type = NewmarkAccelAux
    variable = accel_x
    displacement = disp_x
    velocity = vel_x
    beta = 0.25
    execute_on = timestep_end
  [../]
  [./vel_x]
    type = NewmarkVelAux
    variable = vel_x
    acceleration = accel_x
    gamma = 0.5
    execute_on = timestep_end
  [../]
  [./accel_y]
    type = NewmarkAccelAux
    variable = accel_y
    displacement = disp_y
    velocity = vel_y
    beta = 0.25
    execute_on = timestep_end
  [../]
  [./vel_y]
    type = NewmarkVelAux
    variable = vel_y
    acceleration = accel_y
    gamma = 0.5
    execute_on = timestep_end
  [../]
  [./accel_z]
    type = NewmarkAccelAux
    variable = accel_z
    displacement = disp_z
    velocity = vel_z
    beta = 0.25
    execute_on = timestep_end
  [../]
  [./vel_z]
    type = NewmarkVelAux
    variable = vel_z
    acceleration = accel_z
    gamma = 0.5
    execute_on = timestep_end
  [../]
[]

[Kernels]
  [./damper_disp_x]
    type = StressDivergenceDamper
    block = '3'
    displacements = 'disp_x disp_y disp_z'
    component = 0
    variable = disp_x
    save_in = reaction_x
  [../]
  [./damper_disp_y]
    type = StressDivergenceDamper
    block = '3'
    displacements = 'disp_x disp_y disp_z'
    component = 1
    variable = disp_y
    save_in = reaction_y
  [../]
  [./damper_disp_z]
    type = StressDivergenceDamper
    block = '3'
    displacements = 'disp_x disp_y disp_z'
    component = 2
    variable = disp_z
    save_in = reaction_z
  [../]
[]

[NodalKernels]
  [./x_inertial]
    type = NodalTranslationalInertia
    variable = disp_x
    velocity = vel_x
    acceleration = accel_x
    boundary = '2 3'
    beta = 0.25
    gamma = 0.5
    mass = 25000
    alpha = 0
    eta = 0
  [../]
  [./y_inertial]
    type = NodalTranslationalInertia
    variable = disp_y
    velocity = vel_y
    acceleration = accel_y
    boundary = '2 3'
    beta = 0.25
    gamma = 0.5
    mass = 25000
    alpha = 0
    eta = 0
  [../]
  [./z_inertial]
    type = NodalTranslationalInertia
    variable = disp_y
    velocity = vel_y
    acceleration = accel_y
    boundary = '2 3'
    beta = 0.25
    gamma = 0.5
    mass = 25000
    alpha = 0
    eta = 0
  [../]
[]

[Physics/SolidMechanics/LineElement/QuasiStatic]
  displacements = 'disp_x disp_y disp_z'
  rotations = 'rot_x rot_y rot_z'
  velocities = 'vel_x vel_y vel_z'
  accelerations = 'accel_x accel_y accel_z'
  rotational_velocities = 'rot_vel_x rot_vel_y rot_vel_z'
  rotational_accelerations = 'rot_accel_x rot_accel_y rot_accel_z'
  save_in = 'reaction_x reaction_y reaction_z'

  beta = 0.25 # Newmark time integration parameter
  gamma = 0.5 # Newmark time integration parameter

  # parameters for 5% Rayleigh damping
  zeta = 0  # stiffness proportional damping
  eta = 0   # Mass proportional Rayleigh damping

  [./block_1]
    block = 1
    area = 9e-2
    Iy = 6.75e-4
    Iz = 6.75e-4
    y_orientation = '0.0 0.0 1.0'
  [../]
  [./block_2]
    block = 2
    area = 9e-2
    Iy = 6.75e-1
    Iz = 6.75e-1
    y_orientation = '0.0 1.0 0.0'
  [../]
[]

[Materials]
  [./elasticity_beamcolumn]
    type = ComputeElasticityBeam
    youngs_modulus = 20e9
    poissons_ratio = -0.9
    shear_coefficient = 1
    block = '1 2'
  [../]
  [./stress_beam]
    type = ComputeBeamResultants
    block = '1 2'
  [../]
  [./elasticity_damper]
    type = ComputeFVDamperElasticity
    block = '3'
    y_orientation = '-0.5144 0.85749 0'
    displacements = 'disp_x disp_y disp_z'
    cd = 232764
    alpha = 0.35
    k = 25000000
    gamma = 0.5
    beta = 0.25
  [../]
[]

[BCs]
  [./base_accel_x]
    type = PresetAcceleration
    boundary = 1
    function = acceleration_x
    variable = disp_x
    beta = 0.25
    acceleration = 'accel_x'
    velocity = 'vel_x'
  [../]
  [./fix_y]
    type = DirichletBC
    variable = disp_y
    boundary = 1
    value = 0.0
  [../]
  [./fix_z]
    type = DirichletBC
    variable = disp_z
    boundary = 1
    value = 0.0
  [../]
  [./fix_rot_x]
    type = DirichletBC
    variable = rot_x
    boundary = 1
    value = 0.0
  [../]
  [./fix_rot_y]
    type = DirichletBC
    variable = rot_y
    boundary = 1
    value = 0.0
  [../]
  [./fix_rot_z]
    type = DirichletBC
    variable = rot_z
    boundary = 1
    value = 0.0
  [../]
[]

[Functions]
  [./acceleration_x]
    type = PiecewiseLinear
    data_file = acceleration.csv
    format = columns
  [../]
[]

[Preconditioning]
  [./smp]
    type = SMP
    full = true
  [../]
[]

[Executioner]
  type = Transient
  solve_type = NEWTON
  line_search = none
  nl_rel_tol = 1e-8
  nl_abs_tol = 1e-8
  start_time = 0
  end_time = 0.2
  dt = 0.0005
  dtmin = 0.000001
  timestep_tolerance = 1e-6
[]

[Postprocessors]
  # For basicdisp and basicforce postprocessors, divide the output values
  # by a scale factor equal to 'length of the damper element; (5.83 for this
  # problem). Becasue the object 'ElementIntegralMaterialProperty', calculates
  # the material property value integrating over the element length.
  [./basicdisp]
    type = ElementIntegralMaterialProperty
    mat_prop = basic_deformation
    block = '3'
  [../]
  [./basicforce]
    type = ElementIntegralMaterialProperty
    mat_prop = basic_force
    block = '3'
  [../]
  [./disp_bottom_right]
    type = PointValue
    variable = disp_x
    point = '5 0 0'
  [../]
  [./disp_top_right]
    type = PointValue
    variable = disp_x
    point = '5 3 0'
  [../]
[]

[Outputs]
  csv = true
  exodus = true
  perf_graph = true
  interval = 10
[]

(test/tests/materials/fv_damper/fv_damper_axial.i)

# Test for Damper in Axial (Sinewave loading)

[Mesh]
  type = GeneratedMesh
  xmin = 0
  xmax = 1
  nx = 1
  dim = 1
  displacements = 'disp_x disp_y disp_z'
[]

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

[AuxVariables]
  [./vel_x]
    order = FIRST
    family = LAGRANGE
  [../]
  [./vel_y]
    order = FIRST
    family = LAGRANGE
  [../]
  [./vel_z]
    order = FIRST
    family = LAGRANGE
  [../]
  [./accel_x]
    order = FIRST
    family = LAGRANGE
  [../]
  [./accel_y]
    order = FIRST
    family = LAGRANGE
  [../]
  [./accel_z]
    order = FIRST
    family = LAGRANGE
  [../]
  [./reaction_x]
  [../]
  [./reaction_y]
  [../]
  [./reaction_z]
  [../]
[]

[Kernels]
  [./lr_disp_x]
    type = StressDivergenceDamper
    block = '0'
    displacements = 'disp_x disp_y disp_z'
    component = 0
    variable = disp_x
    save_in = reaction_x
  [../]
  [./lr_disp_y]
    type = StressDivergenceDamper
    block = '0'
    displacements = 'disp_x disp_y disp_z'
    component = 1
    variable = disp_y
    save_in = reaction_y
  [../]
  [./lr_disp_z]
    type = StressDivergenceDamper
    block = '0'
    displacements = 'disp_x disp_y disp_z'
    component = 2
    variable = disp_z
    save_in = reaction_z
  [../]
[]

[AuxKernels]
  [./accel_x]
    type = NewmarkAccelAux
    variable = accel_x
    displacement = disp_x
    velocity = vel_x
    beta = 0.25
    execute_on = timestep_end
  [../]
  [./vel_x]
    type = NewmarkVelAux
    variable = vel_x
    acceleration = accel_x
    gamma = 0.5
    execute_on = timestep_end
  [../]
  [./accel_y]
    type = NewmarkAccelAux
    variable = accel_y
    displacement = disp_y
    velocity = vel_y
    beta = 0.25
    execute_on = timestep_end
  [../]
  [./vel_y]
    type = NewmarkVelAux
    variable = vel_y
    acceleration = accel_y
    gamma = 0.5
    execute_on = timestep_end
  [../]
  [./accel_z]
    type = NewmarkAccelAux
    variable = accel_z
    displacement = disp_z
    velocity = vel_z
    beta = 0.25
    execute_on = timestep_end
  [../]
  [./vel_z]
    type = NewmarkVelAux
    variable = vel_z
    acceleration = accel_z
    gamma = 0.5
    execute_on = timestep_end
  [../]
[]

[Materials]
  [./elasticity]
    type = ComputeFVDamperElasticity
    block = 0
    y_orientation = '0.0 1.0 0.0'
    displacements = 'disp_x disp_y disp_z'
    cd = 2226502
    alpha = 0.3
    k =300000000
    gamma = 0.5
    beta = 0.25
  [../]
[]

[BCs]
  [./fixx0]
    type = DirichletBC
    variable = disp_x
    boundary = left
    value = 0.0
  [../]
  [./fixy0]
    type = DirichletBC
    variable = disp_y
    boundary = left
    value = 0.0
  [../]
  [./fixz0]
    type = DirichletBC
    variable = disp_z
    boundary = left
    value = 0.0
  [../]
  [./disp_x_1]
    type = PresetDisplacement
    boundary = right
    function = displacement
    variable = disp_x
    beta = 0.25
    acceleration = accel_x
    velocity = vel_x
  [../]
[]

[Functions]
  [./displacement]
    type = PiecewiseLinear
    data_file = disp_axial.csv
    format = columns
  [../]
[]

[Preconditioning]
  [./smp]
    type = SMP
    full = true
  [../]
[]

[Executioner]
  type = Transient
  solve_type = NEWTON
  line_search = none
  nl_rel_tol = 1e-8
  nl_abs_tol = 1e-8
  start_time = -0.02
  end_time = 2
  dt = 0.0025
  dtmin = 0.00001
  timestep_tolerance = 1e-6
[]

[Postprocessors]
  [./disp_x]
    type = NodalVariableValue
    nodeid = 1
    variable = disp_x
  [../]
  [./vel_x]
    type = NodalVariableValue
    nodeid = 1
    variable = vel_x
  [../]
  [./accel_x]
    type = NodalVariableValue
    nodeid = 1
    variable = accel_x
  [../]
  [./reaction_x]
    type = NodalSum
    variable = reaction_x
    boundary = left
  [../]
[]

[Outputs]
  csv = true
  exodus = true
  interval = 2
  perf_graph = true
[]

Description
Input Parameters
Input Files
References

Usage

Development

Demonstration