teste_incrementos.py

import time
import pyFEM

# inicia contagem do tempo
start_time = time.time()

# carrega o arquivo de malha do gmsh
mesh = pyFEM.pre.Malha("malha.msh")

# define e aplica o material
mat = pyFEM.materiais.Material(2e+7, 0.3, 0.1)
mesh.defineElementos(mat)

# aplica os apoios
mesh.Apoios(mesh.grupos['apoio'])

# aplica as forças
mesh.Forcas(mesh.linhas['topo'], [0, -100])

# inicializa o solver
solv = pyFEM.solver.Solver(mesh)

# resolve e exporta
solv.calcularInc(10, "malhaA")
umx = solv.U[1*2+1]*1000
print("Deslocamento incremental = %.3f mm\n" % umx)

# resolve e exporta
solv.calcular()

solv.exportar("malhaB")

# tempo de processamento
elapsed_time = time.time() - start_time

# deslocamento vertical do nó inferior da extreminde livre (Nó #1)
umx = solv.U[1*2+1]*1000
print("Deslocamento normal = %.3f mm\n" % umx)