As it is possible to use a buckling deformation as an initial imperfection for bars I would be pleased to see this function extended to shells. But not only:
It would be perfect if it was possible to deform permanently a mesh (bar and shell elements) according to a load case (buckling or not) in order to analyse imperfect models. This function is available in FEMAP.
Have a nice day,
Mathias.