diff --git a/include/constraints/RebarBondSlipConstraint.h b/include/constraints/RebarBondSlipConstraint.h index 2cf43ff8..06a2c0d9 100644 --- a/include/constraints/RebarBondSlipConstraint.h +++ b/include/constraints/RebarBondSlipConstraint.h @@ -99,9 +99,6 @@ class RebarBondSlipConstraint : public EqualValueEmbeddedConstraint std::vector _var_nums; std::vector _vars; - /// flag to turn on printing values for debugging - const bool _debug; - /// maximum bond stress const Real _max_bondstress; diff --git a/src/constraints/RebarBondSlipConstraint.C b/src/constraints/RebarBondSlipConstraint.C index d85416cd..42f7269a 100644 --- a/src/constraints/RebarBondSlipConstraint.C +++ b/src/constraints/RebarBondSlipConstraint.C @@ -33,7 +33,6 @@ RebarBondSlipConstraint::validParams() params.addCoupledVar( "displacements", "The displacements appropriate for the simulation geometry and coordinate system"); - params.addParam("debug", false, "whether to print out debug messages"); params.addParam("max_bondstress", 0.0, "Maximum bond stress"); params.addParam("frictional_bondstress", 0.0, "Bond stress due to friction"); @@ -53,7 +52,6 @@ RebarBondSlipConstraint::RebarBondSlipConstraint(const InputParameters & paramet _mesh_dimension(_mesh.dimension()), _var_nums(_mesh_dimension, libMesh::invalid_uint), _vars(_mesh_dimension, nullptr), - _debug(getParam("debug")), _max_bondstress(getParam("max_bondstress")), _frictional_bondstress(getParam("frictional_bondstress")), _ultimate_slip(getParam("ultimate_slip")), @@ -99,12 +97,6 @@ RebarBondSlipConstraint::timestepSetup() bool RebarBondSlipConstraint::shouldApply() { - if (_debug) - { - std::cout << "===========================================\n"; - std::cout << "node id: " << _current_node->id() << std::endl; - std::cout << "at coord: " << (Point)*_current_node << std::endl; - } auto it = _secondary_to_primary_map.find(_current_node->id()); if (it != _secondary_to_primary_map.end()) @@ -127,7 +119,6 @@ void RebarBondSlipConstraint::computeTangent() { _secondary_tangent = 0.0; - // get normals // get connected elements of the current node const std::map> & node_to_elem_map = _mesh.nodeToElemMap(); @@ -154,9 +145,6 @@ RebarBondSlipConstraint::computeTangent() _secondary_tangent /= _secondary_tangent.norm(); _current_elem_volume /= elems.size(); - - if (_debug) - std::cout << "tangent: " << _secondary_tangent << std::endl; } void @@ -174,9 +162,6 @@ RebarBondSlipConstraint::reinitConstraint() RealVectorValue slip_normal = relative_disp - slip_axial; Real slip_ratio = std::abs(slip) / _transitional_slip[0]; - if (_debug) - std::cout << "Slip = " << slip << ".\n"; - const Node * node = _current_node; auto it = _bondslip.find(node->id()); mooseAssert(it != _bondslip.end(), "Node not found in bond-slip map"); @@ -185,18 +170,6 @@ RebarBondSlipConstraint::reinitConstraint() bond_slip.slip_min = std::min(bond_slip.slip_min_old, slip); bond_slip.slip_max = std::max(bond_slip.slip_max_old, slip); - if (_debug) - { - std::cout << "Slip_min = " << bond_slip.slip_min << ".\n"; - std::cout << "Slip_min_old = " << bond_slip.slip_min_old << ".\n"; - std::cout << "Slip_max = " << bond_slip.slip_max << ".\n"; - std::cout << "Slip_max_old = " << bond_slip.slip_max_old << ".\n"; - std::cout << "Bondstress_min = " << bond_slip.bondstress_min << ".\n"; - std::cout << "Bondstress_min_old = " << bond_slip.bondstress_min_old << ".\n"; - std::cout << "Bondstress_max = " << bond_slip.bondstress_max << ".\n"; - std::cout << "Bondstress_max_old = " << bond_slip.bondstress_max_old << ".\n"; - } - Real slope = 5.0 * _max_bondstress / _transitional_slip[0]; Real plastic_slip_max = bond_slip.slip_max - bond_slip.bondstress_max / slope; Real plastic_slip_min = bond_slip.slip_min - bond_slip.bondstress_min / slope; @@ -238,12 +211,6 @@ RebarBondSlipConstraint::reinitConstraint() else _bond_stress = _frictional_bondstress; - if (_debug) - { - std::cout << "Bondstress = " << _bond_stress << "\n"; - std::cout << "Bondstress Derivative = " << _bond_stress_deriv << "\n"; - } - Real bond_force = 2.0 * libMesh::pi * _bar_radius * _current_elem_volume * _bond_stress; Real bond_force_deriv = 2.0 * libMesh::pi * _bar_radius * _current_elem_volume * _bond_stress_deriv; diff --git a/test/tests/rebar_bondslip/RCBeam_constraint.i b/test/tests/rebar_bondslip/RCBeam_constraint.i index efe874e3..d9a30650 100644 --- a/test/tests/rebar_bondslip/RCBeam_constraint.i +++ b/test/tests/rebar_bondslip/RCBeam_constraint.i @@ -4,7 +4,6 @@ [GlobalParams] displacements = 'disp_x disp_y' - # volumetric_locking_correction = true [] [Modules/TensorMechanics/Master] @@ -17,7 +16,6 @@ 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' [../] [] @@ -29,7 +27,6 @@ area = area displacements = 'disp_x disp_y' save_in = 'resid_x resid_y' - # add_variables = true [../] [] @@ -84,7 +81,6 @@ transitional_slip_values = 0.0005 ultimate_slip = 0.1 rebar_radius = 7.98e-3 - # debug = true [] [rebar_y] type = RebarBondSlipConstraint @@ -98,7 +94,6 @@ transitional_slip_values = 0.0005 ultimate_slip = 0.1 rebar_radius = 7.98e-3 - # debug = true [] [] @@ -113,11 +108,9 @@ [BCs] [./loading] type = FunctionDirichletBC - # type = DirichletBC variable = disp_x boundary = '102' function = loading - # value = 0.00004 preset = true [../] [./left_support_x] @@ -244,28 +237,21 @@ [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' - + nl_max_its = 100 + nl_abs_tol = 1.E-5 + nl_rel_tol = 1E-3 end_time = 30 dtmin = 0.00001 - # num_steps = 5 - + num_steps = 5 dt = 0.1 [] [Outputs] - # print_linear_residuals = false exodus = true csv = true - # file_base = RCBeam_bondslip_test []