The paper presents an approach developed on the basis of the functional voxel method to the geometric representation of the thermal expansion of objects and temperature stresses in a material when exposed to a surface of a heat source. A discrete geometric law of a single temperature stress in an isotropic heat-conducting body is derived, applicable in the concept of functional voxel modeling. Based on this law, functional-voxel models of thermal stress are developed for a single and distributed application of a heat source. Algorithms of functional-voxel modeling of temperature stress and expansion in the case of distributed thermal loading are presented, which make it possible to construct a loading region of a complex configuration, uniformly form a contour (surface) after material expansion and obtain information about changes in the length (volume) of products. The advantages of the proposed functional-voxel approach to modeling thermal expansion and stress over approaches based on the FEM are substantiated.