BondSlip Initial Tests #99

test/tests/rebar_bondslip/RCBeam_constraint.i

@@ -0,0 +1,268 @@
+[Mesh]
+  file = gold/RCBeam_test_mesh.e
+[]
+
+[GlobalParams]
+  displacements = 'disp_x disp_y'
+  # volumetric_locking_correction = true
+[]
+
+[Modules/TensorMechanics/Master]
+  [./Concrete_block]
+    block = 1
+    strain = finite
+    incremental = true
+    generate_output = 'stress_xx stress_xy stress_yy strain_xx strain_xy strain_yy
+    		       max_principal_stress mid_principal_stress min_principal_stress
+    		       secondinv_stress thirdinv_stress vonmises_stress
+    		       secondinv_strain thirdinv_strain
+    		       elastic_strain_xx elastic_strain_xy elastic_strain_yy'
+#    		       plastic_strain_xx plastic_strain_xy plas tic_strain_xz plastic_strain_yy plastic_strain_yz plastic_strain_zz'
+    save_in = 'resid_x resid_y'
+  [../]
+[]
+
+[Modules/TensorMechanics/LineElementMaster]
+  [./Reinforcement_block]
+    block = '2'
+    truss = true
+    area = area
+    displacements = 'disp_x disp_y'
+    save_in = 'resid_x resid_y'
+   # add_variables = true
+  [../]
+[]
+
+[Variables]
+  [./disp_x]
+  [../]
+  [./disp_y]
+  [../]
+[]
+
+[AuxVariables]
+  [./resid_x]
+  [../]
+  [./resid_y]
+  [../]
+  [./area]
+    order = CONSTANT
+    family = MONOMIAL
+  [../]
+  [./axial_stress]
+    order = CONSTANT
+    family = MONOMIAL
+  [../]
+[]
+
+[AuxKernels]
+  [./area]
+    type = ConstantAux
+    block = '2'
+    variable = area
+    value = 2.00e-4 # 509 mm2
+    execute_on = 'initial timestep_begin'
+  [../]
+  [./axial_stress]
+    type = MaterialRealAux
+    block = '2'
+    variable = axial_stress
+    property = axial_stress
+  [../]
+[]
+
+[Constraints]
+  [rebar_x]
+    type = RebarBondSlipConstraint
+    slave = 2
+    master = 1
+    penalty = 1e6
+    variable = 'disp_x'
+    master_variable = 'disp_x'
+    component = 0
+    max_bondstress = 100
+    transitional_slip_values = 0.0005
+    ultimate_slip = 0.1
+    rebar_radius = 7.98e-3
+  []
+  [rebar_y]
+    type = RebarBondSlipConstraint
+    slave = 2
+    master = 1
+    penalty = 1e6
+    variable = 'disp_y'
+    master_variable = 'disp_y'
+    component = 1
+    max_bondstress = 100
+    transitional_slip_values = 0.0005
+    ultimate_slip = 0.1
+    rebar_radius = 7.98e-3
+  []
+[]
+
+[Functions]
+  [./loading]
+    type = PiecewiseLinear
+    x = '0 10       20     30 '
+    y = '0 -0.008 0.003    0.0'
+  [../]
+[]
+
+[BCs]
+  [./loading]
+    type = FunctionDirichletBC
+    # type = DirichletBC
+    variable = disp_x
+    boundary = '102'
+    function = loading
+    # value = 0.00004
+    preset = true
+  [../]
+  [./left_support_x]
+    type = DirichletBC
+    variable = disp_x
+    boundary = '100'
+    value = 0
+  [../]
+  [./left_support_y]
+    type = DirichletBC
+    variable = disp_y
+    boundary = '100'
+    value = 0
+  [../]
+[]
+
+[Postprocessors]
+  [./deformation_x]
+    type = AverageNodalVariableValue
+    variable = disp_x
+    boundary = '101'
+  [../]
+  [./deformation_y]
+    type = AverageNodalVariableValue
+    variable = disp_y
+    boundary = '101'
+  [../]
+  [./react_x]
+    type = AverageNodalVariableValue
+    variable = resid_x
+    boundary = '100'
+  [../]
+  [./react_y]
+    type = AverageNodalVariableValue
+    variable = resid_y
+    boundary = '100'
+  [../]
+  [./react_x2]
+    type = AverageNodalVariableValue
+    variable = resid_x
+    boundary = '101'
+  [../]
+  [./react_y2]
+    type = AverageNodalVariableValue
+    variable = resid_y
+    boundary = '100'
+  [../]
+  [./node1_fx]
+    type = NodalVariableValue
+    variable = resid_x
+    nodeid = 138
+  [../]
+  [./node1_fy]
+    type = NodalVariableValue
+    variable = resid_y
+    nodeid = 138
+  [../]
+  [./node1_dx]
+    type = NodalVariableValue
+    variable = disp_x
+    nodeid = 138
+  [../]
+  [./node1_dy]
+    type = NodalVariableValue
+    variable = disp_y
+    nodeid = 138
+  [../]
+  [./node1_fx2]
+    type = AverageNodalVariableValue
+    variable = resid_x
+    boundary = '102'
+  [../]
+
+  [./stress_xx]
+    type = ElementAverageValue
+    variable = stress_xx
+    block = '1'
+  [../]
+  [./strain_xx]
+    type = ElementAverageValue
+    variable = strain_xx
+    block = '1'
+  [../]
+  [./axial_stress]
+    type = ElementAverageValue
+    variable = axial_stress
+    block = '2'
+  [../]
+[]
+
+[Materials]
+  [Cijkl_concrete]
+    type = ComputeIsotropicElasticityTensor
+    youngs_modulus = 500e6
+    poissons_ratio = 0.2
+    block = 1
+  []
+  [./isotropic_plasticity]
+    type = IsotropicPlasticityStressUpdate
+    yield_stress = 285788383.2488647 # = sqrt(3)*165e6 = sqrt(3) * yield in shear
+    hardening_constant = 0.0
+    block = '1'
+  [../]
+  [./radial_return_stress]
+    type = ComputeMultipleInelasticStress
+    tangent_operator = elastic
+    inelastic_models = 'isotropic_plasticity'
+    block = '1'
+  [../]
+  [truss]
+    type = LinearElasticTruss
+    block = '2'
+    youngs_modulus = 2e11
+  []
+[]
+
+[Preconditioning]
+  [./SMP]
+    type = SMP
+    full = true
+  [../]
+[]
+
+[Executioner]
+  type = Transient
+  solve_type = 'PJFNK'
+  nl_max_its = 100
+  nl_abs_tol = 1.E-5
+  nl_rel_tol = 1E-3
+
+  line_search = none
+
+  petsc_options_iname = '-pc_type'
+  petsc_options_value = 'lu'
+
+  petsc_options = '-snes_converged_reason'
+
+  end_time = 30
+  dtmin = 0.00001
+
+  dt = 0.1
+[]
+
+
+[Outputs]
+  # print_linear_residuals = false
+  exodus = true
+  csv = true
+  # file_base = RCBeam_bondslip_test
+[]