Дано развитие метода [1] численного решения граничной задачи электродинамики для произвольных проводящих поверхностей на случай присутствия произвольных диэлектрических тел. Учет диэлектрических тел осуществлен посредством добавления к исходной системе интегральных уравнений дополнительных членов с неизвестным током поляризации, использовании техники параметрического отображения для представления криволинейного диэлектрического тела и метода Галеркина с граничными и конечными элементами. Объемный ток поляризации внутри диэлектрического тела представлен в виде криволинейных конечных элементов с кусочно-постоянными базисными функциями. Развито параметрическое описание геометрии криволинейного обобщенного шестигранника для представления широкого класса диэлектрических тел. Метод и реализующая его программа апробированы на примере расчета в десятикратной полосе частот входного импеданса поверхностных излучателей с плоской и криволинейной структурами, снабженных диэлектрической подложкой конечных размеров. Показано, что рассчитанные этим методом зависимости входного импеданса от частоты и значений диэлектрической проницаемости подложки хорошо коррелируют с референтными результатами, полученными по широко апробированным программам HFSS и CST Microwave Studio
