The present research intended to propose and evaluate regression models which estimate toughness property as a function of physical, chemical properties of thermally modified hardwood and thermal treatment temperature, using linear, quadratic, cubic, exponential, logarithmic, geometric and multiple linear models. Commercial thermally modified woods were used on the study, being characterized for all referred properties, totalizing 450 experimental determinations. The analyzed models presented a low and moderate coefficient of determination, indicating the impossibility to use such models in the estimation of toughness as a function of physical and chemical factors