Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Деформирование и разрушение композитов

..pdf
Скачиваний:
9
Добавлен:
15.11.2022
Размер:
6.86 Mб
Скачать

типа на

k-u

уровне; N — число типов структурных

элементов на

k-м уровне;

р — число уровней рассмотрения; nt — степень струк­

турного

многообразия элементов i-ro типа

уровня

k с

позиции

m-го уровня.

 

{S, р,,

v},

где S —

Для

каждого уровня задается тройка

совокупность параметров, описывающих топологию уровня; р — совокупность параметров, обусловливающих физические свойства уровня; v — совокупность «соответственных» параметров, связан­ ных с состоянием системы на данном уровне [11].

Пусть (£»<*) (r(fe))e v — математический объект, характеризующий поврежденность Де-го уровня («соответственный» параметр состояния уровня k). Физический смысл понятия «повреждение» (но не поврежденности) для выбранных уровней описания иллюстри­ руется следующей схемой:

Атомно-молекулярные уровни (А—М)

Структурно-композитные уровни (С—К)

Феноменологически-кон- струкционные уровни (Ф -К )

А.Разрыв атомно-молекулярных связей, дисло­ цированные атомы, вакансии, трехмерные дефекты типа дефекта Френкеля и т. п.

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

А.Дефекты типа пор, трещины нормального от­ рыва, полное объемное (когезионное) разрушение элемента уровня и т. п.

Б. Дефекты типа скола или отслоения, сдвиговые трещины, разрушение адгезии структурного эле­ мента и т. п.

А. Отказ конструкции или отдельных элементов, магистральные макротрещины нормального отрыва

ит. п.

Б.Разрушение «связей» между элементами кон­ струкции, макроотслоения и т. п.

Разделение повреждений на два типа Л

и Б связано с выбо­

ром математического объекта для описания

поврежденности ЫкК

Повреждения типа А могут быть заданы концентрацией или

относительным

объемом разрушенных

элементов, количеством

разрывов атомно-молекулярных

связей,

количеством

вакансий

и т. п.

В качестве количественной меры повреждений

типа А

можно

взять

скаляр (шаровая

часть

тензора поврежденности

Sif*). Повреждения типа Б , как величины векторного типа,

41

будем

характеризовать девиаторной частью

тензора поврежден­

ное™

D a \ т. е.

 

 

lo(*)=S,w‘,+ D («.

(2)

Поскольку тензор поврежденности со^) относится к разряду соответственных параметров состояния, то в качестве «сжатых» характеристик (связь 1-го типа между уровнями [12]) поврежден­ ности уровня —1) употребляются функционалы вида

 

 

 

'

(

 

<0(‘)(Л*-Ч, /K‘))du(W,

(3)

 

 

 

Л

Д(А1) ( >

— ! ) )

 

 

 

 

 

где

/Ч*-1) — радиусы-векторы уровней

k

и (k — 1); dv^(А>—

элементарный

объем области

ДГА:>;

— промежуточная

область

(мезоэлемент) уровня k , вводимая с целью установления

связей

между

уровнями

и содержащая

статистический ансамбль струк­

турных элементов

 

Между

координатами

смежных уровней

имеет место соотношение

 

 

 

 

 

 

 

 

 

 

/-(ft-i)=r(*)/e(0<e<l).

 

(4)

В случае, когда

поврежденность

со**) образует случайное эрго-

дичное поле,

вместо (3) можно использовать

 

 

 

 

©<*-1) (/•(*-!)) =

( coi*> (IK*—1),

/■<*)) ) .

(5)

3. Процесс разрушения твердых тел многоуровневый, причем количество уровней рассмотрения р при построении математиче­ ской модели этого процесса определяется двумя факторами: сте­ пенью структурного разнообразия системы & (идеальный кри­ сталл, композиционный образец, конструкция) и степенью «сложности» траектории внешнего воздействия Ж (уровень и ха­ рактер силовой нагрузки, воздействие теплового или нейтрон­ ного облучения), т. е.

Р=*Р{&, Ж).

(6)

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

42

повреждения путем увеличения количества повреждений (про­ исходит локализация повреждений на уровне А — М).

На следующей стадии образуются трещины в элементах структуры композиционного уровня с локализацией главным образом в плоскостях, по которым соприкасаются структурные составляющие. Локализация с энергетической точки зрения — наиболее выгодный способ реагирования на силовой поток: необходимые изменения осуществляются путем перераспределе­ ния силового потока из напрягаемых и ослабленных мест в бо­ лее прочные. Устойчивость сопротивления структурного элемен­ та разрушению («катастрофическому» развитию процесса) зависит от жесткости нагружения элемента [3]. С позиции композиционного уровня этот процесс неравновесный, необрати­ мый. Однако при большом числе уровней р=р(£Р) в результа­ те осреднения «соответственных» параметров [12] с позиции феноменологического уровня он может быть принят за квазиравновесный.

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

4.При многоуровневом подходе возможны два варианта

исследования процесса деформирования — разрушения систе­ мы: с позиций макро- и микроуровней. На макроуровне харак­ теристический элемент уровня &Sk) (мезоэлемент) выступает в роли так называемого «черного ящика» (феноменологический подход). На микроуровне внутренняя структура «черного ящи­ ка» макроуровня в определенной степени раскрывается, однако в этом случае «черный ящик» опускается на микроуровень.

Таким образом, на макроуровне твердое деформируемое тело Q можно рассматривать как глобальную термодинамиче­ скую систему, состоящую из бесконечного числа локальных под­ систем А.

Допустим, что процесс изменения состояния системы задан в пространстве с обобщенными координатами: е — деформация, со — поврежденность, 5 — энтропия. Тогда приращение внутрен­ ней энергии U единицы массы системы можно записать в виде

суммы элементарных

работ

 

d U (е,

со,

S) = o’- • d c o + T d S .

(7)

Здесь а — напряжение;

R — потенциал поврежденности;

Т

абсолютная температура. Уравнение (7) является обобщенным уравнением Гиббса, справедливым как для обратимых, так и не­

43

обратимых процессов [2, 5]. Следует отметить, что второе сла­ гаемое в выражении (7) нашло применение для описания разных по своей природе процессов, где со трактуется как массовая или числовая концентрация [2, 5], или как структурная энтропия [1].

Используя (7), в линейном приближении находим

 

= С • • в—{—/С*- (а — РД7\

(8)

 

©. т

 

где K —d*FldedT, F(e,

со, 7 ) —свободная энергия

единицы массы.

Исходя из того, что нелинейность диаграммы сг~е связана

с накоплением повреждений, в ряде работ [3,

6J обсуждается

вопрос о нахождении

величин типа К, со из макроэксперимента.

Причем, как показано

в [3], уравнение типа (8) справедливо и для

тела с макроповреждением.

Однако определение эффективных коэффициентов С, р и опре­ деление поврежденности для структурно-неоднородных систем

целесообразно вести с позиции микроуровня

с применением свя­

занной

многоуровневой

задачи,

постановка

которой

дана

в

[10—12].

При

 

этом с точки зрения работоспособности конструк­

ции представляет

интерес не столько напряженно-деформирован­

ное,

сколько

поврежденное

состояние системы,

т. е. распределе­

ние

а>(г(°)) в Q.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Л И Т Е Р А Т У Р А

 

 

 

 

 

 

 

1. А л е к с е е в

 

Г. Н. Преобразование энергии. М.: Наука,

1966.

316

с.

2. Б а з а р о в

И. П. Термодинамика. М.:

Высшая

школа,

1983.

344

с.

ков;

3. О моделях

магистральных трещин в деформируемых телах /С. Д. Вол­

УПИ

им. С. М.

Кирова.

Свердловск,

1980.

54

с.

Рукопись

деп. в

ВИНИТИ, № 573-80. Деп.

 

 

 

 

 

 

механика

компо­

 

4. В о л к о в

С. Д„

С т а в р о в В. П. Статистическая

зитных материалов. Минск: Изд-во БГУ, 1978. 206 с.

 

теории

упругости.

5. Г о л ь д е н б л а т

И.

И.

Нелинейные

проблемы

М.: Наука, 1969. 336 с.

 

 

 

 

 

 

 

^

 

 

 

 

 

6. Д з ю б а

В. С. Оценка поврежденности

армированных пластиков при

механическом нагружении.— Проблемы

прочности, 1983,

8,

с. 73—76.

 

 

7. И л ь ю ш и н

А. А., П о б е д р я

Б. Е. Основы

математической

теории

термовязкоупругости. М.: Наука, 1970. 280 с.

 

условиях

ползучести.—

 

8. К а ч а н о в

Д. М. О

времени

разрушения в

Изв. АН СССР. Отделение техн. наук, 1958, № 8 с. 26—31.

 

 

 

 

 

 

9. К о п ь е в

И. М.,

О в ч и н с к и й

А. С.,

Б и л с а ч а е в Н. К. Модели­

рование на ЭЦВМ процессов разрушения композитов с дефектами прочно­ сти связи между компонентами.— В сб.: Разрушение композитных материа­ лов. Рига: Зинатне, 1979, с. 57—61.

10. Элементы механики деформируемых систем на основе многоуровне­

вого подхода: основные

принципы построения моделей /С. В. Мельников;

Ин-т мех. сплошн. сред

УНЦ АН СССР. Пермь, 1982. 12

с. Рукопись деп.

в ВИНИТИ, № 4367-82. Деп.

 

 

11. М е л ь н и к о в С. В. О многоуровневом иерархическом подходе в

структурной механике

композиционных материалов.— В

сб.:

Структурная

механика композиционных материалов. Свердловск: УНЦ

АН

СССР, 1983,

с. 26—30.

 

 

 

44

12. М е л ь н и к о в С. В. К анализу взаимоотношений величин на раз­ ных уровнях в структурной механике деформируемых систем.— В сб.: Струк­ турно-механическое исследование композиционных материалов и конструк­ ций. Свердловск: УНЦ АН СССР, 1984, с. 3—7.

13. Р а б о т н о в Ю. Н. Ползучесть элементов конструкций. М.: Наука,

1966. 752 с.

В. Р., С л у ц к е р А. И. Кинетическая теория прочности.—

14. Р е г е л ь

В сб.: Физика сегодня и завтра. Л.: Наука, 1973, с. 328.

15. Т а м у ж

В. П., К у к с е н к о

В. С. Микромеханика разрушения по­

лимерных материалов. Рига: Зинатне,

1978. 294 с.

АКАДЕМИЯ НАУК СССР

УРАЛЬСКИЙ НАУЧНЫЙ ЦЕНТР

ДЕФОРМИРОВАНИЕ ^РАЗРУШ ЕНИЕ КОМПОЗИТОВ

1985

В. В. МОШЕВ

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

ПРИ ОДНОВРЕМЕННОМ НАКОПЛЕНИИ В НЕЙ ПОВРЕЖДЕННОСТИ:

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

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

эффект Маллинза — Патрикеева

[1, 2], который внешне вос­

принимается

как размягчение — уменьшение

начального моду­

ля упругости

резины. Степень

поврежденности определяется

уровнем деформированного состояния.

геометрической и

Учет совместного проявления

в резинах

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

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

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

46

100 300
500 £,%

Рис. 1. «Стационарные» гистерезисные петли с\Ш1а резины при различных амплитудах.

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

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

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

Эксперимент подтверждает это предположение. Резиновую пластину с круглым отверстием подвергали двукратному растя­ жению. В каждом опыте измеряли максимальную деформацию в месте концентрации у экватора отверстия и на удалении. Отношение первой величины ко второй есть коэффициент кон­ центрации по деформации К.

Цель опыта заключалась в определении зависимости К от внешней окружающей деформации. Рис. 2 показывает, что ко-

А

f

3

 

Рис. 2. Зависимость коэффициента

 

концентрации деформаций от внеш­

2

ней деформации е,,, и деформации

 

у экватора отверстия е* (а) и от

 

внешней деформации при

повтор­

1

ном растяжении образца

(б).

47

Рис. 3. Кривая растяжения резинового образца (см. рис. 2).

Рис. 4. Кривые одноразового и повторных растяжений саженаполненного эластомера.

Рис. 5. Представление результатов на рис. 4 с помощью модели.

эффициент К становится переменной величиной: сначала он увеличивается от значения 3, характерного для малых дефор­ маций, до 3,5, а затем падает до 1,7. В линейной теории малых деформаций К является константой.

Такое поведение обусловлено физической нелинейностью ис­ следованной резины (рис. 3), которая при одноосном растяже­ нии дает кривую вначале с уменьшающимся (в области 80— 150%), а затем возрастающим модулем. Из рис. 2 видно, что рост К сверх трех обусловлен падением модуля в области эква­

тора (ек~100%)

по сравнению

с модулем на удалении (е ~

~ 5 0 % ). Когда

экваториальный

00

участок переходит в область

упрочнения (ек>200 %), деформации окружения попадают в

область уменьшенного модуля (е~100% ): К. начинает умень-

оо

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

На рис. 2,6 приведен график зависимости К от внешней де­ формации при повторном растяжении того же образца. Срав­ нивая эту зависимость с предыдущей, видим, что начальное возрастание К выражено сильнее (4 против 3,5), а последую­ щий спад слабее. Такое различие отражает наличие и неравно­ мерность необратимого изменения свойств, возникших в образ­ це при первом испытании. В наибольшей степени поврежденность накопилась (с соответствующим необратимым падением модуля) у экватора, и поэтому деформация в данной зоне по сравнению с первым опытом возросла.

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

48

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

равновесная степень поврежденности устанавливается не­ медленно и необратимо в соответствии с максимальной дефор­ мацией в данной точке;

степень поврежденности и соответствующий комплекс механических свойств для условий сложного нагружения явля­ ются функциями только максимального главного удлинения

резина несжимаема, т. е. ^1^ 3= 1 .

Для описания свойств резины был принят четырехпарамет­ рический, предложенный Чоглем [3], потенциал, который полу­ чил хорошее экспериментальное подтверждение на нескольких эластомерах до деформаций 500—700 %

W( I E) = (2Gln)IE+BI%,

где п, т, G, В — константы, — инвариант деформации, выра­ женный через главные удлинения:

/ f i = £ М - 1 ) / я . а=1

Расчеты показали, что структура потенциала для исходной резины и для резины с накопленной усталостью сохраняется неизменной, а накопление усталости отражается через изменение численных значений констант п, т, G, В. Так как накопленная усталость является функцией только Яш (по условию), то построе­ ние потенциала сводится к получению зависимости указанных констант от Хт. Упругий потенциал при этом принимает в общем случае вид

w

(Хт) ГЕ < 4

 

п(Ът)

В качестве типового материала для определения констант использовали саженаполненный бутадиенстирольный эластомер, свойства которого при одноосном растяжении иллюстрирует рис. 4 [4], где представлена кривая одноразового растяжения и кривые повторных растяжений, когда первые были доведены со­ ответственно до 200 и 300 %. Этот материал выбирали потому, что он, наряду с сильно выраженной физической нелинейно­ стью, обнаруживает значительную необратимую изменчивость свойств в связи с накоплением поврежденности.

Результаты свидетельствуют о том, что реально наблюдае­ мую изменчивость свойств можно описать, используя зависи­

мость только двух констант G и В от Ят - Числовые

значения

констант получились следующими:

 

/г= 0,8;

m= 8; G= 14— 1,828Я,т ;

 

В = (4,75 - 0.58ЯJ

(17,75 - 5,53Ьт + 0,0578 (Хт -

I)4).

4 Заказ № 39

49

Кривые одноразового и повторных растяжений, вычислен­ ные с помощью указанной модели (рис. 5), удовлетворительно воспроизводят опытные кривые рис. 4.

Для численного исследования выбирали трубу с толщиной стецки, в пять раз превышающей внутренний радиус. Трубу нагружали внутренним давлением. Удлинение вдоль оси трубы исключено, Л2=1 . Канал внутри такой трубы можно рассматри­ вать как полость, производящую осесимметричную концентра­ цию напряжений и деформаций (рис. 6). Рассматриваемая кон­ струкция позволяет выявить закономерности, свойственные не­ однородным силовым полям, без излишних математических трудностей: геометрическая форма и характер нагружения со­ храняют неизменными направления главных осей в ходе де­ формирования. Условие несжимаемости дает возможность вы­ числять деформированное состояние независимо от напряжен­ ного, когда известна деформация хотя бы в одном месте, например на внутреннем или наружном контуре трубы.

Не приводя метода решения, рассмотрим полученные ре­ зультаты. Вначале для учета влияния только геометрической и физической нелинейности исследовали НДС трубы, свойства материала которой описываются потенциалом, не накапливаю­ щим поврежденности, но в остальных отношениях совпадаю­ щим с основным исследуемым потенциалом. Сопротивление де­ формированию у этого материала в области малых деформа­ ций характеризуется модулем упругости около 3,7 МПа. При 30—>100 % наблюдается пологий участок с касательным модулем около 1,0 МПа и, наконец, участок с крутым подъемом (каса­ тельный модуль около 7,0 МПа). Поэтому можно ожидать, что неоднородное поле напряжений будет перестраиваться с изме­ нением нагрузки: при малых деформациях (доли процента) в трубе должно наблюдаться распределение, характерное для модуля Юнга 3,7 МПа. При деформациях на внутренней по­ верхности трубы порядка 50 % общая податливость трубы долж­ на резко возрасти вследствие падения модуля в приканаль­ ном слое, который всегда «принимает» на себя основное сопро­ тивление трубы.

Для числовой характеристики радиальной перестройки НДС с ростом внутреннего давления брали отношение окружного истинного напряжения на внутренней поверхности трубы к той же величине на внешней поверхности трубы (<7е)Вн/(сге)нар.

При малых деформациях (в Гуковой области) это отно-

50