О решении обратных задач численной оптимизации композитных конструкций

О РЕШЕНИИ ОБРАТНЫХ ЗАДАЧ ЧИСЛЕННОЙ ОПТИМИЗАЦИИ КОМПОЗИТНЫХ КОНСТРУКЦИЙ

Е, В, Амелина1, С, К, Голушко1,2

1 Новосибирский государственный университет, 630090, Новосибирск 2Институт вычислительных технологий СО РАН, 630090, Новосибирск

УДК 539.3

Б01: 10.24411/9999-016А-2019-10005

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

Введение

Металлокомпозитные баки высокого давления являются перспективным типом конструкций ракетно-космической техники, обеспечивая их повышенную прочность и минимальный вес. Внутреннее давление является основной нагрузкой для таких баков. Анализ возможностей применения гибридных баков, состоящих из металлических сосудов (лейнеров), обеспечивающих герметичность конструкции, и композитных силовых оболочек, обеспечивающих прочность и жесткость, является весьма актуальным и перспективным направлением исследований [1,2].

К настоящему времени можно выделить несколько подходов к решению задач оптимального и рационального проектирования композитных конструкций [3]- [11]. Наиболее распространенным является подход, при котором используются уравнения безмоментной теории оболочек и "нитяная" модель композиционного материала (КМ) [12,13]. При таком подходе реализация безмоментного состояния конструкции по существу не обеспечивается, а только предполагается. Вопрос о том, будет ли спроектированная оболочка действительно безмоментной при безмоментных краевых условиях остается в таких работах открытым.

Другой подход связан с обоснованным упрощением исходной постановки задачи трехмерной теории упругости и использованием различных вариантов двумерных теорий оболочек [3,14,15].

Именно такой подход использован в данной работе. Для проверки достоверности получаемых решений использованы варианты теории оболочек Тимошенко [16] и Андреева-Немировского [17], учитывающие поперечные сдвиги с различной степенью точности. Достоверность результатов численной оптимизации обеспечена сравнением с известными частными аналитическими решениями.

1 Постановка задачи

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

Работа выполнена при частичной финансовой поддержке Российского фонда фундаментальных исследований (код проекта 18-29-18029).

ISBN 978-5-901548-42-4

V > У0, Р > Ро, М < М0,

где V — внутренний объем, Р — внутреннее давление, М — масса сосуда; Ро> Мо — некоторые заданные значения.

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

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

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

Оболочка задается вращением образующей г = г (в) относительно оси оу (рис. 1), где г — текущий радиус оболочки, в — угол между нормалью к поверхности оболочки и осью вращения, который изменяется в пределах [0о; в\].

Для решения прямых задач расчета напряженно-деформированного состояния композитных баков и верификации решений задач оптимизации будем использовать оболочек Кирхгофа — Лява [18] и уточненные теории Тимошенко [16] и Андреева Физические соотношения опишем с помощью структурных моделей КМ [3].

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

М = 2-к rRihdd [рт

Je о

где рт, Рг — плотности материалов связующего и армирующих волокон, шг — объемное содержание арматуры.

В качестве управляющих функций выберем: главный радиус кривизны R\(9), толщину оболочки h(ß), угол спирального армирования ф(в).

Искомое решение должно удовлетворять ограничению на внутренний объем оболочки

W r2Ri sin OdO = V0

J6o

и требованию ее неразрушения

max|6sr,bsm} < 1,

где bsr, bsm — интенсивности напряжений в связующем материале и армирующих волокнах, определяемым по рассчитанному НДС бака по одной из структурных моделей КМ.

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

0 < ф < 90o, h0 < h < hR*0 < R1 < R¡,

(НДС) многослойных классическую теорию — Немировского [17].

где Л.0, Д| — константы.

При изготовлении композитных оболочек вращения широкое распространение получил метод непрерывной намотки по геодезическим линиям. В этом случае справедлива формула Клеро:

г sin ф(г) = С,

где С — константа, определяемая из условия на экваторе оболочке, а толщина оболочки определяется соотношением

R cos

h(r) = hR-

V cos ф(г)'

которое имеет особенность у кромки отверстия, где угол армирования должен быть равен 90°. Формула (7) применяется при г > ri+rw, где величина гш равна ширине армирующей лепты. Соотношение, определяющее толщину стопки бака в итоге принимает вид:

h

h

R cos фк

гш cos ф(г! + гш) ' R cos фи

г cos ф(г )'

г < т'1 + гш; г > Г1 + гш.

h

2 Прямая задача расчета НДС многослойного бака

Проблема определения НДС композитного бака приводит к необходимости решения краевых задач для жестких систем дифференциальных уравнений. Такие системы являются плохо обусловленными, а в структуре их решений имеются решения с ярко выраженным характером погранслоев. Численные расчеты были выполнены методами сплайн-коллокации и дискретной ортогонализации. реализованными в пакетах СОЬБУБ [19] и (тМ 1)0 [20] [24] соответственно. Эти вычислительные инструменты показали высокую точность и эффективность при решении задач механики композитных пластин и оболочек [3].

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

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

3 Обратная задача оптимизации многослойного бака

Рассмотрим тестовую задачу совместной оптимизации формы и толщины безмоментной осесимметричной изотропной оболочки вращения минимального веса, которая имеет аналитическое решение [29]. Здесь и далее параметры оболочки: У0 = Зм3, г\ = 0.5м, г2 = 1м, д = 108Па. В качестве ограничения по прочности использовано соотношение

где а* — постоянная прочности материала.

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

Тог + пд(г2 - г2)у/Г+^ г~—2 Т\гг"

±1 =---, ±2 = дгу 1 + г'2 + ——2

2

2пг 1+ г

Рассмотрим случай, когда ац = а*. Соотношения для радиуса и толщины оптимальной оболочки в декартовых координатах примут вид [29]:

г

где а, @ зависят от геометрических параметров конструкции.

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

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

Для решения вспомогательных задач безусловной оптимизации проводилась дискретизация управляющих функций на сетке, не связанной с сеткой решения прямой задачи, и, таким образом, искалось решение невыпуклой задачи конечномерной оптимизации [26]. Для решения задачи безусловной оптимизации был использован программный комплекс 0РТС01М-А [27,28].

На рис. 3 приведено сравнение полученных численного и аналитического решений при различных значения параметра а (сплошные кривые — аналитическое решение, пунктирные — численное). В табл. 1 приведены числовые значения погрешностей.

Отличие между расчетным и точным значением массы составило 2 %, т.е. чувствительность функционала массы оболочки к изменению главного радиуса кривизны невелика.

Другая тестовая задача заключалась в поиске проекта бака минимальной массы в рамках полной момент-ной постановки на основе уравнений теории оболочек Кирхгофа - Лява. Значения растягивающих усилий на краях оболочки являлись искомыми параметрами, а перерезывающие усилия и моменты равнялись нулю. Максимальные относительные погрешности по компонентам искомых параметров составили для во = 2.З%, вг = 0.8% к = 6.З%, Кг = 2.9% Тю = 4.0% То2 = 1.2%. Отклонение по значению целевой функции (весу) составило 0.4%.

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

Рис. 3: Влияние параметра геометрии а па вид оптимальных решений

Таблица 1: Максимальные относительные погрешности численных решений

Фупкция Максимальная относительная разность

проектирования а = 1 а = 2 а = 3

h 0.003% 1°/с 2%

Ri 0.003% 4% 17%

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

В таб. 2 приведено сравнение проектов баков, полученных при решении оптимизационных задач. Приняты следующие обозначения: 1 — Ф , 2 — ^,3 — Дх; после 'Т" идут номера функций управления, которые задавались как функции, после "С" функции управления константы. Сокращением "Сеос!" обозначен проект с непрерывной намоткой по геодезическим линиям.

Таблица 2: Сравнение различных проектов многослойных баков

Проект Масса, кг М/Мс 123 ^max ^min hmax /hmin Rimax/Rimin

F123 16.04 72.5% 4 1.08 1.62

F12C3 16.02 72.4% 12 1.55 1.00

F13C2 16.05 72.5% 2 1.00 1.94

F23C1 17.01 76.9% 0 4.74 3.14

F1C23 20.23 91.4% 15 1.00 1.00

F2C13 16.42 74.2% 0 1.83 1.00

F3C12 16.07 72.6% 0 1.00 1.83

Good 19.09 86.2% 85 9.95 5.42

С123 22.13 100.0% 0 1.00 1.00

Масса оптимизированного бака, полученного непрерывной намоткой волокон по геодезическим линиям (Good) составила 19 кг. однако за счет выбора оптимальной формы бака (проект F3C12) можно добиться снижения веса до 16 кг.

Отметим, что проект с геодезической непрорывной намоткой обладает перепадом толщины стенки бака почти в 10 раз с резким градиентом вблизи места закрепления. Проект F123 более предпочтителен, перепад значений параметров не превышает 10% по толщине и 4.5% по углу армирования. Представляют интерес проекты с постоянной толщиной стопки бака (F13C2. F3C12). показывающие возможность получения конструкций с массой близкой к минимальной за счет управляющих функций Дх и ф.

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

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

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

Заключение

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

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

Список литературы

fl] Thesken J. С., Murthy P. L. N., Phoenix S. L. et al. A theoretical investigation of composite overwrapped pressure vessel (COPV) mechanics applied to NASA full scale tests. NASA/TM-2009^215684, 2009.

[2] Амелина E. В., Буров A. E., Голушко С. К., Лепихин А. М., Москвичев В. В., Юрченко А. В. Расчетно-экспериментальная оценка прочности металлокомпозитного бака высокого давления // Вычислительные технологии. 2016. Т. 21, № 5. С. 3-21.

[3] Голушко С. К., Немировский Ю. В. Прямые и обратные задачи механики упругих композитных пластин и оболочек вращения. Москва: ФИЗМАТЛИТ. 2008.

[4] Голушко С. К., Немировский Ю. В. Обзор и анализ подходов к проблеме рационального проектирования армированных оболочек. Красноярск. 1988. — 31 с. (Препринт ВЦ СО АН СССР; № 16).

[5] Голушко С. К., Немировский Ю. В. Об одном подходе к проектированию осесимметричных оболочек с равнонапряженной арматурой / / Численные методы решения задач теории упругости и пластичности / Материалы X Всес. конф. Новосибирск: ИТПМ СО РАН, 1988. С. 58-64.

[6] Голушко С. К., Немировский Ю. В. Построение проектов армированных обол очечных конструкций минимального веса // Вычислительные проблемы механики / Межвуз. сб. Красноярск: Изд-во: КГУ, 1989. С. 117-130.

[7] Голушко С. К. Две задачи рационального проектирования армированных оболочек вращения. Красноярск. 1986. С. 8-14. (Препринт ВЦ СО АН СССР; № 1).

[8] Голушко С.К. Проектирование тонкостенных оболочек с равнонапряженной арматурой. Красноярск. 1987. С. 29-32. (Препринт ВЦ СО АН СССР; № 1).

[9] Golushko S.K. Direct and inverse problems in mechanics of composite shells // Proc. of The Sixth Japan-Russia Joint Symposium on Computational Fluid Dynamics. Nagoya: Nagoya University, 1998. P. 125-130.

[10] Голушко С.К., Немировский Ю.В., Одновал С.В. Расчет и рациональное проектирование композитных оболочек вращения // Динамика сплошной среды. 1998. № 113. С. 39-44.

fll] Golushko S. К. Direct and inverse problems in the mechanics of composite plates and shells // Computational Science and High Performance Computing. Russian-German Advanced Research Workshop, Novosibirsk, Russia. Germany: Springer-Verlag, 2005. P. 205-228.

[12] Stadler W., Krishnan V. Natural structural shapes for shells of revolution in the membrane theory of shells // Structural Optimization. 1989. V. 1. P. 19-27.

[13] Vasiliev V. V. Composite pressure vessels: design, analysis, and manufacturing. Blacksburg, Virginia, USA: Bull Ridge Publishing. 2009.

[14] Liang C.C., Chen H. W., Wang С. H. Optimum design of dome contour for filament-wound composite pressure vessels based on a shape factor // Compos. Struct. 2002. V. 58, iss. 4. P. 469-482.

[15] Zu L., Koussios S., Beukers A. Shape optimization of filament wound articulated pressure vessels based on non-geodesic trajectories // Compos. Struct. 2010. V. 92, iss. 2. P. 339-346.

[16] Григорепко Я. M., Василенко А. Т. Задачи статики анизотропных неоднородных оболочек. Москва: Наука. 1992.

[17] Андреев А. Н., Немировский Ю. В. Многослойные анизотропные оболочки и пластины: Изгиб, устойчивость, колебания. Новосибирск: Наука. 2001.

[18] Новожилов В. В. Теория тонких оболочек. Л.: Судпромгиз. 1951.

[19] Ascher U., Christiansen J., Russel R. D. Collocation software for boundary value ODEs // ACM. Trans, on Math. Software. 1981. V. 7, iss. 2. P. 209-222.

[20] Голушко С. К., Горшков В. В., Юрченко А. В. О двух численных методах решения многоточечных нелинейных краевых задач // Вычислительные технологии. 2002. Т. 7, № 2. С. 24-34.

[21] Голушко С. К., Морозова Е. В., Юрченко А. В. О численном решении краевых задач для жестких систем дифференциальных уравнений // Вестник КазНУ. Серия: Математика, механика, информатика. 2005. № 2. С. 12-26.

[22] Голушко С. К., Горшков В. В., Юрченко А. В. О двух численных методах расчета сопряженных композитных конструкций //В кн.: Современные методы математического моделирования природных и антропогенных катастроф. Тезисы докладов V науч. конференции. 1999. С. 49-55.

[23] Golushko S. К., Yurchenko А. V. Solution of boundary value problems in mechanics of composite plates and shells // Russian Journal of Numerical Analysis and Mathematical Modelling. 2010. V. 25, iss. 1. P. 27-55.

[24] Голушко С. К., Идимешев С. В., Семисалов Б. В. Методы решения краевых задач механики композитных пластин и оболочек. Новосибирск: КТИ ВТ СО РАН, 2014. 131 с.

[25] Бертсекас Д. Условная оптимизация и метод множителей Лагранжа. Москва: Радио и связь. 1987.

[26] Евтушенко Ю. Г. Методы решения экстремальных задач и их применение в системах оптимизации. Москва: Наука. 1982.

[27] Горнов А. Ю. Вычислительные технологии решения задач оптимального управления. Новосибирск: Наука. 2009.

[28] Амелина Е. В., Голушко С. К., Горнов А. Ю., Юрченко А. В. Задачи проектирования многослойных гибридных сосудов давления //В сб.: Безопасность и живучесть технических систем. Материалы и доклады: в 3-х томах. 2015. С. 34-40.

[29] Баничук Н. В. Оптимизация осесимметричных мембранных оболочек // Прикладная математика и механика. 2007. № 1. Р. 19-27.

Амелина Евгения Валерьевна — к.ф.-м.н., ст. науч. сотр. Новосибирского государственного университета;

e-mail: amelina.evgenia@gmail.com;

Голушко Сергей Кузьмич — д.ф.-м.н., профессор Новосибирского государственного университета;

глав. науч. сотр. Института вычислительных технологий СО РАН;

e-mail: s.k.golushko@gmail.com.

Дата поступления — 30 апреля 2019 г.