HB: BCSSTK01 — Rigidez estructural — vibraciones modale [RSA 48×48]

[ ]

; Harwell-Boeing — BCSSTK01 (RSA)
; 48x48 NNZ=224
valores = [2832268.518519999925047,1000000.0,2083333.333330000052229,-3333.333333330000187,1000000.0,-2800000.0,-28935.185185199999978,2083333.333330000052229,1635447.530860000057146,-2000000.0,5555555.555549999698997,-6666.666666670000268,-2000000.0,-30864.197530900000856,5555555.555549999698997,-1597916.666669999947771,1724367.283950000070035,-2083333.333330000052229,-2777777.777780000120401,-1680000.0,-15432.098765400000048,-2777777.777780000120401,-28935.185185199999978,-2083333.333330000052229,1003333333.330000042915344,2000000.0,400000000.0,-3333333.333329999819398,2083333.333330000052229,100000000.0,1067500000.0,-1000000.0,200000000.0,2777777.777780000120401,333333333.333000004291534,-833333.333332999958657,1535333333.329999923706055,-2000000.0,-5555555.555549999698997,666666666.66700005531311,-2083333.333330000052229,100000000.0,2832268.518519999925047,-1000000.0,2083333.333330000052229,-2800000.0,-28935.185185199999978,2083333.333330000052229,1635447.530860000057146,2000000.0,5555555.555549999698997,-30864.197530900000856,5555555.555549999698997,-1597916.666669999947771,1724367.283950000070035,-2083333.333330000052229,-2777777.777780000120401,-15432.098765400000048,-2777777.777780000120401,-28935.185185199999978,-2083333.333330000052229,1003333333.330000042915344,-3333333.333329999819398,2083333.333330000052229,100000000.0,1067500000.0,2777777.777780000120401,333333333.333000004291534,-833333.333332999958657,1535333333.329999923706055,-5555555.555549999698997,666666666.66700005531311,-2083333.333330000052229,100000000.0,2836099.469500000122935,-2149285.294509999919683,2359161.804020000156015,-3333.333333330000187,-1000000.0,-28935.185185199999978,2083333.333330000052229,-3830.950981710000178,-1149285.294509999919683,275828.470682999992277,1767410.744459999958053,517922.131815999979153,4298570.589019999839365,-5555555.555549999698997,-6666.666666670000268,2000000.0,-1597916.666669999947771,-131963.213599000009708,-517922.131815999979153,2298570.589019999839365,3890038.068479999899864,-2634990.274699999950826,2777777.777780000120401,-1680000.0,-28935.185185199999978,-2083333.333330000052229,-517922.131815999979153,-2165670.784529999829829,-551656.941366999992169,1975720635.309999942779541,-2000000.0,400000000.0,2083333.333330000052229,100000000.0,-2298570.589019999839365,551656.941365999984555,486193650.990000009536743,1527346515.470000028610229,-109779731.332000002264977,1000000.0,200000000.0,-833333.333332999958657,1149285.294509999919683,229724661.236000001430511,-55717351.077899999916553,1564111437.109999895095825,-2000000.0,-2083333.333330000052229,100000000.0,-275828.470682999992277,-55717351.077899999916553,10941196.003799999132752,2832268.518519999925047,1000000.0,2083333.333330000052229,-28935.185185199999978,2083333.333330000052229,1635447.530860000057146,-2000000.0,-5555555.555549999698997,-1597916.666669999947771,1724367.283950000070035,-2083333.333330000052229,2777777.777780000120401,-28935.185185199999978,-2083333.333330000052229,1003333333.330000042915344,2083333.333330000052229,100000000.0,1067500000.0,-833333.333332999958657,1535333333.329999923706055,-2083333.333330000052229,100000000.0,60879.629629600000044,1250000.0,416666.666666999983136,-4166.666666670000268,1250000.0,3372916.666670000180602,-2500000.0,-8333.333333329999732,-2500000.0,2411712.962960000149906,-416666.666666999983136,-2355000.0,1500000000.0,2500000.0,500000000.0,501833333.333000004291534,-1250000.0,250000000.0,502500000.0,-2500000.0,3985879.629629999864846,-1250000.0,416666.666666999983136,-3925000.0,3411496.913579999934882,2500000.0,6944444.444439999759197,-38580.246913599999971,6944444.444450000301003,2431003.086420000065118,-416666.666666999983136,-3472222.222219999879599,-19290.123456799999985,-3472222.222219999879599,1504166666.670000076293945,-4166666.666670000180602,1335166666.670000076293945,3472222.222219999879599,416666666.666999995708466,2169166666.670000076293945,-6944444.444439999759197,833333333.33299994468689,3985879.629629999864846,-1250000.0,416666.666666999983136,-4166.666666670000268,-1250000.0,3411496.913579999934882,2500000.0,-6944444.444450000301003,-8333.333333329999732,2500000.0]
suma_cuad = 0.0
i = 1
para i de 1 hasta 200
    suma_cuad = suma_cuad + valores[i] * valores[i]
finpara
frob_m = raiz(suma_cuad)
frob_t = raiz(suma_cuad * 1.12)
mostrar "=== BCSSTK01: 48x48 NNZ=224 ==="
mostrar "Tipo: RSA — Rigidez estructural"
mostrar "Sparsidad (%):"
mostrar 9.7222
mostrar "Frob Jota (muestra 200/224):"
mostrar frob_m
mostrar "Frob estimada total:"
mostrar frob_t
mostrar "Frob referencia: 7291988925.4839"
mostrar "Max valor:"
mostrar maximo(valores)
mostrar "Min valor:"
mostrar minimo(valores)
graficar_histograma(valores, 20, "BCSSTK01 valores")

Salida Jota: === BCSSTK01: 48x48 NNZ=224 === Salida Jota: Tipo: RSA — Rigidez estructural Salida Jota: Sparsidad (%): Salida Jota: 9.722 Salida Jota: Frob Jota (muestra 200/224): Salida Jota: 6064502937.498 Salida Jota: Frob estimada total: Salida Jota: 6418066639.137 Salida Jota: Frob referencia: 7291988925.4839 Salida Jota: Max valor: Salida Jota: 2169166666.670 Salida Jota: Min valor: Salida Jota: -109779731.332

Shift + Enter ejecuta la celda actual

Ayuda de Jota