@open3cl/engine 1.0.0

package/3.2.2_plancher_bas.spec.js

@@ -0,0 +1,88 @@
+import calc_pb from './3.2.2_plancher_bas.js';
+
+describe('Recherche de bugs dans le calcul de déperdition des planchers bas', () => {
+  describe('calcul de déperdition pour les planchers bas de 2287E1043883T', () => {
+    test('Mur  1 Nord', () => {
+      const zc = 3;
+      const pc_id = 1;
+      const ej = 0;
+      const pb1 = {
+        donnee_entree: {
+          description:
+            'Plancher  1 - Plancher lourd type entrevous terre-cuite, poutrelles béton donnant sur un sous-sol non chauffé avec isolation intrinsèque ou en sous-face (5 cm)',
+          reference: '2022_05_11_17_10_23_1067244006611344',
+          tv_coef_reduction_deperdition_id: 6,
+          enum_type_adjacence_id: '6',
+          surface_paroi_opaque: 42.08,
+          tv_upb0_id: 11,
+          enum_type_plancher_bas_id: '11',
+          enum_methode_saisie_u0_id: '2',
+          enum_type_isolation_id: '4',
+          epaisseur_isolation: 5,
+          enum_methode_saisie_u_id: '3',
+          calcul_ue: 1,
+          perimetre_ue: 27,
+          surface_ue: 42.08,
+          ue: 0.37117494
+        },
+        donnee_intermediaire: {
+          b: 1,
+          upb: 0.5915492957746479,
+          upb_final: 0.37117494,
+          upb0: 2
+        }
+      };
+
+      const pb = {
+        donnee_entree: {
+          description: 'Plancher  2 - Dalle béton donnant sur un terre-plein',
+          reference: '2022_05_11_17_10_51_2415353009491032',
+          tv_coef_reduction_deperdition_id: 5,
+          enum_type_adjacence_id: '5', // Terre-Plein
+          surface_paroi_opaque: 24.34,
+          tv_upb0_id: 9,
+          enum_type_plancher_bas_id: '9', // Dalle béton
+          enum_methode_saisie_u0_id: '2', // déterminé selon le matériau et épaisseur à partir de la table de valeur forfaitaire
+          enum_type_isolation_id: '1', // inconnu
+          enum_periode_isolation_id: '1', // avant 1948
+          tv_upb_id: 2,
+          enum_methode_saisie_u_id: '8', // année de construction saisie (table forfaitaire)
+          calcul_ue: 1,
+          perimetre_ue: 22,
+          surface_ue: 24.34,
+          ue: 0.61789474
+        },
+        donnee_intermediaire: {
+          b: 1,
+          upb: 2,
+          upb_final: 0.61789474,
+          upb0: 2
+        }
+      };
+      const pb_list = [pb1, pb];
+
+      /**
+       * Sortie lib (différence)
+       * tv_upb_id: 8       (2 dans dpe origine)
+       * ue: 0.44           (0.61789474 dans dpe origine)
+       * upb: 0.9           (2 dans dpe origine)
+       * upb_final: 0.44    (0.61789474 dans dpe origine)
+       */
+      calc_pb(pb, zc, pc_id, ej, pb_list);
+
+      /**
+       * P = 22
+       * S = 24.34
+       * 2S/P = (2 * 24.34) / 22 = 2,212727273 = 2 (arrondi à l'entier le plus proche)
+       * upb0 = 2
+       * upb = Min(upb0; upbtab) = Min(2; 2) = 2
+       * ue = Umoyen pour tous les planchers du batiment
+       * upb_final = ue
+       */
+      expect(pb.donnee_intermediaire.b).toBe(1);
+      expect(pb.donnee_intermediaire.upb0).toBe(2);
+      expect(pb.donnee_intermediaire.upb).toBe(2);
+      expect(pb.donnee_intermediaire.upb_final).toBe(pb.donnee_entree.ue);
+    });
+  });
+});

package/3.2.3_plancher_haut.js

@@ -0,0 +1,158 @@
+import enums from './enums.js';
+import b from './3.1_b.js';
+import { tv, requestInput, getKeyByValue, bug_for_bug_compat } from './utils.js';
+
+function tv_uph0(di, de, du) {
+  requestInput(de, du, 'type_plancher_haut');
+  const matcher = {
+    enum_type_plancher_haut_id: de.enum_type_plancher_haut_id
+  };
+  const row = tv('uph0', matcher, de);
+  if (row) {
+    di.uph0 = Number(row.uph0);
+    de.tv_uph0_id = Number(row.tv_uph0_id);
+  } else {
+    console.error('!! pas de valeur forfaitaire trouvée pour uph0 !!');
+  }
+}
+
+function tv_uph(di, de, du, pc_id, zc, effetJoule) {
+  const type_adjacence = requestInput(de, du, 'type_adjacence');
+  const type_ph = requestInput(de, du, 'type_plancher_haut');
+  let type_toiture;
+
+  // From "3CL-DPE 2021" page 21 => "Lorsque le local au-dessus du logement est un local non chauffé, ou un local autre que d’habitation., Uph_tab est pris dans la catégorie « Terrasse »."
+  if (type_adjacence === 'locaux non chauffés non accessible') {
+    type_toiture = 'terrasse';
+  } else if (type_adjacence !== 'extérieur') {
+    type_toiture = 'combles';
+  } else {
+    if (type_ph === 'combles aménagés sous rampant') type_toiture = 'combles';
+    else type_toiture = 'terrasse';
+  }
+
+  if (bug_for_bug_compat && de.tv_uph_id) {
+    /**
+     * Pour certains DPE le type d'adjacence du plancher haut n'est pas bien défini.
+     * On va le récupérer depuis les données d'entrée si besoin
+     */
+    const rowUph = tv('uph', {
+      tv_uph_id: de.tv_uph_id
+    });
+
+    if (rowUph) {
+      if (rowUph.type_toiture === 'combles' && rowUph.type_toiture !== type_toiture) {
+        console.error(
+          `Le type d'adjacence pour le plancher haut ${de.description} ne correspond pas aux données saisies (${rowUph.type_toiture})`
+        );
+        type_toiture = rowUph.type_toiture;
+      }
+
+      if (rowUph.effet_joule !== effetJoule) {
+        console.error(
+          `La variable effet_joule utilisée dans le DPE pour le plancher haut '${de.description}' est ${rowUph.effet_joule}.
+        Celle-ci devrait être ${effetJoule}. La valeur ${rowUph.effet_joule} est conservée dans la suite des calculs`
+        );
+
+        effetJoule = rowUph.effet_joule;
+      }
+    }
+  }
+
+  const matcher = {
+    enum_periode_construction_id: pc_id,
+    enum_zone_climatique_id: zc,
+    effet_joule: effetJoule,
+    type_toiture
+  };
+  const row = tv('uph', matcher);
+  if (row) {
+    di.uph = Number(row.uph);
+    de.tv_uph_id = Number(row.tv_uph_id);
+  } else {
+    console.error('!! pas de valeur forfaitaire trouvée pour uph !!');
+  }
+}
+
+function calc_uph0(di, de, du) {
+  const methode_saisie_u0 = requestInput(de, du, 'methode_saisie_u0');
+  switch (methode_saisie_u0) {
+    case 'type de paroi inconnu (valeur par défaut)':
+    case 'déterminé selon le matériau et épaisseur à partir de la table de valeur forfaitaire':
+      tv_uph0(di, de, du);
+      break;
+    case 'saisie direct u0 justifiée à partir des documents justificatifs autorisés':
+    case "saisie direct u0 correspondant à la performance de la paroi avec son isolation antérieure iti (umur_iti) lorsqu'il y a une surisolation ite réalisée":
+      di.uph0 = requestInput(de, du, 'uph0_saisi');
+      break;
+    case 'u0 non saisi car le u est saisi connu et justifié.':
+      break;
+    default:
+      console.warn('methode_saisie_u0 inconnue:', methode_saisie_u0);
+  }
+}
+
+export default function calc_ph(ph, zc, pc_id, effetJoule) {
+  const de = ph.donnee_entree;
+  const du = {};
+  const di = {};
+  di.uph0 = ph.donnee_intermediaire.uph0;
+
+  b(di, de, du, zc);
+
+  const methode_saisie_u = requestInput(de, du, 'methode_saisie_u');
+
+  switch (methode_saisie_u) {
+    case 'non isolé':
+      calc_uph0(di, de, du);
+      di.uph = di.uph0;
+      break;
+    case 'epaisseur isolation saisie justifiée par mesure ou observation':
+    case 'epaisseur isolation saisie justifiée à partir des documents justificatifs autorisés': {
+      const e = requestInput(de, du, 'epaisseur_isolation', 'int') * 0.01;
+      calc_uph0(di, de, du);
+      di.uph = 1 / (1 / di.uph0 + e / 0.04);
+      break;
+    }
+    case "resistance isolation saisie justifiée observation de l'isolant installé et mesure de son épaisseur":
+    case 'resistance isolation saisie justifiée  à partir des documents justificatifs autorisés': {
+      const r = requestInput(de, du, 'resistance_isolation', 'float');
+      calc_uph0(di, de, du);
+      di.uph = 1 / (1 / di.uph0 + r);
+      break;
+    }
+    case 'isolation inconnue  (table forfaitaire)':
+    case "année d'isolation différente de l'année de construction saisie justifiée (table forfaitaire)":
+      tv_uph(di, de, du, de.enum_periode_isolation_id || pc_id, zc, effetJoule);
+      calc_uph0(di, de, du);
+      di.uph = Math.min(di.uph, di.uph0);
+      break;
+    case 'année de construction saisie (table forfaitaire)': {
+      let pi_id = pc_id;
+      if (de.enum_periode_isolation_id) {
+        pi_id = de.enum_periode_isolation_id;
+      } else {
+        const pc = enums.periode_construction[pc_id];
+        switch (pc) {
+          case 'avant 1948':
+          case '1948-1974':
+            pi_id = getKeyByValue(enums.periode_isolation, '1975-1977');
+            break;
+        }
+      }
+      calc_uph0(di, de, du);
+      tv_uph(di, de, du, pi_id, zc, effetJoule);
+      di.uph = Math.min(di.uph, di.uph0);
+      break;
+    }
+    case 'saisie direct u justifiée  (à partir des documents justificatifs autorisés)':
+    case 'saisie direct u depuis rset/rsee( etude rt2012/re2020)':
+      di.uph = requestInput(de, du, 'uph_saisi', 'float');
+      break;
+    default:
+      console.warn('methode_saisie_u inconnue:', methode_saisie_u);
+  }
+
+  ph.donnee_utilisateur = du;
+  ph.donnee_intermediaire = di;
+}

package/3.3.1.4_porte.js

@@ -0,0 +1,32 @@
+import b from './3.1_b.js';
+import { tv, requestInput } from './utils.js';
+
+function tv_uporte(di, de, du) {
+  requestInput(de, du, 'type_porte');
+  const matcher = {
+    enum_type_porte_id: de.enum_type_porte_id
+  };
+  const row = tv('uporte', matcher, de);
+  if (row) {
+    di.uporte = Number(row.uporte);
+    de.tv_uporte_id = Number(row.tv_uporte_id);
+  } else {
+    console.error('!! pas de valeur forfaitaire trouvée pour uporte !!');
+  }
+}
+
+export default function calc_porte(porte, zc) {
+  const de = porte.donnee_entree;
+  const di = {};
+  const du = {};
+
+  requestInput(de, du, 'surface_porte', 'float');
+
+  b(di, de, du, zc);
+  const methode_saisie_uporte = requestInput(de, du, 'methode_saisie_uporte');
+  if (methode_saisie_uporte === 'valeur forfaitaire') tv_uporte(di, de, du);
+  else di.uporte = requestInput(de, du, 'uporte_saisi', 'float');
+
+  porte.donnee_utilisateur = du;
+  porte.donnee_intermediaire = di;
+}

package/3.3_baie_vitree.js

@@ -0,0 +1,308 @@
+import b from './3.1_b.js';
+import { tv, requestInput, requestInputID, bug_for_bug_compat, getRange } from './utils.js';
+import tvs from './tv.js';
+
+function tv_ug(di, de, du) {
+  const matcher = {
+    enum_type_vitrage_id: requestInputID(de, du, 'type_vitrage')
+  };
+
+  if (matcher.enum_type_vitrage_id && matcher.enum_type_vitrage_id !== '1') {
+    // inside if because simple vitrage does not have these fields
+    matcher.enum_type_gaz_lame_id = requestInputID(de, du, 'type_gaz_lame');
+    matcher.enum_inclinaison_vitrage_id = requestInputID(de, du, 'inclinaison_vitrage');
+    matcher.vitrage_vir = requestInput(de, du, 'vitrage_vir', 'bool');
+    matcher.epaisseur_lame = requestInput(de, du, 'epaisseur_lame', 'float');
+  }
+  const row = tv('ug', matcher);
+  if (row) {
+    di.ug = Number(row.ug);
+    de.tv_ug_id = Number(row.tv_ug_id);
+  } else {
+    console.error('!! pas de valeur forfaitaire trouvée pour ug !!');
+  }
+}
+
+function tv_uw(di, de) {
+  const enum_type_baie_id = de.enum_type_baie_id;
+  const matcher = { enum_type_baie_id };
+  let uw;
+
+  /**
+   * Pas de notion de ug pour les parois polycarbonate ou verre
+   * enum_type_baie_id
+   * 1 - paroi en brique de verre pleine
+   * 2 - paroi en brique de verre creuse
+   * 3 - paroi en polycarbonnate
+   */
+  if (matcher.enum_type_baie_id && ['1', '2', '3'].includes(enum_type_baie_id)) {
+    const row = tv('uw', matcher);
+    if (row) {
+      uw = Number(row.uw);
+      de.tv_uw_id = Number(row.tv_uw_id);
+    }
+  } else {
+    const enum_type_materiaux_menuiserie_id = de.enum_type_materiaux_menuiserie_id;
+    matcher.enum_type_materiaux_menuiserie_id = enum_type_materiaux_menuiserie_id;
+
+    // Récupération de toutes les valeurs de Ug présentes dans les tables Uw pour le type de baie et matériaux de la baie vitrée
+    const ugValues = tvs.uw
+      .filter((row) => {
+        return (
+          row.enum_type_baie_id.split('|').includes(enum_type_baie_id) &&
+          row.enum_type_materiaux_menuiserie_id
+            .split('|')
+            .includes(enum_type_materiaux_menuiserie_id)
+        );
+      })
+      .map((row) => parseFloat(row.ug));
+
+    /**
+     * 3.3.2 Coefficients Uw des fenêtres / portes-fenêtres
+     * Les Uw associés à des Ug non présents dans les tableaux peuvent être obtenus par interpolation ou
+     * extrapolation avec les deux Ug tabulés les plus proches.
+     */
+    let ug1, ug2;
+    [ug1, ug2] = getRange(di.ug, ugValues);
+
+    const matcher_1 = { ...matcher, ...{ ug: `^${ug1}$` } };
+    const row_1 = tv('uw', matcher_1);
+    const delta_ug = ug2 - ug1;
+
+    if (delta_ug === 0) {
+      if (row_1) {
+        uw = Number(row_1.uw);
+      }
+    } else {
+      const matcher_2 = { ...matcher, ...{ ug: `^${ug2}$` } };
+      const row_2 = tv('uw', matcher_2);
+
+      if (row_1 && row_2) {
+        const delta_uw = Number(row_2.uw) - Number(row_1.uw);
+        uw = Number(row_1.uw) + (delta_uw * (di.ug - ug1)) / delta_ug;
+      }
+    }
+  }
+
+  if (uw) {
+    di.uw = uw;
+  } else {
+    console.error(`
+      Pas de valeur forfaitaire uw trouvée pour la baie vitrée ${de.description}.
+    `);
+  }
+}
+
+function tv_deltar(di, de, du) {
+  const matcher = {
+    enum_type_fermeture_id: requestInputID(de, du, 'type_fermeture')
+  };
+  const row = tv('deltar', matcher);
+  if (row) {
+    di.deltar = Number(row.deltar);
+    de.tv_deltar_id = Number(row.tv_deltar_id);
+  } else {
+    console.error('!! pas de valeur forfaitaire trouvée pour deltar !!');
+  }
+}
+
+function tv_ujn(di, de, du) {
+  tv_deltar(di, de, du);
+  const matcher = {
+    deltar: di.deltar,
+    uw: Number.isInteger(di.uw) ? `^${Number(di.uw)}$` : `^${Number(di.uw).toPrecision(2)}$`
+  };
+  const row = tv('ujn', matcher);
+  if (row) {
+    di.ujn = Number(row.ujn);
+    de.tv_ujn_id = Number(row.tv_ujn_id);
+  } else {
+    console.error('!! pas de valeur forfaitaire trouvée pour ujn !!');
+  }
+}
+
+function tv_sw(di, de, du) {
+  const matcher = {
+    enum_type_vitrage_id: requestInputID(de, du, 'type_vitrage')
+  };
+
+  matcher.enum_type_baie_id = requestInputID(de, du, 'type_baie');
+  matcher.enum_type_materiaux_menuiserie_id = requestInputID(de, du, 'type_materiaux_menuiserie');
+
+  /**
+   * 6.2.1 Détermination du facteur solaire
+   * Les champs vitrage_vir et enum_type_pose_id ne sont pas présentes pour les matériaux 'brique de verre' et 'polycarbonate'
+   *
+   * enum_type_materiaux_menuiserie_id
+   * 1 - brique de verre
+   * 2 - polycarbonate
+   */
+  if (![1, 2].includes(parseInt(matcher.enum_type_materiaux_menuiserie_id))) {
+    matcher.vitrage_vir = requestInput(de, du, 'vitrage_vir', 'bool');
+    matcher.enum_type_pose_id = requestInputID(de, du, 'type_pose');
+  }
+
+  const row = tv('sw', matcher);
+  if (row) {
+    di.sw = Number(row.sw);
+    de.tv_sw_id = Number(row.tv_sw_id);
+  } else {
+    console.error('!! pas de valeur forfaitaire trouvée pour sw !!');
+  }
+}
+
+function tv_masque_proche(di, de) {
+  if (!de.tv_coef_masque_proche_id) {
+    di.fe1 = 1;
+    return;
+  }
+  const matcher = {
+    tv_coef_masque_proche_id: de.tv_coef_masque_proche_id // TODO remove
+  };
+  const row = tv('coef_masque_proche', matcher);
+  if (row) {
+    di.fe1 = Number(row.fe1);
+    de.tv_coef_masque_proche_id = Number(row.tv_coef_masque_proche_id);
+  } else {
+    console.error('!! pas de valeur forfaitaire trouvée pour coef_masque_proche !!');
+  }
+}
+
+function tv_masque_lointain_homogene(di, de) {
+  if (!de.tv_coef_masque_lointain_homogene_id) return;
+  const matcher = {
+    tv_coef_masque_lointain_homogene_id: de.tv_coef_masque_lointain_homogene_id // TODO remove
+  };
+  const row = tv('coef_masque_lointain_homogene', matcher);
+  if (row) {
+    di.fe2 = Number(row.fe2);
+    de.tv_coef_masque_lointain_homogene_id = Number(row.tv_coef_masque_lointain_homogene_id);
+  } else {
+    console.log('!! pas de valeur forfaitaire trouvée pour coef_masque_lointain_homogene !!');
+  }
+}
+
+function calc_omb(ml) {
+  const matcher = {
+    tv_coef_masque_lointain_non_homogene_id: ml.tv_coef_masque_lointain_non_homogene_id // TODO remove
+  };
+  const row = tv('coef_masque_lointain_non_homoge', matcher);
+  if (row) {
+    const omb = Number(row.omb);
+    /* de.tv_coef_masque_lointain_non_homogene_id = Number(row.tv_coef_masque_lointain_homogene_id) */
+    return omb;
+  } else {
+    console.log('!! pas de valeur forfaitaire trouvée pour coef_masque_lointain_non_homog !!');
+  }
+}
+
+export default function calc_bv(bv, zc) {
+  const de = bv.donnee_entree;
+  const du = {};
+  const di = {};
+
+  b(di, de, du, zc);
+
+  if (de.sw_saisi) di.sw = de.sw_saisi;
+  else tv_sw(di, de, du);
+
+  tv_ug(di, de, du);
+  if (de.uw_saisi) di.uw = de.uw_saisi;
+  else tv_uw(di, de);
+
+  /**
+   * S'il existe une double-fenêtre, calcul des facteurs sw et uw équivalents
+   * 3.3.2 Coefficients Uw des fenêtres / portes-fenêtres - Traitement des doubles fenêtre
+   * 6.2.1 Détermination du facteur solaire
+   */
+  if (de.double_fenetre === 1 && bv.baie_vitree_double_fenetre) {
+    const deDoubleFenetre = bv.baie_vitree_double_fenetre.donnee_entree;
+    const diDoubleFenetre = bv.baie_vitree_double_fenetre.donnee_intermediaire;
+
+    if (deDoubleFenetre.sw_saisi) {
+      diDoubleFenetre.sw = deDoubleFenetre.sw_saisi;
+    } else {
+      tv_sw(diDoubleFenetre, deDoubleFenetre, du);
+    }
+
+    const sw = di.sw * (diDoubleFenetre.sw || 1);
+
+    if (!de.sw_saisi) {
+      di.sw = sw;
+    }
+
+    if (
+      bug_for_bug_compat &&
+      de.sw_saisi &&
+      bv.donnee_intermediaire.sw.toFixed(5) === sw.toFixed(5)
+    ) {
+      console.error(
+        `Le coefficient sw pour la double fenêtre '${de.description}' est saisi mais il est égale au produit des facteurs 
+          sw des deux fenêtres. Il devrait être sw_saisi, le produit des sw est utilisé.`
+      );
+      di.sw = sw;
+    }
+
+    if (deDoubleFenetre.uw_saisi) {
+      diDoubleFenetre.uw = deDoubleFenetre.uw_saisi;
+    } else {
+      tv_uw(diDoubleFenetre, deDoubleFenetre);
+    }
+
+    const uw = 1 / (1 / di.uw + 1 / diDoubleFenetre.uw + 0.07);
+
+    if (!de.uw_saisi) {
+      di.uw = uw;
+    }
+
+    if (
+      bug_for_bug_compat &&
+      de.uw_saisi &&
+      bv.donnee_intermediaire.uw.toFixed(5) === uw.toFixed(5)
+    ) {
+      console.error(
+        `Le coefficient uw pour la double fenêtre '${de.description}' est saisi mais il est égale au uw équivalent
+          des deux fenêtres. Il devrait être uw_saisi, le facteur uw équivalent est utilisé.`
+      );
+      di.uw = uw;
+    }
+  }
+
+  const type_fermeture = requestInput(de, du, 'type_fermeture');
+  if (type_fermeture !== 'abscence de fermeture pour la baie vitrée') {
+    /**
+     * 3.3.3 Coefficients Ujn des fenêtres/portes-fenêtres
+     * Si le Ujn d’une menuiserie est connu et justifié, le saisir directement
+     */
+    if (de.ujn_saisi) {
+      di.ujn = de.ujn_saisi;
+    } else {
+      tv_ujn(di, de, du);
+    }
+
+    di.u_menuiserie = di.ujn;
+  } else {
+    di.u_menuiserie = di.uw;
+  }
+
+  di.fe2 = 1;
+  if (de.masque_lointain_non_homogene_collection) {
+    let mlnh = de.masque_lointain_non_homogene_collection.masque_lointain_non_homogene || [];
+
+    if (!Array.isArray(mlnh)) {
+      mlnh = [mlnh];
+    }
+
+    di.fe2 = Math.max(
+      0,
+      mlnh.reduce((acc, ml) => acc - calc_omb(ml) / 100, 1)
+    );
+  }
+  tv_masque_proche(di, de, du);
+  tv_masque_lointain_homogene(di, de, du);
+
+  delete di.deltar;
+
+  bv.donnee_utilisateur = du;
+  bv.donnee_intermediaire = di;
+}