Different modes of tsunami edge waves can interact through nonlinear resonance. During this process, edge waves that have very small initial amplitude can grow to be as large or larger than the initially dominant edge wave modes. In this study, the effects of dynamical phase are established for a single triad of edge waves that participate in resonant interactions. In previous studies, Jacobi elliptic functions were used to describe the slow variation in amplitude associated with the interaction. This analytical approach assumes that one of the edge waves in the triad has zero initial amplitude and that the combined phase of the three waves φ = θ₁ + θ₂ − θ₃ is constant at the value for maximum energy exchange (φ = 0). To obtain a more general solution, dynamical phase effects and non-zero initial amplitudes for all three waves are incorporated using numerical methods for the governing differential equations. Results were obtained using initial conditions calculated from a subduction zone, inter-plate thrust fault geometry and a stochastic earthquake slip model. The effect of dynamical phase is most apparent when the initial amplitudes and frequencies of the three waves are within an order of magnitude. In this case, non-zero initial phase results in a marked decrease in energy exchange and a slight decrease in the period of the interaction. When there are large differences in frequency and/or initial amplitude, dynamical phase has less of an effect and typically one wave of the triad has very little energy exchange with the other two waves. Results from this study help elucidate under what conditions edge waves might be implicated in late, large-amplitude arrivals.