Recent studies have shown that the preimage set partitions of weakly regular bent functions, particularly those that are vectorial dual–bent, can give rise to association schemes. The first construction of association schemes from non-weakly regular bent functions, namely from ternary generalized Maiorana–McFarland bent functions from \({\mathbb F}_{3^n} \times {\mathbb F}_{3^k} \times {\mathbb F}_{3^k}\) to \({\mathbb F}_3\) for \(n = 1\) and \(n = 2\) , is presented in Özbudak and Pelen (J. Algebraic Combin. 56 (2022), 635–658). This construction was substantially generalized to arbitrary odd primes p and positive integers n in a recent work of Anbar et al. (Finite Fields Appl. 103 (2025), Paper No. 102568), using a variant of generalized Maiorana–McFarland bent functions. In this paper, we refine the construction of Anbar et al. to obtain association schemes on \(\mathbb {F}_{p^{n}} \times \mathbb {F}_{p^{k}} \times \mathbb {F}_{p^{k}}\) with \( \left( \frac{p^{k} - 1}{p - 1}(p + 1) + p \right) , \quad \left( \frac{p^{k} - 1}{p - 1}(p + 1) + p - 1 \right) \; \text {and} \quad \left( \frac{p^{k} - 1}{p - 1}(p + 1) + \frac{p - 1}{2} \right) \) association classes, depending on n and on the choice of bent functions employed in the construction. We further emphasize that the association schemes previously obtained in the works of Özbudak–Pelen and Anbar et al. arise as fusion schemes of those constructed in this paper.
