The thin layer problem is common in some engineering cases, such as a coating target or dielectric radome. If the skin depth of the medium behind the thin layer is much smaller than the media scale, the transmission line impedance boundary can be used to solve. When both sides of the thin layer are free space or the skin depth of the medium behind is greater than or comparable to the scale of the medium, the transmission line impedance boundary condition is not applicable. Such problems require incident and reflected wave data and are also known as penetrable thin-layer problems. When dealing with penetrable thin-layer problems, we need to give the electromagnetic wave data on both sides. In addition, the thin layer problem can be attributed to multi-scale problems. The EM simulation is difficult to implement because of the high resource consumption of accurate modeling. Considering that the accurate modeling is very expensive, and the mesh-size difference between the thin layer region and the surrounding ordinary region is too large will seriously affect the quality of the elements (tetrahedrons). It is difficult to solve quickly in the time-domain discontinuous Galerkin method even by using local time steps or other optimization methods. To get a fast solution, the approximate boundary condition can be used to describe the thin layer, DGTD also has the advantage of the geometric fitting property. If a thin layer model is curved, we can build a conformal geometric model to deal with the thin layer, thus avoiding the complicated interpolation correction scheme. In this paper, A Robust DGTD Method for penetrable thin layers is proposed. In our study, we find that the current excitation scheme is prone to divergence when dealing with large dielectric coefficients and its robustness is poor. To solve this problem, we use the electromagnetic field correction scheme in DGTD to deal with the thin layer problems. In addition, engineering problems may encounter multiple thin layers. In this paper, a simulation scheme for multiple thin layers is proposed and relevant numerical examples are given as verification. The results illustrate the accuracy and effectiveness of the method in this paper, and the method in this paper can be used for the fast calculation of complex thin layer problems.