Weak interactions of solitary waves in the generalized nonlinear Schrödinger (NLS) equations are studied. It is first shown that these interactions exhibit similar fractal dependence on initial conditions for different nonlinearities. Then by using the Karpman-Solov'ev method, a universal system of dynamical equations is derived for the velocities, amplitudes, positions, and phases of interacting solitary waves. These dynamical equations contain a single parameter, which accounts for the different forms of nonlinearity. When this parameter is zero, these dynamical equations are integrable, and the exact analytical solutions are derived. When this parameter is nonzero, the dynamical equations exhibit fractal structures which match those in the original wave equations both qualitatively and quantitatively. Thus the universal nature of fractal structures in the weak interaction of solitary waves is analytically established.