The current study develops a numerical model to investigate the nanoindentation behavior of heat-affected zones (HAZ) on glass material produced via the electrochemical spark machining (ECSM) method. Initially, microchannels were created using the ECSM method on soda–lime glass. Following that, a nanoindentation test was conducted to quantify the Young’s modulus and hardness of the glass sample. After that, a numerical model based on finite elements was created to characterize the changes in mechanical characteristics of HAZ. According to the findings, increasing the electrolyte concentration from 10 to 30% increases the intensity of electrochemical discharges, and thereby decreases the hardness of the work material by 16.29 to 30.58% compared to unmachined glass. The results obtained from the simulation are in close agreement with the experimental values. The maximum error obtained between simulation and experimental results is only 4.18%.