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

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

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

приведенных для 0,001. Рис. 4—б демонстрирует различные обнаружен­ ные при расчетах тенденции. Все три значения Хр/а на рис. 4—б дают напряжение, колеблющееся от растягивающего к сжимающему и в ко­ нечном итоге стремящееся к нулю.

Случай задания граничного условия для температуры. Когда кон­ центрация влаги на поверхностях С* поддерживается постоянной, в то время как температура поверхностей внезапно возрастает от 7\ до Г/ и затем остается постоянной, напряжения могут быть записаны в виде:

 

 

 

Ох = О у = ~1 — Vp { Т f

Т { ) G T \

Та;у = 0>

 

 

 

 

где

ат=

a

[4'(D ,T) - 4 '( |,D ,T)]+

\ l - P + * s ]

J

[ 4 ( D n ) -

 

L

J

 

 

L

a

 

 

 

—т

я 2T)].

Рис.

5 демонстрирует

распределение

напряжений

на

раз­

 

 

 

личных

стадиях

процесса

диффузии

 

 

 

тепла в пластину. Для большинства

 

 

 

случаев

сопутствующий

 

переходный

 

 

 

процесс абсорбции влаги создает коле­

 

 

 

бательный характер изменения

напря­

 

 

 

жения. На рис. 5—а имеется исключе­

 

 

 

ние для Л,р/ос = 0;

в этом случае напря­

 

 

 

жение в центре пластины всегда растя­

 

 

 

гивающее

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

 

 

 

температуры), а

напряжение

на

по­

верхности — всегда сжимающее.

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

Рис. 5. Распределение напряжения по толщине пластины п различные моменты времени. Поверхностная концентрация влаги фиксирована; у = 0>25. и= 0,001 (б, в) и 0,01 (а), Аф/а = 0 (а); 0,1 (о) и 10 000 (б). Цифры у кривых — значения т/,.

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

Приложение: теория диффузии. Количество влаги на единицу объема в некоторой точке тела обозначается С, в то время как f — количество влаги, протекающей (в на­ правлении f) через единичную площадку, перпендикулярную к f, в единицу времени. Закон сохранения массы записывается в виде:

Vf + £ = 0 ,

где C=dC/dt. Аналогично, закон сохранения энергии:

 

Vq + pCpf = 0,

где рСР — количество тепла,

необходимое для поднятия температуры Г на один градус

в единице объема тела; q —

количество тепла, протекающего (в направлении q) через

единичную площадку, перпендикулярную к q, в единицу времени.

Термодинамическое рассмотрение необратимого процесса приводит к обобщениям закона Фика и закона Фурье соответственно для диффузии и теплопроводности:

f= - D mS7C-XDhVT- q = pCp [-D h V r -v D m V C ],

где Dm — коэффициент диффузии влаги; Dh=kjpCp — коэффициент диффузии тепла; k — коэффициент теплопроводности; v характеризует эффект Дюфо: X — эффект Соре термической диффузии. Эффект Дюфо заключается в существовании потока тепла в изо­ термическом теле при наличии неоднородного поля влажности. Отношение потока тепла к потоку влаги Q* = pCPv. В эффекте Соре отношение потока влаги к потоку тепла, когда влага распределена однородно, а температура — неоднородно, равно Х/рСр. Из соотношений Онзагера можно выразить X через Q*:

Х=

D,

1

Q*.

D,

v=-----

 

рСр

 

где R — газовая постоянная. Условия Dm> 0; Д л > 0; М><1 являются следствиями поло­ жительной определенности локального производства энтропии.

Сочетая уравнения для потоков с законами сохранения, приходим к основным диф­ ференциальным уравнениям теории:

DV2C= C-XT\ DW2T= T - \C ,

где D = ( \ — X\)Dm-, v = ( l —Xv)Dh. Подобные уравнения могут быть получены и для других теоретических моделей. Подробности в работе [4].

Авторы искренне благодарят AFOSR за поддержку работ по выпол­ нению контракта F 49620-78-С-0054, а также д-ра Джозефа Д. Моргана за помощь при проведении данных исследований.

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

1.Henry Р. S. Н. Diffusion in absorbing media. — Proc. R. soc. Ind. Sec. A, 1939,

vol. 171, p. 215—241.

2. Hartranft R. J., Sih G. C„ Chen T. S. Interaction of temperature and moisture in diffusion. — Lehigh univ. Inst, of fracture a solid mech. Report IFSM-77-82, August 1977.

3. Shen C.-H., Springer G. S. Moisture absorption and desorption of composite ma­

terials. — J. composite materials, 1976, vol. 10, p. 2—20.

4. Hartranft R. J., Sih G. C. The influence of the Soret and Dufour effects on the diffusion of heat and moisture in solids. — Lehigh univ. Inst, of fracture a solid mech.

Report IFSM-78-92, June 1978.

5. Pipes R. B., Vinson J. R., Chou T.-W On the hydrothermal response of laminated composite systems. — J. composite materials, 1976, vol. 10, p. 129— 148.

Институт разрушения и механики

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

твердого тела Лихайского университета,

 

Бетлехем, Пенсильвания, США

 

Потенциальная и кинетическая энергия деформации

имеют вид:

 

 

 

 

 

 

 

 

 

 

П=—J

{T\Ei + T2E2+ S(ii + M\K\-\-M2K2+ Hx)dxdy\

 

*4

 

 

 

 

 

 

 

1

 

где

г

du

(

w

dv \

d2w

d2w

 

Т

-,

г ,= С п Ж + с „ \

т + 1 ^, ) - к и —

- К а —

Т

г

ди<Г

I w

<dv \

d2w

d2w

;

Г2 =

 

+ С2ц

- +

~

; -

/(,2 —

- к 22 —

, ж

^

 

^

 

д2ш

+ 1< п ^ + К и ( ^ + ^ - ) - ,

All = —D u — —

Di2

<ty2

 

 

C?X2

 

 

 

dx

 

R

dy

M 2= - D 12

д2ДО

^

 

d2ay

 

 

„„du

 

( w

dv \

dx2 _ D ,,_ ^ T + * 12 dx

 

R +-d^l

 

 

 

 

H = - 2 D

 

d2w

 

 

f du

dv

\

 

 

00

dxdy

+ /(66 1

dx

/

 

 

 

 

 

 

' dy

 

 

 

du

 

 

w

 

dv

 

du

dv

 

 

61

dx

t2=J +

W '

 

— h -— ;

 

 

 

 

 

 

dy

dx

 

 

X\ =1

d2w

 

x2= —

 

d2w

т

о

d2w

 

 

M 4

 

 

dy2

 

dxdy

 

 

1

 

 

 

 

 

 

 

 

 

 

Перемещения оболочки представим в виде:

ОО00

оболочки

(1)

(2)

(3)

«(*, У, 0 = X j

X j C?An(t) *'ft(*) ■cos pnt/;

 

71= 0

ft= 1

a h

 

 

 

(*>*/> 0 = X J

X J

Vkn{i]Xk{x) sin p„y;

(4)

 

71

ft

 

 

ffl>(x, */,/)= XJ

X J

^ hn(t) **(*) C0SM ;

 

n

 

ft

 

 

где pn= /z/i? (л — число волн в окружном направлении); Хк(х) Л-я форма собственных изгибных колебаний балки при некоторых гранич­ ных условиях на концах (балочные функции, соответствующие ряду ти­ пичных случаев закрепления, приведены в [14]).

Подставив (4) в (1) —(3) и введя обозначения (штрихом отмечено дифференцирование по х)

L L L

IhmW= J XhXmdx;

Ihmw= | XkX"mdx\

IkmW = ) X \X 'mdx-,

о

о

о

 

L

 

/,„,<4>= 1 Xk"Xm"dx, 0

где

 

 

 

 

 

 

a hnu =

 

~

 

 

( C e e P A f c O ) + C x, / лл(4) ) ;

 

 

 

 

 

 

 

 

 

 

 

 

 

«л*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

ftnl2 =<Xk

 

(

С ia/

ftft<2>—

C 6e/ftfc<31) ;

 

 

 

 

аы13=—

 

[ ( - ~ - + K ^ n 2) I k ^ - 2 K ^ n 2I k k ^ - K u h k ^ \

 

 

 

 

as

1 '

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а„па =С2фпЧккт + С№1ккЮ\

 

 

 

 

 

 

 

a „ n23=

p

 

 

[

(

 

^

 

- + I^kk' "2- KP, 2/„,» .'2Ч)

-

2 /лй(=166/ч]

 

 

 

a*„33=

( - ^ - + 2 « 22- ^ - + D 22p„* ) /»„<•>-2 ( ^ - + D 12p„2)

/ ftli<2>+

 

 

 

 

 

 

 

 

 

+ 4£)ббРп2/ hft(3) +

i/hh(4);

 

 

 

 

 

 

 

 

2

 

n

k

 

 

 

 

 

 

V

 

 

 

 

#

 

“ *

7

 

Работа

Л

 

реакции

заполнителя

q = [^(1), q(2\

^(3)],

действующей

на

внутреннюю

поверхность

 

оболочки,

А = J (qu)dxdy= J‘ (qMu + q№v +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q

 

 

 

Q

 

 

 

+q<?'>w)dxdy. После того, как найдено решение задачи для заполнителя

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

теля на поверхности контакта

 

(этому будет посвящен п. 2), она может

быть записана в виде:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A = - n R

2

J

Z

I

e2i“knt (U ltn2sknu + Vhn2skna + Whn2shn33+

 

 

 

 

 

 

 

n

 

ft

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

Uhn VhnSkn12+

 

u hnWhnShn3+ Vkn Whnskn2* ) ,

 

( 9 )

Де

 

 

 

 

 

 

 

 

1

 

/

 

 

 

 

Oi

41

 

\

 

 

 

 

 

 

 

 

 

sknu = —

 

x

Qknuhh{3)+ - ^ r - I u h {5) I

 

 

 

 

 

 

 

 

 

 

 

 

Як2

 

 

 

 

aft2

 

 

 

 

 

 

 

 

 

 

 

 

 

O,

 

52

 

 

 

 

 

 

 

 

Q.

63

 

 

 

sjn22=Q)m22//i/i(l| + —

 

—//ли ;

sft„33 = Qk„33/ ft„ " 4 -2£V

- / ll),'2';

 

 

 

 

 

 

 

 

 

 

ал2

 

 

 

 

 

 

 

 

 

ал2

 

 

 

Skn'2=

Qkn2'Ikh< "

+

Ofen51

+

Q /m 12

/ „*<3>

Qkn42

м

151;

 

 

 

- ^

 

- / и й

-

^ V

+

- /

 

 

 

 

 

 

 

 

 

 

 

a ft2

 

 

 

aft2

 

 

 

a fe4

 

 

 

 

Sft„13 = Qkn3'Ikktl) + ^ J - h i , ™

+ - ^ V ~

W 3) + - ^ V "

/ftft'5';

 

 

 

 

 

 

 

 

 

 

 

a ft2

 

 

aft2

 

 

a/t4

 

 

 

W

3 =

(Q

im

23+

 

Q im 32)/ftft|l4

 

^

- ( Q

f t n

 

L

 

 

 

 

53+ Q ftnX 'kX62)/fti" \ d x .|2); /ftft|5' = J

 

 

 

 

 

 

 

 

 

 

 

n u*

 

 

 

 

 

 

 

 

 

 

Коэффициенты

Qftnap

 

(B= 1,2,3; a = l , . . . , 6 )

будут

определены в п.

2.

 

Используя

 

(9),

 

 

 

 

 

 

 

 

д(П + /(-Л )

Л

д{П + К - А )

Л

 

 

(10), из условии ------—-------- =0;

-------^

--------=0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dUfm

 

 

 

 

дУкп

 

д(П + К А )

 

 

 

 

 

 

 

 

 

 

 

 

,

п получаем систему трех одно-

■-----— -----— = 0

для фиксированных k,

OWhn

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

J

 

d2ar(1>

1

dar(')

/i2 + 1

 

 

2/i

 

Gir

 

 

dr2

 

г

dr

r2

 

ar(l)+ a2ar<2) H— — аф(1>+ ——ar+ = 0;

 

 

 

 

 

 

 

 

 

 

d2ar{2)

1

dar(2>

/i2 + 1

 

 

2/i

 

0)"

 

 

dr2

 

r

dr

r2

 

ar,2) + a 2ar(14 — г—аф(2) H--- — ar<2>= 0;

 

 

 

 

 

 

 

с2‘

(17)

,

^flcpO

 

1

da,p(l>

n2 + 1

 

 

 

 

 

■+

а Ф(1) +

а 2а

ф(2Ч —

 

СО" аф<‘) = 0;

 

dr2

 

r

dr

~72

 

 

С2г a

r(1) +

 

d2aф(2)

■+

1

da<p<2)

/12 +

1

аФ(2) + а2аф^Ч—

 

СО"

 

 

dr2

 

dr

 

 

аг(2)H---- аФ(2) = 0;

 

 

 

 

 

 

 

 

 

C2"

 

 

 

d2^ 1)

1

daz(1)

-----— az(1) + a2az(2) +

CO"

■flz(,»= 0;

 

 

 

dr2

 

r

dr

 

r2

 

 

c2z

 

(18)

 

 

d2azW

1

daz(2)

 

 

 

 

 

 

 

 

 

TV*

 

 

гл2

 

 

 

 

 

dr2

■+- r

dr

----- — az(2)+ a2az(1>+- с2г

az(2) = 0.

 

Решение (16) —(18) с учетом (14) имеет вид:

 

 

 

 

ф (1) = _ 2 ~ [ Л

1/ 7 1 ( | + r ) A- A J n { ^ ) +

Л 7 У71 ( P i r ) + A \o Y n ( p 2r )

] ;

 

ar(1) —— [-^2/ n+i {y\r) + Л3/ Т1_1(Y+) +-^5/ 71-1 (y2r) +71б/п+1 (угО +

 

+ T18^n+i (Y ir) + A 9 Y n- i

(Yir) + 71цУп_1 (Y 2^) + Л12У74-1 (Y2O] i

 

 

j

 

 

 

 

 

 

 

 

 

 

(19)

 

= — [ —Л2/ n+i (Y+) +71з/71—1(Yir) +715/ 71-1 (Y2O —AQJn+i (Y2^) —

 

—Tig У7i+i (YIO + 7l9У7i—i(YIO + + 11У71-1 (угО —7112 У7i+i (y2r) ];

 

 

0z(1) = - -- [Y13—Л2)/71 (Y+) + Y2(7l5 +б) /тг 2Г) +

 

 

 

 

 

4a

 

 

 

 

 

 

 

 

 

 

 

+ Yi (7I9—7l8) M YIO +Y2(Ли

7112) Утг (Y2^") L

 

где

Pi=

11/

CO2

 

 

P2У — 2 - a 2;

YI =

“I / CD2

 

 

l / 0)2

|^ —2 + a2;

l' —5-+ a ; Y2=

К 7 Г

 

 

Г

C 12

 

 

C i2

 

r C22

 

 

r C22

Выражения для q72), ar<2), аф(2\ az(2) получаются из (19) заменой знака на противоположный перед членами с коэффициентами 714, Т15, Л6, Лю, Ли, Лю.

Окончательные выражения для перемещений и напряжений в запол­ нителе имеют вид:

Ч г = X

J

f

( г ) + X

hl ~ f k n [A) (r) ]

COS /К р егшАп( ;

" о

"

L

a *

a *

 

 

 

 

Xh(z)f/m(2) (r) + X ^

--fhn{5) {r) ]

sinmpe*0*»*; (20)

 

71

/I

a,t

 

Kr= X I [x»(z)f»«<3|(') + X *{— h»,6>(r) ] cosncpe*”»»';

n h