Central Asian Research Journal For Interdisciplinary Studies (CARJIS)
ISSN (online): 2181-2454 Volume 2 | Issue 5 | May, 2022 | SJIF: 5,965 | UIF: 7,6 | ISRA: JIF 1.947 | Google Scholar |
www.carjis.org DOI: 10.24412/2181-2454-2022-5-273-280
ВЫЧИСЛИТЕЛЬНЫЙ ЭКСПЕРИМЕНТ ПО РАСЧЕТУ ОПТИМИЗАЦИИ
АРОК ПО ВЕСУ
Сабир Халмурадович Якубов
Каршинский государственный университет, д.т.н., профессор
Исмаил Уролмахаматович Хушвоков
Денауский институт предпринимательства и педагогики, преподаватель
АННОТАЦИЯ
В данной работе представлены результаты расчета по оптимизации бесшарнирных арок при различных краевых условиях и системах внешних сил. Задача нахождения оптимальной формы арки в общем случае представляет собою определение из условия минимума веса арки при ограничении прочности, устойчивости, а также деформациям.
Ключевые слова: алгоритмы, оптимизация, целевая функция, толщина арки, вес арки.
ВВЕДЕНИЕ
Задачи оптимизации инженерных тонкостенных конструкций, предполагают использование широкого класса методов математического программирования от симплекс-алгоритма до глобальных алгоритмов случайного поиска. Постановка задач оптимизации и обратных задач расчета конкретных конструкций позволяет унифицировать методы их решения на основе применения различных методов. Стоит отметить, что подобные задачи обладают рядом особенностей по сравнению с абстрактными задачами математического программирования, что позволяет разработать новые алгоритмы или модифицировать известные методы с ускоренной сходимостью.
МЕТОДОЛОГИЯ
Из этих особенностей можно выделить следующее. Во-первых, при весовой оптимизации конструкций минимум целевой функции всегда находится на одном или пересечении ограничений по прочности, жесткости, устойчивости рассматриваемых конструкций. Эта особенность позволяет производить параметрическую адаптацию алгоритмов поиска. Во-вторых,
Central Asian Research Journal For Interdisciplinary Studies (CARJIS)
ISSN (online): 2181-2454 Volume 2 | Issue 5 | May, 2022 | SJIF: 5,965 | UIF: 7,6 | ISRA: JIF 1.947 | Google Scholar |
www.carjis.org DOI: 10.24412/2181-2454-2022-5-273-280
задача прямого расчета конструкции, как правило требует на несколько порядков больше затрат машинного времени, чем вычисление целевой функции. Отсюда - возможность структурной адаптации алгоритмов с целью максимально уменьшить количество прямых расчетов конструкций. В -третьих, как прямые расчеты, так и обратные и оптимизационные для достаточно сложных конструкций производится при помощи численных методов. При этом очевидна целесообразность соотношения точности расчетной модели (которая может выражаться в количестве членов ряда координатных функций, узлов разностной сетки, конечных элементов) и положения поисковой системы в области поиска.
При оптимизации арки по весу целевая функция F(X) есть функция, вычисляющая вес арки:
ßo
F (X) = \yRh(ß)dß (1)
0
В качестве оптимизируемых параметров мы принимали параметры, определяющие толщину арки и угол раствора арки [1].
Толщина арки h(ß) может выражаться только гладкой функцией, имеющей ограниченные вторые производные.
Наиболее близко подойти к оптимальной толщине, дающей минимум веса арки при заданной нагрузке и краевых условиях, можно, задавая функцию, определяющую толщину в виде полинома n-ной степени:
n
h ß ) = X nßi
i=0
здесь hi ( i = n) -оптимизируемые параметры. Задавая толщину в виде (2), можно, при достаточно большом n, с любой точностью приблизиться к оптимальному закону h*(ß).
Интересен другой подход к нахождению оптимального закона h*, разработанный Половинкиным А.И. [6]. Введенное им оптимизируемые параметры способно с любой точностью (при введении достаточного количества параметров) подойти к оптимальной толщине h*. При этом у него предусмотрено наложение определенных ограничений на свободу изменения оптимизируемых параметров с тем, чтобы оптимизируемая форма не была
Central Asian Research Journal For Interdisciplinary Studies (CARJIS)
ISSN (online): 2181-2454 Volume 2 | Issue 5 | May, 2022 | SJIF: 5,965 | UIF: 7,6 | ISRA: JIF 1.947 | Google Scholar |
www.carjis.org DOI: 10.24412/2181-2454-2022-5-273-280
"зубчатой", т.е. заранее непригодной к расчету.
Вышеперечисленные подходы решают задачу нахождения оптимальной толщины наиболее полно. Однако при этом может потребоваться большее количество параметров, что резко усложняет оптимизацию.
ОБСУЖДЕНИЕ
При решении инженерных задач, когда особо высокая точность не требуется, могут быть целесообразны другие подходы к отысканию оптимальной толщины, а именно: оптимальная толщину искать в классе функций, зависящих от малого, по возможности, числа параметров, и позволяющих при определенных нагружениях и граничных условиях достаточно близко подойти к оптимальной толщине. При оптимизации арки нами рассмотрены следующие классы функций. Внутренняя и внешняя поверхности арки - эллипсы (рис 1). Уравнения эллипсов в полярных координатах:
а) наружный эллипс:
_ axbx
п ^Jb2+sin2 /? (a ~bi) (3)
б) внутренний эллипс:
n h
Bi B2
Рис. 1. Внутренняя и внешняя поверхности арки - эллипсы
Central Asian Research Journal For Interdisciplinary Studies (CARJIS)
ISSN (online): 2181-2454 Volume 2 | Issue 5 | May, 2022 | SJIF: 5,965 | UIF: 7,6 | ISRA: JIF 1.947 | Google Scholar |
www.carjis.org DOI: 10.24412/2181-2454-2022-5-273-280
Из условия симметрии наружного и внутреннего эллипсов относительно средней линии имеем:
b = 2/ — b - J (5)
Толщина арки в зависимости от угла □ и параметров а2 и Ь2 выражается следующим образом:
h = 2
„ ab
R —
(6)
^ + sin2 j3( a2 -ь22)
при а2 = b2 имеем h = 2(R - а) = const.
Толщина меняется по синусоидальному закону.
h=h_0+h_1 sin npm/p_0 (7) Здесь в оптимизируемые параметры можно включить и m можно изменять в широких пределах.
Можно продолжить описание подобных функций, т.к. класс функций, позволяющих близко подойти к оптимальной толщине при различных нагружениях и граничных условиях, очевидно, довольно большой. Однако мы не ставили перед собой задачу отыскать все виды подобных функций. Вероятно, окажется полезным для каждого класса оптимизируемых конструкций подобрать соответствующие виды функций, позволяющие быстро и с минимальными затратами достаточно близко подойти к оптимальной толщине.
Описанные выше функции, определяющие толщину арки зависят от 2-х, 3-х и более оптимизируемых параметров. Ввиду того, что при выводе уравнений равновесия арки учитывалась гипотеза Кирхгоффа-Лява, на оптимизируемые параметры накладываются ограничения. Поверхности-эллипсы:
х 0,95R < а,b2 < R
Толщина меняется по синусоидальному закону.
h = h0 + ht sin^ < ОДД (9)
Р о
Как было сказано, угол раствора арки р0 также принимался за оптимизируемый параметр. На него также были наложены ограничения.
Central Asian Research Journal For Interdisciplinary Studies (CARJIS)
ISSN (online): 2181-2454 Volume 2 | Issue 5 | May, 2022 | SJIF: 5,965 | UIF: 7,6 | ISRA: JIF 1.947 | Google Scholar |
www.carjis.org DOI: 10.24412/2181-2454-2022-5-273-280
f = A (10)
Так как пролет арки L с изменением угла р0 не должен меняться, то при этом меняется радиус линии R, который зависит от угла р0 и пролета L следующим образом:
R =-(11)
2 sin
2
Таким образом, полная система ограничений на параметры с учетом (811) имеет вид:
0, 9 (12)
8 4 2 Б 2 51Пу
0 < К + /^п ^ < 0,1 *. (13)
2
Как видно из выражений (8-13), ограничения на параметры, определяющие толщину арки, зависят от значения угла ПО. Таким образом, получилась как бы "плавающая" система ограничений. Для того, чтобы избежать неприятностей, надо её "закрепить", т.е. перейти к безмерным ограничениям.
0,95 <а, Ь <1. (14)
f^0*f; (15)
0<Л0 + М \ п^<0, 1 . (16)
Ро
Здесь а, Ь, К, К, соответственно а2, Ь2, К0, К отнесённые к радиусу R.
Ограничения (14, 15) определяют в каждом случае п- мерную область D, в которой производится поиск оптимума. Поиск осложняется тем, что область D "затеняется" дополнительными ограничениями на функцию качества. Это ограничения по прочности, устойчивости, деформациям. Без учета этих
Central Asian Research Journal For Interdisciplinary Studies (CARJIS)
ISSN (online): 2181-2454 Volume 2 | Issue 5 | May, 2022 | SJIF: 5,965 | UIF: 7,6 | ISRA: JIF 1.947 | Google Scholar |
www.carjis.org DOI: 10.24412/2181-2454-2022-5-273-280
ограничений задача имела бы тривиальное решение.
Напряжения, деформации и др. внутренние усилия определяем после решения прямой задачи расчета арки по алгоритмам, описанным в [1,4,5]. Оптимизация производилась при помощи алгоритмов ГП-3 и ГП-4, описанных в [1,2, 3].
Задача 1. Оптимизация бесшарнирной арки постоянной толщины
Оптимизируемые параметры: Угол раствора арки р0 и толщина h=const.
Физические характеристики материала арки:
Модуль упругости Е=2 106 кг/см2;
Коэффициент Пуассона v= 0,3;
Допускаемые напряжения [а]= 2000 кг/см ;
Длина пролета L= 100 см.
Функция цели - объём арки V=R-h- р0 (17)
Ограничения на параметры
-<&<- 0<h <0,05. (18) 8 4
Нагрузка, равномерно распределенная, направленная к центру кривизны арки, интенсивностью q=10 кг/см.
Поиск минимумов производился с точностью в~3%. Результаты расчетов приведены в табл. 1.
Таблица 1. Результаты расчетов оптимизации арки
В V, см3 Ро , рад. h
локальному
минимуме
1 159,63 0,7854 0,011906
2 287,03 0,6729 0,018596
3 182,74 0,46404 0,008329
4 180,07 0,67131 0,011683
На рис. 2 приведены кривые нормальных напряжений а(Р), изгибающих моментов М(Р), прогибов '(Р), отвечающих двум минимумам: глобальному (1-му) и третьему.
Как видно из таблицы 1. и рис. 2, обе конструкции, соответствующие обоим минимумам, работают на пределе прочности, однако в первом минимуме
Central Asian Research Journal For Interdisciplinary Studies (CARJIS)
ISSN (online): 2181-2454 Volume 2 | Issue 5 | May, 2022 | SJIF: 5,965 | UIF: 7,6 | ISRA: JIF 1.947 | Google Scholar |
www.carjis.org DOI: 10.24412/2181-2454-2022-5-273-280
арка имеет меньший вес.
Рис.2. Кривые нормальных напряжений a(P), изгибающих моментов
М(Р), прогибов W(P): 1- глобальный минимум, 3 - третий минимум
Из рис. 2 видно, что у более пологой арки (3-й минимум) изгибающие моменты и прогиб более крутой. Отсюда можно сделать вывод, что при действии внешнего давления рациональны крутые арки постоянной толщины.
Задача 2. Оптимизация бесшарнирной арки постоянной толщины.
Интенсивность нагрузки q = 10 • sin j кг / см.
Остальные параметры те же, что и в задаче 1. Результаты расчетов приведены в табл. 2.
Таблица 2. Результаты расчетов оптимизации арка
№ V, см3 р0 , рад h
1 96,190 0,785398 0,0071743
Кривые а (Р), М(Р), представлены на рис. 3.
Рис. 3. Кривые нормальных напряжений а(Р) и изгибающих моментов
М(Р)
Central Asian Research Journal For Interdisciplinary Studies (CARJIS)
ISSN (online): 2181-2454 Volume 2 | Issue 5 | May, 2022 | SJIF: 5,965 | UIF: 7,6 | ISRA: JIF 1.947 | Google Scholar |
www.carjis.org DOI: 10.24412/2181-2454-2022-5-273-280
ЗАКЛЮЧЕНИЕ
Результаты оптимизации арок показывают, при рассмотренных системах нагрузок рациональны более крутые типы конструкции. Они обладают минимальным весом, а также меньшими значениями изгибающих моментов по сравнению с пологими оболочками [7].
REFERENCES
1. Кабулов В.К., Назиров Ш.А., Якубов С.Х. Алгоритмизация решения оптимизационных задач. - Ташкент: Изд-во «Фан» Академии наук Республики Узбекистан, 2008. - 204 с.
2. Назиров Ш.А., Якубов С.Х. Алгоритмическая система, автоматизирующая процессов оптимизации для проектирования инженерных конструкций и сооружений//Свидетельство Патентного ведомства Республики Узбекистан, DGU 01422, 2007.13.11.
3. Пискорский Л.Ф. Алгоритмы ГП2 и ГП3 поиска глобального экстремума функции многих переменных//Вопросы вычислительной и прикладной математики. Вып. 20, 1973.
4. Пискорский Л.Ф. К выбору оптимальных алгоритмов для решения классов оптимизационных задач//Известия АН УзССР, серия техн. наук. Вып.5, 1980.
5. Пискорский Л.Ф., Вопросы оптимального проектирования арок// Вопросы вычислительной и прикладной математики. Вып. 21, Ташкент, 1974.
6. Половинкин А.И., Бобков Н.К., Буш Г.Я. и др. Автоматизация поискового конструирования. - М.: Радио и связь, 1981. - 344 с.
7.Якубов С.Х. Методы и алгоритмы синтеза и анализа конструкторских и технологических решений в системе автоматизированного проектирования инженерных конструкций и сооружений. - М.: ИНФРА-М, 2019.-164 с.