Главная Развитие электроэнергетической системы [0] [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [ 20 ] [21] [22] [23] [24] [25] [26] [27] [28] [29] [30] [31] [32] [33] [34] [35] [36] [37] [38] [39] [40] [41] [42] [43] [44] [45] [46] [47] [48] [49] [50] [51] [52] [53] [54] [55] [56] [57] [58] [59] [60] [61] [62] [63] [64] [65] [66] [67] [68] [69] [70] [71] [72] [73] [74] [75] [76] [77] [78] [79] [80] [81] [82] [83] [84] рис. 2.18. К расчету максимальной напряженности электрического поля непроволочного галетного силового резистора: а - РЭ; б - контактный узел; / - резистивный материал; 2 - контакт; 3 -изоляция РЗ; в - разбиение границы РЭ на граничные элементы; г - геометрия граничного элемента u(x) =i>il p(f)«o(f. х)с?Г(П); г Эх,- (x) (2.54) (2.55) где \[/ - функция, обратная зависимости электрической проводимости от потенциала, так что \[/(а(и(х))) =и(х). №тегральное уравнение для определения р(х) имеет вид Хз(х) = J р(Г) Г + (Г) X,(r)Uo(f, x)-h Э« (x) rfr(f). (2.56) Xi (?) X2(f) 1 на Ti; 0 на Гг иГз; О на Г,; = < на Гг; 1 на Гз; Хз(?)= < (ф{0) на Тг, ф{и12) на Гз; Ю на Гз- Воспользуемся аксиальной симметрией задачи и снизим ее размерность до двух, для чего перейдем к цилиндрическим координатам Xl =лсо8в; Х2 = rsinO; Х3 = z. Интеграл от первого члена в квадратных скобках уравнения (2.56) можно переписать в виде /i(x)= J (2.57) Интегрирование в (2.56) ведется на поверхности РЭ, а (2.57) пред ставляет собой криволинейный интеграл; К (к) - эллиптический интеграл первого рода, модуль которого равен: где A-jj, л , zj- - цилиндрические координаты точек х и f. При стремлении точки f к х мотуль эллиптического интеграла стремится к единице, а сам он может быть представлен в виде К{к) = 1 In т.е. ядро /1 (х) имеет логарифмическую сингулярность. Для вычисления второго интеграла в квадратных скобках необходимо вычислить производные ядра по ruz. Обозначим ядро через гК(к) и после несложных алгебраических преобразований вычислим ЪЗ/Ъг и Э5/Эг [29]: dS . VCf - ix- f) [Cx--f)- (x- f)] При стремлении f к x производные dS/dr и dS/dz имеют также логарифмические особенности. Для численного решения задачи (2.56) граница РЭ разбивается на ряд граничных элементов М (рис. 2.18, в). Число граничных элементов определяется точностью решения поставленной задачи. Рассмотрим криволинейный граничный элемент, изображенный на рис. 2.18, г. функцию р можно выразить в векторной форме с помощью однородной координаты rj в виде [27] р(?) = («1, ct2, аз) /Р1\ V Рз/ где надстрочный индекс t означает транспонирование ai -0,5j?(j?-1); = 0,5ч(?+1); «3= (1 -»?)(1 + j?)- В свою очередь координаты границы /-го граничного элемента можно записать в виде ,(Л) = а(г1) I \ 22/ ; rf(rj)=a(n) г21 VlJ \"з Интеграл по /-му граничному элементу Iji4) = Pj J «(»?) Jg(?) С (»?) + Э5 Эг 95 95 9г 9»? 9Z 9»? (2.58) где G= y/{dr/dv) + (9z/9j?)/ = 1,2,..., л/; а: = 1,2,..., л/. Вычисление интегралов осуществляется с помощью квадратур Гаусса 2 / = 1 (2.59) В (2.59) под Z (»?) понимается подьштегральная функция в выражении (2.58). В частности, если х не лежит на данном конечном элементе и и = 3, то координаты узлов интегрирования и значение весов будут те, что приведены в табл. 2.1. 5-6319 [0] [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [ 20 ] [21] [22] [23] [24] [25] [26] [27] [28] [29] [30] [31] [32] [33] [34] [35] [36] [37] [38] [39] [40] [41] [42] [43] [44] [45] [46] [47] [48] [49] [50] [51] [52] [53] [54] [55] [56] [57] [58] [59] [60] [61] [62] [63] [64] [65] [66] [67] [68] [69] [70] [71] [72] [73] [74] [75] [76] [77] [78] [79] [80] [81] [82] [83] [84] 0.0011 |