This letter proposes a numerical method for approximating the location of and dynamics on a class of chaotic saddles. In contrast to the conventional strategy of maximizing the escape time, our proposal is to impose a zero-expansion condition along transversely repelling directions of chaotic saddles. This strategy exploits the existence of skeleton-forming unstable periodic orbits embedded in chaotic saddles, and thus can be conveniently implemented as a variant of subspace Newton-type methods. The algorithm is examined through an illustrative and another standard example.