Деформирование и разрушение композитов
..pdfтипа на |
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
Рис. 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, В — константы, 1Е— инвариант деформации, выра женный через главные удлинения:
/ 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