Разработка программного комплекса для расчета с элементами оптимизации конических передач

1.АНАЛИЗ ПРЕДМЕТНОЙ ОБЛАСТИ

1.1 ОПИСАНИЯ ПРЕДМЕТНОЇ ОБЛАСТІ

1.2 АЛГОРИТМ РЕШЕНИЯ ЗАДАЧ НАОПТИМИЗАЦИЮ КОНИЧЕСКИХ ПЕРЕДАЧ

2. ОБОБЩЕННАЯ ПОСТАНОВКА ЗАДАЧИ

2.1 ОПТИМИЗАЦИЯ ПАРАМЕТРОВ

3. РАЗРАБОТКА ПРОГРАММНОГО ПРОДУКТА

3.1 ОСНОВАНИЯ К РАЗРАБОТКЕ ПП ПО ОПТИМИЗАЦИИ КОНИЧЕСКОЙ ПЕРЕДАЧИ

3.2 ТЗ ДЛЯ РЕШЕНИЯ ЗАДАЧИ ОПТИМИЗАЦИИ КОНИЧЕСКОЙ ПЕРЕДАЧИ

3.2.1 ВСТУПЛЕНИЕ

3.2.2 ОСНОВАНИЕ ДЛЯ РАЗРАБОТКИ ПП ПО ОПТИМИЗАЦИИ КОНИЧЕСКОЙ ПЕРЕДАЧИ

3.2.3 НАЗНАЧЕНИЕ РАЗРАБОТКИ ПП ПО ОПТИМИЗАЦИИ КОНИЧЕСКОЙ ПЕРЕДАЧИ

3.2.4 ТРЕБОВАНИЯ К ПП ПО ОПТИМИЗАЦИИ КОНИЧЕСКОЙ ПЕРЕДАЧИ

3.2.5 ТРЕБОВАНИЯ К ФУНКЦИОНАЛЬНЫМ ХАРАКТЕРИСТИКАМ

3.2.6 ТЕХНИКО-ЭКОНОМИЧЕСКИЕ ПОКАЗАТЕЛИ ПП ПО ОПТИМИЗАЦИИ КОНИЧЕСКОЙ ПЕРЕДАЧИ

3.2.7 ПОРЯДОК КОНТРОЛЯ И ПРИНЯТИЯ ПП ПО ОПТИМИЗАЦИИ КОНИЧЕСКОЙ ПЕРЕДАЧИ

3.2.8 ВЫВОДЫ О ВЫБОРЕ МОДЕЛИ СОЗДАНИЯ ПРОГРАММНОГО ПРОДУКТА

3.3 ВЫБОР ЯЗЫКА ПРОГРАММИРОВАНИЯ

3.4 ДИАГРАММА ПРЕЦЕДЕНТОВ

3.5 МОДУЛИ И ИХВЗАИМОДЕЙСТВИЕ МЕЖДУ СОБОЙ

3.6 РУКОВОДСТВО ПОЛЬЗОВАТЕЛЯ

3.6.1 УСЛОВИЯ ВЫПОЛНЕНИЯ ПРОГРАММЫ

3.6.2 ВЫПОЛНЕНИЕ ПРОГРАММЫ

3.7 СПРАВОЧНАЯ СИСТЕМА

4. ПРИМЕР РАСЧЕТА КОНИЧЕСКОЙ ПЕРЕДАЧИ ВЫВОД ИСПОЛЬЗОВАННАЯ ЛИТЕРАТУРА

ВВЕДЕНИЕ

В наше время одним из главных путей повышения эффективности проектно-конструкторских работ есть автоматизация проектирования на основе применения ЭВМ. Одно из важнейших направлений автоматизации проектирования конструкций и технологий — оптимальное проектирование. Для решения задания оптимизации математических моделей объектов или систем предложено достаточно много разных алгоритмов и их модификаций.

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

Основная цель курсового проекта по дисциплине «Математические методы исследования операций» — получение навыков разработки интеллектуального программного продукта для решения задания оптимизации (поддержки решения) в заданной предметной области

1. Анализ предметной области

1.1 Описание предметной области

Конические передачи

Конические зубчатые колеса применяют в передачах, оси валов которых пересекаются под некоторым углом. Обычно (рис 1).

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

Преимущества:

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

· Возможность передачи движения между звеньями с переменным межосевым углом при широком диапазоне его изменения;

· Расширение компоновочных возможностей при разработке сложных зубчатых и комбинированных механизмов.

Недостатки:

· Более сложная технология изготовления и сборки зубчатых колес;

· Большие осевые и изгибные нагрузки на валы, особенно в связи с консольным расположением зубчатых колес.

1.2 Алгоритм решения задачи расчета плоскоконической передачи

1) Выбираем материал и рассчитываем допускаемое напряжение. Допускаемое контактное напряжение:

Для шестерни:

=2ННВ+70;

SH=1,1;

принимается КHL=1;

;

Для колеса:

=2ННВ+70;

SH=1,1;

;

Для прямозубых передач в качестве расчетного контактного напряжения принимается меньше ;

2) Определяем коэффициенты

KH=KHбKHвKHV,

де KHб=1 (прямые зубья), KHв=1,1 (при b2/d1=0,3;колеса прикрепляются), KHV=1,2 (принято приблизительно).

3) Определяем диаметр колеса

4) Определяем модуль

5) Определяем количество зубьев

6) Определяем геометрические параметры зубчатых колес

de1=mZ1; de2=mZ2;

д2=arctgU1−2ф; д1=90є-д2;

Внешние диаметры отступов и впадин зубьев:

da1=de1+2m· cosд1; da2=de2+2m· cosд2;

df1=de1−2,4m· cosд1; df2=de2−2,4m· cosд2;

Внешняя коническая длина

Средняя коническая длина

Rm=Re-0,5b.

Средний модуль зубьев

Средние делительные диаметры шестерни и колеса:

dm1=mZ1; dm2=mZ2;

7) Определяем степень точности зубчатых колес

В зависимости от угловой скорости

2. ОБОБЩЕННАЯ ПОСТАНОВКА ЗАДАЧИ

2.1 Оптимизация параметров

Оптимизация параметров — определение таких номинальных значений внутренних параметров X (…) объекта проектирования, при которых функция f (x) — целевая функция, функция качества, — принимает экстремальное (min и max) значение. В процессе решения задачи оптимизации обычно необходимо найти оптимальные значения некоторых параметров, которые определяют это задание. При решении инженерных задач их принято называть проектными параметрами, а в экономических задачах их обычно называют параметрами плана. В качестве проектных параметров могут быть, в частности, значение линейных размеров объекта, массы, температуры, и тому подобное Число n проектных параметров x1, x2 …, Xn характеризует размерность (и степень сложности) задачи оптимизации

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

Целевую функцию можно записать в виде: U = f (x1, x2 …, xn). Примерами целевой функции, которые встречаются в инженерных и экономических расчетах, является прочность или масса конструкции, мощность установки, объем выпуска продукции, стоимость перевозок грузов, прибыль и тому подобное В случае одного проектного параметра (n = 1) целевая функция является функцией одной переменной, и ее график — некоторая кривая на плоскости. При n = 2 целевая функция является функцией двух переменных, и ее график — поверхность в трехмерном пространстве. Целевая функция не всегда может быть представлена в виде формулы. Иногда она может принимать только некоторые дискретные значения, задаваться в виде таблицы и тому подобное. Во всех случаях она должна быть однозначной функцией проектных параметров. Целевых функций может быть несколько. Например, при проектировании изделий машиностроения одновременно нужно обеспечить максимальную надежность, минимальную материалоемкость, максимальный полезный объем (или грузоподъемность). Некоторые целевые функции могут оказаться несовместимыми. В таких случаях необходимо вводить приоритет той или другой целевой функции.

В этой работе в качестве входных параметров были избраны: P1 — мощность, n — частота вращения, U — передаточное число, P2- мощность. В качестве проектных параметров были избраны: Qсила, которая действует на валы конической передачи, lp — диаметр колеса, Pмодуль зубьев. В качестве целевой функции была избранная функция, которая учитывает эти три критерия (рисунок 2). При этом каждый из критериев имеет свой коэффициент значимости, которые пользователь может задавать (чем более важный параметр, тем больше коэффициент значимости). Целевая функция вычисляется по такой формуле 1.

f (P, Q, lp) = k1 * P + k2 * Q + k3 * lp, (1)

де k1, k2, k3 — коэффициенты приводят исходные параметры к безразмерным, чем более важнее параметр, тем выше коэффициент.

Рисунок 2 — Входные и выходные параметры

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

3. РАЗРАБОТКА ПРОГРАММНОГО ПРОДУКТА

3.1 Основания к разработке ПП по оптимизации конической передачи

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

3.2 ТЗ для решения задачи оптимизации конической передачи

проектирование программа алгоритм передача

3.2.1 Вступление

Программный продукт (ПП) представляет собой модель автоматизированной информационной системы для оптимизации расчетов конической передачи. ПП может использоваться как часть разработки при создании АИС для современных проектов и разработок.

3.2.2 Основание для разработки ПП по оптимизации конической передачи

Разработка ведется на основе индивидуального задания к курсовой работе по дисциплине «Математические методы исследования операций» Организация, которая утвердила настоящий документ: ДГМА

Дата утверждения: 1 апреля 2103 г.

Тема разработки: «ПП по оптимизации конической передачи».

3.2.3 Назначение разработки ПП по оптимизации конической передачи

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

3.2.4 Требования к ПП по оптимизации конической передачи

3.2.4.1 Требования к функциональным характеристикам

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

а) обеспечить наглядное представление входной и расчетной информации;

б) выводить значение оптимизации целевой функции

в) иметь интуитивно понятный интерфейс.

3.2.4.2 Требования к надежности ПП по оптимизации конической передачи

Среди требований к надежности необходимо выделить следующие: а) ПК не должен приводить к сбоям в работе ОС;

б) ПК должен обеспечивать обработку информации, введенной пользователем.

3.2.4.3 Условия эксплуатации ПП по оптимизации конической передачи

Условия эксплуатации ПМК определяются СанПиН 2.2.2 545−96 «Гигиенические требования к видеодисплейным терминалам, персональным вычислительным машинам и организации работы».

3.2.4.4 Требования к составу и параметрам технических средств

CPU: 1GHz или выше ОЗП: 512 MB или выше Видеоадаптер та монитор: VGA (640×480) Super VGA (800×600)

Свободное место на HDD: 50Мб или больше.

Устройства взаимодействия клавиатура и мышь

3.2.4.5 Требования к информационной и программной совместимости Главным требованием для данного продукта является наличие языка программирования Borland Delphi 7.0 и независимость от аппаратной платформы.

3.2.5 Требования к программной документации

Программная документация должна содержать:

а) рабочий проект ПМК;

б) исходных коды ПМК с комментариями;

в) руководство пользователя;

г) контекстно-зависимую помощь;

д) руководство по установлению ПМК.

3.2.6 Технико-экономические показатели ПП по оптимизации конической передачи

Экономическая эффективность от внедрения ПП обеспечивается за счет сокращения расходов времени на ручные расчеты конических передач.

Стоимость разработки ПП становить 300 грн.

Стадии и этапы разработки приведены в таблице 1.

Таблица 1 — Стадии и этапы разработки

3.2.7 Порядок контроля и принятия ПП по оптимизации конической передачи

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

3.2.8 Выводы о выборе модели создания программного продукта

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

3.3 Выбор языка программирования

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

3.4 Диаграмма прецедентов

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

Рисунок 3 — UML-диаграмма «Ввод параметров для расчета и оптимизации конической передачи»

Основные функции управления ПП из оптимизации конической передачи

· задать начальные параметры передачи;

· задать коэффициенты важности составляющих целевой функции;

· рассчитать коническую передачу;

· построить график зависимости целевой функции от ее параметров;

· оптимизировать расчеты конической передачи;

· сохранить полученные данные в текстовом файле.

3.5 Модули и их взаимодействие между собой

Диаграмма переходов наведена на рисунку 6.

Рисунок 4 — Диаграмма переходов состояний (STD)

Рисунок 5 — Контекстная диаграмма программного продукта

3.6 Руководство пользователя

3.6.1 Условия выполнения программы

Для запуска программы убедитесь в соответствие вашего компьютера системным требованиям, изложенным в пункте 3.1. Если Ваш компьютер имеет такие ресурсы, то из диска, прилагаемому к курсовой работе, скопируйте файл Project1. exe ни запустите его. После установления программы на Вашем компьютере приступайте к работе.

3.6.2 Выполнение программы

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

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

В результате проектирования был разработан программный продукт, название которого Project1.exe. Дополнение не требует для своей работы никакого дополнительного ПО.

Для выполнения задания к курсовой работе был выбранный язык Object Pascal (среда разработки Delphi7).

Для поддержки работы данного разработанного продукта нужное наличие следующего ПО: Windows 2000/XP/2003, Vista, Windows 7, Windows 8.

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

Кнопка «Оптимизировать» проводит оптимизацию целевой функции

F (Q, lp, P) = K1*Q + K2*lp + K3*P,

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

Кнопка «Сохранить в файл» позволяет сохранить полученные данные в текстовом файле в каталоге программы (рис. 6).

Кнопки «Очистить» позволяют очистить поле для выведения информации и поле для построения графиков (рис. 6−7).

Рисунок 6 — Главная форма программы

Рисунок 7 — Главная форма программы

3.7 Справочная система

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

Рисунок 8 — Вікно довідки Справочная система состоит из двух разделов: «Общие сведения» и «Информация, по ПП».

4. ПРИМЕР РАСЧЕТА КОНИЧЕСКОЙ ПЕРЕДАЧИ Входные данные:

Р1 = 4,87 кВт Р2 = 4,63 кВт

n1 =1445 мин-1

n1 =802,8 мин-1

Кпер =1,8

U1,2 = 1,8

1) Выбираем материал и рассчитываем допускаемую нагрузку.

Принимается для изготовления шестерни и колеса Сталь 45 с термообработкой — улучшение.

Для шестерни — твердость поверхности зубьев Н1=269.302 НВ (наиболее достоверная твердость Н1=285 НВ), уВ1=890 МПа, уТ1=650 МПа при диаметре заготовки до 80 мм.

Для колеса — твердость поверхности зубьев Н2=235.262 НВ (наиболее достоверная твердость Н2=250 НВ), уВ1=780 МПа, уТ1=540 МПа при диаметре заготовки до 125 мм.

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

Для шестеренки:

=2ННВ+70=2*285+70=640 МПа;

SH=1,1;

принимается КHL=1;

;

ZR=1 (принято Ra=1б25ююю0б63);

ZV=1 (ждет V<5 м/с);

KL=1 (закрытая обильно смазочная передача);

KXH1 (ожидается диаметр колеса меньше 700 мм);

Для колеса:

=2ННВ+70=2*250+70=570 МПа;

SH=1,1;

принимается КHL=1;

;

Для прямозубых передач в качестве расчетного контактного напряжения принимается меньше ;

2) Определяем коэффициенты

KH=KHбKHвKHV,

де KHб=1 (прямые зубья), KHв=1,1 (при b2/d1=0,3;колеса прикрепляются), KHV=1,2 (принято приблизительно).

3) Определяем диаметр колеса Принимаем ближайшее стандартное de2=125мм.

4) Определяем модуль Учитывая, что для силовых конических передач рекомендовано столько зубьев шестерни Z1=17…20 при U=1,0…1,5 предварительно назначаем Z1=20.

Принимаем стандартно me=3,5 мм.

5) Определяем количество зубьев

(подбором стандартных значений модуля не удалось вычислить число зубьев колеса Z2);

6) Определяем геометрические параметры зубчатых колес

de1=mZ1=3,5· 20=70мм;

de2=mZ2=3,5· 36=126мм;

д2=arctgU1−2ф=arctg1,8=60,9454є;

д1=90є-д2=90є-60,9454є=29,0546є;

b=20мм;

Внешние диаметры выступлений и впадин зубьев:

da1=de1+2m· cosд1=70+2·3,5·cos20,0546є=76,129 мм;

da2=de2+2m· cosд2=126+2·3,5·cos60,9454є=129,40 мм;

df1=de1−2,4m· cosд1=70−2,4·3,5·cos20,0546є=62,66 мм;

df2=de2−2,4m· cosд2=126−2,4·3,5·cos60,9454є=121,92 мм;

Внешняя коническая длина Средняя коническая длина

Rm=Re-0,5b=72,0694−0,5· 20=62,0694 мм.

Средний модуль зубьев Средние делительные диаметры шестерни и колеса:

dm1=mZ1=3,014· 20=60,280 мм;

dm2=mZ2=3,014· 36=108,504 мм;

7) Определяем степень точности зубчатых колес В зависимости от угловой скорости

назначаем степень точности 8-В ГОСТ 1758–81.

ВЫВОД

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

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

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

Приложение А

unit Unit1;

interface

uses

Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,

Dialogs, ExtCtrls, StdCtrls, Buttons, TeEngine, Series, TeeProcs, Chart,

Menus, OleCtnrs, Math;

type

TForm1 = class (TForm)

Panel1: TPanel;

Label11: TLabel;

GroupBox1: TGroupBox;

Label1: TLabel;

Label5: TLabel;

Label6: TLabel;

Label7: TLabel;

Edit1: TEdit;

Edit2: TEdit;

Edit6: TEdit;

Edit7: TEdit;

GroupBox2: TGroupBox;

Label8: TLabel;

Label9: TLabel;

Label10: TLabel;

Edit8: TEdit;

Edit9: TEdit;

Edit10: TEdit;

RadioGroup1: TRadioGroup;

RadioGroup2: TRadioGroup;

BitBtn1: TBitBtn;

BitBtn2: TBitBtn;

Label12: TLabel;

Chart1: TChart;

MainMenu1: TMainMenu;

N1: TMenuItem;

N2: TMenuItem;

N3: TMenuItem;

BitBtn3: TBitBtn;

BitBtn4: TBitBtn;

N4: TMenuItem;

OleContainer1: TOleContainer;

Edit4: TEdit;

Label3: TLabel;

RadioGroup4: TRadioGroup;

Memo1: TMemo;

Series1: TLineSeries;

Series2: TLineSeries;

Series3: TLineSeries;

N5: TMenuItem;

procedure BitBtn1Click (Sender: TObject);

procedure FormCreate (Sender: TObject);

procedure BitBtn3Click (Sender: TObject);

procedure BitBtn4Click (Sender: TObject);

procedure BitBtn2Click (Sender: TObject);

procedure N2Click (Sender: TObject);

procedure N4Click (Sender: TObject);

procedure N5Click (Sender: TObject);

private

{ Private declarations }

public

{ Public declarations }

end;

var

Form1: TForm1;

P1, U, n, Q, k1, k2, k3, max_func, min, d1op, v, d2op, lp, P, P0, C0,Ca, Cv, Cp, a1, Ft, b, d1,d2,e, aop, F0, max_Q, max_p, max_lp, func, min2, x, x2: real;

j, i_min, ap: integer;

d_tabl, d_tabl2,mas_raz1, mas_raz2: array [0.11] of real;

implementation

{$R *.dfm}

procedure TForm1. BitBtn1Click (Sender: TObject);

var i: integer;

begin

aop:=0;

min:=10 000;

max_func:=-10 000 000;

d_tabl[0]: =71.0;

d_tabl[1]:=80.0;

d_tabl[2]:=90.0;

d_tabl[3]:=100.0;

d_tabl[4]:=112.0;

d_tabl[5]:=125.0 ;

d_tabl[6]: =140.0 ;

d_tabl[7]: =160.0 ;

d_tabl[8]: =180.0 ;

d_tabl[9]: =200.0 ;

d_tabl[10]: =224.0 ;

d_tabl[11]: =250.0 ;

P1:=StrToFloat (Edit1.Text);

n:=StrToFloat (Edit2.Text);

j:=StrToInt (Edit7.Text);

U:=StrToFloat (Edit6.Text);

P0:=StrToFloat (Edit4.Text); // 72,73)

k1:=StrToFloat (Edit8.Text);

k2:=StrToFloat (Edit9.Text);

k3:=StrToFloat (Edit10.Text);

d1op:=1100*exp (1/3*ln (P1/n)); //

for i:=0 to 11 do //

begin

x:=d1op-d_tabl[i];

mas_raz1[i]:=abs (x);

end;

for i:=0 to 11 do

begin

if (mas_raz1[i]

begin

min:=mas_raz1[i];

i_min:=i;

end;

end;

d1:=d_tabl[i_min]; //

v:=3.14*d1*n/60 000;

case RadioGroup1. ItemIndex of

0: e:=0.02;

1: e:=0.015;

end;

d2op:=U*d1*(1-e); //

min:=1 000 000;

for i:=0 to 11 do //

begin

x2:=d2op-d_tabl[i];

mas_raz2[i]:=abs (x2);

end;

for i:=0 to 11 do

if (mas_raz2[i]

begin

min:=mas_raz2[i];

i_min:=i;

end;

d2:=d_tabl[i_min];

case RadioGroup1. ItemIndex of

0: aop:=2*(d1+d2);

1: aop:=d1+d2;

end;

while aop<=5500 do

begin

lp:=2*aop+0.5*3.14*(d1+d2)+sqr (d2-d2)/4*aop;

case RadioGroup2. ItemIndex of

0: C0:=1; //

1: C0:=0.9;

2: C0:=0.8;

end;

a1:=180-((d2-d1)/aop)*57;

Ca:=1−0.003*(180-a1);

Cv:=1; //

case RadioGroup4. ItemIndex of

0: Cp:=1.6; //

1: Cp:=1.7;

2: Cp:=2;

end;

P:=P0*C0*Ca*Cv/Cp; //

Ft:=1000*P1/v; //

b:=Ft/P;

b:=round (b);

F0:=1.7*b*j; //

Q:=abs (2*sin (a½)*F0);

func:=k1*Q + k2*lp+P*k3;

Chart1.Series[0]. AddXY (Q, func);

Chart1.Series[1]. AddXY (lp, func);

Chart1.Series[2]. AddXY (P, func);

if func > max_func then

begin

max_func:=func;

max_Q:=Q;

max_lp:=lp;

max_p:=P;

end;

Memo1.Lines.Add ('При aop='+FloatToStrF (aop, ffFixed, 5,2)+

' Q='+FloatToStrF (Q, ffFixed, 5,2)+

' lp='+ FloatToStrF (lp, ffFixed, 5,2)+

' P='+FloatToStrF (P, ffFixed, 5,2));

aop:=aop+100;

end;

Memo1.Lines.Add (' ');

Memo1.Lines.Add (' РЕЗУЛЬТАТЫ ОПТИМИЗАЦИИ:');

Memo1.Lines.Add (' ');

Memo1.Lines.Add ('Исходные данные: ');

Memo1.Lines.Add (' ');

Memo1.Lines.Add (' P1, кВт = '+FloatToStrF (strtofloat (edit1.Text), ffFixed, 5,2));

Memo1.Lines.Add (' n1, 1/мин = '+FloatToStrF (strtofloat (edit2.Text), ffFixed, 5,2));

Memo1.Lines.Add (' U = '+FloatToStrF (strtofloat (edit6.Text), ffFixed, 5,2));

Memo1.Lines.Add (' j = '+FloatToStrF (strtofloat (edit7.Text), ffFixed, 5,2));

Memo1.Lines.Add (' P0 = '+FloatToStrF (strtofloat (edit4.Text), ffFixed, 5,2));

Memo1.Lines.Add (' K1 = '+FloatToStrF (strtofloat (edit8.Text), ffFixed, 5,2));

Memo1.Lines.Add (' K2 = '+FloatToStrF (strtofloat (edit9.Text), ffFixed, 5,2));

Memo1.Lines.Add (' K3 = '+FloatToStrF (strtofloat (edit10.Text), ffFixed, 5,2));

Memo1.Lines.Add (' ');

Memo1.Lines.Add (' F (Q, lp, P)=' + FloatToStrF (max_func, ffFixed, 5,2)+ #13#10 +

end;