Neste trabalho, propõe-se um modelo numérico-computacional representativo dos processos de tratamento térmico, que seja uma ferramenta eficiente e forneça meios para um entendimento efetivo do mecanismo de geração de tensões residuais durante a têmpera de aços. Foram investigados os fenômenos térmicos, mecânicos e de transformação de fase observados na têmpera, bem como o acoplamento entre esses três fenômenos. O modelo utiliza o Método dos Elementos Finitos (MEF) e o programa ABAQUS®, além de rotinas numéricas em FORTRAN responsáveis pela resolução do problema termo-mecânico-microestrutural acoplado. A utilização de sub-rotinas, que implementam uma alteração na formulação (matemática e numérica) do programa de Elementos Finitos, permite incluir no modelo as informações presentes em uma curva CRC (curva de resfriamento contínuo) do aço SAE 4140, implementando o cálculo de deformações da peça simulada de modo incremental e cumulativo. Os resultados mostram que a utilização das sub-rotinas desenvolvidas neste trabalho permitiu implementar, conjuntamente com o programa ABAQUS®, o cálculo das frações volumétricas, durezas, distorções e tensões que surgem em um tratamento térmico de têmpera, simulando as transformações martensítica, perlítica, bainítica e ferrítica. Os resultados dos modelos foram equivalentes aos relatados pela literatura, principalmente no que se refere às durezas e tensões associadas a cada transformação de fase. Em particular, os resultados indicam que a transformação martensítica está sempre associada à formação de tensões compressivas. Ensaios experimentais foram realizados a fim de validar os modelos computacionais propostos, utilizando-se um teste Jominy adaptado e instrumentado, de modo a permitir a amostragem da variação de temperaturas no material. Ensaios metalográficos permitiram correlacionar as frações volumétricas transformadas durante a têmpera do corpo de prova Jominy aos valores calculados pelo modelo numérico acoplado.

The objective of this work is to analyze residual strains and stresses and volumetric expansion due to phase transformations that occur during quenching of a steel body, as well as to predict these phase transformations. The coupled thermo-mechanical-phase transformation problem was analyzed, specifically in terms of the quenching process. Different computational models were presented, based on the finite element software ABAQUS® and on the use of FORTRAN subroutines. The continuous-cooling-transformation (CCT) diagrams of SAE 4140 steel are represented differently in each model, depending on the transformed phases and correspondent volumetric expansion. These subroutines include information from the CCT diagrams of SAE 4140 into a FORTRAN code. The subroutine calculates all the microstructures resulting from quenching (ferrite, pearlite, bainite, and martensite), depending on cooling rate. The numerical analysis conducted in this work provided results in terms of the temperature and stresses developed during quenching. The properties determined in this work are hardness, yield strength, volumetric fraction and distortion. Hardness has been predicted by the use of analytical equations. The finite element analyses were able to explain and reproduce phenomena observed during quenching of a steel cylinder. In particular, numerical results indicated that martensite formation is always related to a compressive stress field. The results of the models are in qualitative agreement with data provided by literature, particularly, in relation to the stresses originated by each different phase transformation during quenching process. Experimental testing was conducted, based on the analysis of the quenching of a Jominy probe, in order to validate the computational model developed in this work.