GPa = 1e+9
MPa = 1e+6
KN = 1e3
m = 1.0
#Inicializar modelo
ret = Reticulado()
#Nodos
ret.agregar_nodo(0,0)
ret.agregar_nodo(1,0)
ret.agregar_nodo(1,1)
print(ret)
#Barras
b1 = Barra(0, 1, 20*cm, 4*mm, 200*GPa, 7600*kg/m**3, 420*MPa)
b2 = Barra(1, 2, 20*cm, 4*mm, 200*GPa, 7600*kg/m**3, 420*MPa)
b3 = Barra(0, 2, 20*cm, 4*mm, 200*GPa, 7600*kg/m**3, 420*MPa)
ret.agregar_barra(b1)
ret.agregar_barra(b2)
ret.agregar_barra(b3)
peso_total = ret.calcular_peso_total()
print(f"peso_total = {peso_total}")
ver_reticulado_2d(ret)
ret.agregar_barra(Barra(0, 1, *props))
ret.agregar_barra(Barra(1, 2, *props))
ret.agregar_barra(Barra(3, 4, *props))
ret.agregar_barra(Barra(0, 3, *props))
ret.agregar_barra(Barra(3, 1, *props))
ret.agregar_barra(Barra(1, 4, *props))
ret.agregar_barra(Barra(4, 2, *props))
ret.agregar_restriccion(0, 0, 0)
ret.agregar_restriccion(0, 1, 0)
ret.agregar_restriccion(2, 1, 0)
ret.agregar_fuerza(4, 1, -F)
ret.ensamblar_sistema()
ret.resolver_sistema()
f = ret.recuperar_fuerzas()
print(ret)
ver_reticulado_2d(ret,
opciones_nodos = {
"usar_posicion_deformada": True,
"factor_amplificacion_deformada": 10.,
},
opciones_barras = {
"color_barras_por_fuerza": True,
"ver_numeros_de_barras": False,
"ver_fuerza_en_barras": True
})