Adding latticelink elements and materials

doc/elementlibmanual/elementlibmanual.tex

@@ -373,6 +373,24 @@ \subsubsection{Lattice2dboundary}
 
 The theory pf lattice2dbounday is described in the paper ``P. Grassl and M. Jir\'{a}sek. Meso-scale approach to modelling the fracture process zone of concrete subjected to uniaxial tension. International Journal of Solids and Structures. Volume 47, Issues 7-8, pp. 957-968, 2010.''
 
+\subsubsection{Lattice3d}
+Represents two-node 3d lattice element. Each node has six degrees of freedom.
+The theory of lattice3d is described in the paper ``P. Grassl, J. Bolander. Three-Dimensional Network Model for Coupling of Fracture and Mass Transport in Quasi-Brittle Geomaterials, Materials, 9, 782, 2016.''
+
+\subsubsection{Lattice3dboundary}
+Represents three-noded 3d lattice element for boundary of 3d periodic cells. The first two nodes have 6 degrees of freedom as for the element lattice3d. The third node is used to control the loading of the periodic cell using three normal (xx, yy and zz) and three shear strain (yz, xz, xy) components.
+
+The theory of lattice3dboundary is described in the paper ``I. Athanasiadis, S. Wheeler and P. Grassl. Hydro-mechanical network modelling of particulate composites, International Journal of Solids and Structures. Vol. 130-131, Pages 49-60, 2018.''
+
+\subsubsection{Latticelink3d}
+Represents two-node 3d link element connecting 3d beam and 3d lattice elements. Each node has six degrees of freedom.
+The theory of latticelink3d is described in the paper ``P. Grassl and A. Antonelli. 3D network modelling of fracture processes in fibre-reinforced geomaterials, International Journal of Solids and Structures, vol. 156-157, Pages 234-242, 2019.''
+
+\subsubsection{Latticelink3dboundary}
+Represents three-node 3d boundary link element connecting 3d beam and 3d lattice elements. The first two nodes have the same meaning as for latticelink3d. The third node is used to control the loading of the periodic cell using three normal (xx, yy and zz) and three shear strain (yz, xz, xy) components.
+
+The theory of latticelink3dboundary is described in the paper ``P. Grassl and A. Antonelli. 3D network modelling of fracture processes in fibre-reinforced geomaterials, International Journal of Solids and Structures, vol. 156-157, Pages 234-242, 2019.''
+
 %-----------------------------------------------------------------------------------------------
 \clearpage
 \subsection{Plane Stress Elements}

src/sm/CMakeLists.txt

@@ -270,7 +270,8 @@ set (sm_material
  Materials/latticeviscoelastic.C
  Materials/latticedamageviscoelastic.C
  Materials/latticebondplasticity.C
- Materials/latticeplasticitydamage.C 
+ Materials/latticeplasticitydamage.C
+ Materials/latticeslip.C
  )
 
 if (USE_PYTHON)

src/sm/Elements/lattice3dboundary.C

@@ -200,7 +200,7 @@ Lattice3dBoundary :: computeStiffnessMatrix(FloatMatrix &answer, MatResponseMode
  computeGeometryProperties();
  }
 
- double volume = this->computeVolumeAround( integrationRulesArray [ 0 ]->getIntegrationPoint(0) );
+ double volume = this->computeVolumeAround(integrationRulesArray [ 0 ]->getIntegrationPoint(0) );
 
  this->computeBmatrixAt(integrationRulesArray [ 0 ]->getIntegrationPoint(0), bj);
  this->computeConstitutiveMatrixAt(d, rMode, integrationRulesArray [ 0 ]->getIntegrationPoint(0), tStep);
@@ -212,10 +212,10 @@ Lattice3dBoundary :: computeStiffnessMatrix(FloatMatrix &answer, MatResponseMode
  bjt.beTranspositionOf(bj);
  answerTemp.beProductOf(bjt, dbj);
 
- answer.resize( computeNumberOfDofs(), computeNumberOfDofs() );
+ answer.resize(computeNumberOfDofs(), computeNumberOfDofs() );
  answer.zero();
 
- answerHelp.resize( computeNumberOfDofs(), computeNumberOfDofs() );
+ answerHelp.resize(computeNumberOfDofs(), computeNumberOfDofs() );
  answerHelp.zero();
 
  for ( int m = 1; m <= 12; m++ ) {
@@ -244,13 +244,13 @@ Lattice3dBoundary :: computeStiffnessMatrix(FloatMatrix &answer, MatResponseMode
  IntArray projectionComponentNodeOne(3);
  projectionComponentNodeOne.zero();
  if ( location.at(1) != 0 ) {
- giveSwitches( projectionComponentNodeOne, location.at(1) );
+ giveSwitches(projectionComponentNodeOne, location.at(1) );
  }
 
  IntArray projectionComponentNodeTwo(3);
  projectionComponentNodeTwo.zero();
  if ( location.at(2) != 0 ) {
- giveSwitches( projectionComponentNodeTwo, location.at(2) );
+ giveSwitches(projectionComponentNodeTwo, location.at(2) );
  }
 
  for ( int k = 1; k <= 12; k++ ) {
@@ -317,12 +317,12 @@ Lattice3dBoundary :: giveVTKCoordinates(int nodeNumber, FloatArray &coords) {
  if ( nodeNumber == 1 ) {
  node = this->giveNode(1);
  if ( location.at(1) != 0 ) {
- giveSwitches( projectionComponent, location.at(1) );
+ giveSwitches(projectionComponent, location.at(1) );
  }
  } else if ( nodeNumber == 2 ) {
  node = this->giveNode(2);
  if ( location.at(2) != 0 ) {
- giveSwitches( projectionComponent, location.at(2) );
+ giveSwitches(projectionComponent, location.at(2) );
  }
  } else {
  OOFEM_ERROR("wrong element used in the vtk module");
@@ -359,14 +359,14 @@ Lattice3dBoundary :: computeStrainVector(FloatArray &answer, GaussPoint *gp, Tim
  IntArray projectionComponentNodeOne(3);
  projectionComponentNodeOne.zero();
  if ( location.at(1) != 0 ) {
- giveSwitches( projectionComponentNodeOne, location.at(1) );
+ giveSwitches(projectionComponentNodeOne, location.at(1) );
  }
 
  IntArray projectionComponentNodeTwo(3);
  projectionComponentNodeTwo.zero();
 
  if ( location.at(2) != 0 ) {
- giveSwitches( projectionComponentNodeTwo, location.at(2) );
+ giveSwitches(projectionComponentNodeTwo, location.at(2) );
  }
 
  FloatMatrix rotationMatrix;
@@ -484,18 +484,18 @@ Lattice3dBoundary :: giveInternalForcesVector(FloatArray &answer, TimeStep *tSte
 
  bt.beTranspositionOf(b);
  if ( useUpdatedGpRecord == 1 ) {
- TotalStressVector = ( ( StructuralMaterialStatus * ) mat->giveStatus( integrationRulesArray [ 0 ]->getIntegrationPoint(0) ) )
+ TotalStressVector = ( ( StructuralMaterialStatus * ) mat->giveStatus(integrationRulesArray [ 0 ]->getIntegrationPoint(0) ) )
  ->giveStressVector();
  } else
  if ( !this->isActivated(tStep) ) {
- strain.resize( StructuralMaterial :: giveSizeOfVoigtSymVector( integrationRulesArray [ 0 ]->getIntegrationPoint(0)->giveMaterialMode() ) );
+ strain.resize(StructuralMaterial :: giveSizeOfVoigtSymVector(integrationRulesArray [ 0 ]->getIntegrationPoint(0)->giveMaterialMode() ) );
  strain.zero();
  }
  this->computeStrainVector(strain, integrationRulesArray [ 0 ]->getIntegrationPoint(0), tStep);
 
  this->computeStressVector(TotalStressVector, strain, integrationRulesArray [ 0 ]->getIntegrationPoint(0), tStep);
 
- dV = this->computeVolumeAround( integrationRulesArray [ 0 ]->getIntegrationPoint(0) );
+ dV = this->computeVolumeAround(integrationRulesArray [ 0 ]->getIntegrationPoint(0) );
  bs.beProductOf(bt, TotalStressVector);
  bs.times(dV);
 
@@ -510,14 +510,14 @@ Lattice3dBoundary :: giveInternalForcesVector(FloatArray &answer, TimeStep *tSte
  IntArray projectionComponentNodeOne(3);
  projectionComponentNodeOne.zero();
  if ( location.at(1) != 0 ) {
- giveSwitches( projectionComponentNodeOne, location.at(1) );
+ giveSwitches(projectionComponentNodeOne, location.at(1) );
  }
 
  IntArray projectionComponentNodeTwo(3);
  projectionComponentNodeTwo.zero();
 
  if ( location.at(2) != 0 ) {
- giveSwitches( projectionComponentNodeTwo, location.at(2) );
+ giveSwitches(projectionComponentNodeTwo, location.at(2) );
  }
 
  //Normal stresses
@@ -579,13 +579,13 @@ Lattice3dBoundary :: computeGeometryProperties()
  IntArray projectionComponentNodeOne(3);
  projectionComponentNodeOne.zero();
  if ( location.at(1) != 0 ) {
- giveSwitches( projectionComponentNodeOne, location.at(1) );
+ giveSwitches(projectionComponentNodeOne, location.at(1) );
  }
 
  IntArray projectionComponentNodeTwo(3);
  projectionComponentNodeTwo.zero();
  if ( location.at(2) != 0 ) {
- giveSwitches( projectionComponentNodeTwo, location.at(2) );
+ giveSwitches(projectionComponentNodeTwo, location.at(2) );
  }
 
  for ( int i = 0; i < 3; i++ ) {
@@ -603,7 +603,7 @@ Lattice3dBoundary :: computeGeometryProperties()
  this->normal.at(i + 1) = coordsB.at(i + 1) - coordsA.at(i + 1);
  }
 
- this->length = sqrt( pow(this->normal.at(1), 2.) + pow(this->normal.at(2), 2.) + pow(this->normal.at(3), 2.) );
+ this->length = sqrt(pow(this->normal.at(1), 2.) + pow(this->normal.at(2), 2.) + pow(this->normal.at(3), 2.) );
 
  for ( int i = 0; i < 3; i++ ) {
  this->normal.at(i + 1) /= length;
@@ -688,7 +688,7 @@ void Lattice3dBoundary :: drawRawGeometry(oofegGraphicContext &gc, TimeStep *tSt
  }
 
  EASValsSetLineWidth(OOFEG_RAW_GEOMETRY_WIDTH);
- EASValsSetColor( gc.getActiveCrackColor() );
+ EASValsSetColor(gc.getActiveCrackColor() );
  EASValsSetLayer(OOFEG_RAW_GEOMETRY_LAYER);
 
 
@@ -701,14 +701,14 @@ void Lattice3dBoundary :: drawRawGeometry(oofegGraphicContext &gc, TimeStep *tSt
  IntArray projectionComponentNodeOne(3);
  projectionComponentNodeOne.zero();
  if ( location.at(1) != 0 ) {
- giveSwitches( projectionComponentNodeOne, location.at(1) );
+ giveSwitches(projectionComponentNodeOne, location.at(1) );
  }
 
  IntArray projectionComponentNodeTwo(3);
  projectionComponentNodeTwo.zero();
 
  if ( location.at(2) != 0 ) {
- giveSwitches( projectionComponentNodeTwo, location.at(2) );
+ giveSwitches(projectionComponentNodeTwo, location.at(2) );
  }
 
  p [ 0 ].x = ( FPNum ) this->giveNode(1)->giveCoordinate(1) + projectionComponentNodeOne.at(1) * specimenDimension.at(1);
@@ -781,7 +781,7 @@ void Lattice3dBoundary :: drawDeformedGeometry(oofegGraphicContext &gc, TimeStep
  WCRec p [ 2 ]; /* points */
 
  EASValsSetLineWidth(OOFEG_DEFORMED_GEOMETRY_WIDTH);
- EASValsSetColor( gc.getDeformedElementColor() );
+ EASValsSetColor(gc.getDeformedElementColor() );
  EASValsSetLayer(OOFEG_DEFORMED_GEOMETRY_LAYER);
 
  FloatArray specimenDimension(3);
@@ -816,13 +816,13 @@ void Lattice3dBoundary :: drawDeformedGeometry(oofegGraphicContext &gc, TimeStep
  IntArray projectionComponentNodeOne(3);
  projectionComponentNodeOne.zero();
  if ( location.at(1) != 0 ) {
- giveSwitches( projectionComponentNodeOne, location.at(1) );
+ giveSwitches(projectionComponentNodeOne, location.at(1) );
  }
 
  IntArray projectionComponentNodeTwo(3);
  projectionComponentNodeTwo.zero();
  if ( location.at(2) != 0 ) {
- giveSwitches( projectionComponentNodeTwo, location.at(2) );
+ giveSwitches(projectionComponentNodeTwo, location.at(2) );
  }
 
 

src/sm/Elements/lattice3dboundary.h

@@ -10,7 +10,7 @@
  *
  * OOFEM : Object Oriented Finite Element Code
  *
- * Copyright (C) 1993 - 2013 Borek Patzak
+ * Copyright (C) 1993 - 2019 Borek Patzak
  *
  *
  *
@@ -33,14 +33,14 @@
  */
 
 #ifndef lattice3dboundary_h
- #define lattice3dboundary_h
+#define lattice3dboundary_h
 
- #include "lattice3d.h"
+#include "lattice3d.h"
 
-///@name Input fields for Lattice2d
+///@name Input fields for Lattice3d
 //@{
- #define _IFT_Lattice3dBoundary_Name "lattice3dboundary"
- #define _IFT_Lattice3dBoundary_location "location"
+#define _IFT_Lattice3dBoundary_Name "lattice3dboundary"
+#define _IFT_Lattice3dBoundary_location "location"
 //@}
 
 namespace oofem {
@@ -80,12 +80,12 @@ class Lattice3dBoundary : public Lattice3d
 
  void giveVTKCoordinates(int nodeNumber, FloatArray &coords);
 
- #ifdef __OOFEG
+#ifdef __OOFEG
  virtual void drawYourself(oofegGraphicContext &context, TimeStep *tStep);
  virtual void drawRawGeometry(oofegGraphicContext &, TimeStep *tStep);
  virtual void drawRawCrossSections(oofegGraphicContext &, TimeStep *tStep);
- virtual void drawDeformedGeometry(oofegGraphicContext &, TimeStep * tStep, UnknownType);
- #endif
+ virtual void drawDeformedGeometry(oofegGraphicContext &, TimeStep *tStep, UnknownType);
+#endif
 
 protected:
  void computeBmatrixAt(GaussPoint *, FloatMatrix &, int = 1, int = ALL_STRAINS) override;
@@ -95,7 +95,5 @@ class Lattice3dBoundary : public Lattice3d
  void giveSwitches(IntArray &answer, int location);
  void computeStrainVector(FloatArray &answer, GaussPoint *gp, TimeStep *stepN) override;
 };
-
-
 } // end namespace oofem
 #endif