### ===================================================================================================================
### NLKA tool - NLKA Element
### ===================================================================================================================
# Copyright ©VIIA 2025
### ===================================================================================================================
### 1. Import modules
### ===================================================================================================================
# General imports
from __future__ import annotations
from typing import TYPE_CHECKING, Dict, Optional
from math import pi, ceil
# References for functions and classes in the rhdhv_fem package
from rhdhv_fem.fem_config import Config
# References for functions and classes in the viiaPackage
if TYPE_CHECKING:
from viiapackage.viiaStatus import ViiaProject
# Import module matplotlib
import matplotlib
# Switch to a non-interactive backend to properly import matplotlib.pyplot
matplotlib.use(Config.MPL_NONINTERACTIVE(notify=False))
import matplotlib.pyplot as plt
### ===================================================================================================================
### 2. Class for NLKA element
### ===================================================================================================================
[docs]class NLKA_Element:
""" This is the class used for creating NLKA element and doing the assessment. This class includes functions to
initiate the NLKA assessment, calculate the seismic demand and capacity, unity check, plot the results, etc."""
# Define the class constant used for the assessment
coefficients = {
'one_way_spanning_wall': {
(5, 3): {'c1': -0.0715, 'c2': 2.717, 'c3': 2.788, 'c4': -1.430},
(5, 2): {'c1': -0.1430, 'c2': 5.505, 'c3': 2.788, 'c4': -1.430},
(5, 1): {'c1': -0.2145, 'c2': 8.294, 'c3': 2.788, 'c4': -1.430},
(4, 3): {'c1': -0.1430, 'c2': 5.434, 'c3': 5.577, 'c4': -1.430},
(4, 2): {'c1': -0.2145, 'c2': 8.222, 'c3': 5.577, 'c4': -1.430},
(4, 1): {'c1': -0.2860, 'c2': 11.011, 'c3': 5.577, 'c4': -1.430}},
'cantilever_wall': {}}
def __init__(
self, project: ViiaProject, name: str, height: float, thickness: float, height_of_center_of_gravity: float,
height_of_building: float, overburden_load: float, density: float, pga: float, frequency: float, e_top: int,
e_bottom: int, angle: float, cantilever: bool = False, spectrum_data: Optional[Dict] = None,
npr_version: str = 'NPR9998:2018', importance_factor: Optional[float] = 1.0):
"""
Input:
- project (obj): VIIA project object containing collections of fem objects and project variables.
- name (string): Name contains the name of the element.
- height (float): The height of the element in [m].
- thickness (float): The thickness of the element in [m].
- height_of_center_of_gravity (float): The height of the center of gravity of the element to the top of the
foundation in [m].
- height_of_building (float): The height of the building to the top of the roof in [m].
- overburden_load (float): The overburden load acting on the element in [kN/m].
- density (float): The density of the material of the element in [kg/m3].
- PGA (float): The peak ground acceleration from webtool in [% of g].
- frequency (float): The fundamental frequency of the building in [Hz].
- e_top (int): The number denotes the eccentricity at the top of the element, please refer to Table H.1 in
NPR9998:2018.
- e_bottom (int): The number denotes the eccentricity at the bottom of the element, please refer to Table
H.1 in NPR9998:2018.
- angle (float): The interstorey drift, refer to Figure H.2 in NPR9998:2018.
- cantilever (boolean): This is True if the element is cantilevered.
- npr_version (str): Version of the NPR. Default value 'NPR9998:2018'.
- importance_factor (float): Multiplied by the values of the accelerations that are from the NPR 9998
webtool, to determine the design values of the spectral acceleration at surface level (ag;d).
"""
self.project = project
self.name = name
self.height = height
self.thickness = thickness
self.z_center = height_of_center_of_gravity
self.height_of_building = height_of_building
self.overburden_load = overburden_load
self.density = density
self.pga = pga
self.frequency = frequency
self.e_top = e_top
self.e_bottom = e_bottom
self.angle = angle
self.cantilever = cantilever
self.spectrum_data = spectrum_data
self.npr_version = npr_version
self.importance_factor = importance_factor
# Variables that are calculated
self.sa_r = None
self.one_way_wall_coeff = None
self.uc = None
self.old_uc = None
self.uc_new_npr_h5a = None
self.uc_new_npr_h5b = None
self.h_5a_result = False
self.h_5b_result = False
self.t_b = None
self.t_c = None
self.t_d = None
self.plateau_factor = None
@property
def weight(self):
""" Calculate the weight of the wall by per meter, in [kN]."""
return self.height * self.thickness * self.density * 9.81 / 1000
@property
def fundamental_period_of_element(self):
"""
This is the function to calculate the fundamental period of the element, refer to the formula H.4 and H.5 in
NPR9998:2018.
Input:
- No input is needed.
Output:
- Fundamental period of the element, in [s].
"""
# Calculate the fundamental period based on the boundary conditions of the element
if self.cantilever:
return (0.65 * self.height * (1 + (self.thickness / self.height) ** 2)) ** (1 / 2) # in s
else:
return (0.28 * self.height / (1 + (self.overburden_load / self.weight) * 2)) ** (1 / 2) # in s
@property
def nominal_slenderness(self):
"""
This is the function to calculate slenderness X', refer to the report 'Interpretation, evaluation and
application of appendix H of NPR 9998:2018'.
Input:
- No input is needed.
Output:
- The slenderness X'.
"""
return self.height / self.thickness
[docs] def amplification_factor(self):
"""
This is the function to calculate the amplification factor, refer to the formula H.3 in NPR9998:2018.
..note:: The Amplification factor should be between 1 and 5.5.
Input:
- No input is needed.
Output:
- Amplification factor.
"""
amplification_factor = \
(3 * (1 + self.z_center / self.height_of_building)) / \
(1 + (1 - self.fundamental_period_of_element / (1 / self.frequency)) ** 2) - 0.5
if amplification_factor <= 1:
amplification_factor = 1.0
return amplification_factor
elif amplification_factor <= 5.5:
return amplification_factor
else:
raise Exception("WARNING: Amplification factor should be between 1 and 5.5")
@property
def seismic_demand(self):
"""
This is the function to calculate design seismic coefficient, refer to the formula H.3 in NPR9998:2018.
Input:
- No input is needed.
Output:
- The design seismic coefficient of the element.
"""
# Decide on the behavior factor by the boundary conditions of the element
if self.cantilever:
q_a = 1
else:
q_a = 2
return self.pga * self.importance_factor / q_a * self.amplification_factor()
[docs] def decide_on_return_period(self):
"""
Function to decide which return period the user selected, if the pga matches none of them, the return period
of 2475 will be selected.
"""
allowed_error = 0.001
if abs(self.spectrum_data['return_period_of_475']['agd'] - self.pga) <= allowed_error:
self.t_b = self.spectrum_data['return_period_of_475']['t_b']
self.t_c = self.spectrum_data['return_period_of_475']['t_c']
self.t_d = self.spectrum_data['return_period_of_475']['t_d']
self.plateau_factor = self.spectrum_data['return_period_of_475']['p']
elif abs(self.spectrum_data['return_period_of_975']['agd'] - self.pga) <= allowed_error:
self.t_b = self.spectrum_data['return_period_of_975']['t_b']
self.t_c = self.spectrum_data['return_period_of_975']['t_c']
self.t_d = self.spectrum_data['return_period_of_975']['t_d']
self.plateau_factor = self.spectrum_data['return_period_of_975']['p']
elif abs(self.spectrum_data['return_period_of_2475']['agd'] - self.pga) <= allowed_error:
self.t_b = self.spectrum_data['return_period_of_2475']['t_b']
self.t_c = self.spectrum_data['return_period_of_2475']['t_c']
self.t_d = self.spectrum_data['return_period_of_2475']['t_d']
self.plateau_factor = self.spectrum_data['return_period_of_2475']['p']
else:
self.project.write_log(
'The input for agd is not corresponding to the data entered for the new NPR at the beginning')
self.pga = self.spectrum_data['return_period_of_2475']['agd']
self.t_b = self.spectrum_data['return_period_of_2475']['t_b']
self.t_c = self.spectrum_data['return_period_of_2475']['t_c']
self.t_d = self.spectrum_data['return_period_of_2475']['t_d']
self.plateau_factor = self.spectrum_data['return_period_of_2475']['p']
@property
def s_ed(self):
if self.spectrum_data:
self.decide_on_return_period()
# By using the formulas from NPR 9998_2018 chapter 3.2.2.2.1
if 0 <= self.fundamental_period_of_element <= self.t_b:
return self.pga * (1 + self.fundamental_period_of_element/self.t_b * (self.plateau_factor - 1))
elif self.t_b < self.fundamental_period_of_element <= self.t_c:
return self.pga * self.plateau_factor
elif self.t_c < self.fundamental_period_of_element <= self.t_d:
return self.pga * self.plateau_factor * (self.t_c / self.fundamental_period_of_element)
elif self.t_d < self.fundamental_period_of_element:
return self.pga * self.plateau_factor * (self.t_c*self.t_d / self.fundamental_period_of_element**2)
else:
return None
[docs] def amplification_factor_new_npr_h5a(self):
"""
This is the function to calculate the amplification factor, refer to the formula H.5a in
NPR9998:2018+C1+A1:2020.
..note:: The Amplification factor should not be larger than 5.5.
Input:
- No input is needed.
Output:
- Amplification factor.
"""
amplification_factor = \
(3 * (1 + self.z_center / self.height_of_building)) / \
(1 + (1 - self.fundamental_period_of_element / (1 / self.frequency)) ** 2) - 0.5
if amplification_factor > 5.5:
raise ValueError(
f'ERROR: Amplification factor for calculating the seismic demand of NSCEs should not be larger than 5.5'
f', but {amplification_factor} is calculated. Please check your input.')
else:
return amplification_factor
[docs] def amplification_factor_new_npr_h5b(self):
"""
This is the function to calculate the amplification factor, refer to the formula H.5b in
NPR9998:2018+C1+A1:2020.
Input:
- No input is needed.
Output:
- Amplification factor.
"""
return 2.5 + 3 * (self.z_center/self.height_of_building)
@property
def seismic_demand_new_npr_h5a(self):
"""
This is the function to calculate design seismic coefficient, refer to the formula H.5a in
NPR9998:2018+C1+A1:2020.
Input:
- No input is needed.
Output:
- The design seismic coefficient of the element.
"""
# Decide on the behavior factor by the boundary conditions of the element
q_a = 2
seismic_demand = self.pga * self.importance_factor / q_a * self.amplification_factor_new_npr_h5a()
if self.s_ed:
if seismic_demand > self.s_ed:
return seismic_demand
else:
return self.s_ed
else:
return None
@property
def seismic_demand_new_npr_h5b(self):
"""
This is the function to calculate design seismic coefficient, refer to the formula H.5b in
NPR9998:2018+C1+A1:2020.
Input:
- No input is needed.
Output:
- The design seismic coefficient of the element.
"""
# Decide on the behavior factor by the boundary conditions of the element
q_a = 2
seismic_demand = self.pga * self.importance_factor / q_a * self.amplification_factor_new_npr_h5b()
if self.s_ed:
if seismic_demand > self.s_ed:
return seismic_demand
else:
return self.s_ed
else:
return None
[docs] def seismic_resistence_one_way_spanning(self):
"""
This is the function to calculate seismic resistance of one way spanning element, refer to the report
'Interpretation, evaluation and application of appendix H of NPR 9998:2018'
Input:
- No input is needed.
Output:
- The seismic resistance of the element.
"""
if self.cantilever:
raise Exception('ERROR: This function is only suitable for one way spanning wall')
try:
coefficients = NLKA_Element.coefficients['one_way_spanning_wall'][(self.e_bottom, self.e_top)]
self.one_way_wall_coeff = coefficients
except KeyError:
raise Exception('Boundary conditions are not correct for one way spanning wall')
sa_r = (coefficients['c1'] * (self.overburden_load / self.weight) ** 2
+ coefficients['c2'] * (self.overburden_load / self.weight)
+ coefficients['c3']) / self.nominal_slenderness + coefficients['c4'] * self.angle
self.sa_r = sa_r
return sa_r
[docs] def seismic_resistence_cantilever(self):
"""
This is the function to calculate seismic resistance of cantilever element, refer to the report 'Interpretation,
evaluation and application of appendix H of NPR 9998:2018'
Input:
- No input is needed.
Output:
- The seismic resistance of the element.
"""
if not self.cantilever:
raise Exception('ERROR: This function is only suitable for cantilever wall')
beta = 3.1
# Check the capacity with the old NPR
if not self.spectrum_data:
if self.e_bottom == 8 and self.e_top == 7:
sa_r = 0.6 * (2 * pi / beta) ** 2 * (
(0.5 + self.overburden_load / self.weight) * (0.975 - 0.025 * self.overburden_load / self.weight) /
self.nominal_slenderness)
elif self.e_bottom == 8 and self.e_top == 6:
sa_r = 0.6 * (2 * pi / beta) ** 2 * (
0.5 * (1 + self.overburden_load / self.weight) * (
0.975 - 0.025 * self.overburden_load / self.weight) /
self.nominal_slenderness)
else:
raise Exception('ERROR: Boundary conditions are not correct for cantilever wall')
else:
if self.e_bottom == 9 and self.e_top == 6:
sa_r = 0.6 * (2 * pi / beta) ** 2 * (
0.5 * (0.975 - 0.025 * self.overburden_load / self.weight) /
self.nominal_slenderness)
elif self.e_bottom == 9 and self.e_top == 7:
sa_r = 0.6 * (2 * pi / beta) ** 2 * (
0.5 * (1 + self.overburden_load / self.weight) * (
0.975 - 0.025 * self.overburden_load / self.weight) /
self.nominal_slenderness)
elif self.e_bottom == 9 and self.e_top == 8:
sa_r = 0.6 * (2 * pi / beta) ** 2 * (
(0.5 + self.overburden_load / self.weight) * (
0.975 - 0.025 * self.overburden_load / self.weight) /
self.nominal_slenderness)
elif self.e_bottom == 10 and self.e_top == 6:
sa_r = 0.6 * (2 * pi / beta) ** 2 * (
(-self.overburden_load) * (0.975 - 0.025 * self.overburden_load / self.weight) /
self.nominal_slenderness)
elif self.e_bottom == 10 and self.e_top == 8:
sa_r = 0.6 * (2 * pi / beta) ** 2 * (
(-self.overburden_load * 0.5) * (0.975 - 0.025 * self.overburden_load / self.weight) /
self.nominal_slenderness)
else:
raise Exception('ERROR: Boundary conditions are not correct for cantilever wall')
self.sa_r = sa_r
return sa_r
[docs] def unity_check_old_npr(self):
"""
This is the function to calculate unity check for the NLKA element.
Input:
- No input is needed.
Output:
- The unity check for the NLKA element.
"""
if self.cantilever:
if self.nominal_slenderness >= 20:
# The unity check for slender cantilever walls is set to 1.1. The seismic resistance is still calculated
# the same way and shown in the results
self.seismic_resistence_cantilever()
uc = 1.1
else:
uc = self.seismic_demand / self.seismic_resistence_cantilever()
else:
uc = self.seismic_demand / self.seismic_resistence_one_way_spanning()
self.old_uc = uc
if uc <= 1:
return True
else:
return False
[docs] def unity_check_new_npr_h5a(self):
"""
This is the function to calculate unity check for the NLKA element by H.5a from NPR9998:2018+C1+A1:2020.
Input:
- No input is needed.
Output:
- The unity check for the NLKA element.
"""
if self.cantilever:
if self.nominal_slenderness >= 20:
# The unity check for slender cantilever walls is set to 1.1. The seismic resistance is still calculated
# the same way and shown in the results
self.seismic_resistence_cantilever()
uc = 1.1
else:
uc = self.seismic_demand_new_npr_h5a / self.seismic_resistence_cantilever()
else:
uc = self.seismic_demand_new_npr_h5a / self.seismic_resistence_one_way_spanning()
self.uc_new_npr_h5a = uc
if uc <= 1:
return True
else:
return False
[docs] def unity_check_new_npr_h5b(self):
"""
This is the function to calculate unity check for the NLKA element H.5b from NPR9998:2018+C1+A1:2020.
Input:
- No input is needed.
Output:
- The unity check for the NLKA element.
"""
if self.cantilever:
if self.nominal_slenderness >= 20:
# The unity check for slender cantilever walls is set to 1.1. The seismic resistance is still calculated
# the same way and shown in the results
self.seismic_resistence_cantilever()
uc = 1.1
else:
uc = self.seismic_demand_new_npr_h5b / self.seismic_resistence_cantilever()
else:
uc = self.seismic_demand_new_npr_h5b / self.seismic_resistence_one_way_spanning()
self.uc_new_npr_h5b = uc
if uc <= 1:
return True
else:
return False
[docs] def unity_check(self, select_5a: bool = False):
""" Unity check based on the applied NPR."""
if self.npr_version == 'NPR9998:2018+C1+A1:2020':
# run the unity checks for new NPR
self.unity_check_new_npr_h5a()
self.unity_check_new_npr_h5b()
# Flag to check user input
if 1 / self.frequency > self.fundamental_period_of_element:
# Switch the H.5b to H.5a formula if user wants to overrule
if select_5a:
uc = self.uc_new_npr_h5a
else:
uc = self.uc_new_npr_h5b
else:
uc = self.uc_new_npr_h5a
else:
# run the unity checks for old NPR
self.unity_check_old_npr()
uc = self.old_uc
self.uc = uc
if uc <= 1:
return True
else:
return False
[docs] def plot_result_one_way_spanning(self):
"""
This is the function to plot the results of the one way spannaning element, together with the charts referring
to the chapter H.5 in NPR9998:2018.
Input:
- No input is needed.
Output:
- The results of the element, together with the chart is saved in the output folder.
"""
# Check if results is calculated
try:
if self.uc:
pass
except AttributeError:
self.unity_check()
# Check if coefficients are present
coefficients = None
try:
coefficients = self.one_way_wall_coeff
except AttributeError:
self.seismic_resistence_one_way_spanning()
coefficients = self.one_way_wall_coeff
# Set up the dict to plot
ar_dict = {i: dict() for i in range(6)}
# Set up the functions to plot
x_axis = [x * 0.5 for x in range(1, 2*ceil(self.nominal_slenderness/5)*5+1)]
def calculate_y(a, b):
y = (coefficients['c1'] * a ** 2 + coefficients['c2'] * a + coefficients['c3']) / b \
+ coefficients['c4'] * self.angle
return y
for p_w, data in ar_dict.items():
data['x'] = x_axis
data['y'] = [calculate_y(p_w, x) for x in x_axis]
self.plot_NPR_curves_with_result(ar_dict)
[docs] def plot_result_cantilever(self):
"""
This is the function to plot the results of the cantilever element, together with the charts referring to the
chapter H.5 in NPR9998:2018.
Input:
- No input is needed.
Output:
- The results of the element, together with the chart is saved in the output folder.
"""
try:
if self.uc:
pass
except AttributeError:
self.unity_check()
beta = 3.1
def calculate_y_condition_1(a, b):
y = 0.6 * (2 * pi / beta) ** 2 * ((0.5 + a) * (0.975 - 0.025 * a) / b)
return y
def calculate_y_condition_2(a, b):
y = 0.6 * (2 * pi / beta) ** 2 * (0.5 * (1 + a) * (0.975 - 0.025 * a) / b)
return y
def calculate_y_new_npr(a, b):
y = None
if self.e_bottom == 9 and self.e_top == 6:
y = 0.6 * (2 * pi / beta) ** 2 * (0.5 * (0.975 - 0.025 * a) / b)
elif self.e_bottom == 9 and self.e_top == 7:
y = 0.6 * (2 * pi / beta) ** 2 * (0.5 * (1 + a) * (0.975 - 0.025 * a) / b)
elif self.e_bottom == 9 and self.e_top == 8:
y = 0.6 * (2 * pi / beta) ** 2 * ((0.5 + a) * (0.975 - 0.025 * a) / b)
return y
# Set up the functions to plot
x_axis = [x * 0.5 for x in range(1, 2*ceil(self.nominal_slenderness/5)*5+1)]
# Set up the dict to plot
ar_dict = {i * 0.5: dict() for i in range(3)}
if self.e_bottom == 8 and self.e_top == 7:
calculate_y = calculate_y_condition_1
elif self.e_bottom == 8 and self.e_top == 6:
calculate_y = calculate_y_condition_2
elif self.e_bottom == 9:
calculate_y = calculate_y_new_npr
for p_w, data in ar_dict.items():
data['x'] = x_axis
data['y'] = [calculate_y_new_npr(p_w, x) for x in x_axis]
self.plot_NPR_curves_with_result(ar_dict)
[docs] def plot_NPR_curves_with_result(self, ar_dict: Dict):
"""
This is the function to plot the results of the element, together with the charts referring to the
chapter H.5 in NPR9998:2018.
Input:
- ar_dict: dictionary contains the data needed for the plot of the NPR curves
Output:
- The results of the element, together with the chart is saved in the output folder.
"""
# Create the graph (close existing previous ones)
plt.close()
plt.figure(figsize=(15, 10))
# Plot NPR curves
ax1 = plt.subplot2grid((1, 1), (0, 0))
ax1.set_ylim([0, 1.0])
if self.nominal_slenderness <= 40:
ax1.set_xlim([0, 40])
else:
ax1.set_xlim([0, ceil(self.nominal_slenderness/5)*5])
plt.xticks([i for i in range(0, ceil(self.nominal_slenderness/5)*5+1, 5)])
plt.yticks([i * 0.1 for i in range(11)])
plt.grid(visible=True, which='major', color='#666666', linestyle='-')
ax1.set_title(f"Seismic resistance plot against slenderness X' for eb = {self.e_bottom}, etop = {self.e_top}")
ax1.set_xlabel("Slenderness X' of masonry wall")
ax1.set_ylabel('Seismic resistance [g]')
for p_w, data in ar_dict.items():
ax1.plot(data['x'], data['y'], label=f'P/W = {p_w}')
plt.legend(loc="upper right")
# Plot the results of this element
plt.plot(self.nominal_slenderness, self.sa_r, marker='o', markersize=5, color="green")
ax1.annotate(f'{self.name} Resistance', (self.nominal_slenderness, self.sa_r))
# Plot the result according to the formulas
if self.h_5a_result:
plt.plot(self.nominal_slenderness, self.seismic_demand_new_npr_h5a, marker='o', markersize=5, color="red")
ax1.annotate(f'{self.name} Demand', (self.nominal_slenderness, self.seismic_demand_new_npr_h5a))
elif self.h_5b_result:
plt.plot(self.nominal_slenderness, self.seismic_demand_new_npr_h5b, marker='o', markersize=5, color="red")
ax1.annotate(f'{self.name} Demand', (self.nominal_slenderness, self.seismic_demand_new_npr_h5b))
else:
plt.plot(self.nominal_slenderness, self.seismic_demand, marker='o', markersize=5, color="red")
ax1.annotate(f'{self.name} Demand', (self.nominal_slenderness, self.seismic_demand))
# Save the picture in the workfolder
image_path = self.project.workfolder_location / 'NSCE Assessment' / f'{self.name}' /\
f'NLKA assessment - {self.name}'
plt.savefig(image_path.as_posix())
[docs] def plot_NPR_curves_envelope(self, input_dict: Dict):
"""
This is the function to plot the results of the element, together with the charts referring to the
chapter H.5 in NPR9998:2018.
Input:
- input_dict (Dict): Dictionary contains the data needed for the plot of the NPR curves.
Output:
- The results of the element, together with the chart is saved in the output folder.
"""
# Create the graph for all the possible overburden loads(close existing previous ones)
plt.close()
plt.figure(figsize=(15, 10))
# Plot NPR curves
ax1 = plt.subplot2grid((1, 1), (0, 0))
ax1.set_ylim([0, 1.0])
ax1.set_xlim([0, 10])
plt.xticks([i for i in range(0, 10, 1)])
plt.yticks([i * 0.1 for i in range(11)])
plt.grid(b=True, which='major', color='#666666', linestyle='-')
ax1.set_title(f"Seismic resistance plot against slenderness X' for eb = {self.e_bottom}, etop = {self.e_top}")
ax1.set_xlabel("Height of masonry wall")
ax1.set_ylabel('Seismic resistance [g]')
for p_w, data in input_dict.items():
ax1.plot(data['x'], data['y'])
plt.legend(loc="upper right")
# Plot the results of this element
plt.plot(self.nominal_slenderness, self.sa_r, marker='o', markersize=5, color="green")
ax1.annotate(f'{self.name} Resistance', (self.nominal_slenderness, self.sa_r))
plt.plot(self.nominal_slenderness, self.seismic_demand, marker='o', markersize=5, color="red")
ax1.annotate(f'{self.name} Demand', (self.nominal_slenderness, self.seismic_demand))
# Save the picture in the workfolder
image_path = self.project.workfolder_location / 'NSCE Assessment' / f'{self.name}' / \
f'NLKA assessment - {self.name} - envelope'
plt.savefig(image_path.as_posix())
### ===================================================================================================================
### 3. End of script
### ===================================================================================================================