In order to solve the problem of photonic crystal band structure more conveniently and quickly

based on the COMSOL Multiphysics coefficient partial differential equation（ PDE） module of the finite element numerical simulation software

the mathematical model is reestablished. The partial differential form of the photonic crystal is derived from the Helmholtz equation and combined with the Bloch state. The eigenfunction

and the propagation of Bloch wave vector is taken into full consideration

the corresponding eigenfrequency is solved so as to solve the energy band. In this case

the energy band structure of one and two-dimensional photonic crystals is simulated respectively. The simulation results are compared with the results of traditional methods

and the feasibility of the scheme is verified.