This paper presents a novel procedure based on first-order reliability method (FORM) for structural reliability analysis with hybrid variables, that is, random and interval variables. This method can significantly improve the computational efficiency for the abovementioned hybrid reliability analysis (HRA), while generally providing sufficient precision. In the proposed procedure, the hybrid problem is reduced to standard reliability problem with the polar coordinates, where an n-dimensional limit-state function is defined only in terms of two random variables. Firstly, the linear Taylor series is used to approximate the limit-state function around the design point. Subsequently, with the approximation of the n-dimensional limit-state function, the new bidimensional limit state is established by the polar coordinate transformation. And the probability density functions (PDFs) of the two variables can be obtained by the PDFs of random variables and bounds of interval variables. Then, the interval of failure probability is efficiently calculated by the integral method. At last, one simple problem with explicit expressions and one engineering application of spacecraft docking lock are employed to demonstrate the effectiveness of the proposed methods.