hanging node tests

githubgrasp · githubgrasp · commit c6e5864b905b · 2025-09-08T14:13:06.000+01:00
diff --git a/src/oofemlib/hangingnode.C b/src/oofemlib/hangingnode.C
@@ -40,6 +40,7 @@
 #include "feinterpol.h"
 #include "spatiallocalizer.h"
 #include "classfactory.h"
+#include "elementgeometrytype.h"
 
 namespace oofem {
 REGISTER_DofManager( HangingNode );
@@ -169,11 +170,12 @@ void HangingNode ::postInitialize()
         }
     }
 
+    //Deal with rotations. Works only for linear tetrahedron so far.
+ if (e->giveGeometryType() == EGT_tetra_1) {
+   
     Dof *dof;
     const int nnodes = e->giveNumberOfNodes();
 
-
-    //Deal with rotations.
     // Compute shape function gradients once
     FloatMatrix dNdX;
     fei->evaldNdx(dNdX, lcoords, FEIElementGeometryWrapper(e));
@@ -296,6 +298,10 @@ void HangingNode ::postInitialize()
         }
     }
 
+} else {
+    OOFEM_WARNING("HangingNode: rotational DOFs only supported for EGT_tetra_1 (linear tetrahedron). Skipping rotational constraint.");
+}
+    
 
 }
 } // end namespace oofem
diff --git a/src/sm/Elements/LatticeElements/latticelink3d.C b/src/sm/Elements/LatticeElements/latticelink3d.C
@@ -318,10 +318,7 @@ namespace oofem {
         Node *nodeA, *nodeB;
         FloatArray coordsA(3), coordsB(3);
 
-        //Order of nodes. However, might not matter.
-        //Reinforcement node
         nodeA  = this->giveNode(1);
-        //Lattice node
         nodeB  = this->giveNode(2);
 
         //Calculate components of distance from reinforcement node to lattice node.
@@ -389,7 +386,6 @@ namespace oofem {
         this->globalCentroid.resize(3);
         for ( int i = 1; i <= 3; i++ ) {
             this->globalCentroid.at(i) = nodeA->giveCoordinate(i);
-            ;
         }
 
         this->geometryFlag = 1;
diff --git a/tests/sm/hangingrotation1.in b/tests/sm/hangingrotation1.in
@@ -0,0 +1,42 @@
+hangingrotation1.out
+#
+Test of rotation constrain in tetrahedron. Only rotations (0.1, 0.2, 0.3) are expected at hanging node.
+NonLinearStatic nmsteps 1 nsteps 2 nmodules 1 contextOutputStep 1000
+nsteps 2 rtolf 1.e-3 stiffMode 2 maxiter 1000 controllmode 1 donotfixload
+errorcheck
+domain 3d
+OutputManager
+#
+ndofman 5 nelem 1 ncrosssect 1 nmat 1 nbc 4 nltf 1 nic 0 nset 5
+node 1 coords 3 0. 0. 0.
+node 2 coords 3 1. 0. 0.
+node 3 coords 3 0. 1. 0.
+node 4 coords 3 0. 0. 1.
+hangingnode 5 coords 3 0.25 0.25 0.25 dofType 6 2 2 2 2 2 2 dofidmask 6 1 2 3 4 5 6
+#
+ltrspace 1 nodes 4 1 2 3 4
+#
+SimpleCS 1 material 1 set 1
+#
+con2dpm 1 d 0 E 30.e9 n 0.15 talpha 0. wf 9.3755e-4 fc 3.e6 ft 1.e6 hp 0.01 yieldtol 1.e-5 asoft 5. stype 1 helem 0.1 kinit 0.3
+#
+BoundaryCondition  1 loadTimeFunction 1 dofs 3 1 2 3 values 3 0.025 -0.05 0.025 set 2
+BoundaryCondition  2 loadTimeFunction 1 dofs 3 1 2 3 values 3 0.025 0.25 -0.175 set 3
+BoundaryCondition  3 loadTimeFunction 1 dofs 3 1 2 3 values 3 -0.275 -0.05 0.125 set 4
+BoundaryCondition  4 loadTimeFunction 1 dofs 3 1 2 3 values 3 0.225 -0.150 0.025 set 5
+#
+PiecewiseLinFunction 1 npoints 2 t 2  0. 1. f(t) 2 0. 1.
+set 1 elements 1 1
+Set 2 nodes 1 1
+Set 3 nodes 1 2
+Set 4 nodes 1 3
+Set 5 nodes 1 4
+
+#%BEGIN_CHECK% tolerance 1e-6
+#NODE tStep 2 number 5 dof 1 unknown d value 0.0
+#NODE tStep 2 number 5 dof 2 unknown d value 0.0
+#NODE tStep 2 number 5 dof 3 unknown d value 0.0
+#NODE tStep 2 number 5 dof 4 unknown d value 0.1
+#NODE tStep 2 number 5 dof 5 unknown d value 0.2
+#NODE tStep 2 number 5 dof 6 unknown d value 0.3
+#%END_CHECK%
diff --git a/tests/sm/hangingrotation2.in b/tests/sm/hangingrotation2.in
@@ -0,0 +1,35 @@
+hangingrotation2.out
+Test of rotation constrain in tetrahedron: Only translations applied
+NonLinearStatic nmsteps 1 nsteps 2 nmodules 1 contextOutputStep 1000
+nsteps 2 rtolf 1.e-3 stiffMode 2 maxiter 1000 controllmode 1 donotfixload
+errorcheck
+domain 3d
+OutputManager tstep_all dofman_all element_all
+#
+ndofman 5 nelem 1 ncrosssect 1 nmat 1 nbc 1 nltf 1 nic 0 nset 2
+node 1 coords 3 0. 0. 0.
+node 2 coords 3 1. 0. 0.
+node 3 coords 3 0. 1. 0.
+node 4 coords 3 0. 0. 1.
+hangingnode 5 coords 3 0.25 0.25 0.25 dofType 6 2 2 2 2 2 2 dofidmask 6 1 2 3 4 5 6
+#
+ltrspace 1 nodes 4 1 2 3 4
+#
+SimpleCS 1 material 1 set 1
+#
+con2dpm 1 d 0 E 30.e9 n 0.15 talpha 0. wf 9.3755e-4 fc 3.e6 ft 1.e6 hp 0.01 yieldtol 1.e-5 asoft 5. stype 1 helem 0.1 kinit 0.3
+#
+BoundaryCondition  1 loadTimeFunction 1 dofs 3 1 2 3 values 3 0.025 -0.05 0.025 set 2
+#
+PiecewiseLinFunction 1 npoints 2 t 2  0. 1. f(t) 2 0. 1.
+set 1 elements 1 1
+Set 2 nodes 4 1 2 3 4 
+
+#%BEGIN_CHECK% tolerance 1e-6
+#NODE tStep 2 number 5 dof 1 unknown d value 0.025
+#NODE tStep 2 number 5 dof 2 unknown d value -0.05
+#NODE tStep 2 number 5 dof 3 unknown d value 0.025
+#NODE tStep 2 number 5 dof 4 unknown d value 0.0
+#NODE tStep 2 number 5 dof 5 unknown d value 0.0
+#NODE tStep 2 number 5 dof 6 unknown d value 0.0
+#%END_CHECK%
diff --git a/tests/sm/hangingrotation3.in b/tests/sm/hangingrotation3.in
@@ -0,0 +1,42 @@
+hangingrotation3.out
+Test of rotation constrain in tetrahedron: Rotations (0.1,0.2,0.3) and translations (0.01, 0.01, 0.01) are expected at hanging node. 
+NonLinearStatic nmsteps 1 nsteps 2 nmodules 1 contextOutputStep 1000
+nsteps 2 rtolf 1.e-3 stiffMode 2 maxiter 1000 controllmode 1 donotfixload
+#StaticStructural nsteps 2 solverType "calm" initialguess 0 rtolf 2.e-3 reqIterations 10 stiffMode 2 maxiter 10 stepLength 1.e-4 minsteplength 1.e-4 Psi 0. nmodules 1
+errorcheck
+domain 3d
+OutputManager
+#
+ndofman 5 nelem 1 ncrosssect 1 nmat 1 nbc 4 nltf 1 nic 0 nset 5
+node 1 coords 3 0. 0. 0.
+node 2 coords 3 1. 0. 0.
+node 3 coords 3 0. 1. 0.
+node 4 coords 3 0. 0. 1.
+hangingnode 5 coords 3 0.25 0.25 0.25 dofType 6 2 2 2 2 2 2 dofidmask 6 1 2 3 4 5 6
+#
+ltrspace 1 nodes 4 1 2 3 4
+#
+SimpleCS 1 material 1 set 1
+#
+con2dpm 1 d 0 E 30.e9 n 0.15 talpha 0. wf 9.3755e-4 fc 3.e6 ft 1.e6 hp 0.01 yieldtol 1.e-5 asoft 5. stype 1 helem 0.1 kinit 0.3
+#
+BoundaryCondition  1 loadTimeFunction 1 dofs 3 1 2 3 values 3 0.035 -0.04 0.035 set 2
+BoundaryCondition  2 loadTimeFunction 1 dofs 3 1 2 3 values 3 0.035 0.26 -0.165 set 3
+BoundaryCondition  3 loadTimeFunction 1 dofs 3 1 2 3 values 3 -0.265 -0.04 0.135 set 4
+BoundaryCondition  4 loadTimeFunction 1 dofs 3 1 2 3 values 3 0.235 -0.140 0.035 set 5
+#
+PiecewiseLinFunction 1 npoints 2 t 2  0. 1. f(t) 2 0. 1.
+set 1 elements 1 1
+Set 2 nodes 1 1 
+Set 3 nodes 1 2
+Set 4 nodes 1 3
+Set 5 nodes 1 4
+
+#%BEGIN_CHECK% tolerance 1e-6
+#NODE tStep 2 number 5 dof 1 unknown d value 0.01
+#NODE tStep 2 number 5 dof 2 unknown d value 0.01
+#NODE tStep 2 number 5 dof 3 unknown d value 0.01
+#NODE tStep 2 number 5 dof 4 unknown d value 0.1
+#NODE tStep 2 number 5 dof 5 unknown d value 0.2
+#NODE tStep 2 number 5 dof 6 unknown d value 0.3
+#%END_CHECK%

Original file line number	Diff line number	Diff line change
`@@ -40,6 +40,7 @@`
`40`	`40`	`#include "feinterpol.h"`
`41`	`41`	`#include "spatiallocalizer.h"`
`42`	`42`	`#include "classfactory.h"`
	`43`	`+#include "elementgeometrytype.h"`
`43`	`44`
`44`	`45`	`namespace oofem {`
`45`	`46`	`REGISTER_DofManager( HangingNode );`
`@@ -169,11 +170,12 @@ void HangingNode ::postInitialize()`
`169`	`170`	`}`
`170`	`171`	`}`
`171`	`172`
	`173`	`+ //Deal with rotations. Works only for linear tetrahedron so far.`
	`174`	`+ if (e->giveGeometryType() == EGT_tetra_1) {`
	`175`	`+`
`172`	`176`	`Dof *dof;`
`173`	`177`	`const int nnodes = e->giveNumberOfNodes();`
`174`	`178`
`175`		`-`
`176`		`- //Deal with rotations.`
`177`	`179`	`// Compute shape function gradients once`
`178`	`180`	`FloatMatrix dNdX;`
`179`	`181`	`fei->evaldNdx(dNdX, lcoords, FEIElementGeometryWrapper(e));`
`@@ -296,6 +298,10 @@ void HangingNode ::postInitialize()`
`296`	`298`	`}`
`297`	`299`	`}`
`298`	`300`
	`301`	`+} else {`
	`302`	`+ OOFEM_WARNING("HangingNode: rotational DOFs only supported for EGT_tetra_1 (linear tetrahedron). Skipping rotational constraint.");`
	`303`	`+}`
	`304`	`+`
`299`	`305`
`300`	`306`	`}`
`301`	`307`	`} // end namespace oofem`