In this article, we investigate the structures of the pseudoscalar mesons (π, K, D, D s , B and B s ) in the framework of the coupled rainbow Schwinger-Dyson equation and ladder Bethe-Salpeter equation with the confining effective potential (infrared modified flat bottom potential). The Schwinger-Dyson functions for the u, d and s quarks are greatly renormalized at small momentum region and the curves are steep at about q 2 = 1 GeV 2 which indicates an explicitly dynamical symmetry breaking. The Euclidean time Fourier transformed quark propagators have no mass poles in the time-like region which naturally implements confinement. As for the c and b quarks, the current masses are very large, the renormalization are more tender, however, mass poles in the time-like region are also absent. The Bethe-Salpeter wavefunctions for those mesons have the same type (Gaussian type) momentum dependence and center around small momentum which indicate that the bound states exist in the infrared region. The decay constants for those pseudoscalar mesons are compatible with the values of experimental extractions and theoretical calculations, such as lattice simulations and QCD sum rules.