avaliação da rigidez à flexão de toras de madeira por meio de ...

290
e com os graus dos polinômios, o programa calcula todos os termos da Matriz [A] da Equação
2.89. Na seqüência, são calculados os coeficientes do vetor {b} da mesma equação.
O programa calcula os termos do vetor {x} utilizando uma função para a resolução de
sistemas lineares denominada “SIS_LIN”. Essa função busca a solução de um sistema linear
por fatoração direta. Os termos dessa solução correspondem aos coeficientes dos polinômios
da Equação 2.77. Os parâmetros de entrada da função “SIS_LIN” são: a matriz [A], o vetor
{b}, a ordem da matriz, a matriz [Z] e o vetor {w}. Estes dois últimos parâmetros são
auxiliares nos cálculos. A matriz [A] é quadrada de ordem m+1+n.
Com os coeficientes do polinômio denominador da Equação 2.77, o programa calcula
as raízes ou pólos do polinômio usando a função “CNEWTON”. Esta função emprega o
método de Newton para a obtenção das raízes complexas do polinômio. Os parâmetros de
entrada são: o grau do polinômio (n) e os coeficientes do polinômio (bk).
A partir das raízes do polinômio denominador, o programa Modal-Id calcula o
amortecimento modal e as frequências naturais utilizando a Equação 2.92. O programa
desenha também o gráfico da FRF ajustada para a verificação visual da qualidade do ajuste.
Para tanto, o programa resolve a Equação 2.77 para cada valor de frequência do arquivo de
entrada.
Para auxiliar o usuário na verificação da qualidade do ajuste da FRF foi implementado
o cálculo do FDAC e do FRSF, respectivamente conforme Equações 2.96 e 2.99.
No botão “Salvar resultados” foi implementada a gravação dos resultados em um
arquivo texto como mostra a Figura C.4.
Vale destacar que não foi implementada nenhuma estratégia de eliminação dos modos
computacionais dos resultados. Portanto, cabe ao usuário verificar e descartar esses modos.
Entretanto, para facilitar a eliminação dos modos computacionais, foi inserida logo abaixo do
botão “Salvar resultados” uma caixa de texto na qual o usuário deverá informar qual é o
máximo amortecimento esperado para o sistema dinâmico testado. O programa irá ler esse
valor e eliminará todos os modos com amortecimentos superiores. O valor padrão adotado
para o amortecimento máximo foi de 2%.