diff --git a/test/tests/rebar_bondslip/RCBeam_constraint.i b/test/tests/rebar_bondslip/RCBeam_constraint.i new file mode 100644 index 000000000..c968f9b22 --- /dev/null +++ b/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 +[] diff --git a/test/tests/rebar_bondslip/RCBeam_constraint_equal.i b/test/tests/rebar_bondslip/RCBeam_constraint_equal.i new file mode 100644 index 000000000..ef557f4ae --- /dev/null +++ b/test/tests/rebar_bondslip/RCBeam_constraint_equal.i @@ -0,0 +1,267 @@ +[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 = small + strain = finite + incremental = true + # add_variables = 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_x2] + type = EqualValueEmbeddedConstraint + slave = 2 + master = 1 + penalty = 1e15 + variable = 'disp_x' + master_variable = 'disp_x' + formulation = penalty + [../] + [./rebar_y2] + type = EqualValueEmbeddedConstraint + slave = 2 + master = 1 + penalty = 1e15 + variable = 'disp_y' + master_variable = 'disp_y' + formulation = penalty + [../] +[] + +[Functions] + [./loading] + type = PiecewiseLinear + x = '0 10 20 30' + y = '0 0.0001 -0.0001 0.0' + [../] +[] + +[BCs] + [./loading] + type = FunctionDirichletBC + variable = disp_x + boundary = '102' + function = loading + 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 + # [] + [./truss] + type = PlasticTruss + youngs_modulus = 2.0e11 + yield_stress = 500e5 + hardening_constant = 100 + block = '2' + outputs = exodus + [../] +[] + +[Preconditioning] + [./SMP] + type = SMP + full = true + [../] +[] + +[Executioner] + type = Transient + solve_type = 'PJFNK' + nl_max_its = 20 + nl_abs_tol = 1e-6 + nl_rel_tol = 1e-05 + l_tol = 1e-03 + + 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] + exodus = true + csv = true +[] diff --git a/test/tests/rebar_bondslip/RCBeam_constraint_test.i b/test/tests/rebar_bondslip/RCBeam_constraint_test.i new file mode 100644 index 000000000..33c06a800 --- /dev/null +++ b/test/tests/rebar_bondslip/RCBeam_constraint_test.i @@ -0,0 +1,279 @@ +[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 = small + strain = finite + incremental = true + # add_variables = 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 + 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 = 1e5 + transitional_slip_values = 0.00005 + ultimate_slip = 0.5 + rebar_radius = 7.98e-3 + frictional_bondstress = 100 + [] + [rebar_y] + type = RebarBondSlipConstraint + slave = 2 + master = 1 + penalty = 1e6 + variable = 'disp_y' + master_variable = 'disp_y' + component = 1 + max_bondstress = 1e5 + transitional_slip_values = 0.00005 + ultimate_slip = 0.5 + rebar_radius = 7.98e-3 + frictional_bondstress = 100 + [] +[] + +[Functions] + [./loading] + type = PiecewiseLinear + x = '0 10 20 30 ' + y = '0 0.0001 -0.0001 0.0' + [../] +[] + +[BCs] + [./loading] + type = FunctionDirichletBC + variable = disp_x + boundary = '102' + function = loading + 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 + # boundary = '1001' + 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 + # [] + [./truss] + type = PlasticTruss + youngs_modulus = 2.0e11 + yield_stress = 500e5 + hardening_constant = 0 + block = '2' + outputs = exodus + [../] +[] + +[Preconditioning] + [./SMP] + type = SMP + full = true + [../] +[] + +[Executioner] + type = Transient + solve_type = 'PJFNK' + nl_max_its = 20 + nl_abs_tol = 1e-6 + nl_rel_tol = 1e-05 + l_tol = 1e-03 + + 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 +[] diff --git a/test/tests/rebar_bondslip/concrete_rebar_bond_test.i b/test/tests/rebar_bondslip/concrete_rebar_bond_test.i new file mode 100644 index 000000000..d9c97e5a0 --- /dev/null +++ b/test/tests/rebar_bondslip/concrete_rebar_bond_test.i @@ -0,0 +1,318 @@ +[GlobalParams] + displacements = 'disp_x disp_y' +[] + +[Mesh] + file = gold/geo.e +[] + +[Variables] + [disp_x] + [] + [disp_y] + [] +[] + +# [Functions] +# [E] +# type = PiecewiseMultilinear +# data_file = 'gold/E.txt' +# [] +# [] + +[AuxVariables] + [stress_xx] + order = CONSTANT + family = MONOMIAL + [] + [stress_yy] + order = CONSTANT + family = MONOMIAL + [] + [strain_xx] + order = CONSTANT + family = MONOMIAL + [] + [strain_yy] + order = CONSTANT + family = MONOMIAL + [] + [area] + order = CONSTANT + family = MONOMIAL + [] + [damage_index] + order = CONSTANT + family = MONOMIAL + [] + [./axial_stress] + order = CONSTANT + family = MONOMIAL + [../] +[] + +[AuxKernels] + [stress_xx] + type = RankTwoAux + index_i = 0 + index_j = 0 + variable = stress_xx + execute_on = 'TIMESTEP_END' + rank_two_tensor = stress + block = 1 + [] + [stress_yy] + type = RankTwoAux + index_i = 1 + index_j = 1 + variable = stress_yy + execute_on = 'TIMESTEP_END' + rank_two_tensor = stress + block = 1 + [] + [strain_xx] + type = RankTwoAux + index_i = 0 + index_j = 0 + variable = strain_xx + execute_on = 'TIMESTEP_END' + rank_two_tensor = total_strain + block = 1 + [] + [strain_yy] + type = RankTwoAux + index_i = 1 + index_j = 1 + variable = strain_yy + execute_on = 'TIMESTEP_END' + rank_two_tensor = total_strain + block = 1 + [] + [area] + type = ConstantAux + block = 2 + variable = area + value = 7.07e-4 + execute_on = 'initial' + [] + [damage_index] + type = MaterialRealAux + variable = damage_index + property = damage_index + block = 1 + execute_on = timestep_end + [] + [./axial_stress] + type = MaterialRealAux + block = '2' + variable = axial_stress + property = axial_stress + [../] +[] + +[Materials] + [concrete_strain] + type = ComputeIncrementalSmallStrain + block = 1 + [] + [Cijkl_concrete] + type = ComputeIsotropicElasticityTensor + youngs_modulus = 500e6 + poissons_ratio = 0.2 + block = 1 + [] + [damage] + type = MazarsDamage + use_old_damage = true + tensile_strength = 1e6 + a_t = 0.87 + a_c = 0.65 + b_t = 20000 + b_c = 2150 + block = 1 + [] + [concrete_stress] + type = ComputeDamageStress + damage_model = damage + block = 1 + [] + # [concrete_stress] + # type = ComputeFiniteStrainElasticStress + # block = 1 + # [] + # [truss] + # type = LinearElasticTruss + # block = 2 + # youngs_modulus = 2e11 + # [] + [./truss] + type = PlasticTruss + youngs_modulus = 2.0e11 + yield_stress = 500e5 + hardening_constant = 100 + block = '2' + outputs = exodus + [../] +[] + +[Postprocessors] + [./stress_xx] + type = ElementAverageValue + variable = stress_xx + block = '1' + [../] + [./strain_xx] + type = ElementAverageValue + variable = strain_xx + block = '1' + [../] + [./stress_yy] + type = ElementAverageValue + variable = stress_yy + block = '1' + [../] + [./strain_yy] + type = ElementAverageValue + variable = strain_yy + block = '1' + [../] + [./axial_stress] + type = ElementAverageValue + variable = axial_stress + block = '2' + [../] +[] + +[Preconditioning] + [./SMP] + type = SMP + full = true + [../] +[] + +[Executioner] + type = Transient + solve_type = PJFNK + line_search = basic + # petsc_options_iname = '-pc_type -pc_hypre_type -ksp_gmres_restart -mat_mffd_type' + # petsc_options_value = 'hypre boomeramg 101 ds' + + # petsc_options_iname = '-pc_type -ksp_gmres_restart' + # petsc_options_value = 'jacobi 101' + petsc_options_iname = '-pc_type' + petsc_options_value = 'lu' + nl_rel_tol = 1e-5 + nl_abs_tol = 1e-6 + l_tol = 1e-03 + l_max_its = 100 + nl_max_its = 30 + dt = 2e-5 + end_time = 1 +[] + +[Outputs] + exodus = true + print_linear_residuals = false + file_base = concrete_rebar_bond_test +[] + +[Kernels] + [solid_x_concrete] + type = StressDivergenceTensors + block = 1 + component = 0 + variable = disp_x + use_displaced_mesh = true + [] + [solid_y_concrete] + type = StressDivergenceTensors + block = 1 + component = 1 + variable = disp_y + use_displaced_mesh = true + [] + [solid_x_truss] + type = StressDivergenceTensorsTruss + block = 2 + component = 0 + variable = disp_x + area = area + [] + [solid_y_truss] + type = StressDivergenceTensorsTruss + block = 2 + component = 1 + variable = disp_y + area = area + [] +[] + +[BCs] + [concrete_top_ydisp] + type = FunctionDirichletBC + variable = disp_y + boundary = '1' + function = '-t' + [] + [concrete_bottom_yfix] + type = DirichletBC + variable = disp_y + boundary = '2 3' + value = 0 + [] + [concrete_bottom_xfix] + type = DirichletBC + variable = disp_x + boundary = '2' + value = 0 + [] +[] + +[Constraints] + [rebar_x] + type = RebarBondSlipConstraint + slave = 2 + master = 1 + penalty = 1e6 + variable = 'disp_x' + master_variable = 'disp_x' + component = 0 + max_bondstress = 3.5e-10 + transitional_slip_values = 0.05e-3 + ultimate_slip = 1.0 + rebar_radius = 0.015 + frictional_bondstress = 1e-15 + [] + [rebar_y] + type = RebarBondSlipConstraint + slave = 2 + master = 1 + penalty = 1e6 + variable = 'disp_y' + master_variable = 'disp_y' + component = 1 + max_bondstress = 3.5e-10 + transitional_slip_values = 0.05e-3 + ultimate_slip = 1.0 + rebar_radius = 0.015 + frictional_bondstress = 1e-15 + [] + # [rebar_x] + # type = EqualValueEmbeddedConstraint + # slave = 2 + # master = 1 + # penalty = 1e6 + # variable = 'disp_x' + # master_variable = 'disp_x' + # formulation = penalty + # [] + # [rebar_y] + # type = EqualValueEmbeddedConstraint + # slave = 2 + # master = 1 + # penalty = 1e6 + # variable = 'disp_y' + # master_variable = 'disp_y' + # formulation = penalty + # [] +[] diff --git a/test/tests/rebar_bondslip/gold/RCBeam_test_mesh.e b/test/tests/rebar_bondslip/gold/RCBeam_test_mesh.e new file mode 100644 index 000000000..dcfb3119f Binary files /dev/null and b/test/tests/rebar_bondslip/gold/RCBeam_test_mesh.e differ diff --git a/test/tests/rebar_bondslip/gold/geo.e b/test/tests/rebar_bondslip/gold/geo.e new file mode 100644 index 000000000..876b78ff2 Binary files /dev/null and b/test/tests/rebar_bondslip/gold/geo.e differ