The thermodynamic topology optimization for hardening materials [1] is a novel approach to consider hardening material behavior into the optimization process in a resource-efficient manner due to a surrogate model. Here, we present our numerical implementation in Julia programming language [2]. Using this Julia code, we optimized a quasi 2D clamped beam, a quasi 2D classical Messerschmitt-Bölkow-Blohm (MBB) beam and a 3D cantilever as numerical results in [1]. We also provide the used job files for the three examples here. Different material parameters of hardening can be investigated due to the parameters in the configuration file. The given zip folder includes the Julia code organized as Julia project and additionally a directory with the job files (configuration, mesh and boundary value problem) for the mentioned examples. [1] Miriam Kick and Philipp Junker. Thermodynamic topology optimization for hardening materials. Submitted. (2022) [2] Jeff Bezanson, Alan Edelman, Stefan Karpinski and Viral B Shah. Julia: A fresh approach to numerical computing, www.julialang.org. SIAM Review, 59(1):65–98. (2017)