Combinaciones_de_acciones/Comb_acciones.py

import os
import sys
import comtypes.client
import pandas as pd
import itertools
import tkinter as tk
from tkinter import *
from tkinter.filedialog import askopenfilename
import math
import numpy as np



# Print iterations progress
def printProgressBar (iteration, total, prefix = '', suffix = '', decimals = 1, length = 100, fill = '█', printEnd = "\r"):
    """
    Call in a loop to create terminal progress bar
    @params:
        iteration   - Required  : current iteration (Int)
        total       - Required  : total iterations (Int)
        prefix      - Optional  : prefix string (Str)
        suffix      - Optional  : suffix string (Str)
        decimals    - Optional  : positive number of decimals in percent complete (Int)
        length      - Optional  : character length of bar (Int)
        fill        - Optional  : bar fill character (Str)
        printEnd    - Optional  : end character (e.g. "\r", "\r\n") (Str)
    """
    percent = ("{0:." + str(decimals) + "f}").format(100 * (iteration / float(total)))
    filledLength = int(length * iteration // total)
    bar = fill * filledLength + '-' * (length - filledLength)
    print(f'\r{prefix} |{bar}| {percent}% {suffix}', end = printEnd)
    # Print New Line on Complete
    if iteration == total: 
        print()


#diccionario donde se encuentran todos los coef de comb acciones de las normas
normas = {
    "IAPF" : {
        "CoefVar" : {
            "Trafico" : {
                "psi0" : 0.8,
                "psi1" : {
                    "1" : 0.8,
                    "2" : 0.6,
                    "3" : 0.4
                },
                "psi2" : 0
            },
            "Resto" : {
                "psi0" : 0.6,
                "psi1" : 0.5,
                "psi2" : 0.2
            }
        },
        "ELU" : {
            "Favorable" : {
                "Persistente" : {
                    "Cte" : 1,
                    "NoCte" : 1,
                    "Variable" : 0
                },
                "Accidental" : {
                    "Cte" : 1,
                    "NoCte" : 1,
                    "Variable" : 0,
                    "Accidental" : 1
                }
            },
            "Desfavorable" : {
                "Persistente" : {
                    "Cte" : 1.35,
                    "NoCte" : 1.5,
                    "Variable" : 1.5
                },
                "Accidental" : {
                    "Cte" : 1,
                    "NoCte" : 1,
                    "Variable" : 1,
                    "Accidental" : 1
                }
            }
        },
        "ELS" : {
            "Favorable" : {
                "Persistente" : {
                    "Cte" : 1,
                    "NoCte" : 1,
                    "Variable" : 0
                }
            },
            "Desfavorable" : {
                "Persistente" : {
                    "Cte" : 1,
                    "NoCte" : 1,
                    "Variable" : 1
                }
            }
        }
    },
    "IAP" : {
        "CoefVar" : {
            "SCUso" : {
                "CVerticales" : {
                    "psi0" : 0.4, #cuidado que para vehiculos pesados es 0.75
                    "psi1" : 0.4, #cuidado que para vehiculos pesados es 0.75
                    "psi2" : 0
                },
                "CHorizontales" : {
                    "psi0" : 0, 
                    "psi1" : 0, 
                    "psi2" : 0
                },
            },
            "SCConstruccion" : {
                "psi0" : 1, 
                "psi1" : 0, 
                "psi2" : 1
            },
        },
        "ELU" : {
            "Favorable" : {
                "Persistente" : {
                    "Cte" : 1,
                    "NoCte" : 1,
                    "Variable" : 0
                }
            },
            "Desfavorable" : {
                "Persistente" : {
                    "Cte" : 1.35,
                    "NoCte" : 1.5,
                    "Variable" : 1.5
                }
            }
        },
        "ELS" : {
            "Favorable" : {
                "Persistente" : {
                    "Cte" : 1,
                    "NoCte" : 1,
                    "Variable" : 0
                }
            },
            "Desfavorable" : {
                "Persistente" : {
                    "Cte" : 1,
                    "NoCte" : 1,
                    "Variable" : 1
                }
            }
        }
    }
}


#Conexion con SAP2000

helper = comtypes.client.CreateObject('SAP2000v1.Helper')
helper = helper.QueryInterface(comtypes.gen.SAP2000v1.cHelper)
mySapObject = helper.GetObject("CSI.SAP2000.API.SapObject")
SapModel = mySapObject.SapModel
SapModel.SetModelisLocked(False)

#se pregunta si se desea borrar todo lo que hay en el SAP

def borrar_combos_existentes():
    comboNames = SapModel.RespCombo.GetNameList()
    if comboNames[0] > 0:
        for name in comboNames[1]:
            ret = SapModel.RespCombo.Delete(name)
            printProgressBar(comboNames[1].index(name)+1, comboNames[0], prefix = 'Borrando combinaciones existentes:', suffix = 'Completado', length = 50)

def borrar_load_cases():
    case_names = SapModel.LoadCases.GetNameList()
    if case_names[0] > 0:
        for name in case_names[1]:
            ret = SapModel.LoadCases.Delete(name)
            printProgressBar(case_names[1].index(name)+1, case_names[0], prefix = 'Borrando casos de carga existentes:', suffix = 'Completado', length = 50)


def borrar_patrones_carga():
    patternNames = SapModel.LoadPatterns.GetNameList()
    if patternNames[0] > 0:
        for name in patternNames[1]:
            ret = SapModel.LoadPatterns.Delete(name)
            printProgressBar(patternNames[1].index(name)+1, patternNames[0], prefix = 'Borrando patrones de carga existentes:', suffix = 'Completado', length = 50)


respuesta = input("¿Deseas borrar todos las combinaciones existentes? (Y/n): ")
if respuesta.lower() == 'y' or len(respuesta) == 0:
    borrar_combos_existentes()

respuesta = input("¿Deseas borrar todos las cargas existentes? (Y/n): ")
if respuesta.lower() == 'y' or len(respuesta) == 0:
    borrar_load_cases()
    borrar_patrones_carga()


#se selecciona la norma a utilizar para la comb acciones
respuesta = input("¿La norma a seguir es la IAPF - 07? (Y / n [IAP - 11]): ")
if respuesta.lower() == 'y' or len(respuesta) == 0:
    norma_proyecto = "IAPF"
    vias_cargadas = input("¿Cuantas vias hay cargadas a la vez 1, 2 o 3 (o más)?")

else: "IAP"



#Extraccion de datos del excel


tk.Tk().withdraw()
ruta_excel = askopenfilename()

#ruta_excel = "C:/Users/Daniel.p/Documents/14. Ay-A 0042 Foso ataque empuje tubos y cimentaciones grua ( PALENCIA)/1. MEMORIA/Comb_acciones.xlsx"

compatibilidades_df = pd.read_excel(ruta_excel,
                              sheet_name='Compatibilidades')

restricciones_df = pd.read_excel(ruta_excel,
                              sheet_name='Restricciones')


#Limpieza de filas vacias

x, y = compatibilidades_df.shape

for i in range(x):
    if (compatibilidades_df.loc[i, compatibilidades_df.columns.values[0]] == 0):
        compatibilidades_df=compatibilidades_df.drop(i)



for i in range(y-1, 0, -1):
    col_name = compatibilidades_df.columns.values[i]
    if (isinstance(compatibilidades_df[col_name].name, int)):
        compatibilidades_df = compatibilidades_df.drop(col_name, axis=1)
    elif isinstance(col_name, str) and col_name.startswith("0"):
        compatibilidades_df = compatibilidades_df.drop(col_name, axis=1) 


x, y = compatibilidades_df.shape

for i in range(1, x+1):
    for j in range(i-1):
        if not isinstance(compatibilidades_df.loc[j, compatibilidades_df.columns.values[i]], str):
            if math.isnan(compatibilidades_df.loc[j, compatibilidades_df.columns.values[i]]):
                compatibilidades_df.loc[j, compatibilidades_df.columns.values[i]] = 'r'



#cambio de los valores Nan por 0
compatibilidades_df = compatibilidades_df.fillna(0)
print(compatibilidades_df)

restricciones_df = restricciones_df.fillna(0)


 
#creacion de load patterns

for i in restricciones_df.loc[:, "Definición de 'Load Cases'"]:
    ret = SapModel.LoadPatterns.Add(i, 3) #el numero indica el tipo de carga: 1 dead, 3 live ..


#Se almacena en esta variable un diccionario indicando los estado en los que se puede encontrar la carga

valores_por_patron = {}
for index, row in restricciones_df.iterrows():
    patron = restricciones_df.loc[:,restricciones_df.columns[0]]

    valores = [0, 1]

    if row["Bidireccional?"]:
        valores.append(-1)
    
    valores_por_patron[patron[index]] = list(set(valores))


#creacion del array donde se almancenan los nombres de las cargas
patrones = list(valores_por_patron.keys())



#se hace con itertools todas las combinaciones de carga posibles sin ninguna restriccion
todas_combinaciones = list(itertools.product(*[valores_por_patron[pat] for pat in patrones]))


#funcion para comprobar que la combiacion de cargas es posible debido a las restricciones

def validar_combinacion(comb, patrones, compatibilidad_df, restricciones_df):
    #se comprueba la matriz de compatibilidad buscando por cada caga que este activa en este combo
    #otra targa que este activa y se comprueba en la matriz buscando una 'r'
    for i in range(len(patrones)):
        if comb[i]:
            for j in range(i+1, len(patrones)):
                if comb[j]:
                    if compatibilidad_df.iloc[i,j+1] == 'r':
                        return False
    #se comprueba que la carga que esta activa esta con otra carga especificada como condicional
    #para que la primera exista
    restriccion = restricciones_df.loc[:,"Si y solo si"]           
    for iter, res in enumerate(restriccion):
        if res:
            idx = patrones.index(res)
            if comb[iter] and not comb[idx]:
                    return False
    #comprueba que las cargas especificadas como permanentes se encuentran en este caso de carga
    permanentes = restricciones_df.loc[:,"Permanentes"]
    for iter, per in enumerate(permanentes):
        if per != 0:
            if not(per=='x' and comb[iter]):
                return False
    #comprueba que no hay dos acciones accidentales en la misma combinacion
    accidentales = restricciones_df.loc[:,"Tipo"]
    num_acc = 0
    for iter, acc in enumerate(accidentales):
        if acc == "Accidental":
            if comb[iter]:
                num_acc += 1
                if num_acc > 1:
                    return False
    #comprueba que con sismo solo hay una carga variable
    sismo = restricciones_df.loc[:,"Accion IAPF"]
    num_var = 0
    for iter, sis in enumerate(sismo):
        if sis == "Sismo":
            if comb[iter]:
                for jter, _ in enumerate(comb):
                    if (restricciones_df.loc[jter,"Tipo"] == "Variable") and comb[jter] != 0:
                        num_var += 1
                        if num_var > 1:
                            return False
    #si no cumple ninguna de las restricciones se considera válida la combinacion
    return True


#comprobacion de la validez de las combinaciones
combinaciones_validas = []
for comb in todas_combinaciones:
    if validar_combinacion(comb, patrones, compatibilidades_df, restricciones_df):
        combinaciones_validas.append(comb)


#funcion para hacer la suma en valor absoluto de un array
def suma_abs(lista):

    for _ in range(len(lista)): total += abs(lista[_])

    return total

num_permanentes = 0

#se obtiene el numero de cargas que son permanentes
for iter, per in enumerate(restricciones_df.loc[:,"Permanentes"]):
    if per == 'x' : num_permanentes += 1



def CrearCombo(nombre, numero, combo, coef):

    global SapModel
    global patrones

    nombre_combo = nombre + "{:04}".format(numero)

    ret = SapModel.RespCombo.Add(nombre_combo, 0)

    for idx, patron in enumerate(patrones):
        coeficiente = combo[idx]
        if coeficiente != 0:
            ret = SapModel.RespCombo.SetCaseList(nombre_combo, 0, patron, coeficiente*coef[idx])



combos_ELUP = []
combos_ELUA = []
combos_ELUS = []
combos_ELSC = []
combos_ELSF = []
combos_ELSCP = []
minorados_ELUP = []
minorados_ELSC = []
minorados_ELSF = []
minorados_ELSCP = []
minorados_ELUA = []
minorados_ELUS = []


for combo in combinaciones_validas:

    coeficientes_ELUP = np.zeros(len(patrones))
    coeficientes_ELUA = np.zeros(len(patrones))
    coeficientes_ELUS = np.zeros(len(patrones))
    coeficientes_ELSC = np.zeros(len(patrones))
    coeficientes_ELSF = np.zeros(len(patrones))
    coeficientes_ELSCP = np.zeros(len(patrones))
    n_var = 0
    n_acc = 0

    for i in range(len(combo)):


        if norma_proyecto == "IAPF":
            #primero se comprueba si la carga es permanente o no
            if restricciones_df.loc[i,"Permanentes"]:
                #luego constante o no
                if restricciones_df.loc[i,"Tipo"] == "Constante":
                    if restricciones_df.loc[i,"Direccion"] == "Desfavorable":
                        #ELU
                        coeficientes_ELUP[i] = normas["IAPF"]["ELU"]["Desfavorable"]["Persistente"]["Cte"]
                        coeficientes_ELUA[i] = normas["IAPF"]["ELU"]["Desfavorable"]["Accidental"]["Cte"]
                        coeficientes_ELUS[i] = normas["IAPF"]["ELU"]["Desfavorable"]["Accidental"]["Cte"]
                        #ELS
                        coeficientes_ELSC[i] = normas["IAPF"]["ELS"]["Desfavorable"]["Persistente"]["Cte"]
                        coeficientes_ELSF[i] = normas["IAPF"]["ELS"]["Desfavorable"]["Persistente"]["Cte"]
                        coeficientes_ELSCP[i] = normas["IAPF"]["ELS"]["Desfavorable"]["Persistente"]["Cte"]
                    elif restricciones_df.loc[i,"Direccion"] == "Favorable":
                        #ELU
                        coeficientes_ELUP[i] = normas["IAPF"]["ELU"]["Favorable"]["Persistente"]["Cte"]
                        coeficientes_ELUA[i] = normas["IAPF"]["ELU"]["Favorable"]["Accidental"]["Cte"]
                        coeficientes_ELUS[i] = normas["IAPF"]["ELU"]["Favorable"]["Accidental"]["Cte"]
                        #ELS
                        coeficientes_ELSC[i] = normas["IAPF"]["ELS"]["Favorable"]["Persistente"]["Cte"]
                        coeficientes_ELSF[i] = normas["IAPF"]["ELS"]["Favorable"]["Persistente"]["Cte"]
                        coeficientes_ELSCP[i] = normas["IAPF"]["ELS"]["Favorable"]["Persistente"]["Cte"]
                elif restricciones_df.loc[i,"Tipo"] == "No Constante":
                    if restricciones_df.loc[i,"Direccion"] == "Desfavorable":
                        #ELU
                        coeficientes_ELUP[i] = normas["IAPF"]["ELU"]["Desfavorable"]["Persistente"]["NoCte"]
                        coeficientes_ELUA[i] = normas["IAPF"]["ELU"]["Desfavorable"]["Accidental"]["NoCte"]
                        coeficientes_ELUS[i] = normas["IAPF"]["ELU"]["Desfavorable"]["Accidental"]["NoCte"]
                        #ELS
                        coeficientes_ELSC[i] = normas["IAPF"]["ELS"]["Desfavorable"]["Persistente"]["NoCte"]
                        coeficientes_ELSF[i] = normas["IAPF"]["ELS"]["Desfavorable"]["Persistente"]["NoCte"]
                        coeficientes_ELSCP[i] = normas["IAPF"]["ELS"]["Desfavorable"]["Persistente"]["NoCte"]
                    elif restricciones_df.loc[i,"Direccion"] == "Favorable":
                        coeficientes_ELUP[i] = normas["IAPF"]["ELU"]["Favorable"]["Persistente"]["NoCte"]
                        coeficientes_ELUA[i] = normas["IAPF"]["ELU"]["Favorable"]["Accidental"]["NoCte"]
                        coeficientes_ELUS[i] = normas["IAPF"]["ELU"]["Favorable"]["Accidental"]["NoCte"]
                        #ELS
                        coeficientes_ELSC[i] = normas["IAPF"]["ELS"]["Favorable"]["Persistente"]["NoCte"]
                        coeficientes_ELSF[i] = normas["IAPF"]["ELS"]["Favorable"]["Persistente"]["NoCte"]
                        coeficientes_ELSCP[i] = normas["IAPF"]["ELS"]["Favorable"]["Persistente"]["NoCte"]
            elif restricciones_df.loc[i,"Tipo"] == "Accidental":
                if combo[i] != 0: n_acc += 1
                if restricciones_df.loc[i,"Direccion"] == "Desfavorable":
                    #ELU
                    coeficientes_ELUP[i] = 0
                    coeficientes_ELUA[i] = normas["IAPF"]["ELU"]["Desfavorable"]["Accidental"]["Accidental"]
                    coeficientes_ELUS[i] = normas["IAPF"]["ELU"]["Desfavorable"]["Accidental"]["Accidental"]
                    #ELS
                    coeficientes_ELSC[i] =  0
                    coeficientes_ELSF[i] =  0
                    coeficientes_ELSCP[i] = 0
                elif restricciones_df.loc[i,"Direccion"] == "Favorable":
                    #ELU
                    coeficientes_ELUP[i] = 0
                    coeficientes_ELUA[i] = normas["IAPF"]["ELU"]["Favorable"]["Accidental"]["Accidental"]
                    coeficientes_ELUS[i] = normas["IAPF"]["ELU"]["Favorable"]["Accidental"]["Accidental"]
                    #ELS
                    coeficientes_ELSC[i] =  0
                    coeficientes_ELSF[i] =  0
                    coeficientes_ELSCP[i] = 0
            elif restricciones_df.loc[i,"Tipo"] == "Variable":
                if combo[i] != 0: n_var += 1
                if restricciones_df.loc[i,"Direccion"] == "Desfavorable":
                    #ELU
                    coeficientes_ELUP[i] = normas["IAPF"]["ELU"]["Desfavorable"]["Persistente"]["Variable"]
                    coeficientes_ELUA[i] = normas["IAPF"]["ELU"]["Desfavorable"]["Accidental"]["Variable"]
                    coeficientes_ELUS[i] = normas["IAPF"]["ELU"]["Desfavorable"]["Accidental"]["Variable"]
                    #ELS
                    coeficientes_ELSC[i] = normas["IAPF"]["ELS"]["Desfavorable"]["Persistente"]["Variable"]
                    coeficientes_ELSF[i] = normas["IAPF"]["ELS"]["Desfavorable"]["Persistente"]["Variable"]
                    coeficientes_ELSCP[i] = normas["IAPF"]["ELS"]["Desfavorable"]["Persistente"]["Variable"]
                elif restricciones_df.loc[i,"Direccion"] == "Favorable":
                    #ELU
                    coeficientes_ELUP[i] = normas["IAPF"]["ELU"]["Favorable"]["Persistente"]["Variable"]
                    coeficientes_ELUA[i] = normas["IAPF"]["ELU"]["Favorable"]["Accidental"]["Variable"]
                    coeficientes_ELUS[i] = normas["IAPF"]["ELU"]["Favorable"]["Accidental"]["Variable"]
                    #ELS
                    coeficientes_ELSC[i] = normas["IAPF"]["ELS"]["Favorable"]["Persistente"]["Variable"]
                    coeficientes_ELSF[i] = normas["IAPF"]["ELS"]["Favorable"]["Persistente"]["Variable"]
                    coeficientes_ELSCP[i] = normas["IAPF"]["ELS"]["Favorable"]["Persistente"]["Variable"]
        elif norma_proyecto == "IAP":
            pass    
        
        
    minorados_ELUP.clear()
    minorados_ELUA.clear()
    minorados_ELUS.clear()
    minorados_ELSC.clear()
    minorados_ELSF.clear()
    minorados_ELSCP.clear()


    if n_acc:
        coef_var_ELUS = np.ones(len(patrones))
        coef_var_ELUA = np.ones(len(patrones))
        for i in range(len(combo)):
            if (restricciones_df.loc[i,"Accion IAPF"] == "Sismo") and combo[i] != 0:
                for j in range(len(combo)):
                    if (restricciones_df.loc[j,"Tipo"] == "Variable") and j:
                        if restricciones_df.loc[j,"Accion IAPF"] == "Trafico":
                            coef_var_ELUS[j] = normas["IAPF"]["CoefVar"]["Trafico"]["psi2"]
                        if restricciones_df.loc[j,"Accion IAPF"] == "Resto":
                            coef_var_ELUS[j] = normas["IAPF"]["CoefVar"]["Resto"]["psi2"]
                minorados_ELUS.append(coef_var_ELUS)
            if (restricciones_df.loc[i,"Accion IAPF"] != "Sismo") and (restricciones_df.loc[i,"Tipo"] == "Accidental") and combo[i] != 0:
                for k in range(len(combo)):
                    coef_var_ELUA = np.ones(len(patrones))
                    if restricciones_df.loc[k,"Tipo"] == "Variable":
                        if restricciones_df.loc[k,"Accion IAPF"] == "Trafico":
                            coef_var_ELUA[k] = normas["IAPF"]["CoefVar"]["Trafico"]["psi1"][vias_cargadas if int(vias_cargadas) < 4 else "3"]
                        else:
                            coef_var_ELUA[k] = normas["IAPF"]["CoefVar"]["Resto"]["psi1"]
                        for j in range(len(combo)):
                            if (restricciones_df.loc[j,"Tipo"] == "Variable") and (j != k):
                                    #se asigna el coeficiente de la variable al resto de variables
                                if restricciones_df.loc[j,"Accion IAPF"] == "Trafico":
                                    coef_var_ELUA[j] = normas["IAPF"]["CoefVar"]["Trafico"]["psi2"]
                                else:
                                    coef_var_ELUA[j] = normas["IAPF"]["CoefVar"]["Resto"]["psi2"]
                        minorados_ELUP.append(coef_var_ELUA)

    else:
        if n_var == 1:
            coef_var_ELSF = np.ones(len(patrones))
            coef_var_ELSCP = np.ones(len(patrones))
            for i in range(len(combo)):
                if (restricciones_df.loc[i,"Tipo"] == "Variable") and combo[i] != 0:
                    if restricciones_df.loc[i,"Accion IAPF"] == "Trafico":
                        coef_var_ELSF[i] = normas["IAPF"]["CoefVar"]["Trafico"]["psi1"][vias_cargadas if int(vias_cargadas) < 4 else "3"]
                        coef_var_ELSCP[i] = normas["IAPF"]["CoefVar"]["Trafico"]["psi2"]
                    else:
                        coef_var_ELSF[i] = normas["IAPF"]["CoefVar"]["Resto"]["psi1"]
                        coef_var_ELSCP[i] = normas["IAPF"]["CoefVar"]["Resto"]["psi2"]


        while n_var > 1:
            n_var -= 1 
            #se busca la primera variable en el combo y se asigna el coeficiente
            for i in range(len(combo)):
                coef_var_ELUP = np.ones(len(patrones))
                coef_var_ELSC = np.ones(len(patrones))
                coef_var_ELSF = np.ones(len(patrones))
                coef_var_ELSCP = np.ones(len(patrones))
                if (restricciones_df.loc[i,"Tipo"] == "Variable") and combo[i] != 0:
                    if restricciones_df.loc[i,"Accion IAPF"] == "Trafico":
                        coef_var_ELSF[i] = normas["IAPF"]["CoefVar"]["Trafico"]["psi1"][vias_cargadas if int(vias_cargadas) < 4 else "3"]
                        coef_var_ELSCP[i] = normas["IAPF"]["CoefVar"]["Trafico"]["psi2"]
                    else:
                        coef_var_ELSF[i] = normas["IAPF"]["CoefVar"]["Resto"]["psi1"]
                        coef_var_ELSCP[i] = normas["IAPF"]["CoefVar"]["Resto"]["psi2"]
                    for j in range(len(combo)):
                        if (restricciones_df.loc[j,"Tipo"] == "Variable") and (j != i):
                                #se asigna el coeficiente de la variable al resto de variables
                            if restricciones_df.loc[j,"Accion IAPF"] == "Trafico":
                                coef_var_ELUP[j] = normas["IAPF"]["CoefVar"]["Trafico"]["psi0"]
                                coef_var_ELSC[j] = normas["IAPF"]["CoefVar"]["Trafico"]["psi0"]
                                coef_var_ELSF[j] = normas["IAPF"]["CoefVar"]["Trafico"]["psi2"]
                                coef_var_ELSCP[j] = normas["IAPF"]["CoefVar"]["Trafico"]["psi2"]
                            else:
                                coef_var_ELUP[j] = normas["IAPF"]["CoefVar"]["Resto"]["psi0"]
                                coef_var_ELSC[j] = normas["IAPF"]["CoefVar"]["Resto"]["psi0"]
                                coef_var_ELSF[j] = normas["IAPF"]["CoefVar"]["Resto"]["psi2"]
                                coef_var_ELSCP[j] = normas["IAPF"]["CoefVar"]["Resto"]["psi2"]
                    minorados_ELUP.append(coef_var_ELUP)
                    minorados_ELSC.append(coef_var_ELSC)
                    minorados_ELSF.append(coef_var_ELSF)
                    minorados_ELSCP.append(coef_var_ELSCP)

    
    if minorados_ELUP:
        for minorado in minorados_ELUP: 
            combos_ELUP.append(np.multiply(np.multiply(coeficientes_ELUP, combo),minorado))
        for minorado in minorados_ELSC: 
            combos_ELSC.append(np.multiply(np.multiply(coeficientes_ELSC, combo),minorado))
        for minorado in minorados_ELSF: 
            combos_ELSF.append(np.multiply(np.multiply(coeficientes_ELSF, combo),minorado))
        for minorado in minorados_ELSCP: 
            combos_ELSCP.append(np.multiply(np.multiply(coeficientes_ELSCP, combo),minorado))
    else:
        combos_ELUP.append(np.multiply(coeficientes_ELUP, combo))
        combos_ELSC.append(np.multiply(coeficientes_ELSC, combo))
        combos_ELSF.append(np.multiply(coeficientes_ELSF, combo))if not n_var else combos_ELSF.append(np.multiply(np.multiply(coeficientes_ELSF, combo),coef_var_ELSF))
        combos_ELSCP.append(np.multiply(coeficientes_ELSCP, combo)) if not n_var else combos_ELSCP.append(np.multiply(np.multiply(coeficientes_ELSCP, combo),coef_var_ELSCP))
    if minorados_ELUS:
        for minorado in minorados_ELUS:
            combos_ELUS.append(np.multiply(np.multiply(coeficientes_ELUS, combo),minorado))
    if minorados_ELUA:
        for minorado in minorados_ELUA:
            combos_ELUA.append(np.multiply(np.multiply(coeficientes_ELUA, combo),minorado))
    


combos_ELUP = np.unique(np.array(combos_ELUP), axis=0)
combos_ELUA =  np.unique(np.array(combos_ELUA), axis=0)
combos_ELUS =  np.unique(np.array(combos_ELUS), axis=0)
combos_ELSC =  np.unique(np.array(combos_ELSC), axis=0)
combos_ELSF =  np.unique(np.array(combos_ELSF), axis=0)
combos_ELSCP =  np.unique(np.array(combos_ELSCP), axis=0)

print("------------------------------------------------")
print(combos_ELUP)
print(len(combos_ELUP))
print("------------------------------------------------")
print(combos_ELUA)
print(len(combos_ELUA))
print("------------------------------------------------")
print(combos_ELUS)
print(len(combos_ELUS))
print("------------------------------------------------")
print(combos_ELSC)
print(len(combos_ELSC))
print("------------------------------------------------")
print(combos_ELSF)
print(len(combos_ELSF))
print("------------------------------------------------")
print(combos_ELSCP)
print(len(combos_ELSCP))

total_combos = len(combos_ELUP) + len(combos_ELSC) + len(combos_ELSF) + len(combos_ELSCP) + len(combos_ELUA) + len(combos_ELUS)
#creacion de los combos en SAP2000

#creacion de los combos ELU Permantentes

combo_num = 0 

for combo in combos_ELUP:
    combo_name = "ELU_P_{:03d}".format(combo_num)
    ret = SapModel.RespCombo.Add(combo_name, 0)
    for idx, patron in enumerate(patrones):
        coeficiente = combo[idx]
        if coeficiente != 0:
            ret = SapModel.RespCombo.SetCaseList(combo_name, 0, patron, coeficiente)
    printProgressBar(combo_num, total_combos, prefix = 'Creando combos ELU Persistente:', suffix = 'Completado', length = 50)
    combo_num += 1

#creacion de los combos ELU Accidentales
tot = combo_num
combo_num = 0 

for combo in combos_ELUA:
    combo_name = "ELU_ACC_{:03d}".format(combo_num)
    ret = SapModel.RespCombo.Add(combo_name, 0)
    for idx, patron in enumerate(patrones):
        coeficiente = combo[idx]
        if coeficiente != 0:
            ret = SapModel.RespCombo.SetCaseList(combo_name, 0, patron, coeficiente)
    printProgressBar(combo_num, total_combos, prefix = 'Creando combos ELU Accidental:', suffix = 'Completado', length = 50)
    combo_num += 1

#creacion de los combos ELU Sismo
tot += combo_num
combo_num = 0
for combo in combos_ELUS:
    combo_name = "ELU_SIS_{:03d}".format(combo_num)
    ret = SapModel.RespCombo.Add(combo_name, 0)
    for idx, patron in enumerate(patrones):
        coeficiente = combo[idx]
        if coeficiente != 0:
            ret = SapModel.RespCombo.SetCaseList(combo_name, 0, patron, coeficiente)
    printProgressBar(tot+combo_num, total_combos, prefix = 'Creando combos ELU Sismo:', suffix = 'Completado', length = 50) 
    combo_num += 1

#creacion de los combos ELS Caracteristico
tot += combo_num
combo_num = 0 

for combo in combos_ELSC:
    combo_name = "ELS_C_{:03d}".format(combo_num)
    ret = SapModel.RespCombo.Add(combo_name, 0)
    for idx, patron in enumerate(patrones):
        coeficiente = combo[idx]
        if coeficiente != 0:
            ret = SapModel.RespCombo.SetCaseList(combo_name, 0, patron, coeficiente)
    printProgressBar(tot+combo_num, total_combos, prefix = 'Creando combos ELS Caracteristico:', suffix = 'Completado', length = 50) 
    combo_num += 1

#creacion de los combos ELS Frecuente
tot += combo_num
combo_num = 0 

for combo in combos_ELSF:
    combo_name = "ELS_F_{:03d}".format(combo_num)
    ret = SapModel.RespCombo.Add(combo_name, 0)
    for idx, patron in enumerate(patrones):
        coeficiente = combo[idx]
        if coeficiente != 0:
            ret = SapModel.RespCombo.SetCaseList(combo_name, 0, patron, coeficiente)
    printProgressBar(tot+combo_num, total_combos, prefix = 'Creando combos ELS Frecuente:', suffix = 'Completado', length = 50)
    combo_num += 1

#creacion de los combos ELS Cuasi Permanente
tot += combo_num
combo_num = 0 

for combo in combos_ELSCP:
    combo_name = "ELS_CP_{:03d}".format(combo_num)
    ret = SapModel.RespCombo.Add(combo_name, 0)
    for idx, patron in enumerate(patrones):
        coeficiente = combo[idx]
        if coeficiente != 0:
            ret = SapModel.RespCombo.SetCaseList(combo_name, 0, patron, coeficiente)
    printProgressBar(tot+combo_num, total_combos, prefix = 'Creando combos ELS Cuasipermanente:', suffix = 'Completado', length = 50)
    combo_num += 1
printProgressBar(total_combos, total_combos, prefix = 'Creando combos ELS Cuasipermanente:', suffix = 'Completado', length = 50)
#creacion del combo envolvente 
print()


combo_name = "1. ENV ELU P"
ret = SapModel.RespCombo.Add(combo_name, 1)
comboNames = SapModel.RespCombo.GetNameList()
patron = []
for i in range(comboNames[0]):
    if "ELU_P" in comboNames[1][i]: patron.append(comboNames[1][i])  
for i, name in enumerate(patron):
    ret = SapModel.RespCombo.SetCaseList(combo_name, 1, name, 1)
    printProgressBar(i+1, comboNames[0], prefix = 'Creando combo envolvente ELU P:', suffix = 'Completado', length = 50)
n=i

combo_name = "2. ENV ELU ACC"
ret = SapModel.RespCombo.Add(combo_name, 1)
comboNames = SapModel.RespCombo.GetNameList()
patron = []
for i in range(comboNames[0]):
    if "ELU_ACC" in comboNames[1][i]: patron.append(comboNames[1][i])  
for i, name in enumerate(patron):
    ret = SapModel.RespCombo.SetCaseList(combo_name, 1, name, 1)
    printProgressBar(n+i+1, comboNames[0], prefix = 'Creando combo envolvente ELU ACC:', suffix = 'Completado', length = 50)
n+=i

combo_name = "3. ENV ELU SIS"
ret = SapModel.RespCombo.Add(combo_name, 1)
comboNames = SapModel.RespCombo.GetNameList()
patron = []
for i in range(comboNames[0]):
    if "ELU_SIS" in comboNames[1][i]: patron.append(comboNames[1][i])  
for i, name in enumerate(patron):
    ret = SapModel.RespCombo.SetCaseList(combo_name, 1, name, 1)
    printProgressBar(n+i+1, comboNames[0], prefix = 'Creando combo envolvente ELU ACC:', suffix = 'Completado', length = 50)
n+=i

combo_name = "4. ENV ELS C"
ret = SapModel.RespCombo.Add(combo_name, 1)
comboNames = SapModel.RespCombo.GetNameList()
patron = []
for i in range(comboNames[0]):
    if "ELS_C_" in comboNames[1][i]: patron.append(comboNames[1][i])  
for i, name in enumerate(patron):
    ret = SapModel.RespCombo.SetCaseList(combo_name, 1, name, 1)
    printProgressBar(n+i+1, comboNames[0], prefix = 'Creando combo envolvente ELS C:', suffix = 'Completado', length = 50)
n+=i   


combo_name = "5. ENV ELS F"
ret = SapModel.RespCombo.Add(combo_name, 1)
comboNames = SapModel.RespCombo.GetNameList()
patron = []
for i in range(comboNames[0]):
    if "ELS_F_" in comboNames[1][i]: patron.append(comboNames[1][i])  
for i, name in enumerate(patron):
    ret = SapModel.RespCombo.SetCaseList(combo_name, 1, name, 1)
    printProgressBar(n+i+1, comboNames[0], prefix = 'Creando combo envolvente ELS F:', suffix = 'Completado', length = 50)
n+=i
    

combo_name = "6. ENV ELS CP"
ret = SapModel.RespCombo.Add(combo_name, 1)
comboNames = SapModel.RespCombo.GetNameList()
patron = []
for i in range(comboNames[0]):
    if "ELS_CP" in comboNames[1][i]: patron.append(comboNames[1][i])  
for i, name in enumerate(patron):
    ret = SapModel.RespCombo.SetCaseList(combo_name, 1, name, 1)
    printProgressBar(n+i+1, comboNames[0], prefix = 'Creando combo envolvente ELS CP:', suffix = 'Completado', length = 50)
printProgressBar(comboNames[0], comboNames[0], prefix = 'Creando combo envolvente ELS CP:', suffix = 'Completado', length = 50)

print()
print("Combinaciones de acciones creadas con éxito.")