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

УД/< 539.3.001:678

И. Н. Преображенский, Ю. М. Коляно, О. И. Борисенко

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

Одной из основных проблем механики деформируемого твердого тела является термомеханика тел неоднородной структуры [1]. Основные урав­ нения классической термоупругости тел двухмерной и одномерной струк­ туры выведены в [2, 3] с использованием аппарата обобщенных функций, и на их основе получены решения квазистатических и динамических за­ дач несвязанной термоупругости, единые для всей области их определе­ ния. При этом1предполагалось, что скорость распространения тепла — величина бесконечно большая. Однако учет конечной скорости распро­ странения тепла оказывает существенное влияние на динамические про­ цессы в твердых телах [4]. Приведем основные уравнения неклассической (обобщенной) термоупругости тел трехмерной, двухмерной и одномер­ ной структуры, используя при этом развитую в работах [5, 6] алгебру асимметричных обобщенных функций.

Пусть в изотропном теле, занимающем область Q, ограниченную по­ верхностью 5, имеется М3 слоев изотропных параллелепипедных включе­ ний по М2 рядов в каждом слое, состоящих из М\ включений в каждом ряду. Грани включений параллельны координатным плоскостям декарто­ вой прямоугольной системы координат **, начало которой разместим в геометрическом центре системы включений. Размер включений по оси Xk обозначим через 2 ЪкУа расстояние между центрами включений по оси хь.

через 2dk. Тогда akr= (2r- 1 - M k)dk- b ky $kr= ( 2 r - \ - M k)dk+ bhy где

г — номер слоя включений по оси Хи\ анг, |3ьг — координаты начала и конца этого слоя соответственно.

Система дифференциальных уравнений взаимосвязанной термоупру­ гости неоднородного изотропного тела имеет вид [7]

[p(Xi, Х2, *з)

+

+ Я(Хь X2J X3) U/i'kSij —

 

 

“ P { X u

*2, * 3) *6ij] , j = p ( * b * 2, X 3 )

 

(1 )

[h (xu x2, X3) t,i] = l[cv(xu x2y X3) i +p {xu x2y xz) tQeu-

wt],

(2)

где Я, |x — коэффициенты, учитывающие упругие свойства

материала;

Р — коэффициенты, учитывающие механические и теплофизические свойства материала; щ — компоненты вектора перемещений; р — плот­ ность; cv — объемная теплоемкость при постоянной деформаций; Я* — коэффициент теплопроводности; /0 — температура тела в ненапряжен­ ном состоянии; t — приращение температуры тела; wt — плотность внут­

ренних источников тепла; еZj = (Uitj + Ujyi) — компоненты тензора дефор­

маций; т — время; / = 1+ T*J^; T* = const — время релаксации теплового

потока; i yj y k = 1,3; точки над величинами обозначают производные по времени т; 6„ — символ Кронекера.

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

температуры принимаем в виде [4]

 

tiiXt (Р) [t (Р, т) ] ,1 + las (Р) [*(Р, т) - to*(Р, т) ] = 0,

P eS ;

t(P, 0) = t0(P)-, i(P, 0) =f0(P); P ^Q ,

(3)

где tii — компоненты вектора внешней нормали п к поверхности тела; а8 — коэффициент теплопередачи с поверхности 5 тела; tc8 — темпера­ тура среды, омывающей эту поверхность.

Термомеханические характеристики тела как единого целого пред­

ставим в виде

з

p(xi,X2,xa) = p i + -(po-Pi) Ц м(Хк),

(4)

А-1

 

где ро, р1 — термомеханические характеристики материалов включения и основного материала;

мк

 

 

N(Xh)= ^

[S_(*ft —<хаг) —S+(XA—£аг)] —характеристическиефункции;

Г=1

 

 

 

1,

6> 0,

5±(6) =

0,5+ 0,5,

| = 0, — асимметричные единичные функции.

 

о,

К О

Подставляя в (1), (2) выражения для термомеханических характе­ ристик тела в виде (4), используя при этом свойства единичных асиммет­ ричных функций [5, 6], а также соотношение

Р ( Х 1,Х 2,Х з) Pi \ ро Р\ >

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

“«+[-дг+1+(-^--7г) П

%.«-

- [ - Ы £ - £ ) П * ( * М £ + ( £ - £ ) х

Pi

 

ро Pi

3

 

 

X П N (xk) 1

(Xi)N(xm)N (xi)t— ( — — 1) X

X TiN'(x})N(x,)N(xn) (Uij+uj,t) - ■— %l- N ' ( Xi)N(xm)N (x,)ujl};

(5)

i-1

 

3

 

[ т г + ( - ^ - 7 г ) П * < * > ] i+

x

где m=<o(0 ; / = со2(0 ; *=*>(/); n=co2(/);

ш = ( * * ®) ;

a = — — коэффициент температуропроводности; c9= ]/-^

— скорость

r T

 

распространения тепла; p= at(3A,-t-2p); at — температурный коэффи­ циент линейного расширения; X, р, — постоянные Ламе.

Если поверхность 5 и поверхности включений не имеют общих точек, то граничное условие (3) в этом случае имеет вид

ЪХи [t (Р, т) ] .i+ a.i/[f (Р, т) - tc‘ (Р, т) ] = 0; Р е 6\

Если часть граничной поверхности S ic S

есть плоскость, совпадаю­

щая с боковыми гранями включений, то (3) запишется в виде

2

 

п\ [ Яп+ (hto—Xu) J J N (Xk) j

X)] fi+

2

+ [a si+ (aso -asi) Д ] У ( ^ ) ] l[t(P, r) - t c8(P, T)] =0;

P e S b 6 = 1, 2.

Если влияние деформации на изменение температурного поля прене­ брежимо мало, то уравнения (5), (6) становятся несвязанными и урав­ нение теплопроводности (6) принимает вид

[w+(

~v)ПN{Xh)1/+[ с<г1'2+

 

з

 

 

+ (Сдо~2—Сд1~г) Ц

N (Xh) ]

i —

- ( у 2- - 1 ) £ * ' ( * ) X

4

<_1

 

3

 

 

 

X Ы(Хт)М(Х1) и - [ - i - + ( J -

- - i - )

П**{хк) ] lwt. (7)

Рассмотрим некоторые частные случаи армирования. Положив в урав­ нениях (5), (6) и (5), (7) Мг=\, приходим к уравнениям обобщенной взаимосвязанной и несвязанной задачи термоупругости для тела, арми­ рованного параллелепипедными включениями, расположенными в слое, ограниченном плоскостями х3= ± 63. Из тех же уравнений при M i=\, Мз= 1 вытекают уравнения для тела, армированного рядом параллелепипедных включений.

Положив в уравнениях (5), (6) и (5), (7) N(x3)== 1, приходим к уравнениям обобщенной взаимосвязанной и несвязанной задачи термо­ упругости для тела со сквозными включениями, поперечное сечение кото­ рых плоскостью х3 = 0 — прямоугольники размера 2 Ь\Х262-

Уравнения для слоистого тела с плоскопараллельными равностоя­ щими границами раздела получаются из (5), (6) и (5), (7) в случае

N(xг) = 1, N(x3) ==1. Например, из (5), (6) получим

 

 

 

|*о

м

N(xi)]Uj,ji

[SL +\[ К . _1L X

М-1

!*1 '

 

хЦ-о

II

 

 

 

*

 

 

[ - -+ (—

1N(Xl)]1

 

 

 

LНЧ

х|Хо I*I'

 

 

' z J L l L j V ' f x J t -

- l )

N'(xI)(Ui,i + ul,i)+ - ? ? - ^ - N '(X i)uj;;

Ц1

'

Hi

 

 

 

 

l*i

 

 

 

t,u= Г — +

( - ------- - )

/ + [ ^ Г 2+ (с д0- 2- с 9г 2) X

L CL\

' CLQ

CL\

J

 

 

 

[ - ^ - +

 

^ (* i)]

M u -

-

 

 

[

i i r + ( i r - т г )

A,(J" ) ] * " •

Рассмотрим тело, состоящее из М изотропных слоев, краевые поверх­ ности 5 1и SM которого совпадают с плоскостями Х\ = 0, х х= gM. Термоме­ ханические характеристики тела как единого целого представим в виде

М -1

 

р (*0 = P l + J j (Pr+l—Pr)S-(Xi—lr),

(8 )

где рг — термомеханические характеристики г-го слоя, — координата границы раздела r-го и г + 1-го слоев.

Подставляя в (1), (2) представление (8) термомеханических харак­ теристик тела и используя при этом свойства асимметричных единичных функций [5, 6], а также соотношение

т т

(я ж “ Яг)5±(л:—£r) J

(br+\ — Ьг)8±{х—£>г) =

Т=1

г-1

 

m

 

 

= ^ Л ±'6±(* -£ г),

 

Г=1

 

 

где A+r=ar+\ {br+\ — br) ; А-?=аг(Ьг+\ — Ьг) ; 6±(g)

, получим сис­

тему дифференциальных уравнений с коэффициентами импульсного типа обобщенной взаимосвязанной термоупругости изотропных многослойных тел с плоскопараллельными границами раздела:

М -1

 

иг,Н+ [ — - + 1+

(

— —--------- S-(*i —Er)]

 

 

 

 

L

Щ

~

\

Pr+l

\ir

1

 

 

J

 

 

 

 

 

M -\

 

 

 

 

 

 

 

 

 

 

- [ A + X

( i S L — t

)

s

.

(

[ - *- +

 

M - l

L

Hi

"

x F^-H

 

Hr7

 

M -l

J

 

1 F*1

 

 

 

 

 

 

 

 

 

 

 

 

 

'

Jirf|

 

Jlr >

 

 

J

 

l

L

|*T

 

- ^

7

r

Luui :

 

 

L

 

 

 

 

 

(9>

 

 

^ T + §

( а

7

~ ^ г ) 5- ,*|~ Ег)]

i+

[ c ,r ’ +

 

 

M - l

 

 

 

 

 

 

 

 

 

M - l

 

 

+

^

(C q r+ r2— Cqr~2) S _ ( X i — gr) 1

t +

F

-----------------------------------

 

 

r-1

 

 

 

 

 

ЛГ-1

J

 

L

"

 

V

 

 

 

 

 

 

 

 

£

( - £ f - l )

 

 

6- t a - i , ) -

 

 

 

 

 

 

M - l

 

 

 

 

 

 

 

 

 

 

-

[ Х7Г +

s (

i

b -

- i

r )

s -<*‘- S ' > ] t o ‘-

(10)

Обобщенное граничное условие теплообмена третьего рода (3) в данном случае принимает вид

h d t(P , т)] ,1 - a sll[t(P,'c)-tca(P ,r)]= 0\

P<=S,;

^•tM\t(P,t)\,i + o>ml\t{P, т) —/с8(Р,т)] = 0;

P ^ S M-

Температурные напряжения в теле определим из соотношений Дюа- меля—Неймана [7]

<т,,=2 ц(хи х2, х3)e{j+ [Я,(^ь х2,х з)е - р (х и х 2,х з)ф ц ,

(11)

где e=Uk,h, а величины р, К, р выбираем в виде (4) или (8) в зависимости от геометрии неоднородностей тела.

Рассмотрим частный

случай,

когда v = const;

Е = Е (хi), at = at{xi).

В этом случае с использованием

(11)

уравнения

(9),

(10)

запишутся в

виде

 

 

 

 

 

 

 

 

 

 

М-1

 

 

 

 

 

 

 

1

 

 

 

2(l+v)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

“MJ+ l - 2v Uj,ji

1—2v

■[«fl+

r-1

(0C(r+l —air)

X

 

 

 

 

 

 

 

 

 

 

 

M-l

 

 

 

 

 

 

X S~ (JCI

£r) \ t ,r —2(1 H-v) [ -£L+

t i

(

 

Г

 

Er

)

 

 

 

 

 

 

 

 

 

Ei

r-1

X£r+1

 

 

 

 

 

 

 

 

 

 

 

M-l

 

Er+\

 

Er f

 

 

 

 

X S -(xi~ ■Ir) | « i= 2 (l+v) ^

 

 

6nX

 

 

 

 

Er2

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v ^ [ ( a j r + i — l ) £ r + i — { a t r — 1 )£ V ] . 1

1

 

 

с /

t 4 .

П о \

Х ---------(£ж

 

- & ) ( l - 2v)--------- <12>

 

 

 

 

 

 

 

 

 

 

 

М-1

 

 

 

 

 

 

 

 

 

 

 

 

 

t.ii— [ ~

 

НX!l

(~

------- Т- )

 

 

С*1

 

1г) ]

 

£+

 

 

 

L

ai

 

 

\ar+i

 

йг /

 

 

 

 

 

J

 

 

 

 

 

 

 

М-1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ [

c 9 l“"2+

^ j

( c gr+l“ 2 — ^qrr~ 2) S _ (л:1— £r) j

 

t +

 

 

 

 

 

 

r-1

 

M-l

 

 

 

 

 

 

 

 

 

 

 

 

1

 

Г

at\E\

/

a/r+i^r+i

 

atrEr

 

 

 

 

^

 

 

 

1—2v

 

L

Xti

 

^

'

Xir+i

 

 

 

O x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M-l

 

 

 

 

 

 

 

 

 

 

 

 

X S ~ (* i-g r) ]

w<«

£

(

Y

 

'

0

1

 

,

6

 

 

 

 

 

 

 

Xl=lr

 

 

 

 

 

 

 

M-l

r-1

 

Air

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

" U

r

+ 5

^

 

 

 

 

xir )

 

s ~{Xl - ? r ) ]

 

 

(13)

 

 

 

 

 

 

 

 

Уравнения связанной динамической термоупругости многослойных тел получается из (9), (10) или (12), (13) в предположении, что ско­ рость распространения тепла есть величина бесконечно большая (т* = 0). В этом слуЧЪЪ (10), (13) принимают вид

М-1

* , « « [ —

+ £

М

------ J - ) s

_ ( x 1- £ r) ] / +

[ ^

- +

L Oi

'

От+1

ar '

J

L

At!

M-l

 

 

 

 

M-l

 

 

X/.:

hjr+£

•“

 

 

 

t v +

%

) 5 -<jci- £' ) ] ' + т ^ г [ ' л г 1 _ +

+ £

( “ 7 g w ■-

 

s - f e - ы ] м „ -

г-1

\ Atr+ 1

W r

9'

М -1

М -1

Пренебрегая влиянием деформации на изменение температурного поля, а также инерционными членами в уравнениях (9), (12), приходим

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

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

1.Ильюшин А. А. Две проблемы механики деформируемого твердого тела (МДТТ). — В кн.: Тр. Всесоюз. конф. по уравнениям с части, производными, посвящ.

75-летию со дня рождения И. Г. Петровского. М., 1978, с. 127—130.

2.Коляно Ю. М. Термоупругость тел неоднородной структуры. — В кн.: Всесоюз. семинар по теории упругости неоднород. тела. Ереван, 1981, с. 22—30.

3.Коляно Ю. М. Применение обобщенных функций в термомеханике кусочно-од­ нородных тел. — В кн.: Математические методы и физико-механические поля. 1978,

вып. 7, с. 7—11 (Киев).

4.Подстригая Я. С., Коляно Ю. М. Обобщенная термомеханика. Киев, 1976. 312 с.

5.Коляно Ю. М., Попович В. С. Об одном эффективном методе решения задач термоупругости для кусочно-однородных тел, нагреваемых внешней средой. — Физ.- хим. механика материалов, 1976, № 2, с. 108—112.

6.Коляно Ю. М., Процюк Б. В. Термоупругость неоднородных и кусочно-однород­ ных пластин, обладающих цилиндрической анизотропией. — В кн.: Обобщенные функ­ ции в термоупругости. Киев, 1980, с. 19—34.

7.Коляно Ю. М., Штер 3. И. Термоупругость неоднородных сред. — Инж.-физ. журн., 1980, т. 38, № 6, с. 1111—1114.

Институт прикладных проблем механики и математики Поступило в редакцию 10.02.82 АН Украинской ССР, Львов