Решение инженерно-технических задач в среде Mathcad

Математический пакет Mathcad [23, 42, 44] позволяет специалистам в конкретной научно-технической области быстро освоить работу на компьютере и реализовать на них математические модели, не вдаваясь в тонкости программирования. Mathcad — это зарегистрированная торговая марка фирмы MathSoft, Inc. (http://www.mathsoft.com), представляющая собой сокращение английской фразы: Mathematical Computer Aided Design — математическое проектирование с помощью компьютера.

Ниже приведены три примера решения типовых задач теплоэнергетики [68], которые дают представления об особенностях работы в среде Mathcad.

Расчет термодинамических свойств воды и водяного пара.

Ниже представлено формирование функции пользователя vw (p), предназначенной для определения удельного объема кипящей воды на линии насыщения (vw) от давления (p), и функции hss (T, p), предназначенной для определения удельной энтальпии перегретого пара (hss) от температуры (T) и давления (p).

Расчет удельного объема кипящей воды на линии насыщения сводится к интерполяции (здесь задействована встроенная Mathcad-функция interp) кубическими сплайнами (cspline) табличных данных, хранящихся в двух векторах: P — табличные значения давления и V — табличные значения удельного объема. Элементы векторов — величины размерные: давление измеряется в мегапаскалях (10⁶Pa), а удельный объем — в литрах на килограмм (m³/1000kg). В векторах P и V (они транспонируются из векторов-строк в вектора-столбцы) показаны только крайние элементы (далее будут показаны экраны дисплея, отображающие работу в среде Mathcad):

Интерполяция сплайном для определения удельного объема кипящей воды ведется в диапазоне давлений от 611 Па до 22,088 МПа. Функция error отслеживает значение аргумента p, прерывает выполнение вычислительной процедуры и выдает пользовательское сообщение об ошибке («Давление низкое» или «Давление высокое»), если значение аргумента не попадает в оговоренный диапазон (эту работу выполняет оператор if).

Созданную функцию с одним аргументом можно вызвать в традиционном виде (функция), а также в виде префиксного и постфиксного операторов.

При этом в расчетах возможны различные размерности давления, объема и массы (здесь kgf — килограмм силы, galUK — британский галлон, lb — фунт и lbf — фунт силы и ft — фут). Далее представлено решение обратной задачи через встроенную Mathcad-функцию root (поиск корня уравнения) — определение давления при заданном значении удельного объема кипящей воды.

Расчет удельной энтальпии перегретого пара ведется также через кубический сплайн, но уже двухмерной интерполяцией. Табличные значения энтропии (килоджоули на килограмм) заносятся в матрицу M, «боковик» которой (без первого элемента) значения температуры (градусы Цельсия), а «шапка» — давления (мегапаскали). Далее эти значения с помощью встроенной функции submatrix разносятся по матрицам TP (у нее два столбца и 18 строк) и M (квадратная матрица с порядком 18). Так достигается компромисс между требованиями наглядного представления табличных данных и спецификой двумерной сплан-интерполяции в среде Mathcad.

К пользовательским сообщениям об ошибках ввода неверных значений аргументов (температура и давление пара) добавлено новое — «Ошибка: измените p и/или T». Ошибка возникает в том случае, если рассчитанное значение энтальпии окажется меньше наименьшего табличного значения (1930,8 кДж/кг), что может случиться на границах матрицы H, где вместо отсутствующих термодинамических табличных данных записано: -1.

Пользовательскую функцию h_ss(T, p) с двумя аргументами можно вызывать в традиционной форме (функция), а также инфиксным и древовидным оператором.

Аналогичным образом построены и другие функции термодинамических свойств воды и водяного пара: удельная энтальпия кипящей воды на линии насыщения h_w(p), удельная энтальпия сухого пара на линии насыщения h_s(p), удельная энтропия кипящей воды на линии насыщения s_w(p), удельная энтропия сухого пара на линии насыщенияs_s(p), удельная энтропия перегретого пара s_ss(T, p) и др. Соответствующий самораскрывающийся архивный файл с именем WaterSteamPro. exe хранится на FTP-cервере по Internet-адресу: http://twt.mpei.ac.ru/orlov/watersteampro/. Этот файл можно по сети Internet скопировать на компьютер конкретного пользователя, разархивировать файл и задействовать его в термодинамических расчетах, что и будет сделано при решении второй задачи.

Особенности Mathcad, проиллюстрированные примерами:

1. Интерфейс пользователя в среде Mathcad опирается на операторы ввода и операторы вывода значений переменных и выражений.

В среде Mathcad три оператора ввода:

A ¦ - ввод значения локальной переменной, область видимости которой ограничена местом ввода переменной и концом вычислительной процедуры, отмеченной вертикальной чертой, фиксирующей начало и конец программного блока (см. выше программу формирования функции v_w(p));

A := ¦ - ввод значения полуглобальной переменной, область видимости которой — это часть Mathcad-документа, расположенная правее и ниже места ввода переменной;

A ¦ - ввод значения глобальной переменной, область видимости — весь Mathcad-документ.

Кроме того, в Mathcad-документе возможны ссылки на другой Mathcad-документ (на соответствующий Mathcad-файл на диске — файл с расширением mcd), переменные и пользовательские функции которого, становятся видимыми в текущем Mathcad-документе.

В среде Mathcad два оператора вывода:

A = ¦ ¦ - вывод численного значения переменной или выражения: на месте первого операнда появляется числовая константа, а на месте второго соответствующая размерность (см. табл. 5.17), если константа размерная; пользователь вправе менять размерность по умолчанию на любую другую, определенную к данному моменту;

A ¦ - вывод символьного значения переменной или выражения.

Переменные в среде Mathcad могут хранить скалярное, векторное и матричное значения.

Скаляром, а также элементом вектора (V) и матрицы (M) может быть число (целое, вещественное или комплексное — см. элементы вышеприведенного вектора V), текст, а также новый вектор или матрицу (механизм составных массивов). Булевы значения «True» (Истина) и «False» (Ложь) в среде Mathcad имеют значения 1 и 0, соответственно. Рассчитанные значения могут быть выведены в десятеричном (умолчание), двоичном, восьмеричном или шестнадцатеричном формате.

Переменные в среде Mathcad подразделяются на пользовательские и встроенные (предопределенные табл. 5.16).

интерфейс мathcad сплайн Таблица 5.16. Встроенные переменные Mathcad.

Имена переменных в среде Mathcad могут содержать пробелы, буквы латиницы, кириллицы и (чего нет в традиционных языках программирования) греческого алфавита. Кроме того, в именах переменных допустимы нижние индексы и штрихи. Все это позволяет делать Mathcad-документ максимально «читабельным» — давать переменным и функциям те имена, которые закрепились за ними в процессе формирования той или иной научной дисциплины задолго до появления компьютеров.

2. В работе с Mathcad есть особенность, позволяющая называть Mathcad не просто математическим, а физико-математическим пакетом. Часто решая физическую задачу, пользователь делает ошибки не в формулах и не в счете, а в размерностях физических величин. Пакет Mathcad поддерживает основные системы физических величин (СИ, килограмм-метр-секунда, грамм-сантиметр-секунда и британскую систему единиц) и ведет контроль за соответствием размерностей (табл. 5.17).

Для присваивания величине размерности за числовой константой ставится знак умножения (но его можно и не ставить) и вводится название соответствующей размерности. — L := 5 m (или L := 5 m, что более соответствует общепринятой норме записи размерностей).

Таблица 5.3.2. Список единиц измерения физических величин, встроенных в Mathcad.

Для ввода размерностей предусмотрено диалоговое окно «Ввод размерности» («Insert Unit»), определяющее физическую величину (Dimension), размерность (Unit) и систему единиц (System).

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

функции Бесселя (Bessel);

функции комплексных чисел (Complex Numbers);

функции решения дифференциальных уравнений и систем (задача Коши, краевая задача, уравнения в частных производных — Differential Equation Solving);

функции типа выражения (Expression Type);

функции работы с файлами (File Access);

функции преобразований Фурье (Fourier Transform);

гиперболические функции (Hyperbolic);

функции обработки образов (Image Processing);

функции интерполяции и экстраполяции (Interpolation and Prediction);

логарифмические и экспоненциальные функции (Log and Exponential);

функции теории чисел и комбинаторики (Numbers Theory/Combinatorics);

функции ступенек и условия (Piecewise Continuous);

функции плотности вероятности (Probably Density);

функции распределения вероятности (Probably Distribution);

функции случайных чисел (Random Numbers);

функции регрессии и сглаживания (Regression and Smoothing);

функции обработки сигналов (Signal Processing);

функции решения алгебраических уравнений и систем, а также решения оптимизационных задач (Solving);

функции сортировки (Sorting);

специальные функции (Special);

статистические функции (Statistics);

текстовые функции (String);

тригонометрические функции (Trigonometric);

функции округления и работы с частью числа (Truncation and Round-Off);

функции работы с векторами и матрицами (Vector and Matrix);

функции волнового преобразования (Wavelet Transform).

Для функции предусмотрено диалоговое окно «Ввод функции» («Insert Function»), показывающее группы функций (Function Category), имя функции (Function Name), ее синтаксис и краткую аннотацию.

В систему Mathcad, кроме того, интегрированы средства символьной математики, что позволяет решать поставленные задачи (этап задачи) не только численно, но и аналитически см. третью задачу.

Кроме математических функций в Mathcad-документе можно работать с математическими операторами, которые вводятся через нажатие соответствующих кнопок панелей инструментов.

Одна из причин популярности Mathcad заключается в том, что пользователь вправе вставлять в документы либо функцию, либо оператор в зависимости от того, к чему он привык, изучая математику в школе или в институте. Благодаря этому Mathcad-документ максимально похож на лист с математическими выкладками, написанными от руки или созданными в среде какого-либо текстового процессора (MS Word, ScientificWord, ChiWriter и др.).

• 4. Математические выражения в среде Mathcad записываются в их общепринятой нотации: числитель находится сверху, а знаменатель — внизу, в интеграле пределы интегрирования также расположены на своих привычных местах. Казалось бы, это все мелочи, никак не влияющие на вычислительный процесс. Но!.. Программа должна быть понятной не только для компьютера, но и для человека. Пользователь, анализируя Mathcad-документ на экране дисплея или на бумаге принтера, видит, что данная величина записана в числителе и ее рост приводит к возрастанию всего выражения. А это очень важно при анализе математических моделей, форма и содержание которых едины.
• 5. В среде Mathcad процесс создания программы идет параллельно с ее отладкой и оптимизацией. Пользователь, введя в Mathcad-документ новое выражение, может не только сразу подсчитать, чему оно равно при определенных значениях переменных и в выбранной системе размерностей, но и построить график или поверхность, беглый взгляд на которые может безошибочно показать, где кроется ошибка, если она была допущена при вводе формул или при создании самой математической модели. Отладочные фрагменты можно оставить в готовом документе для того, чтобы, например, еще раз убедить воображаемого или реального оппонента в правильности модели. Система Mathcad оборудована средствами анимации, что позволяет реализовать созданные модели не только в статике (числа, таблицы, графики), но и в динамике (анимационные клипы).
• 6. Пакет Mathcad дополнен справочником по основным математическим и физико-химическим формулам и константам, которые можно автоматически переносить в документ без опасения внести в них искажения, нередкие при ручной работе. К пакету Mathcad можно приобрести те или иные электронные учебники по различным дисциплинам: решение обыкновенных дифференциальных уравнений, статистика, термодинамика, теория управления, сопротивление материалов и т. д. Прежде чем решать возникшую проблему, пользователь может изучить электронный учебник и перенести из него в свой документ нужные фрагменты, отдельные формулы и константы.
• 7. Не выходя из среды Mathcad, возможно открывать новые документы на других серверах и пользоваться теми преимуществами информационных технологий, предоставляемых Internet. Ниже приведено диалоговое окно, открываемое из среды Mathcad, с помощью которого пользователи Mathcad могут обмениваться информацией и вести совместные проекты по таким разделам: математика и статистика, астрономия и навигация, электроэнергетика, физика, химия и химическая промышленность и т. д.

Пользователи пакета Mathcad (или пакета Mathcad Explorer свободно распространяемой в Internet укороченной версии Mathcad) могут послать на сервер свое сообщение (кнопка Post…) или вступить в дискуссию (Post Reply…) по интересующему его вопросу.