A stochastic finite element method is developed for the buckling analysis of frames with random initial imperfections, uncertain sectional and material properties. The random geometrical imperfections of the frames are described by member initial crookednesses which are modeled as given initial displacement functions with amplitudes treated as random variables. The effects of the random initial geometric imperfections are formulated as a set of equivalent random nodal coordinates in the finite element discretization of the members. The mean-centered second-order perturbation technique is used to formulate the stochastic finite element method for the buckling analysis of the imperfect frames. Use of the present method is illustrated by several examples of buckling analysis of random frames. Results derived from the Monte Carlo method are also obtained for comparison.