In this paper, the critical clearing time, tcc in power system transient stability analysis is modeled as a random variable due to the randomness nature of power system load. A linear approximation approach that mainly involves the sensitivity calculation of the derivative of tcc to system load is developed to obtain the probability density function (pdf) of tcc based on the pdf of system load. This approach is verified to be accurate under the condition of small load deviation by Monte Carlo simulation method. The probability of system being transiently stable for particular disturbances are calculated based on the pdf of tcc.