diff --git a/test/test_model/test_solid_mechanics_model/test_contact/test_contact_rigid/CMakeLists.txt b/test/test_model/test_solid_mechanics_model/test_contact/test_contact_rigid/CMakeLists.txt index 15c1060b9..71656fecd 100644 --- a/test/test_model/test_solid_mechanics_model/test_contact/test_contact_rigid/CMakeLists.txt +++ b/test/test_model/test_solid_mechanics_model/test_contact/test_contact_rigid/CMakeLists.txt @@ -1,52 +1,54 @@ #=============================================================================== # @file CMakeLists.txt # @author David Kammer # @date Wed Jan 19 12:37:24 2011 # # @section LICENSE # # Copyright (©) 2010-2011 EPFL (Ecole Polytechnique Fédérale de Lausanne) # Laboratory (LSMS - Laboratoire de Simulation en Mécanique des Solides) # # Akantu is free software: you can redistribute it and/or modify it under the # terms of the GNU Lesser General Public License as published by the Free # Software Foundation, either version 3 of the License, or (at your option) any # later version. # # Akantu is distributed in the hope that it will be useful, but WITHOUT ANY # WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR # A PARTICULAR PURPOSE. See the GNU Lesser General Public License for more # details. # # You should have received a copy of the GNU Lesser General Public License # along with Akantu. If not, see . # # @section DESCRIPTION # #=============================================================================== add_subdirectory("patch_tests") #=============================================================================== set(LIST_TYPES triangle_3 quadrangle_4 tetrahedron_4 #hexahedron_8 ) #=============================================================================== foreach(_type ${LIST_TYPES}) register_test(test_contact_rigid_no_friction_${_type} test_contact_rigid_no_friction.cc io_helper_tools.cc) set_target_properties(test_contact_rigid_no_friction_${_type} PROPERTIES COMPILE_DEFINITIONS TYPE=_${_type}) file(COPY ./data/${_type}.msh DESTINATION .) + file(COPY ./ref_files/test_contact_rigid_no_friction_${_type}.ref DESTINATION .) endforeach() #=============================================================================== -register_test(test_contact_rigid_restart_triangle_3 test_contact_rigid_restart_triangle_3.cc) +register_test(test_contact_rigid_restart test_contact_rigid_restart.cc) +file(COPY ./ref_files/test_contact_rigid_restart.ref DESTINATION .) #=============================================================================== file(COPY material.dat DESTINATION .) file(MAKE_DIRECTORY ${CMAKE_CURRENT_BINARY_DIR}/paraview) file(MAKE_DIRECTORY ${CMAKE_CURRENT_BINARY_DIR}/restart) diff --git a/test/test_model/test_solid_mechanics_model/test_contact/test_contact_rigid/ref_files/test_contact_rigid_no_friction_quadrangle_4.ref b/test/test_model/test_solid_mechanics_model/test_contact/test_contact_rigid/ref_files/test_contact_rigid_no_friction_quadrangle_4.ref new file mode 100644 index 000000000..cb38f8cc6 --- /dev/null +++ b/test/test_model/test_solid_mechanics_model/test_contact/test_contact_rigid/ref_files/test_contact_rigid_no_friction_quadrangle_4.ref @@ -0,0 +1,33 @@ + +passing step 1/3 +Before solveContact +we have 0 penetrating nodes: + +After solveContact +we have 0 penetrating nodes: +we have 0 active impactor nodes: + +passing step 2/3 +Before solveContact +we have 3 penetrating nodes: +node 4 with disp: -0.02 0 +node 7 with disp: -0.02 0 +node 16 with disp: -0.02 0 + +After solveContact +we have 0 penetrating nodes: +we have 3 active impactor nodes: +node 4 with disp: -0.005 0 and force: -2120192308 -908653846.2 +node 7 with disp: -0.005 0 and force: -2120192308 908653846.2 +node 16 with disp: -0.005 0 and force: -4240384615 1.192092896e-07 + +passing step 3/3 +Before solveContact +we have 0 penetrating nodes: + +After solveContact +we have 0 penetrating nodes: +we have 3 active impactor nodes: +node 4 with disp: -0.005 -0.0001730769231 and force: -2082040496 34948224.85 +node 7 with disp: -0.005 0.0001730769231 and force: -2082040496 -34948224.85 +node 16 with disp: -0.005 1.135326567e-20 and force: -4171070636 -1.192092896e-07 diff --git a/test/test_model/test_solid_mechanics_model/test_contact/test_contact_rigid/ref_files/test_contact_rigid_no_friction_tetrahedron_4.ref b/test/test_model/test_solid_mechanics_model/test_contact/test_contact_rigid/ref_files/test_contact_rigid_no_friction_tetrahedron_4.ref new file mode 100644 index 000000000..5daa2f8e8 --- /dev/null +++ b/test/test_model/test_solid_mechanics_model/test_contact/test_contact_rigid/ref_files/test_contact_rigid_no_friction_tetrahedron_4.ref @@ -0,0 +1,60 @@ + +passing step 1/3 +Before solveContact +we have 0 penetrating nodes: + +After solveContact +we have 0 penetrating nodes: +we have 0 active impactor nodes: + +passing step 2/3 +Before solveContact +we have 12 penetrating nodes: +node 8 with disp: -0.02 0 0 +node 11 with disp: -0.02 0 0 +node 12 with disp: -0.02 0 0 +node 15 with disp: -0.02 0 0 +node 44 with disp: -0.02 0 0 +node 51 with disp: -0.02 0 0 +node 52 with disp: -0.02 0 0 +node 55 with disp: -0.02 0 0 +node 126 with disp: -0.02 0 0 +node 127 with disp: -0.02 0 0 +node 128 with disp: -0.02 0 0 +node 129 with disp: -0.02 0 0 + +After solveContact +we have 0 penetrating nodes: +we have 12 active impactor nodes: +node 8 with disp: -0.005 0 0 and force: -237061729.3 -96454326.92 56790865.38 +node 11 with disp: -0.005 0 0 and force: -198497596.2 96454326.92 56790865.38 +node 12 with disp: -0.005 0 0 and force: -178410836.5 -83247862.81 -61736575.66 +node 15 with disp: -0.005 0 0 and force: -223669327.9 85740916.47 -59102833.53 +node 44 with disp: -0.005 0 0 and force: -301826674.9 694629.7749 122952820.2 +node 51 with disp: -0.005 0 0 and force: -259624073.6 3773540.109 -106576393.5 +node 52 with disp: -0.005 0 0 and force: -150361356 -135784899.9 5854896.109 +node 55 with disp: -0.005 0 0 and force: -161149085.5 137604093.9 8680563.982 +node 126 with disp: -0.005 0 0 and force: -345952305.3 -29886873.34 27566584.21 +node 127 with disp: -0.005 0 0 and force: -339103405.6 6801584.431 -36123541.65 +node 128 with disp: -0.005 0 0 and force: -287336292.2 29109532.75 31886839.75 +node 129 with disp: -0.005 0 0 and force: -356695350.7 -8455718.438 -37684537.28 + +passing step 3/3 +Before solveContact +we have 0 penetrating nodes: + +After solveContact +we have 0 penetrating nodes: +we have 12 active impactor nodes: +node 8 with disp: -0.005 -5.252148627e-06 3.092386575e-06 and force: -236906752.3 129012.8535 -69846.0467 +node 11 with disp: -0.005 5.645512477e-06 3.323993328e-06 and force: -198361709.3 -137713.7039 -78940.8853 +node 12 with disp: -0.005 -5.266916173e-06 -3.905942542e-06 and force: -178287070.3 118602.6904 79920.40268 +node 15 with disp: -0.005 4.551947184e-06 -3.137743188e-06 and force: -223524083.3 -108325.5697 60182.96939 +node 44 with disp: -0.005 2.57596845e-08 4.559588391e-06 and force: -301656046.8 24098.50197 -187317.2234 +node 51 with disp: -0.005 1.135979314e-07 -3.208355412e-06 and force: -259487112.5 -15987.83741 100522.7372 +node 52 with disp: -0.005 -3.195806575e-06 1.37799678e-07 and force: -150326656.1 112318.9415 4945.56448 +node 55 with disp: -0.005 3.296970855e-06 2.079848472e-07 and force: -161114758 -107967.1748 499.8248915 +node 126 with disp: -0.005 -4.113827814e-07 3.794447802e-07 and force: -345878463.4 114639.8229 -116519.9735 +node 127 with disp: -0.005 1.336221257e-07 -7.096735285e-07 and force: -339013371.3 -73880.43674 159358.0287 +node 128 with disp: -0.005 6.226939363e-07 6.821044476e-07 and force: -287251508.1 -130220.7594 -110290.0135 +node 129 with disp: -0.005 -2.989869283e-07 -1.332492813e-06 and force: -356549334.8 65161.27992 110922.0719 diff --git a/test/test_model/test_solid_mechanics_model/test_contact/test_contact_rigid/ref_files/test_contact_rigid_no_friction_triangle_3.ref b/test/test_model/test_solid_mechanics_model/test_contact/test_contact_rigid/ref_files/test_contact_rigid_no_friction_triangle_3.ref new file mode 100644 index 000000000..4223c6451 --- /dev/null +++ b/test/test_model/test_solid_mechanics_model/test_contact/test_contact_rigid/ref_files/test_contact_rigid_no_friction_triangle_3.ref @@ -0,0 +1,33 @@ + +passing step 1/3 +Before solveContact +we have 0 penetrating nodes: + +After solveContact +we have 0 penetrating nodes: +we have 0 active impactor nodes: + +passing step 2/3 +Before solveContact +we have 3 penetrating nodes: +node 4 with disp: -0.02 0 +node 7 with disp: -0.02 0 +node 16 with disp: -0.02 0 + +After solveContact +we have 0 penetrating nodes: +we have 3 active impactor nodes: +node 4 with disp: -0.005 0 and force: -3483173077 -1665865385 +node 7 with disp: -0.005 0 and force: -3483173077 1665865385 +node 16 with disp: -0.005 0 and force: -6966346154 0 + +passing step 3/3 +Before solveContact +we have 0 penetrating nodes: + +After solveContact +we have 0 penetrating nodes: +we have 3 active impactor nodes: +node 4 with disp: -0.005 -2.041552242e-05 and force: -3474797181 7176662.884 +node 7 with disp: -0.005 2.041552242e-05 and force: -3474797181 -7176662.884 +node 16 with disp: -0.005 0 and force: -6953323228 -2.384185791e-07 diff --git a/test/test_model/test_solid_mechanics_model/test_contact/test_contact_rigid/ref_files/test_contact_rigid_restart.ref b/test/test_model/test_solid_mechanics_model/test_contact/test_contact_rigid/ref_files/test_contact_rigid_restart.ref new file mode 100644 index 000000000..fac652022 --- /dev/null +++ b/test/test_model/test_solid_mechanics_model/test_contact/test_contact_rigid/ref_files/test_contact_rigid_restart.ref @@ -0,0 +1,204 @@ +CONTACT +Active impactor nodes (contact): + +node: 4, master element type: segment_2 +Direction 0 + master normal = 1 + friction force = -0 + stick position = 0.5 + residual force = -3.47523e+09 + previous velocity = 0 +Direction 1 + master normal = 0 + friction force = 3.47523e+08 + stick position = 0.375 + residual force = -3.41916e+08 + previous velocity = -48.4752 +stick 0: 0 +stick 1: 0 + +node: 7, master element type: segment_2 +Direction 0 + master normal = 1 + friction force = -0 + stick position = 0.5 + residual force = -3.47523e+09 + previous velocity = 0 +Direction 1 + master normal = 0 + friction force = -3.47523e+08 + stick position = 0.625 + residual force = 3.41916e+08 + previous velocity = 48.4752 +stick 0: 0 +stick 1: 0 + +node: 16, master element type: segment_2 +Direction 0 + master normal = 1 + friction force = -0 + stick position = 0.5 + residual force = -6.95341e+09 + previous velocity = 0 +Direction 1 + master normal = 0 + friction force = -0 + stick position = 0.5 + residual force = 0 + previous velocity = 0 +stick 0: 1 +stick 1: 1 + +MAP +Active impactor nodes (map): + +node: 4, master element type: segment_2 +Direction 0 + master normal = 1 + friction force = -0 + stick position = 0.5 + residual force = -3.47523e+09 + previous velocity = 0 +Direction 1 + master normal = 0 + friction force = 3.47523e+08 + stick position = 0.375 + residual force = -3.41916e+08 + previous velocity = -48.4752 +stick 0: 0 +stick 1: 0 + +node: 7, master element type: segment_2 +Direction 0 + master normal = 1 + friction force = -0 + stick position = 0.5 + residual force = -3.47523e+09 + previous velocity = 0 +Direction 1 + master normal = 0 + friction force = -3.47523e+08 + stick position = 0.625 + residual force = 3.41916e+08 + previous velocity = 48.4752 +stick 0: 0 +stick 1: 0 + +node: 16, master element type: segment_2 +Direction 0 + master normal = 1 + friction force = -0 + stick position = 0.5 + residual force = -6.95341e+09 + previous velocity = 0 +Direction 1 + master normal = 0 + friction force = -0 + stick position = 0.5 + residual force = 0 + previous velocity = 0 +stick 0: 1 +stick 1: 1 + +RESTART MAP +Active impactor nodes (map): + +node: 4, master element type: segment_2 +Direction 0 + master normal = 1 + friction force = -0 + stick position = 0.5 + residual force = -3.47523e+09 + previous velocity = 0 +Direction 1 + master normal = 0 + friction force = 3.47523e+08 + stick position = 0.375 + residual force = -3.41916e+08 + previous velocity = -48.4752 +stick 0: 0 +stick 1: 0 + +node: 7, master element type: segment_2 +Direction 0 + master normal = 1 + friction force = -0 + stick position = 0.5 + residual force = -3.47523e+09 + previous velocity = 0 +Direction 1 + master normal = 0 + friction force = -3.47523e+08 + stick position = 0.625 + residual force = 3.41916e+08 + previous velocity = 48.4752 +stick 0: 0 +stick 1: 0 + +node: 16, master element type: segment_2 +Direction 0 + master normal = 1 + friction force = -0 + stick position = 0.5 + residual force = -6.95341e+09 + previous velocity = 0 +Direction 1 + master normal = 0 + friction force = -0 + stick position = 0.5 + residual force = 0 + previous velocity = 0 +stick 0: 1 +stick 1: 1 + +RESTART CONTACT +Active impactor nodes (contact): + +node: 4, master element type: segment_2 +Direction 0 + master normal = 1 + friction force = -0 + stick position = 0.5 + residual force = -3.47523e+09 + previous velocity = 0 +Direction 1 + master normal = 0 + friction force = 3.47523e+08 + stick position = 0.375 + residual force = -3.41916e+08 + previous velocity = -48.4752 +stick 0: 0 +stick 1: 0 + +node: 7, master element type: segment_2 +Direction 0 + master normal = 1 + friction force = -0 + stick position = 0.5 + residual force = -3.47523e+09 + previous velocity = 0 +Direction 1 + master normal = 0 + friction force = -3.47523e+08 + stick position = 0.625 + residual force = 3.41916e+08 + previous velocity = 48.4752 +stick 0: 0 +stick 1: 0 + +node: 16, master element type: segment_2 +Direction 0 + master normal = 1 + friction force = -0 + stick position = 0.5 + residual force = -6.95341e+09 + previous velocity = 0 +Direction 1 + master normal = 0 + friction force = -0 + stick position = 0.5 + residual force = 0 + previous velocity = 0 +stick 0: 1 +stick 1: 1 + diff --git a/test/test_model/test_solid_mechanics_model/test_contact/test_contact_rigid/test_contact_rigid_restart_triangle_3.cc b/test/test_model/test_solid_mechanics_model/test_contact/test_contact_rigid/test_contact_rigid_restart.cc similarity index 100% rename from test/test_model/test_solid_mechanics_model/test_contact/test_contact_rigid/test_contact_rigid_restart_triangle_3.cc rename to test/test_model/test_solid_mechanics_model/test_contact/test_contact_rigid/test_contact_rigid_restart.cc