CC BY
45
ВИАМ/2013-Тр-01-03
УДК 678.027.7
К ВОПРОСУ УДАЛЕНИЯ ИЗЛИШКОВ СВЯЗУЮЩЕГО ПРИ АВТОКЛАВНОМ ФОРМОВАНИИ ИЗДЕЛИЙ ИЗ ПОЛИМЕРНЫХ КОМПОЗИЦИОННЫХ МАТЕРИАЛОВ
М. И. Душин
кандидат технических наук
А.В. Хрульков
А.Е. Раскутин кандидат технических наук
Январь 2013
Всероссийский институт авиационных материалов (ФГУП «ВИАМ» ГНЦ) - крупнейшее российское государственное материаловедческое предприятие, на протяжении 80 лет разрабатывающее и производящее материалы, определяющие облик современной авиационно-космической техники. 1700 сотрудников ВИАМ трудятся в более чем тридцати научно-исследовательских лабораториях, отделах, производственных цехах и испытательном центре, а также в четырех филиалах института. ВИАМ выполняет заказы на разработку и поставку металлических и неметаллических материалов, покрытий, технологических процессов и оборудования, методов защиты от коррозии, а также средств контроля исходных продуктов, полуфабрикатов и изделий на их основе. Работы ведутся как по государственным программам РФ, так и по заказам ведущих предприятий авиационнокосмического комплекса России и мира.
В 1994 г. ВИАМ присвоен статус Государственного научного центра РФ, многократно затем им подтвержденный.
За разработку и создание материалов для авиационно-космической и других видов специальной техники 233 сотрудникам ВИАМ присуждены звания лауреатов различных государственных премий. Изобретения ВИАМ отмечены наградами на выставках и международных салонах в Женеве и Брюсселе. ВИАМ награжден 4 золотыми, 9 серебряными и 3 бронзовыми медалями, получено 15 дипломов.
Возглавляет институт лауреат государственных премий СССР и РФ, академик РАН, профессор Е.Н. Каблов.
Статья подготовлена для опубликования в научно-техническом сборнике «Авиационные материалы и технологии», №2, 2005 г.
УДК 678.027.7
М.И. Душин, А.В. Хрульков, А.Е. Раскутин
К ВОПРОСУ УДАЛЕНИЯ ИЗЛИШКОВ СВЯЗУЮЩЕГО ПРИ АВТОКЛАВНОМ ФОРМОВАНИИ ИЗДЕЛИЙ ИЗ ПОЛИМЕРНЫХ КОМПОЗИЦИОННЫХ МАТЕРИАЛОВ
Описан подход к определению времени удаления излишков связующего из пакета заготовки препрегов изделий из ПКМ на основе закона Дарси с предварительным определением экспериментальным путем коэффициента проницаемости армирующего наполнителя, а также вязкости и жизнеспособности связующего при температуре формования в автоклаве. Для определения коэффициента проницаемости использован расчетный метод сочетания теории протекания (пер-коляции) и приведения к элементарной ячейке. Представлены результаты расчета проницаемости пористой системы, состоящей из волокон диаметром 7 мкм, в зависимости от пористости структуры, а также экспериментальные данные по изменению вязкости и жизнеспособности эпоксидного связующего в зависимости от температуры.
Ключевые слова: автоклавное формование, содержание связующего, проницаемость, давление, продолжительность.
M.I. Dushin, A.V. Khrulkov, A.E. Raskutin
ABOUT REMOVAL RESIN EXCESS UNDER AUTOCLAVE FORMING OF PRODUCTS FROM POLYMER COMPOSITE MATERIALS
Method of approach to a definition of the time for resin excess removal from the prepreg semiproduct package for PCM products on the base of Darsi’s law with preliminary experimental definition of permeability coefficient of reinforcing filler , viscosity and durability of the resin at autoclave forming temperature also, was described. For definition of permeability coefficient the calculated method of combination of leakage
(percolation) and reduction to elementary cell was used. The results of calculation of porous system permeability consisting of 7 /im 0 fibers depending on porosity of the structure were shown, but also experimental data on changing of viscosity and lifetime of the epoxy resin depending on temperature.
Key words: autoclave forming, content of resin, permeability, pressure, duration.
Соотношение армирующего наполнителя и связующего при прочих равных условиях определяет физико-механические характеристики полимерных композиционных материалов (ПКМ) и изделий на их основе. Поэтому при изготовлении изделий из ПКМ необходимо обеспечивать заданное соотношение компонентов. При автоклавном формовании используются препреги, в которых содержание связующего колеблется в пределах ±2,0^±2,5% от среднего номинала. В этом случае возникает необходимость удалить излишки связующего из собранного пакета препрегов. Удаление может происходить по двум направлениям: вдоль слоев препрегов к краям формуемых изделий и в поперечном к ним направлении, т. е. в направлении толщины изделия.
Удалить излишки связующего в направлении укладки слоев в случае, когда габаритные размеры изделия значительно превосходят его толщину, нереально, так как удаление произойдет только вблизи кромок изделия, что может привести затем к его короблению. Поэтому, если необходимо обеспечить равномерное распределение армирующего наполнителя и связующего по всему объему изделия, надо осуществлять удаление излишков в поперечном направлении. При этом, для исключения вытекания связующего в направлении укладки слоев, по периметру изделия укладываются ограничители типа герметизирующего жгута. А при использовании однонаправленных препре-гов на основе жгутовых наполнителей эти ограничители обязательны для предотвращения расползания слоев вместе с потоком отжимаемых излишков связующего.
Для сбора излишков связующего используются так называемые впитывающие слои (блидеры) (Проспект материалов для автоклавного формования фирмы «Airtech», 1993), обладающие большой впитывающей способностью. Фирма «Hexcel» (Проспект фирмы «Hex-cel», «Prepreg Technology», 1997) предлагает формулу расчета количества слоев блидера, необходимого для получения заданного содержания волокна Vf в пластике:
N =
r Wf • pr • V ^
f
(1)
где Ыъ - количество слоев блидера; А - поглощающая способность блидера, г/м2; Жг -поверхностная плотность смолы в препреге, г/м2; Ж/ - поверхностная плотность волокна в препреге, г/м2; V/ - заданное (выбранное) объемное содержание волокна, %; V -объемное содержание смолы, %; рг - плотность смолы, г/см3; р/ - плотность волокна, г/см3; Ыр - количество слоев в пакете ламината.
Там же (в проспекте) приводится предлагаемый порядок проведения опытов по оценке впитывающей способности материала блидера, поверхностной плотности смолы и наполнителя в препреге.
Как показывает практика, впитывающий слой играет важную, но не главную роль в удалении излишков связующего, так как процесс удаления требует значительных усилий и времени и зависит от многих факторов: структуры армирующего наполнителя, количества слоев препрега, количества излишков связующего, его жизнеспособности, величины давления формования и т. д.
Для установления необходимых времени и давления формования для удаления излишков связующего, необходимо провести предварительные экспериментальные исследования на небольших образцах, имитирующих реальные изделия по количеству слоев препрега из того же материала, что и изделие, при разных давлениях и скоростях нагрева пакета. Таким путем можно подобрать режимы формования, обеспечивающие получение изделий с заданным содержанием связующего. Но для изделий с другим соотношением связующего и наполнителя или же с другим количеством слоев препрега придется снова прибегать к трудоемким экспериментальным исследованиям.
Однако если знать так называемые фильтрационные характеристики наполнителя и связующего, то можно предварительно рассчитать время, которое потребуется для удаления излишков связующего из пакета препрега, состоящего из определенного количества слоев. Дело в том, что процесс удаления излишков связующего подобен процессу пропитки под давлением, используемому в авиации для изготовления антенных обтекателей. Процесс пропитки стеклонаполнителей подчиняется закону фильтрации, называемому законом Дарси - в честь французского инженера, заложившего основы теории фильтрации, - и выражается функцией:
V =— ср / йп, П
где V - скорость фильтрации; К - коэффициент проницаемости, характеризующий пористую среду, Д (дарси); ^ - коэффициент динамической вязкости жидкости; йр/йп -градиент давления в направлении «п».
Этот закон, как показано в работе [1], хорошо подтверждается экспериментальными данными, полученными при изготовлении изделий различных геометрических форм на основе стеклонаполнителей.
Время пропитки (в нашем случае время отжима) определяется зависимостью:
где т - продолжительность пропитки, с; т - пористость наполнителя; Н - протяженность пути в направлении пропитки, см; АР - перепад давления, действующего на связующее, МПа; К - коэффициент проницаемости, Д (дарси).
Как показывает формула (3), скорость пропитки зависит от фильтрационных характеристик наполнителя, к которым относятся коэффициент проницаемости К и величина пористости т, а также - от вязкости связующего ^. Кроме того, продолжительность пропитки (в нашем случае - отжима связующего) пропорциональна квадрату пути течения и обратно пропорциональна перепаду давления.
Проницаемость представляет собой проводимость пористой среды под влиянием перепада давления и является мерилом подачи жидкости через эту среду под действием проталкивающего давления. Коэффициент проницаемости определяется экспериментальным путем из соотношения:
т =---------,
2 К-АР
(3)
ь
(4)
где Q - объемный расход жидкости, см3/с; £ - площадь сечения пропитываемого образца, см ; АР - перепад давления, действующего на жидкость, МПа; Ь - длина пути в направлении пропитки образца, см.
Проницаемость определяется в единицах «дарси» (Д). За единицу «дарси» принимается проницаемость такого материала, в котором перепад давления в 1 ат (~98 кПа)
поддерживает расход жидкости в 1 см3/с при вязкости, равной 1 сантипуазу (сП), через куб с длиной ребра в 1 см. Таким образом в системе СИ
[к ]=^ •
с
м3 Н - с - м м2 - Н
.2
■ = м
м
м
Следовательно, 1 Д = 9,87-10"9 см2 = 9,87-10'13 м2 ~ 10'12 м2.
По Б.В. Дерягину [2], определение коэффициента проницаемости можно производить и по методике, основанной на измерениях объема жидкости, впитавшейся в образец за единицу времени:
К = (в- -«■) , (5)
28 • т -АР
V? V,2
где вх = —^; в2 = —^; У1 - объем жидкости, впитавшейся в пористое тело за время t1 ^ t2
под воздействием только капиллярного давления; У2 - объем жидкости, впитавшейся в пористое тело за время Ь при наличии добавочного давления АР; 8 - поперечное сечение образца; т - пористость (отношение объема пор ко всему объему).
Из уравнения (5) следует, что скорость пропитывания (впитывания блидером в нашем случае) тем больше, чем меньше давление воздуха внутри пористого тела, и что пропитывание (впитывание) с наложением внешнего давления должно протекать быстрее и тем лучше, чем больше пористость материала и выше давление. Наибольшая скорость отжима (впитывания в блидер) должна осуществляться при одновременном наложении вакуума и внешнего давления. Кроме того, ясно, что крупнопористое тело (впитывающий слой) пропитывается быстрее (быстрее и отдает жидкость), чем мелкопористое. Пористость, определяющая насыщенность наполнителя связующим, есть отношение объема пор ко всему объему тела.
Армирующие наполнители, особенно тканые, имеют два типа пористости: крупную (макро-) и мелкую (микро-) пористость. Макропористость - это поры между отдельными нитями или жгутами, а микропоры (капилляры) - между элементарными волокнами в нитях и жгутах. Как макро-, так и микропористость носят направленный характер вдоль нитей и жгутов, причем размер как макро-, так и микропор и их доля от общей пористости различны в зависимости от направления (вдоль основы, утка, по высоте).
Следовательно, армирующие наполнители представляют собой сложную пористую структуру с выраженной направленностью пор (капилляров), размеры которых колеблются от нескольких микрометров (внутри нитей и жгутов) до долей миллиметра (между нитями и жгутами). Следовательно, величина проницаемости К определяется как структурной пористостью материала, так и ее направленностью.
Не существует и не следует ожидать общей закономерности между пористостью и проницаемостью. Отсюда вполне очевидно, что любая деформация армирующего наполнителя, которая создает уменьшение пористости, влечет за собой в результате и уменьшение проницаемости. Уменьшение пористости заключается в уменьшении размеров пор, а отсюда в значительно большей степени уменьшается проницаемость, что было показано автором работы [3] для разных структур стеклотканей. В работе [3] определение проницаемости проводилось в направлении плоскости тканей вдоль основы и утка. При формовании изделий из ПКМ в автоклаве отжим связующего производится в основном в направлении, перпендикулярном плоскости композита, т. е. течение происходит поперек плоскости укладки препрега - в направлении убывания градиента давления. В этом направлении упаковка волокон имеет наибольшее уплотнение, особенно это относится к структурам из однонаправленных жгутов и лент на их основе или близких к ним (углеродных лент типа «Элур», УОЛ, кордных стеклоструктур и др.). Движение связующего должно преодолеть сопротивление плотно упакованных моноволокон, диаметр которых составляет всего несколько микрометров (7-10 мкм). Конечно, связующее будет перемещаться и в направлении плоскости укладки, отчего может произойти смещение слоев друг относительно друга, поэтому при формовании изделий из однонаправленных структур по краям обязательно выкладываются ограничительные жгуты. Однако следует заметить, что движение связующего в направлении торца изделия, при облегании его герметизирующей пленкой, в процессе вакуумирова-ния и создания давления будет остановлено, так как давление на композит со всех сторон одинаково. Таким образом, удаление излишков связующего может быть осуществлено в этом случае в основном только в направлении впитывающего слоя (блидера), укладываемого на пакет наполнителя (препрега). Для исключения прилипания блидера к изделию между ними прокладывается фторопластовая пористая ткань или перфорированная антиадгезионная пленка.
В работе [4] показано, что процесс вытекания смолы в перпендикулярном к плоскости оснастки направлении происходит по следующему механизму. После созда-
ния давления первый (верхний) слой препрега двигается по направлению ко второму слою. При этом связующее вытесняется из пространства между этими слоями и просачивается между волокнами первого слоя. Когда волокна первого слоя подходят к волокнам второго слоя, то уже два слоя двигаются в сторону третьего слоя, вытесняя связующее из пространства между вторым и третьим слоями. Эта последовательность событий повторяется для последующих слоев. Таким образом, взаимодействие волокон распространяется вниз по слоям препрегов в волнообразной форме. Причем при создании давления на препрег слои, из которых удалены излишки связующего, уплотняются и объем пористости и размеры пор уменьшаются, так же, но в еще большей степени уменьшается проницаемость, о чем говорилось выше, и течение вязкой композиции (связующего) проходит через поперечно расположенные капиллярные микропоры. При этом с увеличением количества уплотненных слоев время истечения, согласно формуле (3), растет пропорционально квадрату расстояния. Поэтому, при прочих равных условиях, сократить время отжима связующего можно и, очевидно, нужно как путем использования препрега с незначительным превышением заданного содержания его, или же при значительном увеличении перепада давления, или же за счет того и другого одновременно. Если взять, например, ленту углеродную ЭЛУР-П марки «А», имеющую по техническим условиям (ГОСТ 28006-88) поверхностную плотность 30±5 г/м при ширине 245±30 мм, плотность нити 1,71±0,04 г/см3, толщину монослоя по пластику 0,11^0,13 мм (в исходном состоянии при замере микрометром лента такой марки имеет толщину 0,14^0,15 мм), можно подсчитать, какова исходная и конечная (по пластику) пористость наполнителя. При толщине 0,15 мм и массе 30 г/м, средней ширине 245 мм, пористость составляет 52,3%, а при толщине 0,11 она будет равна уже 34,9%. Таким образом, при уплотнении ленты происходит уменьшение пористости на 17,4%, что, в свою очередь, влечет за собой изменение проницаемости. Если предположить, что этот объем занят связующим, то в процессе формования его необходимо удалить во впитывающий слой, чтобы получить толщину монослоя пластика 0,11 мм. Если плотность отвержденного связующего составляет 1,237 г/см3, как у связующего ЭДТ-69Н(М), то из каждого квадратного сантиметра поверхности одного слоя препрега необходимо удалить
0,004x1,237=0,004948-0,00495 г.
Для того чтобы можно было предварительно рассчитать время удаления излишков связующего из заданного пакета препрегов, необходимо, исходя из формулы (3), знать коэффициент проницаемости К, который определяется, как правило, экспериментальным путем [5].
Однако в работе [6] предлагается определение коэффициента проницаемости К производить путем расчета структуры со взаимопроникающими компонентами, используя метод сочетания теории протекания (теории перколяции и приведения к элементарной ячейке). Как следует из теории протекания [6], при малых значениях пористости т1, в материале появляются одиночные поры или в виде небольших изолированных скоплений или «кластеров». По мере возрастания т1 появляются и большие скопления пор, которые вместе с малыми скоплениями образуют так называемые изолированные кластеры (ИК). Когда пористость становится равной критической тс, ИК сливаются и образуют так называемый бесконечный кластер (БК), - система становится проводящей. При дальнейшем возрастании т1>тс БК увеличивается и поры пронизывают всю систему, образуя структуру с взаимопроникающими компонентами. Значение т1=тс называется порогом протекания. Рассматривая параллелепипед длиной Ь, шириной I и высотой к (рис. 1), к торцам которого приложена разность давления, а боковые стенки непроницаемы для потока жидкости, описывают сопротивление фильтрации Яф, которое может быть представлено по аналогии с электрическим сопротивлением в виде
Уравнение (6) использовалось авторами при расчете фильтрационных сопротивлений различных пористых систем. Рис. 1 представляет собой модель структуры со взаимопроникающими ком
(6)
Рис. 1. Структура со взаимопроникающими компонентами (К1 и К2 - коэффициент проницаемости поры и материала соответственно; q - поток жидкости)
Для расчета коэффициента проницаемости К пористых материалов, в которых каркас и поры образуют структуру со взаимопроникающими компонентами (к которым относятся тканые структуры), рекомендуется использовать формулу
К/К = с2 + у(1 - с)2 + 2у-с(1 - с)(у-с +1 - с)-1; V = К2/К1 , (7)
где с - положительный корень кубического уравнения
2с3-3с2 + 1= шх , (8)
с = 0,5+ВеоБф/3; 270°<ф<360°; К и К2 - проницаемость поры и каркаса соответственно; ш1 -пористость. При 0<ші<0,5 В=-1 и ф=агсов(1-2ш1); при 0,5<ші<1 В=1 и ф=агеов(2ш1-1).
На рис. 2 приведено графическое выражение зависимости (8): с от ш1, которой удобно пользоваться при расчетах. По данным [6], сопоставление расчетов по формуле (7) с результатами многочисленных экспериментов для систем с v^0 приводит к хорошему результату (расхождение не более 6,5%, т. е. сопоставимо с погрешностью численных расчетов).
0 0,2 0,4 0,6 0,8 1 ш1
Рис. 2. Зависимость с=/(тх): с - см. уравнение (8); шх - пористость материала со взаимопроникающими компонентами
При условии, что проницаемость каркаса (волокон) К2=0, а проницаемость поры К\, формула (7) примет вид
К=К г с2. (9)
Неизвестной величиной в этой формуле остается проницаемость поры К1, которую определим из следующих соображений. Расход жидкости через капилляр дается законом Гагена-Пуазейля
^-52 -АР (10)
я = т , (10)
32^- Т1
где 5, 5Т=/2 - диаметр и площадь поперечного сечения капилляра (БК); Ь1 - длина капилляра (БК) или длина траектории жидкой части в элементарной ячейке (см. рис. 1,
где / берется равной Ъ-Т).
С другой стороны, расход жидкости через элементарную ячейку может быть найден из закона Дарси:
К ^ АР (11)
Я = К------(11)
П т
В этом случае Ъ, 8=/ 2 - толщина и площадь поперечного сечения элементарной ячейки (см. рис. 1, где Ъ=1). Приравнивая (10) и (11), получим выражение для К элементарной ячейки:
К=(52/32) - (Т/Т1)- (^ / 8 )=(52/32) - (/2 /Т2) (Т/Т1 )=52/32^-с2, (12)
где ^=Т1/Т - извилистость поры (БК).
Из сравнения уравнений (9) и (12) видно, что проницаемость поры
К1=52/32^. (13)
Проницаемость определяется геометрией пористой структуры материала. На коэффициент проницаемости пористых материалов существенное влияние оказывает извилистость пор £, равная отношению средней длины траектории жидкой частицы в образце к длине образца. В этом случае формула (13) запишется в виде:
52 2
К = — с2, (14)
32-\
где £ - отношение длины траектории жидкости в образце к длине образца.
Для расчета проницаемости К предлагаются зависимости:
при шс<ш1<0,5 (шс - пористость границы проницаемости) с=0,5
при 0,5<ш1<1,0 2с3 -3с2 +ш1=0.
Г Л0,8
ш1 -шс
V0,5-шс у
Извилистость пор оценивается по формуле
/ \ 0,4
( 1—т 1
£, = -----— при да1>0,2. (15)
I т1 —тс )
Графическая зависимость £ от пористости представлена на рис. 3.
Наибольшую неопределенность при расчетах по формуле (15) вносит выбор диаметра пор 5. Для волокнистых структур этот размер предлагается оценивать по формуле:
5=0,8941/(с-1)]'1, (16)
где ё - диаметр волокна, а с определяется из рис. 2.
Графическая зависимость 5 от ш1 представлена на рис. 4.
На основании предлагаемых методик, авторами проведен расчет изменения коэффициента проницаемости в зависимости от пористости ш1 жгутового наполнителя в пределах 30-55%. В расчете диаметр волокна принят 7 мкм. Результаты расчета представлены в виде графика (рис. 5).
Полученная зависимость относится только к структурам типа однонаправленных жгутовых наполнителей, в которых при качественной укладке (без пропусков) жгутов отсутствуют межжгутовые зазоры, а пористость определяется межволоконным пространством. Пористость же, а следовательно и проницаемость, тканых структур определяется межнитяными зазорами, размеры которых на 2 и более порядка больше меж-волоконных зазоров.
Пользуясь полученными данными, можно проводить предварительный расчет времени, необходимого для удаления излишков связующего из пакета препрега при автоклавном формовании с использованием впитывающего материала (блидера).
Второй составляющей в формуле (3) расчета времени отжима связующего является вязкость связующего. В любой момент времени исходя из требования динамического равновесия скорость вытекания связующего из композиции равна скорости поступления ее во впитывающий материал. Температура, а следовательно, и вязкость связующего внутри впитывающего слоя предполагаются не зависящими от положения, но не от времени. Время как раз имеет большое значение, так как связующее при повышенных температурах обладает ограниченной жизнеспособностью, причем чем температура выше, тем меньше жизнеспособность. Поэтому для достижения максимальной скорости удаления связующего из композита необходимо найти оптимальные его свойства, к которым относятся наименьшая вязкость ^ при максимальной жизнеспособности
[7].
0,3 0,35 0,4 0,45 0,5 0,55 т1 0,3 0,35 0,4 0,45 0,5 0,55 т1 0,3 0,35 0,4 0,45 0,5 0,55т1
Рис. 3. Связь коэффициента Рис. 4. Зависимость 5 от Рис. 5. Зависимость прони-
ИЗВИЛИСТОСТИ £, С ПОРИСТОСТЬЮ Ш\
ПОРИСТОСТИ Ш\
(5 - дИамеТр поры)
цаемоСТИ от пористости матерИала (волокнИСтая Структура С дИаметром волокна 7 мкм)
Формула (3) показывает, что время удаления излишков связующего пропорционально вязкости связующего. Поэтому желательно, чтобы в процессе формования (отжима) связующее, во-первых, имело минимальную вязкость, во-вторых, она (вязкость) сохранялась бы постоянной в период отжима излишков.
Жизнеспособность связующего, мин
Рис. 6. Динамическая вязкость эпоксидного связующего при температуре 60 (•), 80 (о), 90 (■), 100 (□), 110 (▲) и 120°С (Л)
Определение вязкости и жизнеспособности эпоксидного связующего в зависимости от времени проводилось на реоспектрометре фирмы «Инстрон» модели 3250. Определялись динамические характеристики связующего*. Для этих целей использовались образцы связующего, из которого был удален растворитель путем вакуумирования при температуре 40-50°С. Скорость нагрева при определении вязкости и жизнеспособности связующего составляла 2-3°С в минуту. Оценка проводилась в интервале температур 80-100°С. Результаты оценки представлены на рис. 6. На рис. 6 видно, что минимальная вязкость связующего при жизнеспособности его —110 мин приходится на температуру 80°С, которая и рекомендуется для расчета времени отжима связующего из пакета препрегов. Пользуясь данными, полученными при определении оптимальных показателей вязкости связующего при повышенных температурах, и коэффициентами проницаемости армирующих наполнителей, можно произвести расчет времени отжима связующего в зависимости от толщины пакета наполнителей при постоянном или изменяемом давлении.
* Исследования проводились С.И. Литовченко.
ЛИТЕРАТУРА
1. Душин М.И., Хрульков А.В., Мухаметов Р.Р., Чурсова Л.В. Особенности изготовления изделий из ПКМ методом пропитки под давлением //Авиационные материалы и технологии. 2012. №1. C. 18-26.
2. Дерягин Б.В., Захаваева Н.Н., Лопатина А.М. Исследование фильтрации растворов электролитов в высокодисперсных порошках //В кн.: Исследование в области поверхностных сил. М.: АН СССР. 1961.
3. Душин М.И. Автореф. дис. на соиск. степени канд. техн. наук, ВИАМ, 1971, 32 с.
4. Springer G.S. Resin Flow During the Cure of Fiber Reinforced Composites //Journal of Composite Materials. 1982. V. 16. Р. 400-410.
5. Душин М.И., Хрульков А.В., Мухаметов Р.Р. Выбор технологических параметров автоклавного формования деталей из полимерных композиционных материалов //Авиационные материалы и технологии. 2011. №3. С. 20-26.
6. Дульнев Г.Н., Новиков В.В. Процессы переноса в неоднородных средах. М.: Энер-гоатомиздат. 1991. 224 с.
7. Мухаметов Р.Р., Ахмадива К.Р., Ким М. А., Бабин А.Н. Расплавные связующие для перспективных методов изготовления ПКМ нового поколения /В сб. Авиационные материалы и технологии: Юбилейный науч.-технич. сб. (приложение к журналу «Авиационные материалы и технологии»). М.: ВИАМ. 2012. С.260-265.
