irie 0.0.16__py3-none-any.whl → 0.0.17__py3-none-any.whl

irie/apps/prediction/runners/hazus.py

@@ -15,6 +15,8 @@ The seismic design of a bridge is taken into account in terms of the
 (iii) drift limits, and 
 (iv) the longitudinal reinforcement ratio. 
 
+REFERENCES
+
 [1] Hazus earthquake technical manual
     https://www.fema.gov/sites/default/files/2020-10/fema_hazus_earthquake_technical_manual_4-2.pdf
 
@@ -24,148 +26,108 @@ The seismic design of a bridge is taken into account in terms of the
 [3] Mander and Basoz
     https://www.researchgate.net/profile/Jb-Mander/publication/292691534_Seismic_fragility_curve_theory_for_highway_bridges/links/5a7346d7aca2720bc0dbb653/Seismic-fragility-curve-theory-for-highway-bridges.pdf
 
--------------------------------------------
-
-Table 9-6 in [2] (also 7-1 in [1])
-
-HWB1    All    Non-CA <  1990 N/A > 150 N/A EQ1 0 Conventional Major Bridge - Length >  150 meters 
-        All    CA     <  1975 N/A > 150 N/A EQ1 0 Conventional Major Bridge - Length >  150 meters 
-HWB2    All    Non-CA >= 1990 N/A > 150 N/A EQ1 0 Seismic      Major Bridge - Length >  150 meters 
-        All    CA     >= 1975 N/A > 150 N/A EQ1 0 Seismic      Major Bridge - Length >  150 meters 
-HWB3    All    Non-CA <  1990   1   N/A N/A EQ1 1 Conventional Single Span 
-        All    CA     <  1975   1   N/A N/A EQ1 1 Conventional Single Span 
-HWB4    All    Non-CA >= 1990   1   N/A N/A EQ1 1 Seismic      Single Span 
-        All    CA     >= 1975   1   N/A N/A EQ1 1 Seismic      Single Span 
-HWB5  101 106  Non-CA <  1990 N/A   N/A N/A EQ1 0 Conventional Multi-Col. Bent, Simple  Support - Concrete 
-HWB6  101 106  CA     <  1975 N/A   N/A N/A EQ1 0 Conventional Multi-Col. Bent, Simple  Support - Concrete 
-HWB7  101 106  Non-CA >= 1990 N/A   N/A N/A EQ1 0 Seismic      Multi-Col. Bent, Simple  Support - Concrete 
-      101 106  CA     >= 1975 N/A   N/A N/A EQ1 0 Seismic      Multi-Col. Bent, Simple  Support - Concrete
-HWB8  205 206  CA     <  1975 N/A   N/A N/A EQ2 0 Conventional Single Col., Box Girder -  Continuous Concrete 
-HWB9  205 206  CA     >= 1975 N/A   N/A N/A EQ3 0 Seismic      Single Col., Box Girder -  Continuous Concrete 
-HWB10 201 206  Non-CA <  1990 N/A   N/A N/A EQ2 1 Conventional Continuous Concrete 
-      201 206  CA     <  1975 N/A   N/A N/A EQ2 1 Conventional Continuous Concrete 
-HWB11 201 206  Non-CA >= 1990 N/A   N/A N/A EQ3 1 Seismic      Continuous Concrete 
-      201 206  CA     >= 1975 N/A   N/A N/A EQ3 1 Seismic      Continuous Concrete 
-HWB12 301 306  Non-CA <  1990 N/A   N/A No  EQ4 0 Conventional Multi-Col. Bent, Simple  Support - Steel 
-HWB13 301 306  CA     <  1975 N/A   N/A No  EQ4 0 Conventional Multi-Col. Bent, Simple  Support - Steel 
-HWB14 301 306  Non-CA >= 1990 N/A   N/A N/A EQ1 0 Seismic      Multi-Col. Bent, Simple  Support - Steel 
-      301 306  CA     >= 1975 N/A   N/A N/A EQ1 0 Seismic      Multi-Col. Bent, Simple  Support - Steel 
-HWB15 402 410  Non-CA <  1990 N/A   N/A No  EQ5 1 Conventional Continuous Steel 
-      402 410  CA     <  1975 N/A   N/A No  EQ5 1 Conventional Continuous Steel 
-HWB16 402 410  Non-CA >= 1990 N/A   N/A N/A EQ3 1 Seismic      Continuous Steel 
-      402 410  CA     >= 1975 N/A   N/A N/A EQ3 1 Seismic      Continuous Steel 
-HWB17 501 506  Non-CA <  1990 N/A   N/A N/A EQ1 0 Conventional Multi-Col. Bent, Simple  Support - Prestressed  Concrete 
-HWB18 501 506  CA     <  1975 N/A   N/A N/A EQ1 0 Conventional Multi-Col. Bent, Simple  Support - Prestressed  Concrete 
-HWB19 501 506  Non-CA >= 1990 N/A   N/A N/A EQ1 0 Seismic      Multi-Col. Bent, Simple  Support - Prestressed  Concrete
-      501 506  CA     >= 1975 N/A   N/A N/A EQ1 0 Seismic      Multi-Col. Bent, Simple  Support - Prestressed  Concrete 
-HWB20 605 606  CA     <  1975 N/A   N/A N/A EQ2 0 Conventional Single Col., Box Girder -  Prestressed Continuous  Concrete 
-HWB21 605 606  CA     >= 1975 N/A   N/A N/A EQ3 0 Seismic      Single Col., Box Girder -  Prestressed Continuous  Concrete 
-HWB22 601 607  Non-CA <  1990 N/A   N/A N/A EQ2 1 Conventional Continuous Concrete
-      601 607  CA     <  1975 N/A   N/A N/A EQ2 1 Conventional Continuous Concrete
-HWB23 601 607  Non-CA >= 1990 N/A   N/A N/A EQ3 1 Seismic      Continuous Concrete
-      601 607  CA     >= 1975 N/A   N/A N/A EQ3 1 Seismic      Continuous Concrete 
-HWB24 301 306  Non-CA <  1990 N/A   N/A Yes EQ6 0 Conventional Multi-Col. Bent, Simple  Support - Steel 
-HWB25 301 306  CA     <  1975 N/A   N/A Yes EQ6 0 Conventional Multi-Col. Bent, Simple  Support - Steel 
-HWB26 402 410  Non-CA <  1990 N/A   N/A Yes EQ7 1 Conventional Continuous Steel 
-HWB27 402 410  CA     <  1975 N/A   N/A Yes EQ7 1 Conventional Continuous Steel 
-HWB28 N/A N/A All other bridges that are not classified 
-
-
---------------------------
-
-Table 9-7: Hazus Highway System Classification
-
-HWB1  Major Bridge - Length > 150 meters (Conventional Design) 
-HWB2  Major Bridge - Length > 150 meters (Seismic Design) 
-HWB3  Single Span - (Not HWB1 or HWB2) (Conventional Design) 
-HWB4  Single Span - (Not HWB1 or HWB2) (Seismic Design) 
-HWB5  Concrete, Multi-Column Bent, Simple Support (Conventional Design), Non-California (Non CA) 
-HWB6  Concrete, Multi-Column Bent, Simple Support (Conventional Design), California (CA) 
-HWB7  Concrete, Multi-Column Bent, Simple Support (Seismic Design) 
-HWB8  Continuous Concrete, Single Column, Box Girder (Conventional Design) 
-HWB9  Continuous Concrete, Single Column, Box Girder (Seismic Design) 
-HWB10 Continuous Concrete, (Not HWB8 or HWB9) (Conventional Design) 
-HWB11 Continuous Concrete, (Not HWB8 or HWB9) (Seismic Design) 
-HWB12 Steel, Multi-Column Bent, Simple Support (Conventional Design), Non-California (Non-CA) 
-HWB13 Steel, Multi-Column Bent, Simple Support (Conventional Design), California (CA) 
-HWB14 Steel, Multi-Column Bent, Simple Support (Seismic Design) 
-HWB15 Continuous Steel (Conventional Design) 
-HWB16 Continuous Steel (Seismic Design) 
-HWB17 PS Concrete Multi-Column Bent, Simple Support (Conventional Design), Non-California  
-HWB18 PS Concrete, Multi-Column Bent, Simple Support (Conventional Design), California (CA) 
-HWB19 PS Concrete, Multi-Column Bent, Simple Support (Seismic Design) 
-HWB20 PS Concrete, Single Column, Box Girder (Conventional Design) 
-HWB21 PS Concrete, Single Column, Box Girder (Seismic Design)
-HWB22 Continuous Concrete, (Not HWB20/HWB21) (Conventional Design) 
-HWB23 Continuous Concrete, (Not HWB20/HWB21) (Seismic Design) 
-HWB24 Same definition as HWB12 except the bridge length is less than 20 meters 
-HWB25 Same definition as HWB13 except the bridge length is less than 20 meters 
-HWB26 Same definition as HWB15 except the bridge length is less than 20 meters and Non-CA 
-HWB27 Same definition as HWB15 except the bridge length is less than 20 meters and in CA 
-HWB28 All other bridges that are not classified (including wooden bridges)
-
-
------------------------------
-
-Table 7-6 Fragility Function Median Values for Highway Bridges 
-
-HWB1  0.40 0.50 0.70 0.90 3.9 3.9 3.9 13.8 
-HWB2  0.60 0.90 1.10 1.70 3.9 3.9 3.9 13.8 
-HWB3  0.80 1.00 1.20 1.70 3.9 3.9 3.9 13.8 
-HWB4  0.80 1.00 1.20 1.70 3.9 3.9 3.9 13.8 
-HWB5  0.25 0.35 0.45 0.70 3.9 3.9 3.9 13.8 
-HWB6  0.30 0.50 0.60 0.90 3.9 3.9 3.9 13.8 
-HWB7  0.50 0.80 1.10 1.70 3.9 3.9 3.9 13.8 
-HWB8  0.35 0.45 0.55 0.80 3.9 3.9 3.9 13.8 
-HWB9  0.60 0.90 1.30 1.60 3.9 3.9 3.9 13.8 
-HWB10 0.60 0.90 1.10 1.50 3.9 3.9 3.9 13.8 
-HWB11 0.90 0.90 1.10 1.50 3.9 3.9 3.9 13.8 
-HWB12 0.25 0.35 0.45 0.70 3.9 3.9 3.9 13.8 
-HWB13 0.30 0.50 0.60 0.90 3.9 3.9 3.9 13.8 
-HWB14 0.50 0.80 1.10 1.70 3.9 3.9 3.9 13.8 
-HWB15 0.75 0.75 0.75 1.10 3.9 3.9 3.9 13.8 
-HWB16 0.90 0.90 1.10 1.50 3.9 3.9 3.9 13.8 
-HWB17 0.25 0.35 0.45 0.70 3.9 3.9 3.9 13.8 
-HWB18 0.30 0.50 0.60 0.90 3.9 3.9 3.9 13.8 
-HWB19 0.50 0.80 1.10 1.70 3.9 3.9 3.9 13.8 
-HWB20 0.35 0.45 0.55 0.80 3.9 3.9 3.9 13.8 
-HWB21 0.60 0.90 1.30 1.60 3.9 3.9 3.9 13.8 
-HWB22 0.60 0.90 1.10 1.50 3.9 3.9 3.9 13.8 
-HWB23 0.90 0.90 1.10 1.50 3.9 3.9 3.9 13.8 
-HWB24 0.25 0.35 0.45 0.70 3.9 3.9 3.9 13.8 
-HWB25 0.30 0.50 0.60 0.90 3.9 3.9 3.9 13.8 
-HWB26 0.75 0.75 0.75 1.10 3.9 3.9 3.9 13.8 
-HWB27 0.75 0.75 0.75 1.10 3.9 3.9 3.9 13.8 
-HWB28 0.80 1.00 1.20 1.70 3.9 3.9 3.9 13.8 
+
 """
 import math
+import json
 from scipy.stats import norm
+from irie.apps.events.models import EventRecord
+from irie.apps.prediction.runners import Runner, RunID
+from irie.apps.prediction.models import PredictorModel
 
 Slight, Moderate, Extensive, Complete = range(4)
+LEVELS = ["Slight", "Moderate", "Extensive", "Complete"]
 
 # State codes; We'll add more later, right now we assume
 # everything is in California
 class StateCodes:
     California = 22
 
+class HazusRunner(Runner):
+    platform = "mdof"
+
+    schema = {
+      "title": "System ID",
+      "name": "P2",
+      "type": "object",
+      "required": [
+        "name",
+        "decimation",
+        "method",
+        "channels"
+      ],
+      "properties": {
+        "name": {
+          "type": "string",
+          "title": "Name",
+          "description": "Predictor name",
+          "minLength": 2,
+          # "default": "S1"
+        },
+        "method": {
+          "type": "string",
+          "title": "Method",
+          "enum": ["Fourier Spectrum","Response Spectrum","SRIM","OKID"]
+        },
+        "damping": {
+          "type": "number",
+          "title": "Damping",
+          "default": 0.02,
+          "options": {"dependencies": {"method": ["Response Spectrum"]}},
+          "description": "assumed damping ratio"
+        },
+        "channel": {
+          "type": "integer",
+          "title": "Channel",
+          "description": ""
+        }
+      }
+    }
+
+    def render(self):
+        return 
+
+    @classmethod
+    def create(cls, asset, request):
+        predictor = PredictorModel()
+        data = json.loads(request.body)
+
+        data["metrics"] = [""]
+
+        predictor.name = data.pop("name")
+        predictor.config = data
+        predictor.asset = asset
+        predictor.protocol = "IRIE_PREDICTOR_T3"
+        predictor.active = True
+        return predictor
+
+    def newPrediction(self, event):
+        self.event = event
+        return RunID(1)
+
+    def runPrediction(self, run_id: RunID) -> bool:
+        try:
+            self.metric_data = ...
+            return True
+        except Exception as e:
+            self.metric_data = {"error": str(e)}
+            return False
+
+    def getMetricData(self, run, metric):
+        if not hasattr(self, "metric_data"):
+            raise Exception(f"Error {self.name}({id(self)}), {run}")
+        return self.metric_data
+
+
 def hazus_fragility(
-    nbi_data: dict,
-    pga: float = 0.5,  # Peak Ground Acceleration (g)
-    sa_03s: float = 1.1,  # Spectral Acceleration at 0.3 seconds (g)
-    sa_10s: float = 1.4,  # Spectral Acceleration at 1.0 seconds (g)
-    soil_type: str = "B",  # Soil classification ("A", "B", "C", "D", "E")
-    level: int = None,  # Optional: Specify a damage state (0 = Slight, 1 = Moderate, etc.)
-    generate_plot: bool = False,  # Generate fragility curve plot
-    return_data: bool = False
-) -> dict:
+        nbi_data: dict,
+        soil_type: str = "B",  # Soil classification ("A", "B", "C", "D", "E")
+    ) -> dict:
     """
     Compute fragility probabilities for a given bridge using the Hazus methodology.
 
     Args:
     - nbi_data (dict): NBI data containing bridge-specific properties.
-    - pga (float): Peak Ground Acceleration (g).
-    - sa_03s (float): Spectral Acceleration at 0.3 seconds (g).
-    - sa_10s (float): Spectral Acceleration at 1.0 seconds (g).
     - soil_type (str): Soil classification ("A", "B", "C", "D", "E").
     - level (int): Specify a damage state (0 = Slight, 1 = Moderate, etc.) (optional).
     - generate_plot (bool): Whether to generate and return fragility curve plots.
@@ -177,43 +139,53 @@ def hazus_fragility(
     # Step 0: Extract relevant bridge properties
     properties = _bridge_info(nbi_data)
 
+    # Step 1: Determine Hazus bridge type
+    hazus_type: int = _hazus_type(properties)
+    if hazus_type == -1:
+        raise ValueError("Bridge type not found in Hazus classification")
+
+    # Step 3: Generate fragility curve 
+
+    # Adjust sa_range to start from 0
+    sa_range = [0.0] + [0.1 * i for i in range(1, 21)]  # Include 0 explicitly
+
+    curves = {state: [] for state in LEVELS}
+
+    # Generate fragility values, handling the case for Sa = 0
+    for sa in sa_range:
+        for state in LEVELS:
+            median = _get_old_medians(hazus_type)[state]
+            dispersion = 0.6  # β (dispersion factor)
+            curves[state].append(
+                norm.cdf((math.log(sa / median)) / dispersion) if sa > 0 else 0
+            )
+
+    return dict(sa_range=sa_range, curves=curves)
+
+def hazus_prediction(
+        nbi_data: dict,
+        soil_type: str = "B",  # Soil classification ("A", "B", "C", "D", "E")
+        sa_03s: float = 1.1,   # Spectral Acceleration at 0.3 seconds (g)
+        sa_10s: float = 1.4,   # Spectral Acceleration at 1.0 seconds (g)
+        level: int = None,     # Optional: Specify a damage state (0 = Slight, 1 = Moderate, etc.)
+    ):
+    # Step 0: Extract relevant bridge properties
+    properties = _bridge_info(nbi_data)
+
     # Step 1: Determine Hazus bridge type
     hazus_type: int = _hazus_type(properties)
     if hazus_type == -1:
         raise ValueError("Bridge type not found in Hazus classification")
 
     # Step 2: Call _hazus_curve to compute fragility probabilities
-    fragility_probs = _hazus_curve(hazus_type, properties, pga, sa_03s, sa_10s, soil_type)
-
-    # Step 3: Generate fragility curve plot if requested
-    if generate_plot or return_data:
-        # Adjust sa_range to start from 0
-        sa_range = [0.0] + [0.1 * i for i in range(1, 21)]  # Include 0 explicitly
-
-        curves = {state: [] for state in fragility_probs.keys()}
-
-        # Generate fragility values, handling the case for Sa = 0
-        for sa in sa_range:
-            for state in fragility_probs.keys():
-                median = get_old_medians(hazus_type)[state]
-                dispersion = 0.6  # β (dispersion factor)
-                curves[state].append(
-                    norm.cdf((math.log(sa / median)) / dispersion) if sa > 0 else 0
-                )
-        
-        if return_data:
-            fragility_probs["curves"] = curves 
-            fragility_probs["sa_range"] = sa_range
+    fragility_probs = _hazus_curve(hazus_type, properties, sa_03s, sa_10s, soil_type)
 
     # Step 4: Handle `level` as an integer
-    damage_states = ["Slight", "Moderate", "Extensive", "Complete"]
     if level is not None:
         if level not in range(4):
             raise ValueError(f"Invalid level: {level}. Must be an integer in range(4).")
-        damage_state = damage_states[level]
-        return fragility_probs[damage_state]
+        return fragility_probs[LEVELS.index(level)]
 
-    # Step 5: Return all probabilities and optionally the plot
     return fragility_probs
 
 
@@ -254,10 +226,11 @@ def _bridge_info(nbi: dict) -> dict:
 
 
 
-def _hazus_curve(type: int, properties: dict, pga: float, sa_03s: float, sa_10s: float, soil_type: str) -> dict:
+def _hazus_curve(type: int, properties: dict, sa_03s: float, sa_10s: float, soil_type: str) -> dict:
     """
     Compute fragility probabilities for the four damage states and optionally generate fragility curves.
 
+    See page 7-14 of [1]
     Args:
     - type (int): Bridge classification (integer from 1 to 28).
     - properties (dict): Dictionary containing bridge-specific properties.
@@ -265,7 +238,6 @@ def _hazus_curve(type: int, properties: dict, pga: float, sa_03s: float, sa_10s:
     - sa_03s (float): Spectral Acceleration at 0.3 seconds (g).
     - sa_10s (float): Spectral Acceleration at 1.0 seconds (g).
     - soil_type (str): Soil classification ("A", "B", "C", "D", "E").
-    - generate_curve (bool): If True, generate and display fragility curves.
 
     Returns:
     - dict: Fragility probabilities for Slight, Moderate, Extensive, and Complete damage states.
@@ -283,31 +255,32 @@ def _hazus_curve(type: int, properties: dict, pga: float, sa_03s: float, sa_10s:
     span_count = properties["span_count"]
     skew_angle = properties["skew_angle"]
 
-    # Call modify_ground_motion to get the modified values
-    modified_values = modify_ground_motion(pga, sa_03s, sa_10s, soil_type)
-    modified_pga = modified_values['modified_pga']
-    modified_sa_03s = modified_values['modified_sa_03s']
-    modified_sa_10s = modified_values['modified_sa_10s']
+    # Step 2: Get soil-amplified shaking parameters
+    modified_values = modify_ground_motion(soil_type, None, sa_03s, sa_10s)
+    modified_sa_03s = modified_values['sa_03s']
+    modified_sa_10s = modified_values['sa_10s']
 
+    # Step 3: Compute modification factors
     # Compute K_skew, K_shape, K3D
     K_skew = math.sqrt(math.sin(math.radians(90 - skew_angle)))
     if modified_sa_03s == 0:
         raise ValueError("Modified Sa(0.3 sec) cannot be zero.")
+
     K_shape = (2.5 * modified_sa_10s) / modified_sa_03s
-    A, B = get_a_b(type)
+    A, B = _get_a_b(type)
     if span_count - B == 0:
         raise ValueError("Invalid span count resulting in division by zero.")
     K3D = 1 + A / (span_count - B)
 
+    # Step 4: Modify shaking medians
     # Retrieve old medians and compute new medians
-    old_medians = get_old_medians(type)
-    I_shape = get_i_shape(type)
+    old_medians = _get_old_medians(type)
+    I_shape = _get_i_shape(type)
     factor_slight = 1 if I_shape == 0 else min(1, K_shape)
 
     new_medians = {"Slight": old_medians["Slight"] * factor_slight}
-    damage_states = ["Moderate", "Extensive", "Complete"]
     new_medians.update(
-        {state: old_medians[state] * K_skew * K3D for state in damage_states}
+        {state: old_medians[state] * K_skew * K3D for state in ["Moderate", "Extensive", "Complete"]}
     )
 
     # Compute fragility probabilities
@@ -316,6 +289,7 @@ def _hazus_curve(type: int, properties: dict, pga: float, sa_03s: float, sa_10s:
 
     if modified_sa_10s <= 0:
         raise ValueError("Modified Sa(1.0 sec) must be positive.")
+
     fragility_probs = {
         state: compute_probability(modified_sa_10s, median)
         for state, median in new_medians.items()
@@ -324,7 +298,7 @@ def _hazus_curve(type: int, properties: dict, pga: float, sa_03s: float, sa_10s:
     return fragility_probs
 
 
-def modify_ground_motion(pga: float, sa_03s: float, sa_10s: float, soil_type: str) -> dict:
+def modify_ground_motion(soil_type: str, pga: float, sa_03s: float, sa_10s: float) -> dict:
     """
     Modify PGA, Sa(0.3 sec), and Sa(1.0 sec) based on soil amplification factors (Table 4.7).
 
@@ -386,16 +360,18 @@ def modify_ground_motion(pga: float, sa_03s: float, sa_10s: float, soil_type: st
         raise ValueError("Interpolation failed unexpectedly.")
 
     # Calculate modified values
-    modified_pga = pga * get_factor(pga, "FPGA")
-    modified_sa_03s = sa_03s * get_factor(sa_03s, "FA")
-    modified_sa_10s = sa_10s * get_factor(sa_10s, "FV")
+    outputs = {}
+    if pga is not None:
+        outputs["pga"]    =    pga * get_factor(pga, "FPGA")
+    if sa_03s is not None:
+        outputs["sa_03s"] = sa_03s * get_factor(sa_03s, "FA")
+    if sa_10s is not None:
+        outputs["sa_10s"] = sa_10s * get_factor(sa_10s, "FV")
 
-    return {"modified_pga": modified_pga,
-        "modified_sa_03s": modified_sa_03s,
-        "modified_sa_10s": modified_sa_10s,
-    }
+    return outputs
 
-def get_a_b(bridge_type: int) -> tuple:
+
+def _get_a_b(bridge_type: int) -> tuple:
     """
     Retrieve coefficients A and B for K3D calculation based on bridge type.
 
@@ -431,7 +407,7 @@ def get_a_b(bridge_type: int) -> tuple:
     # Retrieve and return A and B values
     return equation_to_ab[equation]
 
-def get_i_shape(bridge_type: int) -> int:
+def _get_i_shape(bridge_type: int) -> int:
     """
     Retrieve I_shape (indicator for skew effects) for the given bridge type
     from Table 7-1.
@@ -456,7 +432,7 @@ def get_i_shape(bridge_type: int) -> int:
 
     return i_shape_mapping[bridge_type]
 
-def get_old_medians(bridge_type: int) -> dict:
+def _get_old_medians(bridge_type: int) -> dict:
     """
     Retrieve the old medians for Slight, Moderate, Extensive, and Complete damage states
     from Table 7-6 based on the bridge type.
@@ -466,17 +442,50 @@ def get_old_medians(bridge_type: int) -> dict:
 
     Returns:
     - dict: Old median values for each damage state.
+
+
+    Table 7-6 Fragility Function Median Values for Highway Bridges 
+
+    HWB1  0.40 0.50 0.70 0.90 3.9 3.9 3.9 13.8 
+    HWB2  0.60 0.90 1.10 1.70 3.9 3.9 3.9 13.8 
+    HWB3  0.80 1.00 1.20 1.70 3.9 3.9 3.9 13.8 
+    HWB4  0.80 1.00 1.20 1.70 3.9 3.9 3.9 13.8 
+    HWB5  0.25 0.35 0.45 0.70 3.9 3.9 3.9 13.8 
+    HWB6  0.30 0.50 0.60 0.90 3.9 3.9 3.9 13.8 
+    HWB7  0.50 0.80 1.10 1.70 3.9 3.9 3.9 13.8 
+    HWB8  0.35 0.45 0.55 0.80 3.9 3.9 3.9 13.8 
+    HWB9  0.60 0.90 1.30 1.60 3.9 3.9 3.9 13.8 
+    HWB10 0.60 0.90 1.10 1.50 3.9 3.9 3.9 13.8 
+    HWB11 0.90 0.90 1.10 1.50 3.9 3.9 3.9 13.8 
+    HWB12 0.25 0.35 0.45 0.70 3.9 3.9 3.9 13.8 
+    HWB13 0.30 0.50 0.60 0.90 3.9 3.9 3.9 13.8 
+    HWB14 0.50 0.80 1.10 1.70 3.9 3.9 3.9 13.8 
+    HWB15 0.75 0.75 0.75 1.10 3.9 3.9 3.9 13.8 
+    HWB16 0.90 0.90 1.10 1.50 3.9 3.9 3.9 13.8 
+    HWB17 0.25 0.35 0.45 0.70 3.9 3.9 3.9 13.8 
+    HWB18 0.30 0.50 0.60 0.90 3.9 3.9 3.9 13.8 
+    HWB19 0.50 0.80 1.10 1.70 3.9 3.9 3.9 13.8 
+    HWB20 0.35 0.45 0.55 0.80 3.9 3.9 3.9 13.8 
+    HWB21 0.60 0.90 1.30 1.60 3.9 3.9 3.9 13.8 
+    HWB22 0.60 0.90 1.10 1.50 3.9 3.9 3.9 13.8 
+    HWB23 0.90 0.90 1.10 1.50 3.9 3.9 3.9 13.8 
+    HWB24 0.25 0.35 0.45 0.70 3.9 3.9 3.9 13.8 
+    HWB25 0.30 0.50 0.60 0.90 3.9 3.9 3.9 13.8 
+    HWB26 0.75 0.75 0.75 1.10 3.9 3.9 3.9 13.8 
+    HWB27 0.75 0.75 0.75 1.10 3.9 3.9 3.9 13.8 
+    HWB28 0.80 1.00 1.20 1.70 3.9 3.9 3.9 13.8 
     """
+
     old_medians = {
-        1: {"Slight": 0.40, "Moderate": 0.50, "Extensive": 0.70, "Complete": 0.90},
-        2: {"Slight": 0.60, "Moderate": 0.90, "Extensive": 1.10, "Complete": 1.70},
-        3: {"Slight": 0.80, "Moderate": 1.00, "Extensive": 1.20, "Complete": 1.70},
-        4: {"Slight": 0.80, "Moderate": 1.00, "Extensive": 1.20, "Complete": 1.70},
-        5: {"Slight": 0.25, "Moderate": 0.35, "Extensive": 0.45, "Complete": 0.70},
-        6: {"Slight": 0.30, "Moderate": 0.50, "Extensive": 0.60, "Complete": 0.90},
-        7: {"Slight": 0.50, "Moderate": 0.80, "Extensive": 1.10, "Complete": 1.70},
-        8: {"Slight": 0.35, "Moderate": 0.45, "Extensive": 0.55, "Complete": 0.80},
-        9: {"Slight": 0.60, "Moderate": 0.90, "Extensive": 1.30, "Complete": 1.60},
+        1:  {"Slight": 0.40, "Moderate": 0.50, "Extensive": 0.70, "Complete": 0.90},
+        2:  {"Slight": 0.60, "Moderate": 0.90, "Extensive": 1.10, "Complete": 1.70},
+        3:  {"Slight": 0.80, "Moderate": 1.00, "Extensive": 1.20, "Complete": 1.70},
+        4:  {"Slight": 0.80, "Moderate": 1.00, "Extensive": 1.20, "Complete": 1.70},
+        5:  {"Slight": 0.25, "Moderate": 0.35, "Extensive": 0.45, "Complete": 0.70},
+        6:  {"Slight": 0.30, "Moderate": 0.50, "Extensive": 0.60, "Complete": 0.90},
+        7:  {"Slight": 0.50, "Moderate": 0.80, "Extensive": 1.10, "Complete": 1.70},
+        8:  {"Slight": 0.35, "Moderate": 0.45, "Extensive": 0.55, "Complete": 0.80},
+        9:  {"Slight": 0.60, "Moderate": 0.90, "Extensive": 1.30, "Complete": 1.60},
         10: {"Slight": 0.60, "Moderate": 0.90, "Extensive": 1.10, "Complete": 1.50},
         11: {"Slight": 0.90, "Moderate": 0.90, "Extensive": 1.10, "Complete": 1.50},
         12: {"Slight": 0.25, "Moderate": 0.35, "Extensive": 0.45, "Complete": 0.70},
@@ -513,6 +522,86 @@ def _hazus_type(properties: dict) -> int:
         
     Returns:
         int: Hazus type classification (1-28) or -1 if no match is found
+
+
+    Table 9-6 in [2] (also 7-1 in [1])
+
+    HWB1    All    Non-CA <  1990 N/A > 150 N/A EQ1 0 Conventional Major Bridge - Length >  150 meters 
+            All    CA     <  1975 N/A > 150 N/A EQ1 0 Conventional Major Bridge - Length >  150 meters 
+    HWB2    All    Non-CA >= 1990 N/A > 150 N/A EQ1 0 Seismic      Major Bridge - Length >  150 meters 
+            All    CA     >= 1975 N/A > 150 N/A EQ1 0 Seismic      Major Bridge - Length >  150 meters 
+    HWB3    All    Non-CA <  1990   1   N/A N/A EQ1 1 Conventional Single Span 
+            All    CA     <  1975   1   N/A N/A EQ1 1 Conventional Single Span 
+    HWB4    All    Non-CA >= 1990   1   N/A N/A EQ1 1 Seismic      Single Span 
+            All    CA     >= 1975   1   N/A N/A EQ1 1 Seismic      Single Span 
+    HWB5  101 106  Non-CA <  1990 N/A   N/A N/A EQ1 0 Conventional Multi-Col. Bent, Simple  Support - Concrete 
+    HWB6  101 106  CA     <  1975 N/A   N/A N/A EQ1 0 Conventional Multi-Col. Bent, Simple  Support - Concrete 
+    HWB7  101 106  Non-CA >= 1990 N/A   N/A N/A EQ1 0 Seismic      Multi-Col. Bent, Simple  Support - Concrete 
+        101 106  CA     >= 1975 N/A   N/A N/A EQ1 0 Seismic      Multi-Col. Bent, Simple  Support - Concrete
+    HWB8  205 206  CA     <  1975 N/A   N/A N/A EQ2 0 Conventional Single Col., Box Girder -  Continuous Concrete 
+    HWB9  205 206  CA     >= 1975 N/A   N/A N/A EQ3 0 Seismic      Single Col., Box Girder -  Continuous Concrete 
+    HWB10 201 206  Non-CA <  1990 N/A   N/A N/A EQ2 1 Conventional Continuous Concrete 
+        201 206  CA     <  1975 N/A   N/A N/A EQ2 1 Conventional Continuous Concrete 
+    HWB11 201 206  Non-CA >= 1990 N/A   N/A N/A EQ3 1 Seismic      Continuous Concrete 
+        201 206  CA     >= 1975 N/A   N/A N/A EQ3 1 Seismic      Continuous Concrete 
+    HWB12 301 306  Non-CA <  1990 N/A   N/A No  EQ4 0 Conventional Multi-Col. Bent, Simple  Support - Steel 
+    HWB13 301 306  CA     <  1975 N/A   N/A No  EQ4 0 Conventional Multi-Col. Bent, Simple  Support - Steel 
+    HWB14 301 306  Non-CA >= 1990 N/A   N/A N/A EQ1 0 Seismic      Multi-Col. Bent, Simple  Support - Steel 
+        301 306  CA     >= 1975 N/A   N/A N/A EQ1 0 Seismic      Multi-Col. Bent, Simple  Support - Steel 
+    HWB15 402 410  Non-CA <  1990 N/A   N/A No  EQ5 1 Conventional Continuous Steel 
+        402 410  CA     <  1975 N/A   N/A No  EQ5 1 Conventional Continuous Steel 
+    HWB16 402 410  Non-CA >= 1990 N/A   N/A N/A EQ3 1 Seismic      Continuous Steel 
+        402 410  CA     >= 1975 N/A   N/A N/A EQ3 1 Seismic      Continuous Steel 
+    HWB17 501 506  Non-CA <  1990 N/A   N/A N/A EQ1 0 Conventional Multi-Col. Bent, Simple  Support - Prestressed  Concrete 
+    HWB18 501 506  CA     <  1975 N/A   N/A N/A EQ1 0 Conventional Multi-Col. Bent, Simple  Support - Prestressed  Concrete 
+    HWB19 501 506  Non-CA >= 1990 N/A   N/A N/A EQ1 0 Seismic      Multi-Col. Bent, Simple  Support - Prestressed  Concrete
+        501 506  CA     >= 1975 N/A   N/A N/A EQ1 0 Seismic      Multi-Col. Bent, Simple  Support - Prestressed  Concrete 
+    HWB20 605 606  CA     <  1975 N/A   N/A N/A EQ2 0 Conventional Single Col., Box Girder -  Prestressed Continuous  Concrete 
+    HWB21 605 606  CA     >= 1975 N/A   N/A N/A EQ3 0 Seismic      Single Col., Box Girder -  Prestressed Continuous  Concrete 
+    HWB22 601 607  Non-CA <  1990 N/A   N/A N/A EQ2 1 Conventional Continuous Concrete
+        601 607  CA     <  1975 N/A   N/A N/A EQ2 1 Conventional Continuous Concrete
+    HWB23 601 607  Non-CA >= 1990 N/A   N/A N/A EQ3 1 Seismic      Continuous Concrete
+        601 607  CA     >= 1975 N/A   N/A N/A EQ3 1 Seismic      Continuous Concrete 
+    HWB24 301 306  Non-CA <  1990 N/A   N/A Yes EQ6 0 Conventional Multi-Col. Bent, Simple  Support - Steel 
+    HWB25 301 306  CA     <  1975 N/A   N/A Yes EQ6 0 Conventional Multi-Col. Bent, Simple  Support - Steel 
+    HWB26 402 410  Non-CA <  1990 N/A   N/A Yes EQ7 1 Conventional Continuous Steel 
+    HWB27 402 410  CA     <  1975 N/A   N/A Yes EQ7 1 Conventional Continuous Steel 
+    HWB28 N/A N/A All other bridges that are not classified 
+
+
+    --------------------------
+
+    Table 9-7: Hazus Highway System Classification
+
+    HWB1  Major Bridge - Length > 150 meters (Conventional Design) 
+    HWB2  Major Bridge - Length > 150 meters (Seismic Design) 
+    HWB3  Single Span - (Not HWB1 or HWB2) (Conventional Design) 
+    HWB4  Single Span - (Not HWB1 or HWB2) (Seismic Design) 
+    HWB5  Concrete, Multi-Column Bent, Simple Support (Conventional Design), Non-California (Non CA) 
+    HWB6  Concrete, Multi-Column Bent, Simple Support (Conventional Design), California (CA) 
+    HWB7  Concrete, Multi-Column Bent, Simple Support (Seismic Design) 
+    HWB8  Continuous Concrete, Single Column, Box Girder (Conventional Design) 
+    HWB9  Continuous Concrete, Single Column, Box Girder (Seismic Design) 
+    HWB10 Continuous Concrete, (Not HWB8 or HWB9) (Conventional Design) 
+    HWB11 Continuous Concrete, (Not HWB8 or HWB9) (Seismic Design) 
+    HWB12 Steel, Multi-Column Bent, Simple Support (Conventional Design), Non-California (Non-CA) 
+    HWB13 Steel, Multi-Column Bent, Simple Support (Conventional Design), California (CA) 
+    HWB14 Steel, Multi-Column Bent, Simple Support (Seismic Design) 
+    HWB15 Continuous Steel (Conventional Design) 
+    HWB16 Continuous Steel (Seismic Design) 
+    HWB17 PS Concrete Multi-Column Bent, Simple Support (Conventional Design), Non-California  
+    HWB18 PS Concrete, Multi-Column Bent, Simple Support (Conventional Design), California (CA) 
+    HWB19 PS Concrete, Multi-Column Bent, Simple Support (Seismic Design) 
+    HWB20 PS Concrete, Single Column, Box Girder (Conventional Design) 
+    HWB21 PS Concrete, Single Column, Box Girder (Seismic Design)
+    HWB22 Continuous Concrete, (Not HWB20/HWB21) (Conventional Design) 
+    HWB23 Continuous Concrete, (Not HWB20/HWB21) (Seismic Design) 
+    HWB24 Same definition as HWB12 except the bridge length is less than 20 meters 
+    HWB25 Same definition as HWB13 except the bridge length is less than 20 meters 
+    HWB26 Same definition as HWB15 except the bridge length is less than 20 meters and Non-CA 
+    HWB27 Same definition as HWB15 except the bridge length is less than 20 meters and in CA 
+    HWB28 All other bridges that are not classified (including wooden bridges)
+
     """
     year_built = properties["year_built"]
     span_count = properties["span_count"]