CC BY
93
НАУЧНОЕ ИЗДАНИЕ МЕТУ ИМ. Н. Э. БАУМАНА
НАУКА и ОБРАЗОВАНИЕ
Эл № ФС77 - 48211. Государственная регистрация №0421200025. КБМ 1994-0408
электронный научно-технический журнал
Метод решения задачи оптимизации структуры дисперсно-армированных композитов при ограничениях на тепловые и прочностные свойства
# 11, ноябрь 2013 Б01: 10.7463/1113.0621065
Димитриенко Ю. И., Дроголюб А. Н., Соколов А. П., Шпакова Ю. В.
УДК 539.3
Россия, МГТУ им. Н.Э. Баумана dimit.bmstu@gmail.com adrogolub@gmail.com alsokolo@bmstu.ru ayulya@yandex.ru
Введение. Полимерные композиционные материалы находят широкое применение при создании тепловой защиты аэрокосмических летательных аппаратов и возвращаемых аппаратов. При проектировании теплозащитных материалов к ним предъявляется комплекс требований: обладания относительно высокими теплопрочностными
свойствами, пониженной теплопроводностью, пониженной плотностью и другими. В [1] был предложен аналитический подход к проектированию оптимальной микроструктуры теплозащитных композитов для обеспечения выполнения комплекса требуемых эксплуатационных характеристик материалов. В данной работе предложено развитие этого метода с учетом более детального рассмотрения микронапряжений и тепловых полей в композитах, и с применением метода конечного элемента для решения локальных задач для расчета микрохарактеристик композитов, возникающих при использовании метода асимптотического осреднения композитов с периодической структурой [2-7].
Постановка задачи оптимального проектирования композиционного материала. Рассмотрим композиционный материал, состоящий из изотропной матрицы и N-1 изотропного наполнителя со следующими характеристиками: ра - плотность; Е -модуль Юнга; V1 - коэффициент Пуассона; Оа - модуль сдвига; а^ ,а“ - пределы
прочности при растяжении, сжатии и сдвиге, соответственно; X -коэффициент теплопроводности. Здесь индекс а = 1, N — 1 соответствует наполнителю с номером а, а
индекс а=Ы соответствует матрице. Объемные концентрации компонент обозначим фа,
N
для них выполняется очевидное соотношение: '^Фа = 1.
а=1
Введем обозначения для эффективных характеристик композиционного материала: р - плотность; Ськ1 - компоненты тензора эффективных модулей упругости; ст, се -эффективные пределы прочности композита при одноосных растяжениях и сжатиях вдоль трех ортогональных осей, соответственно; Л - компоненты эффективного тензора
теплопроводности. Все эти эффективные характеристики зависят от объемных концентраций компонентов, а также от их геометрической формы. Будем далее полагать, что геометрическая форма может быть параметризована только с помощью одного числового параметра. Это допущение является достаточно сильным, ему удовлетворяет ограниченное число геометрических форм армирующих компонентов композитов, например, дисперсные частицы с шаровой формой, частицы в виде микросфер с постоянной толщиной стенок микросфер, цилиндрические волокна и некоторые другие типы, в которых меняется только один геометрический параметр. Тогда для такого типа однопараметрических композитов их эффективные характеристики могут быть представлены как функции только концентраций компонентов:
_ N _ _ __ __
РаРа ; СуЫ (р...р) ; Л (р...р); ст (<Pl,..., р) ; се (р...р). (1)
а=1
Аналитическое выражение для функций (1) при любых структурах армирования возможно только для плотности (линейная функция), для остальных функций эти зависимости имеют сложный нелинейный вид, аналитически они могут быть вычислены только приближенно ([2]), или для очень простых микроструктур , например, слоистых ([1]). Методика численного конечно-элементного вычисления функций (1) представлена в следующих разделах статьи.
Зададим ограничения на эффективные характеристики композита:
Ли < Л, ст > , се > сте, г = 1,3, (2)
где Лг, ст.т, ст.е - некоторые числовые константы - ограничения, задаваемые
проектировщиком композиционного материала.
Задачей оптимального проектирования микроструктуры композита будем называть задачу минимизации плотности композита
___ N
Р(>1,...,^г ) = 1 РаРа^ т1п (3)
а=1
при выполнении ограничений (2). Задача (3) с ограничениями (2) представляет собой задачу поиска минимума линейной функции в условиях нелинейных ограничений.
Методика решения вспомогательной задачи расчета эффективных характеристик композита. Для нахождения функций С]к1 (р...р) и Л (Р—Р) -
эффективных характеристик упругости и теплопроводности композита применим метод асимптотического осреднения [5-12], согласно которому для нахождения этих функций необходимо решить специальные краевые трехмерные задачи теории упругости и теории теплопроводности на 1/8 части ячейки периодичности композита. Эти задачи называются локальными задачами Ьрд и Ьд , соответственно.
В соответствии с [7] локальные задачи Ьрч теории упругости имеют вид:
где Ъа^ - поверхности контакта компонентов композита; - координатные плоскости, £ж ^ = 0}; £'ж = (<^ : ^ = 1 / 2} -торцевые поверхности ячейки периодичности,
имеющей форму единичного куба, £=1,2,3; - безразмерные координаты, изменяющиеся
в единичном кубе (локальные координаты), дд) - перемещения,
напряжения (микронапряжения) и деформации в компонентах композита Уа, С°‘к1 -
тензоры модулей упругости компонентов, вычисляемые через упругие константы Е , Vе1 ,
Оа изотропных компонентов стандартным образом [13], / - производная по локальной
координате ^, П] - компоненты вектора нормали. Функции ) удовлетворяют
следующим граничным условиям на Е :
(4)
=0,икрЧ) = 0,£ е£к,г*]*к*г пРиг,]*[p,q},р*ъ
0,£е£'., при г, ] = [р, д}, р * д
0,^е£,г * ] * к * г при р = д
(5)
где 8рд - заданные средние деформации (макродеформациии) композита.
Локальная задача Ьр теории теплопроводности имеет вид:
Яг(р)/,■ = 0 6 У,,
(р) =-А^; (р) 6 у, (6)
§ г (р) =%, )/г 6 Ув, ’ ()
в{Р) =6'((р), (Я«р) - А) ) Пг = 0 На ЪаМ , а = 1, Я, Р = 1,3 ,
где qi( ) - компоненты вектора теплового потока, §.( ) -компоненты вектора градиента температуры, #(“) - температура. Граничные условия для задачи Ьр теории
теплопроводности имеют вид:
в = 2 §р на Е’р, в/р = 0 на Ер
в,г = 0 на Е;, Ег, Р = 1,3; i = 1,3.
§ - заданные компоненты осредненного вектора градиента температуры.
Локальные задачи (4),(5) и (6),(7) для каждых значений р,я={1,2,3} решаем методом конечного элемента [7], в результате решения этих задач вычисляем поля микронапряжений а„(р,) и тепловых микропотоков ) для всей 1/8 ячейки
N
периодичности V = Уа. С помощью этих полей вычисляем эффективные
а=1
характеристики композита по формулам:
_ 1 N _ 1 N
с= — 2 I а„,Л-, Лр, = - —2 | С,,dV, (8)
£рд а=1 Vа §, а=1 Vа
Осуществляя описанным выше методом решение локальных задач Ьр, и Ьр при различных значениях концентраций компонентов сра, с помощью формул (8) находим искомые
функции (1) - зависимости эффективных характеристик композита от концентраций сра
его структурных компонентов.
Методика расчета эффективных констант прочности композита. Для вычисления
пределов прочности ат и ас композита используем следующую методику, предложенную в [8]. Вычислим поля компонент тензоров концентрации напряжений В./^) в компонентах композита, которые связывают микронапряжения со средними напряжениями
а„ (& ) = в;к1 & а, (9)
а=1
где а=2 а р,) . Компоненты этих тензоров вычисляются по формулам:
в;а &)=сI, & )П р,*, (10)
где Пр,к1 - компоненты тензора эффективных упругих податливостей композита, обратного к эффективному тензору модулей упругости Сг„р, (8), а поля компонент тензоров С„р, (О) - локальных модулей упругости, вычисляются с помощью соотношения
„ о)=а> |(О). (11)
5 р,
В [8] был предложен метод расчета первичного разрушения композитов на основе критериев прочности Малмейстера-Ву, а в [14] был исследован критерий прочности Писаренко-Лебедева [15]. В данной работе использован критерий прочности Малмейстера-Ву, согласно которому будем считать, что происходит разрушение в некоторой точке материала, принадлежащей компоненте а, если выполняется критерий прочности (условие разрушения)
I(а;(£)) = вааа(£) + ва(а а(О))2 + Ваъ (аиа(О))2 * 1, (12)
здесь обозначены аа - среднее напряжение и а^ -интенсивность напряжений
_ —ос . а . _а /1л\
а = а11 +а22 +аз3 , (13)
_а2 1 а а\2 . /-_а а\2 / а а\2 а 2 . а 2 а 2ч\
а = ^ ((а11 —а22) + (а11 —а33) + (а22 —а33) +6(а12 + а13 + а23 )) ,
а вг - параметры прочности материала компоненты а , определяемые по формулам
а а а а 11 1 ___
В« = ^^г_, ва = —------1_, Ваъ =-0- а = 1, N. (14)
1 а а 5 2 а а о а 2 5 3 ~_______а 2 5 4 '
ат ас ат ас 3а5 3а5
Рассмотрим процесс простого (пропорционального) нагружения в 6-ти мерном пространстве макронапряжений , при котором компоненты макронапряжений
изменяются пропорционально одному параметру х, т.е. = о„х . Если осуществляется
одноосное нагружение, то среди компонент тензора о„ имеется только одна ненулевая
компонента<а°ву, т.е. = 1(8в8„г+8„р8^)авг. Будем называть точку Ох ячейки
периодичности композита точкой первоначального разрушения, а напряжение ав = X * ав - первоначальным пределом прочности композита при нагружении по
направлению в 6-ти мерном пространстве, если для некоторого значения х* > 0 выполняется условие
I (в;к1 оа) * 1, (15)
а для значений х < X * выполняются условия /(в„к^1 (О )а*к1) < 1 для всех значений О е V .
Определение (15) дает конструктивный способ нахождения первоначального предела прочности. Для этого необходимо решить задачи Ьр, теории упругости, в соответствии с формулами (8)-(9) вычислить тензор концентрации напряжений, задать малое начальное значение параметра нагружения х, и последовательно увеличивать его, пока при некотором его значении х * в некоторой точке ячейки периодичности композита не будет выполнен критерий (15).
Для определения пределов не первоначального, а полного разрушения композита, будем использовать следующий алгоритм. Назовем трещиной непрерывную область композита, во всех точках которой выполняется критерий (15), а полным разрушением композита в целом будем называть ситуацию, когда существует по крайней мере одна трещина, соединяющая противоположные грани ячейки периодичности. Минимальное
_** . . __0
значение напряжение а ру = х**ару, приводящее к полному разрушению ячейки периодичности, будем называть эффективным пределом прочности композита при нагружении по лучу <у0ру .
Алгоритм расчета пределов прочности композита был следующим.
1) Выбираем произвольные значения макродеформаций 5р,. Решая задачи Ьрч (4),(5) для р,,=1,2,3 методом конечных элементов, по формулам (10),(11) вычислим тензор концентрации напряжений в„/(^).
2) Выбираем луч а°Ру в 6-мерном пространстве макронапряжений, по которому будет осуществляться нагружение композита. Выбираем начальное малое значение параметра нагружения х{0}, при котором критерий (15) первоначального разрушения композита заведомо не выполняется. Далее осуществляем итерационный цикл по параметру нагружения, на каждом т-м шаге последовательно увеличивая значения параметра х{т}. На каждом т-м шаге решения локальных задач Ьрч не требуется, а проверяется только выполнение условия (15). Значение х*, при котором впервые выполняется условие (15), соответствует первоначальному пределу прочности 0Ру .
3) Далее увеличиваем значение параметра нагружения х{т}. Конечные элементы, в
которых выполнен критерий (15), будем трактовать как новую компоненту композита с номером N+1 и называть разрушенной компонентой. При этом модуль упругости EN+1 разрушенной компоненты полагается существенно меньшим, чем модуль упругости всех остальных компонентов: при афЫ+1.
4) При х{т} > X* на каждом шаге итерации пересчитываем значение эффективного модуля упругости с{и по формуле (8), в которой напряжения и средние
—{т-1}
деформации а рд берутся с предыдущего шага итерации :
_, , 1 N
Ст==т* I !<т$'• (16)
8рд а=1 уа
—{т}
затем вычисляется новые средние значения деформации ард на т-м шаге
{т} _ ТГ{т} ^{т} ТГ{т} _ //^{т} \-1 ^{т} _ „,{т} — 0 {Л Н\
а рд —11 рдк1 ск1 , 11 рдк1 = (С рдк1) , ск1 =х ск1 . (17)
„ _{т} //1Ч ...
С этими значениями а рд решается локальная задача (4), (5) на т-м шаге итерации, по ее
а{ т} *-»
результатам вычисляются компоненты оц тензоров микронапряжений на т-м шаге.
5) Далее осуществляется проверка выполнения критерия (15) в каждом не разрушенном конечном элементе для шага итерации т. В случае выполнения условия (15) для дополнительных КЭ, разрушенная область композита увеличивается за счет этих КЭ и далее осуществляется переход к следующему шагу итерации.
6) Эффективный предел прочности композита к**,* будет соответствовать
ситуации, когда разрушенные конечные элементы сформируют трещину, соединяющую противоположные грани ячейки периодичности.
7) Описанную процедуру проведем для 6 различных лучей нагружения <7°^, соответствующих одноосному растяжению и сжатию вдоль 3-х координатных осей ^ , т.е. выбирая значения <7°^ равными: к°а > 0 и каа < 0, а —1,2,3. В результате получим
эффективные пределы прочности композита сат —к*а, о**а > 0, а —1,3 при
растяжении и сжатии, с ат — |, к**а < 0, а — 1,3, соответственно.
Алгоритм метода численного решения задачи оптимального проектирования.
Рассмотрим далее случай, когда имеется только одна армирующая компонента композита с концентрацией (р1 — ф. Для численного решения задачи оптимального
проектирования (1)-(3) композиционного материала с одной армирующей компонентой применим следующий алгоритм.
Задачу условной минимизации (1)-(3) представим в виде задачи безусловной минимизации с помощью метода барьерных функций [16]. С учетом ограничений (2) введем барьерную функцию
3 _ 3 ___ 3 _
Л(ф) — I тах{0, Яи (ф) Я} + тах{0, с^ скт (ф)} + тах{0, с с скс (ф)}. (18)
I—1 I—1 I—1
Тогда задача (1)-(3) сводится к задаче безусловной минимизации
P(r) = inf(p(^) + rR(p)) ^ min, (19)
ф
которая, в свою очередь, сводится к последовательности следующих задач минимизации
фk+1) = argmin(p(^) + rkR(p)), rk = Prk_1, в = const > 1 (20)
с критерием остановки
rkR(<Pk+i)<£ (21)
для заранее определенного s.
Барьерная функция (18) имеет нелинейный характер из-за присутствия функций
Git (ф), <7ю (ф), для решения каждой из задач (20) будем пользоваться методом Хука-
Дживса. В качестве начальной точки для метода барьерных функций требуется точка,
заведомо не удовлетворяющая ограничениям. В качестве такой точки выберем решение
задачи (3) без ограничений (2): ф(0) = argminр(ф). В случае, если ф(0) удовлетворяет
ограничениям (2), то, очевидно, является точным решением исходной задачи (1)-(3). В
(0)
противном случае ф можно использовать в качестве начальной точки процесса оптимизации. Таким образом, решение последовательности задач (13) с начальной точкой ф^0) и критерием остановки (21) является решением исходной задачи (1)-(3).
Итерационную процедуру вычисления эффективных функций (2), описанную в предыдущих разделах, основанную на применении методов МКЭ, используем на каждом шаге оптимизации структуры.
Результаты численного моделирования. Численные расчеты проводились для дисперсно-армированного композита с изотропной матрицей и одним изотропным наполнителем, представляющим собой тонкие микросферы (рис. 1).
Рис. 1. Изображение 1/8 ячейки периодичности дисперсно-армированного композита на
основе микросфер.
В процессе расчета для каждой концентрации наполнителя ячейка периодичности разбивалась приблизительно на 400 тысяч конечных элементов. Разбиение проводилось с помощью свободно распространяемого программного продукта Netgen 4.4. Локальные задачи для поиска эффективных характеристик решались в программном комплексе GCD, разрабатываемом на кафедре ФН-11 МГТУ им. Н.Э. Баумана.
В качестве материала тонкой оболочки микросфер было использовано стекло с характеристиками: р1=2600 кг/м3, Е1=70 ГПа, v1=0,2, А1 =0,7 Вт/(мК), а1т — 0,5 ГПа ,
сС — 0,6 ГПа , с1 — 0,4 ГПа; в качестве матрицы рассматривалось кремнийорганическое
2 3 2 2
связующее со следующими характеристиками: р =1100 кг/м, Е =2 ГПа, V =0,35,
Х2=0,2 Вт/(мК), — 0,01 ГПа, а2с — 0,02 ГПа , <у25 — 0,005 ГПа . Вследствие симметрии
структуры рассмотренного типа композита относительно всех 3-х координатных направлений, он будет обладать изотропными эффективными свойствами, т.е. все эффективные модули упругости Е (ф) композита совпадают, совпадают также все 3
коэффициента теплопроводности Я1 (ф) — Я (ф), а также пределы прочности с т (ф) — ст (ф) .
В процессе численного моделирования были рассчитаны зависимости эффективных характеристик от концентрации стеклянной компоненты микросфер, эти зависимости показаны на рисунках 2-5. В процессе численного моделирования были рассчитаны зависимости эффективных характеристик от концентрации стеклянной компоненты микросфер, эти зависимости показаны на рисунках 2-5. Графики зависимости модуля упругости и коэффициента теплопроводности от коэффициента армирования являются монотонными, как это и следует из теоретических соображений [3]. Коэффициент Пуассона у(ф) композита меняется немонотонным образом - до значения концентрации ф — 0.27 он уменьшается, а затем снова возрастает, однако эти изменения невелики по абсолютной величине (от 0.35 до 0.337), вцелом, значения коэффициента Пуассона композита близки по своим значениям к коэффициенту Пуассона матрицы. График зависимости прочности на растяжение ст (ф) не является монотонным: при увеличении концентрации микросфер
прочность композита вначале резко падает почти в 2 раза по сравнению с прочностью неармированной матрицы, далее при увеличении содержания микросфер она немного возрастает (примерно на 15%), а затем при предельном содержании микросфер - снова падает. Небольшие осцилляции прочности, показанные на рисунке 5, обусловлены численной погрешностью вычислений прочности, которая вызвана сеточными эффектами (при различных концентрациях микросфер конечно-элементная сетка была различной и это вносило определенную погрешность в результат расчета прочности).
1300 I--------------------1-------------------1--------------------1-------------------1-------------------1-------------------1-------------------1--------------------1-------------------1-------------------
О 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
Ф
Рисунок 2. Зависимость эффективного модуля упругости Ei (ф) дисперсно-армированного
композита от концентрации микросфер ф.
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
Ф
Рис. 3. Зависимость эффективного коэффициента Пуассона дисперсно-армированного
композита от концентрации микросфер ф .
1(4)
Ф
Рис. 4. Зависимость эффективного коэффициента теплопроводности Я (ф) дисперсно-армированного композита от концентрации микросфер
Рис. 5. Зависимость прочности на растяжение ат (ф) дисперсно-армированного композита
от концентрации микросфер
Для задачи оптимизации были выбраны следующие ограничения (2): Д = 0,3Вт / (м ■ К) и (гт = ас = 5,6МПа.
При этом значение оптимума было достигнуто при концентрации стеклянной фазы ф = 0,084, что соответствует концентрации микросфер ф0 = 0,31. Соответствующее
оптимальное значение плотности композита составляет р = 977кг/м3.
Выводы. Разработан численный метод решения задачи оптимизации микроструктуры композиционного материала, с одновременными ограничениями на теплофизические и механические характеристики. Метод основан на конечно-элементном расчете эффективных характеристик композита: модулей упругости, коэффициентов теплопроводности и прочностей, а также на итерационном решении задач оптимизации функционала плотности. Предложенный метод реализован на примере дисперсно-армированного композиционного материала с наполнителем в виде стеклянных микросфер. Численно рассчитаны зависимости эффективных характеристик композита от коэффициента армирования, а также найдено оптимальное содержание наполнителя композита. Предложенный метод может использоваться для определения оптимальных структур композитов в задачах проектирования тепловой защиты.
Исследование выполнено при поддержке Министерства образования и науки Российской Федерации (Соглашения Xq14.B37.21.0448, 14.B37.21.1869, государственные контракты №14.514.11.4052, 14.516.11.0083, НИР № 1.5433.2011), грантов Президента РФ (МК-6421-2012-9, МК-765.2012.8, МК-6573.2013.3) и гранта РФФИ№12-08-00998.
Список литературы
1. Yu.I. Dimitrienko, G.A. Efremov, S.A. Chernyavsky. Optimum design of erosion-stable heatshield composite materials // Applied Composite Materials. 1997. № 4. 35-52.
2. Бахвалов Н.С., Панасенко Г.П. Осреднение процессов в периодических средах. М.: Наука, 1984. 352 с.
3. Победря Б.Е. Механика композиционных материалов. М.: МГУ, 1984. 336 с.
4. Санчес-Паленсия Э. Неоднородные среды и теория колебаний: пер. с англ. М.: Мир, 1984. 336 с.
5. Димитриенко Ю.И., Соколов А.П. Разработка автоматизированной технологии вычисления эффективных упругих характеристик композитов методом асимптотического осреднения // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2008. № 2. С. 57-67.
6. Димитриенко Ю.И., Соколов А.П. Автоматизация прогнозирования свойств композиционных материалов на основе метода асимптотического осреднения // Информационные технологии. 2008. № 8. С. 31-38.
7. Димитриенко Ю.И., Кашкаров А.И. Конечно-элементный метод для вычисления эффективных характеристик пространственно-армированных композитов // Вестник МГТУ им. Баумана. Сер. Естественные науки. 2002. № 2. С. 95-108.
8. Димитриенко Ю.И., Соколов А.П., Шпакова Ю.В., Юрин Ю.В. Моделирование поверхностей прочности композитов на основе микроструктурного конечно-элементного анализа // Наука и образование. МГТУ им. Н.Э. Баумана. Электрон. журн. 2012. № 11. DOI: http://dx.doi.org/10.7463/1112.0496336.
9. Димитриенко Ю.И., Соколов А.П. Многомасштабное моделирование упругих композиционных материалов// Математическое моделирование.- 2012. - Т. 24. - № 5. -С. 3-20.
10. Димитриенко Ю.И., Соколов А.П. Современный численный анализ механических свойств композиционных материалов // Известия РАН. Физическая серия - Т. 75, № 11. -
2011. - с. 1551-1556.
11. Димитриенко Ю.И., Морозов А.Н., Соколов А.П., Ничеговский Е.С. Моделирование эффективных пьезоэлектроупругих свойств композиционных материалов// Вестник МГТУ им.Н.Э. Баумана. Сер. Естественные науки.-№ 3. -2010.- с. 86-97.
12. Димитриенко Ю.И., Ничеговский Е.С. Численное моделирование магнитных свойств композиционных материалов// Вестник МГТУ им.Н.Э. Баумана. Сер. Естественные науки.-№ 1. -2010.- с. 3-12.
13. Димитриенко Ю.И. Основы механики твердого тела/ Механика сплошной среды.Т.4.-Изд-во МГТУ им.Н.Э. Баумана.-2013.-580 с.
14. Ю.И. Димитриенко, А.П. Соколов, С.В. Сборщиков Моделирование микроразрушения тканевых композитов// Вестник МГТУ им.Н.Э. Баумана. Сер. Естественные науки. Спец.выпуск № 3 «Математическое моделирование».-2012.-с. 5-19.
15. Лебедев А.А. Развитие теорий прочности в механике материалов // Проблемы прочности. 2010. № 5. С. 127-146.
16. Банди Б. Методы оптимизации. Вводный курс. М.: Радио и связь, 1988.
SCIENTIFIC PERIODICAL OF THE BAUMAN MSTU
SCIENCE and EDUCATION
EL № FS77 - 48211. №0421200025. ISSN 1994-0408
electronic scientific and technical journal
Method for solving an optimization problem of the structure of fiber reinforced composites under restrictions on thermal and strength properties
# 11, November 2013
DOI: 10.7463/1113.0621065
Dimitrienko Yu.I., Drogolyub A.N., Sokolov A.P., Shpakova Yu.V.
Bauman Moscow State Technical University, 105005, Moscow, Russian Federation
dimit.bmstu@gmail.com
adrogolub@gmail.com
alsokolo@bmstu.ru
ayulya@yandex.ru
In this paper the authors formulated a problem of optimising the composite microstructure, which is a problem of minimisation of the composite density under given restrictions on the components of effective elasticity and heat-conduction tensors, and effective limits of the composite strength. The authors proposed a method for solving an optimal design problem of a composite material reinforced with micro-spheres and designed for thermal shielding. The method of asymptotical averaging (homogenisation method) was used for calculating effective composite characteristics such as elastic and heat-conduction modulus tensors and strength limits; the finite element method was used for solving local 3D heat conduction and elasticity theory problems on the composite’s periodicity cells. The finite element method was implemented with the use of GCD software developed at “FN-11” department of Bauman MSTU. The Hooke Jeeves algorithm was used for solving optimisation problems. Numerical implementation of the proposed method of micro-structure optimisation was developed by the example of a composite based on glass micro-spheres.
Publications with keywords: strength, elasticity, composites, finite element method, optimization methods, optimization, thermal conductivity, asymptotic averaging method, micro-sphere
Publications with words: strength, elasticity, composites, finite element method, optimization methods, optimization, thermal conductivity, asymptotic averaging method, micro-sphere
References
1. Dimitrienko Yu.I., Efremov G.A., Chernyavsky S.A. Optimum design of erosion-stable heatshield composite materials. Applied Composite Materials, 1997, no. 4, pp. 35-52. DOI: 10.1007/BF02481387
2. Bakhvalov N.S., Panasenko G.P. Osrednenieprotsessov vperiodicheskikh sredakh [Averaging processes in periodic media]. Moscow, Nauka, 1984. 352 p.
3. Pobedria B.E. Mekhanika kompozitsionnykh materialov [Mechanics of composite materials]. Moscow, MSU Publ., 1984. 336 p.
4. Sanchez-Palencia E. Non-Homogeneous Media and Vibration Theory. Springer Verlag, Berlin, 1980. (Ser. Lecture Notes in Physics, vol. 127.). DOI: 10.1007/3-540-10000-8 (Russ. ed.: Sanchez-Palencia E. Neodnorodnye sredy i teoriia kolebanii. Moscow, Mir, 1984. 471 p.).
5. Dimitrienko Iu.I., Sokolov A.P. Razrabotka avtomatizirovannoi tekhnologii vychisleniia effektivnykh uprugikh kharakteristik kompozitov metodom asimptoticheskogo osredneniia [Development of automated technology of calculation of effective elastic characteristics of composites by method of asymptotic averaging]. VestnikMGTUim. N.E.Baumana. Ser. Estestvennye nauki [Herald of the Bauman MSTU. Ser. Natural science], 2008, no. 2, pp. 57-67.
6. Dimitrienko Iu.I., Sokolov A.P. Avtomatizatsiia prognozirovaniia svoistv kompozitsionnykh materialov na osnove metoda asimptoticheskogo osredneniia [Automation of predicting the properties of composite materials on basis of the method of asymptotic averaging]. Informatsionnye tekhnologii, 2008, no. 8, pp. 31-38.
7. Dimitrienko Iu.I., Kashkarov A.I. Konechno-elementnyi metod dlia vychisleniia
effektivnykh kharakteristik prostranstvenno-armirovannykh kompozitov [Finite-element method to calculate the effective characteristics of spatially reinforced composites]. Vestnik MGTU im. N.E.Baumana. Ser. Estestvennye nauki [Herald of the Bauman MSTU. Ser. Natural science], 2002, no. 2, pp. 95-108.
8. Dimitrienko Ju.I., Sokolov A.P., Shpakova Ju.V., Jurin Ju.V. Modelirovanie poverhnostej
prochnosti kompozitov na osnove mikrostrukturnogo konechno-jelementnogo analiza [Modeling of surfaces of composite strength on the basis of the microstructural finite element analysis]. Nauka i obrazovanie MGTU im. N.E. Baumana [Science and Education of the Bauman MSTU],
2012, no. 11. DOI: 10.7463/1112.0496336
9. Dimitrienko Ju.I., Sokolov A.P. Mnogomasshtabnoe modelirovanie uprugih
kompozicionnyh [Multiscale modeling of elastic composite materials]. Matematicheskoe modelirovanie, 2012, vol. 24, no. 5, pp. 3-20.
10. Dimitrienko Ju.I., Sokolov A.P. Sovremennyj chislennyj analiz mehanicheskih svojstv kompozicionnyh materialov [Modern numerical analysis of mechanical properties of composite materials]. IzvestijaRAN. Ser. Fizicheskaja, 2011, vol. 75, no. 11, pp. 1551-1556.
11. Dimitrienko Ju.I., Morozov A.N., Sokolov A.P., Nichegovskij E.S. Modelirovanie jeffektivnyh p'ezojelektrouprugih svojstv kompozicionnyh materialov [Simulation of Efficient Piezoelectric-Elastic Composite Materials]. Vestnik MGTU im. N.E.Baumana. Ser. Estestvennye nauki [Herald of the Bauman MSTU. Ser. Natural science], 2010, no. 3, pp. 86-97.
12. Dimitrienko Ju.I., Nichegovskij E.S. Chislennoe modelirovanie magnitnyh svojstv
kompozicionnyh materialov [Numerical Simulation of Magnetic Properties of Composites]. Vestnik MGTU im. N.E.Baumana. Ser. Estestvennye nauki [Herald of the Bauman MSTU. Ser. Natural science], 2010, no. 1, pp. 3-12.
13. Dimitrienko Ju.I. Mehanika sploshnoj sredy. V 4 t. T. 4. Osnovy mehaniki tverdyh sred [Continuum Mechanics. 4 vols. Vol. 4. Fundamentals of mechanics of solid media]. Moscow, Bauman MSTU Publ., 2013. З80 з.
14. Dimitrienko Ju.I., Sokolov A.P., Sborshhikov S.V. Modelirovanie mikro-razrushenija tkanevyh kompozitov [Simulation of Microfracture of Textile Composites]. Inzhenernyy zhurnal: nauka i innovatsii [Engineering Journal: Science and Innovation], 2012, no. 2. Available at; http;//engjournal.ru/catalog/mathmodel/material/32.html , accessed 01.10.2013.
1З. Lebedev A.A. Razvitie teorij prochnosti v mehanike materialov [Development of the theories of strength in the mechanics of materials]. Problemyprochnosti, 2010, no. З, pp. 127146. (English Translation; Strength of Materials, October 2010, vol. 42, iss. З, pp. З78-З92. DOI; 10.1007/s11223-010-9246-9 ).
16. Bunday B. Basic optimization methods. London, Edward Arnold, 1984. (Russ. ed.; Bunday B. Metody optimizacii. Vvodnyj kurs. Moscow, Radio i svjaz', 1988. 128 p.).
