This paper introduces an improved response surface-based fuzzy finite element analysis of structural dynamics. The free vibration of structure is established using superposition method, so that fuzzy displacement responses can be presented as functions of fuzzy mode shapes and fuzzy natural frequencies. Instead of direct determination of these fuzzy quantities by modal analysis which will involve the calculation of the whole finite element model, the paper proposes a felicitous approach to design the response surface as surrogate model for the problem. In the design of the surrogate model, complete quadratic polynomials are selected with all fuzzy variables are transformed to standardized fuzzy variables. This methodology allows accurate determination of the fuzzy dynamic outputs, which is the major issue in response surface based techniques. The effectiveness of the proposed fuzzy finite element algorithm is illustrated through a numerical analysis of a linear two-storey shear frame structure.