Программа может быть применена для любых размеров и числа витков катушки. В ней реализован метод суммирования. Приводятся результаты расчёта индуктивности катушки для полного осевого размера Ик и для расчётного hKp, равного эквивалентной глубине проникновения электромагнитного поля в катушку. Глубина проникновения поля рассчитана в п. 1.2.З.1. Задача 3.
PROGRAM IndKat IMPLICIT NONE.
INTEGER: nr, nWK, m, nrr, nrz, i INTEGER, ALLOCATABLE: NACHW (:), KONW (:).
REAL (8), ALLOCATABLE: RAD (:), Z (:).
REAL (8): rindkat, rKvn, tWr, tWz, shrzb, shr, shz, tizwwr NAMELIST/LIST/ nWK, rKvn, tWr, tWz, tizwwr, shrzb.
OPEN (7,FILE='lndKat.dat'); OPEN (8,FILE='lndKat.ch').
READ (7,list); WRITE)*, list); WRITE (8,list);
CALL RAZBIVKAftWr, tWz, shrzb, shr, shz, nrr, nrz, i) N-кол.площадок в витке nr=i*nWK !nrполное число разбиений; размер массивов RAD и Z ALLOCATEfRAD (nr), Z (nr), NACHW (nWK), KONW (nWK)) !Распр.массивов RAD (l)=rKvn+shr/2dO; Z (l)=shz/2d0 !Коорд-ты центра нач. площадки 1-го витка CALL KRD_PLOSK_KAT (RAD, Z, nr, 1, nrr, nrz, shz, shr, tizwwr, .false., nr).
NACHW (1)=1; KONW (l)=i do m=2, nWK.
NACHW (m)=NACHW (m-l)+i; KONW (m)=KONW (m-l)+i !Запомин. нач. и кон.номеров.
enddo
CALL lndKat_plosk_q (RAD, Z, nr, nWK, NACHW, KONW, nrr, nrz, shr, shz, rindkat) и/гКе (8,*)'Индуктивность катушки rindkat=', rindkat, Th'.
END PROGRAM IndKat.
Если в исходных данных задан осевой размер катушки, равный эквивалентной глубине проникновения поля SK= 2.26 мм (размер получен в п. 1.2.3.1. Задача 3; в файле данных в этом случае осевой размер катушки обозначается tWz=2.26d-3), то весь файл данных с именем lndKat.dat может быть записан с комментариями в виде:
&LIST
nWK=50 (число Витков Катушки rKvn=20d-3 (радиус Катушки внутренний tWr=ld-3 (толщина Витка радиальная tWz=2.26d-3 (толщина Витка осевая
tizwwr=ld-4 (толщина изоляции виток-виток радиальная shrzb=ld-3 (шаг разбиения /
Файл результатов IndKat.ch содержит следующие записи:
&LIST
NWK = 50,
RKVN = 2.0Е-002,
TWR = 1.0Е-003,
TWZ = 2.260 000 000 000 000Е-003,
TIZWWR = 1.0Е-004,
SHRZB = 1.0Е-003
Индуктивность катушки rindkat= 2.178 921 472 493 825Е-004 Гн
Для расчёта с осевым размером катушки 15 мм с исходными данными из п. 1.2.2.6. Задача 1. файл данных с именем lndKat.dat записан без комментариев и имеет вид:
&LIST nWk=50,rKvn=20d-3,tWr=ld-3,tWz=15d-3,tizwwr=ld-4,shrzb=ld-3 /
При работе программы в файл IndKat.ch помещаются прочитанные исходные данные (для уменьшения объёма пособия эти данные удалены) и результаты рассчёта в виде записи:
Индуктивность катушки rindkat= 1.900 918 679 588 754Е-004 Гн
Следует обратить внимание на тот факт, что чем меньше осевой размер катушки (при прочих равных размерах) гем болше её индуктивность.
