Optional segment length in IntElPoint element. Bondceb02 test case

src/sm/Elements/Interfaces/intelpoint.C

@@ -213,9 +213,14 @@ IntElPoint :: computeGlobalCoordinates(FloatArray &answer, const FloatArray &lco
 double
 IntElPoint :: computeAreaAround(GaussPoint *gp)
 {
+ if ( (this->giveCrossSection()->hasProperty(CS_Thickness)) ) {
+ double thickness = this->giveCrossSection()->give(CS_Thickness, gp);
+ return thickness*this->length;
+ } else {
  // The modeled area/extension around the connected nodes. 
  // Compare with the cs area of a bar. ///@todo replace with cs-property? /JB
  return this->area;
+ }
 }
 
 
@@ -253,6 +258,9 @@ IntElPoint :: initializeFrom(InputRecord *ir)
  this->area = 1.0; // Default area ///@todo Make non-optional? /JB
  IR_GIVE_OPTIONAL_FIELD(ir, this->area, _IFT_IntElPoint_area);
 
+ this->length = 1.0;
+ IR_GIVE_OPTIONAL_FIELD(ir, this->length, _IFT_IntElPoint_length);
+
  this->computeLocalSlipDir(); 
  return IRRT_OK;
 }

src/sm/Elements/Interfaces/intelpoint.h

@@ -46,6 +46,7 @@
 #define _IFT_IntElPoint_refnode "refnode"
 #define _IFT_IntElPoint_normal "normal"
 #define _IFT_IntElPoint_area "area"
+#define _IFT_IntElPoint_length "length"
 //@}
 
 namespace oofem {
@@ -66,6 +67,7 @@ class IntElPoint : public StructuralInterfaceElement
  int referenceNode = 0;
  FloatArrayF<3> normal;
  double area = 0.;
+ double length = 0.;
 
 public:
  IntElPoint(int n, Domain *d);

tests/sm/bondceb02.in

@@ -0,0 +1,56 @@
+bondceb02.out
+Test case for BondCEB interface material - rebar pull-out in 3D. Added by Adam Sciegaj
+StaticStructural nsteps 20 deltat 1.0 rtolf 1.0e-3 MaxIter 50 initialguess 1 nmodules 1
+errorcheck
+#vtkxml tstep_all domain_all primvars 1 1
+#vtkxml tstep_all domain_all ipvars 2 98 99 regionsets 1 3
+#matlab tstep_all integrationpoints internalvars 2 98 99
+domain 3d
+OutputManager tstep_all dofman_all element_all
+ndofman 14 nelem 5 ncrosssect 3 nmat 3 nbc 4 nic 0 nltf 1 nset 6
+node 1 coords 3 -0.1 -0.1 0.1
+node 2 coords 3 -0.1 0.0 0.1
+node 3 coords 3 0.1 0.0 0.1
+node 4 coords 3 0.1 -0.1 0.1
+node 5 coords 3 -0.1 -0.1 -0.1
+node 6 coords 3 -0.1 0.0 -0.1
+node 7 coords 3 0.1 0.0 -0.1
+node 8 coords 3 0.1 -0.1 -0.1
+node 9 coords 3 -0.1 0.1 0.1
+node 10 coords 3 0.1 0.1 0.1
+node 11 coords 3 -0.1 0.1 -0.1
+node 12 coords 3 0.1 0.1 -0.1
+node 13 coords 3 -0.1 0.0 0.1 dofidmask 6 1 2 3 4 5 6
+node 14 coords 3 0.1 0.0 0.1 dofidmask 6 1 2 3 4 5 6
+lspace 1 nodes 8 1 2 3 4 5 6 7 8
+lspace 2 nodes 8 2 9 10 3 6 11 12 7
+libeam3d 3 nodes 2 13 14 refnode 9
+IntELPoint 4 nodes 2 2 13 normal 3 0 0 1 length 0.1
+IntElPoint 5 nodes 2 3 14 normal 3 0 0 1 length 0.1
+Set 1 elements 2 1 2
+Set 2 elements 1 3
+Set 3 elements 2 4 5
+Set 4 nodes 1 1
+Set 5 nodes 2 1 4
+Set 6 nodes 1 14
+SimpleCS 1 material 1 set 1
+SimpleCS 2 area 3.1416e-4 material 2 set 2
+InterfaceCS 3 thickness 6.28319e-2 material 3 set 3
+IsoLE 1 d 1.0 E 3.36e10 n 0.2 tAlpha 1.0
+IsoLE 2 d 1.0 E 2e11 n 0.3 tAlpha 1.0
+bondceb 3 kn 6e12 ks 6.135e10 s1 0.001 s2 0.002 s3 0.0065 taumax 1.541e7 tauf 6.164e6
+BoundaryCondition 1 loadTimeFunction 1 dofs 3 1 2 3 values 3 0 0 0 set 4
+BoundaryCondition 2 loadTimeFunction 1 dofs 2 2 3 values 2 0 0 set 5
+BoundaryCondition 3 loadTimeFunction 1 dofs 3 4 5 6 values 3 0 0 0 set 2
+BoundaryCondition 4 loadTimeFunction 1 dofs 1 1 values 1 8e-3 set 6
+PiecewiseLinFunction 1 t 2 0.0 20.0 f(t) 2 0.0 1.0
+
+#%BEGIN_CHECK% tolerance 1.e1
+##Bond stress values check
+#ELEMENT tStep 1 number 5 gp 1 keyword 99 component 2 value 9.01726573e+06
+#ELEMENT tStep 5 number 5 gp 1 keyword 99 component 2 value 1.54100000e+07
+#ELEMENT tStep 8 number 5 gp 1 keyword 99 component 2 value 1.342479552e+07
+#ELEMENT tStep 15 number 5 gp 1 keyword 99 component 2 value 7.45817170e+06
+#ELEMENT tStep 19 number 5 gp 1 keyword 99 component 2 value 6.16400000e+06
+##
+#%END_CHECK%