CC BY
261
УДК 621.1.36.7 (035.5) В.Ф. ОчКов
Задачи по физике: новый подход к решению
В статье рассмотрены вопросы использования современных программных средств для решения классических задач физики.
Ключевые слова: физика, Mathcad, анимация.
PHYSIC PROBLEMS: A NEW WAY TO SOLUTIONS
The paper deals with the use of advanced software solutions for the classical problems ofphysics.
Keywords: physics, Mathcad, animation.
История развития методов решения задач по математике, физике, химии и другим школьным и вузовским учебным дисциплинам - это, помимо прочего, и история борьбы с... вычислительными средствами. Сначала (на уроках устного счета, например) запрещали пользоваться бумагой и ручкой, затем (при изучении счета «столбиком» ручкой на бумаге) - калькулятором и, наконец, при решении более сложных задач - компьютером с современными математическими программами: МаШетайса, Маріє, МайаЬ, Mathcad, SMath и др. Вернее, прямо не запрещали и не запрещают, а говорят в том плане, что, мол, решение школьной задачи по физике в среде, например, программы Mathcad равносильно решению задачи для устного счета (65 помножить на 9, например) на калькуляторе. Этот запрет касается не самих вычислительных средств, которые можно и нужно осваивать с помощью специально подобранных примеров на занятиях по информатике, а применения этих средств для решения задач, придуманных для уроков «математики, физики, химии» задолго до появления компьютеров и калькуляторов. Но современные школьники и студенты этого не понимают. Более того, они уже не мыслят учебы, а иногда, к сожалению, и досуга без компьютера. Еще один «способ» решения задач с использованием современных информационных технологий
таков. Задача отправляется по SMS эксперту или через интернет вывешивается на каком-либо форуме. Слово «способ» тут можно писать и без кавычек. Почему?! Преподаватель может задать школьникам или студентам задачу и попросить решить ее всеми возможными способами, включая и тот, какой мы отметили выше, - с использованием современных средств связи: «звонок другу». На форуме вышеотме-ченных математических программ www.exponenta.ru часто можно увидеть такие «звонки».
С устным здесь счетом все более-менее ясно. Эти упражнения - прекрасная гимнастика для ума. Использование на таких занятиях калькулятора равносильно дооборудованию спортивного тренажера (гимнастика для тела)... гидроусилителями. Счет «столбиком» также можно рассматривать как гимнастику для ума. Но тут возникает еще один довод: считать в уме нужно уметь, если под рукой не окажется карандаша и бумаги, считать карандашом на бумаге нужно уметь, если под рукой не окажется калькулятора, и т.д. (на. манжете рубашки или на бумажной салфетке в ресторане, если вспомнить многочисленные легенды о спонтанном доказательстве различных теорем.). И далее. добывать огонь трением «нужно уметь, если под рукой не окажется» спичек или зажигалки, определять стороны света по деревьям в лесу «нужно уметь, если под
рукой не окажется» компаса и т.д. и т.п. К сожалению (или к счастью), с развитием цивилизации мы разучились считать в уме, добывать огонь трением, ориентироваться на местности и т.д. и т.п. Когда-то считали и на счетах, щелкая костяшками, на логарифмической линейке, двигая ее ползунок, на арифмометре, крутя его ручку. В античные времена был в ходу абак. Всего и не припомнишь! Школьный учитель автора не уставал повторять на занятиях по арифметике, что если мы, его ученики, не освоим счет в уме или столбиком на бумаге, то нас будут обвешивать и обсчитывать в магазинах. Сейчас этот довод не работает. В современных супермаркетах практикуют современные методы «обвешивания и обсчитывания». Устный счет тут «отдыхает».
Противники использования современных компьютерных средств решения школьных и вузовских задач по «математике, физике, химии...» опираются также на ряд других доводов, о которых они, правда, открыто не говорят.
Во-первых, многие школьные учителя и преподаватели вузов, к сожалению, просто-напросто не умеют работать с такими программами. Компьютер они освоили на уровне офисных программ (текстовый и табличный редакторы, электронная почта, работа в интернете) и азов операционной системы, но дальше идти не хотят или не могут, оправдывая это и тем, что, мол,
Валерий Федорович Очков,
д.т.н., профессор Национального исследовательского университета «Московский энергетический институт Тел.: 8 (495) 362-71-71 Эл. почта: [email protected]
Valery F Ochkov,
Dr., Professor, National research university "Moscow Power Engineering Institute " Tel.: 8 (495) 362-71-71 E-mail: [email protected]
такие программы вредны для учащихся (см. выше).
Во-вторых, внедрение этих программ в учебный процесс требует кардинального пересмотра содержания и методик преподавания, а также переписывания всех учебников и задачников или, по крайней мере, существенной их переработки. Примеры в задачниках, конечно, переписываются. Прежде были «фунты и аршины» - сейчас метры и килограммы. Раньше в задаче мы читали: «Землекоп выкопал столько-то метров канавы». Сейчас -«память компьютера имеет объем столько-то гигабайт». Но суть задач при этом, увы, не меняется.
В-третьих, вышеупомянутые компьютерные программы довольно дороги. Их не в состоянии купить наши бедные школы и вузы. Бедные в смысле «небогатые». Но и слово «беда» тут будет уместно, если вспомнить, какие беды свалились на нашу систему образования в конце прошлого века. Но эта проблема, сразу скажем, решаема. Фирмы-разработчики и их дилеры предоставляют учебным заведениям существенные скидки, а в ряде случаев передают программы «бесплатно, т.е. даром». Но с корыстной, конечно, мотивацией. Школьники и студенты, освоив бесплат-
ную программу, в конце концов, став инженерами, купят сами или попросят купить программу своего работодателя. Преподаватели должны не сетовать на дороговизну программ или на невозможность по разным причинам работы с пиратскими копиями [1], а должны искать пути решения этой проблемы, связываться, например, с разработчиками программ и их дилерами. Программа, подобная Mathcad, -SMath (www.smath.info), которую можно видеть в ранее приведенном списке («Mathematica, Maple, Matlab, Mathcad, SMath и др.») -бесплатная. Но в ней пока нет инструментов для решения приведенных в статье задач.
Рассмотрим несколько типичных задач по физике, иллюстрирующих вышесказанное.
Задача 1
Два тела бросили одновременно из одной точки: одно вертикально вверх, другое под углом в = 60° к горизонту. Начальная скорость каждого тела v0 = 25 м/с. Найти расстояние между телами через t = 1.70 с. (Эта задача взята с сайта http://exir.ru - бесплатный решатель по физике.)
Такую задачу современный школьник или студент, приняв по
Начальные условия:
Расстояние первого тела от начальной точки при t
Расстояние второго тела по горизонтали от начальной точки при I
Расстояние второго тела по вертикали от начальной точки при t
Расстояние между телами при t
Расстояние между телами при t=1.7 Б
#:=60 * г?0:=25
3
Ух(*) -t—1.0.£2
Х2 W V0 * COS (0) * *
y2(t):=V0-sin(e)-t- *T'9-t2
S(t) := ^/Х2(*)2+(У1(*)-У2(*))2 S(l.r s) = 22 m
Рис. 1. Задача о свободном полете двух тел - традиционное решение
ТО
з
1.7 s
Solve
Первое тело
Solve
Второе тело
Уі(0 з) = 0 m Уі'(Оз) = У0
y1"(t)-m1 = -m1-g Уі '■= odesolve (Vi (t), tend)
^(OsjzOm x2'(0 s) = Vq-cos^)
y2(0i) = 0m y2'(0 s) = Vo-»in(0)
x2"(t)*m2 = 0JV у2"(і)-іщ = -т2-д
У2
:= odesolve
Хз(4)
У2 (О
01807 s-~\j(0 m-x2(tend)) +(Уі (tend)-У2 (tend)) —22
Рис. 2. Задача о свободном полете двух тел - решение в среде Mathcad
умолчанию ряд допущений, скорее всего, будет решать по алгоритму, показанному на рис. 1, вспоминая или черпая из учебников и справочников по физике формулы свободного движения тел под действием силы тяжести. При этом школьник или студент часто толком не понимает, откуда взялись сами эти формулы.
Приводимые в статье расчеты привязаны к инерциальной системе отсчета, в которой справедлив первый закон Ньютона (закон инерции). Деление систем на инерци-альные и неинерциальные, как кажется автору, в первую очередь связано не с физикой той или иной задачи, а с устаревшими методами расчетов, связанными с упрощением задач так, чтобы их можно было решить теми средствами, какие были доступны во времена формирования этих задач, - без привлечения компьютера.
Формулы, показанные на рис. 1, являются решением дифференциальных уравнений движения тел: сумма сил, действующих на тело, равно произведению массы тела на его ускорение. Когда автор учился в школе, произведение массы тела на его ускорение называли на уроках физики силой инерции: тело толкают и этому действию есть противодействие в виде этой самой силы инерции. Но современные «физики» уверили автора, что такой силы (больше) нет и надо формировать задачу так, как показано выше. Ускорение же - это вторая производная пути по времени. Отсюда и возникают эти самые дифференциальные уравнения, связывающие функции и их производные. Если говорить о допущениях, то в первую очередь нужно отметить, что не учитывается сопротивление воздуха. Кроме того, по умолчанию имеется в виду, что тело бросается на земле, где ускорение свободного падения равно 9.8 м/с2. Иногда в задачниках по физике рекомендуют принимать ускорение свободного падения равным 10 м/с2. Сути задачи это не меняет, но упрощает математические выкладки при их решении. (На сайте http://twt.mpei. ас. ru/MCS/Worksheets/g-h-psi.xmcd можно рассчитать значение ускоре-
ния свободного падения в зависимости от высоты над уровнем моря и географической широты.)
Просить школьников или студентов составить такие дифференциальные уравнения еще можно, но заставить их решить эти уравнения без использования современных компьютерных математических программ нереально. Реально - это заставить школьников запомнить (вызубрить) частные решения этих уравнений, что и делается с разной степенью успешности последние несколько веков со времен, когда Ньютону «на голову упало яблоко».
Но в настоящее время стали вполне доступны школьникам, не говоря уже о студентах и преподавателях, компьютерные средства решения дифференциальных уравнений аналитическими (символьными) и численными (приближенными) методами.
На рис. 2 показано численное решение задачи о двух подброшенных телах в среде пакета Mathcad [2, 3] через использование блока Solve. В этом блоке пользователь в области Ограничения должен записать начальные условия (положение тел и их скорости в начальный момент времени) и само дифференциальное уравнение (для первого тела -см. левый блок Solve на рис. 2) или систему двух дифференциальных уравнений (для второго тела - правый блок Solve), а также (в области Решатель) встроенную в Mathcad функцию odesolve, возвращающую пользовательские функции (координаты тел во времени) с именами
Ух, х2 и у2, которые описывают движение наших двух тел под действием силы тяжести. Эту задачу через несложные преобразования можно решить в среде Mathcad и аналитически (символьно), получив формулы, записанные на рис. 1. Задача легко решается и в уме, если вспомнить следующее: если вторая производная функции (ускорение движения тела) равна константе, то сама функция - это парабола вида а + Ь • t+с • ґ2. Константы а, Ь и с этого квадратного многочлена несложно найти после анализа исходных данных при t = 0: начального положения тела и его начальной скорости. Но если вторая производная функции не является константой, а есть сложная функция от ґ (от времени), то такую задачу в уме уже не решишь. Тут нужно прибегать к помощи «калькулятора» - к современным компьютерным средствам решения дифференциальных уравнений, благо они сейчас, повторяем, вполне доступны.
Решение, зафиксированное на рис. 2, легко расширить, сняв с задачи ряд ограничений. Можно рассматривать движение тел с учетом сопротивления воздуха, приняв его, например, пропорциональным площади поперечного сечения тела, умноженной на плотность среды (воздуха) и на квадрат скорости. Можно учесть изменение по высоте значений ускорения свободного падения, плотности воздуха и др. Можно также отказаться от декартовой системы координат и перейти к полярным координатам, т.е. учитывать кривизну поверхности
Рис. З. Анимация задачи о свободном полете двух тел
земли и непараллельность векторов сил тяжести. Движение второго тела в задаче, показанной на рис. 2, мы раскладывали на две составляющие - горизонтальную (ось х) и вертикальную (ось у): получали и решали систему двух дифференциальных уравнений.
Пакет Mathcad предоставляет еще одну очень интересную возможность. Можно не только построить на графиках траектории движения тел, но и создать анимацию [4, 5]. Она отображена на рис. 3. (Сама анимация доступна на сайте http://communities.ptc.com/ videos/3221.)
Анимация упомянута здесь неслучайно. Дело в том, что школьники и студенты тянутся к компьютеру, ожидая от него в первую очередь занимательности, игры. Эта тяга часто ведет к компьютерным играм («стрелялки», «бегалки» и проч.) и к. болезненной зависимости от них. Но эту тягу можно направить по более плодотворному и здоровому руслу. Изучил школьник законы движения тел и компьютерные методы их моделирования, прочитав, например, эту статью?! Хорошо! Пусть школьник создаст анимацию движения снаряда к цели - несложную, а может быть и довольно сложную «стрелялку»!
Основное допущение при решении задач о свободном падении тела, как уже было отмечено, - это отсутствие сопротивления воздуха. В двух последующих задачах это ограничение снято.
Задача 2
Стальной шар радиусом 5 см был подброшен вертикально вверх со скоростью 100 м/с и упал на землю через 15 с. Определить значение коэффициента k в формуле сопротивления воздуха: k • £ • р • V2, где £ - площадь поперечного сечения шара, р - плотность воздуха и V - скорость шара.
На рис. 4 показано решение этой задачи в среде Mathcad. Блок Solve формирует функцию y(k) -высоту полета шара в зависимости от коэффициента k в формуле сопротивления воздуха (см. выше). Далее к этой функции добавляется еще один аргумент t и c помощью функции root (поиск нуля функции) определяется значение k, равное О.877, при котором у = О. Функция root ищет нуль функции у(15 s, k) методом половинного деления на отрезке значений k от О до 1. При k = О.877 наш шар упадет на землю (у = О) ровно через 15 с, что показано на графике. В безвоздушном пространстве (k = О) шар упал бы на землю через примерно 2О с, что следует из решения квадратного уравнения v0 • t -- V g • t2 = О, которое также вставлено в рис. 4 и решается аналитически с помощью оператора (не путать с блоком!) Solve. На графиках рис. 4 двумя горизонтальными маркерами зафиксированы и другие особые точки - максимальные высоты полета шара в вакууме (513 м) и в воздухе (275 м).
d==100 тт
Psteel
■- 7800
кд
Pair
„==1.225
кд
S-
7r-d~
= 78.54 cm2 V~—тг- — =523.599 cm
m:=V-p4teel = 4.084 kg
vO —100— vO-t——g-t = 0m-
0
200 -m
s-g
0
20.394
t:=0 s,0.1 S..30 s
Solve
y(t,0) (m) у(t,k) (m)
. 600 ' 20
400 2751 300
200 МО А "Х
3 5 8 10 13 1 > 16 20 23 25 t (s)
Рис. 4. Задача о свободном полете тела с сопротивлением воздуха
Рис. 5. Задача о парашютисте
Коэффициент k - это «головная боль» всех автомобильных и авиационных дизайнеров, которые должны проектировать автомобили и самолеты с минимальным сопротивлением воздуху. Принятая нами формула - весьма упрощенная. В реальной жизни показатель степени у скорости меняется. При ламинарном обтекании тела (без завихрений воздуха, что имеет место при невысоких скоростях, вернее, при небольших значений числа Рейнольдса, связывающего скорость тела, его размеры и вязкость среды) он близок к единице, а при турбулентном обтекании (с завихрениями) - к двойке. В формуле k • 51 • р • V2 на рис. 4 записано в виде - k • £ • р • V • VI. Это позволяет учитывать направление вектора силы сопротивления воздуха: этот вектор имеет противоположное направление вектору скорости. На рис. 5 можно видеть и другую, эквивалентную запись: k • £ • р • V2 • sign(v). где sign(v) - это встроенная в Mathcad функция, возвращающая «знак» (1, 0 или минус 1) своего аргумента.
Задача 3
Парашютист прыгает с высоты 3 км (}ц) и на высоте 500 м (Н2) раскрывает парашют. Определить, как долго он будет лететь до земли и с какой скоростью он приземлится.
Для решения этой задачи мы сформулируем такую упрощенную модель. Свободно летящий парашютист - это шар с радиусом 30 см и плотностью 950 кг/м3, а он же под парашютом - это шар радиусом 2.7 м и плотностью 1.303 кг/м3 (чуть больше плотности воздуха). При таком раскладе масса такого «парашютиста» равна 107 кг. Коэффициент в формуле сопротивления воздуха мы возьмем из предыдущей задачи: k = 0.877 (см. рис. 4). В новой задаче уже нельзя пренебрегать изменением плотности воздуха по высоте - нужно будет «порыться» в справочниках или интернете и найти формулы, связывающие температуру и давление воздуха с высотой над уровнем моря. Эти формулы можно видеть на рис.
5 над графиком изменения плотности воздуха по высоте.
Плотность и радиус нашего «парашютиста» будут также меняться по высоте. Но они будут не гладкими, а ступенчатыми функциями от высоты - см. две конструкции с if под первым графиком на рис. 5. Наше дифференциальное уравнение полета «парашютиста», связывающее силы тяжести и сопротивления воздуха, мы разобьем на два: скорость - это производная пути по времени (первое уравнение) и баланс сил, действующих на «парашютиста», где ускорение - это производная скорости по времени (второе уравнение). Функция odesolve в решении, показанном на рис. 5, возвращает нам две функции: высоту полета парашютиста у и скорость парашютиста V. Эти функции на рис. 5 отображены графически.
На графике вертикальными линиями отмечено время свободного полета парашютиста (без раскрытого парашюта) - 0.738 минут и общее время полета 5.04 минут. (Скачать файлы задачи о парашютисте можно с сайта http://communities.ptc. com/message/187928.)
14 октября австрийский парашютист Феликс Баумгартнер совершил прыжок из стратосферы с высоты 39 километров и удачно приземлился в штате Нью-Мексико. Данные о скорости его полета еще уточняются, но, скорее всего, он стал первым человеком, преодолевшим звуковой барьер без специальных устройств. Кроме того 43-летний австриец побил три мировых рекорда: самая большая высота прыжка с парашютом, самая боль-
v(t) (*pft)
Рис. 5а. Полет парашютиста из стратосферы
шая скорость падения и самый высокий пилотируемый полет на воздушном шаре. (См. http://lenta.ru/ articles/2012/10/15/skydiver.)
На рис. 5а показана попытка смоделировать этот прыжок в среде Mathcad Prime.
Задача 4
Спортсмен прыгает с десятиметровой вышки. Определить параметры прыжка: время полета до касания воды, глубину погружения в воду и др. Эта задача схожа с предыдущей. Только там резко менялись параметры «прыгуна», а в новой задаче резко меняются параметры окружающей «прыгуна» среды.
Здесь, как и в задаче о парашютисте, нужно будет принять ряд допущений. Мы примем, что прыгун в воду - это шар радиусом 30 см и плотностью 950 кг/м3. В этой задаче плотность среды так же, как и в задаче о парашютисте, будет зависеть от высоты, но не гладко, а ступенькой: при у > 0 (над водой) -это 1.225 кг/м3 (воздух), а при у < 0 (под водой) - это 1000 кг/м3 (вода). В предыдущих задачах мы игнорировали Архимедову силу. В задаче о прыгуне в воду этого делать уже нельзя. На рис. 6 дано решение задачи в среде Mathcad. Ответы отмечены маркерами на графике: через 1.42 с прыгун коснется воды, через 2.76 с он погрузится на максимальную глубину примерно 1.5 м, а через 6.62 с выньгрнет на поверхность воды и будет как маятник «колебаться на ней с затуханием». (На сайте http://communities.ptc.com/ videos/3248 можно видеть анимацию прыжка в воду и скачать расчетные файлы.)
На рис. 6 колебания прыгуна в воду на поверхности воды видны слабо. Но на вышеотмеченном сайте можно найти график изменения скорости прыгуна в воду от его высоты. Это некий аттрактор - график закручивается в некую спираль у начала координат.
Задача 5
А теперь рассмотрим еще одну популярную задачу по физике - задачу о колебании маятника и решим ее в еще одной математической программе - в среде Мар1е.
Любой мало-мальски успевающий школьник знает, что период колебания маятника зависит от его длины и значения ускорения свободного падения. Если к маятнику мысленно прикрепить карандаш, а под него подсунуть лист бумаги и равномерно двигать ее, то будет нарисована синусоида.
Многие их «многих школьников» вспомнят также формулу, по которой можно вычислить период колебания маятника, и даже упомянут, что эта формула справедлива только по отношению к математическому маятнику - к точечному маятнику, колеблющемуся в вакууме на идеально гибкой нерастя-
жимой и невесомой нити. При этом угол отклонения нити маятника от вертикали не должен превышать семи градусов. Но мало кто (не только среди школьников и студентов, но даже среди преподавателей) знает, откуда взялась эта семерка и как она определяет «математич-ность» маятника. Другие два условия (вакуум и свойства нити) касаются больше физики задачи. Многие вполне серьезно полагают, что эта семерка находится в одном ряду с другими «физическими и метафизическими» великолепными семерками: цветами радуги, нотами звукоряда, чудесами Света, днями недели, древними мудрецами, смертными грехами и т.д. и т.п. На самом деле все намного проще (или намного сложнее), но главное, интереснее.
На рис. 7 и 8 показано аналитическое (символьное) решение дифференциального уравнения движения маятника, выполненное в среде математической программы Maple.
На маятник в вакууме кроме силы тяжести (см. предыдущие задачи) действует еще одна сила -сила натяжения нити. Анализ баланса этих сил дает дифференциальное уравнение, введенное нами
h==10 т
г ==30 cm
4 а з
volume•——• 7г• г =0.113 т 3
4 ,
rnass~Pmaп’ — "п"г =107.442 kg
kg m3
area-~n-r2 =0.283 m;
p(y)==if
Solve
y>0 m,1.225 *£,1000 —j mJ m°
k ==0.877
17 з
£
О
= y(0s) = /l v(0s):0 ---
3
-mass • v' (t) = g • volume • (pman - p (y (t))) +k • p(y (t)) • area • v (t) • | v (t) | (\
odesolve
y(‘)
v(t)
> tend
•«.«*: t==0 3,----..t.
30000
end
t (.)
Г 9 9 10 U 12 13 14 IS
Рис. 6. Задача о прыгуне в воду
Рис. 7. Задача о колебании маятника - полное решение
в расчет первым оператором, показанным на рис. 7. Пакет Маріє решил это уравнение в общем виде (оператор 2 на рис. 7) и выдал частное решение (оператор 4) после ввода начальных условий (оператор 3) - начального угла отклонения нити маятника от вертикали 0^ и начальной скорости маятника, равной нулю. Решение, показанное на рис. 7, как говорится, ничего не дает ни уму, ни сердцу. Но его можно намного упростить, если. не очень сильно отклонять нить маятника от вертикали. Если угол отклонения по модулю меньше семи градусов, то синус такого угла примерно равен самому углу. В решении задачи о маятнике, показанном на рис. 8, синус угла был заменен на сам угол. Эту задачу также можно решить в уме, если вспомнить, что производная синуса равна косинусу, а производная от косинуса равна минус синусу. Отсюда и получился косинус в решении, показанном на рис. 8. Если же вторая производная функции пропорциональна не самой функции (рис. 8), а синусу функции (рис. 7), то такая задача не имеет аналитического решения, вернее, имеет «решения», показанные на рис. 7, с которыми нужно еще работать численно. Отсюда вывод - такую задачу лучше сразу решать численно.
Такое упрощение задачи позволило пакету Маріє выдать нам традиционное решение, в котором фигурирует упомянутая нами форму-
ла периода колебания маятника и его закон - синусоида.
В задаче о колебании маятника также можно использовать численные методы решения. Это позволит отойти от идеальности маятника - отклонять его на угол, больший семи градусов, учитывать сопротивление среды, эластичность нити и другие реальные факторы.
На рис. 9 указан сайт, на котором после ввода исходных данных и нажатия кнопки Recalculate (Пересчитать) выдается траектории колебания маятника в трех случаях [4]:
• нить нерастяжимая (кн = да) и сопротивления среды нет (к = 0);
• нить нерастяжимая (кн = да) и есть сопротивление среды (к > 0);
• нить растяжимая (кн < да) и есть сопротивление среды (к > 0). Здесь кн - это коэффициент в
формуле закона Гука (Hooke), гласящего, что удлинение эластичной нити пропорционально силе, к ней приложенной, а к - это коэффициент в формуле сопротивления воздуха (см. рис. 4).
Методы решения дифференциальных уравнений в сочетании со средствами анимации позволяют в среде Mathcad создавать очень интересные анимации - колебания, например, трех связанных маятников в вязкой среде (рис. 10 и сайт http:// communities.ptc.com/groups/dynamic-models-in-mathcad). В среде Mathcad такую задачу решить довольно просто - нужно только правильно составить дифференциальные уравнения.
Маятники, показанные на рис. 10, оставляют за собой цветные следы, формирующие замысловатые узоры. На синий (последний в цепочке) грузик действуют сила тяжести и сила натяжения синей нити, а точка крепления (красный грузик) синего (среднего) маятника подвижна; на красный маятник кроме силы тяжести действу-
1. Ввод дифференциального уравнения
> ode ■= =-у ot(f)
ode := —— a(t) ■ Аґ
2. Решение в общем виде
> ско1\'е{ос1е)
а(*) = ^1111 | + ^С05|
3. Ввод начальных условий
> 1С5 >- а(0) = а)(0) = О
/с* := а( 0) = Б( а) (0) = О
4. Решение с учетом начальных условий
> ско!\е( {ос!е, гсз})
/7 <
a(f) = ос cos
/г
Рис. 8. Задача о колебании маятника - упрощенное решение
Рис. 9. Задача о колебании маятника на Mathcad Server^
■ Play Animation
vblack(*) = 0007 m/s vred(t) = 0.038 m/s Vb|Ug(t) = 0.094 m/s 1 J / / / / / / / 0black = -3-Q4°
9red = -0.81° 0blue = 1-O4°
_ J t = 9.80s
■ ш і и— і
Рис. 10. Анимация задачи о колебании маятника
ют силы натяжения синей и красной нитей и т.д. Только точка крепления черного (верхнего) маятника неподвижна.
В [5] дано описание решения с помощью современных компьютерных средств еще одной интересной задачи - задачи о движении тел под действием сил, определяемым законом всемирного тяготения. На сайте http://communities. ptc.com/groups/dynamic-models-in-mathcad можно увидеть анимации
об описании динамических моделей, рассмотренных выше, и многих других. Модели движущихся физических конструкций, без учета сил, действующих на них, собраны на сайте http://communities. ptc.com/groups/kinematic-models-in-mathcad. При их создании решались не дифференциальные уравнения и системы, а алгебраические уравнения и системы. Описание этих решений дано на сайте http:// twt.mpei.ac.ru/ochkov/Mathcad-15/ kinematic.html.
Современные информационные технологии позволяют преобразовать не только решение задач на семинарских занятиях в вузе. Кардинально можно изменить и чтение лекций. На лекции теперь представляется возможным, во-первых, показать реальный физический опыт (колебание связанных маятников, например), потом на большом экране (замедленно или, наоборот, ускоренно) - соответствующий фильм, так, чтобы детали движения были всем студентам ясно видны, а затем - анимацию этого движения, рассказав, как она создавалась, какие в нее заложены физические законы, допущения и ограничения.
Выводы
1. Современные компьютерные средства решения задач позволяют по-новому поставить преподавание физики в школе и вузе, учитывая тягу школьников и студентов к компьютерам.
2. С задач по физике можно снять многие искусственные ограничения и допущения, приблизив их к реальности.
3. Современные информационные технологии позволяют преобразовать традиционные решения задач (поиск формул и расчет по ним) в написание некоего реферата, где численные и/или аналитические решения будут дополнены рассуждениями об истории задачи, о развитии методов ее решения, о допустимости тех или иных ограничений и т.д. в стиле, например, данной статьи.
Продолжение данной статьи можно найти на сайте http:// twt.mpei.ac.ru/ochkov/Mathcad-15/ Physic.pdf. Автор благодарит Алана Стеванса (Alan Stevens - http:// communities.ptc.com/docs/DOC-3019) за подготовку решений задач на форуме PlanetPTC (http:// c o mmunities .ptc.co m/community/ mathcad).
Литература
1. Очков В.Ф. Открытое письмо пирата в адрес компьютерных изданий // Глава из книги «Mathcad PLUS 6.0 для студентов и инженеров». - М.: КомпьютерПресс, 1996. - http://twt.mpei.ac.ru/ochkov/mc8Pro.book/6_text. htm#_Toc525449695.
2. Очков В. Mathcad 14 для студентов и инженеров: русская версия. - СПб.: BHV, 2009. - http://twt.mpei.ac.ru/ ochkov/Mathcad_14/RusIndex.html
3. Гурский Д., Турбина Е. Mathcad для студентов и школьников. Популярный самоучитель. - СПб.: Питер, 2005.
4. Очков В. MCS (Mathcad Calculation Server) на занятиях по математике, физике, информатике. // Компьютерные учебные программы и инновации. - 2008. - № 3. - С. 187-194. - http://twt.mpei.ac.ru/ochkov/Pendulum
5. Очков В. MAS (Mathcad Application Server) на занятиях по математике, физике, информатике. // Компьютерные учебные программы и инновации. - 2006. - № 2. - http://twt.mpei.ac.ru/ochkov/Mathcad_12/Planet
