Вычислительной программы

Анализ показывает, что нули главных миноров матрицы Н (Я) строго разделяются, а упорядоченная совокупность главных миноров этой матрицы обладает свойством последовательности Штурма. Указанное служит основой эффективной вычислительной процедуры для локализации собственных значений Т7^1'-моделей [2]. Для многомерных моделей эта процедура по быстродействию и затратам оперативной памяти ЭВМ существенно превосходит наиболее прогрессивные современные вычислительные схемы, базирующиеся на методах К. Якоби, В. Гивенса, А. Хаус-холдера [3]. Помимо эффективного определения собственных зна-

В настоящее время наиболее эффективный в вычислительном отношении алгоритм решения проблемы собственных значений симметричных матриц произвольной структуры базируется на методе Хаусхолдера ортогонального подобного приведения анализируемой матрицы к трехдиагональному виду. Трехдиаго-нализация (п X п)-матрицы А осуществляется на основе неитерационной вычислительной процедуры, состоящей из п — 2 шагов последовательных преобразований подобия исходной матрицы А. На каждом шаге в качестве матриц преобразования используются ортогональные матрицы отражения Р следующего вида [95]:

упрощения вычислительной процедуры и снижения ее погрешностей начальный радиус полости выбирался несколько завышенным (Ю-5 м) по сравнению с рекомендацией /32/, где указано значение (0.5-3)-10'6 м. Используя указанные допущения, в /12/ найдены распределения а,(г)к

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

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

При расчете вязкости газовой смеси весьма эффективно (в отношении точности расчета и компактности вычислительной процедуры) использовать формулу [103]

А <г- ¦ — для значения локальной переменной, область видимости которой ограничена местом ее ввода и концом вычислительной процедуры, отмеченной вертикальной чертой, фиксирующей начало и конец программного блока;

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

А <— П — ввод значения локальной переменной, область видимости которой ограничена местом ввода переменной и концом вычислительной процедуры, отмеченной вертикальной чертой, фиксирующей начало и конец программного блока (см. выше программы формирования функций v w и hss);

Универсальность вычислительной процедуры обеспечивается тремя следующими моментами:

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

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

Кавана [105] предложил общий способ для описания свойств композиционных материалов на базе экспериментальных данных и вычислительной программы, основанной па методе конечных элементов.

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

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

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

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

Не менее удовлетворительная сходимость между расчетными данными и результатами эксперимента получена при холодных проливках смесительных головок с 72 двухструйными од-нокомпонентными и 8 двухструйными двухкомпонентными форсунками при установке пробоотборника на расстоянии 3,5 и 4 см от смесительной головки (рис. 84). Однако опыты с пробоотборником, установленным на расстоянии 7—8 см от смесительной головки, показывают, что столкновения капель от разных форсунок и жгутование могут приводить к неравномерности расходонапряженности. Влияние этих столкновений учитывается в системе уравнений сохранения количества движения для капель того же компонента и группы диаметров, но не учитывается для капель из разных групп. Таким образом, использование вычислительной программы LISP применимо лишь! для зоны, отстоящей от смесительной головки не более -на 5 см.

Рис. 90. Подпрограммы вычислительной программы TDK [134].

Не менее удовлетворительная сходимость между расчетными данными и результатами эксперимента получена при холодных проливках смесительных головок с 72 двухструйными од-нокомпонентными и 8 двухструйными двухкомпонентными форсунками при установке пробоотборника на расстоянии 3,5 и 4 см от смесительной головки (рис. 84). Однако опыты с пробоотборником, установленным на расстоянии 7—8 см от смесительной головки, показывают, что столкновения капель от разных форсунок и жгутование могут приводить к неравномерности расходонапряженности. Влияние этих столкновений учитывается в системе уравнений сохранения количества движения для капель того же компонента и группы диаметров, но не учитывается для капель из разных групп. Таким образом, использование вычислительной программы LISP применимо лишь! для зоны, отстоящей от смесительной головки не более -на 5 см.

Рис. 90. Подпрограммы вычислительной программы TDK [134].

Изложенный алгоритм явился основой для создания вычислительной программы на языке ФОРТРАН-IV, которая реализована на ЭВМ ЕС-1022.

Изложенный алгоритм явился основой для создания вычислительной программы на языке ФОРТРАН-IV.