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

1130

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

1.

Москвитин В. В. Об одной нелинейной модели вязко-упругой среды, учитыва­

ющей

влияние вида напряженного состояния. — Механика

полимеров

1969 № 6

с. 994— 1001.

 

 

 

2.

Победря Б. Е. Математическая теория нелинейной вязко-упругости. __ В кн.:

Упругость и неупругость. М., 1973, с. 95— 173.

 

 

3..Findley

W., Reed A., Stern Р. Hydorstatic creep of solid

plastics. —

Mod. Plast.,

1968, vol. 45,

N 8, p. 144.

 

 

4.Findley W., Reed A. Effect of cross-linking on hydrostatic creep of epoxy — Polym. Eng. and Sci., 1977, vol. 17, N 12, p. 837—841.

5.Гольдман A. Цыганков С. А., Григорян Э. С. Прогнозирование объемной ползучести сетчатых полимеров с различной степенью отверждения. — Механика композитных материалов, 1979, № 4, с. 586—593.

6. Гольдман А. Я., Цыганков С. А. Об учете влияния вида напряженного состояния при описании нелинейных вязкоупругих свойств частично кристаллических полимеров (политетрафторэтилен). — Проблемы прочности, 1978, № 8, с. 60—64.

Охтинское научно-производственное

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

объединение «Пластполимер» , Ленинград

Механика композитных материалов,

 

 

1980, № 4, с. 733—737

УДК 539.3.001:678.067.5

Л. П. Хорошун, Н. С. Солтанов

СТАЦИОНАРНОЕ ТЕМПЕРАТУРНОЕ ПОЛЕ ДВУХКОМПОНЕНТНОЙ ПОЛОСЫ СЛОИСТОЙ СТРУКТУРЫ (КРАЕВОЙ ЭФФЕКТ)

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

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

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

 

 

 

 

{kTti) ti= 0

 

 

(1)

при

граничных

условиях на

границе L области

S = { — °°< X i< + оо; 0^л.'2^ а } : Т\ь =

=

Здесь k(x),

Т(х),

являются

случайными функциями

координат

плоскости

л'= (JCI, лг2) , причем Х\ перпендикулярна к слоям, х2 — параллельна слоям. Случайные

функции k(x),

Т(х)

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

суммы

математических

ожидании

и флюкту­

аций:

 

 

О

О

 

 

 

 

 

 

 

(2)

 

 

 

 

k= (k )+k;

Т = (Г )+ Т .

 

 

 

 

(k)T,u = — (kT,i),i.

 

(3)

Применяя к (3) известные формулы Грина

[6],

получим интегральное уравнение отно­

сительно Т {х):

 

 

 

 

 

 

 

 

Т(х) = - \q>(x')G ,i'{x,x')n'idl'--------- f Gti' (х, x')k(x')T ,i' (x')dS',

(4)

 

7

 

 

-<A>-

я

 

 

где tii

является внешней

нормалью к границе L,

a G (x,x') определяется формулой [6]

 

 

00

 

 

 

 

 

 

 

1

*

П Л Х 2

плх' 2

 

 

 

G(х, х') ——

— sin---------- sin ---------- exp

 

 

 

п

П=1 п

а

и

 

 

 

Для исследования краевого эффекта достаточно найти градиенты температур фаз.

Решение этой задачи будем строить методом условных моментов [7].

 

 

2.

Принимаем,

что k(x), Т(х), <p(L) являются статистически однородными функ­

циями координаты Xi и удовлетворяют свойству эргодичности, причем k(x) и

 

принимают постоянные значения

ki, ф0г, q>ai

(i= l, 2) в пределах

i-й фазы. Средние

температуры Т{ и Т2 в фазах 1

и 2 соответственно зависят только

от х2, т. е.:

 

 

 

Ti = Ti(X2)\

Т2=Т2(Х2).

 

(5)

Учитывая сказанное, умножим уравнение (4) на условную плотность f(T, T',k,k'lvx) (плотность распределения температуры и коэффициента теплопроводности в точке х' при условии, что в точке х находится фаза v) и проведем статистическое усреднение. В результате получим систему интегральных уравнений относительно условных мате­ матических ожиданий температуры и градиента температуры:

2

 

T v (x ) = -

J

Pvm{x, ^ ) < P m v (X, x')G,i' U , х')п'idl' -

 

 

 

 

L

m = 1

 

 

 

 

 

 

,

 

 

 

( 6)

------ —

f G ,i'(x,x ')

(km- (k ))T ,i'mv(x,x') pvm(x,n')dS'

{v =1,2}.

 

(k>

t

 

 

 

 

 

 

s

 

m = I

 

 

 

/

 

X^X

)

— математическое ожидание температуры в точке а: 'при

Здесь Tm*(x, x') = (T{x') |

условии, что в этой

точке находится фаза т , а в точке х фаза v; pVm —х')

)~

вероятность нахождения в точке х' фазы т , при условии, что в

точке х находится

фаза v; T,imv(x, х'), <p”lv(дс, х')

определяются

аналогично.

 

 

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

двухточечные моменты T , i mv,

<p,nv вырождаются в одноточечные (Г,jm), (фт ) и

урав­

нения (6) приводятся к виду

 

 

 

 

 

 

 

 

2

 

 

 

 

 

Tv ( x ) = ~

j G , i , ( x , x ' )

Pvm(x, x')(pm(x')n,,dr —

 

 

 

 

L

m=1

 

 

 

 

 

 

 

2

 

 

(7)

 

 

 

 

 

 

 

“ T IT

S G>{'( X>X')

У .

>x')dS'

(v= 1,2).

 

'

s

 

m = 1

 

 

 

Условные вероятности pvm(x, x') возьмем в виде

 

 

 

 

Pvm = P m + (—l) v+m( l —pv) -exp( — Y |A:I X 'I |) ,

 

( 8 )

Зависимость модуля упругости в продольном Е ц (а) и поперечном Е ± (б) направлениях одно­ направленного композитного материала от объемного соотношения волокно—связующее.

ннях. Упругие и другие деформационные свойства высокомодульных органических

волокон в продольном направлении

описаны в ряде работ [ 1—3], а

данные

по этим

же свойствам в трансверсальном

направлении крайне отрывочны.

Имеется

только

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

Упругие свойства однонаправленных армированных материалов в первом прибли­ жении можно определять по аддитивному вкладу компонентов [5, 6]. Однако в более общем случае математический аппарат для определения модуля упругости армирован­ ного материала на основе свойств компонентов достаточно сложен [7, 8] и включает не только характеристики свойств и объемную долю компонентов, но и ориентацию, размеры и другие характеристики армирующего наполнителя.

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

Вкачестве объекта исследования было взято высокомодульное волокно па основепараполиамида, описанное в работе [2], и эпоксидное связующее типа ЭДТ-10. Изме­ рение модуля упругости проводилось акустическим методом на приборе УК-22П на частоте 60 кГц при температуре 20 °С.

На рисунке показаны зависимости величин модулей упругости армированных ма­ териалов с различным соотношением компонентов. Как видно из этих данных, модули упругости материала в продольном и поперечном направлениях хорошо отвечают ад­ дитивной схеме. Линейная зависимость величины продольного модуля упругости компо­ зиции при экстраполяции до содержания волокна 100% хорошо соответствует экспери­

ментально определенному модулю волокна

= 125 •109 Н/м2. Аналогичная экстрапо­

ляция поперечного модуля упругости композиции дает величину поперечного модуля упругости волокна £ ^ = (6±0,4) •10э Н/м2. Для сравнения по методике, изложенной в работе [9], был проведен расчет £j_ из величины теоретического модуля волокна Е т имодуля специально полученной изотропной пленки Ех. Полученная в результате вели­ чина £j_ = 5,5- 109 Н/м2 находится в хорошем соответствии с экспериментально найден­ ной величиной.

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

 

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

 

1. Рысюк В. Д.,

Носов М. //. Механическая анизотропия

полимеров. Киев, 1978.

232 с.

 

 

2. Курземниекс

А. X. Деформативные свойства структуры

органических волокон

на основе параполиамидов. — Механика композитных материалов, 1979, № 1, с. 10— 14. 3. Перепелкин К. Е. Предельные механические свойства ориентированных поли­ мерных структур как армирующих наполнителей. — В кн.: Волокнистые и дисперси­

онно-изотропные композиционные материалы. М., 1976, с. 165— 171.

4.

Phoenix S. L., Skelton G. Transverse compressive moduli and yield behavior of

some

orthotropic, high-modulus filaments. — Text. Res. J., 1974, vol. 44, N 12, p. 934—940

5. Скудра A. M., Булаве Ф. Я. Структурная теория армированных материалов

Рига,

1978. 192 с.

6. Котосонова В. Я., Кузьмин А. М. Исследование упругих свойств композицион­ ных материалов медь—углеродное волокно. — Физика и химия обработки материалов

1976, № 2, с. 110— 112.

7. Kudrna М., Kamtner H.-W., Berger W. Zur Ausbreitung von Schallimpulsen in binaren Polymerfolien aus Hochdruckpolyathylene und Polyathylenterephthalat bzw. Polycapro-

lactam. — Faserforsch. u. Textiltechnik, 1978, Bd 129,

N 11/12,

S. 669—675.

8. Шермергор T. Д. Теория упругости

микронеоднородных

сред. М., 1977. 400 с.

9. Утевский Л. Е. Оценка ориентации

проходных

и держащих нагрузку участков

в аморфных прослойках кристаллизующихся полимеров. — Высокомолекуляр. соедине­ ния. Сер. А, 1974, т. 16, № 10, с. 2339—2344.

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

Механика композитных материалов 1980, № 4, с. 740-742

УДК 539.4:678.067.5

В.В. Беляев, И. А. Буяков

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

ВЗОНЕ ФЛАНЦЕВ*

Сосуды давления, изготовляемые методом непрерывной намотки ленты из компо­ зитного материала, имеют в качестве конструктивных элементов полярные фланцы [1]. Опыт эксплуатации таких сосудов показал, что под действием внутреннего давления фланцы смещаются относительно оболочки. Существующие методики расчета цельнонамотанных сосудов этот факт не учитывают; при этом предполагается, что контактное давление равномерно распределено по первоначальной поверхности контакта фланца с оболочкой [2].

Предлагаемая работа посвящена определению контактного взаимодействия фланца с оболочкой. В первом приближении фланец считается абсолютно жестким телом (штампом). Таким образом, ставится задача о контакте оболочки со штампом.

1.

Известно, что корректное решение контактной задачи для тонкостенных элемен­

тов конструкций можно получить на основе теории типа Тимошенко [3]. Для

однослой­

ной оболочки вращения при осесимметричном нагружении разрешающая система нели­

нейных уравнений теории типа Тимошенко принимает вид [4]:

 

 

Ац,1=ф1 (N22— N ц) — £iQi3—pi+thp;

(hi)

 

Q i 3 , i = —

+

+

^ з + Р ;

(1-2)

 

A f 11,1 = ф 1 ( A f 22— М п ) +

Q l3 + К \\^\ —к2У\М22—ГП\\

U\li = —k\Ui+8,\\ — 1/2'9,l2i

 

 

U 3,l = klU\ —‘&l]

 

 

 

 

Y I ,I

= X I I + Vz^ieis2;

 

(1.3)

 

е22= ф | Ц 1 + ^2 и з;

X 2 2 = "ф iY i — */2 A2Y i 2I

 

 

 

N = /(e ;

^13=^44613.

 

( 1.4)

Здесь

N= |yVii/V22MiiM22|/; е=|ецв2 2 Х 11Х 2 2 Г

(штрих означает транспонирование);

 

Вп В i2 A\i

A i2

 

 

 

 

 

В22 А 12 А22

 

 

 

 

 

К--

— симметричная матрица жесткостей;

 

 

Du D\2

 

 

 

 

 

 

D22

 

 

 

 

 

Ni3=Qi3+Nnfti+k\Ei3Mn-,

(1.5)

Yi ='в,1+ ei3;

( 1.6)

р — нормальная составляющая контактного давления.

* Д о к л а д , п р е д с т а в л е н н ы й н а

I I I В с е с о ю з н ы й

с и м п о з и у м п о м е х а н и к е к о н с т р у к ц и й н з компо­

з и ц и о н н ы х м а т е р и а л о в ( Л с н и н а к а н ,

с е н т я б р ь 1979

г . ) .

По определению 'бч = дсю—*i — угол поворота касательной

к поверхности приве­

дения в момент контакта со штампом; тогда

 

fti,i = k \ o - k i = Kiu

(1.7)

где *ю, Аю — угловая координата и кривизна штампа.

Чтобы получить разрешающую систему уравнений контактной задачи, необходимо

разрешить (1.4) относительно ei3 с

учетом (1.5);

всюду

исключить у, с помощью

(1.6); исключить контактное давление

р из ( 1.1) с

помощью

( 1.2).

После подстановки (1.6) в (1.3) с учетом (1.7) получим соотношение, содержащее величину Qi3,i. Вследствие этого уравнение (1.2) исключается из системы дифферен­ циальных уравнений и служит для определения контактного давления.

Разрешающая система уравнений контактной задачи приобретает следующий вид:

10'lQl3,I= — ( & 1 ~ Ql3+ (Ф1 —^21Э’|)Л^22~^1+‘0'19з;

ftiWn.i-t- Q13,1+ &1613^ 11,1= (В 44k\Mn) (xu + VaAieo2—Хц) — ХцЛ^п+ ( 6 4 4 ,1 —^иМп) 613;

Мц,1=ф1 ( М 22 —Mn)+Qi3+'6,i(iVii —k2M22) к 2&\2М 22— т й Ы|,i = —&1 U3 +E1 1 —1/г19,12;

U*A =

k\U{ - f t u

Z2 2 = ^\U\ +

k 2Uz\

Х2 2

= Ф1 (^l + 6 1 3) ~ ’/2 ^ 2 (Oi + 6l3)2;

6 1 3 = ( 6 4 4 —

k \ M \ \ ) _ 1

( Q i 3 + ^ i ^ i i ) ;

N =

K E ',

p = q3 + Q\3,i + t y iQ i3 ~ k iN i\ — k 2N 22.

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

2. При решении задачи о нагружении баллона внутренним давлением с учетом смещения фланца возник вопрос о корректном задании условий перехода от области контакта к смежному участку оболочки. Если кромка фланца скруглена, то обеспечи­ вается непрерывность вектора решения контактной задачи; к этому добавляется непре­ рывность угла Yi. В противном случае на границе области контакта обеспечивается непрерывность компонент Nu, Ми, а также условие сплошности среды, т. е. непре­ рывность ии м3, Yi.

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

То, что дело обстоит именно таким образом, следует из решения контактной задачи теории упругости. Известно, что при вдавливании в упругое тело штампа, граница области контакта которого не изменяется в процессе нагружения, по пери­ метру области контакта возникают напряжения, устремляющиеся к бесконечности [5]; однако интеграл от напряжений по площади контакта — величина ограниченная. В данном случае со стороны штампа на оболочку по границе контакта действует сосре­

доточенная кольцевая

сила, нормальная к деформированной поверхности. По своей

сути эта сила является

эквивалентом «бесконечных» напряжений.

У отверстия граничные условия записываются в виде Мц = Мц = 0; 8221“ 1+ ^2«з=

= 822°,

где 622°

— заданная кольцевая деформация у отверстия.

Возможны

другие

варианты граничных условий.

 

 

На

другом

краю оболочки задаются три соответствующих

граничных

условия.

3. Задача решается методом последовательного уточнения зоны контакта, которая сначала устанавливается произвольно. Контроль точности решения можно вести, на­ пример, по величине силы Nn на правой границе оболочки (рис. 1 ).

Для определения напряженно-деформированного состояния сосуда давления с учетом смещения фланца авторами разработана универсальная программа на языке АЛГОЛ-60 применительно к ЭВМ БЭСМ-6. Уравнения линеаризованы по Ньютону, система линейных уравнений решается методом прогонки Абрамова [6].

На рис. 1 показана расчетная схема сосуда давления с фланцем, который плотно прилегает к оболочке до деформирования. Линейные размеры даны в миллиметрах;

р а с т я ­

I мм

Рис. 2.

Рис. 1. Рис. 3.

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

упругости и сдвига и коэффициент Пуассона

равны

соответственно: £ = 2 -1 0 4

кгс/мм2;

G13 = 7 . 102 кгс/мм2; v= 0,3. На оболочку и

штамп

действует внутреннее

давление

<7з= 0,1 кгс/мм2.

На рис. 2 дана эпюра контактного давления; на границе контакта приложена

сосредоточенная нормальная сила Q= — Qi3K-f-Qi3= — 0,053

кгс/мм2, действующая в том

же направлении, что и контактное давление (здесь Qi3K —

перерезывающая сила слева

от границы контакта, Q13 — справа от границы). Положительное направление кон­ тактного давления совпадает с направлением внутренней нормали к оболочке. Пункти­ ром показано статически эквивалентное «контактное» давление.

На рис. 3 приведены графики изменения кольцевой деформации при учете сме­ щения фланца (кривая 1) и без учета при условии, что контактное давление — равно­ мерное (кривая 2). Смещение штампа в осевом направлении составило 0,17 мм.

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

1.Росато Д. В., Грове К. С. Намотка стеклонитью. М., 1969. 310 с.

2.Образцов И. Ф., Васильев В. В., Бунаков В. А. Оптимальное армирование обо­ лочек вращения из композиционных материалов. М., 1977. 144 с.

3.Филоненко-Бородич М. М., Изюмов С. Н., Олисов Б. А., Кудрявцев И. Н., Мальгинов Л. И. Курс сопротивления материалов. Ч. 2. 4-е изд. М., 1956. 540 с.

4.Буянов И. А. Нелинейные уравнения теории типа Тимошенко многослойных анизотропных оболочек. — Механика композитных материалов, 1979, № 3, с. 501—507.

5.Тимошенко С. П., Гудьер Дж. Теория упругости. М., 1975. 575 с.

6.Абрамов А. А. О переносе граничных условий для систем линейных обыкно­ венных дифференциальных уравнении (вариант метода прогонки). — Журн. вычисл. математики и мат. физики, 1961, т. 1, № 3, с. 542—545.

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

Механика композитных материалов, 1980, Л® 4. С. 742-744

УДК 539.4:678.01.001

С. П. Левицкий, Н. В. Шакиров

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

Поведение газовых включений в полимерных и композитных материалах в динами­

ческих условиях

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

прочности

при

деформации

растяжения [1, 2], некоторыми проблемами технологии

[3] и т. п.

Как

известно [4,

5], одной из характерных особенностей многих текучих

д и с п е р с н ы х

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

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