We investigate the structure of gravitational self interactions of coherently oscillating axions in the general relativistic framework. A generic action for a massive scalar field in the Friedmann-Robertson-Walker background is first introduced based on the effective field theory approach to cosmological perturbations. Using the obtained setup, we evaluate the effective gravitational interaction of the massive scalar field, i.e. scalar quartic interactions mediated by metric perturbations. Applying the results to the system of dark matter axions, we estimate the self interaction rate of them, and discuss its implications for axion Bose-Einstein condensate dark matter scenario. Leading contributions for gravitational interactions of axions are given by the process mediated by the dynamical graviton field, which is essentially the Newtonian potential induced by fluctuations of the background fluids. We find that it leads to the same order of magnitude for the interaction rate of dark matter axions in the condensed regime, compared with the result of previous studies using the Newtonian approximation.