update v0.2.8 --> v0.3.0

Add pressure load and cyclic symm. boudary condicions
avgcampos · Jun 17, 2024 · d3e1488 · d3e1488
1 parent 246cb8c
commit d3e1488
Show file tree

Hide file tree

Showing 19 changed files with 783 additions and 104 deletions.
diff --git a/docs/history.md b/docs/history.md
@@ -51,16 +51,17 @@ Esta nova versão permite a inclusão de módulos e códigos escritos em C/Cytho
 		1. Multi-material
 			- [X] Interface de multi-material (multi E, v, ...)
 		2. Acoplamento multi-físico
-			- [ ] Fluid-Structure Interaction 
-			- [ ] Thermal-Structure Interaction
-			- [ ] Acoustic-Structure Interaction
+			- [ ] Fluid-Structure interaction
+			- [ ] Thermal-Structure interaction
+			- [ ] Acoustic-Structure interaction
     5. Comportamento mecânico e material:
         - [X] Plane stress
         - [X] Plane strain
-        - [X] Solid isotropic
+        - [ ] Solid isotropic
         - [ ] Plate Kirchhoff
         - [ ] Plate Reissner-Mindlin
-        - [ ] Cyclic symmetry boundary
+        - [X] Cyclic symmetry boundary
+		- [X] Pressure load application
         - [ ] Homogeinização/ micro escala (tensor)
 
 8. Validar a análise com a implementação seguintes tipos de elementos isoparametrico:
@@ -78,6 +79,5 @@ Esta nova versão permite a inclusão de módulos e códigos escritos em C/Cytho
 
 10. [X] Desenvolver a saida dos resultados por meio de arquivos .vtk;
 
-## 
-
+##
 Após a nova implementação, foi feito um _merge_ no repositório **myfempy_dev** (repositório oficial para desenvolvimento e testes), e após passar por todos os teste foi feito o upload para o repositório principal **myfempy** do projeto.
diff --git a/docs/index.md b/docs/index.md
@@ -206,69 +206,143 @@ The changelog is available
 
 ``` bash
 /myfempy
-│
-├───core
-│      alglin.py
-│      harmonicforced.py
-│      mesh.py
-│      modallinear.py
-│      solver.py
-│      staticlinear.py
-│      utils.py
-│   
-│
-├───felib
-│   │
-│   ├───elements
-│   │     element.py
-│   │     plane.py
-│   │     platekc.py
-│   │     solid.py
-│   │
-│   ├───geometry
-│   │     geometry.py
-│   │     rectangle.py
-│   │     thickness.py
-│   │   
-│   │
-│   ├───material
-│   │     material.py
-│   │     planestrain.py
-│   │     planestress.py
-│   │     solid.py
-│   │
-│   ├───physic
-│   │     bcstruct.py
-│   │     loadstruct.py
-│   │     structural.py
-│   │
-│   ├───shapes
-│         hexa8.py
-│         line.py
-│         line2.py
-│         quad4.py
-│         shape.py
-│         tetr4.py
-│         tria3.py
-│
-├───io
-│      iocsv.py
-│      iogmsh.py
-│      iovtk.py
-│   
-│
-├───misc
-│      logo.png
-│      logo.txt
-│      path.py
-│      utils.py
-│
-├───plots
-        meshquality.py
-        physics.py
-        plotmesh.py
-        plotxy.py
-        postplot.py
-
+|   __about__.py
+|   __init__.py
+|
++---core
+|   |   utilities.py
+|   |   __init__.py
+|   |
+|   +---elements
+|   |   |   acustic2DWP.py
+|   |   |   beamEB.py
+|   |   |   beamNL.py
+|   |   |   element.py
+|   |   |   fluid1DLF.py
+|   |   |   fluid2DLF.py
+|   |   |   heat2DSS.py
+|   |   |   heat2DTR.py
+|   |   |   plane.py
+|   |   |   plateKC.py
+|   |   |   solid.py
+|   |   |   specials.py
+|   |   |   __init__.py
+|   |
+|   +---geometry
+|   |   |   geometry.py
+|   |   |   rectangle.py
+|   |   |   thickness.py
+|   |   |   __init__.py
+|   |
+|   +---material
+|   |   |   fluidplane.py
+|   |   |   material.py
+|   |   |   microscale.py
+|   |   |   planestrain.py
+|   |   |   planestress.py
+|   |   |   solid.py
+|   |   |   __init__.py
+|   |
+|   +---mesh
+|   |   |   gmsh.py
+|   |   |   legacyquad4.py
+|   |   |   legacytria3.py
+|   |   |   mesh.py
+|   |   |   __init__.py
+|   |
+|   +---physic
+|   |   |   acustic.py
+|   |   |   bcstruct.py
+|   |   |   coupling.py
+|   |   |   fluidflow.py
+|   |   |   loadstruct.py
+|   |   |   structural.py
+|   |   |   thermal.py
+|   |
+|   +---shapes
+|   |   |   hexa20.py
+|   |   |   hexa8.py
+|   |   |   line.py
+|   |   |   line2.py
+|   |   |   line3.py
+|   |   |   quad4.py
+|   |   |   quad4_tasks.c
+|   |   |   quad4_tasks.cp311-win_amd64.pyd
+|   |   |   quad4_tasks.html
+|   |   |   quad4_tasks.pyx
+|   |   |   quad8.py
+|   |   |   shape.py
+|   |   |   tetr10.py
+|   |   |   tetr4.py
+|   |   |   tria3.py
+|   |   |   tria3_tasks.c
+|   |   |   tria3_tasks.cp311-win_amd64.pyd
+|   |   |   tria3_tasks.html
+|   |   |   tria3_tasks.pyx
+|   |   |   tria6.py
+|   |   |   __init__.py
+|   |
+|   +---solver
+|   |   |   acustic.py
+|   |   |   assembler.py
+|   |   |   assemblerfull.py
+|   |   |   assemblerfull_cython_v4.c
+|   |   |   assemblerfull_cython_v4.cp311-win_amd64.pyd
+|   |   |   assemblerfull_cython_v4.html
+|   |   |   assemblerfull_cython_v4.pyx
+|   |   |   assemblerfull_numpy_v1.py
+|   |   |   assemblersymm.py
+|   |   |   assemblersymm_cython_v4.c
+|   |   |   assemblersymm_cython_v4.cp311-win_amd64.pyd
+|   |   |   assemblersymm_cython_v4.html
+|   |   |   assemblersymm_cython_v4.pyx
+|   |   |   assemblersymm_numpy_v1.py
+|   |   |   buckling.py
+|   |   |   cyclicsymm.py
+|   |   |   dynamic.py
+|   |   |   fluid.py
+|   |   |   harmoniclinear.py
+|   |   |   harmonicmodal.py
+|   |   |   modallinear.py
+|   |   |   multiphysic.py
+|   |   |   solver.py
+|   |   |   staticlinear.py
+|   |   |   staticlineariterative.py
+|   |   |   staticnonlinear.py
+|   |   |   thermal.py
+|   |   |   __init__.py
+|
+|
++---io
+|   |   iocsv.py
+|   |   iogmsh.py
+|   |   iovtk.py
+|   |   __init__.py
+|
++---plots
+|   |   meshquality.py
+|   |   physics.py
+|   |   plotmesh.py
+|   |   plotxy.py
+|   |   postplot.py
+|   |   prevplot.py
+|   |   __init__.py
+|
++---setup
+|   |   asi.py
+|   |   fea.py
+|   |   fsi.py
+|   |   model.py
+|   |   physics.py
+|   |   results.py
+|   |   topopt.py
+|   |   tsi.py
+|   |   __init__.py
+|
++---utils
+|   |   logo.png
+|   |   logo.txt
+|   |   utils.py
+|   |   __init__.py
 
 ```
diff --git a/myfempy/__about__.py b/myfempy/__about__.py
@@ -1 +1 @@
-__version__ = "0.2.9"
+__version__ = "0.3.0"
diff --git a/myfempy/__init__.py b/myfempy/__init__.py
@@ -1,6 +1,7 @@
 from __future__ import absolute_import
 
 from .__about__ import __version__
+from .core.solver.cyclicsymm import StaticLinearCyclicSymm
 from .core.solver.harmoniclinear import HarmonicLinear
 from .core.solver.modallinear import ModalLinear
 from .core.solver.staticlinear import StaticLinear
@@ -11,6 +12,7 @@
     "newAnalysis",
     "StaticLinear",
     "StaticLinearIterative",
+    "StaticLinearCyclicSymm",
     "ModalLinear",
     "HarmonicLinear",
     "__version__",

diff --git a/myfempy/core/physic/bcstruct.py b/myfempy/core/physic/bcstruct.py
@@ -20,6 +20,10 @@ def getBCApply(modelinfo, bclist):
                 bcl = bclist[bc_index]
                 bcapp = BoundCondStruct.__BCDispl(modelinfo, bcl)
                 boncdnodeaply = np.append(boncdnodeaply, bcapp, axis=0)
+            elif bclist[bc_index][0] == "csymm":
+                bcl = bclist[bc_index]
+                bcapp = BoundCondStruct.__BCCS(modelinfo, bcl)
+                boncdnodeaply = np.append(boncdnodeaply, bcapp, axis=0)
             else:
                 pass
         boncdnodeaply = boncdnodeaply[1::][::]
@@ -63,3 +67,25 @@ def __BCDispl(modelinfo, bclist):
 
         boncdnodeaply = boncdnodeaply[1::][::]
         return boncdnodeaply
+
+    def __BCCS(modelinfo, bclist):
+        boncdnodeaply = np.zeros((1, 4))
+
+        nodelist = bclist[2:]
+        node_list_bc, dir_fc = get_nodes_from_list(
+            nodelist, modelinfo["coord"], modelinfo["regions"]
+        )
+
+        if bclist[1] == "left":
+            bcdof = 11
+        elif bclist[1] == "right":
+            bcdof = 12
+        else:
+            bcdof = 0
+
+        for j in range(len(node_list_bc)):
+            bcapp = np.array([[int(node_list_bc[j]), bcdof, 0.0, int(bclist[8])]])
+            boncdnodeaply = np.append(boncdnodeaply, bcapp, axis=0)
+
+        boncdnodeaply = boncdnodeaply[1::][::]
+        return boncdnodeaply
diff --git a/myfempy/core/solver/assembler.py b/myfempy/core/solver/assembler.py
@@ -54,3 +54,7 @@ def getConstrains():
     @abstractmethod
     def getDirichletNH():
         pass
+
+    @abstractmethod
+    def getRotationMatrix():
+        pass
diff --git a/myfempy/core/solver/assemblerfull.py b/myfempy/core/solver/assemblerfull.py
@@ -126,14 +126,7 @@ def getConstrains(constrains, nodetot, nodedof):
         elif nodedof == 2:
             for ii in range(ntbc):
                 no = int(constrains[ii, 0])
-                if int(constrains[ii, 1]) == 0:
-                    if constrains[ii, 2] == 0.0:
-                        fixedof[nodedof * no - 2, 0] = nodedof * no - 1
-                        fixedof[nodedof * no - 1, 0] = nodedof * no
-                    else:
-                        constdof[0, nodedof * no - 2, 0] = nodedof * no - 1
-                        constdof[0, nodedof * no - 1, 0] = nodedof * no
-                elif int(constrains[ii, 1]) == 1:
+                if int(constrains[ii, 1]) == 1:
                     if constrains[ii, 2] == 0.0:
                         fixedof[nodedof * no - 2, 0] = nodedof * no - 1
                     else:
@@ -143,20 +136,18 @@ def getConstrains(constrains, nodetot, nodedof):
                         fixedof[nodedof * no - 1, 0] = nodedof * no
                     else:
                         constdof[nodedof * no - 1, 0] = nodedof * no
-
-        elif nodedof == 3:
-            for ii in range(ntbc):
-                no = int(constrains[ii, 0])
-                if int(constrains[ii, 1]) == 0:
+                else:  # int(constrains[ii, 1]) == 0:
                     if constrains[ii, 2] == 0.0:
-                        fixedof[nodedof * no - 3, 0] = nodedof * no - 2
                         fixedof[nodedof * no - 2, 0] = nodedof * no - 1
                         fixedof[nodedof * no - 1, 0] = nodedof * no
                     else:
-                        constdof[nodedof * no - 3, 0] = nodedof * no - 2
-                        constdof[nodedof * no - 2, 0] = nodedof * no - 1
-                        constdof[nodedof * no - 1, 0] = nodedof * no
-                elif int(constrains[ii, 1]) == 1:
+                        constdof[0, nodedof * no - 2, 0] = nodedof * no - 1
+                        constdof[0, nodedof * no - 1, 0] = nodedof * no
+
+        elif nodedof == 3:
+            for ii in range(ntbc):
+                no = int(constrains[ii, 0])
+                if int(constrains[ii, 1]) == 1:
                     if constrains[ii, 2] == 0.0:
                         fixedof[nodedof * no - 3, 0] = nodedof * no - 2
                     else:
@@ -171,6 +162,15 @@ def getConstrains(constrains, nodetot, nodedof):
                         fixedof[nodedof * no - 1, 0] = nodedof * no
                     else:
                         constdof[nodedof * no - 1, 0] = nodedof * no
+                else:  # int(constrains[ii, 1]) == 0:
+                    if constrains[ii, 2] == 0.0:
+                        fixedof[nodedof * no - 3, 0] = nodedof * no - 2
+                        fixedof[nodedof * no - 2, 0] = nodedof * no - 1
+                        fixedof[nodedof * no - 1, 0] = nodedof * no
+                    else:
+                        constdof[nodedof * no - 3, 0] = nodedof * no - 2
+                        constdof[nodedof * no - 2, 0] = nodedof * no - 1
+                        constdof[nodedof * no - 1, 0] = nodedof * no
 
         fixedof = fixedof[np.nonzero(fixedof)]
         fixedof = fixedof - np.ones_like(fixedof)
@@ -236,3 +236,23 @@ def getDirichletNH(constrains, nodetot, nodedof):
             else:
                 pass
         return Uc
+
+    # https://en.wikipedia.org/wiki/Rotation_matrix
+    def getRotationMatrix(node_list, coord, ndof):
+        # Initialize RM as a sparse matrix
+        RM = sparse.lil_matrix((ndof, ndof))
+        for n in range(node_list.shape[0]):
+            nol = int(node_list[n] - 1)
+            RonX = coord[nol, 1]  # - Og[0]
+            RonY = coord[nol, 2]  # - Og[1]
+            Ron = np.sqrt(RonX**2 + RonY**2)
+            S_the = RonY / Ron
+            C_the = RonX / Ron
+            # Assign values to the sparse matrix
+            RM[2 * n, 2 * n] = C_the
+            RM[2 * n, 2 * n + 1] = -S_the
+            RM[2 * n + 1, 2 * n] = S_the
+            RM[2 * n + 1, 2 * n + 1] = C_the
+
+        # Convert to CSR format for more efficient arithmetic and matrix-vector operations
+        return RM.tocsr()
diff --git a/myfempy/core/solver/assemblersymm.py b/myfempy/core/solver/assemblersymm.py
@@ -48,3 +48,6 @@ def getConstrains(constrains, nodetot, nodedof, Uc):
 
     def getDirichletNH(constrains, nodetot, nodedof):
         return AssemblerFULL.getDirichletNH(constrains, nodetot, nodedof)
+
+    def getRotationMatrix(node_list, coord, ndof):
+        return AssemblerFULL.getRotationMatrix(node_list, coord, ndof)