The linear spectrum and corresponding Bloch modes of shallow honeycomb lattices near Dirac points are investigated. Via perturbation theory, the dispersion relation is found to have threefold degeneracy at leading order with eigenvalue splitting at the following two orders; i.e., the threefold eigenvalue splits into single and double values. Multiscale perturbation methods are
employed to describe the nonlinear dynamics of the associated wave envelopes. The dynamics of the envelope depends on different asymptotic balances whereupon a three-level nonlinear Dirac-type equation or a two-level nonlinear Dirac equation is derived. The analysis agrees well with direct numerical simulations.