Статья 'Алгоритм определения оптимальных параметров тонкостенных пространственных конструкций' - журнал 'Кибернетика и программирование' - NotaBene.ru
по
Journal Menu
> Issues > Rubrics > About journal > Authors > About the Journal > Requirements for publication > Council of Editors > Peer-review process > Policy of publication. Aims & Scope. > Article retraction > Ethics > Online First Pre-Publication > Copyright & Licensing Policy > Digital archiving policy > Open Access Policy > Article Processing Charge > Article Identification Policy > Plagiarism check policy
Journals in science databases
About the Journal

MAIN PAGE > Back to contents
Cybernetics and programming
Reference:

The algorithm for determining the optimal parameters of thin-walled spatial structures

Kolesnikov Aleksandr Georgievich

PhD in Technical Science

Associate Professor, Department of Urban, Road Construction and Structural Mechanics, Southwestern State University

305054, Russia, Kurskaya Oblast' oblast', g. Kursk, ul. 50 Let Oktyabrya, 94

ag-kolesnikov@mail.ru
Other publications by this author
 

 

DOI:

10.7256/2306-4196.2016.5.20500

Received:

22-09-2016


Published:

29-01-2017


Abstract: The subject of this study is an algorithm for determining the optimal parameters of shallow thin-walled spatial geometrically nonlinear shells on the elastic base. A main attention is paid to the selection algorithm optimization of such structures and its adaptation to the problem of finding the optimal form of shallow shell construction taking into account the material nonlinearity. The paper also demonstrates the possibility of achieving significant savings in volume (weight) of structures such as thin shallow shells on elastic base by changing its shape and thickness distribution. The optimization algorithm is based on a modification of a method of random search, including a combination of random and gradient search, as well as a method of "ravines". The algorithm for determining the optimal parameters of thin-walled spatial structures can be used to determine critical force and stresses for geometrically nonlinear shallow shells on an elastic base with a variable middle surface shape under various constraints. The novelty of the research is to use the combined method of finding the extremum of nonlinear functions with various restrictions.


Keywords:

algorithm, gradient descent, random search, geometric nonlinearity, shallow shells, optimization, elastic base, optimal shape, objective function, construction simulation

This article written in Russian. You can find original text of the article here .
Введение

В строительстве вопросы снижения стоимости несущих конструкций и повышения их эксплуатационных характеристик выходят в настоящее время на первый план.

Существенный вклад в решение этих задач вносит использование в конструкторских решениях элементов типа пологих оболочек [1-3], которые уже нашли широкое применение в строительстве, машиностроении и других областях техники (рисунок 1).

obol_1

Рисунок 1 – Оболочка на упругом основании.

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

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

Известны два основных вида задач оптимизации: задачи с ограничениями 1-го рода и 2-го рода. Под ограничениями 1-го рода будем понимать физические ограничения, наложенные на параметры проектирования: ограничения на толщину оболочки, ограничения на показатель степени ξ, который характеризуют форму оболочки, ограничения на параметр формы толщины t, на параметр k, характеризующий отношение толщины в центре оболочки и на краях [4, 12, 13].

Алгоритм определения оптимальных параметров при ограничении 1-го рода.

Необходимо решить следующие задачи оптимизации: определение формы срединной поверхности и толщины оболочки минимального объема (V(ξ, t) → min), оболочки, воспринимающей максимальную критическую нагрузк (p(ξ, t) → max), оболочки, имеющей минимальные напряжения (σ(ξ, t) → min) на всем множестве допустимых форм срединных поверхностей и постоянных толщин оболочек.

Алгоритм определения оптимальных форм конструкций типа пологих оболочек при ограничениях 1-го рода рассмотрим на примере поиска минимума объема оболочки на всем множестве допустимых форм срединных поверхностей оболочки.

В основу алгоритма определения оптимальных форм положен алгоритм определения локального минимума функции. Ввиду того, что функция объема оболочки V(xi) монотонная, вогнутая и гладкая на всей области изменения параметров x = (ξ, t) [5-6], то можно утверждать, что найденный локальный экстремум будет являться и глобальным.

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

Перед началом поиска производится нормировка, устраняющая "простые" овраги, возникающие из-за разницы в размерности параметров или разной степени зависимости функций V(xi) и p(xi) от аргументов. Перейдем к новым безразмерным параметрам:

xi=Xi/xхар, i=1,n (1)

где xi - новое пространство безразмерных параметров, Xi - старое пространство параметров, xхар - характерные размеры.

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

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

На следующем этапе поиска экстремума используется метод "наказания случайностью". Многомерный поиск сводится к последовательности одномерных. Направления для них выбираются случайным образом. Использование в ходе поиска элемента случайности делает итерационный процесс нахождения минимума более устойчивым и нечувствительным к мелкому рельефу поверхности V(xi).

Шаг поиска в методе "наказания случайностью" делается следующим образом:

∆xik+1= α∆xik, Vx<Vxk-1, xk Ω ; -∆xik+hς, Vx≥Vxk-1, xk `!in` Ω ; (2)

Здесь h – параметр, определяющий степень локализации случайного поиска,

ς - случайная величина, равномерно распределенная на отрезке (-1, 1) .

Если при очередном шаге критерий качества V(xi) уменьшился и точка не вышла за границы области Ω, следующий шаг делается в том же направлении, что и предыдущий, но в α раз больший по величине. В случае увеличения V(xi) или выхода точки за границу Ω, поиск возвращается в предыдущую точку, после чего следует "наказание случайностью": экстремум ищется на множестве точек, равномерно распределенных внутри куба со стороной 2h и центром в точке, имеющей наименьшее значение критерия. После определения точки с еще меньшим значением V(x), центр куба перемещается в эту точку. В каждой очередной серии из N проб происходит уменьшение размеров куба

hk = chk-1, (3)

где с - коэффициент локализации (0<c<1). Таким образом, случайный дрейф куба к минимуму сочетается с локализацией поиска.

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

∆xik+1 = α[xik- xik-N]h/N, i=1,n (4)

После проведения заданного числа серий из N случайных проб поиск прекращается.

Максимизация функции f(x), очевидно, эквивалентна минимизации функции - f(x). Следовательно, с помощью описанного алгоритма можно решать задачу об определении как минимума объема оболочки V(xi) и минимума значений напряжений σ(xi), так и максимума критической нагрузки p(xi) на всем множестве допустимых форм оболочек.

Алгоритм определения оптимальных форм оболочек при ограничении 2-го рода

Алгоритм определения оптимальной формы оболочки рассмотрим на примере отыскания минимума объема оболочки при ограничении на величину критической нагрузки Поставленная задача относится к весьма сложному для решения типу задач нелинейного программирования, когда целевая функция и ограничения нелинейны. Выпуклость [8] функций критической нагрузки p(xi), напряжений σ(xi) и объема V(xi) позволяет применить для ее решения алгоритмы метода штрафных функций или метод скользящего допуска [7, с. 50]. Функция, описывающая величины критических нагрузок p(xi) не может быть вычислена вне области поиска П поэтому для реше­ния задачи используется метод внутренних штрафных функций. В этом методе поиск экстремума ведется внутри области поиска. В случае нарушения ограничений штрафная функция возвращает поиск внутрь области.

Допустим, ограничения выполняются с точностью до некоторого малого εR, то есть

p(xi)≤ εR, (5)

Тогда критерий качества можно представить в виде:

V(xi) = V(x)+ `sum_(k=1)^m` LR(pk)Vk(x), (6)

где LR(VR) - штрафные коэффициенты, которые выбираются сле­дующим образом:

LR(VR) = LR, если VR>0; 0, если VR≤0. (7)

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

obol_2

Рисунок 2 – Иллюстрация нахождения оптимальной формы и толщины пологой оболочки

Алгоритм опробован для решения задач об определении формы срединной поверхности и толщины пологой оболочки минимального объема при ограничении на критическую нагрузку и величину напряжений [9-11]:

V(ξ, t)→Vmin; p(ξ, t) – p0≥0, (8)

V(ξ, t)→Vmin; σ(ξ, t) – σ0≤0, (9)

а также, оболочки, воспринимающей максимальную критическую нагрузку или напряжения в которой не превышают допустимые значения при ограниче­нии на величину объема:

p(ξ, t) →pmax; V(ξ, t) – V0≤0, (10)

σ(ξ, t)→σmin; V(ξ, t) – V0≤0. (11)

В задачах оптимизации формы оболочек минимального объема при заданной критической нагрузке, экономия объема (веса) составляет 6%-12% по сравнению с традиционной сферической формой.

В задачах оптимизации формы оболочек минимального объема при заданном значении напряжений, экономия объема (веса) составляет 4%-9%.

В задачах оптимизации формы оболочек, воспринимающих максимальную критическую нагрузку при заданной величине объема, возрастание критической нагрузки составляет 5-10%.

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

Заключение

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

Работа выполнена при поддержке гранта Президента РФ № МК-9203.2016.8

References
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
13.
Link to this article

You can simply select and copy link from below text field.


Other our sites:
Official Website of NOTA BENE / Aurora Group s.r.o.