Analisis Estructural Dinamico Avanzado

[ ]

rigidez_piso = 45000.0
mostrar "=== Analisis Estructural Dinamico Avanzado ==="
mostrar "Simulacion real en el tiempo (no carga estatica equivalente):"
mostrar "integracion directa de las ecuaciones modales acopladas bajo"
mostrar "un sismo sintetico, con 5% de amortiguamiento critico -- lo"
mostrar "que de verdad se hace en un analisis tiempo-historia profesional."
mostrar ""
mostrar "Rigidez de entrepiso actual (kN/m):"
mostrar rigidez_piso

masa_piso = 80.0
altura_piso = 3.2

K = [[0,0,0,0,0],[0,0,0,0,0],[0,0,0,0,0],[0,0,0,0,0],[0,0,0,0,0]]
K[1,1] = rigidez_piso * 2
K[1,2] = 0 - rigidez_piso
K[2,1] = 0 - rigidez_piso
K[2,2] = rigidez_piso * 2
K[2,3] = 0 - rigidez_piso
K[3,2] = 0 - rigidez_piso
K[3,3] = rigidez_piso * 2
K[3,4] = 0 - rigidez_piso
K[4,3] = 0 - rigidez_piso
K[4,4] = rigidez_piso * 2
K[4,5] = 0 - rigidez_piso
K[5,4] = 0 - rigidez_piso
K[5,5] = rigidez_piso
M = [[masa_piso,0,0,0,0],[0,masa_piso,0,0,0],[0,0,masa_piso,0,0],[0,0,0,masa_piso,0],[0,0,0,0,masa_piso]]

mostrar ""
mostrar "--- Motor espectral: frecuencias naturales ---"
analisis_modal(K, M, 100, frecuencias, modos)
mostrar frecuencias[1]
mostrar frecuencias[2]
mostrar frecuencias[3]
mostrar frecuencias[4]
mostrar frecuencias[5]
graficar_barras(frecuencias, "Frecuencias naturales (rad/s)")

mostrar ""
mostrar "--- Factor de participacion modal ---"
Gamma = [0,0,0,0,0]
para i = 1 hasta 5
    suma_num = 0.0
    suma_den = 0.0
    para j = 1 hasta 5
        suma_num = suma_num + modos[j,i]
        suma_den = suma_den + modos[j,i]*modos[j,i]
    finpara
    Gamma[i] = suma_num / suma_den
finpara

mostrar ""
mostrar "--- Simulacion tiempo-historia (integracion directa, 400 pasos) ---"
mostrar "Sismo sintetico: envolvente gaussiana modulando una onda de 2Hz"
zeta_amort = 0.05
dt_paso = 0.005
n_pasos = 400
paso_muestreo = 4

q_modal = [0,0,0,0,0]
qp_modal = [0,0,0,0,0]

dpiso1xx = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]
dpiso2xx = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]
dpiso3xx = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]
dpiso4xx = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]
dpiso5xx = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]
tiempo_muestra = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]

idx_muestra = 0
contador_sub = 0
qfund_hx = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]
qpfund_h = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]

para k = 1 hasta n_pasos
    t_val = k * dt_paso
    diff_t = t_val - 1.0
    exponente = 0 - (diff_t*diff_t)/(2*0.5*0.5)
    envolvente = exp(exponente)
    ag_val = 2.5 * envolvente * seno(2*3.14159265*2.0*t_val)

para i = 1 hasta 5
        qpp = 0 - 2*zeta_amort*frecuencias[i]*qp_modal[i] - frecuencias[i]*frecuencias[i]*q_modal[i] - Gamma[i]*ag_val
        qp_modal[i] = qp_modal[i] + qpp*dt_paso
        q_modal[i] = q_modal[i] + qp_modal[i]*dt_paso
    finpara

contador_sub = contador_sub + 1
    si contador_sub == paso_muestreo
        contador_sub = 0
        idx_muestra = idx_muestra + 1

u1 = 0.0
        u2 = 0.0
        u3 = 0.0
        u4 = 0.0
        u5 = 0.0
        para i = 1 hasta 5
            u1 = u1 + modos[1,i]*q_modal[i]
            u2 = u2 + modos[2,i]*q_modal[i]
            u3 = u3 + modos[3,i]*q_modal[i]
            u4 = u4 + modos[4,i]*q_modal[i]
            u5 = u5 + modos[5,i]*q_modal[i]
        finpara

dpiso1xx[idx_muestra] = u1
        dpiso2xx[idx_muestra] = u2
        dpiso3xx[idx_muestra] = u3
        dpiso4xx[idx_muestra] = u4
        dpiso5xx[idx_muestra] = u5
        tiempo_muestra[idx_muestra] = t_val
        qfund_hx[idx_muestra] = q_modal[5]
        qpfund_h[idx_muestra] = qp_modal[5]
    finsi
finpara

mostrar "Desplazamiento maximo del techo durante el sismo (m):"
maximo_techo = 0.0
para i = 1 hasta 100
    val_abs = valor_abs(dpiso5xx[i])
    si val_abs > maximo_techo
        maximo_techo = val_abs
    finsi
finpara
mostrar maximo_techo

Historia_pisos = [[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0],[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0],[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0],[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0],[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]]
para j = 1 hasta 100
    Historia_pisos[1,j] = dpiso1xx[j]
    Historia_pisos[2,j] = dpiso2xx[j]
    Historia_pisos[3,j] = dpiso3xx[j]
    Historia_pisos[4,j] = dpiso4xx[j]
    Historia_pisos[5,j] = dpiso5xx[j]
finpara
graficar_multi_matriz(Historia_pisos, "Respuesta dinamica de cada piso durante el sismo (m)")

dispersion_3d(qfund_hx, qpfund_h, tiempo_muestra, "Retrato de fase del modo fundamental (posicion, velocidad, tiempo) -- la espiral hacia el reposo por el amortiguamiento")

Salida Jota: === Analisis Estructural Dinamico Avanzado === Salida Jota: Simulacion real en el tiempo (no carga estatica equivalente): Salida Jota: integracion directa de las ecuaciones modales acopladas bajo Salida Jota: un sismo sintetico, con 5% de amortiguamiento critico -- lo Salida Jota: que de verdad se hace en un analisis tiempo-historia profesional. Salida Jota: Salida Jota: Rigidez de entrepiso actual (kN/m): Salida Jota: 45000.000 Salida Jota: Salida Jota: --- Motor espectral: frecuencias naturales --- Salida Jota: 45.513 Salida Jota: 39.904 Salida Jota: 31.063 Salida Jota: 19.705 Salida Jota: 6.751 Salida Jota: Salida Jota: --- Factor de participacion modal --- Salida Jota: Salida Jota: --- Simulacion tiempo-historia (integracion directa, 400 pasos) --- Salida Jota: Sismo sintetico: envolvente gaussiana modulando una onda de 2Hz Salida Jota: Desplazamiento maximo del techo durante el sismo (m): Salida Jota: 0.035

Shift + Enter ejecuta la celda actual

Ayuda de Jota