The generalized Polynomial Chaos Expansion Method (gPCEM), which is a random uncertainty analysis method by employing the orthogonal polynomial bases from the Askey scheme to represent the random space, has been widely used in engineering applications due to its good performance in both computational efficiency and accuracy. But in gPCEM, a nonlinear transformation of random variables should always be used to adapt the generalized Polynomial Chaos theory for the analysis of random problems with complicated probability distributions, which may introduce nonlinearity in the procedure of random uncertainty propagation as well as leading to approximation errors on the probability distribution function (PDF) of random variables. This paper aims to develop a flexible polynomial expansion method for response analysis of the finite element system with bounded random variables following arbitrary probability distributions. Based on the large family of Jacobi polynomials, an Improved Jacobi Chaos Expansion Method (IJCEM) is proposed. In IJCEM, the response of random system is approximated by the Jacobi expansion with the Jacobi polynomial basis whose weight function is the closest to the probability density distribution (PDF) of the random variable. Subsequently, the moments of the response can be efficiently calculated though the Jacobi expansion. As the IJCEM avoids the necessity that the PDF should be represented in terms of the weight function of polynomial basis by using the variant transformation, neither the nonlinearity nor the errors on random models will be introduced in IJCEM. Numerical examples on two random problems show that compared with gPCEM, the IJCEM can achieve better efficiency and accuracy for random problems with complex probability distributions.