Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Механика композитных материалов 1 1982..pdf
Скачиваний:
13
Добавлен:
15.11.2022
Размер:
10.09 Mб
Скачать

Yn1= Pn1[ (a21,nai5,n—|0&25,7iCCil,n OC21,7i] (a 23,nOCi1,71

0&13,7i X

X a 2i,n)_1; Y^2 = PTl2[(a21»7iai5»n"“ a25»nail,n^ ^"ail»7tl (a23.n X

Xan.n-ais.naai.n)"1;

Yn3 = Pn3(a21inai5ln-a25lnanfn) X

^2 11)

X (a23,nCXnln‘-ai3,na21,n)“ 1; 6n1= (1 —OCis.nYn^'^lS.nPn1

X

X (aii,n)_1i

6n2=

( - a i 3,n»Yn2- 'a i5,nPn2) (an.n)

>

 

6n3=

( —ai3,nYn3—аю.пРп3) (an.n)"1,

 

 

где

a,j\,n = aj\,n + cij\,n{{)d4a4\,n+ cij\,n{])dsa5itn + aj2tn(l)C4Ci4\tn +

+ tfj2,7i(1 *£5^51,71 + #j3,71^^4^41,71+ #j3,71^ ^ 5051,71;

aj3,n = 0j3,n+ 0ji,7i(l,d4043,n + 0ji,7i(1,d5a53,n+ 0ji,7i(1 d6063,7i +

+ 0j2,7i(1 ^4043,7i+ 0j2,n(1 C5053,7i+ 0j2,7i(1^6063,n+ 0j3,7i(1 b4 X

^ j 2)

X 043,71 + Oj3,7i(1 65053,7!+ Oj3,n(1 J1ЬбОбЗ.п;

|Otj5,7i= Oj5,Ti+

 

+ 0л,п(1 ^4045,71+ 0ji,7i(1 ^ 5055,71+ Ол,71(1 ^бОб5,71+ Oj2,7i(1 X

 

X £4045,71 + Оj2,71^ £ 5055,71 + Оj2,71^^£бОб5,71, /

1» 2, 3.

 

Заметим, что выражения для aij.n1 и бп1^,

Ртг1Л, уп1к (/—1,3;

6 = 1,3;

у= 1,6) можно получить из выражений ац,п и

6тЛ

РтЛ УтЛ заменяя в

оj1,71, cij2,n и Oji>?i(1 знак на противоположный.

Обозначим определитель бесконечной системы алгебраических урав­ нений (2.7) через Д(А,). Тогда, пользуясь обычной процедурой, получим

характеристическое уравнение в форме

 

А(Х)=0.

(2.13)

В результате решения уравнения (2.13) получаем зависимость Х=Я(х). Ясно, что при решении уравнения (2.13) необходимо заменить бесконеч­ ный определитель конечным. Чтобы оправдать такую замену, докажем, что этот определитель является определителем нормального типа [п]. Следуя [И], для этого нужно доказать, что следующий ряд сходится:

+00 -foo

-J-oo

3

3

(Anvplk + B nvplh + Cnvplk+ D nvplh) j

 

М =

j

 

 

(2.14)

7i=0 p=0

v=0

fe=l

i=l

 

 

 

Асимптотический

анализ

выражений

(2.6), (2.9) — (2.13)

позволяет

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

 

 

 

 

+оо

+оо

+оо

 

 

M<consti

71=0

 

П21п(If*) /v (£(«) Км+п (S/Кр) X

 

 

 

р=0 v=0

 

(2.15)

X

[Av{Z,tx)]-'Kv+p(£txp)[ (р+ const2) КР(£«ч)]-';

 

 

consti>0;

const2> 0 ;

p = R l2R~l>2 .

 

При произвольных конечных значениях аргумента всегда справедливо [10} неравенство

| М Ь * ) | «

const3

(2.16)

 

п\

При конечном значении аргумента всегда можно так подобрать

const4, что будет выполняться соотношение

 

 

 

|Кп(£<х) | «const4

2

Y

(2.17)

 

 

| Ь * | /

 

 

 

 

 

 

 

Табл. 2

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

 

 

пространственного армирования

 

 

 

Схема

 

Количество

а х х хх

а ххуу

 

a Vzyz

a xxxx a xxyy

армиро­

 

направлений

avvvv

~~a xxzz

 

a xzxz

вания

. армирования

 

^a xyxy

по рис. 2

Н

^ N ____

a zzzz

a yyzz

 

a xyxy

1

 

3

0,0714

0,0084

0,3730

0,11

2

 

3

0,4046

0,3415

0,3730

9,35

 

 

 

0,4046

0,0084

0,3730

3

 

4

0,0714

0,0084

0,0399

6,57

 

0,5156

0,2305

0,0568

4

 

6

0,1793

0,0630

0,0721

1,68

5

 

6

0,1328

0,0391

0,0859

1,00

1 + 3

 

7

0,1269

0,0361

0,0892

0,91

1 + 4

 

9

0,1143

0,0298

0,0986

0,73

3 + 4

 

10

0,2390

0,0922

0,0651

2,54

1 + 3 + 4

 

13

0,1461

0,0457

0,0804

1,19

Примечания. £ а/£ с = Ю0; Е с= 1; va=0,2;

Vc = 0,4; |Хх=0,35;

=

х, у, z — оси

симметрии механических свойств композита.

чение со можно установить экспериментально из условия минимума сред­ неквадратичного относительного отклонения.

Были обработаны экспериментальные данные для шести компози­ тов [1]. Для этих материалов (трехмерно ортогонально армированные и гибридные материалы) было обнаружено, что численные значения со при использовании в (7) компонент тензора податливости больше, чем при использовании компонент тензора жесткости: соа>сол.

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

 

Л .» ,.. £ J

A&U

,

(8)

1V

7Ь— 1 L /

\

 

 

<п)

 

 

где J ] |

J F(n)dUn)=V) A a$y6=Ai%lialftlhyli6]

N — количество траекто-

" - 1 L C„)

рий армирования, различающихся уравнением траектории или материалом волокна; F^n) — поперечное сечение расчетного структурного эле­ мента в математической модели материала; d/(n) — элементарная длина расчетного криволинейного структурного элемента; V — общий объем повторяющегося элемента композита; Лг- ^ п) — компоненты тензора жесткости однонаправленно армированного расчетного структурного элемента в его главных осях симметрии свойств; Аа$уб{п) — те же компо­ ненты, но приведенные к выбранным осям композита. Величина F ^ d l^ определяет элементарный объем отрезка структурного элемента.

2 - 1939

17

Умножая знаменатель и числитель в (8) на общую длину данной траектории волокон Ь п\ в пределах повторяющегося элемента и исполь­

зуя зависимости ^(п)£(п)= ^(п)’> V(n)IV=V а(п)/У а! V a = Va(l)+ Va(2)+ • • •

---- Ь Уа(п), зависимость [8] удобно привести к следующему виду:

 

N

 

 

i4<xfiv6—

 

У а(п)

f -<») Л (n)

71=1

~ v T }

 

 

 

V )

где Уа(11 — объем арматуры n-й траектории армирования; Уа — общий объем арматуры в элементе композита.

Предложенный математический аппарат позволяет оценить деформа-

тивные характеристики композитов с весьма сложными пространственно­

криволинейными схемами армирования [5].

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

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

 

Принимая во внимание вязкоупругие свойства структурного элемента,

компоненты тензора деформаций Хг/П =е^<п> при x,ij =

Oij представим в

виде полинома возрастающей кратности:

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

6 ij(n) ( 0

= d ijk lG h l{t)

+ CLijhlm pGkl(t) Gmp ( t ) +

+ J

 

T ) Ghl (t) d% +

 

 

 

+ J J(iij& m p (f Ti,

T2) Ghl (TI )

(Т2)ЙТ1^Т2 +

 

(9)

 

 

 

0

0

 

 

 

 

 

 

 

rno

n(П*

ДП

 

Я(П

Л(Л

 

 

с

использованием

где

aijkiy

uijhimp , . . . , G'Hhu

aijMmp , . . . определяются

теорий

армирования

для однонаправленно

армированного композита

[3, 4]. Подставляя (1) и (9) в (3), получаем

 

 

 

 

 

 

 

 

Сар ( 0

[

a i j h l G y b ( t ) lia lj $ lh y l l6 [ l ( n)

+

 

 

 

 

 

 

 

M'S

n

 

 

 

 

 

 

 

 

+

ЯijhlmpGyb [ t ) CT|ri ( t ) ^ia^jplh y lldlm ^lpr]^71^ +

 

+

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

+

J CLijlil ( t

т) Оуб (T ) lialj$lh yll6C lT \l(n )

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

t

t

 

 

 

 

 

 

 

 

]

+

0

0

 

 

t — %2 ) Gyb (T I ) (Tfcn (T2) h (iljfilh y ll6 lm llp x \d x id x 2 \i^ +

J

J Q i]k lm p (t — Ti,

 

( 10)

Считая, что (10) пригодно для любого пути нагружения, в том числе и для постоянных во времени напряжений, имеем следующие зависимости для определения компонент тензора податливости пространственно ар­ мированного композита (при ava=const):

с

°apve(0. “

-----J dijkl(t — х) lialjfil-hykbdx'l |x<n>;

 

 

^

n

0

J

 

#aPv6£ri (^) —

n I Q'ijlilmplicLlj$lhyll6lm%,Lpri~{~

 

t

t

 

 

"HJ

J ^ ijh lm p ( t —Tj, t —T2) liali$W yU blm \lpx\dx\d% 2^

(X^n^.

0

 

0

 

 

Вязкоупругие свойства структурного элемента композита для вари­ анта усреднения (6) можно представить аналогичным образом [4].

В качестве примера рассмотрим ползучесть композитов с полимерной эпоксидной матрицей и борными волокнами при однонаправленном и пространственно-хаотическом армировании с объемным коэффициентом армирования piz= 0,4. Вязкоупругие свойства связующего описываются с использованием дробно-экспоненциального ядра Работнова, а для струк­ турного элемента было принято, что деформирование композита во вре­ мени по разным направлениям описывается одной и той же нормирован­ ной функцией времени. Кривые изменения некоторых компонент тензора податливости двух структур армирования во времени, с учетом того, что со в (7) одинаковы для мгновенного и длительного модуля упругости (соа= 0,8), показаны на рис. 3.

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

Примем, что пластическое деформирование структурного элемента композита в виде цилиндра можно описать инкрементальной ассоцииро­ ванной теорией пластичности с линейным изотропным упрочнением в следующем виде [7]:

h ± -( J L - + - ? ! —) df — если

0 и d f ^ 0;

2 ' d%ij

дкн '

 

deij^P

 

 

0 — в остальных случаях,

 

где dei/n)p — тензор пластических деформаций в структурном элементе;

d f= -J^ —d%hi\ f= p 2—R2 — функция нагружения для структурного эле-

OXftl

мента; р — квадратичная функ­ ция тензора (напряжений или деформаций); i?=ri + c; с — на­ чальный радиус поверхности

Рис. 3. Изменение во времени некоторых компонент тензора податливости одно­ направленно (Д 4У5) и пространственно хаотически армированного (2, 3) боропластика при p,s = 0,4 и o)a= 0,8. 1

#111Ь 2 Qxyxy] 3 й хххх\ 4 — #1212;

5 — #2222-

пластичности для структурного элемента; ti — его приращение при на­

гружении; h=lJ{zop2) — функция упрочнения; z0=const.

В случае пропорционального активного нагружения, когда xa = kijt,

kij = const,

0, приходим к теории

деформационного типа:

 

 

si}MP = Cijhi{n)kitiin)(t - t*)

=Сцы(п){кы-к*ы).

(11)

Если хц = CTij, то полные деформации ег/ п) для структурного элемента на­

ходим, как сумму его упругих и

пластических деформаций:

8i/”)=

= ец(п)е (ом) +ец(п)р. Если же хц=гц, то напряжения

в структурном

элементе определяются функцией

 

 

 

сгij(n) = Oij(n>(еы{п) —гиМр)

 

(12)

В случае линейной упругости вместо (11) и (12) имеем

 

 

вц(п)= а«ы<»)ам<п) + 8««п>*’;

оц<п)= А т (пЦеы'п)-еы<п)р).

 

После усреднения этих величин по всем структурным элементам, при­ ходим к искомым соотношениям композита:

e a p = J 8ij<n>/ia/ipM !)4

(13)

15

 

 

ИЛИ

 

 

(Tap=J

O ijm iJiM D d l,

(14)

s

 

 

где по-прежнему g= g(0,<p =

— единичный вектор по оси

1; Ц£) —

функция плотности распределения арматуры по направлениям.

В случае, когда число направлений армирования конечно, интегралы

(13), (14) вырождаются в конечные суммы:

 

eap= "7T“

СГаР=—

Oij^nHialj^V^n\

n

Va

n

Уравнение поверхности пластичности для композита имеет вид шах р2 с2 = 0.

Рассматриваемая теория пластичности композитов содержит, таким образом, следующие задаваемые величины: ацм{п) (или Лда<п>), Сцм{п\ К(1) и с. Функция Щ ) определяется структурой армирования, а значе­ ние константы с, регулирующей величину упругой области, должно вы-

Рис. 4. Экспериментальная кривая растяжения матрицы (У) и ее аппроксимация (2); экспериментальная (3) и расчетная (4) кривые поперечного растяжения композита при (Оа = 0,5.

Рис. 5. Сечения поверхностей пластичности: 1 — первичная, 2 — после простого на­ гружения; 3 — после сложного нагружения для трехмерно ортогонально армирован­ ного композита.

бираться на основе экспериментальных данных. Компоненты тензоров a,ijhi{n\ Aijki{n) и можно вычислить по известным свойствам арма­ туры и матрицы с использованием теории армирования [6].

В качестве примера на рис. 4 показаны аппроксимация эксперимен­ тальной кривой растяжения матрицы из алюминия и экспериментальная кривая растяжения композита (алюминий—волокно борсик) поперек во­ локон, а также ее теоретическая аппроксимация при соае = о)ар = 0,5.

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

Оптимизация структуры ортотропных цилиндрических оболочек при действии внешнего давления показала, что для сравнительно толстых оболочек (h/R 0,05) при р,2= 0,4 применение пространственной струк­ туры армирования по сравнению с плоской приводит к существенному увеличению критической нагрузки, причем наибольшие значения крити­ ческих нагрузок дают схемы, хорошо сопротивляющиеся сдвиговым нап­ ряжениям [8].

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

Третий ранг имеет тензор пьезоэлектрических модулей [10]. Напри­ мер, в [ 11] метод усреднения использован для определения пьезоэлек­ трических свойств костной ткани по известным свойствам коллагеновых волокон и их пространственного распределения в кости.

Второй ранг имеют такие тензоры, как тензор удельной электропро­ водности, тензор коэффициентов теплопроводности и др. [10]. В однона­ правленно армированном (в направлении 1) трансверсально-изотропном композите скорость теплового потока h прямо пропорциональна гра­ диенту температуры:

(15)

где i, / = 1,2,3; — тензор коэффициентов теплопроводности, который по известным теплотехническим свойствам арматуры и матрицы опреде­ ляется зависимостями работы [ 12].

Отметим, что kn = kn и в данном случае &i3 = &i2. Если (15) разрешить относительно градиента температуры, то получим

дТ

d x i

где rfj — тензор коэффициентов теплового сопротивления, определяется обращением тензора кц.

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

в целом:

N

(16)

С другой стороны, в предположении одинаковости скорости теплового потока во всех расчетных элементах тензор гарг композита в целом опре­

деляется как

-4-Е J

 

 

 

 

N

 

 

d lfi

 

 

 

 

 

Га&

 

 

 

ri}lnUial№Va™

(17)

Ка

п -1

L ( п )

Цп)

 

 

 

 

 

 

 

Отметим, что raPft= {ka^ )~ \ но та^ Ф г арг. Верхний индекс (k или г) ука­ зывает на объект, который подвергался усреднению. Для слоистого пластика суммарный коэффициент теплопроводности в направлении парал­ лельно слоям определяется по (16), а перпендикулярно слоям — по (17).

В общем случае использование описанного математического аппарата связано с большим объемом вычислительных работ и немыслимо без применения ЭВМ. Поэтому был разработан ряд алгоритмов, основные из которых переданы в Госфонд алгоритмов и программ [13, 14].

С П И С О К Л И Т Е Р А Т У Р Ы

1. Жигун И. Г П о л я к о в В. А. Сврйства пространственно армированных пласти­

ков. Рига, 1978.2 15 с.

2. Малмейстер А. К. Пластичность квазиизотропного тела. — В кн.: Вопросы ди­ намики и динамической прочности, 1956, N° 4, с. 37—48 (Рига).

3. Крегер А. Ф., Мелбардис Ю. Г. Определение деформируемости пространственно армированных композитов методом усреднения жесткостей. — Механика полимеров, 1978, N° 1, с. 3—8.

4.Крегер А. Ф., Тетере Г. А. Применение методов усреднения для определения вязкоупругих свойств пространственно армированных композитов. — Механика компо­ зитных материалов, 1979, N° 4, с. 617—624.

5.Крегер А. Ф. Определения деформативных свойств композитного материала, ар­ мированного пространственно-криволинейной арматурой. — Механика композитных материалов, 1979, N° 5, с. 790—793.

6.Крегерс А. ФТе т е ре Г. А. Определение упругопластических свойств пространст­ венно армированных композитов методом усреднения. — Механика композитных мате­ риалов, 1981, N° 1, с. 30—36.

7. Зилауц А. Ф., Крегерс А. Ф., Лагздиныи А. Ж Т е т е р е Г. А. Расчет упруго­ пластических деформаций композита при сложном нагружении. — Механика композит­

ных материалов,

1981, N° 6, с. 987—992.

 

8. Крегер А. Ф., Тетере Г. А. Оптимизация структуры пространственно армирован­

ных композитов

в задачах устойчивости. — Механика композитных материалов,

1979

№ 1, с. 79—85.

F

 

9. Мелбардис Ю. Г., Крегерс А. Ф. Аппроксимация поверхностей прочности транс­

версально-изотропного материала. — Механика композитных материалов, 1980,

N° 3

с.436—443.

10.Най Дж. Физические свойства кристаллов. М., 1960. 385 с.

11.Martin R. В. Theoretical analysis of the piezoelectric effect in bone. — J. Biomech., 1979, vol. 12, N 1, p. 55—63.

12.Composite Materials Workbook/Ed. by S. W. Tsai. 1977, N 3. 728 p. (AFML).

13. Крегер А. Ф., Мелбардис Ю. Г. Определение характеристик деформируемости

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

композитов. — Алгоритмы и программы,

1979

N° 5 (31), с. 74.

г г »

,

14. Крегер А. Ф. Алгоритм определения деформативных свойств гибридного компо­

зита, армированного пространственно криволинейной арматурой. — Алгоритмы и прог­ раммы, 1980, № 6 (38), с. 32. * *

Институт механики полимеров

Поступило в редакцию 15.09.81

АН Латвийской ССР, Рига

 

УДК 539.582:539.377:678.067

Ю. П. Зезин

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

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

Наполненные полимерные материалы нашли широкое применение в современной технике. Многофазные материалы с высокой степенью на­ полнения (до 70—80% по объему), содержащие высокомолекулярные соединения в качестве модификаторов, получили название высоконаполненных полимерных систем (ВПС).

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

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

a (ei + е2) = а (ei) + а (е2) ,

где а — напряжение; еь е2 — деформации, в общем виде зависящие от времени.

Второе условие линейности вязкоупругого тела при этом может быть справедливым (в пределах погрешности эксперимента и естественного разброса механических свойств материала):

сг(ае) = аа(е)

или выполняется в более общем виде [3] —

а(ае) = а па(е).

(1)

В работе [1] предложена, а в работах [3—5] получила дальнейшее развитие простая модель деформирования и разрушения ВПС, в которой предполагается, что релаксационные процессы в подобных материалах обусловлены последовательным разрушением микрообъемов материала (области с относительно однородным распределением деформаций поли­ мерных цепей). Для микрообъемов принимается справедливым единый критерий длительной прочности (критерий Бейли в работе [1], критерий Ильюшина в работах [3—5]). Связь между напряжениями и деформа­ циями в рамках этой модели устанавливается на основе простейших представлений о соотношениях между деформациями микрообъемов и деформацией материала в целом. При этом используется предположение об упругости (не обязательно линейной) полимерного связующего.

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

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

ведены в [5].

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

боте [3]:

 

сг(е) =E & 2nAJF[e(t)].

(2)

Здесь Е2, А, п — постоянные материала; функционал Не (01 определя­ ется соотношением

F [e(0 ] = [ I

]'/h

(3)

О

 

 

t — время; т — переменная итерирования; п, т — константы. Условие разрушения материала имеет вид [3]

AIF[B (Q ] = 1,

(4)

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

Все приведенные в работе экспериментальные данные получены при испытаниях на растяжение образцов в виде двойной лопатки с длиной рабочей части 35 мм и сечением 7X7 мм. Материал — ВПС на основе полиуретана НП-1 [5]. Испытания при различных температурах (от 217 до 313 К) проводили на модернизированной испытательной машине WPM-250 [6], а испытания при различных гидростатических давлениях (до 15,0 МПа) рабочей жидкости (силиконовое масло) — на установке МИД-160 [7].

Методика определения параметров определяющих соотношений по экспериментальным данным подробно изложена в работе [5], поэтому здесь приведем лишь их значения, принятые для расчета зависимостей между напряжениями и деформациями в образце при различных режи­ мах нагружения: А =0,69; £ 2= 1,47 МПа; т = 0,04; лг = 0,6; k = 2S.

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

ления Р при описании механических свойств ВПС с помощью соотноше- t

ний (2) — (4) введением модифицированного времени f: f = f a Ttp~l (Т, Р).

о

Из определения функции временного сдвига следует, что ат,р(Т8,Р 8) =1, где Т81 Р3 — температура и давление приведения. Очевидно, что при Р = Р8должно выполняться равенство ат,р(Т,Р3) = а т (Т), а при Т = Т3 ат,р{Т8у Р) = а р (£ ), где ат (Т) — функция температурно-временного сдвига [8], аР(Р) — функция баро-временного сдвига [9].

Отметим, что методы временных аналогий широко применяются для прогнозирования поведения полимерных материалов в различных усло­ виях нагружения [10]. Примеры использования баро-временного сдвига для оценки влияния гидростатического давления на ползучесть ряда по­ лимерных материалов приведены в монографии [И]. Известно также ис­

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

Для определяющих соотношений (2)— (4) вид функции ат>р (Г ,Р) можно установить по данным испытаний материала с постоянной ско-ж

ростью деформации. В

этом случае

соотношение

(2)

между напряже­

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

и временем

при условиях

приведения

(Т= Т8,

P = PS, а г ,р = 1) можно записать в виде

 

 

 

lg a (e )= lg {А0Е2) + п lg e+m lg f,

 

(5)

где AQ= A /(nkB) ; B = B(mk + 1, nk) — бета-функция [13].

Р ф Р 8,

При условиях, отличных от условий приведения

(Т ф Т 8,

ат,рф 1), соотношение

(2) для нагружения с постоянной скоростью де­

формации дает

 

 

 

 

 

lg [су(е) ] г,р= lg (А0Е2) + ttlg e + m lg f- m lg аг,р(7\ Р ).

(6)

При выводе (6) предполагалось, что температура и давление не меня­ лись в процессе испытания. Из (5) и (6) легко определить значение ат.р при фиксированных в опыте Т и Р:

lg ат,р$Т, Р) = m - 4 g {[от(е)] т ,р / су(е)}.

(7)

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

На рис. 1 представлены кривые растяжения материала НП-1, полу­ ченные с постоянной скоростью деформации 5,8-10~~4 с-1 при различных температурах (217, 233, 253, 273, 293 и 313 К). Результаты определения зависимости ат (Т) по соотношению (7) в исследованном интервале тем-

Puc. 7. Кривые растяжения материала НП-1 при постоянной скорости деформации 8=5,8- 10-4 с - 1 и при Г=313 (7); 293 (2) \ 273 (5); 253 (4); 233 (5) и 217 К (6).

Рис. 2. Кривые растяжения материала НП-1 при различных значениях гидростатиче­

ского давления Р = 0,1 (7); 1 (2); 2 (5); 4

(4); 8 (5) и 15,7 МПа (б). --------

экспери­

мент; ---------

расчет.

 

Т. к

(Год. МПа

lg ат

Р, МПа

CTQ.I * МПа

Igflp

313

0,220

-0,4826

0,0981

0,245

0

293

0,230

0

0,9810

0,310

2,5549

273

0,245

0,6859

1,9620

0,355

4,0266

253

0,265

1,5379

2,9430

0,380

4,7654

233

0,320

3,5855

3,9240

0,400

5,3224

217

0,440

7,0431

7,8480

0,470

7,0733

 

 

 

15,6960

0,480

7,3019

ператур приведены в табл. 1. Здесь <7o,i — напряжение при величине де­ формации 0,1. Температура приведения была выбрана равной 293 К. Экс­ периментальная зависимость йт(Т) аппроксимирована известным соот­ ношением Вильямса—Ланделла—Ферри (ВЛФ) [8] с коэффициентами Ci = —2,1 и С2 = 98: lg a T(T) = —2,1 —Ts)/(98 + T—Т8). Сопоставление экспериментальной и расчетной зависимостей ат(Т) представлены на

вставке рис. 1.

Функция аР(Р) при Т =Т 8 определена аналогично по серии кривых растяжения, полученных при скорости 0,895 с- 1 и при различных гидро­ статических давлениях (до 15,7 МПа), приведенных на рис. 2 сплош­ ными линиями. Результаты определения аР(Р) по соотношению (7) даны в табл. 2. Давление приведения выбрано равным атмосферному. Зави­ симость аР(Р) аппроксимировалась уравнением, аналогичным уравне­ нию ВЛФ: \gaP(P) = 8 ,5 (P —Ps)/ (22,4 + Р —Р8). Следует заметить, что функция баро-временного сдвига введена лишь для описания влияния внешнего гидростатического давления среды. Влияние вида напряжен­ ного состояния в рассматриваемой модели предполагается учесть введе­ нием зависимости параметра А от параметра Лоде.

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

На рис. 3 и 4 представлены примеры сопоставления эксперименталь­ ных и расчетных данных по влиянию температуры на механические свойства материала НП-1 в опытах на растяжение со сложной времен-

и

-НЫе (------

} "

Расчетные

(-------

) зависимости

между напря­

жением и

деформацией в

опыте

с разгрузкой и

повторным нагружением. 7=233 К-

 

 

е=5,2 • Ю-з

с-ь Л/7[е( /,) }=

1,17.

 

Рис. 4. Экспериментальная

(-------

) „

расчетная (---------

) зависимости

между напря-

жением и деформацией в опыте с остановкой. 7=253 К; ё = 5,2* 10-3 с-1; A/F[e(t )}=1 01.

ной программой нагружения при давлении приведения. На рис. 3 пока­ заны экспериментальное и расчетные зависимости между напряжением и деформацией в опыте с разгрузкой и повторным нагружением при тем­ пературе 233 К. Образец растягивали с постоянной скоростью деформа­ ции ё= 5,2 • 10-3 с-1 в течение 30 с, затем полностью ^разгружали и вновь нагружали до .разрушения с той же скоростью деформации. Как видно из рис. 3, принятые определяющие соотношения удовлетворительно опи­ сывают эффект Маллинза, проявляющийся в искажении кривой растя­ жения при повторном нагружении. Условие разрушения (4) выполняется при этом с точностью до 17%.

На рис. 4 приведены экспериментальная и расчетная зависимости нап­ ряжения от деформации в опыте с остановкой (см. вставку) при темпе­ ратуре 253 К. Сопоставление показывает, что и при таком режиме нагру­ жения принятые определяющие соотношения удовлетворительно описы­ вают зависимость напряжения ОТ -деформации. Условие разрушения в этом опыте удовлетворяется с точностью до 1 %. г

Использование функции баро-временной аналогии и условия разру­ шения (4) позволяет рассчитать зависимость предельной деформации е* от гидростатического давления для режима нагружения с постоянной

скоростью деформации:

 

в Л ^ ) - И о ^ а Р- ( Р ) ] 1/п..

(8)

На рис. 5 показана расчетная зависимость, полученная по соотноше­ нию (8). При давлениях более 3,92 МПа образцы не разрушались в про­ цессе испытания из-за ограниченного хода рабочего штока машины МИД-160.

Использование температурно-временной аналогии позволяет устано­ вить диапазон скоростей нагружения, в котором удовлетворяется усло­ вие разрушения ВПС (4). С этой целью на рис. 6 представлена экспери­ ментальная зависимость величины отношения A/F[b(t*)] от логарифма приведенной скорости деформации lg (eaT). Как видно из графика, в широком диапазоне скоростей деформаций, составляющем шесть поряд­ ков, это отношение равно единице с точностью ±15% , что говорит о вполне приемлемом для инженерных оценок удовлетворении критерию разрушения ВПС. Отметим, что параметры функционала, входящего в это соотношение, определены по данным исследований деформационных свойств материала.

Удовлетворительное совпадение экспериментальных и расчетных данных в значительном диапазоне температур и давлений может слу­ жить косвенным подтверждением справедливости гипотезы о едином ме-

Рис. 5. Зависимость предельной деформации

при

растяжении материала НП-1 от

гидростатического давления: • — эксперимент;

------- расчет; ё = 0,895 с-1.

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

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

С П И С О К Л И Т Е Р А Т У Р Ы

1. Fitzgerald J. Е., Farris R. J. Deficiencies of viscoelastic theories as applied to solid propellants. — Bull, of the first joint army-navy-NASA-air force meeting of the working group on mech. beh., 1970, publ. N 160.

2.Mullins L. Effect of stretching on the properties of rubber. — J. Rubber Res., 1947, vol. 16, N 12, p. 275—289.

3.Зезин Ю. П. О механическом поведении материалов, обладающих свойствами

незатухающей памяти. — В кн.: Физика структуры и свойств твердых тел, 1977, вып. 2,

с.136—141 (Куйбышев).

4.Трифонов В. П., Малинин Н. И. О связи между напряжениями и деформациями для полимерных материалов, проявляющих свойства незатухающей памяти. — В кн.:

Научн. тр. Ин-та механики МГУ. М., 1975, с. 77—85.

5.Зезин Ю. П., Малинин Н. И. О методах описания деформационных и прочност­ ных свойств высоконаполненных полимерных систем. — Механика композитных мате­ риалов, 1980, № 4, с. 592—600.

6.Зезин Ю. П., Малинин Н. И. Экспериментальная проверка концепции Фицдже­ ральда о незатухающей памяти наполненных полимеров. — Изв. АН СССР. Механика твердого тела, 1977, N° 3, с. 125—129.

7. Зезин Ю. П., Козырев Ю. И., Малинин Н. И., Никифорова В. М., Смир­ нов В. С. Машины для испытаний полимерных материалов, созданные в отделе пластич­ ности НИИМ при МГУ. — В кн.: Методы и приборы для механических испытаний по­ лимерных материалов. Ростов-на-Дону, 1979, с. 11()— 115.

8.Ферри Дж. Вязкоупругие свойства полимеров. ИЛ, М., 1963. 536 с.

9.Айнбиндер С. Б., Алксне К. И., Тюнина Э. Л. Свойства полимеров при высоких давлениях. М., 1973. 192 с.

10.Уржумцев Ю. С., Максимов Р. Д. Прогностика деформативности полимерных материалов. Рига, 1975. 416 с.

11.Голъдман А. Я. Прочность конструкционных пластмасс. Л., 1979. 320 с.

12.Filler R. W., Tschoegl N. W. Time-temperature pressure superposition in elasto­ mers. — In: Proc. 7th Int. Congr. Rheol. Gothenburg, 1976, p. 170—171.

13.Бейтмен Г., Эрдейн А. Высшие трансцендентные функции. М., 1965. 296 с.

Институт механики

Поступило в редакцию 31.03.81

Московского государственного университета

 

им. М. В. Ломоносова