kib-lap 0.5__cp313-cp313-win_amd64.whl → 0.7.7__cp313-cp313-win_amd64.whl

KIB_LAP/FACHWERKEBEN/main.py

@@ -0,0 +1,27 @@
+from Iteration import IterationClass
+from Plotting import PlottingStructure
+
+
+Iteration = IterationClass()
+try:
+    Iteration.Summarize()
+except:
+    Iteration.SummarizeLinear()
+
+plotter = PlottingStructure(
+    Iteration.UG_FINAL,
+    Iteration.FI_FINAL,
+    Iteration.EXTFORCES,
+    Iteration.MBRFORCES,
+    Iteration.Inp.nodes,
+    Iteration.Inp.members,
+    Iteration.Mat.P0,
+    Iteration.Inp.restrainedIndex,
+    _reactionsFlag=True,
+)
+
+plotter.Kspring_diag = Iteration.Kspring_diag  # (nDoF,)
+# optional (noch besser, damit nur echte Federn geplottet werden)
+# plotter.spring_dofs = Iteration.spring_dofs
+
+plotter.create_structure_plot()

KIB_LAP/Plattentragwerke/PlateBendingKirchhoff.py

@@ -33,10 +33,10 @@ class PlateBendingKirchhoffClass:
         loading="Regular",
         support="hhhh",
         reihen=8,
-        x_s_positions = [],
-        I_s_values= [],
-        y_s_positions =  [],
-        I_t_values = []
+        x_s_positions=[],
+        I_s_values=[],
+        y_s_positions=[],
+        I_t_values=[],
     ):
         """
         Initialisierung der Klasse PlateBendingKirchhoff.
@@ -97,15 +97,12 @@ class PlateBendingKirchhoffClass:
 
         self.load = np.zeros((self.reihen**2))
 
-
         self.x_s_array = np.array(x_s_positions, dtype=np.float64)
         self.I_s_array = np.array(I_s_values, dtype=np.float64)
 
-
         self.y_s_array = np.array(y_s_positions, dtype=np.float64)
         self.I_t_array = np.array(I_t_values, dtype=np.float64)
 
-
     def CalculateAll(self):
         """
         Führt alle Berechnungsschritte aus.
@@ -142,7 +139,7 @@ class PlateBendingKirchhoffClass:
             self.x_s_array,
             self.I_s_array,
             self.y_s_array,
-            self.I_t_array
+            self.I_t_array,
         )
 
         self.mat = np.array(self.mat)  # Konvertierung in ein NumPy-Array
@@ -779,49 +776,59 @@ class PlateBendingKirchhoffClass:
         plt.close()
 
 
-
 a = 15
 b = 6
 
 Plate = PlateBendingKirchhoffClass(
-    35000, 0.30, a, b, 1, 0.5, 0.5, 0.5, 0.5, 0.0, 0, 0, n_inte=50, loading="Liste", support = "hhhh",reihen = 8)
+    35000,
+    0.30,
+    a,
+    b,
+    1,
+    0.5,
+    0.5,
+    0.5,
+    0.5,
+    0.0,
+    0,
+    0,
+    n_inte=50,
+    loading="Liste",
+    support="hhhh",
+    reihen=8,
+)
 
 
 btire = 0.37
-p_tire = 120/(4*btire**2) / 1000
+p_tire = 120 / (4 * btire**2) / 1000
 
 print("tire pressure [MN/m²]: ", p_tire)
 
 # Daten als Liste von Listen
 new_data = [
-    ["No.","x0[m]","x1[m]","y0[m]","y1[m]","p0[MN/m**2]"],
-    [1,3+2-btire,3+2+btire,0.5,0.5+2*btire,p_tire *1.5],
-    [2,3+3.2-btire,3+3.2+btire,0.5,0.5+2*btire,p_tire *1.5],
-    [3,3+3.2-btire,3+3.2+btire,2.5,2.5+2*btire,p_tire *1.5],
-    [4,3+2-btire,3+2+btire,2.5,2.5+2*btire,p_tire *1.5], 
-
-    [5,0,15,0,3,6.5/1000*1.5],
-
-    [6,0,15,0,6,2.5/1000*1.5],
-    [7,0,15,0,6,0.35*25/1000*1.35],
-
+    ["No.", "x0[m]", "x1[m]", "y0[m]", "y1[m]", "p0[MN/m**2]"],
+    [1, 3 + 2 - btire, 3 + 2 + btire, 0.5, 0.5 + 2 * btire, p_tire * 1.5],
+    [2, 3 + 3.2 - btire, 3 + 3.2 + btire, 0.5, 0.5 + 2 * btire, p_tire * 1.5],
+    [3, 3 + 3.2 - btire, 3 + 3.2 + btire, 2.5, 2.5 + 2 * btire, p_tire * 1.5],
+    [4, 3 + 2 - btire, 3 + 2 + btire, 2.5, 2.5 + 2 * btire, p_tire * 1.5],
+    [5, 0, 15, 0, 3, 6.5 / 1000 * 1.5],
+    [6, 0, 15, 0, 6, 2.5 / 1000 * 1.5],
+    [7, 0, 15, 0, 6, 0.35 * 25 / 1000 * 1.35],
     # [6,2-0.2,2+0.2,3.5,3.9,100/0.4**2/1000],
     # [7,3.2-0.2,3.2+0.2,3.5,3.9,100/0.4**2/1000],
     # [8,3.2-0.2,3.2+0.2,5.5,5.9,100/0.4**2/1000],
     # [9,2-0.2,2+0.2,5.5,5.9,100/0.4**2/1000],
-
     # [10,2-0.2,2+0.2,6.5,6.9,100/0.4**2/1000],
     # [11,3.2-0.2,3.2+0.2,6.5,6.9,100/0.4**2/1000],
     # [12,3.2-0.2,3.2+0.2,8.5,8.9,100/0.4**2/1000],
     # [13,2-0.2,2+0.2,8.5,8.9,100/0.4**2/1000]
-
 ]
 
 # Konvertiere die Daten in ein DataFrame
 df = pd.DataFrame(new_data[1:], columns=new_data[0])
 
 # Dateipfad zur CSV-Datei
-file_path = 'Loading/Constant_Loading.csv'
+file_path = "Loading/Constant_Loading.csv"
 
 # Schreibe das DataFrame in eine CSV-Datei
 df.to_csv(file_path, index=False)
@@ -829,9 +836,9 @@ df.to_csv(file_path, index=False)
 
 Plate.CalculateAll()
 
-Plate.SolutionPointDisp(0.5,0.5)
-Plate.SolutionPointMomentx(0.5,0.5)
-Plate.SolutionPointMomenty(0.5,0.5)
+Plate.SolutionPointDisp(0.5, 0.5)
+Plate.SolutionPointMomentx(0.5, 0.5)
+Plate.SolutionPointMomenty(0.5, 0.5)
 
 
 Plate.PlotLoad()
@@ -840,4 +847,4 @@ Plate.PlotMomentGrid()
 print("SHEAR FORCES")
 Plate.PlotShearForceGrid()
 print("DRILL MOMENTS")
-Plate.PlotTorsionalMomentGrid()
+Plate.PlotTorsionalMomentGrid()

KIB_LAP/STABRAUM/InputData.py

@@ -20,6 +20,7 @@ class Input:
         self.ElementLoads()
         self.MaterialData()
         self.CrossSectionData()
+        self.Springs()
 
     def MaterialData(self):
         try:
@@ -47,7 +48,7 @@ class Input:
             df_numeric = df[numeric_columns]
 
             self.nodes = df_numeric
-
+            self.nDoF = (len(self.nodes))*7
             print("1. 🟢 Vertices.csv imported")
         else:
             self.forceLocationData = []
@@ -60,7 +61,7 @@ class Input:
             numeric_columns = ['na', 'ne',"cs"]
             df_numeric = df[numeric_columns]
             self.members = df_numeric
-            self.nDoF = (len(self.members)+1)*7
+            
             print("2. 🟢 Edges.csv imported")
         else:
             print("2. 🛑 STOP: Edges.csv not found")
@@ -72,6 +73,16 @@ class Input:
         df_numeric = df[numeric_columns]
         self.RestraintData = df_numeric
 
+    def Springs(self):
+        if not glob("data/Springs.csv"):
+            self.SpringsData = pd.DataFrame(columns=["node_a","node_e","dof","cp[MN]"])
+            print("⚠️ Springs.csv not found (optional)")
+            return
+
+        df = pd.read_csv("data/Springs.csv")
+        self.SpringsData = df[["node_a","node_e","dof","cp[MN]"]].copy()
+        print("🟢 Springs.csv imported")
+
     def NodalForces(self):
         df = pd.read_csv("data/NodalForces.csv")
         self.NodalForces = df

KIB_LAP/STABRAUM/Output_Data.py

@@ -0,0 +1,61 @@
+import pandas as pd
+import numpy as np
+
+
+class OutputData:
+    def __init__(self):
+        print("Init OutputData")
+
+    def support_reactions_from_springs_table(self, res, u0_default=0.0):
+        # --- robuster Typ-/Attribut-Check ---
+        needed = ("Inp", "GesMat", "u_ges", "FGes")
+        missing = [a for a in needed if not hasattr(res, a)]
+        if missing:
+            raise TypeError(
+                "Du musst ein AnalysisResults-Objekt übergeben (res = mainloop().run()). "
+                f"Übergeben wurde: {type(res)}; fehlende Attribute: {missing}"
+            )
+
+        dfR = res.Inp.RestraintData.copy()
+
+        # u als 1D
+        u = np.asarray(res.u_ges, dtype=float).reshape(-1)
+
+        # DOF Mapping (7 dof/node)
+        dof_name = {
+            0: ("Fx", "MN"),
+            1: ("Fy", "MN"),
+            2: ("Mz", "MNm"),
+            3: ("Fz", "MN"),
+            4: ("My", "MNm"),
+            5: ("Mx", "MNm"),
+            6: ("W", "-"),
+        }
+
+        rows = []
+        for _, row in dfR.iterrows():
+            node = int(row["Node"])
+            dof = int(row["Dof"])
+            k = float(row["Cp[MN/m]/[MNm/m]"])
+
+            gdof = 7 * (node - 1) + dof
+            ui = float(u[gdof])
+            u0 = float(u0_default)
+
+            Ri = k * (ui - u0)
+
+            name, unit = dof_name.get(dof, (f"DOF{dof}", ""))
+
+            rows.append(
+                {
+                    "Node": node,
+                    "DOF": dof,
+                    "Type": name,
+                    "k": k,
+                    "u": ui,
+                    "R_spring": Ri,
+                    "Unit": unit,
+                }
+            )
+
+        return pd.DataFrame(rows).sort_values(["Node", "DOF"]).reset_index(drop=True)