Using a sample of (2712.4±14.3)×106 ψ(3686) events collected by the BESIII detector at the BEPCII collider, we measure the branching fractions of the decays χcJ→φφη, φφη′, and~φK+K−η (J=0,1,2). The obtained branching fractions are B(χc0→φφη)=(7.40±0.23±0.55)×10−4, B(χc1→φφη)=(3.33±0.14±0.25)×10−4, B(χc2→φφη)=(5.46±0.17±0.40)×10−4, B(χc0→φφη′)=(2.96±0.23±0.29)×10−4, B(χc1→φφη′)=(0.69±0.10±0.08)×10−4, B(χc2→φφη′)=(0.65±0.09±0.07)×10−4, B(χc0→φK+K−η)=(1.23±0.08±0.10)×10−4, B(χc1→φK+K−η)=(1.00±0.07±0.07)×10−4, and B(χc2→φK+K−η)=(1.82±0.09±0.14)×10−4, where K+K− is not from the decay of a φ meson, the first uncertainties are statistical and the second systematic. The branching fractions of χcJ→φφη are measured with precision improved by factors of 1.5−1.9, and those of χcJ→φφη′ and φK+K−η are measured for the first time.