In this paper, a new framework for accurate reliability analysis is proposed based on improving the directional simulation by using metaheuristic algorithms. Usually for highly nonlinear and complex performance functions, finding the unit vector direction requires very high calculations or impossible practically. Hence, the novel improved version incorporates the Harris Hawks Optimization algorithm, where the unit vector of direction is formulated as a constrained optimization problem and estimated using metaheuristic algorithms. Given that metaheuristic algorithms have been introduced to solve unconstrained problems, the penalty function method is used to convert the constrained problem into an unconstrained problem. The applicability of the proposed framework is firstly tested on five highly nonlinear benchmark functions and then applied to solve four high-dimensional engineering problems. The performance of six simulations-based reliability analysis methods and the first-order reliability method were compared with the proposed method. Besides the feasibility of other metaheuristic algorithms were investigated. The results show high-performance abilities of the improved version of the directional simulation for solving highly nonlinear engineering problems.