Оглавление

1 АРХИТЕКТУРА И ПОСЛЕДОВАТЕЛЬНОСТЬ РАБОТЫ С ПРОГРАММНОЙ СИСТЕМОЙ
 1.1 Архитектура программной системы
 1.2 Установка программной системы
 1.3 Организация расчетов в системе
 1.4 Последовательность работы с системой
  1.4.1 Создание модели и программы ее решения (шаг МОДЕЛЬ)
  1.4.2 Интегрирование уравнений движения (шаг РАСЧЕТ)
  1.4.3 Обработка и вывод результатов решения
 1.5 Использование макро в оболочке
2 ФОРМАТ ОПИСАНИЯ РАСЧЕТНОЙ СХЕМЫ
 2.1 Общие замечания
 2.2 Синтаксис языка описания моделей
  2.2.1 Допустимые знаки
  2.2.2 Определения
  2.2.3 Структура файла исходных данных
  2.2.4 Структура блока данных
  2.2.5 Типы строк в файле исходных данных
 2.3 Формат описания структуры расчетной схемы модели
  2.3.1 Блок описания подсистем (тел) модели
  2.3.2 Блок описания характеристик соединительных элементов
  2.3.3 Блок описания соединительных элементов (с. э.)
  2.3.4 Блок описания специальных сил (с.с.)
  2.3.5 Блок задания узлов действия специальных сил
 2.4 Формат описания параметров расчетной схемы модели
  2.4.1 Блок описания инерционных параметров подсистемы твердое тело
  2.4.2 Блок геометрических параметров точек крепления подсистемы твердое тело
  2.4.3 Блок описания параметров функций сил в механической характеристике
  2.4.4 Блок описания параметров функции специальных сил
  2.4.5 Блок описания коэффициентов демпфирования по формам подсистемы упругое тело
  2.4.6 Блок задания расстояния между первым и последующим мостом
  2.4.7 Блок описания параметров колеи
 2.5 Размерность вводимых данных
 2.6 Препроцессор ввода описания модели
  2.6.1 Глобальные переменные
  2.6.2 Вставки
  2.6.3 Условный оператор
  2.6.4 Макроопределения
3 УПРАВЛЯЮЩАЯ ИНФОРМАЦИЯ ДЛЯ ПРОГРАММЫ ИНТЕГРИРОВАНИЯ УРАВНЕНИЙ ДВИЖЕНИЯ
 3.1 Общие положения
 3.2 Управляющие переменные и формат их задания
  3.2.1 Задание типа процедуры интегрирования
 3.3 Управляющая информация для расчета собственных частот и форм
 3.4 Управляющая информация для расчета моделей с кинематическими парами
 3.5 Управляющая информация для сборки моделей с кинематическими парами
 3.6 Задание начальных условий
 3.7 Управляющая информация для определения статического положения системы
4 ОБРАБОТКА И ВЫВОД РЕЗУЛЬТАТОВ РЕШЕНИЯ
 4.1 Задание контрольной точки
 4.2 Описание классов вывода.
 4.3 Описание кодов вывода
 4.4 Визуализация графических зависимостей
  4.4.1 Вывод графиков на экран
  4.4.2 Формат представления данных для построения графика
 4.5 Анимация результатов расчетов
  4.5.1 Входные данные программы.
  4.5.2 Управление картинкой при визуализации.
  4.5.3 Формат представления координат тел.
  4.5.4 Формат представления геометрии тел.
  4.5.5 Формат представления сил, действующих на точки тел
  4.5.6 Другая дополнительная информация для аниматора
  4.5.7 Режим работы ”кино”
  4.5.8 Изображение соединительных элементов при анимации
5 ВСПОМОГАТЕЛЬНЫЕ ПОДПРОГРАММЫ - ПОДГОТОВКА МИКРОПРОФИЛЯ И Т.Д.
 5.1 Подготовка микропрофиля
 5.2 Нормирование собственных форм
 5.3 Связь с конечно-элементными комплексами и передача форм колебаний упругих тел
  5.3.1 Взаимодействие с комплексом NASTRAN
  5.3.2 Взаимодействие с комплексом BASYS
  5.3.3 Формат файла с собственными частотами и формами
 5.4 Подготовка данных по рельсовому пути
6 ТИПЫ ФУНКЦИОНАЛЬНЫХ ЗАВИСИМОСТЕЙ СИЛ В СОЕДИ НИТЕЛЬНЫХ ЭЛЕМЕНТАХ И СИЛ НАГРУЗОК
 6.1 Структура сил в соединительных элементах
 6.2 Типы функций, описывающих упругое слагаемое
 6.3 Типы функций, описывающих демпфирующее слагаемое
 6.4 Структура сил нагрузок (специальных сил)
 6.5 Типы функций спец. сил
 6.6 Свободные типы функций
  6.6.1 Свободные типы сил в соединительных элементах
  6.6.2 Свободные типы сил нагрузок
  6.6.3 Включение подпрограмм свободных типов в систему
7 ОПИСАНИЕ МОДЕЛЕЙ СИСТЕМ ТЕЛ С КИНЕМАТИЧЕСКИМИ ПАРАМИ
 7.1 Задание кинематических пар
  7.1.1 Одномерный случай
  7.1.2 Плоский случай
  7.1.3 Пространственный случай
8 ОПИСАНИЕ И РЕДАКТИРОВАНИЕ ГЕОМЕТРИИ ТЕЛ И МОДЕЛИ
 8.1 Cистемы координат
  8.1.1 Глобальная система координат
  8.1.2 Локальная система координат
 8.2 Занесение тел в библиотеку
  8.2.1 Формат данных для занесения тела в библиотеку
 8.3 Редактирование геометрии
9 БОЛЬШИЕ ДВИЖЕНИЯ СИСТЕМ ТЕЛ, ДИНАМИКА МЕХАНИЗМОВ
 9.1 Общие положения
 9.2 Задание кинематических пар
 9.3 Дополнительная информация для расчета больших движений
  9.3.1 Дополнительная информация для шага МОДЕЛЬ
  9.3.2 Дополнительная информация для шага РАСЧЕТ
 9.4 Задание степеней свободы
 9.5 Задание начальных условий
 9.6 Задание нагрузок
 9.7 Податливые связи
  9.7.1 Пружина в шарнирах
  9.7.2 Многомерный податливый элемент
  9.7.3 Податливая кинематическая пара.
  9.7.4 Задание шины
 9.8 Большие движения упругих тел
 9.9 Вывод результатов
  9.9.1 Вывод кинематических параметров тел
  9.9.2 Вывод реакций в кинематических парах
  9.9.3 Вывод смещений в кинематических парах
 9.10 Интегрирование уравнений, устойчивость решения
10 МОДЕЛИРОВАНИЕ СИСТЕМ С КОНТУРАМИ УПРАВЛЕНИЯ
 10.1 Описание связей в контурах управления
  10.1.1 Описание входа от соединительных элементов
  10.1.2 Описание входа от специальных сил
 10.2 Типы элементов контуров управления
  10.2.1 Двухполюсники
  10.2.2 Многополюсники
  10.2.3 Свободно программируемые многополюсники
  10.2.4 Многополюсники для задания траекторий движения точек тел
 10.3 Вывод параметров элементов контуров управления
11 РАБОТА ФРУНДА СОВМЕСТНО С ПАКЕТОМ ОПТИМИЗАЦИИ И ИДЕНТИФИКАЦИИ MOVI
 11.1 Особенности работы пакета MOVI.
 11.2 Входные данные для оптимизации.
  11.2.1 Описание варьируемых параметров
  11.2.2 Файл описания критериев оптимизации OPKRI
 11.3 Контрольные точки для вычисления критериев.
12 МОДЕЛИРОВАНИЕ СПЕЦИФИЧЕСКИХ ВЗАИМОДЕЙСТВИЙ - ШИНА АВТОМОБИЛЯ, КАЧЕНИЕ ЖЕЛЕЗНОДОРОЖНОГО КОЛЕСА И Т.Д.
 12.1 Модель качения эластичного колеса
 12.2 Модель рулевых воздействий водителя
 12.3 Модель распределения тяги по колесам с учетом дифференциалов
 12.4 Модель качения железнодорожного колеса
 12.5 Описание движения железнодорожного экипажа в кривых
 12.6 Модель виброизолирующей гидроопоры

Глава 1
АРХИТЕКТУРА И ПОСЛЕДОВАТЕЛЬНОСТЬ РАБОТЫ С ПРОГРАММНОЙ СИСТЕМОЙ

1.1 Архитектура программной системы

Система моделирования ФРУНД состоит из нескольких программных модулей, выполняющих различные функции при расчете динамики исследуемого объекта.

Список имен выполняемых модулей:

1.
avmodel* - ввод тел в библиотеку,
2.
crmodel* - ввод структуры и параметров расчетной схемы,
3.
vizual* - визуализация и редактирование геометрии и расчетной схемы,
4.
umnform - нормирование коэффициентов форм упругих тел,
5.
obpron - сглаживание случайного микропрофиля,
6.
model32 - генерация уравнений движения модели,
7.
sbor32 - сборка программы интегрирования уравнений движения,
8.
rashet32 - интегрирование уравнений движения,
9.
anaz - обработка результатов расчета (интегрирования),
10.
viv - вывод результатов обработки,
11.
benafin* - визуализация графиков,
12.
multic* - анимация движений,
13.
frund - сервисная оболочка для работы с системой,
14.
hf - справочник по системе.

Програмы системы работают в операционной системе MS DOS, версии 5.0 и выше.

В систему входит также библиотека объектных модулей rashet32.lib. Кроме того, система для своей работы использует компилятор с языка ФОРТРАН и сборщик программ. В частности используется Microsoft FORTRAN PowerStation Version 1.0. Для моделирования колебаний транспортных средств могут поставляться файлы с значениями ординат микропрофиля некоторых дорожных покрытий, полученные прямыми замерами.

1.2 Установка программной системы

Все системные файлы ФРУНДа должны лежать в директории с именем UTILF. Для установки системных файлов, находясь в этой директории нужно набрать строку arj x -v utilf.arj Директория UTILF должна быть включена в PATH. В autoexec.bat требуется добавить строку «SET FRUND = dirname», где dirname - полное имя директории UTILF.

Перед первым запуском программы следует также поставить полное имя библиотеки rashet32.lib, которая находится в директории UTILF, в файле rashet32.bld.

Оболочка системы - bat файл с именем frund использует в работе редактор текстов, в качестве которого может применяться любой из имеющихся в распоряжении пользователя редакторов. Для настройки на требуемый редактор нужно внести имя его вызываемого модуля в файл editor.dat директории UTILF. По умолчанию используется редактор Multi-Edit - в файле editor.dat занесено "me".

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

1.
*.mhl - modhelp;
2.
*.upr - uprf;
3.
*.ico - initcond.inp;
4.
*.int - epsilon.dat;
5.
*.rhl - rashelp.

Программы, помеченные в п.п. 1.1 звездочкой, могут работать в двух режимах ввода данных - запрос входного файла с терминала и считывание имени входного файла из вспомогательного файла. Режим работы определяется управляющей переменной хранящейся в файле с именем *.cnt , где * - имя программы. Значение управляющей переменной 0 - ввод по запросу , 1 -считывание из файла *.tmp, где * имя программы.

1.3 Организация расчетов в системе

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

Ниже рассматривается работа каждого программного модуля совместно с пунктом меню оболочки, который ему соответствует. Работу с моделями рекомендуется проводить в индивидуальных директориях для каждой модели. Для вызова оболочки набирается команда frund, и далее работа ведется по пунктам меню.

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

1.4 Последовательность работы с системой

Программная система функционально разбита на две части (два шага) - шаг МОДЕЛЬ и шаг РАСЧЕТ. На шаге МОДЕЛЬ создается численная модель заданной структуры. На шаге РАСЧЕТ производится интегрирование уравнений движения, обработка и вывод результатов решения. На шаге РАСЧЕТ может проводится многократная модификация численных значений параметров модели и повторные расчеты без изменения структуры модели.

1.4.1 Создание модели и программы ее решения (шаг МОДЕЛЬ)

Описание расчетной схемы

Описание расчетной схемы производится на специальном языке - см. главу 2. Расчетная схема задается двумя частями - структурой и параметрами.Описание структуры расчетной схемы заносится в произвольный файл с расширением mdl, а описание параметров в файл с расширением ras. Эти файлы могут непосредственно набираться пользователем с помощью редактора текстов, или формироваться в из программы vizual - пункт меню оболочки - ’Редактирование геометрии и структуры’. При описании структуры расчетной схемы модели используются типы функций, определяющих силовую характеристику связей (соединительных элементов) и силы нагрузки (спец. силы) приведенные в главе 6.

Для описания тел используются существующие в библиотеке типы подсистем (тел) или, при необходимости, создаются новые - см. главу 7

Для упрощения описания расчетной схемы можно использовать существующие и создавать новые макроопределения - п.п. 2.6.

Ввод структуры расчетной схемы

Ввод расчетной схемы проводится программой crmodel - пункт меню оболочки - ‘Ввод структуры расчетной схемы’. При инициировании этого пункта меню требуется указать файл с параметрами модели - файл с расширением mdl. В файл crmodel.ls дублируется содержимое входного файла с диагностикой ошибок. При работе из оболочки перед запуском программы ввода crmodel входной файл обрабатывается препроцессором, раскрывающим макроопределения. Результат работы препроцессора на шаге МОДЕЛЬ сбрасывается в файл с постоянным именем modelp.mdl. Если в расчетную схему входят упругие подсистемы, то после ввода расчетной схемы нужно обязательно нормировать их собственные формы в соответствии с системами единиц, которые использовались при расчете конечно-элементной модели и которые будут использоваться при расчете данной модели. Такая нормировка производится программой umnform - пункт меню оболочки - ‘Нормировка коэффициентов форм упругих тел’, для каждой подсистемы отдельно - см. главу 5.

После ввода структуры расчетной схемы образуются или обновляются следующие файлы:

1.
********.gmr - файл с геометрией всей модели, где ******** - код задачи;
2.
******NP.gmr - файлы с геометрией тел в формате аниматора, где ****** - первые шесть символов кода задачи, NP - номер подсистемы. Указанные файлы в дальнейшем используются аниматором;
3.
******NP.suz - файлы соответствия номеров узлов в конечно-элементной модели и номеров в изображении тела (для упругих тел), где ****** - первые шесть символов кода задачи, NP - номер подсистемы;
4.
********.anc - управляющий файл для аниматора, где ******** - код задачи;
5.
default.upr - файл с управляющей информацией для решения и обработки, в котором сгенерированы контрольные точки по коду 19 для всех тел (вывод данных для анимации движения). Если для какого либо тела не задана геометрия, то в этом файле тело будет пропущено;
6.
default.ico - файл с начальными условиями для интегрирования, в котором за начальное положения тел приняты их положения, заданные геометрией;
7.
defaulm.upr - файл описания контрольных точек для обработки результатов решения, описанные в макроопределениях - п.п. 2.6;

default.mhl - файл с дополнительными данными о расчетной схеме, необходимый, например, при расчете больших движений - п.п. 9.3.

После ввода структуры расчетной схемы можно посмотреть геометрию всей модели программой vizual. Геометрия модели строится объединением изображений отдельных тел, которые хранятся в библиотеке тел. Файл с геометрией всей модели - ********.gmr.

Формирование уравнений движения и генерация программы их интегрирования

Формирование уравнений и генерация программы их решения производится программой model32 - пункт меню оболочки - ‘Генерация модели - шаг МОДЕЛЬ’. При запуске этого пункта запрашивается файл с дополнительными данными для шага МОДЕЛЬ, как правило это файл default.mhl. После отработки этого пункта образуются следующие файлы:

1.
idant.inf - в файл выводятся для справок список параметров, полученных уравнений движения, и символический вид самих уравнений;
2.
fpravf, fpravff - подпрограммы вычисления правых частей уравнений движения, в первом файле содержится заголовок подпрограммы, во втором - исполнимые операторы;
3.
fpravfl, fpravffl - подпрограммы вычисления линейных составляющих правых частей уравнений движения, в первом файле содержится заголовок подпрограммы, во втором - исполнимые операторы;
4.
fas1, fas2.. - файлы подпрограмм вычисления коэффициентов уравнений связей;
5.
fsila.for - файл подпрограмм вычисления сил в соединительных элементах и от нагрузок;
6.
defaults.upr - файл описания контрольных точек для вывода графиков сил в соединительных элементах. Все контрольные точки в этом файле закомментированы;
7.
animasil.upr - файл описания контрольных точек для показа сил при анимации. Все контрольные точки в этом файле закомментированы;
8.
pointsil - файл указания точек приложения сил при анимации. Все описания точек в этом файле закомментированы. Порядок следования точек в этом файле соответствует порядку контрольных точек в файле animasil.upr;
9.
********.sel - файл с описанием изображений соединительных элементов при анимации. Все описания в этом файле закомментированы.

Для программы model32 существует несколько файлов с ключами, задающими особенности генерирования уравнений движения. К таким файлам относятся - weight и signv. Учет изменения направления силы веса. При использовании уравнений малых перемещений не учитывается эффект ’маятника’, определяемый проекцией силы веса на плоскость XY (горизонтальную плоскость). В уравнениях малых движений это эквивалентно добавлению в уравнения движения тела по X и Y проекций ускорения свободного падения на эти оси от углов поворота тела вокруг этих же осей. Для включения во ФРУНДе этой возможности перед генерацией модели в файл WEIGHT нужно ввести 1. При отсутствии этого файла, эффект маятника не учитывается.

Указание направления перемещения профиля кинематического возмущения. По умолчанию предполагается, что микропрофиль под колесом двигается в отрицательном направлении оси X. Для задания движения в положительном направлении, в файл с именем signv заносится значение 2.

В случае моделирования больших движений тел перед запуском программы model32 готовится соответствующий файл modhelp - см. п.п. 9.3.1.

Сборка программы интегрирования уравнений движения

Сборка программы интегрирования производится пакетным файлом sbor32.bat - пункт меню оболочки - ‘Сборка программы решения уравнений модели’. Процедура производит компиляцию программ вычисления правых частей уравнений и сборку их с постоянными модулями процедур интегрирования (библиотека rashet32.lib). После окончания работы пакетного файла образуется программа (загрузочный модуль) с именем rashet32.exe, которая в дальнейшем используется для конкретных расчетов по данной расчетной схеме. После окончания работы процедуры sbor32 все объектные модули программы вычисления правых частей для данной модели автоматически уничтожаются.

1.4.2 Интегрирование уравнений движения (шаг РАСЧЕТ)

Описание параметров уравнений движения

Описание параметров уравнений движения производится на специальном языке - см. главу 2.

Описание параметров модели заносится в произвольный файл с расширением ras. Этот файл может непосредственно набираться пользователем с помощью редактора текстов, или формироваться в из программы vizual - пункт меню оболочки - .’Редактирование геометрии и структуры’. В случае необходимости готовятся также файлы со случайным микропрофилем - см. п.п. 5.1.

Ввод параметров математической модели.

Ввод параметров математической модели производится программой crmodel аналогично шагу МОДЕЛЬ. В оболочке системы соответствующий пункт меню - ‘Ввод параметров расчетной схемы’. При инициировании этого пункта меню требуется указать файл с параметрами модели - файл с расширением ras. При работе из оболочки перед запуском программы ввода crmodel входной файл обрабатывается препроцессором, раскрывающим макроопределения. Результат работы препроцессора на шаге РАСЧЕТ сбрасывается в файл с постоянным именем rashetp.ras. После ввода структуры расчетной схемы образуются или обновляются следующие файлы:

1.
default.rhl - файл с дополнительными данными для шага РАСЧЕТ, необходимый, например, при расчете больших движений - п.п. 9.3.2.
2.
crmodel.ls - протокол ввода параметров с сообщениями об ошибках.

Все пропущенные в файле описания параметров модели инерционные и геометрические параметры берутся из библиотеки. Другими словами блоки BL31 и BL32 могут быть опущены в описании параметров модели. Если геометрия тела была откорректирована и тело заново занесено в библиотеку со старым номером типа, то для учета этих изменений перед расчетом параметры модели вводятся заново.

Интегрирование уравнений движения

Интегрирование уравнений движения производится программой rashet32, полученной на шаге МОДЕЛЬ. В оболочке системы соответствующий пункт меню - ‘Решение уравнений модели - шаг РАСЧЕТ’. При инициировании этого пункта меню требуется указать следующие файлы:

1.
файл с дополнительными данными для шага Расчет, расширение rhl;
2.
файл с начальными условиями интегрирования, расширение ico. При первом запуске расчета это может быть файл default.ico;
3.
файл с параметрами расчета кинематических пар, расширение int. Описание структуры файла - п.п. 9.10.
4.
файл с управляющей информацией для интегрирования уравнений движения - п.п. 3.2.

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

При расчете собственных частот и форм колебаний системы дополнительная информация задается в файле с именем epsilst - см. п.п. 3.3.

После отработки программы интегрирования уравнений образуются следующие файлы:

1.
rezr - рабочий файл с результатами интегрирования, которые используются затем для обработки и вывода;
2.
default.ico - файл с конечным положением интегрируемой системы, который в дальнейшем можно использовать как начальные условия для других расчетов.

1.4.3 Обработка и вывод результатов решения

Описание вида обработки и вывода

Описание вывода требуемых параметров движения конструкции и их характеристик производится в файле с постоянным именем uprf - см. главу 4.

Обработка результатов интегрирования

Обработка результатов интегрирования производится программой anaz по заданному в п.п. 1.4.3.описанию обработки. При работе в оболочке соответствующий пункт меню - ‘Обработка результатов решения’. При инициализации этого пункта меню требуется указать файл с описанием контрольных точек для обработки - расширение upr - п.п. 4.1.

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

Вывод результатов обработки

Вывод результатов обработки производится программой viv. При работе в оболочке соответствующий пункт меню -‘Вывод результатов решения’. После отработки программы вывода результатов образуются следующие файлы:

1.
benafi.bnf - файл с графиками для программы просмотра;
2.
rezob - файл с результатами имеющими табличное представление (алфавитно - цифровая информация).

Обработка и вывод результатов может задаваться и производиться многократно без повторного интегрирования уравнений.

Просмотр графиков

Просмотр графиков производится программой benafin. При работе в оболочке соответствующий пункт меню -‘Просмотр графиков’. При инициализации этого пункта меню требуется указать файл с графиками - расширение bnf - п.п. 4.4.

Анимация движений

Анимация движений производится программой multic. При работе в оболочке соответствующий пункт меню - ‘Анимация движений’. При инициализации этого пункта меню требуется указать файл для аниматора - расширение anc - п.п. 4.5.

1.5 Использование макро в оболочке

При работе в оболочке существует возможность запоминать последовательность вызываемых пунктов меню и их аргументов, именовать эти последовательности и использовать в дальнейшем как одну команду - макро.

Для записи макро следует нажать клавишу ‘w’, внизу экрана появляется надпись ’Запись макро’. После этого можно вызвать несколько пунктов меню с указанием файлов данных, требуемых для работы с этими пунктами. При этом сами пункты не выполняются, а происходит только формирование bat файла с командами, необходимыми для выполнения этого пункта. Для окончания набора макро снова нажимается клавиша ‘w’.

Сформированное таким образом макро содержится в файле с именем frundmac.bat.

Для запуска макро нажимается клавиша ‘m’. Если макро предполагается использовать неоднократно, то файл frundmac.bat имеет смысл переименовать (оставив расширение bat), во избежание затирания его при записи другого макро.

Макро, содержащее пункт меню ‘Решение уравнений модели - шаг РАСЧЕТ’ может вводить необходимые исходные данные из файла rashet.inp, для чего эти данные нужно предварительно внести в указанный файл. К таким данным относятся имена файлов с микропрофилем и признак учета скорости движения транспортного средства в начальных условиях.

Кроме того, такое макро при запуске запрашивает варианты файлов с параметрами расчетной схемы (расширение ras). При нажатии клавиши Esc, происходит отказ от ввода вариантов и данные вводятся из файла, указанного при записи макро. В противном случае нужно указать на предварительно переименованные файлы с вариантами параметров, для которых макро проведет последовательные расчеты. При обработке результатов файлы с графиками будут иметь имена файлов с вариантами параметров с расширением bnf.

Глава 2
ФОРМАТ ОПИСАНИЯ РАСЧЕТНОЙ СХЕМЫ

2.1 Общие замечания

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

2.2 Синтаксис языка описания моделей

2.2.1 Допустимые знаки

0123456789+-. ABCDEFGHIJKLMNOPQRSTUVWXYZ,()*/;=
abcdefghijklmnopqrstuvwxyz :|‘йцукенгшщзхъфыва
проджлэячсмитьбюЙЦУКЕНГШЩЗХЪФЫВАПРОЛДЖЭЯЧСМИТЬБЮ

2.2.2 Определения

1.
Алфавитно-цифровая цепочка (цепочка) - набор допустимых знаков, исключая знак ’пробел’.
2.
Цепочка данных - цепочка, состоящая из знаков 0123456789+-. (то есть цепочка есть целое или вещественное число).
3.
Разделительный знак (делиметр) - разделительный знак применяется для разделения цепочек между собой. Он состоит из одного или нескольких знаков "пробел".
4.
Строка - набор цепочек, разделенных делиметром.

2.2.3 Структура файла исходных данных

1.
строка ’имя задачи’
2.
строка ’код задачи’
3.
блок данных
4.
............
5.
блок данных
6.
строка ’конец файла данных’

Примечание: Последовательность блоков данных определена в п.п. 2.32.4.

2.2.4 Структура блока данных

1.
Строка ’имя блока’ ( первая и единственная в блоке данных )
2.
Строка ’комментарий’ ( строка может быть опущена )
3.
Строка ’комментарий для запоминания’ ( строка может быть опущена )
4.
Строка ’данные’
5.
Строка ’комментарий’ ( строка может быть опущена )
6.
Строка ’комментарий для запоминания’ ( строка может быть опущена )
7.
Строка ’данные’

Примечание:

1.
Строка ’имя блока’ в каждом блоке данных различна.
2.
Строка ’комментарий’ в блоке данных может стоять перед любой строкой.
3.
Строка ’комментарий для запоминания’ в блоке данных требует после себя наличия строки данных.

2.2.5 Типы строк в файле исходных данных

Здесь и далее в квадратных скобках указано то, чем может быть дополнена строка. Длина строки не должна превышать 80-ти знаков. Файл исходных данных является дисковым файлом, состоящим из строк (типы строк описаны нижеше).

Строка ’имя задачи’

Формат строки: LIN1 XXXX имя, где XXXX - код данных. Для шага МОДЕЛЬ XXXX=STPM, для шага РАСЧЕТ XXXX=STPR. Имя (имя задачи)- набор цепочек, заключенных в апострофы (количество знаков между апострофами не более 40). Эта строка должна быть первой в файле исходных данных. Строка для одной модели на шаге МОДЕЛЬ и на шаге РАСЧЕТ должна совпадать.

Пример: LIN1 STPR ’ВАЗ 2123, модель с упругим кузовом’

Строка ’код задачи’

Формат строки: LIN2 имя [цепочки], где имя (имя кода задачи) - цепочка, заключенная в апострофы ( количество знаков между апострофами - 8 ); Эта строка должна быть второй в файле исходных данных. Строка для одной модели на шаге МОДЕЛЬ и на шаге РАСЧЕТ должна совпадать. Имя кода задачи используется для формирования некоторых вспомогательных имен файлов модели, поэтому в нем следует применять символы, допустимые в именах файлов для операционной системы.

Пример: LIN2 ’vaz2123u’

Строка ’имя блока’

Формат строки: [делиметр] BLXX [цепочки], где XX - номер блока. Эта строка стоит перед строками данных.

Строка ’комментарий’

Формат строки: $ [делиметр] [цепочки], знак ’$’ должен быть первым в строке. Эта строка может стоять в любом месте файла исходных данных после строк ’имя задачи’, ’код задачи’. Строка можетотсутствовать.

Строка ’комментарий для запоминания’

Формат строки: @ [делиметр] имя, знак ’@’ должен быть первым в строке, имя - набор цепочек, заключенных в апострофы (количество знаков между апострофами не более 40). Комментарии такого типа используются для формирования комментариев в программах, генерируемых системой. Строка может отсутствовать. В дальнейшем такой комментарий может называться ‘@’- комментарием.

Строка ’данные’

Формат строки: A1 ... AN, где Ai - i-ая цепочка данных. Перед строкой ‘данные’ может стоять строка ‘комментарий для запоминания’ ).

Строка ’конец файла’

Формат строки: ENDE [делиметр][цепочки], эта строка должна быть последней в файле данных.

2.3 Формат описания структуры расчетной схемы модели

Последовательность описания блоков данных, представленная ниже, определяет последовательность их расположения в файле данных.

2.3.1 Блок описания подсистем (тел) модели

Присутствие блока обязательно. Имя блока - BL14. Строка данных - состоит из 8 чисел :

1.
Номер подсистемы (целое число (ц. ч.)), номер подсистемы может принимать значения от 1 до 5000,
2.
Тип подсистемы ( ц. ч. ), см. п.п. 8.2.
3.
Признак наличия степени свободы по оси X (ц. ч.), 0 - степень подавлена, 1 - степень разрешена, N - номер узла упругого тела. Задаваемые таким образом три узла упругого тела, используются для пересчета начальных условий и расчета матрицы поворота. В качестве таких узлов выбираются любые узлы не лежащие на одной прямой и, желательно, максимально удаленные друг от друга.
4.
Признак наличия степени свободы по оси Y (ц. ч.),
5.
Признак наличия степени свободы по оси Z (ц. ч.),
6.
Признак наличия поворота вокруг оси X (ц. ч.), 0 - поворот подавлен, 1 - поворот разрешен, малые движения, 3 - поворот разрешен, большие движения,
7.
Признак наличия поворота вокруг оси Y (ц. ч.),
8.
Признак наличия поворота вокруг оси Z (ц. ч.).

Строк данных столько, сколько подсистем.

Номера подсистем 49 и 50 имеют специальный смысл - 49 - для задания кинематического возмущения от дороги, 50 - резервный. Поэтому их использование для других целей запрещено. При использовании подсистемы с номером 49 для задания кинематического возмущения от двухколейного пути, считается, что первые десять узлов крепления (1-10) принадлежат левой колее, вторые десять (11-20) правой. Сама подсистема номер 49 должна быть твердым телом со всеми запрещенными степенями свободы.

Если подсистемы заданного типа нет в библиотеке тел, то выдается предупреждение и считается, что такое тело имеет 20 узлов крепления с номерами 1-20.

Если в блоке были указаны твердые тела с большими поворотами, то для этих тел автоматически генерируется необходимая дополнительная информация (см. п.п. 9.3.1), которая записывается в файл default.mhl.

2.3.2 Блок описания характеристик соединительных элементов

Присутствие блока обязательно. Имя блока - BL16. Строка данных - состоит из 6-и чисел:

1.
Номер характеристики (ц. ч.),
2.
Номер (направление) составляющей (ц. ч.), 1 - вдоль оси X, 2 - вдоль оси Y, 3 - вдоль оси Z, 4 - поворот вокруг оси X, 5 - поворот вокруг оси Y, 6 - поворот вокруг оси Z, (если нет составляющей - нет и соответствующей строки),
3.
Номер типа упругой составляющей функции (0 - упругая составляющая отсутствует),
4.
Число параметров упругой составляющей функции,
5.
Номер типа демпфирующей составляющей функции (0 - демпфирующая составляющая отсутствует),
6.
Число параметров демпфирующей составляющей.

Типы функций слагаемых - см. описание типов функций - п.п. 6.26.3

2.3.3 Блок описания соединительных элементов (с. э.)

Присутствие блока обязательно. Имя блока - BL17. Строка данных - состоит из 6-и чисел:

1.
Номер соединительного элемента (ц. ч.),
2.
Номер первой подсистемы (ц. ч.),
3.
Номер узла крепления первой подсистемы (ц. ч.),
4.
Номер второй подсистемы (ц. ч.),
5.
Номер узла крепления второй подсистемы (ц. ч.),
6.
Номер механической характеристики.

Строк данных столько сколько с. э..

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

2.3.4 Блок описания специальных сил (с.с.)

Блок может отсутствовать, если имеются мосты. Имя блока - BL18. Строка данных - состоит из 4-х чисел:

1.
Номер специальной силы (ц. ч.),
2.
Номер (направление) составляющей (ц. ч.), 1 - вдоль оси X, 2 - вдоль оси Y, 3 - вдоль оси Z, 4 - поворот вокруг оси X, 5 - поворот вокруг оси Y, 6 - поворот вокруг оси Z,
3.
Номер типа функции составляющей (ц. ч.),
4.
Число параметров функции.

Спец. сила может содержать до шести составляющих. Спец. сил может быть не более 1000.

2.3.5 Блок задания узлов действия специальных сил

Блок может отсутствовать, если отсутствует BL18. Имя блока - BL19. Строка данных - состоит из 3-х чисел:

1.
Номер специальной силы (ц. ч.),
2.
Номер подсистемы (ц. ч.),
3.
Номер узла крепления подсистемы (ц. ч.).

Строк столько, сколько спец. сил описано в BL18. Спец. сила может быть приложена только в узле крепления. Точка центр масс (Ц.М.) - всегда точка крепления. Спец. сила, имитирующая вес, может быть приложена только к подсистеме ’твердое тело’ в Ц.М.

2.4 Формат описания параметров расчетной схемы модели

Последовательность описания блоков данных, представленная ниже, определяет последовательность ввода.

2.4.1 Блок описания инерционных параметров подсистемы твердое тело

Блок может отсутствовать, в этом случае данные будут взяты из библиотеки тел (подсистем). Имя блока - BL31. Строка данных - состоит из 5-ти чисел:

1.
Номер подсистемы (ц. ч.),
2.
Масса тела (вещественное число (в. ч.)),
3.
Главный центральный момент инерции относительно оси X,
4.
Главный центральный момент инерции относительно оси Y,
5.
Главный центральный момент инерции относительно оси Z.

Строк данных столько, сколько подсистем.

2.4.2 Блок геометрических параметров точек крепления подсистемы твердое тело

Блок необходим только при изменении координат точек тел. При отсутствии описания координат точки тела в блоке данные будут взяты из библиотеки тел (подсистем). Имя блока - BL32. Строка данных - состоит из 5-ти чисел:

1.
Номер подсистемы (ц. ч.),
2.
Номер узла крепления в подсистеме (ц. ч.),
3.
Координата X точки крепления (в. ч.),
4.
Координата Y точки крепления (в. ч.),
5.
Координата Z точки крепления (в. ч.).

Строк данных столько, сколько точек крепления. Координаты задаются в локальной системе координат в метрах. Вводятся не сами координаты точек а смещения точек относительно ц. м. тела.

2.4.3 Блок описания параметров функций сил в механической характеристике

Присутствие блока обязательно. Имя блока - BL33. Строка данных - состоит из 4-х чисел:

1.
Номер механической характеристики (ц. ч.),
2.
Номер (направление) составляющей (ц. ч.), 1 - вдоль оси X, 2 - вдоль оси Y, 3 - вдоль оси Z, 4 - поворот вокруг оси X, 5 - поворот вокруг оси Y, 6 - поворот вокруг оси Z,
3.
Признак составляющей функции (ц. ч.), 1- упругая составляющая, 2- демпфирующая составляющая
4.
Количество строк по 4 в. ч. в строке, которые вводятся после описанной строки. В этих строках задаются параметры функции. Если число параметров не кратно 4 , то последняя строка дополняется нулями.

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

2.4.4 Блок описания параметров функции специальных сил

Блок может отсутствовать. Имя блока - BL34. Строка данных - состоит из 3-х чисел:

1.
Номер специальной силы (ц. ч.),
2.
Номер (направление) составляющей (ц. ч.), 1 - вдоль оси X, 2 - вдоль оси Y, 3 - вдоль оси Z, 4 - поворот вокруг оси X, 5 - поворот вокруг оси Y, 6 - поворот вокруг оси Z,
3.
Количество строк по 3 в. ч. в строке, которые вводятся после описанной строки. В этих строках заданы параметры функции. Если число параметров не кратно трем, то последняя строка дополняется нулями.

В строках должны быть описаны параметры всех специальных сил, заданных в блоке 19.

2.4.5 Блок описания коэффициентов демпфирования по формам подсистемы упругое тело

Блок может отсутствовать если в модели нет упругих тел. Имя блока - BL35. Строка данных - состоит из 5-ти чисел:

1.
Номер подсистемы (ц. ч.),
2.
Количество строк по 5 в. ч. в строке ,которые вводятся после описанной строки. В этих строках заданы коэффициенты демпфирования по формам колебаний тела. Если число коэффициентов демпфирования не кратно 5, то последняя строка дополняется нулями,
3.
Число вводимых коэффициентов демпфирования (ц. ч.),
4.
Любое целое число,
5.
Любое целое число.

Строк данных столько, сколько упругих подсистем.

2.4.6 Блок задания расстояния между первым и последующим мостом

Блок может отсутствовать. Имя блока - BL36. Строка данных - состоит из 5-ти чисел:

1.
Номер подсистемы передний мост (ц. ч.),
2.
Номер узла контакта с дорогой переднего моста (ц. ч.),
3.
Номер подсистемы следующий мост (ц. ч.),
4.
Номер узла контакта с дорогой следующего моста (ц. ч.),
5.
Расстояние между этими мостами (в. ч.).

Строк данных столько сколько мостов, не считая передний мост.

2.4.7 Блок описания параметров колеи

Блок может отсутствовать. Имя блока - BL37. Строка данных - состоит из 6-ти чисел:

1.
Признак колеи ( ц. ч. ), 1 - левая колея, 2 - правая колея.

Следующие пять чисел в зависимости от вида возбуждения имеют следующий смысл:

1.
случайный профиль.
(a)
Шаг дискретизации реализации микропрофиля (в метрах),
(b)
Величина радиального биения шины (в метрах),
(c)
Радиус качения шины (в метрах),
(d)
Фаза радиального биения шин передних колес,
(e)
Фаза радиального биения шин задних колес.
2.
периодический профиль (все параметры задаются в метрах):
(a)
Амплитуда волны,
(b)
Длина волны,
(c)
Фаза волны,
(d)
Любое вещественное число,
(e)
Любое вещественное число.
3.
единичная неровность (все параметры задаются в метрах).
(a)
Высота неровности,
(b)
Длина подъема неровности,
(c)
Длина спуска неровности,
(d)
Длина пути до наезда на неровность,
(e)
Длина горизонтального участка неровности.

Строк данных 2 - левая и правая колеи.

2.5 Размерность вводимых данных

Система не предопределяет жестко пользование только одной системой единиц. Однако, вся выводимая информация сопровождается комментариями в системе СИ, поэтому предпочтительнее пользоваться ею.

МАССА - кг (килограмм), ДЛИНА - м (метр), ВРЕМЯ - с (секунда), СИЛА - Н (ньютон).

2.6 Препроцессор ввода описания модели

Для упрощения описания структуры и параметров моделей, вводятся дополнительные возможности, которые обеспечиваются препроцессором. Программа препроцессора имеет имя crmodelm.exe для обработки данных шага МОДЕЛЬ и crmodelr.exe для обработки данных шага РАСЧЕТ. Имя входного файла для препроцессора содержится в файле с именем crmodelp.tmp. Выходной файл препроцессора имеет постоянное имя modelp.mdl для шага МОДЕЛЬ или rashetp.ras для шага РАСЧЕТ. Эти имена заносятся в файл crmodel.tmp из которого они считываютя при вводе. Рассмотрим основные возможности препроцессора.

2.6.1 Глобальные переменные

Численные значения параметров уравнений движения могут задаваться с помощью имен-идентификаторов, которым предварительно присвоено значение. Имя может встречаться после своего определения произвольное количество раз. Синтаксис строки определения константы: имя = значение, где ‘имя’ - символьная строка без пробелов, ’значение’ - арифметическое выражение из констант и ранее определенных имен. В качестве символов для имени могут использоваться только буквы латинского алфавита. В выражении допускаются операции сложения, вычитания, умножения, деления, а также левая и правая круглые скобки. Кроме этого предусмотрены функции SIN, COS, SQRT (квадратный корень). Определенная константа может использоваться в строках блоков описания, например:

 
.................................................  
$ параметры оптимизации  
 ad=2  
 massa_kol_pary  = ad*1650.  
 jxz_kol_pary    = 1800.  
.................................................  
BL31 NPDS  MAS           JNX     JNY     JNZ  
$   ’1 колесная пара’  
  1 massa_kol_pary  jxz_kol_pary 200. jxz_kol_pary  
$   ’2 колесная пара’  
  2 massa_kol_pary  jxz_kol_pary 200. jxz_kol_pary*0.5  
.................................................  

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

2.6.2 Вставки

В файле допускаются вставки. Для этого в первую позицию строки заносится символ #, за которым без пробела следует имя файла. Например:

#fn1.ras

Подразумевается, что в файле fn1.ras содержится какая-то информация для препроцессора. Если имя файла макро имеет расширение mhl (файл с дополнительной информацией для шага МОДЕЛЬ), то его содержимое после замены констант выводится в файл с постоянным именем modhelpp.mhl. Если имя файла макро имеет расширение rhl (файл с дополнительной информацией для шага РАСЧЕТ), то его содержимое после замены констант выводится в файл с постоянным именем rashelpp.rhl.

2.6.3 Условный оператор

Препроцессор содержит условный оператор типа

if-then-else

Формат оператора:

if условие then  
  (then  часть)  
else  
  (else часть)  
endif

Условие записывается в форме: выражение1 операция_сравнения выражение2.

Допускаются следующие операции сравнения:
== - равно,
!= - не равно,
> - больше,
< - меньше.

Пример:

a=1.  
b=2.  
if a*2==b then  
  c=10.  
else  
  c=20.  
endif

Если выражение не определено, то его значение считается равным 0.

2.6.4 Макроопределения

Макроопределения или просто макро - параметризованные фрагменты описания отдельных частей модели, позволяющие упростить описание сложных моделей. Синтаксис макро - аналогично вставке, только имя файла должно иметь расширение *.mac, за которым через пробелы должны быть перечислены фактические параметры макро. Например:

#rail.mac 1 a1 49 1 massa

Тексты макро хранятся в в одноименных файлах, которые собраны в директории
../UTILF/MACRO. Количество параметров макро и вложенность не ограничиваются. Формальные параметры в тексте макро именуются идентификаторами вида pXX, где XX порядковый номер формального параметра в обращении ( например p1, p12 и т.д.). Перед строкой вызова макро должен стоять ’@’ комментарий, задающий имя уникальное имя данного экземпляра макро, и который может использоваться в макро специальным образом - см. ниже. Вызов макро должен стоять после описания всех блоков. Макро может содержать описания блоков как для шага МОДЕЛЬ, так и РАСЧЕТ, смысл формальных параметров при этом различен для блоков разных шагов. Все идентификаторы, объявленные в макро носят локальный характер и теряют свои значения после выхода.

Системные переменные макро

Для согласования нумерации основных типов элементов, применяемых для описания моделей во ФРУНДе, в макро и в основной программе (файле описания рсчетной схемы) используются системные переменные, значения которых автоматически вычисляются препроцессором. Такие системные переменные обозначаются идентификаторами, начинающимися со следующих символов:

npt - номера подсистем твердых или упругих тел; npu - номера подсистем объектов контуров управления; nmx - номера механических характеристик; nse - номера соединительных элементов; nss - номера сил нагрузок - спец. сил.

Например npt1, nmx22 и т.д.

Системные функции макро

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

w_mhl2 - запись в файл default.mhl строчки указания, что соединительный элемент является "пружиной в шарнирах". Функция использует 1 параметр: номер соединительного элемента.

w_mhl3 - запись в файл default.mhl строчки указания, о направлении оси вращения колеса и точки его контакта с доргой. Использует 4 параметра:

1.
номер подсистемы тела колесо,
2.
номер узла контакта колеса с дорогой,
3.
номер узла на оси колеса в плоскости колеса,
4.
номер узла на оси колеса.

w_rhl3 - запись в файл default.rhl величины свободной длины пружины в шарнире. Функция использует 2 параметра:

1.
номер соединительного элемента,
2.
длина пружины.

w_ktcom - запись в файл defaulm.upr комментария к контрольной точке. Текст комментария должен задаваться как ’@’ комментарий в строчке, предшествующей вызову макро.

w_copyf - копирование файлов. Имена файлов задаются последовательно в ’@’ комментарии в строчке, предшествующей вызову макро. Первым задается имя файла который нужно скопировать, вторым имя файла в который копируется.

w_kt - запись в файл defaulm.upr контрольной точки описания выводимого параметра. Функция использует следующие параметры:

1.
класс вывода;
2.
код вывода;
3.
номер узла тела или соединительного элемента;
4.
номер тела;
5.
признак направления (1-6);
6.
признак величины (1-3).

Стандартные макро

В стандартную конфигурацию ФРУНДа входят макро, необходимые для полного описания модели из файлов *.mdl и *.ras. К таким макро относятся:

tspring.mac - задание признака пружины в шарнире для шага МОДЕЛЬ. Макро имеет один параметр: номер соединительного элемента.

lspring.mac - задание свободной длины пружины в шарнире для шага РАСЧЕТ. Макро имеет два параметра:

1.
номер соединительного элемента,
2.
длина пружины в свободном состоянии.

Если данное макро вызывается из другого макро, то номер соединительного элемента определяется автоматически, а на месте фактического параметра может стоять любое значение.

wkt.mac - генерация контрольной точки для вывода результатов. Параметры - см. параметры системных функции w_kt и w_ktcom.

copyf.mac - копирование файлов. Имена файлов задаются последовательно в ’@’ комментарии в строчке, предшествующей вызову макро. Первым задается имя файла который нужно скопировать, вторым - имя файла в который копируется.

kkoleso.mac - запись в файл default.mhl строчки указания, о направлении оси вращения колеса и точки его контакта с дорогой. Использует 4 параметра:

1.
номер подсистемы тела колесо,
2.
номер узла контакта колеса с дорогой,
3.
номер узла на оси колеса в плоскости колеса,
4.
номер узла на оси колеса.

Кроме указанных набор макро, хранящийся в директории ../UTILF/MACRO может содержать параметризованные описания специфических элементов - взаимодействие колесной пары с рельсом, листовую рессору и т.д. Описание которых содержится в комментариях, хранящихся в самих макро. Часть таких макро описана в главе 11.

Глава 3
УПРАВЛЯЮЩАЯ ИНФОРМАЦИЯ ДЛЯ ПРОГРАММЫ ИНТЕГРИРОВАНИЯ УРАВНЕНИЙ ДВИЖЕНИЯ

3.1 Общие положения

Управляющая информация, необходимая для интегрирования уравнений математической модели, задается в файле с постоянным именем uprf, при работе в оболочке файлы с расширением upr. Первая строка указанного файла содержит значения управляющих переменных в свободном формате, вторая строка - обязательный комментарий (опускаться не может!) содержит, например, символическое обозначение управляющих переменных. Начиная с третьей строки могут задаваться описания контрольных точек для вывода результатов - см. главу 4.

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

3.2 Управляющие переменные и формат их задания

Управляющие переменные интегрирования задаются в одной первой строке в следующей последовательности:

1.
Количество вариантов расчета - целое число. Вариантом расчета считается расчет при определенной скорости движения транспортной машины или частоте возмущающих сил. Число вариантов не ограничивается.
2.
Начальная скорость (частота) варианта - вещественное число. Скорость задается в м/с, частота в герцах.
3.
Шаг изменения скорости (частоты) - в тех же единицах, что и в предыдущем пункте.
4.
Шаг интегрирования уравнений - вещественное число. Шаг задается в секундах.
5.
Шаг записи результатов интегрирования в промежуточный файл хранения - вещественное число. Шаг записи результатов должен быть не меньше шага интегрирования.
6.
Код расчета - целое число. Код расчета определяет тип возмущения, действующего на конструкцию. Установлены следующие коды расчета:
(a)
код 1 - случайное воздействие со стороны дороги.
(b)
код 2 - периодическое воздействие со стороны дороги.
(c)
код 3 - воздействие со стороны дороги от единичной неровности
(d)
код 4 - воздействие периодической силой, непосредственно приложенной к какой либо точке конструкции, кинематическое воздействие со стороны дороги при этом обнуляется.
(e)
код 5 - возмущение специальными силами приложенными в точках конструкции.
(f)
код 7 - линеаризация и расчет собственных частот и форм колебаний. Расчет собственных частот производится для конечного момента времени интегрирования.
(g)
код 8 - сборка механизмов.
(h)
код 9 - определение статического положения системы.
(i)
код 10 - расчет управляемого движения системы методом обратной задачи.
7.
Длина пути машины, на котором интегрируются уравнения движения, или время интегрирования (см следующий пункт) - действительное число. Задается в метрах или в секундах. Время интегрирования при задании пути определяется отношением пути к скорости.
8.
Длина пути или время, с которого начинается запись результатов интегрирования в промежуточный файл хранения. По кодам 1 - 3 в п.п. 7,8 задается путь, по кодам 4,5 - задается время. По кодам 1 - 3 в п.п. 2, 3 задается скорость, по кодам 4,5 - частота.
9.
Количество вариантов по начальным условиям - целое число. В данной версии всегда задается единица.
10.
Код накопления результатов интегрирования - целое число. Код накопления указывает, какие параметры, получаемые при интегрировании, необходимо запомнить для дальнейшей обработки.

Установлены следующие коды накопления:
1 - накапливаются все результаты интегрирования,
2 - накапливаются перемещения, скорости, ускорения обобщенных координат,
3 - накапливаются только перемещения обобщенных координат,
4 - накапливаются только скорости обобщенных координат,
5 - накапливаются только ускорения обобщенных координат,
6 - накапливаются относительные перемещения, скорости в с.э. и силы в них,
7 - накапливаются только те переменные, которые описаны в контрольных точках в файле uprf. Этот способ требует минимального объема дисковой памяти.

3.2.1 Задание типа процедуры интегрирования

В системе моделирование подключено две процедуры интегрирования - метод Рунге - Кутта четвертого порядка и метод Чанга - метод первого порядка. Оба метода явные. Неявные методы (метод Гира) не включены в стандартную конфигурацию пакета, поскольку при тестовых расчетах для больших моделей не показали преимуществ.

Тип процедуры интегрирования задается в коде расчета. При обычном задании кода расчета используется метод Рунге - Кутта. Для задания метода Чанга необходимо код расчета увеличить на 10. Например: код расчета 1 - метод Рунге Кутта, код расчета 11 - метод Чанга.

3.3 Управляющая информация для расчета собственных частот и форм

При расчете собственных частот и форм используется дополнительная информация, которая заносится в файл epsilst. В указанном файле в первой строке идут последовательно:

1.
отклонение координаты по перемещению при линеаризации, по умолчанию - 0.0001,
2.
отклонение координаты по скорости при линеаризации, по умолчанию - 0,
3.
коэффициент кратности условной жесткости кинематических пар. Условная жесткость в кинематических парах определяется произведением максимального диагонального коэффициента матрицы жесткости упругой части конструкции на коэффициент кратности. Значение по умолчанию - 100.

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

Укажем, что линеаризация нелинейности типа рессоры НАМИ (п.п. 6.2) производится по внешней характеристике на текущем участке. При построении матрицы жесткости программа запрашивает наличие в модели форм перемещений как твердого тела, и в случае утвердительного ответа запрашивает корректировочный множитель, который используется для корректировки матрицы. Корректировочный множитель имеет величину порядка низшей собственной ненулевой частоты колебаний. При выводе частоты приводятся к частотам исходной матрицы. Корректировка матрицы может производиться так же при наличии нулевых диагональных членов, такие члены в этом случае заменяются значениями равным минимальному диагональному коэффициенту, деленному на 100.

При наличии в системе кинематических связей, такие связи заменяются упругими связями, с жесткостями, равными максимальному диагональному коэффициенту, умноженному на коэффициент кратности (по умолчанию 100). При выводе частот в этом случае выводятся только частоты, соответствующие реальным степеням свободы системы. Если при расчете получается отрицательная собственная частота, то выдается соответствующее сообщение, а самой частоте присваивается малое положительное значение - 0.00008. Значения частот и соответствующие им формы можно посмотреть, построив анимацию для форм колебаний.

3.4 Управляющая информация для расчета моделей с кинематическими парами

При интегрировании уравнений с кинематическими парами используется информация из файла epsilon.dat - см. главу 9, при работе в оболочке такие файлы имеют расширение int. Если не используется схема стабилизации решения (первые три числа в файле равны 0), то для получения допустимого по точности решения, следует задать только погрешность смещения в шарнирах - пятый параметр в файле. Рекомендуемое значение погрешности - 0.0001. В случае превышения смещения в шарнирах заданного расчет будет останавливаться с сообщением "раскрытие шарниров". В таких случаях следует использовать стабилизацию решения в соответствии с изложенным в главе 9.

3.5 Управляющая информация для сборки моделей с кинематическими парами

При работе с моделями с кинематическими парами для начала решения необходимо знать начальное положение механизма - начальные условия. Для их определения служит код расчета 8 - сборка механизмов. Для задания параметров сборки используется информация из файла epsilon.dat - см. главу 9. При этом существенное значение имеют следующие параметры:
1 - коэффициент стабилизирующих консервативных сил,
2 - коэффициент стабилизирующих диссипативных сил,
5 - максимальная допустимая погрешность в шарнирах,
6 - код схемы стабилизации, 1-улучшенная, 0-обычная,
9 - код интегрирования уравнений Эйлера,
10 - коэффициент условия прекращения сборки из-за несходимости (целое число >= 1),
12 - шаг времени пересчета погрешности в шарнирах,
13 - шаг подсоединения уравнений связей при сборке,
14 - условие печати вычисления погрешности в шарнирах (0 есть печать, 1 нет).

Алгоритмически сборка производится путем последовательного присоединения звеньев механизма или отдельных связей , что определяется 13-ым параметром файла epsilon.dat. Процесс сборки можно наблюдать, построив анимационное изображение движения системы при сборке. Сборка может не сходиться по двум причинам - из-за несобираемости самого механизма и из-за того, что в процессе сборки механизм попадает в особое положение, из которого не может выйти. Причина эта может быть установлена при просмотре процесса сборки на аниматоре. В первом случае исправляется расчетная схема модели или ее параметры. Во втором - строится модифицированная модель, обеспечивающая достижение одного из возможных положений реальной модели, но содержащая измененные кинематические пары, возможно и с избыточными связями. Например, при использовании только сферических шарниров сборка будет сходиться всегда (при условии собираемости самого механизма). Файл начальных условий, получаемый в конце сборки содержит одно из возможных начальных положений механизма, с которого можно начинать решение динамической задачи. Сборка может проводиться поэтапно, параллельно с усложнением модели, путем добавления в нее новых тел. При этом каждое конечное положение предыдущей модели служит начальным для последующей.

3.6 Задание начальных условий

Начальные условия по умолчанию считаются нулевыми. Специально начальные условия могут задаваться из файлов с начальными условиями, при работе в оболочке - файлы с расширением ico.

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

При каждом запуске rashet32 программа проверяет наличие файла с именем initcond.inp, и, если он существует, то считывает из него начальные условия, считая, что они записаны в формате файла initcond.ico. Другими словами, в файл initcond.inp копируется содержимое нужного файла initcond.ico. Файл initcond.ico содержит записи с номерами тел и начальными условиями для этих тел. Файлы указанного формата могут использоваться для разных моделей, если в этих моделях нумерация идентичных тел совпадает. То есть, начальные условия для более простой модели могут использоваться для более сложной модели, если эта простая модель целиком входит в сложную и имеет одинаковую нумерацию общих тел и совпадающие степени свободы у таких тел. Начальные условия для моделей с упругими телами находятся в два этапа. Сначала определяются начальные условия для модели в которой упругие тела описаны как твердые. Затем в полученном файле с начальными условиями модифицируются номера упругих тел - номера записываются со знаком минус. Такой модифицированный файл однократно используется для расчета модели с упругими телами. Полученные после такого расчета начальные условия можно в дальнейшем использовать при расчетах с упругими телами.

3.7 Управляющая информация для определения статического положения системы

При расчетах систем с упругими элементами и нагруженных постоянными силами (силами веса) для начала расчета необходимо определять статическое положение системы. Для этого используется код расчета 9. Условием нахождения статического положения является уменьшение кинетической энергии системы в максимуме ниже пороговой величины, определяемой параметром 5 в файле epsilon.dat - п.п. 3.5. Статическое положение системы выводится в файл начальных условий - initcond.ico.

Глава 4
ОБРАБОТКА И ВЫВОД РЕЗУЛЬТАТОВ РЕШЕНИЯ

4.1 Задание контрольной точки

Описание вида обработки задается в постоянном файле uprf (файлы с расширением upr в оболочке), в котором содержится и управляющая информация для интегрирования уравнений. Задание вида обработки производится парами строк, начиная с третьей строки файла. Каждая пара строк идентифицирует один объект - контрольную точку. В первой строке каждой пары задается имя контрольной точки, содержащее не более сорока символов, во второй строке пары задаются последовательно следующие параметры:

1.
1. Класс вывода: 1 - 8 (см. описание классов вывода).
2.
2. Код вывода (см описание кодов вывода.)
3.
3. В зависимости от класса вывода задаются: 1) класс 1,5,6,7,8. Номер узла подсистемы твердое тело или упругое тело, описанные в исходных данных шага МОДЕЛЬ. 2) класс 2,3. Номер соединительного элемента. 3) класс 4. Порядковый номер обобщенной координаты упругой подсистемы.
4.
4. Для классов вывода 1,4,5,6,7,8 задается номер подсистемы. Для классов вывода 2 и 3 - любое целое число.
5.
5 - 10. В этих позициях задаются признаки вывода параметров по нужным направлениям, соответственно по X,Y,Z и по угловым направлениям UX,UY,UZ. Признаком вывода является внесение в соответствующую позицию целого числа, отличного от нуля. Для класса вывода 4, нужно задавать один признак в любой из позиций 5 - 10.
6.
11 - 13. В этих позициях задаются следующие признаки: 1) классы вывода 1,6. Задаются признаки вывода соответственно перемещений, скоростей, ускорений - признаком является целое число отличное от нуля. 2) класс вывода 2. Задаются признаки вывода соответственно относительных перемещений и скоростей в соединительных элементах. Позиция 13 всегда игнорируется. 3) класс вывода 3. Задаются признаки вывода соответственно упругого, демпфирующего слагаемых сил в соединительном элементе и суммарной силы. 4) класс вывода 4. Аналогично классу вывода 1. 5) класс вывода 5. Задается признак вывода только в позиции 13.
7.
14. Задается признак масштабирования осей графиков при графическом выводе результатов: 0 - линейный масштаб, 1 - логарифмический масштаб. (логарифмический масштаб корректно работает для кодов вывода, работающих со спектрами и их производными).
8.
15 - 17. Задаются параметры только для класса вывода 1: 1) подсистема твердое тело - задаются координаты контрольной точки относительно ц.м. подсистемы. 2) подсистема упругое тело - задаются номер первой учитываемой формы колебаний и количество форм, которые нужно учесть. 17 - ая позиция всегда игнорируется. Если задаются все нули, то учитываются все формы.
9.
18 - 20. Задаются вспомогательные параметры для обработки в зависимости от указанного кода вывода - см. описание кодов вывода. При задании вида обработки и вывода следует помещать пары строк, описывающие контрольные точки с одним классом и кодом вывода, последовательно, друг за другом. Символ ’$’ в первой позиции строки файла задания вида обработки является признаком комментария.

4.2 Описание классов вывода.

Класс вывода определяет качественный вид параметра, характеристики которого следует обработать и вывести. Установлены следующие классы вывода:

1.
Абсолютные перемещения, скорости, ускорения точек подсистем упругое тело и твердое тело.
2.
Относительные перемещения и скорости в соединительных элементах.
3.
Упругие, демпфирующие и суммарные силы в соединительных элементах.
4.
Перемещения, скорости, ускорения обобщенных координат упругих подсистем.
5.
Третья производная перемещений точек подсистем.
6.
Абсолютные перемещения, скорости, ускорения точек упругого тела, совершающего большие движения.
7.
Скорости и ускорения тела, совершающего большие движения в проекции на подвижную систему координат. Такой класс вывода необходим, например, для определения скорости движения автомобиля и бокового ускорения при расчете управляемости.
8.
Радиус поворота автомобиля.

Примечание: 7 и 8 классы действуют только для кодов вывода 1 и 2.

4.3 Описание кодов вывода

Код вывода определяет вид обработки тех или иных параметров движения системы. Установлены следующие коды вывода:
Код вывода 1 - графики реализаций процессов. В позициях 18-20 задаются начальное, конечное значение интервала времени, в котором необходимо построить график, и шаг абсциссы графика в секундах.
Код вывода 2 - графики спектров мощности процессов в интересующем диапазоне частот. В позициях 18-20 задаются начальная, конечная частота спектра и шаг (разрешение) спектра в герцах. Размерность спектра - (размерность параметра)^2/Гц. Спектр вычисляется с помощью БПФ (быстрого преобразования Фурье). Реализация разбивается на блоки с длиной кратной степени двойки и наложением 75%. Результирующий спектр усредняется по блокам. Минимальное разрешение спектра определяется длиной блока. Если заданный шаг разрешения меньше минимального, то при выводе он назначается равным минимальному. Если заданный шаг разрешения больше минимального, то он округляется до наименьшего целого. При вычислении спектров используются временные окна. Длина блока под реализации и вид окна определяются в файле windef в котором последовательно задаются код вида окна и степень двойки длины под реализации. Определены следующие коды вида окон:

1.
- трапециевидное окно,
2.
- окно Хэннинга,
3.
- окно Кайзера - Бесселя,
4.
- прямоугольное окно,
5.
- окно Flat - Top.

По умолчанию устанавливается окно Хэннинга и длина под реализации 1024 точки - степень двойки 10. Настройки, задаваемые в файле windef действуют на все коды вывода, использующие БПФ.
Код вывода 3 - таблицы среднеквадратичных значений процессов во всем диапазоне частот. В позициях 18-20 задаются нули.
Код вывода 4 - таблицы среднеквадратических значений скоростей и ускорений в треть октавных полосах частот по сравнению с нормами. В позиции 18 задается число октавных полос(максимум 7), в позициях 19, 20 - нули.
Код вывода 5 - график функции когерентности двух процессов. В позициях 18 - 20 задаются начальное, конечное значение частотного интервала и шаг по частоте. По данному коду каждая функция когерентности описывается двумя парами строк, причем в каждой строке в позициях 5 - 10 и 11 - 13 должно стоять только по одному признаку. Комментарий будет выводиться только из первой пары строк. Каждая пара определяет один из процессов, функции когерентности которых определяются. Код вывода 6 - таблицы характеристик выбросов процессов за интересующий уровень. В позиции 18 задается значение уровня, в позиции 19 - признак направления выброса - положительное число - выброс вверх, отрицательное - вниз. В позиции 20 задается нуль.
Код вывода 7 - графики спектров мощности процессов в интересующем диапазоне частот со сглаживанием по треть октавам. Позиции 18-20 задаются аналогично коду 2.
Код вывода 8 - гистограммы распределения процессов. В позициях 18-20 задаются нули.
Код вывода 9 - таблицы мощностей, рассеиваемых в диссипативных элементах системы (только для класса вывода 3). В позициях 18-20 задаются нули, в позициях 11 - 12 задаются всегда нули, в позиции 13 - отличное от нуля число.
Код вывода 10 - графики мнимых и действительных частей амплитуд процессов при периодическом возбуждении. В позициях 18-20 задаются нули. Этот код вывода следует применять только при коде расчета 4.
Код вывода 11 - графики мнимых и действительных частей амплитуд процессов при периодическом возбуждении. В позициях 18-20 задаются нули.
Код вывода 12 - фазовые портреты процессов. В позициях 18-20 задаются нули.
Код вывода 13 - графики среднеквадратических значений скоростей и ускорений в треть октавных полосах частот по сравнению с нормами. В позициях 18 - 20 задаются значения, аналогичные коду 4.
Код вывода 14 - график двумерной плотности распределения вероятности двух процессов. Возможна одновременная обработка только одной пары процессов. Описание пар - аналогично коду 5.
Код вывода 15 - таблицы квантилей процессов. В позициях 18-20 задаются нули.
Код вывода 16 - таблицы спектров мощности скоростей и ускорений в треть октавных полосах частот по сравнению с нормами. В позиции 18 задается число октавных полос, в позициях 19, 20 - нули. Код вывода 17 - графики спектров мощности процессов в треть октавных полосах частот. В позициях 18-20 задаются значения, аналогичные коду 16.
Код вывода 19 - подготовка данных для аниматора. В позициях 18-20 задаются начальное, конечное значение интервала времени, и шаг по времени. Код работает только с классом вывода 1.
Код вывода 20 - подготовка данных для анимации сил. В позициях 18-20 задаются начальное, конечное значение интервала времени, и шаг по времени. Код работает только с классом вывода 3.
Код вывода 21 - сброс необработанных реализаций в ASCII файлы для дальнейшей обработки другими программами. Данные сбрасываются в файлы - fcN.dat, frN.dat, где N - номер класса вывода. fcN.dat содержит имена процессов и названия размерностей. frN.dat - шаг дискретизации (первое число), затем группы значений реализаций. Код работает для классов вывода 1,2,3.
Код вывода 25 - графики передаточной функции. Задание контрольных точек и параметров аналогично коду вывода 5.
Код вывода 26 - графики фазы двух процессов. Задание контрольных точек и параметров аналогично коду вывода 5.
Код вывода 27 - графики взаимных спектров. Задание контрольных точек и параметров аналогично коду вывода 5.
Код вывода 33 - построение СКЗ процессов в зависимости от скорости.
Код вывода 38 - построение максимальных значений процессов в зависимости от скорости.
Код вывода 49 - подготовка данных для анимации собственных форм колебаний. В позиции 18 задается количество выводимых форм (по умолчанию все формы модели), в позиции 19-количество выводимых положений для одной формы (по умолчанию 2). Код вывода применяется только для кода рас чета 7.

4.4 Визуализация графических зависимостей

Программа визуализации benafin предназначена для построения графиков. Считается, что графики предварительно рассчитаны с помощью программы viv и соответствующие данные записаны в файле в требуемом формате - файлы с расширением bnf. Программа позволяет изображать одновременно один график из файла. В файле может быть до 50 графиков. К программе можно подключить до 10 файлов с данными (имена файлов могут повторяться).

4.4.1 Вывод графиков на экран

В начале работы программа запрашивает имя файла с данными для построения графиков ( по умолчанию файл имеет имя benafi.bnf). Если необходимо отказаться от ввода имени файла то нужно ввести ^ Z. После ввода имени файла программа открывает его и полностью прочитывает (операция чтения файла занимает некоторое время), при этом происходит проверка правильности задания данных и, если будет обнаружена ошибка в данных, то весь файл закрывается как непригодный для работы. При просмотре графиков возможен переход к последующему графику и к предыдущему (PgDn и PgUp), а также установка на первый и последний графики в файле (Home и End). Для подключения еще одного файла с данными необходимо нажать клавишу [O/o], и ввести (выбрать) имя файла. Переключение между графиками осуществляется по вводу номера канала, к которому подсоединен файл (номера каналов 0...9), путем нажатия клавиши с цифрой номера канала. Информацию о том, к какому каналу подсоединен файл можно получить из таблицы подсказки, которая вызывается при нажатии клавиши [H/h].

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

При просмотре графиков можно производить следующие операции:


клавиша

операция

0...9

переключение на канал номер 0 - 9,

PgUp/PgDn

следующий/предыдущий график,

Home/End

первый/последний график,

A/a

наложить на текущий график другой график,

W/w

построить зависимость функции первого графика от функции второго,

P

сбросить график в файл формата HPGL, дополнительно в файл text.txt при этом, сбрасываются массивы абсцисс и ординат графика, которые можно использовать для редактирования другими программами, например Excel,

p

сбросить график в файл в формате матричного принтера,

R/r

записать график в файл с другим именем (расширением bnf),

E/e

редактировать надписи комментариев графика,

S/s

редактировать надписи осей графика,

M/m

редактировать границы графика (масштаб),

L/l

выдать координаты указанной точки графика,

K/k

выдать кратные частоты указанной точки графика (команда применяется после команды L/l),

H/h

выдать список команд,

O/o

открыть новый файл с графиками по следующему каналу,

G/g

логарифмический масштаб по обеим осям,

C/c

линейный масштаб по обеим осям,

D/d

интеграл/СКО графика,

B/b

интеграл/СКО графика в указанном диапазоне,

Z/z

нормированная свертка спектра вибраций по оси Z,

x/y

нормированная свертка спектра вибраций по осям X/Y,

N/n

определение параметров затухания свободных колебаний,

^p

сжать поле графика,

^n

разжать поле графика,

^w

черно - белое/цветное изображение,

+-*/

арифметические операции с графиками,

Alt q

квадратный корень из суммы квадратов двух графиков,

Alt + Alt *

сложение/умножение графиков на константы,

Alt c

умножение первого графика на косинус второго. При умножении передаточной функции на косинус фазы получаем коэффициент жесткости,

Alt s

умножение первого графика на синус второго. При умножении передаточной функции на синус фазы получаем коэффициент вязкого сопротивления,

F/f

умножение/деление на аргумент*6.2832,

V/v

дифференцирование/интегрирование спектра,

I/i

двойное дифференцирование/интегрирование спектра,

Q/q

квадратный корень из графика,

^o

закрыть текущий канал,

^z

выйти из программы,

Alt x

выйти из программы.


4.4.2 Формат представления данных для построения графика

Данные, необходимые для построения графика записаны в файл в следующем порядке (бесформатный вывод программой на языке FORTRAN):

Данные по спецификации поля графика

Первая запись: символ ’%’, обозначающий начало каждого графика.

Вторая запись:

1.
код поля - 1 - линейный масштаб, 2 - логарифмический масштаб;
2.
код масштаба (всегда 1);
3.
минимальное значение по оси X;
4.
максимальное значение по оси X;
5.
минимальное значение по оси Y;
6.
максимальное значение по оси Y.

Название графика

Третья запись: имя расчетной задачи (char *80).

Четвертая запись: имя контрольной точки (char *80).

Имена осей графика

Пятая запись: имя оси X (char *32).

Шестая запись: имя оси Y (char *32).

Обозначение кривой графика

Первая запись: символ ’@’, обозначающий начало каждой кривой данного графика.

Вторая запись: имя кривой (char *80).

Данные по спецификации кривой

Третья запись:

1.
тип кривой (1 - равномерный шаг по оси X, 2 - неравномерный шаг по оси X);
2.
количество точек в кривой;
3.
минимальное значение по оси X;
4.
максимальное значение по оси X;
5.
минимальное значение по оси Y;
6.
максимальное значение по оси Y;
7.
шаг по оси X (если тип кривой - 2, значение игнорируется).

Значения координат точек кривой

В каждой записи файла задается 1 точка кривой. Формат зависит от типа кривой:

1.
тип кривой - 1 - координата Y (значения координат точек по оси X вычисляются по шагу);
2.
тип кривой - 2 - координаты точек по оси X и Y.

Число записей с координатами точек равно числу точек графика.

Расположение данных нескольких графиков

Данные для нескольких графиков располагаются в файле последовательно (для каждого графика задаются свои данные). В конце всего файла с графиками производится запись символа ‘#’.

4.5 Анимация результатов расчетов

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

4.5.1 Входные данные программы.

Входные данные для аниматора вводятся из единого файла, имя которого запрашивается на терминале. При работе программы crmodel для шага МОДЕЛЬ входной файл для аниматора и все необходимые для него файлы (кроме файлов с координатами тел и значениями сил) образуются автоматически. Входной файл в этом случае имеет расширение anc. Управление картинкой при просмотре осуществляется с помощью "горячих"клавиш. Последовательность расположения данных для аниматора во входном файле приводится на нижеследующем примере:

4 количество тел в системе,
sdvigd.dat файл с начальным взаимным расположением тел,
1 признак наличия фона (неподвижного тела), 0 -нет фона,
5 число показываемых сил ( 0 - нет сил ),
dsil.dat файл с описанием изображаемых сил,
sil.dat файл со значениями сил,
stanina.datфайл с геометрией фона ( фон считается телом 1),
geod4.dat файл с геометрией второго тела,
mov4.dat файл с координатами второго тела,
geod1.dat файл с геометрией третьего тела,
mov1.dat файл с координатами третьего тела,
geod3.dat файл с геометрией четвертого тела,
mov3.dat файл с координатами четвертого тела.

При визуализации одного тела файл с начальным расположением тел не вводится. При отсутствии изображения сил не вводятся файлы с описанием и значениями сил. Если какое-то тело необходимо исключить из показа, то перед именами файлов с геометрией и координатами этого тела ставится знак ’-’ (минус). Если для данной модели требуется автоматический запуск ”кино”, то имя файла ”кино” указывается последним. В качестве признака конца данных в файле для аниматора используется в случае необходимости символ ’&’.

4.5.2 Управление картинкой при визуализации.

При анимации можно производить следующие манипуляции с изображением:

клавиша

операция,

x/X

поворот картинки вокруг оси x вперед/назад,

y/Y

поворот картинки вокруг оси y вперед/назад,

z/Z

поворот картинки вокруг оси z вперед/назад,

l/L

уменьшение/увеличение масштаба линейных перемещений,

u/U

уменьшение/увеличение масштаба угловых перемещений,

v/V

уменьшение/увеличение скорости движения,

s

остановка движения тел,

b/B

возобновление движения тел,

p/P

сброс картинки в файл для печати матричным принтером,

^z

завершение работы,

a/A

добавление шлейфа в изображение,

d/D

снятие шлейфа в изображении,

?

смещение картинки вверх,

?

смещение картинки влево,

?

смещение картинки вправо,

?

смещение картинки вниз,

r

обращение направления прокрутки времени,

w/W

уменьшение/увеличение изображения,

e/E

уменьшение/увеличение масштаба векторов сил,

f/F

выкл/вкл изображения векторов сил,

c

смена цвета фона экрана,

C

установка режима сохранения следа,

O/o

вкл/выкл режима раскраски тел различными цветами,

K/k

вкл/выкл режима записи кино,

G/g

вкл/выкл режима демонстрации кино,

t/T

центровка изображение,

n/N

переход к пред/след форме колебаний,

^f

измерение изображаемых сил,

M/m

перейти к следующей/предыдущей форме колебаний при показе перемещений упругого тела только от одной формы,

^m

выключить режим показа перемещений упругого тела от одной формы,

Q/q

перейти к следующей/предыдущей форме колебаний при показе форм колебаний упругого тела,

^q

выключить режим показа форм колебаний упругого тела,

i/I

перейти/выйти в систему координат тела,

S

увеличить фрагмент изображения,

T

измерить расстояние в плоскости экрана,

D

измерить расстояние между точками,

B

выдать координаты точки,

A

измерить углы направляющих косинусов прямой,

R

измерить угол развала колеса,

P

измерить углы ориентации оси поворота колеса,

1

подготовить контрольные точки для вывода,

^l/Alt l

увеличить/уменьшить масштаб показа упругих перемещений тел,

^h

выкл/вкл показ комментариев,

^i

выкл/вкл показ комментариев.

4.5.3 Формат представления координат тел.

Координаты положения тела задаются шестью величинами - перемещениями центра масс по трем неподвижным осям, и тремя углами поворота в подвижной системе последовательно относительно осей X, Y и Z. Координаты положения каждого тела записываются в отдельном файле в свободном формате следующим образом: первая строка - имя задачи - char *40, вторая строка - имя тела - char *40. Третья строка содержит вспомогательную информацию - признак варианта по частоте - отрицательное вещественное число, значение переменной варианта (частота), шаг приращения переменной варианта, четыре вещественных нуля(для твердого тела). После строки вспомогательной информации идут строки с координатами - первое число - время, затем шесть координат. Если по какой-то координате нет движения, то на ее место пишутся нули. При больших пространственных движениях тела после строки с координатами тела идет строка из 9 чисел с матрицей поворота тела. В одном файле может быть любое число вариантов по частоте.

Файлы описанного формата готовятся системой моделирования ФРУНД при обработке результатов решения по коду 19 и 49, классу вывода 1. Пример фрагмента файла с координатами твердого тела:

 issledovanie  
   ЦЕНТР МАСС КРИВОШИПА 1  
  -.30000E+01 .250E+01 .500E-01 .000E+00 .000E+00 .000E+00 .000E+00  
   .00000E+00 .000E+00 .000E+00 .000E+00 .000E+00 .000E+00 .000E+00  
   .25000E-01 .301E-05 .000E+00 .174E-02 .000E+00 .347E-02 .000E+00  
   .40000E-01 .198E-04 .000E+00 .444E-02 .000E+00 .889E-02 .000E+00  
   .60000E-01 .100E-03 .000E+00 .100E-01 .000E+00 .200E-01 .000E+00  
   .80000E-01 .316E-03 .000E+00 .178E-01 .000E+00 .355E-01 .000E+00  
   ............................................................................

В данном файле заданы координаты движения тела для плоского случая.

4.5.4 Формат представления геометрии тел.

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

Файл с геометрией имеет следующую структуру:

1.
число узлов в геометрии;
2.
описание координат узлов;
3.
описание элементов геометрии.

Число узлов - целое число, стоит в первой строке файла геометрии.Описание координат узлов - тройки чисел, задающие координаты узла в пространстве. Номер узла равен порядковому номеру тройки координат. Число троек координат должно строго равняться заданному в первой строке числу узлов. Описание элементов геометрии. Элементы геометрии задаются после описания координат узлов. Общий формат описания - тип геометрического элемента, параметры элемента. Ниже представлены существующие элементы геометрии.

Отрезок

Тип элемента - 2. Параметры:

1.
номер первого узла,
2.
номер второго узла,
3.
цвет отрезка (1 - 15).

Пример строки описания отрезка:

2      2     3      4

Окружность

Тип элемента - 4. Параметры:

1.
номер узла, в котором расположен центр окружности,
2.
номер узла на оси окружности,
3.
радиус окружности,
4.
количество отрезков ломанной кривой, изображающей окружность. Рекомендуемое значение - 8 - 20,
5.
начальный угол рисования окружности в градусах - 0 - 360,
6.
конечный угол рисования окружности в градусах - 0 - 360,

Пример строки описания окружности:

4    2    3   0.4   18    0.0    360

Центр масс

Тип элемента - 100. Параметры: номер узла геометрии, в котором расположен центр масс данного тела. Считается, что тело вращается относительно этого узла.

Пример строки описания центра масс:

100           2

Признак конца задания геометрии

Тип элемента - 300 или 999. Параметров нет. Пример строки описания конца задания геометрии:

999

Текст

Тип элемента - 6. Параметры:

1.
номер узла к которому привязан текст,
2.
длина текста в символах,
3.
вспомогательное вещественное число,
4.
высота шрифта текста - вещественное число,
5.
цвет текста - целое число,
6.
строка текста - задается в отдельной строке, следующей за первыми пятью параметрами.

Для аниматора существенными являются первый и шестой параметры.

Пример строки описания конца задания геометрии:

6      3     13     1.      3.     4  
строка_текста

Файл с собственными частотами и формами

Тип элемента - 20. Параметры:

1.
имя файла с собственными формами и частотами. Считается, что указанный файл содержит частоты и формы в формате, который используется при вводе форм в библиотеку тел,
2.
количество собственных частот и форм в файле,
3.
нормировочный множитель для коэффициентов форм по поступательным степеням свободы узла.Данная информация используется при анимации движения упругих тел. Если формы колебаний вычислены не для всех узлов геометрии тела, то при анимации смещения таких узлов будут равны нулю.

Пример строки описания конца задания геометрии:

20      rama.txt     20     0.0319

Пример задания геометрии тела
        3  
 0.0    0.0    0.0  
-0.2    0.0    0.0  
0.2     0.0    0.0  
2    2    3    6  
3    2    2   0.4   18    0.0    360.  
3    2    3   0.4   18    0.0    360.  
100   1  
20   RAMA.TXT  20  0.0319  
300

В данном примере описаны две окружности и отрезок, соединяющий их центры.

4.5.5 Формат представления сил, действующих на точки тел

Все силы, действующие на точки тел и которые необходимо визуализировать, записываются в одном файле, в одной строке для каждого момента времени. Первая строка такого файла содержит вспомогательную информацию аналогично п.п. 4.5.3. После строки вспомогательной информации идут строки со значениями сил для каждого момента времени. Строки со значениями сил должны быть записаны для тех же моментов времени что и координаты тел. В одном файле сил может быть любое число вариантов по скорости (частоте). Файлы такого формата готовятся системой моделирования ФРУНД при обработке результатов решения по коду 20 и классу вывода 3. После работы шага МОДЕЛЬ образуются файлы ANIMASIL.UPR с контрольными точками по коду 20. В этом файле все строки закомментированы символом ’$’ в первой позиции. Для показа нужных сил необходимо убрать комментарии с соответствующих контрольных точек в файле ANIMASIL.UPR. При обработке результатов решения по этому коду запрашивается имя файла для записи значений сил, и образуется файл с указанием точек приложения сил. Максимальное количество одновременно показываемых в аниматоре сил - 20. Пример фрагмента файла со значениями сил.

  -.30000E+01  .250E+01  .500E-01  .000E+00  .000E+00  .000E+00  .000E+00...  
   .00000E+00  .222E+02  .111E+02  .000E+00  .113E-14  
   .75228E+00  .222E+02  .111E+02  .598E+00 -.904E-03  
   ......................................................  
   .76708E+01  .221E+02  .113E+02  .609E+01 -.944E-01  
   .13104E+02  .219E+02  .118E+02  .104E+02 -.277E+00

В данном файле записаны значения для пяти сил. Длина реализаций сил может быть любая.

4.5.6 Другая дополнительная информация для аниматора

Начальное взаимное расположение тел

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

Описание изображаемых сил.

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

1.
номером тела, куда она приложена,
2.
номером узла тела,
3.
направлением действия - 1 - по оси x, 2 - по оси y, 3 - по оси z.

Номер тела может быть указан со знаком минус, что означает, что данная сила будет изображаться в системе координат связанной с телом, при положительном номере тела сила изображается в неподвижной системе координат. Порядок описания изображаемых сил соответствует порядку задания сил в файле представления значений сил (п.п. 4.5.5).

4.5.7 Режим работы ”кино”

Аниматор дает возможность запоминания произвольной последовательности команд и последующего их воспроизведения, что позволяет создавать компьютерные фильмы. Режим работы ”кино” содержит две группы команд - запись кино - K/k и демонстрации кино - G/g. При включенном режиме записи кино все команды с клавиатуры и временные задержки сбрасываются в файл с именем kino.dat. После записи кино и выхода из программы аниматора, при необходимости следует переименовать этот файл, иначе после следующего сеанса работы программы он может быть изменен. Для автоматического запуска кино в конце файла данных для аниматора следует поставить имя файла кино. Клавиша G вызывает файл кино с именем kino.dat. В режиме демонстрации кино блокируются все клавиши кроме HELPа и выхода из режима демонстрации - g. Файл кино допускает редактирование для вставки комментариев и объединения файлов кино разных моделей. Структура файла кино - строки управляющей информации и строки команд. Строки управляющей информации начинаются с отрицательного целого числа -1,-2,-3,-4. Строки с кодами -1,-2 образуются автоматически. Строки с кодами -3,-4 служат для вставки комментариев и новых моделей. Код -3 задает новую модель, после кода строки идет имя файла для аниматора. После этой строки содержимое файла кино для данной модели копируется в сборный файл. Код -4 задает комментарии, которые появляются на экране во время демонстрации фильма. После кода строки задается количество строк комментария и в последующих строках сам текст комментария.Строки с кодами -3,-4 не могут стоять между строками с кодами -1,-2.

4.5.8 Изображение соединительных элементов при анимации

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

Файл с расширением sel состоит из закомментированных пар строк. Каждая пара относится к одному соединительному элементу. Первая строка в паре - имя соединительного элемента, служит для его указания. Вторая строка в паре содержит описание положения и схемы изображения.

Во второй строке пары последовательно перечисляются:

1.
номер тела 1 точки крепления;
2.
номер узла тела 1 точки крепления;
3.
номер тела 2 точки крепления;
4.
номер узла тела 2 точки крепления;
5.
тип схемы - целое число: 1 - винтовая пружина, 2 - мышца, 3 - амортизатор;
6.
первый параметр схемы, вещественное число;
7.
второй параметр схемы, вещественное число.

В зависимости от типа схемы изображения соединительного элемента параметры имеют различный смысл.

Винтовая пружина - первый параметр - диаметр пружины, второй - число витков.

Мышца - первый параметр - максимальный поперечный размер, второй свободный.

Амортизатор - первый параметр - минимальная длина амортизатора, второй - максимальная длина.

По умолчанию, при генерации файла ********.sel, все элементы показываются в виде пружин диаметром 55 мм с четырьмя витками.

Для показа соединительных элементов при анимации в файле *.sel нужно раскомментировать вторые строки пар нужных соединительных элементов и модифицировать, при необходимости, параметры этих строк.

Глава 5
ВСПОМОГАТЕЛЬНЫЕ ПОДПРОГРАММЫ - ПОДГОТОВКА МИКРОПРОФИЛЯ И Т.Д.

5.1 Подготовка микропрофиля

Для расчетов динамики движения транспортных средств по неровностям дороги используются записи реализаций случайного микропрофиля, которые могут быть получены экспериментально, например, с помощью геодезической съемки, или сгенерированы с помощью численных алгоритмов. В состав системы входит программа подготовки микропрофиля, работающая с файлами экспериментальных ординат неровностей дорог трех типов автополигона НАМИ. Указанные экспериментальные файлы лежат в системной директории ФРУНДа - ../UTILF/DOROGI, а именно: bulglad - ровный булыжник; bwb1 - булыжник с выбоинами и буграми; dmd - динамометрическая дорога.

Подготовка микропрофиля производится специальной программой - obpron, которая центрирует экспериментальные ординаты высот неровностей и записывает их в необходимом виде в файл, который используется затем на шаге РАСЧЕТ. При работе из оболочки - пункт меню ‘Вспомогательные функции’, подпункт ‘Подготовка автомобильного микропрофиля’. Программа obpron работает с двумя файлами, имена которых запрашиваются с терминала. Первый файл содержит экспериментальные ординаты микропрофиля, записанные в свободном формате по 2 ординаты в строке (левая и правая колея). Первые 8 строк данного файла содержат обязательные комментарии. В качестве примера можно обратиться к поставляемым файлам bulglad, bwb1, dmd. Второй файл содержит результаты предварительной обработки микропрофиля и подается по запросу программы rashet32. Программа запрашивает также некоторые вспомогательные данные, смысл которых определяется в запросах. Более точное описание сглаживающей способности шины получается при подготовке микропрофиля под заданную скорость движения, на которой будет производиться расчет. Для пространственной модели готовятся файлы с микропрофилем, различные для каждой колеи. Для плоской модели на каждую колею можно подавать одинаковые файлы.

5.2 Нормирование собственных форм

Если собственные формы упругих подсистем были рассчитаны в системе единиц отличной от той в которой в дальнейшем предполагается вести расчеты, то перед генерацией модели (программа model32), нужно произвести нормировку форм с учетом выбранной системы единиц. Нормировка форм производится программой umnform - пункт меню оболочки ‘Нормировка коэффициентов форм упругих тел’. Указанная программа за один вызов производит нормировку форм для одной упругой подсистемы, номер которой она запрашивает с терминала. Программа использует файл с постоянным именем mnogit, в котором перечислены нормировочные множители последовательно для всех форм. В первой строке файла mnogit указывается число форм в теле (может быть больше реального количества форм), затем в каждой строке два числа - коэффициент нормировки для перемещений и для углов, как правило эти строки одинаковые.

5.3 Связь с конечно-элементными комплексами и передача форм колебаний упругих тел

При моделировании упругих тел используются конечно - элементные модели таких тел, построенные с помощью существующих программ расчета МКЭ. При этом нужно учитывать два момента - нумерация узлов в КЭМ может не совпадать с нумерацией узлов в геометрии ФРУНДа и формат выдачи форм колебаний не соответствует входному формату ФРУНДа. Поэтому для связи с программами КЭМ могут использоваться специальные программы для согласования геометрии и формата передачи собственных форм. Такие программы разработаны для взаимодействия с МКЭ комплексом NASTRAN и BASYS. Рассмотрим эти программы.

5.3.1 Взаимодействие с комплексом NASTRAN

Для перезаписи собственных частот и форм из формата NASTRAN в формат ФРУНД разработана программа nasf, которая осуществляет требуемое переформатирование. Для работы программы nasf требуются следующие файлы:

1.
файл описания модели NASTRANа (расширение *.bdf);
2.
файл с частотами и формами колебаний, рассчитанными в NASTRANе (расширение *.out);
3.
файл геометрии в формате ФРУНДа (расширение *.geo).

Программа nasf запрашивает перечисленные файлы и записывает собственные частоты и формы в формате ФРУНДа. Имя для выходного файла можно давать произвольное, рекомендуемое расширение - txt.

Программа nasf работает следующим образом. Из файла *.bdf считываются координаты узлов модели NASTRANа, а из файла *.geo считываются координаты узлов тела в модели ФРУНДа. Затем программа для каждого узла тела из геометрии ФРУНДа находит ближайший или совпадающий узел в геометрии NASTRANa. Далее записывает коэффициенты форм для ближайших узлов в выходной файл, сохраняя при этом нумерацию узлов геометрии ФРУНДа.

Геометрия для ФРУНДа может быть получена из геометрии конечно-элементной модели. Так как для расчета динамики во ФРУНДе не требуется столь точная геометрия, которая используется в NASTRANe, то необходимые узлы и контуры тела должны быть заданы в NASTRANе с помощью фиктивных стержневых элементов типа CROD. Для получения стержневой геометрии из NASTRANа в формате ФРУНДа разработана программа nrods. Для работы программы nrods требуется файл описания модели NASTRANа (расширение *.bdf) в котором обязательно описана упрощенная геометрия с помощью фиктивных стержневых элементов. После работы nrods образуется файл в формате *.geo, соответствующий геометрии фиктивных стержневых элементов.

При использовании в моделях упругих тел, как правило складывается ситуация, когда первоначально строится полностью твердотельная модель, а потом то или иное тело заменяется на упругое. В этом случае геометрия, полученная из NASTRANа с помощью программы nrods, будет обладать другой нумерацией узлов и, скорее всего, не точно совпадать со «старой» геометрией ФРУНДа. Поэтому для использования геометрии, полученной из NASTRANа, ее следует перенумеровать.

Для перенумерации геометрии служит программа renuzl. Эта программа в качестве исходных данных использует файл с геометрией тела в твердотельной модели (старую геометрию) и файл с геометрий тела, полученный из NASTRANа с помощью программы nrods (новую геометрию). Renuzl ищет в новой геометрии узлы, совпадающие (с заданной точностью)с узлами старой геометрии, и присваивает этим узлам номера старой геометрии.

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

Поскольку комплекс NASTRAN может быть установлен не только на IBM совместимом компьютере, то тексты источников вышеперечисленных программ (язык C) приложены к системе ФРУНД и находятся в директории ../UTILF/NASTRAN.

5.3.2 Взаимодействие с комплексом BASYS

basgeom - программа готовит по входному файлу системы BASYS файл геометрии в формате ФРУНДа (файл *.geo), а также файл соответствия номеров узлов в КЭМ и файле геометрии для ФРУНДа. При этом учитывается, что нумерация узлов во ФРУНДе сплошная, начиная с 1. Поэтому нумерация КЭМ изменяется в соответствии с этим а изменения записываются в файл соответствия номеров узлов.

renbas - программа готовит файл собственных форм с учетом новой нумерации узлов в геометрии ФРУНДа. Входными данными для нее являются файл собственных частот и форм, полученный из BASYSа, а также файл соответствия номеров узлов в КЭМ и ФРУНДе. После работы программы образуется файл собственных частот, который можно вводить в библиотеку ФРУНДа.

5.3.3 Формат файла с собственными частотами и формами

Собственные частоты и формы колебаний упругого тела могут получаться с помощью любой МКЭ программы, но для ввода в библиотеку тел они должны быть записаны в нужном формате. В файле для каждой частоты указывается ее номер, значение частоты, номера узлов и коэффициенты форм по шести направлениям. Такой формат проиллюстрируем примером:

 nomer ~astoty    1   ~astota       0.00  
 
 
 1     1 +1.9437e+00 -8.7443e-02 -7.9858e-01 -6.0283e-05 +0.0000e+00 -1.6815e-05  
 1     2 +1.9168e+00 -8.3994e-02 -8.0930e-01 -6.0283e-05 -5.2256e-04 -1.6815e-05  
 1     3 +1.9674e+00 -8.8939e-02 -8.3695e-01 -6.0283e-05 -5.2256e-04 -1.6815e-05  
 1     4 +1.9249e+00 -1.1773e-01 +2.1016e-01 -6.0285e-05 -5.2256e-04 -1.6815e-05  
 1     5 +1.9070e+00 -1.1773e-01 +2.7406e-01 -6.0283e-05 -5.2256e-04 -1.6815e-05  
 1   228 +1.9663e+00 -8.8939e-02 -8.3312e-01 -6.0283e-05 -5.2256e-04 -1.6815e-05  
 1   229 +1.9851e+00 -1.1120e-01 -2.0860e-01 -6.0283e-05 -5.2256e-04 -1.6815e-05  
 6     0 +0.0000e+00 +0.0000e+00 +0.0000e+00 +0.0000e+00 +0.0000e+00 +0.0000e+00  
 
 nomer ~astoty    2   ~astota       0.00  
 
 
 1     1 +1.3993e-01 +2.1221e+00 +4.1607e-01 +5.3991e-04 +0.0000e+00 -1.9497e-04  
 1     2 +1.5037e-01 +2.0983e+00 +4.2023e-01 +5.3991e-04 +2.0283e-04 -1.9497e-04  
 1     3 +1.3074e-01 +2.1608e+00 +4.3096e-01 +5.3991e-04 +2.0283e-04 -1.9497e-04  
 1     4 +3.2823e-01 +1.6904e+00 +4.8846e-01 +5.3991e-04 +2.0283e-04 -1.9497e-04  
 1     5 +1.2157e-01 +1.6904e+00 -8.3845e-02 +5.3991e-04 +2.0283e-04 -1.9497e-04  
 1   228 +1.1834e-01 +2.1608e+00 +3.9663e-01 +5.3991e-04 +2.0283e-04 -1.9497e-04  
 1   229 +2.0593e-01 +1.9284e+00 +3.9739e-01 +5.3991e-04 +2.0283e-04 -1.9497e-04  
 6     0 +0.0000e+00 +0.0000e+00 +0.0000e+00 +0.0000e+00 +0.0000e+00 +0.0000e+00

5.4 Подготовка данных по рельсовому пути

Моделирование качения железнодорожного колеса по рельсовому пути с помощью многополюсника типа 420 - п.п. 12.4, использует данные по неровности рельсового полотна и профилю колеса и рельса. Подготовка профиля производится с помощью программы railway.exe. Программа использует три файла с постоянными именами:

1.
spcdataz.txt - значения ординат вертикальных неровностей - одинаковые для левой и правой колеи;
2.
spcdatac.txt - разница неровностей левой и правой колей по вертикали;
3.
spcdatay.txt - поперечное смещение левой и правой колей.

В каждом из этих файлов данные хранятся в формате файлов с параметрами микропрофиля п.п. 5.1, только для левой и правой колеи заносятся одинаковые числа. Файлы с такой информацией могут быть получены из экспериментальных данных или случайные реализации в них сгенерированы численными методами. Файл spcdataz.txt должен присутствовать обязательно. При отсутствии spcdatac.txt левая и правая колея считаются одинаковыми. При отсутствии spcdatay.txt, в поперечном направлении профиль считается гладким. Выходной бесформатный файл имеет постоянное имя railway.dat и содержит суммарную информацию о неровностях.

Все реализации микропрофиля колей должны быть сгенерированы с одинаковым шагом, который задается в блоке 37 файла задания параметров расчетной схемы ФРУНДа, аналогично заданию шага дискретизации микропрофиля автомобильных дорог

Подготовка данных по профилям сечений рельса и колеса производится программой spectr.exe. Профили задаются плоскими кривыми с помощью встроенного графического редактора и хранятся в файлах с именами, задаваемыми пользователем. В настоящий момент подготовлены два файла:

1.
wheel.dat - данные по профилю нового колеса;
2.
r50.dat - данные по профилю неизношенного рельса r50 в.

С помощью соответствующей функции программы Spectr.exe подготавливаются текстовые файлы профилей в формате файлов wheelpt.dat, railpt.dat, необходимых для работы многополюсника типа 420. Рекомендуемый шаг дисретизации данных профилей - 0.5 мм.

Глава 6
ТИПЫ ФУНКЦИОНАЛЬНЫХ ЗАВИСИМОСТЕЙ СИЛ В СОЕДИ НИТЕЛЬНЫХ ЭЛЕМЕНТАХ И СИЛ НАГРУЗОК

6.1 Структура сил в соединительных элементах

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

Рисунок 6. 1

На рис. 6.1. изображены две подсистемы 1 и 2, узлы 10 и 11 которых связаны соединительным элементом 20. Если в блоке BL17 первым был описан узел 10, т.е. строка блока имела вид

20    1    10    2    11    1,

то, в этом случае, физическому сжатию соединительного элемента 20 будет соответствовать положительное значение его относительного смещения. При описании первым узла 11, его сжатию будет соответствовать отрицательное значение относительного смещения.

Такое правило справедливо при рассмотрении с. э. в случае малых движений тел и смещений в кинематических парах при больших движениях. Для с. э. типа «пружина в шарнирах» сжатию всегда соответствует положительное значение деформации, независимо от очередности описания узлов. Функция, описывающая каждое слагаемое, зависит от определенного числа параметров, задаваемых пользователем в соответствии с приводимыми ниже описаниями возможных в данной системе видов функций. Пользователь описывает вид функции силы в с. э. в исходных данных шага МOДЕЛЬ в блоке BL16, для каждого слагаемого отдельно. Затем значения параметров этих функций задаются в блоке BL33 шага РАСЧЕТ. Число параметров функции может быть больше числа вводимых параметров (в блоке BL33), поскольку некоторые параметры могут быть зависимыми и вычисляться внутри системы.

6.2 Типы функций, описывающих упругое слагаемое

Тип номер 1 - линейная функция - рис. 6.2. Число параметров - 1. Число вводимых параметров - 1. Вводимые параметры: 1. коэффициент жесткости.

Рисунок 6. 2

Тип номер 2 - кусочно-линейная функция из двух участков, причем первый участок проходит через начало координат - Рис. 6.3. Число параметров - 4. Число вводимых параметров - 3. Вводимые параметры:

1.
коэффициент жесткости на первом участке,
2.
коэффициент жесткости на втором участке,
3.
абсцисса границы участков - всегда положительная величина.

Здесь и далее участки номеруются слева направо.

Рисунок 6. 3

Тип номер 3 - кусочно-линейная функция из трех участков, причем функция на среднем участке проходит через начало координат рис. 6.4. Число параметров - 7. Число вводимых параметров - 5. Вводимые параметры:

1.
коэффициент жесткости на первом участке,
2.
коэффициент жесткости на втором участке,
3.
коэффициент жесткости на третьем участке,
4.
абсцисса границы первого и второго участков - всегда меньше 0,
5.
абсцисса границы второго и третьего участков - всегда больше 0.

Рисунок 6. 4

Тип номер 4 - кусочно-линейная функция из N участков - рис. 6.5. Число параметров - (N+1)*2+1. Число вводимых параметров - (N+1)*2+1.Вводимые параметры:
1. число участков,
2. абсцисса первой точки,
3. ордината первой точки,
........................................,
(n+1)*2. - абсцисса N+1 точки,
((n+1)*2+1). - ордината N+1 точки.

Рисунок 6. 5

Тип номер 5 - симметричная кубическая парабола с линейным членом рис. 6.6. Число параметров - 2. Число вводимых параметров - 2. Вводимые параметры:

1.
коэффициент при линейном члене,
2.
коэффициент при кубическом члене.

Рисунок 6. 6

Тип номер 7 - линейная функция, проходящая через начало координат с участком постоянного значения - рис. 6.7. Число параметров - 3. Число вводимых параметров - 2. Вводимые параметры:

1.
коэффициент жесткости линейного участка,
2.
величина на постоянном участке- всегда < 0.

Рисунок 6. 7

Тип номер 8 - зазор с односторонним линейным упором - рис. 6.8. Число параметров - 3. Число вводимых параметров - 2. Вводимые параметры:

1.
коэффициент жесткости упора,
2.
величина одностороннего зазора положительная величина.

Рисунок 6. 8

Тип номер 10 - кусочно - линейная функция с постоянным участком - рис. 6.9. Число параметров - 6. Число вводимых параметров - 4. Вводимые параметры:

1.
коэффициент жесткости первого участка,
2.
коэффициент жесткости второго участка,
3.
величина постоянной силы - отрицательная величина,
4.
абсцисса границы второго и треть его участков - всегда > 0.

Рисунок 6. 9

Тип номер 11 - кусочно - линейная функция с постоянным участком - рис. 6.10. Число параметров - 6. Число вводимых параметров - 4. Вводимые параметры:

1.
коэффициент жесткости первого участка,
2.
коэффициент жесткости второго участка,
3.
величина постоянной силы - отрицательная величина,
4.
абсцисса границы второго и третьего участков - всегда меньше нуля.

Рисунок 6. 10

Тип номер 13 - линейная функция с постоянным смещением - рис. 6.11. Число параметров - 2. Число вводимых параметров - 2. Вводимые параметры:

1.
коэффициент жесткости,
2.
величина смещения.

Рисунок 6. 11

Тип номер 30 - модель рессоры НАМИ - рис. 6.12. Число параметров - 20. Число вводимых параметров - 20. Вводимые параметры:

1.
величина жесткости на первом участке ветви разгрузки - нижней ветви (участки нумеруются и задаются слева направо),.
2.
величина жесткости на втором участке ветви разгрузки,
3.
величина жесткости на третьем участке ветви разгрузки,
4.
величины жесткости на первом участке ветви нагрузки - верхней ветви,
5.
величины жесткости на втором участке ветви нагрузки,
6.
величины жесткости на третьем участке ветви нагрузки,
7.
любое число,
8.
величина постоянного члена ветви разгрузки на втором участке,
9.
любое число,
10.
любое число,
11.
величина постоянного члена ветви нагрузки на втором участке,
12.
любое число,
13.
начало первого участка характеристики - отрицательная величина,
14.
начало второго участка характеристики - всегда равно нулю,
15.
начало третьего участка характеристики - положительная величина (включение подрессорника),
16.
показатель первой экспоненты,
17.
показатель второй экспоненты,
18.
весовой коэффициент первой экспоненты,
19.
весовой коэффициент второй экспоненты (сумма весовых коэффициентов равна единице),
20.
ширина зоны нечувствительности рессоры - не более 0.1 мм.

Рисунок 6. 12

Тип номер 31 - линейный матричный упругий элемент с матрицей 3*3. Число параметров - 11. Число вводимых параметров - 11. Вводимые параметры:

1-9. значения элементов матрицы жесткости, задаваемые последовательно по столбцам,
10. угол поворота элемента вокруг оси X в градусах,
11. угол поворота элемента вокруг оси Y в градусах.

Матричный упругий элемент может применяться только для случая наличия трех составляющих сил в механической характеристике - по осям X,Y,Z, хотя реальная жесткость может присутствовать только по одной оси. Матрица жесткости задается в неподвижной системе координат, а реальная ее ориентация определяется двумя углами.

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

6.3 Типы функций, описывающих демпфирующее слагаемое

Тип номер 50 - линейная функция - линейное вязкое демпфирование - рис. 6.13. Число параметров - 1. Число вводимых параметров - 1. Вводимые параметры: 1. коэффициент демпфирования.

Рисунок 6. 13

Тип номер 51 - линейное вязкое демпфирование с сухим трением - рис. 6.14. Число параметров - 2. Число вводимых параметров - 2. Вводимые параметры:

1.
коэффициент вязкого демпфирования,
2.
величина силы сухого трения.

Рисунок 6. 14

Тип номер 52 - кусочно-линейная функция из N участков. Число параметров - (N+1)*2+1. Число вводимых параметров - (N+1)*2+1. Вводимые параметры:
1. число участков (N),
2. абсцисса первой точки,
3. ордината первой точки,
........................................,
(n+1)*2. - абсцисса N+1 точки,
((n+1)*2+1). - ордината N+1 точки.

Тип номер 53 - сухое трение с начальным вязким участком - рис. 6.15. Число параметров - 3. Число вводимых параметров - 2. Вводимые параметры:

1.
величина силы сухого трения,
2.
величина скорости, до которой сопротивление вязкое.

Рисунок 6. 15

Тип номер 54 - кусочно - линейная функция из двух участков, проходящая через ноль - рис. 6.16.Число параметров - 2. Число вводимых параметров - 2. Вводимые параметры:

1.
коэффициент вязкого демпфирования при отрицательном аргументе,
2.
коэффициент вязкого демпфирования при положительном аргументе.

Рисунок 6. 16

Тип номер 55 - кусочно - линейная функция из четырех участков, проходящая через начало координат (автомобильный амортизатор) - рис.6.17. Число параметров - 8. Число вводимых параметров - 6. Вводимые параметры:

1.
коэффициент наклона прямой на первом участке, любое число,
2.
коэффициент вязкого демпфирования на втором участке, положительное число,
3.
коэффициент вязкого демпфирования на третьем участке, положительное число,
4.
коэффициент наклона прямой на четвертом участке, любое число,
5.
абсцисса границы первого и второго участков, отрицательное число,
6.
абсцисса границы третьего и четвертого участков, положительное число.

Рисунок 6. 17

Тип номер 56 - сухое трение с начальным вязким и ниспадающим основным участком - рис. 6.18.Число параметров - 5. Число вводимых параметров - 3. Вводимые параметры:

1.
величина силы сухого трения,
2.
коэффициент жесткости ниспадающего участка,
3.
величина скорости ,до которой сопротивление вязкое.

Рисунок 6. 18

Тип номер 58 - вязкое демпфирование, включающееся при ударе в односторонний ограничитель. Число параметров - 2. Число вводимых параметров -2. Вводимые параметры:

1.
коэффициент вязкого сопротивления,
2.
величина хода, при которой включается демпфирование.

Данный тип является функцией двух переменных - относительной скорости и относительного смещения.

Тип номер 61 - матричный демпфирующий элемент с матрицей 3*3. Число параметров - 11. Число вводимых параметров - 11. Вводимые параметры:
1 - 9.значения элементов матрицы демпфирования. Задаются последовательно по столбцам.
10.угол поворота элемента вокруг оси X в градусах.
11.угол поворота элемента вокруг оси Y в градусах.

Замечания - см. замечания к типу 31 упругой составляющей.

6.4 Структура сил нагрузок (специальных сил)

Специальная сила (с. с.) -это сила, зависящая от времени. Пользователь указывает вид функциональной зависимости с. с. в исходных данных шага МОДЕЛЬ. Различаются следующие типы функций с. с.:

1.
Спец. силы с аргументом T. Номера типов функций 1 - 30.
2.
Спец. силы с аргументом V,T. Номера типов функций 31-40 и 201 - 209.

Указанные функции применяются для моделирования сил, вызванных дисбалансом вращающихся масс.Функции, не требующие аргументов, могут относиться к любому типу.

6.5 Типы функций спец. сил

Тип номер 1 - периодическое возбуждение. Число параметров - 2. Число вводимых параметров - 2. Вводимые параметры:

1.
амплитуда силы,
2.
фазовый сдвиг силы.

Тип номер 2 - постоянная сила. Число параметров - 1. Число вводимых параметров - 1. Вводимые параметры: величина силы.

Тип номер 3 - сила веса. Число параметров - 1. Число вводимых параметров - 1. Данный вид с.с. применяется для моделирования весовой нагрузки. Величина силы вычисляется автоматически по массе тела. Сила имеет отрицательный знак. Может действовать только на твердое тело по оси Z.

Тип номер 4 - периодическое импульсное воздействие синусообразной формы. Число параметров - 5. Число вводимых параметров - 4. Вводимые параметры:

1.
период силы,
2.
фазовый сдвиг силы,
3.
продолжительность импульса,
4.
амплитуда импульса силы.

Тип номер 12 - периодическое воздействие с переменной частотой. Число параметров - 4. Число вводимых параметров - 4. Вводимые параметры:

1.
амплитуда силы,
2.
начальная частота,
3.
коэффициент приращения частоты,
4.
фазовый сдвиг силы.

Тип номер 13 - Прямоугольный импульс. Число параметров - 3. Число вводимых параметров - 3. Вводимые параметры:

1.
величина импульса,
2.
момент начала действия,
3.
момент окончания действия.

Тип номер 31 - дисбаланс вращающихся масс. Число параметров - 3. Число вводимых параметров - 3. Вводимые параметры:

1.
постоянная, связывающая скорость движения с угловой скоростью вращающихся масс, находится из выражения C1=OMEGA*V,
2.
величина дисбаланс вращающихся масс,
3.
фазовый сдвиг силы неуравновешенности.

Тип номер 32 - периодическое возмущение, с амплитудой, пропорциональной квадрату частоты. Число параметров - 2. Число вводимых параметров - 2. Вводимые параметры:

1.
амплитуда силы (дисбаланс),
2.
фазовый сдвиг силы.

Примечание. Данный тип функции предназначен для моделирования сил от дисбаланса при коде расчета 4 (резонансные испытания). Тип номер 201 – скачок с косинусоидальными участками возрастания и убывания. Число параметров - 5. Число вводимых параметров - 5. Вводимые параметры:

1.
Время начала действия скачка,
2.
время выхода на постоянный участок силы,
3.
время начала убывания силы,
4.
время окончания действия слы,
5.
величина скачка.

Тип номер 202 – табличная периодическая функция. Число параметров - 3. Число вводимых параметров - 3. Вводимые параметры:

1.
Номер функции – 1 - 5,
2.
частота функции (Гц),
3.
фаза (радианы).

Данный тип функции использует таблицу значений ординат, хранящихся в файле с именем ss202N, где N – номер функции. Например имя файла для функции с номером 3 будет – ss2023. В файле с ординатами перечислены значения функции для равномерного шага за один период.

6.6 Свободные типы функций

Система моделирования ФРУНД позволяет пользователю создавать свои произвольные типы зависимостей сил в податливых связях и нагрузках и использовать их в моделях ФРУНДа. Такая возможность реализуется путем написания пользователем своей подпрограммы вычисления силы и стандарта подсоединения такой подпрограммы в систему моделирования.

6.6.1 Свободные типы сил в соединительных элементах

Свободно программируемые типы сил в соединительных элементах имеют зарезервированные номера типов начиная с 60. Такой номер типа относится к упругому слагаемому в механической характеристики. При описании механической характеристики в блоке 16, упругому слагаемому нужной составляющей присваивается свободный, выбранный пользователем тип, например 60. Количество параметров упругого слагаемого и их смысл определяется пользователем. Демпфирующее слагаемое в таком случае можно указывать любое (например линейное), в модели оно будет игнорироваться. Например:

BL16 NMX PRSOST TUF KPAR TDF KPAR  
@ ’гидравлический буфер’  
             1           3         60      5       50        1

Для вычисления значения силы пользователю необходимо написать подпрограмму на языке ФОРТРАН со стандартным именем и входом. Имя подпрограммы состоит из шести символов, первые три - FUN, последние три - номер типа функции. Так имя подпрограммы типа 60 имеет вид FUN060. Пользователь пишет подпрограмму :

SUBROUTINE FUN060(X,Y,XO,YO,P,FUN)  
REAL P(1)  
........................................  
тело программы  
........................................  
XO=............  
YO=............  
FUN=XO+YO  
RETURN  
END

Здесь X и Y - входные переменные - относительная деформация и скорость относительной деформации соединительного элемента. XO,YO,FUN - выходные переменные, значение упругой, демпфирующей и суммарной силы в связи. Значения XO и YO передаются для обеспечения дальнейшего вывода сил. В элементах массива P(1),P(2), ....., передаются параметры подпрограммы, смысл и значения которых определяются пользователем. Сами значения параметров задаются обычным путем - в блоке 33.

6.6.2 Свободные типы сил нагрузок

Свободно программируемые типы специальных сил (нагрузок) имеют зарезервированные номера типов начиная с 210. При описании спец. силы в блоке 18, нужной составляющей присваивается свободный, выбранный пользователем тип, например 210. Количество параметров и их смысл определяется пользователем. Например:

BL18  NSS   RPSOST   TFUN   KPAR  
@    ’импульсная сила’  
       1      3      210      2

Для вычисления значения силы пользователю необходимо написать подпрограмму на языке ФОРТРАН со стандартным именем и входом. Имя подпрограммы состоит из шести символов, первые три - FUN, последние три - номер типа функции. Так имя подпрограммы типа 210 имеет вид FUN210. Пользователь пишет подпрограмму:

SUBROUTINE FUN210(T,V,P,FUN)  
REAL P(1)  
........................................  
тело программы  
........................................  
FUN=........  
RETURN  
END

Здесь T и V - входные переменные - время и скорость движения, FUN - выходная переменная - значение силы. В элементах массива P(1),P(2), ....., передаются параметры подпрограммы, смысл и значения которых определяются пользователем. Сами значения параметров задаются обычным путем - в блоке 34.

6.6.3 Включение подпрограмм свободных типов в систему

Модели, содержащие свободные типы функций, должны при сборке исполнимого модуля подключать объектные модули пользовательских функций. Это достигается путем внесения имен файлов объектных модулей в сборочный файл RASHET32.BLD.

Для подключения, например, функции типа 60, в конец этого файла добавляется строка с именем объектного модуля: fun060.obj

Глава 7
ОПИСАНИЕ МОДЕЛЕЙ СИСТЕМ ТЕЛ С КИНЕМАТИЧЕСКИМИ ПАРАМИ

7.1 Задание кинематических пар

Кинематическая пара - кинематическая связь задается соединительным элементом как обычная податливая связь. Указание на то, что связь является кинематической, производится в описании механической характеристике. Для этого служит специальный тип функции упругого слагаемого составляющей в механической характеристики - тип 15. Возможное количество параметров - 1,2,3,4, смысл которых раскрывается в главе 9. В случае задания у слагаемого силы такой упругой составляющей демпфирующее слагаемое игнорируется. При использовании кинематических пар параметры механической характеристики должны быть равны 0, если специально не моделируется податливая пара.

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

7.1.1 Одномерный случай

Жесткую связь узлов 3 и 4 тел 1 и 2 - рис. 7.1 можно описать механической характеристикой.

BL16  NMX PRSOST TUF KPAR TDF KPAR  
   1           1        15       1       50       1

Считая, что каждое тело имеет только одну степень свободы - перемещение вдоль X, такая связь обеспечит одинаковые перемещения тел 1 и 2.

Рисунок 7. 1

7.1.2 Плоский случай

Если каждое тело на рис. 7.1 обладает тремя степенями свободы (совершает плоское движение), то связь 3 и 4 узлов с помощью следующей механической характеристики

BL16  NMX PRSOST TUF KPAR TDF KPAR  
  1           1         15       1      50        1  
  1           2         15       1      50        1

обеспечивает вращательное движение одного тела относительно другого, т.е. моделируется плоский цилиндрический шарнир. На рис. 7.2 показан случай поступательной кинематической пары. Тело 1 может перемещаться вдоль поверхности тела 2 без отрыва по направлению X. Задание такой связи возможно с помощью двух соединительных элементов со следующей механической характеристикой:

BL16  NMX PRSOST TUF KPAR TDF KPAR  
   1          2         15       1       50       1

Описание соединительных элементов выглядит при этом следующим образом

BL17  NSE  NT1  NU1  NT2  NU2  NMX1  
   1        1        3        2      10        1  
   2        1        4        2      10        1

Рисунок 7. 2

7.1.3 Пространственный случай

Сферический шарнир в пространственном случае - рис. 7.3 моделируется с помощью трехмерной механической характеристики

BL16  NMX PRSOST TUF KPAR TDF KPAR  
   1          1         15       1        50      1  
   1          2         15       1        50      1  
   1          3         15       1        50      1

Рисунок 7. 3

Цилиндрический шарнир в пространственном случае - рис. 7.4 моделируется двумя соединительными элементами с разными механическими характеристиками

Рисунок 7. 4

BL16  NMX PRSOST TUF KPAR TDF KPAR  
   1          1          15      1       50      1  
   1          2          15      1       50      1  
   1          3          15      1       50      1  
   2          1          15      1       50      1  
   2          2          15      1       50      1

На рис. 7.4 приняты одинаковые номера узлов для обоих тел в точках цилиндрического шарнира. Соединительные элементы будут задаваться:

BL17  NSE  NT1  NU1  NT2  NU2  NMX1  
   1        1        3        2       3         1  
   2        1        4        2       40       2

На аналогичном принципе может быть построено моделирование и поступательных пар в пространственном случае.

Глава 8
ОПИСАНИЕ И РЕДАКТИРОВАНИЕ ГЕОМЕТРИИ ТЕЛ И МОДЕЛИ

8.1 Cистемы координат

8.1.1 Глобальная система координат

Геометрия модели задается в правой декартовой системе координат.

Для транспортных средств рекомендуется использовать ориентацию машины относительно системы координат, принятую в конструкторской документации. Для автомобилей:

1.
ось ’X’ - горизонтальная продольная, параллельная плоскости симметрии машины, положительное направление - против хода машины,
2.
ось ’Y’ - горизонтальная поперечная, перпендикулярная оси ’X’, положительное направление справа налево, если смотреть спереди ,
3.
ось ’Z’ - вертикальная, положительное направление вверх, начало координат находится в произвольной точке машины.

Для железнодорожного транспорта:

1.
ось ’X’ - горизонтальная продольная, параллельная плоскости симметрии машины, положительное направление - по ходу движения,
2.
ось ’Y’ - горизонтальная поперечная, перпендикулярная оси ’X’, положительное направление слева напрво, если смотреть спереди ,
3.
ось ’Z’ - вертикальная, положительное направление вверх, начало координат находится в произвольной точке машины.

Примечание. Данные условия необязательны, но целесообразны при работе с моделями транспортных средств.

8.1.2 Локальная система координат

Локальная система совпадает с глобальной с точностью до параллельного переноса, каждая подсистема (тело) может описываться в своей локальной системе координат. Начало системы координат выбирается произвольно (учитывается удобство задания координат узлов подсистем).

8.2 Занесение тел в библиотеку

Программы ФРУНДа могут получать доступ к геометрии тел только через библиотеку тел. Поэтому для использования того или иного тела в модели нужно предварительно описать его геометрию и занести ее в библиотеку.

Описание геометрии можно производить двумя способами.

Можно задавать геометрию непосредственно в формате входного файла библиотеки и затем вносить ее в библиотеку с помощью программы avmodel, при работе из оболочки пункт меню ‘Ввод геометрии тел в библиотеку’. При инициализации этого пункта меню запрашивается имя фала с геометрией тела в формате библиотеки - расширение lba.

Более удобно описывать геометрию с помощью программы редактирования геометрии и модели - vizual, пункт меню в оболочке - ‘Редактирование геометрии и структуры’, из которой, затем ,можно ввести тело в библиотеку.

8.2.1 Формат данных для занесения тела в библиотеку

Геометрия каждого тела записывается в отдельный файл в следующем формате:

1.
1-я строка - код вида работ целое число (ц. ч.), может принимать следующие значения
(a)
1 - инициализация библиотеки,
(b)
2 - запись тела в библиотеку,
(c)
3 - стирание тела из библиотеки,
(d)
4 - контроль параметров тела из библиотеки (выдает информацию о теле в файл ’avmodel.ls’ в форме файла исходных данных ),
(e)
6 - контроль состояния всей библиотеки,
(f)
7 - результат контроля смотри в ’avmodel.ls’,
(g)
8 - очистка библиотеки ( после выполнения этого вида работ библиотека будет в состоянии инициализации ).
2.
2-я строка - тип тела ц. ч. (ITIP), может принимать следующие значения
(a)
ITIP< -1000 - упругое плоское тело,
(b)
-1000 < = ITIP < -100 - твердое плоское тело,
(c)
-100 <= ITIP <= -1 - твердые плоские мосты,
(d)
1 <= ITIP <= 100 - твердые объемные мосты,
(e)
100 < ITIP <= 1000 - твердое объемное тело,
(f)
1000 < ITIP упругое объемное тело. Примечание. Рекомендуется при решении плоских задач пользоваться объемными телами, для большей универсальности визуализации геометрии.
3.
3-я строка - масса тела (в. ч.),
4.
4-я строка - имя тела (алфавитно-цифровая информация, не боле 40 символов),
5.
5-я строка - 6 чисел - моменты инерции тела относительно осей X,Y,Z (считается, что это главные центральные оси), остальные - три вещественных нуля.
6.
6-я строка - 6 чисел
(a)
общее число узлов в описании тела (ц. ч.),
(b)
число узлов крепления в теле без учета узла центр масс (ц. ч.),
(c)
число узлов контакта с дорогой, (ц. ч.)
(d)
число стержней в описании геометрии тела (ц. ч.),
(e)
число узлов в конечно-элементной модели для упругого тела (для твердого тела 0),
(f)
число учитываемых форм колебаний упругого тела (для твердого тела 0),
7.
7-я строка - 4 числа,
(a)
номер узла ц. м. подсистемы,
(b)
координаты узла ц. м. в локальной системе, координат - 3 числа.
8.
C 8 строки вводятся координаты узлов подсистемы, причем первыми вводятся узлы крепления, затем узлы контакта с дорогой (считается , что узлы контакта входят в число узлов крепления),узел привязки системы координат, потом остальные узлы (координаты вводятся по осям X Y Z, в мм).

После описания геометрии узлов вводятся стержни подсистемы - парами узлов. Число стержней должно быть равно их числу, заданному в 6 строке. В конце файла можно вводить любое количество строк комментариев. Информация об инерционных параметрах и геометрии тела автоматически используется при расчете, в случае отсутствия ее в явном виде в файле описания параметров модели.

8.3 Редактирование геометрии

Редактирование геометрии тел производится с помощью программы vizual, пункт меню оболочки ‘Редактирование геометрии и структуры’. В том случае файлы геометрии тел имеют расширение geo. Программа vizual предлагает выбрать входной файл из списка или образовать новый.

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

клавиша

операция

f1

подсказка к программе vizual,

x/X

поворот картинки вокруг оси x вперед/назад,

y/Y

поворот картинки вокруг оси y вперед/назад,

z/Z

поворот картинки вокруг оси z вперед/назад,

?

сместить картинку вверх,

?

сместить картинку влево,

?

сместить картинку вправо,

?

сместить картинку вниз,

n/N

показать/убрать номера узлов,

l/L

тонкие/толстые линии,

w/W

уменьшить/увеличить изображение,

c

сменить цвет линий,

C

сменить фон экрана,

a/A

добавить новое тело, программа позволяет одновременно отображать несколько тел, последнее введенное тело считается текущим,

o/O

выкл/вкл раскраску тел разными цветами,

t/T

поместить изображение в центр,

p

сбросить картинку в файл для печати,

^Z

завершить работу,

Esc

выход из подсказки,

f/F

создать файл тела для библиотеки и ввести тело в библиотеку,

e/E

редактировать файл с геометрией тела, по этой команде текущий файл с геометрией помещается в редактор текстов, в котором можно его редактировать,

D

удалить последнее тело,

S

выбрать и увеличить фрагмент изображения,

P

приложить к телу весовую нагрузку (изменение структуры модели),

u/U

выкл/вкл показ центров масс тел,

M

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

m

отсоединить файл структуры расчетной схемы,

k/K

добавить в модель соединительный элемент,

v/V

выкл/вкл показ соединительных элементов,

R

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

r

отсоединить файл параметров модели,

b/B

добавить тело в файл структуры модели,

Alt d

удалить тело из файла структуры,

^d

вернуть тело в файл структуры,

^k

сформировать контрольную точку для вывода кинематических параметров движения тела,

^f

сформировать контрольную точку для вывода параметров сил и деформаций соединительных элементов,

Alt s

сместить геометрию тела параллельно,

Alt a

добавить новый узел в геометрию,

Alt m

модифицировать координаты узла в геометрии,

Alt l

добавить новый стержень в геометрию,

Alt u

удалить существующий стержень,

g/G

добавить макро в модель,

Q/q

вкл/выкл показ макро в модели,

Alt r

добавить окружность.

Программа vizual позволяет одновременно визуализировать несколько тел, входящих в модель и редактировать структуру и параметры модели для изображаемых тел.

Глава 9
БОЛЬШИЕ ДВИЖЕНИЯ СИСТЕМ ТЕЛ, ДИНАМИКА МЕХАНИЗМОВ

9.1 Общие положения

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

Как везде, во ФРУНДе используется правая система координат. Система углов Эйлера - последовательный поворот вокруг осей X,Y,Z. Движение упругого тела описывается в нормальных координатах с учетом форм движения как твердого тела. Для вычисления матрицы поворота задается взаимное расположение трех любых узлов тела, не лежащих на одной прямой. Для моделирования динамики больших движений тел используются дополнительные типы элементов расчетных схем - кинематические пары и пружины в шарнирах - упруго - демпфирующие элементы, направление действия сил в которых совпадает с линией, соединяющей узлы двух тел.

9.2 Задание кинематических пар

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

Тип кинематической пары задается в механической характеристике. Код типа функции упругого слагаемого составляющей механической характеристики в этом случае должен быть равен 15 - кинематическая связь. Число параметров может принимать значения 1,2,3 или 4, смысл которых расписывается ниже.

В случае задания в механической характеристике по некоторому направлению такого упругого слагаемого демпфирующее слагаемое игнорируется.Принципы моделирования различных кинематических пар приведены в главе 7

Положение кинематической пары в модели задается путем описания соответствующего соединительного элемента со ссылкой на необходимую механическую характеристику. Для каждого соединительного элемента можно задавать свою вспомогательную систему координат - п.п. 9.3.19.3.1,т.е. каждая кинематическая пара одного типа в модели может быть ориентирована произвольно. Считается, что такая вспомогательная система координат связана с первым телом, указываемом в соединительном элементе. Число параметров в типе упругого слагаемого 15 имеет следующий смысл:

1.
координаты точек первого и второго тела , входящих в кинематическую пару, фиксированы, уравнения связей строятся в проекции на систему координат, связанную с первым телом, совпадающую с подвижными осями координат первого тела или повернутую относительно ее на постоянные.
2.
координаты точки первого тела фиксированы, а координаты точки второго тела вычисляются из условия ее совпадения с точкой первого тела в неподвижной системе координат (данный вид связи не применяется),
3.
Уравнения связей строятся в проекции на неподвижную инерциальную систему координат,
4.
координаты точек первого и второго тела , входящих в кинематическую пару, вычисляются из условия совпадения с некоторой точкой на третьем теле, относительно которого эти два тела могут перемещаться. Такой точкой может быть, например, геометрический центр зубчатого зацепления, для моделирования которого предназначен данный способ построения уравнений связи.

9.3 Дополнительная информация для расчета больших движений

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

Вспомогательная информация для шага МОДЕЛЬ задается в файле с именем modhelp, для шага РАСЧЕТ - rashelp (в оболочке файлы с расширениями mhl и rhl). В указанных файлах должна содержаться информация об элементах расчетной схемы, необходимых для расчета больших движений. Если часть тел модели совершает малые движения, то соответствующие элементы не указываются в дополнительной информации . Если указанные вспомогательные файлы отсутствуют, то движения тел модели считаются малыми.

9.3.1 Дополнительная информация для шага МОДЕЛЬ

Описание тел, совершающих большие движения

Тип строки - 1. Далее идут по порядку 8 целых чисел:

1.
номер подсистемы;
2.
код типа движения тела, вводятся следующие коды:
(a)
1 - малые движения с единичной матрицей, поворота,
(b)
3 - большие движения с точной матрицей поворота,
3.
код типа вращательной степени свободы вокруг оси X, вводятся следующие коды
(a)
0 - степень свободы подавлена,
(b)
1 - малые повороты, с линеаризованной матрицей поворотов,
(c)
3 - большие повороты, с точной матрицей поворотов,.
4.
код типа поворота вокруг оси Y,
5.
код типа поворота вокруг оси Z,
6.
номер узла 1 (для упругого тела),
7.
номер узла 2 (для упругого тела),
8.
номер узла 3 (для упругого тела).

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

Строки этого типа автоматически формируются в файле default.mhl при задании в блоке 14 признаков больших поворотов и номеров узлов в позициях признаков поступательных степеней свободы.

Описание соединительных элементов

Тип строки - 2. Далее идут по порядку 5 чисел:

1.
номер соединительного элемента,
2.
тип соединительного элемента, вводятся следующие типы:
(a)
0 - обычный для ФРУНДа, или кинематическая пара,
(b)
1 - "пружина в шарнирах считается, что механическая характеристика такого элемента имеет только одну составляющую по оси X,
3.
признак ориентации относительно оси X (0 или 1),
4.
признак ориентации относительно оси Y (0 или 1),
5.
признак ориентации относительно оси Z (0 или 1).

Признаки ориентации относительно осей могут задаваться только для кинематической пары. Признак ориентации означает, что система координат кинематической пары повернута относительно подвижной системы координат тела.

Строки с такой информацией для пружины в шарнирах автоматически заносятся в файл default.mhl при использовании макро tspring.mac.

Описание точек контакта и ориентации осей колес транспортных средств

Тип строки - 3. Далее идут по порядку 4 числа:

1.
номер тела «колесо»,
2.
номер узла на колесе, касающегося дороги, в плоскости вращения,
3.
номер узла на оси колеса в плоскости вращения,
4.
номер узла на оси вращения колеса в направлении к центру машины.

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

Строки с такой информацией для точек контакта автоматически заносятся в файл default.mhl при использовании макро kkoleso.mac.

9.3.2 Дополнительная информация для шага РАСЧЕТ

Задание координат вспомогательных узлов упругих тел

Тип строки -1. Далее идут по порядку 7 чисел:

1.
номер подсистемы,
2.
координата X второго узла относительно первого - п.п. 9.3.1,
3.
координата Y второго узла относительно первого,
4.
координата Z второго узла относительно первого,
5.
координата X третьего узла относительно первого,
6.
координата Y третьего узла относительно первого,
7.
координата Z третьего узла относительно первого.

Данная информация является не обязательной (берется автоматически из геометрии).

Задание начального положения тел.

Тип строки -2. Далее идут по порядку 7 чисел:

1.
номер подсистемы,
2.
X координата ц. м. тела в неподвижной (абсолютной) системе координат,
3.
Y координата ц. м. тела в неподвижной системе координат,
4.
Z координата ц. м. тела в неподвижной системе координат,
5.
начальный угол поворота относительно оси X,
6.
начальный угол поворота относительно оси Y,
7.
начальный угол поворота относительно оси Z.

Предназначено для приближенного задания первоначального положения тел при моделировании механизмов (больших движений).

Задание специальных характеристик соединительных элементов

Тип строки -3. Далее идут по порядку 7 чисел:

1.
номер соединительного элемента,
2.
угол поворота относительно оси X - всегда 0,
3.
угол поворота относительно оси Y - всегда 0,
4.
угол поворота относительно оси Z - всегда 0,
5.
свободная длина соединительного элемента,
6.
смещение по оси Y - всегда 0,
7.
смещение по оси Z - всегда 0.

Данная информация может задаваться для кинематической пары или с. э. типа "пружина в шарнирах".

Для кинематической пары имеют значение углы поворота относительно осей. Такие углы определяют ориентацию местной системы координат кинематической пары и задаются в случае, если были заданы признаки поворота в п. п. 9.3.1

Для пружины в шарнирах имеет значение свободная длина соединительного элемента - которая равна длине пружины в несжатом состоянии. Если для пружины такой параметр не будет задан, то ее свободная длина считается равной нулю.

Строки с такой информацией для пружин в шарнирах автоматически заносятся в файл default.rhl при использовании макро lspring.mac.

9.4 Задание степеней свободы

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

9.5 Задание начальных условий

При моделировании с учетом больших движений является существенным взаимное начальное положение тел. В простейших случаях положение центра масс тела в неподвижной системе координат и начальные углы поворота согласно можно задать с помощью дополнительной информации - п.п. 9.3.2.Наиболее универсальный способ определения начальных условий - применение процедуры сборки см. главу 3.Для упругого тела начальные условия находятся так, как описано в п. п. 3.6. При расчете больших движений, программа расчета может сама корректировать начальные условия в соответствии с заданной в управляющей информации скоростью движения. Такая корректировка заключается в добавке к переменным скоростейтел по направлению X скорости движения со знаком минус, и производится по запросу программы. Необходимость в этом возникает, например, при моделировании движения автомобиля с различными скоростями. Автоматическое добавление к начальным условиям скорости, позволяет пользоваться одними начальными условиями, определенными в покое, для разных скоростей движения.

9.6 Задание нагрузок

В случае больших движений тел предполагается, что силы от нагрузок - специальные силы в терминологии ФРУНДа, могут быть также приложены в любой точке тела. При этом считается, что силы действуют всегда в направлении неподвижной системы координат, а моменты - в системе координат, связанной с телом. Такое задание внешних сил удобно для задания весовых и инерционных нагрузок.

9.7 Податливые связи

При моделировании динамики больших движений тел могут быть использованы, кроме кинематических пар, податливые связи двух типов. К таким связям относятся "пружина в шарнирах"и обычные для ФРУНДа многомерные податливые связи. Кроме того есть способ задания упругости в самой кинематической паре по выбранным направлениям.

9.7.1 Пружина в шарнирах

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

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

9.7.2 Многомерный податливый элемент

При больших движениях тел связи между телами можно задавать обычным образом - соединительным элементом, передающим силы и моменты по направлениям подвижной системы координат одного из тел.В таком случае принимается, что взаимные смещения тел друг относительно друга малы, и соединительный элемент считается "вмороженным"в первое тело, указанное в описании соединительного элемента. Если первое тело совершает малые движения или неподвижно, то программа считает первым