Assessing structural reliability problem with high-dimensional random inputs is still challenging due to the “curse of dimensionality”. In this paper, this challenge is addressed by extending the Generalized Distribution Reconstruction Method via Characteristic Function Inversion (GDRM-CFI). Specifically, a clustering-based partially stratified sampling method is proposed for selecting high-dimensional points to numerically evaluate the characteristic function (CF) curve of complex high-dimensional problems. An improved number-theoretical method (i-NTM) is used to establish a uniform, efficient point set, ensuring determinism and reducing variability. Subsequently, a partial stratification approach partitions the high-dimensional space into orthogonal two-dimensional subspaces. The fundamental point set is projected into each subspace, and the k-means clustering algorithm identifies centroids within each, acting as representative points. The complete set of representative points from all subspaces formulates the high-dimensional point set. Numerical examples are investigated, which demonstrate the proposed method is effective for high-dimensional structural reliability assessment. •The proposed method can generate high-dimensional samples with great uniformity.•The proposed method is applicable to high-dimensional reliability assessment.•The proposed method outperforms some widely-used high-dimensional sampling methods.•Numerical examples validate the effectiveness of the proposed method.