We present new micromagnetic simulations, where an influence of the wire shape on the spatial distribution of internal magnetoelastic stress is taken into account. Local magnetoelastic anisotropy is approximated to be uniaxial, with spatial dependence due to the stress distribution. The hysteresis loop is calculated with the mumax 3 software for the wire volume and separately for selected parts of the wire, including the surface. The bistability of the calculated surface hysteresis loops depends on the behaviour of the neighboring part of the inner core, most close to the surface element where the surface loop is measured.