bifurcation buckling

Hello every one.

I'm doing work on buckling  analysis using abaqus. I want to continue my work on postbuckling of plates also. I'm not getting any tutorials on how to perform Postbuckling analysis of plates using abaqus. When i searched the abaqus help documentation, i got .inp file of a plate, but i didnt get the steps on how to perform the analysis. 

I read in some books that after buckling of a plate, it still can take the load. How to do this in abaqus interface ?how to perform the postbuckling analysis in a proper way ?

Plz help

Does anyone know how to modify Bathe's subspace iteration eigensolver to compute the highest eigenvalue instead of the smallest ?

Anyone knows an efficient implementation of the subspace iteration method to compute a few lower eigenvalues and eigenvectors of a generalized problem KG*u=lambda*KS*u or any similar method? i.e., to compute the buckling loads and modes.

I'd like to use something optimized to Intel or amd64 processors or similar, like MKL, etc.

This is for bifurcation buckling analysis. KG is the stiffness matrix. KS is the geometric stress matrix. Both constructed from a FEM discretization.