Translate this page:
Please select your language to translate the article


You can just close the window to don't translate
Library
Your profile

Back to contents

Cybernetics and programming
Reference:

On the use of computational experiments in the study of physics

Maier Robert Valerievich

Doctor of Pedagogy

Professor, the department of Physics and Didactics of Physics, Glazov State Pedagogical Training Institute

427628 Russia, The Udmurt Republic, Glazov, Kalinina Street 8A, unit #79

robert_maier@mail.ru
Other publications by this author
 

 

DOI:

10.7256/2306-4196.2014.6.13483

Received:

02-11-2014


Published:

16-11-2014


Abstract: The author reviews a problem of use of educational computational experiments in the study of physical phenomena. By educational computational experiment author means experiment on the mathematical model of the object carried out with the aid of a computer for the purpose of learning. A set of simplified versions of educational computational experiments adapted to the conditions of learning form a System of educational computational experiments. The article analyzes the examples of the use of educational computer experiment for: 1) for the study of the magnetization of a ferromagnet, calculation of the magnetization curve and hysteresis loop; 2) study of chaotic oscillations of a Dafing pendulum, occurrence of bifurcation in changes in the profile of the potential well, study of Poincaré section and the evolution of the phase volume. The authors apply mathematical and computational methods (simulation) modeling, which involves the construction of a mathematical model and building a software simulating the phenomenon under study based on the numerical solution of the corresponding to its’ system of equations. The novelty of the work is in the proposed four simple computer programs in Pascal, allowing: 1) obtain the magnetization curve and hysteresis loop for a ferromagnet in a changing magnetic field; 2) simulate the oscillation of the Dafing pendulum; 3) study the transition of oscillator into the chaotic regime at a change of the potential well; 4) get a Poincaré section and to study the evolution of the phase volume for Dafing pendulum.


Keywords:

information technology, computer simulation, programming, computer simulations, physical phenomena, education, magnetization of a ferromagnet, magnetic hysteresis, Dafing pendulum, dynamic chaos


1. Введение

Развитие информационных технологий создало предпосылки для использования компьютерных моделей при изучении физических явлений [1-6]. Имеет смысл говорить об учебном вычислительном эксперименте, под которым следует понимать упрощенные варианты научного вычислительного эксперимента, адаптированные к условиям обучения. Фактически это эксперимент над математической моделью объекта, проводимый с помощью ЭВМ с целью обучения. В процессе выполнения УВЭ учащиеся изменяют параметры модели исследуемой системы, характер и величину внешних воздействий, начальные условия и "наблюдают", как при этом изменяется отклик системы, скорость и направление протекающих процессов. Это позволяет изучить динамику изменения различных величин, характеризующих изучаемое явление, сформировать его наглядный образ, повысить интерес студентов к физике. УВЭ следует рассматривать как дополнение к учебной теории и учебному натурному эксперименту и использовать в сочетании с ними. Рассмотрим два примера использования УВЭ для изучения физических систем.

2. Компьютерное моделирование намагничивания ферромагнетика

При изучении явления намагничивания ферромагнетика обычно рассматривают следующий опыт. Ферромагнитный стержень вставляют в обмотку из n витков, которая подключена к источнику тока. Увеличивают силу тока I в обмотке, измеряют индукцию B магнитного поля. Затем уменьшают силу тока, пропускают ток в противоположном направлении, увеличивая и затем уменьшая его до нуля. Используя результаты измерения H и B, строят кривую намагничивания B(H). Для более глубокого понимания физической сущности этого явление можно решить следующую задачу.


Задача 1. Ферромагнитный стержень находится внутри обмотки, через которую пропускают переменный ток. Методами компьютерного моделирования изучите процесс намагничивания ферромагнетика, получите кривую намагниченности и петлю гистерезиса. Выясните, как изменятся ее форма при уменьшении амплитуды колебаний тока через обмотку.


Рассмотрим двумерную модель ферромагнетика, состоящую из прямоугольной сетки N x M (кристаллической решетки), в узлах которой находятся атомы с магнитным моментом p. Сначала магнитное поле отсутствует, ферромагнетик полностью размагничен, магнитные моменты атомов ориентированы произвольно, образуя с осью Ox случайные углы alfa[i,j]. Если ось Ox сонаправлена с силовой линией магнитного поля, то среднее значение индукции:

fo1

где H=nI – напряженность магнитного поля. Индукция магнитного поля Bij в данном узле сетки (i, j) обусловлена внешним полем обмотки H и магнитным полем, создаваемым соседними атомами. В проекции на ось Ox:

fo2_01
Со стороны магнитного поля на атом, находящийся в узле (i, j), действует вращающий момент MВР=Bij*p*sin(alfa[i,j]). При повороте атома возникает момент “упругих сил”, препятствующий повороту. В результате действия этих моментов, а также сил внутреннего трения, атом поворачивается на k*MВР, и угол alfa[i,j] становится равным alfa[i,j]+k*MВР+x. Здесь x – небольшая случайная величина в интервале [–0,1; 0,1], учитывающая хаотические колебания атомов в узлах кристаллической решетки. Давая случайное приращение x углу alfa[i,j], мы как бы “встряхиваем” нашу модель, предоставляем возможность найти более “удачные” значения alfa[i,j]. Это важно еще и потому, что в анализируемой модели число атомов существенно меньше, чем в реальном ферромагнитном стержне.


Суть вычислительного эксперимента состоит в плавном увеличении напряженности магнитного поля и пересчете alfa[i,j] у всех атомов с нахождением соответствующих средних значений индукции магнитного поля BSR. В результате на плоскости B от H получится кривая намагниченности с характерным горизонтальным участком, соответствующим магнитному насыщению, когда магнитные моменты всех атомы ориентированы по полю. Для получения петли гистерезиса необходимо, чтобы напряженность поля H совершила одно гармоническое колебание.

risu1_01


Представленная ниже программа ПР–1 содержит цикл по времени t, в котором пересчитываются H, alfa[i,j], BSR и строится зависимость B(H). Шаг по времени должен быть достаточно мал. После ее запуска на экране рисуется петля гистерезиса (рис. 1.1). При уменьшении амплитуды колебаний напряженности поля H получается петля гистерезиса, изображенная на рис. 1.2. Аналогичным образом можно промоделировать поляризацию сегнетоэлектриков в изменяющемся электрическом поле.

pro1_02

Задача 2. С помощью компьютерной модели исследуйте зависимость формы петли гистерезиса от величины магнитного момента p атомов ферромагнетика.


Для решения этой задачи также используется программа ПР–1. Если уменьшить магнитный момент атомов p в 2 - 3 раза (от 0,15 до 0,05), то уменьшится индукция магнитного поля, соответствующая насыщению, петля гистерезиса сожмется по вертикали (рис. 1.3).

3. Моделирование хаотических колебаний маятника Дафинга

Примером простой механической системы, обнаруживающей способность к хаотическому движению, является осциллятор Дафинга. Он представляет собой частицу, движущуюся в потенциальной яме с двумя углублениями. Рассмотрим несколько задач, которые могут быть решены при изучении хаотического движения маятника Дафинга.

Задача 3. Промоделируйте вынужденные колебания маятника Дафинга, на который действует внешняя гармонически изменяющаяся сила.


Рассмотрим колебания шарика внутри потенциальной ямы с двумя углублениями под действием внешней периодически изменяющейся силы F(t). Пусть потенциальная энергия внутри ямы задается функцией U(x). Потенциальное поле создает возвращающую силу F'(t). Профиль потенциальной ямы, уравнение вынуждающей и возвращающей сил, а также дифференциальное уравнение колебаний маятника Дафинга представлены ниже:

fo3
Результат решения задачи представлен на рис. 2. Из графиков x(t) и v(t) и фазовой кривой видно, что система совершает нерегулярные стохастические колебания относительно двух положений равновесия. Движение шарика является трудно предсказуемым, бесконечно малые изменения параметров системы, начальных условий и внешней силы приводят к тому, что система эволюционирует по иному пути.

risu2_01

Задача 4. Получите сечения Пуанкаре для вынужденных колебаний маятника Дафинга, происходящих под действием внешней периодически изменяющейся силы. Для этого рассмотрите фазовую кривую в пространстве, образованном осями x, p и F, и получите множество точек ее пересечения с плоскостью F=const.


Результаты вычислений приведены на рис. 3. Сечения Пуанкаре для маятника Дафинга получены так: на ЭВМ моделируются колебания системы и определяются ее состояния в моменты времени, когда фаза вынуждающей силы равна 0, 3,14/4, 3,14/2, 3*3,14/4, 3,14. Так как режим хаотический, точки сечения имеют фрактальную структуру. Программа ПР-2 представлена ниже.

pro2

risu3_03

Задача 5. Промоделируйте переход к хаосу в случае, когда на маятник Дафинга действует гармоническая сила при изменяющемся профиле потенциальной ямы. Получите точку бифуркации, в которой происходит раздвоение пути эволюции системы.


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


На рис. 4.1 изображен профиль потенциальной ямы U(x)=0,25*x4-b*x2 при различных значениях параметра b1=0,1, b2=0,5, b3=1. Видно, что при увеличении b в центре ямы растет бугорок и образуются два симметричных углубления.

risu4_01

Необходимо многократно (1000 раз) моделировать вынужденные колебания маятника Дафинга при случайной начальной фазе вынуждающей силы, а на экран компьютера выводить координаты x через заданное время после начала колебаний, когда они уже установились. Одновременно с этим должен изменяться бифуркационный параметр b, определяющий профиль потенциальной ямы (рис. 4.1). При удачном подборе параметров модели когда параметр b начинает превышать критическое значение bk=0, вместо одной потенциальной ямы получается две, и система начинает совершать колебания относительно одного из двух положений равновесия. На экране ставится точка, показывающая координату x в заданный момент t, и все повторяется снова. Начальная фаза силы изменяется случайным образом, поэтому шарик оказывается то в левом, то в правом углублении. В результате получается бифуркация типа вил (рис. 4.2).

pro3

Задача 6. Промоделируйте перемешивание фазового объема в случае свободных незатухающих колебаний маятника Дафинга.


Представим себе совокупность из 160 одинаковых маятников Дафинга, совершающих свободные незатухающие колебания (r=0), отличающиеся только начальными условиями x, p. Пусть в момент t=0 фазовые точки, характеризующие начальное состояние маятников, находятся внутри прямоугольника [x, x+dx, p, p+dp]. Предлагаемая программа рассчитывает состояние каждого маятника в момент времени t1, и ставит соответствующую точку на фазовой плоскости. Результаты вычислительного эксперимента представлены на рис. 5. Видно, что происходит расползание фазового объема, его перемешивание в фазовом пространстве, что свидетельствует о хаотичности колебаний. Так как колебания незатухающие, то в соответствии с теоремой Лиувилля величина фазового объема остается постоянной.

pro4

risu5_01

4. Заключение

В настоящей статье предложены примеры использования учебного вычислительного эксперимента для изучения намагничивания ферромагнетика и хаотических колебаний маятника Дафинга. Рассмотренные программы позволяют: 1) получить кривую намагниченности и петлю гистерезиса для ферромагнетка в изменяющемся магнитном поле; 2) промоделировать колебания маятника Дафинга; 3) изучить переход осциллятора в хаотический режим при изменении профиля потенциальной ямы; 4) получить сечение Пуанкаре и изучить эволюцию фазового объема для маятника Дафинга. Рассмотренные компьютерные модели могут быть использованы в качестве лекционных демонстраций, на практических занятиях и лабораторных работах при изучении курса физики и основ компьютерного моделирования. Кроме того, они могут быть частью выпускной квалификационной работы или исследовательского проекта. Их применение способствует повышению интереса студентов к численным методам, компьютерным моделям и информационным технологиям в целом.

References
1. Guld Kh. Tobochnik Ya. Komp'yuternoe modelirovanie v fizike: V 2 ch. / Kh. Guld, Ya. Tobochnik. – M.: Mir, 1990. – Ch.2. – 400 s.
2. Kunin S. Vychislitel'naya fizika. –– M.: Mir, 1992. –– 518 s.
3. Maier R.V. Zadachi, algoritmy, programmy. [Elektronnyi resurs] –– Glazov, 2011.-URL: http://maier-rv.glazov.net.
4. Maier R.V. Ispol'zovanie vychislitel'nykh eksperimentov pri izuchenii volnovykh protsessov v lineinykh i nelineinykh sredakh // NB: Kibernetika i programmirovanie.-2014.-№4.-C. 57-65. DOI: 10.7256/2306-4196.2014.4.12683. URL: http://www.e-notabene.ru/kp/article_12683.html.
5. Maier R.V. Komp'yuternoe modelirovanie fizicheskikh yavlenii. –– Glazov: GGPI, 2009. –– 112 s. (Rezhim dostupa: maier-rv.glazov.net).
6. Samarskii A.A., Mikhailov A.P. Matematicheskoe modelirovanie: Idei. Metody. Primery. –– M.: Fizmatlit, 2001. –– 320 s.