Изменил номиналы R-C на 120 кОм и 47наноФарад, токовый резистор в первичке 60 Ом, обмотки - по 100 витков, поставил разделительный конденсатор на выходе УНЧ . Вроде более-менее стало рисовать картинку
http://savepic.org/4249575.jpg
http://savepic.org/4242407.jpg

В качестве конденсатора интегратора можно использовать только плёночные конденсаторы

Вторичную обмотку образца нужно сделать с таким числом витков, чтобы оптимально "вписаться" в полный размах индукции на экране осциллографа при соответствующих "милливольтах" на интеграторе макета.
( вроде так по логике нужно сделать ? )

Подумал, по поводу снятия характеристик В-Н .
Насколько правильно понимаю саму "философию" этого вопроса.
Считаю так:
1. Обязательным условием является синусоидальность напряжения на обмотке.

2. При синусоидальном напряжении известно, что индукция в магнитопроводе тоже синусоидальна и сдвинута на 90 градусов относительно ЭДС.
3. Таким образом, для получения сигнала, подобного ! индукции ( но не являющегося непосредственным отражением индукции в магнитопроводе !!!! ) , используется ФНЧ в виде R-C цепочки, для получения сдвига сигнала ЭДС на 90 градусов.

4. Правильно для случая 10кГц , рассчитать ФНЧ на частоту не выше 100 Гц . При этом выполняется условие приближения фазового сдвига к 90 градусам, а амплитуда сигнала на интеграторе составляет порядка 0,01 от входного напряжения. Далее, для настройки , можно понижать частоту среза ниже 100Гц изменением номиналов R-C , при этом условие "фазового" сдвига будет лучше, но существенно снижается амплитуда сигнала наблюдаемой "индукции" .

5. Необходимо учитывать, что при малых значениях R в интеграторе будет излишне большой ток, что будет вносить искажения , а при высоких значениях R скажется влияние входного сопротивления осциллографа.

Попробовал в модели посравнивать снимаемые B-Н при условии синусоидального сигнала, и при условии существенных искажений сигнала.
Модели рисуют одинаковые "петли" при условии подбора напряжения генератора, для создания одинаковой "индукции" . Это даёт основание полагать, что исследовать несоответствие реальной индукции и напряжения интегратора - задача достаточно сложная. Но те переплетения восходящей и убывающей ветви на В-Н, которые обозначил рисунком на прошлой странице, не спроста.... "неточность " интегрирования ???

система работает для любой формы периодического переменного напряжения, не имеющего постоянной составляющей.
Учитывая то, что производная индукции пропорциональна приложенному напряжению, делённому на сечение и витки
dB/dt=V/(A*N)
Получить текущее значение индукции (при учёте её начального значения) можно интегрируя функцию напряжения.

Расчётное значение индукции - порядка 0,36 Тл . ( *** надо ещё подумать и перепроверить, по формуле которая фигурировала выше получается около 0,251 Тл , - результат неверный ).

Возможно чего-то не учёл в формуле. Это не долго поправить. Однако с выводами торопиться не стоит. Надо насытить материал, т.е. довести напряженность до 1200 А/м.

Я так понял, надо отрисовать семейство кривых, чтобы получить "начальную" кривую по точкам вершин ? Это нам потом понадобится.