Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Программа и тезисы докладов VI Всероссийской научно-технической конфе..pdf
Скачиваний:
11
Добавлен:
15.11.2022
Размер:
13.42 Mб
Скачать

НЕЛИНЕЙНЫЕ МОДЕЛИ ПРОГНОЗИРОВАНИЯ ЭФФЕКТИВНЫХ УПРУГИХ СВОЙСТВ ВЫСОКОПОРИСТЫХ композитов

Макарова Е.Ю.

Пермский государственный технический университет

В работе [1] рассматриваются линейные модели прогнозирования макроскопических свойств квазиизотропных композитов. Экспериментальные и теоретические исследования свидетельствуют о том, что линейные модели имеют достаточную для практики точность лишь для пористости, не превышающей 8%, что имеет важное значение для спеченных порошковых структурно-неоднородных материалов. Для нового класса высокопористых композитов эта теория не дает удовлетворительной точности. Поэтому в данной работе поставлена задача разработать нелинейные структурно­ феноменологические модели механики композитов, которые имеют необходимую для инженерной практики точность прогнозирования эффективных свойств. Отметим при этом, что впервые нелинейные структурно-феноменологические модели применены к прогнозированию эффективных свойств гранулированных композитов с полимерными матрицами с жесткими сферическими включениями [1, 2]. Эти модели дали удовлетворительную сходимость теории и эксперимента. В качестве структурно­ феноменологической модели для описания высокопористых композитов предлагаются следующие соотношения:

%m (r) =C"m„ {\-p ll-X {r)), (1)

где дф„„(г) - структурные модули упругости композита, cfjmn- модули упругости матрицы, р - объемная пористость композита, Х(г) - случайная индикаторная функция:

f 1

с вероятностью

р

(2)

=

.

с вероятностью (\ - р ) ,

 

Из соотношения (1) с учетом (2) находим эффективные модули упругости композита в приближении Фойгта

с;тп= { % Л г ))= с 1 Л - р )г (з)

Из соотношения (3) видно, что в приближении Фойгта эффективные модули определяются нелинейной

зависимостью от пористости.

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

c;jm„= с " „(1 - р)2+ J^ C'^C'I A 1 - р)ъ

(4 )

Расчеты эффективного модуля сдвига по соотношению (4) показывают, что корреляционная поправка невелика и значение модуля сдвига попадает в вилку Фойгта-Рейсса по нелинейной модели, что свидетельствует о необходимости применения нелинейных моделей к расчету эффективных упругих

свойств высокопористых композитов.

Работа выполнена при поддержке Российского фонда фундаментальных исследований (грант РФФИУрал № 02-01-96-403)

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

1.Соколкин Ю.В., Ташкинов А.А. Механика деформирования и разрушения структурно неоднородных

2.1а в ! Ташттов А.А. Макроскопические упругие модули гранулиропаннык эластомерных

композитов. - Механика эластомеров, 1981, вып. 5, с. 46-53.

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

Макарова Е.Ю.

Пермский государственный технический университет

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

eijm„(r) = efjmn{l-pXl-Xir)),

( 1)

где где вутп(г) - структурные модули упругости композита, Qjjmn - случайные модули упругости матрицы,

р - объемная пористость композита, Л(г) - случайная индикаторная функция:

Г 1

с вероятностью

р

\( г ) = \

с вероятностью (\ -

(2)

р )

Предполагается, что поля упругих свойств 0 ^ и случайной индикаторной функции Л(г) являются

собственными, то есть отсутствует взаимная корреляция между полями. Осредняя соотношение (1), находим эффективные модули упругости для модели Фойгта:

с>1т» = ^ А г ) ) = с ;тп{’ 1 - р )г

(з )

Аналогичные соотношения получаем для модели Рейсса:

(4)

где П ijm„(r) - структурные коэффициенты податливости композита, П “.тп- коэффициенты податливости матрицы.

На основе решения краевой задачи микромеханики композитов, аналогично тому, как это сделано в работе [1], находим формулу для расчета макроскопических модулей упругости высокопористого композита с учетом естественного разброса упругих свойств матрицы:

^ijmn ~ C (jmn (l р ) + * А х у В 6 Р )

И II

/ о//

 

 

(5)

(/аР^убтл + (

Как видно из формулы (5), последнее слагаемое учитывает естественный разброс упругих свойств матрицы. Изменяя коэффициенты вариации модулей сдвига и объемного сжатия в диапазоне от 0 до 40%, приводятся оценки макроскопических модулей высокопористых композитов.

Работа выполнена при поддержке Российского фонда фундаментальных исследований (грант РФФИУрал№ 02-01-96-403)

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

1.Макарова Е.Ю. Нелинейные модели прогнозирования эффективных упругих свойств высокопористых композитов// Аэрокосмическая техника и высокие технологии - 2003. Материалы VI Всероссийской конференции. - Пермь: ПГТУ, 2003. - С. 99.

АНАЛИЗ ПОЖАРООПАСНОСТИ ИЗДЕЛИЙ С ГОРЮЧИМИ МАТЕРИАЛАМИ ПРИ ВНЕШНЕМ ПОЖАРЕ

Мальцев И.С, Одинцов Ю. I

ФГУП «Научно-исследовательский институт полимерных материалов», г. Пермь

Рассматривается осесимметричный многослойный полый цилиндр, состоящий из металлической оболочки (диаметром 500 мм, толщиной 3 мм), теплозащитно-крепящего слоя (резины толщиной 1 мм) и полимерного блока из горючего материала (диаметр канала 100 мм), подвергаемый в аварийной ситуации внешнему пожару с температурой пламени Тс = 1000 °С. Требуется оценить через какое время после возникновения пожара возможно воспламенение полимерного блока.

Для расчета теплового состояния изделия с учетом кинетики тепловыделения при термическом разложении горючего материала используется система дифференциальных уравнений [1]:

 

PiC, dT,/dT = X.i d2T,/dr2 +Vr-dT,/dr + Q- d^/dx;

(1.1)

 

PiCi-dTj/dx = ki d^i/dr2 + Vr-dT/dr;

(1.2)

 

dr\l dx = Z • exp(-Efl/(R(T+273))) x (1 - TI)(TI0 + rj),

(1.3)

при r0< r < R ,

0 < r j < l ,

1 = 2,3.

 

Здесь Z = 6,86T07 1/c - предэкспоненциальный множитель; Efl = 26218,3 кал/моль - эффективная энергия

активации; Q =

1065,7 кал/г - тепловой эффект реакции термического разложения; х)а = 0,2

- фактор

автокаталитичности; R =

1,987 кал/(моль град) -

газовая постоянная; Т

-

температура

в

градусах

Кельвина.

 

 

 

 

 

 

 

 

 

Система уравнений (1.1 - 1.3) дополняется начальными

 

 

 

 

 

 

Т(г, 0) = То = 20°С;

л0*- 0) = ло = 0;

 

 

 

(1.4)

и граничными условиями на поверхности оболочки

 

 

 

 

 

 

Х3 dT3(R, т)/ dr

=

а(Тс -

Т0).

 

 

 

(1.5)

Приведенная

система

уравнений (1.1

-

1.3)

с начальными (1.4)

и

граничными

условиями

(1.5)аппроксимирована системой уравнений в конечных разностях с использованием четырехточечной неявкой схемы и решена методом «прогонки». Для ее решения использовались следующие

теплофизические характеристики материалов: pj = 1,8 г/см3; Cj = 0,28 кал/(г*град);

Х\

1,36*10

кал/(см*с*град) - для полимерного горючего;

р2 = 1,195 г/см3;

с2 = 0,334

кал/(г*град);

Х2 = 7,8810 4

кал/(см с град) — для резины;

р3

= 7,8 г/см3;

с3 = 0,105

кал/(гград);

Х.3

— 6,05*10

кал/(см*с*град) - для оболочки.

 

 

 

 

 

 

Коэффициент теплообмена, вычисленный с помощью критерия Нуссельта Nu = ot*2RA.m= f(Prm*Grm) по формуле

а= 0,135(р-(2R)3gPrm(Tc - T0)/v2m)0'333-^/(2R)

сиспользованием теплофизических характеристик дымовых газов (индекс ш), для рассматриваемой

задачи составил а = 20,35 -Ю"4 кал/(см2-с-град)

В результате расчета установлено, что после возникновения внешнего пожара при заданных условиях время начала термического разложения горючего материала (до воспламенения) составляет х = 3,4 мин.

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

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

1 Амосов А.П. Тепловая теория воспламенения. Издание КуАИ, г. Куйбышев, 1982. - 96с.

УСТРОЙСТВО ЗАПОРНОЕ ОГНЕСТОЙКОЕ - УЗО (СПОСОБ ПЕРЕКРЫТИЯ ТРУБОПРОВОДА)

Марихин А.В.

Пермский государственный технический университет

Предлагается способ для перекрытия газового потока в трубопроводах например, при их ремонте. По предлагаемому способу в отверстие стенки трубопровода вводят гибкий мешок на раздвижной

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

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

Предлагаемый способ позволяет надёжно удерживать надутый мешок внутри трубопровода в неподвижном положении при возможном увеличении давления газа в перекрытой полости. Так как часть нагрузки ложиться на рамку, надуваемый мешок делается короче, чем мешок, используемый в прототипе. Это позволяет значительно уменьшить его вес, стоимость и время, необходимое для перекрытия и ремонта трубопровода.

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

В настоящее время изготавливается опытный образец на предлагаемое устройство и стенд для проведения испытаний.

Предлагаемое устройство имеет несколько типоразмеров на газопровод диаметром от 530 до 1420 мм.

Основные технические параметры устройства

приведены в табл 1.

 

 

 

 

 

 

 

 

 

Табл. 1

Наименование параметра и единица измерения

 

Тип, величина

 

УЗО

 

УЗО

УЗО

УЗО

УЗО

 

531

721

1021

1221

1420

Диаметр перекрываемого трубопровода, мм

530

720

1020

1220

1420

Длина камеры, мм, не более

350

350

450

450

450

Давление в камере, мм вод.ст.

5000

5000

5000

5000

5000

Масса камеры, кг, не более

3.5

4.5

8.5

10.7

13.2

Масса рамки, кг, не более

4.7

5.0

5.7

6.0

6.4

КВОПРОСУ ПОСТРОЕНИИ НЕСТАНДАРТНЫХ АВТОМАТИЗИРОВАННЫХ СИСТЕМ

СПРИМЕНЕНИЕМ СТАНДАРТНЫХ ПРОГРАММНЫХ И АППАРАТНЫХ СРЕДСТВ

Масков А.В., Сальников А. Ф.

Пермский государственный технический университет

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

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

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

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

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

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

ЛЕВЫЙ ИНВАРИАНТ РИМАНА В ГРАНИЧНЫХ УСЛОВИЯХ м е т о д а КРУПНЫХ ЧАСТИЦ

Матюнин В.П.

Пермский государственный технический университет

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

инвариант Римана R_ (м/с), рассчитываемый в /-ой входной частице. В этом случае скорость С”х (м/с) в

фиктивных частицах определяется исходя из условия постоянства температуры (энтальпии) по рассчитанному значению инварианта Римана:

С", = — R

k + 1Г

А

.•Л

1 - J 1 -

1 - — л

 

" * + 1

к - 1

к - \ R lj

 

При сверхзвуковой скорости противотока в фиктивной частице С"_х < 0 (отрицательный дискриминант в

(1)) в них следует записывать отрицательные значения критической скорости звука. Статические температура Т (К) и давление р (Па) рассчитывают из условия постоянства параметров торможения. Условие протекания при использовании инварианта Римана обеспечивается автоматически.

Рисунок 1. Изменение скорости при прохождении прямой и отраженной ударных волн через первую расчетную частицу

На Рис. 1 представлено решение модельной одномерной задачи движения ударной волны и отражения ее от препятствия (для скорости за прямой ударной волной Л/=0,47). Лагранжев этап имеет первый порядок точности. Без релаксации на этом этапе (использование скоростей с текущего временного слоя - штриховая линия на Рис. 1) применение левого инварианта Римана обеспечивает хорошую воспроизводимость и регулярность решения для многократных циклов движения. При этом в отличие от обычной постановки граничных условий не наблюдается затухания колебаний, обуславливаемых вязкостными эффектами разностной схемы. При наличии релаксации на лагранжевом этапе амплитуда колебаний хотя и уменьшается (сплошная линия), но остается выше, чем при постоянстве только параметров торможения (пунктирная линия).

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

Модорский В.Я., Соколкин Ю.В., Кашина И.А. Пермский государственный технический университет

Изучается проблема непрогнозируемого возникновения колебательных газодинамических процессов в РДТТ.

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

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

Численно-аналитические модели:

-описывают конструкцию РДТТ с заполнителем канального типа как бесконечный полый цилиндр, подкрепленный снаружи оболочкой;

-описывают конструкцию РДТТ торцевого типа как защемленный стержень, скрепленный с газовой полостью - трубкой.

Численные модели:

-описывают конструкцию РДТТ торцевого типа как защемленный стержень, скрепленный с газовой полостью - трубкой.

-описывают конструкцию РДТТ в осесимметричной двумерной постановке.

-описывают конструкцию РДТТ в трехмерной постановке.

Материал конструкции - упругий; нагрузка - газодинамическая.

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

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

Использование классических форм записи граничных условий не позволяет вскрыть причин аномального поведения РДТТ. Предлагается подход, реализующий уточненную постановку граничных условий с учетом жесткости нагружающей системы.

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

Разработана методика определения приведенной жесткости системы.

Вводится критерий, характеризующий степень влияния газа на напряженно-деформированное состояние конструкции S, который определяется как отношение а гтахжнсК а гтвх, где ст гтюжнсмаксимальные нормальные напряжения в направлении перпендикулярном контактной границе, определяемые в рамках унифицированного подхода, а ст Гшах" максимальные нормальные напряжения в направлении перпендикулярном контактной границе, определяемые без учета взаимовлияния в системе « газ-конструкция».

При этом выявлено, что:

1.Снижение жесткости конструкции приведет к усилению влияния газа на величину компонент динамического напряженно-деформированного состояния конструкции.

2.Снижение приведенной жесткости системы ведет к усилению колебаний газа и конструкции в РДТТ.

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

4.Совместный анализ динамического поведения газа и конструкции позволяет выявить условия

изменения амплитудно-частотных характеристик РДТТ.

_____

5.Сформулированы условия изменения амплитудно-частотных характеристик РДТТ.

6.Выявлены области изменения параметров газа и конструкции при которых возникает усиление амплитуд колебаний газа и конструкции в РДТТ.

7.Показаны области применения традиционных подходов.

Нихамкин М.Ш Пермский государственный технический университет

Развитие трещин в дисках газотурбинных двигателей (ГТД) в реальных эксплуатационных условиях необходимо рассматривать как случайный процесс. Этот процесс определяется действием нескольких групп случайных факторов. Первая группа - рассеяние свойств материала, в частности - характеристик трещиностойкости, чувствительных к влиянию состояния поверхности, остаточных напряжений, макро- и микроструктуры материала. Вследствие причин металлургического и технологического характера эти факторы могут быть различными у различных экземпляров дисков. Вторая группа - вариация от полета к полету условий нагружения диска, обусловленная влиянием метеоусловий, высоты, продолжительность полета, загрузки самолета. Третья группа рассеяние размеров дисков в пределах допусков на изготовление, определяющее рассеяние полей напряжений. Наконец последняя группа - форма, размеры и расположение исходного дефекта, инициализирующего процесс разрушения. В связи с этим, характеристики живучести диска также являются случайными величинами и их определение должно сопровождаться анализом закона распределения.

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

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

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

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

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

1.Нихамкин М.Ш. Закономерности развития поверхностных трещин в дисках газотурбинных двигателей. Техника машиностроения. - 2000. - № 5.- С. 15-17.

2.Нихамкин М.Ш, Нихамкина Е.М., Аторкин И.А. Статистическое моделирование живучести дисков газотурбинного двигателя. Вестник ПГТУ. Аэрокосмическая техника. Вып. 6, Пермь, 2000. С. 6269.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]