Numerical simulation of the process manufacture of large-scale composite shell taking into account thermo viscoelastic.pdf Введение Полимерные композиционные материалы занимают одно из лидирующих мест среди конструкционных материалов. Улучшенные физико-механические и эксплуатационные характеристики позволили значительно расширить применение полимерных композитов в машиностроении [1], строительстве [2], авиастроении [3], недропользовании [4] и других отраслях промышленности. Полимерные 166 Сахабутдинова Л.Р., Сметанников О.Ю., Ильиных Г.В. Численное моделирование процесса материалы применяются в качестве защитных и защитно-упрочняющих покрытий и прослоек [5], для создания элементов конструкций [6], а также для создания изделий с использованием метода намотки [7-8]. Намотка - это процесс изготовления высокопрочных армированных изделий, форма которых определяется вращением произвольных образующих (оправок). Для связующей матрицы выбирают эпоксидные или полиэфирные смолы, полимеры сложных виниловых эфиров. Выбор связующего зависит от типа армирующего материала [9-10]. Оправки для намотки должны иметь форму и размеры, точно соответствующие внутренней поверхности наматываемых изделий, эффективно воспринимать силовые и температурные воздействия при намотке и отверждении связующего. При производстве крупногабаритных изделий преимущественно используются разрушаемые оправки, для их изготовления возможно применение песчанополимерных смесей или легкоплавких материалов [11]. В процессе роста объемов производства изделий из композитных материалов, в том числе и способом мокрой намотки, актуальной становится проблема точного прогнозирования несущей способности конструкций, что тесно связано с разработкой эффективных математических моделей поведения как композиционных материалов в целом, так и матрицы, связующего и армирующих элементов в частности. Качественные модели поведения композитов позволят оценить влияние наиболее значимых конструктивных и технологических факторов на напряженно-деформированное состояние и прочность конструкций [12]. На протяжении последних десятилетий многие группы и отечественных, и зарубежных ученых занимались исследованием проблем, связанных с изготовлением намоточных изделий. Основные исследуемые вопросы можно разделить на три группы: определение эффективных механических характеристик [13]; определение усилий, реализуемых в материале в процессе намотки армирующих волокон на формообразующую оправку [14]; определение остаточных напряжений после изготовления [15]. Отдельно стоит проблема изучения термовязкоупругого поведения материалов, используемых в процессе производства конструкций [16]. В силу особенностей процесса изготовления крупногабаритных намоточных изделий, таких как длительность процесса намотки (от 10 до 15 суток) и высокотемпературный режим полимеризации связующего (до 170°С), целесообразно учитывать реологические процессы, протекающие в материале песчано-полимерных оправок и наматываемых оболочек, поскольку они оказывают существенное влияние на напряженно-деформированное состояние изделия. Экспериментальные методы определения характеристик достаточно дорогостоящие и имеют ряд физических ограничений, не позволяющих установить весь набор необходимых для проектирования характеристик [17]. Поэтому развитие численных моделей является перспективным и быстроразвивающимся направлением в вопросах изучения и применения полимерных композиционных материалов. Объединение в единый расчетный модуль задач по определению эффективных характеристик материалов, усилий, реализуемых в процессе укладки армирующих материалов, температурных полей и деформаций в процессе отверждения позволит рассмотреть поведение конструкции с технологической оснасткой на всех ключевых этапах процесса производства, а также локализовать зоны возникновения дефектов. 167 Механика / Mechanics Объект исследования Объектом данного исследования является крупногабаритное намоточное изделие, представляющее собой пластиковую оболочку типа «кокон» с технологической оснасткой, состоящей из вала, песчано-полимерной оправки и прижимающего устройства, обеспечивающего фиксацию оправки. На рис. 1 представлены схема рассматриваемой конструкции и модель для решения тестовых задач в виде трехслойного сегмента цилиндра. Поскольку поставлена задача разработки комплексной методики прогнозирования напряженно-деформированного состояния конструкции в процессе изготовления, то в качестве предмета исследования выбраны два этапа технологического процесса - намотка и отверждение связующего. ххх; - формообразующая оправка модель прижимном оснастки пластиковая оболочка раздели тельный слои соорочныи ват Рис. 1. Расчетная схема конструкции: a - схема реальной конструкции; b - тестовая модель Fig. 1. Scheme of the construction: (a) real construction and (b) simplified model Реальная конструкция не осесимметрична, так как секции формообразующей оправки укреплены ребрами жесткости и имеют технологические отверстия. Поэтому при численном моделировании разработан трехмерный конечно -элементный аналог, учитывающий не только габариты конструкции, но и ее технические особенности. Конечно-элементная модель параметризирована по размерам, углам, схеме армирования оболочки, усилиям и режиму полимеризации, что значительно расширяет ее применимость на производстве. Подобные крупногабаритные изделия, как правило, малосерийные. К готовой оболочке предъявляются повышенные требования по прочности и точности геометрической формы. Особенно важно, что готовое изделие не должно иметь геометрических дефектов на внутренней поверхности. На появление таких отклонений, прежде всего, влияют свойства и поведение оправки в процессе производства. Поэтому при проектировании изделий данного класса важно вычислить зависимость нормального давления на поверхности оправки от условий техноло-168 Сахабутдинова Л.Р., Сметанников О.Ю., Ильиных Г.В. Численное моделирование процесса гического процесса. Для этого применяют упрощенные модели в упругой постановке, поскольку экспериментально измерить значения контактного давления на производстве трудоемко, нетехнологично и финансовозатратно [18]. Использование численных моделей позволяет как оценить эволюцию нормального давления в конкретных точках на поверхности оправки, так и получить интегральные оценки. В рамках исследования для первого этапа валидации модели выбрано три ключевых точки на поверхности оправки, обозначенные на рис. 1, для количественной и качественной оценки эволюции нормального давления. Точки 1 и 3 расположены на полусферических донных участках, точка 2 - на цилиндрическом. Таким образом, разработка трехмерного численного аналога - новый этап для повышения точности расчетов и численных исследований напряженно -деформированного состояния систем «оправка-оболочка». Основываясь на особенностях технологического процесса изготовления пластиковых оболочек, построены расчетные схемы этапов технологического процесса. Определена система внешних силовых и температурных воздействий. Постановка задачи Поскольку технологический процесс сопряжен с температурным воздействием, необходимо последовательно решать два типа задач. Для определения изменения температурных полей в конструкции на протяжении процесса полимеризации решается задача нестационарной теплопроводности в условиях конвективного теплообмена с окружающей средой. Общая постановка представлена в [19]. Для определения эволюции напряженно-деформированного состояния конструкции в процессе изготовления решается квазистатическая краевая задача. В общем случае постановка задачи без учета объемных сил описана в [20]. В качестве исходных данных для описания поведения материалов выступали экспериментальные кривые релаксации материалов оправки [21], однонаправленного пластика и эпоксидного связующего, полученные при различных температурах. Согласно [22] рассматриваемые материалы можно принять линейновязкоупругими и термореологически простыми. В качестве определяющих соотношений выбрана вязкоупругая модель Prony, использующая в качестве релаксационного ядра сумму экспонент с постоянным модулем объемного сжатия. В выбранной модели предполагается, что материал испытывает только сдвиговую релаксацию при постоянном значении модуля объемного сжатия. Для получения экспериментальных параметров удобно использовать данные испытаний на одноосное растяжение-сжатие. Модуль релаксации для одноосного растяжения-сжатия имеет вид: E (t ) = E0 Ne f +Z ci exp '=1 v Pi (1) здесь c - относительные модули растяжения-сжатия для времен релаксации , N - количество времен релаксации растяжения-сжатия. Определяющее соотношение для одноосного случая: t j(t w о Г Ne f t-tYI Q! E» + Eo Z ci exP '=1 V Pi )J afe(x), (2) 169 Механика / Mechanics где Р Р, A (T) приведенное время, Л(т) - функция сдвига. Поскольку матери ал предполагается термореалогически простым, то используется функция сдвига Вильямса-Ландела-Ферри Q (Г - T) C +(Г - Г) (3) lg(A(T)) где Г - текущая температура, Гг - постоянная базовая температура, C , С2 -эмпирические постоянные для материала. Механическое поведение материалов в среде ANSYS Mechanical можно описать средствами модели Prony, предполагающей для среды наличие одного вязкоупругого оператора. При этом данная модель допускает описание поведения материала в виде сочетания ортотропно-упругих характеристик и релаксации компонент пропорционально упругим значениям [23]. Параметрами модели Prony являются ct, Р, - коэффициенты аппроксимации обобщенной функции релаксации (1). Учет влияния температуры на скорость протекания релаксационных процессов описывается моделью Shift (в документации описана как Williams-Landel-Ferry Shift Function) с параметрами модели Вильямса-Ланделла-Ферри (3). Особенности численного решения задачи Для разработки численного аналога исследуемой конструкции предложено несколько новых подходов к конечно-элементному моделированию ключевых технологических этапов изготовления оболочки вращения, получаемой методом непрерывной намотки. Одной из рекомендаций является разделение процесса изготовления оболочки на два основных этапа: намотку высокомодульных лент на формообразующую оправку и полимеризацию связующей матрицы оболочки в печи. Такой подход позволяет проводить последовательные расчеты и оперативно варьировать параметры технологических процессов - усилия натяжения армирующих лент, усилия в прижимной оснастке, а также режим полимеризации. При численном моделировании процесса изготовления композиционной оболочки необходимо определение эффективных вязкоупругих характеристик материала. Исходя из особенностей армирования, материал оболочки на макроуровне конструкции принято считать ортотропным. Для решения проблемы прогнозирования характеристик оболочки при изменении схемы армирования была разработана и предложена комплексная экспериментально-расчетная методика, основанная на решении краевых квазистатических задач термомеханики и принципах многоуровневого моделирования [24]. Моделирование технологического процесса изготовления осуществляется в четыре последовательных этапа. На первом этапе определяются усилия, возникающие в конструкции при натяжении высокомодульных армирующих лент на формообразующую оправку. Для реализации первого этапа моделирования на тестовой схеме (см. рис. 1, b) решена квазистатическая задача механики деформируемого тела в упругой постановке, имитирующая послойное наращивание толщины оболочки с применением технологии «оживления» конечных элемен-170 Сахабутдинова Л.Р., Сметанников О.Ю., Ильиных Г.В. Численное моделирование процесса тов. В каждом элементе слоя оболочки в качестве начальных условий вносились значения компонент тензора напряжений, возникающих от усилий наматывания армирующих лент. Для оценки возможности рационализации процедур численного моделирования проведено исследование влияния послойного и одномоментного оживления наматываемых слоев на изменение значения контактного давления по толщине оболочки и на поверхности оправки во время процесса намотки (рис. 2). Рис. 2. Динамика контактного давления: a - по толщине; b - в процессе намотки; 1 - послойное оживление; 2 - одномоментное оживление Fig. 2. Contact pressure curve: (a) by the thickness and (b) on the surface of the mandrel during winding; (1) layer-by-layer activation; (2) single-step activation При одномоментном оживлении всех слоев значение максимального давления получилось на 10% ниже эталонного. Для коррекции данного значения было предложено ввести поправочный коэффициент для значений компонент тензора напряжений. На основе полученных данных сделан вывод о целесообразности объединения отдельных слоев в один слой с эффективными характеристиками и эффективными начальными напряжениями, которые обеспечивают необходимый уровень давления от оболочки на оправу. Для исследования реальной конструкции реализован параметризированный модуль для вычисления начальных напряжений в многослойной оболочке от действия силы натяжения армирующих волокон. На втором этапе решается нестационарная задача теплопроводности для определения эволюции температурных полей в конструкции «оправка-оболочка» на протяжении полимеризации. При решении задачи теплопроводности реализуются граничные условия третьего рода, определяющие нагрев оболочки путем конвективного теплообмена с окружающей средой. Значения температуры внешней среды задаются на внешней поверхности оболочки и технологической оснастки, начальное распределение температур в конструкции постоянно. Коэффициенты конвективного теплообмена определены из серии вычислительных экспериментов, в рамках которых проводилось сравнение получаемого решения с термометрией реальной конструкции в процессе полимеризации. Третий этап включает моделирование начального напряженного состояния в оболочке, возникающего от усилия натяжения ленты в процессе намотки, и расчет релаксации напряжений в материалах оправки и оболочки на протяжении времени, соответствующего процессу намотки. Вычисленные начальные деформации учитываются в материале оболочки как преднапряженное состояние. 171 Механика / Mechanics Затем проводится решение квазистатической задачи. На четвертом этапе проводится расчет по определению эволюции напряженно-деформированного состояния конструкции в процессе полимеризации с использованием итерационной процедуры, в которой для последовательных моментов времени подгружается соответствующее решение задачи нестационарной теплопроводности в качестве объемной нагрузки. Полученные таким образом результаты, основанные на экспериментальных данных и предложенных приемах конечно-элементного моделирования, позволяют оценить вклад релаксации напряжений в оправке и оболочке на эволюцию напряженно-деформированного состояния в процессе изготовления реальной конструкции. О многоуровневом подходе к прогнозированию вязкоупругих характеристик материала оболочки В рамках численного моделирования процесса изготовления исследуемой конструкции при изменении схемы армирования материала многослойной оболочки каждый раз необходимо решать проблему определения эффективных вязкоупругих характеристик. Выбранный метод многоуровневого моделирования позволяет выделить этот этап решения в отдельную расчетную процедуру в рамках решения общей задачи. В исследовании выбрано два масштабных уровня. На первом уровне определяются свойства ячейки периодичности однонаправле-ного армированного слоя. На втором уровне определяются механические характеристики пластиковой оболочки, рассмотренной на представительном объеме наматываемого композита в виде набора отдельных слоев с постоянными углами армирования. В композиционных оболочках слои в зависимости от угла укладки арматуры разделяют на два вида - кольцевые и спиральные. Кольцевые слои - слои с арматурой, уложенной перпендикулярно оси изделия. Спиральные слои представляют собой переплетенные ленты арматуры с симметричным относительно оси оболочки углом. Традиционно [13, 25] при решении задач по определению эффективных характеристик спиральные слои принято разделять на два слоя с симметричными углами армирования. Направление армирования в каждом слое соответствует заранее определенной схеме укладки волокон. Данный подход успешно применяется при прогнозировании эффективных упругих характеристик материалов. Численная процедура прогнозирования вязкоупругих характеристик материала реализована в ANSYS Mechanical APDL. На каждом масштабном уровне решался набор квазистатических задач итерационным методом [26]. Метод основан на решении задач линейной вязкоупругости с помощью квазикон-статных операторов. При численной реализации задачи методом конечных элементов выбранный подход вырождается в решении упругой задачи на каждой итерации (каждом шаге) интегрирования по времени. При рассмотрении первого масштабного уровня предполагаем, что волокна уложены параллельно друг другу и материал имеет однонаправленную структуру. Такой материал описывается в литературе [13, 25] как трансверсально-изотропный в плоскостях, перпендикулярных направлению армирования. В упругой постановке поведение такого композита характеризуется пятью независимыми кон-172 Сахабутдинова Л.Р., Сметанников О.Ю., Ильиных Г.В. Численное моделирование процесса стантами. Для определения необходимых констант формулируется три вычислительных эксперимента для трехмерной ячейки периодичности: исследование продольной, поперечной и сдвиговой продольной деформаций ячейки. Для определения коэффициентов линейного температурного расширения формулируется вычислительный эксперимент по вычислению температурной деформации. Для обобщения метода квазиконстантных операторов на данный случай упругие константы заменяются соответствующими вязкоупругими операторами. В таком случае ставятся вычислительные эксперименты на исследование релаксации ячейки периодичности при продольном, поперечном и сдвиговом продольном деформировании. Исходными данными для моделирования ячейки периодичности выступают результаты аппроксимации экспериментальных кривых продольной релаксации однонаправленного пластика и эпоксидного связующего, полученные при различных температурах. Для материала связующего функция релаксации была построена ранее на основе экспериментальных данных [27]. Однако для численной реализации первого уровня необходимо выделить функцию релаксации волокна. При моделировании механических свойств материал волокна считается изотропным и термореологически простым. Для получения функции релаксации материала используются соотношения для определения эффективного продольного модуля однонаправлено-армированного композиционного материала в операторной форме [13, 25]. Выделенная в результате кривая функции релаксации волокна представлена на рис. 3, a. Результаты реализации предложенного алгоритма решения задачи по прогнозированию эффективных характеристик в виде кривой ядра продольной релаксации для однонаправленного органопластика представлены на рис. 3, b. Значение мгновенного модуля E* в поперечном направлении однонаправленного пластика получено равным 33 ГПа, а длительного - 2 ГПа. Рис. 3. Кривые ядра релаксации при различных температурах (справа налево) 20°С, 40°С, 60°С, 80°С, 100°С, 120°С, 140°С, 160°С: а - волокно, б - уровень монослоя Fig. 3. Relaxation kernel curves at different temperatures (from right to left) 20°С, 40°С, 60°С, 80°С, 100°С, 120°С, 140°С, 160°С: (a) fiber and (b) laminar monolayer Определение эффективных вязкоупругих характеристик материала оболочки в виде представительного объема многослойного пакета проведено на основе вычислительных экспериментов, в рамках которых исследуется эволюция напряженно-деформированного состояния с использованием метода конечных элементов. Поскольку исследуемый материал на макроуровне в глобальной системе ко-173 Механика / Mechanics ординат конструкции принято считать ортотропным, то для определения всех искомых операторов требуется проведение следующих вычислительных экспериментов: релаксация ячейки при деформировании в направлении кольцевого армирования; релаксация ячейки при деформировании в направлении, перпендикулярном кольцевому армированию; релаксация ячейки при деформировании в направлении, перпендикулярном плоскости слоев; сдвиговая релаксация ячейки в трех плоскостях. Результаты реализации предложенного алгоритма решения задачи по прогнозированию эффективных вязкоупругих характеристик на уровне многослойного * пакета представлены на рис. 4. Для операторного продольного модуля мгновенное и длительное значения на 40% ниже, чем у однонаправленного пластика. * Мгновенное значение операторного поперечного модуля Еу получено на уровне 55 ГПа, а длительное значение равно 26 ГПа, что выше значений для монослоя на порядок. Таким образом, рассматриваемая схема армирования приводит к более равномерному распределению прочностных характеристик и повышению прочности при повышенных температурах и / или длительных периодах эксплуатации. Рис. 4. Операторные ядра релаксации оболочки при различных температурах (справа налево) 20°С, 40°С, 60°С, 80°С, 100°С, 120°С, 140°С, 160°С: a - модуль Е*, b - модуль Е* Fig. 4. Relaxation kernels at different temperatures (from right to left) 20°С, 40°С, 60°С, 80°С, 100°С, 120°С, 140°С, 160°С: Young's moduli (a) Е* and (b) E* Для адаптации полученных кривых релаксации в вычислительном комплексе конечно-элементного анализа выбрана методика обработки данных методами нелинейного программирования. Параметры коэффициента температурно-временного сдвига (3) определяются решением системы линейных алгебраических уравнений, составляемых из сравнения кривых релаксации при трех различных температурах. Коэффициенты аппроксимации (1) определяются минимизацией среднеквадратичной невязки. Поскольку в ANSYS Mechanical APDL модель вязкоупругого поведения материалов Prony предполагает наличие для среды одного ядра релаксации, исходя из особенностей деформирования исследуемой конструкции сделано предположение, что наиболее существенным на эволюцию напряженно-деформированного состояния будет влияние релаксационных процессов в направлении кольцевого армирования. Параметры модели Prony определялись аппроксимацией ядра операторного модуля Е* . Для вычисления коэффициентов аппроксимации выбрана температура приведения Tb = 80°С, и рассмотрен временной интервал, на котором 174 Сахабутдинова Л.Р., Сметанников О.Ю., Ильиных Г.В. Численное моделирование процесса происходит достижение длительного модуля. Параметры аппроксимации определялись из решения задачи о поиске минимума среднеквадратичной невязки, относительная погрешность аппроксимации не превысила 0,6%. Анализ влияния модели поведения оболочки на значение нормального давления на поверхности оправки Для оценки влияния модели Prony с полученными параметрами на эволюцию напряженно-деформированного состояния конструкции проведено сравнение релаксации напряжений, возникающих от усилия укладки лент, на поверхности оправки при температуре 80°C. Серия вычисленных экспериментов с тремя различными вариантами представления материала оболочки выполнена на тестовой модели (рис. 1, b). Методом переменных параметров упругости были реализованы расчеты для послойного представления оболочки и оболочки с эффективными вязкоупругими характеристиками слоев. Третьим вариантом было однородное представление материала оболочки с набором упругих характеристик, описывающих ортотропность материала, и одним ядром релаксации посредством модели Prony с учетом влияния температуры на скорость протекания релаксационных процессов. Результаты решения задач представлены на рис. 5 в виде средних напряжений на поверхности оправки в трех направлениях. Обозначения соответствуют цилиндрической системе координат. ““ - оболочка с эффективными характеристиками и аппроксимацией одного ядра релаксации моделью Prony Рис. 5. Эволюция средних напряжений при температуре 80°C Fig. 5. Evolution of mean stresses at temperature of 80°C - - оболочка с послойным представлением материала, с учетом всех ядер релаксации ■Ж- - оболочка с эффективными характеристиками, с учетом всех ядер релаксации 175 Механика / Mechanics Решения тестовой задачи для многослойного представления и для осредненного методом переменных параметров упругости полностью идентичны. При сравнении с решением квазистатический задачи с моделью Prony получено, что отклонения длительных значений средних напряжений модели с аппроксимацией составляют: (az) - 3,7%, (a^j - 23%, (az) - 34%. На интервале времени до 107 с максимальные погрешности для разных направлений равны 3,7, 8,7 и 6,5% соответственно. Сходимость решений по радиальному напряжению наиболее значима, поскольку оправка подвергается сжатию со стороны оболочки. Таким образом, выбранная модель описания поведения материала оболочки в исследуемой конструкции применима при решении задач по оценке прочности формообразующей оправки в процессе изготовления композиционной оболочки. Сравнение нормального давления на поверхности оправки для двух моделей поведения оболочки Используя вышеизложенные приемы конечно-элементного моделирования, реализовано решение комплексной задачи прогнозирования эволюции напряженно-деформированного состояния формообразующей оправки в процессе изготовления композиционной оболочки. Для оценки влияния релаксационных процессов в материале оболочки на значение нормального давления на поверхности оправки решено провести сравнение с упругой моделью, применяемой сегодня на производстве. Для этого моделирование этапов намотки и полимеризации оболочки реализовано с двумя моделями поведения материала оболочки - упругой и термовязкоупругой. Материал оправки описан как термовязкоупругий, параметры модели Prony и Shift определены по аналогичному алгоритму, что и для оболочки. В качестве исходных данных выступали результаты натурных экспериментов на сжатие при различных температурах [21]. Сравнение результатов расчетов проводится по значениям нормального давления в трех контрольных точках. Поскольку результаты, полученные для точек 1 и 3, идентичны, для качественного сравнения взято их среднее значение. Для наглядности и интерпретации полученных результатов оценка проведена для двух временных интервалов -этапа намотки оболочки и процесса полимеризации связующего. Результаты сравнения на этапе намотки представлены на рис. 6, а, на этапе полимеризации -на рис. 6, b. На этапе намотки влияние модели поведения материала оболочки несущественно, поскольку наибольший вклад в снижение контактного давления на поверхности оправки вносит вязкоупругость материала оправки. Это подтверждается и реологическим поведением материала при нормальной температуре. По результатам экспериментальных данных при температуре 22°С снижение напряжений в образцах к моменту времени 105 с составило порядка 70%, с 8,5 до 2,5 ГПа. На этапе полимеризации значения контактного давления на поверхности оправки, полученные для вязкоупругой модели оболочки, ниже на 10-15%, что также соответствует снижению продольного модуля материала многослойной оболочки. Например, для продольного модуля оболочки падение составляет 14,5% при всех рассмотренных температурах, с 70 до 60 ГПа. 176 Сахабутдинова Л.Р., Сметанников О.Ю., Ильиных Г.В. Численное моделирование процесса Р, МПа о Р, МПа 175 150 100 *"'~ѳ-ё 300 240 250 время, час режим полимеризации время, час контрольные точки на поверхности оправки Э - оболочка с учетом реологии - - линейно упругая оболочка Рис. 6. Сравнение эволюции нормального давления на поверхности оправки: a - в процессе намотки, b - в процессе полимеризации Fig. 6. Comparison of evolution of normal pressure on the surface of the mandrel: (a) winding process and (b) polymerization process На основе полученных данных можно сделать вывод о целесообразности и важности учета термовязкоупругого поведения в композиционном материале многослойной оболочки при прогнозировании эволюции напряженно-деформированного состояния конструкции в целом и оценке прочности формообразующей оправки в процессе изготовления изделий. Заключение В результате численного моделирования для реальной конструкции получены количественные и качественные закономерности эволюции нормального давления на поверхности формообразующей оправки на протяжении процесса намотки и полимеризации с учетом реологического поведения ключевых элементов конструкции - многослойной оболочки и оправки. Предложенные методы и разработанные программные модули позволят в дальнейшем количественно оценивать влияние на напряженно-деформированное состояние конструкции таких технологических параметров, как схема армирования, усилия натяжения армирующего материала, усилия, возникающие от действия технологической оснастки, режима термообработки. Планируемые в дальнейшем исследования по совершенствованию и верификации разрабатываемой модели позволят оптимизировать технологические процессы и повысить качество получаемых изделий. Предложенные методы и подходы к конечно-элементному моделированию дадут возможность существенно расширить области применения изделий из полимерных композиционных материалов, получаемых методом непрерывной намотки. Вклад авторов: все авторы сделали эквивалентный вклад в подготовку публикации. Авторы заявляют об отсутствии конфликта интересов.

Скачать электронную версию публикации