Several useful Green's functions are derived for the quasi-static analysis of shielded planar transmission lines by the boundary element method. These newly employed Green's functions satisfy forced boundary conditions in a rectangular region. The integral equation does not have a singular point and the integrand contains only the normal derivative of the electric potential. The present method is proposed to characterize multilayered and multi-conductor structures. Numerical results are presented for a microstrip, a suspended line, and a coplanar waveguide. For the coplanar waveguide, a combined Green's function is also studied. This combined Green's function further reduces the memory size in computation. All of these Green's functions are represented in infinite series. The resulting matrix equation has slowly convergent matrix elements. To reduce the computation time of matrix elements, we split the original series in two parts. The geometric series method is employed to convert one part into a fast convergent form (4 terms). For the other part, only a few terms (less than 20) require computation.