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

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

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

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

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

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

1985

А. Н. МАШКИН, И. Н. ШАРДАКОВ, Н. А. ШЕВЕЛЕВ

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

С УЧЕТОМ ВРАЩЕНИЯ

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

В настоящей работе предлагается метод численного анализа пространственных форм и частот колебаний вращающихся осе­ симметричных тел. Учет всех инерционных сил, участвующих в колебательном движении, приводит к комплексной проблеме собственных значений, что уже само по себе представляет опре­ деленные трудности. Однако они оправданы — постановка за­ дачи лишена упрощений и позволяет получить ряд интересных результатов. Например, при определенных сочетаниях жестко­ сти конструкции, скорости вращения, размеров собственная частота колебаний для отдельных гармоник может отличаться до 40 % по отношению к невращающемуся диску. Наблюдается эффект перераспределения частот по сравнению со случаем (о = 0.

Для математической постановки задачи используем соотно­ шения линейной теории упругости в цилиндрической системе

координат и принцип возможных перемещений

[3], имеющий

вид

 

6 А = Ь А а +6ЛИ+6ЛК+6ЛС.

( 1)

В вариационном уравнении (1) учитываются работа внутренних напряжений

работа силы инерции

6ЛИ= — JfJpUfiUdu,

91

работа центробежных сил, связанная с полем перемещений II,

6АС=—j*JJ* рюх (wXU) бШи,

и работа сил Кариолиса

6ЛК= — Jj'J 2р (wxU) 6Ши.

В выражениях для работ обозначено: аи, ei7-— компоненты тен» зоров напряжений и деформаций; U — вектор перемещений в ци­ линдрической системе координат; со — вектор угловой скорости вращения (совпадает с осью z); р — плотность материала конструк­ ции. Геометрия конструкции позволяет применить для отыскания пространственных собственных форм колебаний следующее разло­ жение по угловой координате:

 

U = (Uc+U a)ef‘,

(2)

где Uc и Ua — комплексные, симметричные и

антисимметричные

относительно угловой

координаты формы колебаний; р — ком­

плексная собственная

частота колебаний.

 

Применение к вариационному уравнению (1) метода конечных элементов приводит к алгебраической проблеме отыскания собст­ венных значений следующей системы уравнений (согласно обозна­ чениям, принятым в [2]):

([К]+Ра [М]+ш2 [C]+(i)p [С,]) {6}=0,

(3)

причем матрицы [К] — жесткости, [М] —масс, [С] — обусловлен­ ная центробежными силами, симметричны. Матрица [Ск], обуслов­ ленная силами Кариолиса (в [2] аналогичную по природе матри­ цу называют матрицей демпфирования), несимметрична.

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

Указанная схема решения была реализована в пакете про­ грамм для ЭВМ на языке ФОРТРАН с использованием ком­ плексной арифметики. Необходимо отметить, что наличие матри­

цы демпфирования не позволяет найти приведенные частоты,

как

в случае консервативной задачи, полагая отношение р/£ = 1 .

Это

приводит к

возникновению больших порядков при вычислении

частотного определителя по методу

конечного элемента и различ­

ному вкладу

составляющих матриц

в (3), что делает решение не-

устойчивым.

Однако если ввести

о

замену переменных t = tQrx,

92

тогда d/df = fld/d? и d2/d/ 2 = Q2d2/d/2. Если теперь положить Й= ©, то в (3) получим

([К]+[С]+р2 [М]+р [CJ) {6} = 0.

(4 )

Матрицы [С], [М],

[Ск] имеют множитель Q= р©2/Е. Уравнение [4]

можно переписать

в виде

 

где

([К0]+Р2 (М]+?[Ск]){6}=0,

(5>

 

 

[K0]= [R ]+ [q .

Форма записи (5) позволяет исключить указанные сложности при численной реализации. При определении реальной частоты необ­ ходимо учесть, что р = рй.

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

1 . Определены собственные частоты и формы колебаний диска для случая консервативной задачи, содержащие все три компо­ ненты вектора перемещений (это возможно только при рассмот­ рении пространственной задачи).

2.Силы Кариолиса оказывают влияние именно на эти собст­ венные формы и частоты (см. п. 1 ).

3.При учете сил Кариолиса происходит разделение частот, соответствующих симметричным и антисимметричным формам при

со = 0 , может наблюдаться

 

перераспределение частот по номерам.

4.

Для диска

из стали

£ = 2-10и Па

при ©=10 000 об/мин,

изменение

частот составляет

0н - 3%

(по сравнению с

© = 0 )

и 0 -^13%

при

© = 50 000 об/мин.

Если

материал

диска

имеет

модуль порядка

~ 109 Па,

то соответствующие изменения состав­

ляют уже

-V. 26 и 40 %. Чтобы реально представить

эти измене­

ния, приведем следующие значения (все для

первой

гармоники):

© = 0

об/мин,

рс = ра = 3514,9

рад/с;

©=10000

об/мин,

р1 =

= 2526,6

рад/с;

р2 = 4435,9

рад/с;

© = 50 000

об/мин,

рх=

= 2093,8 рад/с,

р, = 5014,8

рад/с.

 

 

 

 

 

 

 

 

 

 

 

Л И Т Е Р А Т У Р А

 

 

 

 

 

 

1. Б е р е з и н

И. С.., Ж и д к о в

Н.

П.

Методы

вычислений. Т.

1.

М.:

Наука,

1966. 632 с.

О. Метод

конечных элементов

в технике. М.: Мир,

1975.

2. З е н к е в и ч

541 с.

Л у р ь е А.

И. Теория

упругости. М.: Наука,

1970.

490 с.

 

 

3.

 

 

93

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

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

ЦЕНТР

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

1985

Л. А. ГОЛОТИНА

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

Определенный класс применяемых в промышленности ком­ позиционных материалов характеризуется высоким объемным наполнением матрицы твердой фазой при большом различна в жесткости материалов матрицы и наполнителя. Специфические свойства и особенности в механическом поведении таких мате­ риалов обусловлены появлением и развитием микроповрежде­ ний при деформировании. Необратимые повреждения структуры происходят даже при малых деформациях. В работах [3, 4] приводятся модели механического поведения высоконаполненных композитов и исследуются особенности деформирования конструкций из таких материалов на примере кругового ци­ линдра при различных программах нагружения.

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

В настоящей работе исследовались процессы квазистатиче- •ского деформирования конструкций, представляющих собой длинный полый цилиндр с криволинейным каналом, скреплен­ ный с оболочкой, под действием циклических температурных неоднородных нестационарных полей 1 3 = 1 3 (*, у, t).

Для описания свойств композиционного материала цилиндра использовали физические уравнения в виде [4]

t

11 / р

о

•94

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

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

Л И Т Е Р А Т У Р А

1. Г о л о т и н а Л. А. Численный расчет напряженно-деформированного состояния вязкоупругого цилиндра с криволинейным каналом.— В сб.: Упру­

гое и вязко-упругое поведение

материалов и конструкций. Свердловск:

УНЦ АН СССР, 1980, с. 41—44.

напряженно-деформированного состояния

2. Г о л о т и н а Л. А. Расчет

конструкций из высоконаполненного композиционного материала.— Методы расчета изделий из высокоэластичных материалов: Тезисы докл. Всесоюзн. конф. Рига: РПИ, 1983, с. 65—66.

3. 3 е з и н Ю. П. Особенности расчета напряжений и деформаций в из­ делиях из наполненных эластомеров.— Там же, с. 61—62.

4. Т у р к и н В. Т„ Г о л о в а н о в В. П. Расчет напряженно-деформиро­ ванного состояния двухслойного цилиндра из полимерного материала с не­ линейными вязко-упругими свойствами.— Там же, с. 63—64.

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

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

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

1985

Ю. В. СОКОЛКЯН, Ю. и. СИРИН

УСТАЛОСТЬ КОМПОЗИТОВ: ОСОБЕННОСТИ ПРОЦЕССА РАЗРУШЕНИЯ И МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ

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

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

Однако до сих пор не существует комплексных моделей уста­ лостной долговечности, позволяющих прогнозировать, кроме пределов выносливости, изменение жесткости и предела проч­ ности после произвольного числа циклов нагружения при слож­ ном напряженном состоянии. Это, по-видимому, объясняется тем, что, как указывал академик Л. И. Седов [5], адекватное описание деформирования и разрушения материалов при цикли­ ческих нагрузках непосредственно связано с решением пробле­ мы сложного нагружения. Наиболее полные достигнутые ре­ зультаты, отражающие существующее положение в области теоретических исследований усталости однородных деформируе­ мых тел, приведены в монографии [3].

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

чески на основе минимального

числа опытов предсказывать

7 Заказ № 39

97

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

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

С самого начала должно быть признано, что волокнистые композиты обладают уникальным сочетанием свойств, главными из которых являются неоднородность, анизотропия, упрочнение за счет передачи нагрузки на волокна, наличие поверхностей раздела. Отметим также, что усталостное поведение композита во многом зависит от его типа. Мы выделим три типа: стекло­ пластики, металлокомпозиты (композиты с металлической матрицей) и углепластики, и рассмотрим особенности усталост­ ного разрушения материалов на основе обобщения данных, со­ держащихся в обзорных монографиях [2, 4, 7], в разделах книги [1], написанных М. Оуэном и Дж. Хэнкоком, и в ориги­ нальной статье [8], где приведены результаты испытаний сереб­ ра, упрочненного стальной проволокой. Там, где это не будет специально оговорено, подразумевается одноосное растяже­ ние— сжатие в направлении армирования композита.

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

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

При любой усталостной выносливости отслоение волокон от матрицы начинается при напряжениях, составляющих 25—30 % от разрушающих, а растрескивание в матрице — при напряже­ ниях 60—70 % от разрушающих. В стеклопластиках при пульси­ рующем растяжении расслаивание близко к насыщению на ран­ ней стадии испытаний, но оно снова увеличивается перед кон* цом испытаний. Растрескивание матрицы растет медленнее и следует квадратичной функции от отношения текущего числа циклов до разрушения.

98

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

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

Принято, что уменьшение модуля Юнга на 2 % свидетель­ ствует о начале процесса расслаивания в стеклопластиках при пульсирующем растяжении. При возникновении растрескивания смолы уменьшение модуля лежит в пределах 10 % от началь­

ного значения.

При

растрескивании смолы обнаружено на­

растающее снижение

прочности (максимальное снижение —

около 13 % независимо от формы цикла).

Интересно,

что наибольшее изменение модулей' упругости

наблюдается при кручении полимерных композиционных мате­ риалов на основе высокомодульных углеродных волокон. Жест­ кость композитов уменьшается до 50—60 % исходного значения после десяти циклов при уровне напряжений, составляющем около 65 % статической прочности при сдвиге.

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

кследующему заключению. Отношение усталостной прочности

кпределу прочности при статическом нагружении зависит от конфигурации и не зависит от объемного содержания упрочняю­ щих волокон. Это отношение при пульсирующем растяжении находится в пределах 0,25-^0,37 для стеклопластиков, 0,3-^0,72 для металлокомпозитов и составляет 80 % и более для совре­ менных углепластиков.

Хотя введение волокон, например в металлическую матрицу,

увеличивает усталостную выносливость материала, отношение

7*

99

предела усталости к пределу прочности для матрицы характе­ ризуется несколько большей величиной, чем для композита в целом. Так, отношение предела усталости при пульсирующем растяжении к пределу прочности при растяжении уменьшается от 0,68 для чистого серебра до 0,50 для серебра с 30 об.% воль­ фрама. Отсюда следует, что однонаправленное армирование более эффективно для одноосного статического растяжения, чем для циклического нагружения.

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

Влияние вида напряженного состояния на предел усталости.

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

При пульсирующем растяжении прочность материала падает с увеличением угла 0, образованного основным направлением волокна и направлением приложения нагрузки. Например, для стеклопластиков предел выносливости при N = 106 циклов умень­ шается в два-семь раз при 0 = 45° и в четыре-десять раз и более при 0= 90°

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

Естественно, что математическое описание такого многогран­ ного и сложного процесса, как усталость композитов, сопряжено с большими трудностями на этапе постановки задачи. Нами предлагается подход, связанный с расширением области приме­ нения структурно-феноменологических моделей и поставленных для них краевых задач механики деформирования и разрушения структурно-неоднородных тел [6], в частности для прогнозиро­ вания усталостного поведения материалов.

100