This paper is based on the hypothesis that introducing distribution of mechanical properties is beneficial for modeling all kinds of mechanical behavior, even of ordinary metallic materials. To bring proof of its admissibility, it has to be first shown that modeling based on this assertion is able to efficiently describe standard mechanical behavior of materials. Searching for typical study case, it has been assessed that at a low scale, yield stresses could be strongly distributed in ultrathin aluminum foils used in packaging industry, offering opportunities to identifying their distribution and showing its role on the mechanical properties. Considering initially reduced modeling allow to establish a valuable connection between the hardening curve and the distribution of local yield stresses. This serves for finding initial value of distribution parameters in a more sophisticated identification procedure. With finally limited number of representative classes of local yield stresses, concretely 3 is enough, it is shown that a 3D finite element simulation involving limited numbers of elements returns realistic behavior of an ultrathin aluminum foil exerted to tensile test, in reference to experimental results. This gives way to large possibilities in modeling in order to give back complex experimental evidence. © 2018 Author(s).