A simulação computacional em biomecânica permite analisar o comportamento dos movimentos do corpo humano, diminuindo, e inclusive evitando ensaios experimentais invasivos. A locomoção humana resulta das forças desenvolvidas pelos músculos esqueléticos. Os mecanismos que produzem essas forças ainda são um tema de investigação aberto. O pouco entendimento desse fenômeno tem levado a subestimar propriedades importantes nos modelos mecânicos, as quais são essenciais para a simulação do comportamento do músculo. O objetivo desta tese foi desenvolver um código computacional que permita obter de maneira precisa e exata, a representação numérica do comportamento mecânico dos músculos esqueléticos. O código visa compilar diversas pesquisas numéricas de tal forma que a simulação possa considerar os fenômenos essenciais no comportamento mecânico do músculo e posteriormente avaliar sua influência na geração de força muscular. A formulação utilizada é baseada no Método dos Elementos Finitos (MEF), que é escrito em função das posições nodais. Os músculos esqueléticos foram discretizados por elementos planos e sólidos e uma análise não linear geométrica foi realizada. O programa considera fibras longas colocadas dentro de um domínio contínuo (passivo) sem adicionar graus de liberdade ao sistema). Um modelo transversalmente isotrópico, hiperelástico quase incompressível foi utilizado para simular o tecido muscular. A energia livre de Helmholtz foi usada para modelar o comportamento muscular ativo e passivo do músculo. Os resultados da pesquisa mostram que o código computacional é adequado para representar um modelo hiperelástico quase incompreensível no modelo transversalmente isotrópico. Permitindo considerar o músculo esquelético em duas partes distintas: intramuscular (matriz) e extracelular (fibras) utilizando a energia livre de Helmholtz e com ativação uniaxial, tanto em modelos estáticos como dinâmicos não lineares. Os resultados numéricos demonstraram que o algoritmo implementado é adequado para realizar análises não lineares geométricas de músculos esqueléticos via MEF. A condição de incompressibilidade foi comprovada nos problemas com materiais hiperelásticos. Também, foi demostrada a necessidade de realizar uma análise de convergência para as fibras. Finalmente, foi notada a complexidade na construção e na análise estrutural dos músculos esqueléticos, sendo necessário continuar desenvolvendo estratégias numéricas para maior aprofundamento.

Computational Modeling in Biomechanics allows analyzing of human body's movements, decreasing and some cases avoiding invasive experimental tests. The human locomotion is the result of forces developed by skeletal muscles. The mechanisms that produce this force are still an open research topic. The little knowledge of this phenomenon has led to underestimating important properties in mechanical models. The goal of this thesis was developed a computer code to obtain, in a precise and exact manner, the numerical representation of the mechanical behavior of skeletal muscles. The code aims to compile several numerical research, such that the simulation can consider the essential phenomena in mechanical behavior and then evaluate their influence in the muscle strength development. The used formulation is based on the Finite Element Method (FEM), which is written as a function of nodal positions. The skeletal muscles were discretized by plane and solid elements, and a geometrically nonlinear analysis was performed. The program considers long fibers placed inside a continuum domain (passive) without adding degrees of freedom to the system. A transversely isotropic model almost incompressible hyperelastic model was used to simulate the muscle tissue. The Helmholtz free energy was used to model the active and passive muscle behavior of muscle. The findings from the research indicate that the computer code is adequate to represent a transversely isotropic model almost incompressible hyperelastic model. The code allows skeletal muscle to be considered in two parts: intramuscular (matrix) and extracellular (fibers) using the Helmholtz free energy and with uniaxial activation, in nonlinear statical and dynamical models. The results support the model implemented for nonlinear geometrical analyzes of skeletal muscle using FEM. The almost incompressibility condition was tested in problems with hyperelastic materials. Also, numerical simulations confirm that a convergence analyzes for fibers is always required. Finally, it was noted the complexity in the construction and the structural analyzes of skeletal muscles, being necessary to continue developing numerical strategies for further deepening.