Merge pull request #672 from MassimoPetracca/add-asd-absorbing-boundary

Adding 2D and 3D Absorbing Boundary Elements
zhuminjie · Oct 6, 2021 · 09bda8b · 09bda8b
2 parents 4ff0314 + f63ae97
commit 09bda8b
Show file tree

Hide file tree

Showing 18 changed files with 4,555 additions and 3 deletions.
diff --git a/SRC/Makefile b/SRC/Makefile
@@ -647,6 +647,8 @@ ELE_LIBS   =  $(FE)/element/Element.o \
 	$(FE)/element/catenaryCable/CatenaryCable.o \
 	$(FE)/element/tetrahedron/FourNodeTetrahedron.o \
 	$(FE)/element/absorbentBoundaries/LysmerTriangle.o \
+	$(FE)/element/absorbentBoundaries/ASDAbsorbingBoundary2D.o \
+	$(FE)/element/absorbentBoundaries/ASDAbsorbingBoundary3D.o \
 	$(FE)/element/UWelements/EmbeddedBeamInterfaceL.o \
 	$(FE)/element/UWelements/EmbeddedBeamInterfaceP.o \
 	$(FE)/element/UWelements/EmbeddedEPBeamInterface.o \

diff --git a/SRC/actor/objectBroker/FEM_ObjectBrokerAllClasses.cpp b/SRC/actor/objectBroker/FEM_ObjectBrokerAllClasses.cpp
@@ -307,6 +307,9 @@
 
 #include "CEqElement/ASDEmbeddedNodeElement.h"
 
+#include "absorbentBoundaries/ASDAbsorbingBoundary2D.h"
+#include "absorbentBoundaries/ASDAbsorbingBoundary3D.h"
+
 #include "LinearCrdTransf2d.h"
 #include "LinearCrdTransf3d.h"
 #include "PDeltaCrdTransf2d.h"
@@ -904,6 +907,12 @@ FEM_ObjectBrokerAllClasses::getNewElement(int classTag)
     case ELE_TAG_ASDEmbeddedNodeElement:
       return new ASDEmbeddedNodeElement();
 
+    case ELE_TAG_ASDAbsorbingBoundary2D:
+      return new ASDAbsorbingBoundary2D();
+
+    case ELE_TAG_ASDAbsorbingBoundary3D:
+      return new ASDAbsorbingBoundary3D();
+
     default:
       opserr << "FEM_ObjectBrokerAllClasses::getNewElement - ";
       opserr << " - no Element type exists for class tag " ;

diff --git a/SRC/classTags.h b/SRC/classTags.h
@@ -795,6 +795,8 @@
 #define ELE_TAG_MasonPan3D                    216
 #define ELE_TAG_ASDEmbeddedNodeElement             217  // Massimo Petracca (ASDEA)
 #define ELE_TAG_InertiaTruss              218	//Added by Xiaodong Ji, Yuhao Cheng, Yue Yu
+#define ELE_TAG_ASDAbsorbingBoundary2D    219  // Massimo Petracca (ASDEA)
+#define ELE_TAG_ASDAbsorbingBoundary3D    220  // Massimo Petracca (ASDEA)
 #define ELE_TAG_ExternalElement           99990
 
 

diff --git a/SRC/element/TclElementCommands.cpp b/SRC/element/TclElementCommands.cpp
@@ -166,6 +166,8 @@ extern void *OPS_ASDEmbeddedNodeElement(void); // Massimo Petracca (ASDEA)
 extern void *OPS_ShellANDeS(void);
 extern void *OPS_FourNodeTetrahedron(void);
 extern void *OPS_LysmerTriangle(void);
+extern void *OPS_ASDAbsorbingBoundary2D(void); // Massimo Petracca (ASDEA)
+extern void *OPS_ASDAbsorbingBoundary3D(void); // Massimo Petracca (ASDEA)
 extern void *OPS_TwoNodeLink(void);
 extern void *OPS_LinearElasticSpring(void);
 extern void *OPS_Inerter(void);
@@ -1279,6 +1281,28 @@ TclModelBuilderElementCommand(ClientData clientData, Tcl_Interp *interp,
       }
   }
 
+  else if (strcmp(argv[1], "ASDAbsorbingBoundary2D") == 0) {
+      void *theEle = OPS_ASDAbsorbingBoundary2D();
+      if (theEle != 0) {
+    theElement = (Element*)theEle;
+      } else {
+    opserr<<"tclelementcommand -- unable to create element of type : "
+    <<argv[1]<<endln;
+    return TCL_ERROR;
+      }
+  }
+
+  else if (strcmp(argv[1], "ASDAbsorbingBoundary3D") == 0) {
+      void *theEle = OPS_ASDAbsorbingBoundary3D();
+      if (theEle != 0) {
+    theElement = (Element*)theEle;
+      } else {
+    opserr<<"tclelementcommand -- unable to create element of type : "
+    <<argv[1]<<endln;
+    return TCL_ERROR;
+      }
+  }
+
   else if (strcmp(argv[1], "FourNodeTetrahedron") == 0) {
       void *theEle = OPS_FourNodeTetrahedron();
       if (theEle != 0)