ourocode 2.1.4__tar.gz → 2.1.6__tar.gz

{ourocode-2.1.4/ourocode.egg-info → ourocode-2.1.6}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: ourocode
-Version: 2.1.4
+Version: 2.1.6
 Summary: Ceci est un catalogue de fonction permettant une utilisation rapide pour la réalisation de note de calcul personnalisée.
 Author-email: Anthony PARISOT <contact@ourea-structure.fr>
 License-Expression: Apache-2.0

{ourocode-2.1.4 → ourocode-2.1.6}/ourocode/eurocode/core/model_result.py

@@ -834,9 +834,10 @@ class Model_result(Projet):
             Un glisseur retournera une force nulle dans la direction libre.
         """
         reaction = {}
-        supports = self._model_generator.get_all_supports()
-        for support_id, support in supports.items():
-            node_id = support["Noeud"]
-            reaction[support_id] = self._get_node_reaction(node_id, combination)
+        all_supports = self._model_generator.get_all_supports()
+        for support_type, supports in all_supports.items():
+            for support_id, support in supports.items():
+                node_id = support["Noeud"]
+                reaction[support_id] = self._get_node_reaction(node_id, combination)
         return reaction
 

{ourocode-2.1.4 → ourocode-2.1.6}/ourocode/eurocode/ec1/neige.py

@@ -455,7 +455,7 @@ class Neige(Batiment):
 		"""Retourne la charge de neige en kN/m sur un obstacle ou une barre à neige
 
 		Args:
-			Sn (float): la charge de neige en kN/m²
+			Sn (float): la charge de neige normale au sol en kN/m²
 			entraxe (int): entraxe horizontal entre deux obstacle à la neige en m
 			alpha (float): angle du versant considéré.
 
@@ -470,6 +470,46 @@ class Neige(Batiment):
 			return F_s
 		return val()
 
+	def poussee_hydrostatique(self, h: si.m):
+		"""Retourne la poussée horizontale de neige "hydrostatique" selon l'EN 1991-1-3 §6.2.
+
+		Modélise l'action horizontale due à une accumulation de neige contre un obstacle
+		vertical (parapet, écran, mur...) par une charge triangulaire :
+		pression nulle au sommet de la couche de neige et maximale à la base.
+
+		Formules :
+
+			- hauteur de neige : h_n = S_n / γ
+			- pression maximale à la base : q_h = γ × h_n
+			- résultante horizontale par mètre linéaire : F_hl = 0.5 × q_h × h
+			- hauteur équivalente d'application de la résultante: h_eq = 1/3 * h_n
+
+		Args:
+			h (float): Hauteur de la paroi verticale en m.
+
+		Returns:
+			tuple: (latex_string, dict) avec :
+				- "q_h" : pression horizontale maximale à la base en kN/m².
+				- "F_hl" : résultante horizontale par unité de largeur en kN/m.
+				- "h_neige" : hauteur de neige en m.
+		"""
+		h = h * si.m
+		gamma = 2 * si.kN/si.m**3
+		Sk_alt = self.Sk_alt[1]
+		C_e = self.C_e
+		C_t = self.C_t
+
+		@handcalc(override="short", precision=2, jupyter_display=self.JUPYTER_DISPLAY, left="\\[", right="\\]")
+		def val():
+			mu_i = min(max(gamma * h / Sk_alt, 0.8), 2)
+			S_n = Sk_alt * mu_i * C_e * C_t
+			h_neige = S_n / gamma # hauteur de neige
+			h_eq = h_neige * 1/3 # hauteur équivalente d'application de la résultante
+			q_h = gamma * h_neige # pression maximale à la base
+			F_hl = 0.5 * q_h * h_neige # résultante horizontale par mètre linéaire
+			return q_h, F_hl, h_neige, h_eq
+		latex, (q_h, F_hl, h_neige, h_eq) = val()
+		return latex, {"q_h": q_h, "F_hl": F_hl, "h_neige": h_neige, "h_equivalent": h_eq}
 
 
 if __name__ == "__main__":

{ourocode-2.1.4 → ourocode-2.1.6}/ourocode/eurocode/ec5/cvt.py

@@ -379,11 +379,11 @@ class MOB(Batiment):
         if couturage_ext_rive:
             wall['Couturage ext. rive'] = couturage_ext_rive * si.mm
             wall['Couturage ext. intermédiaire'] = couturage_ext_inter * si.mm
-            if Ff_Rk_ext is None:
+            if Ff_Rk_ext is None or Ff_Rk_ext == 0 * si.N:
                 wall['Ff,Rk ext'] = self._data_connecteurs['Extérieur']['Ff,Rk']
             else:
                 wall['Ff,Rk ext'] = Ff_Rk_ext * si.N
-            if Kser_ext is None:
+            if Kser_ext is None or Kser_ext == 0 * si.N/si.mm:
                 wall['Kser ext'] = self._data_connecteurs['Extérieur']['Kser']
             else:
                 wall['Kser ext'] = Kser_ext * si.N / si.mm
@@ -391,11 +391,11 @@ class MOB(Batiment):
         if couturage_int_rive:
             wall['Couturage int. rive'] = couturage_int_rive * si.mm
             wall['Couturage int. intermédiaire'] = couturage_int_inter * si.mm
-            if Ff_Rk_int is None:
+            if Ff_Rk_int is None or Ff_Rk_int == 0 * si.N:
                 wall['Ff,Rk int'] = self._data_connecteurs['Intérieur']['Ff,Rk']
             else:
                 wall['Ff,Rk int'] = Ff_Rk_int * si.N
-            if Kser_int is None:
+            if Kser_int is None or Kser_int == 0 * si.N/si.mm:
                 wall['Kser int'] = self._data_connecteurs['Intérieur']['Kser']
             else:
                 wall['Kser int'] = Kser_int * si.N / si.mm
@@ -529,6 +529,10 @@ class MOB(Batiment):
             k_ser_r (float): raideur du panneau.
             k_ser_i (float): raideur du panneau.
         """ 
+        if b == 0 * si.mm or h == 0 * si.mm:
+            raise ValueError("Dimensions de panneau invalides: b et h doivent être non nulles.")
+        if k_ser_r == 0 * si.N / si.mm**2:
+            raise ValueError("Raideur linéique de rive nulle (k_ser_r): vérifier Kser et les couturages de rive.")
         K_p = (2 / (k_ser_r * (b + h/3) + k_ser_i * h/6) + 2*h**2 / (k_ser_r * b**2 * (h + b/3)))**-1
         return K_p
 
@@ -556,11 +560,11 @@ class MOB(Batiment):
             for sys_name, sys_wall in sys.items():
                 h_sys = sys_wall['Hauteur du système de mur']
                 for wall_name, wall in sys_wall['Murs'].items():
-                    if wall.get('Kser ext', None):
+                    if wall.get('Kser ext', None) is not None:
                         kser_ext = wall['Kser ext']
                         couturage_ext_rive = wall['Couturage ext. rive']
                         couturage_ext_inter = wall['Couturage ext. intermédiaire']
-                    if wall.get('Kser int', None):
+                    if wall.get('Kser int', None) is not None:
                         kser_int = wall['Kser int']
                         couturage_int_rive = wall['Couturage int. rive']
                         couturage_int_inter = wall['Couturage int. intermédiaire']
@@ -867,6 +871,16 @@ class MOB(Batiment):
             long_wall = self.data_walls[wall_name]['Longueur']
             nb_anc_lb = floor(long_wall / e_anc_lb)
             nb_anc_li = floor(long_wall / e_anc_li)
+            if nb_anc_lb == 0:
+                raise ValueError(
+                    f"Mur '{wall_name}' (L={long_wall}) : e_anc_lb={e_anc_lb} est supérieur "
+                    f"à la longueur du mur — aucun ancrage lisse basse possible (nb_anc_lb=0)."
+                )
+            if nb_anc_li == 0:
+                raise ValueError(
+                    f"Mur '{wall_name}' (L={long_wall}) : e_anc_li={e_anc_li} est supérieur "
+                    f"à la longueur du mur — aucun ancrage lisse d'implantation possible (nb_anc_li=0)."
+                )
             Ft_Rd_wood = (Ft_Rk_wood * K_mod_montant / Assemblage.GAMMA_M_ASS)
             Fv_Rd_anc_lb = (Fv_Rk_anc_lb * K_mod_lb / Assemblage.GAMMA_M_ASS) * nb_anc_lb
             Fv_Rd_anc_li = (Fv_Rk_anc_li * K_mod_li / Assemblage.GAMMA_M_ASS) * nb_anc_li

{ourocode-2.1.4 → ourocode-2.1.6/ourocode.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: ourocode
-Version: 2.1.4
+Version: 2.1.6
 Summary: Ceci est un catalogue de fonction permettant une utilisation rapide pour la réalisation de note de calcul personnalisée.
 Author-email: Anthony PARISOT <contact@ourea-structure.fr>
 License-Expression: Apache-2.0

{ourocode-2.1.4 → ourocode-2.1.6}/pyproject.toml

@@ -8,7 +8,7 @@ include = ["ourocode*"]
 
 [project]
 name = "ourocode"
-version = "2.1.4"
+version = "2.1.6"
 description = "Ceci est un catalogue de fonction permettant une utilisation rapide pour la réalisation de note de calcul personnalisée."
 readme = "README.md"
 authors = [{name = "Anthony PARISOT", email = "contact@ourea-structure.fr"}]

{ourocode-2.1.4 → ourocode-2.1.6}/tests/test_EC1_Neige.py

@@ -166,3 +166,31 @@ def test_fs():
     # F_s = S_n * entraxe * sin(alpha)
     expected = 1.0 * si.kN/si.m**2 * 1.5 * si.m * math.sin(math.radians(30))
     assert pytest.approx(out[1].value, rel=1e-6) == expected.value
+
+
+def test_poussee_hydrostatique():
+    n = Neige(
+        exposition="Normal",
+        type_toit="1 versant",
+        blocage=False,
+        h_bat=5, d_bat=10, b_bat=10,
+        alpha_toit=10,
+        code_INSEE=1001,
+        alt=100,
+    )
+    h = 1 * si.m
+    latex, res = n.poussee_hydrostatique(h=h)
+    gamma = 2 * si.kN/si.m**3
+    Sk_alt = n.Sk_alt[1]
+    mu_i = min(max(gamma * h / Sk_alt, 0.8), 2)
+    S_n = Sk_alt * mu_i * n.C_e * n.C_t
+    h_neige_expected = S_n / gamma
+    q_expected = gamma * h_neige_expected
+    F_expected = 0.5 * q_expected * h_neige_expected
+    h_equivalent_expected = h_neige_expected / 3
+
+    assert set(res.keys()) == {"q_h", "F_hl", "h_neige", "h_equivalent"}
+    assert pytest.approx(res["h_neige"].value, rel=1e-6) == h_neige_expected.value
+    assert pytest.approx(res["h_equivalent"].value, rel=1e-6) == h_equivalent_expected.value
+    assert pytest.approx(res["q_h"].value, rel=1e-6) == q_expected.value
+    assert pytest.approx(res["F_hl"].value, rel=1e-6) == F_expected.value