Дифференциальные уравнения связанной задачи

advertisement
УДК 539.3
А.Н. Андреев
Дифференциальные уравнения связанной задачи термоупругого
деформирования слоистой композитной оболочки
Аннотация
А.Н. Андреев
Дифференциальные уравнения связанной задачи термоупругого
деформирования слоистой композитной оболочки
Исследуется связанная задача термоупругого деформирования слоистой
оболочки. Дифференциальные уравнения задачи, позволяющие обеспечить
сопряжение полей деформаций и температур, учесть явление поперечных
сдвиговых деформаций и нелинейный характер распределения температуры по
толщине, получены из обобщенного принципа виртуальных перемещений.
A.N. Andreev
Differential Equations of the Coupled Thermoelastic of a Layered Composite Shell
Abstract
The Coupled Thermoelastic Problem for a Layered Composite Shell Is Investigated.
Differential Equations of this Problem Considering Conjugation of Fields of
Deformations and Temperatures, and Taking Into Account the Phenomenon of
Transverse Shear Deformations and Nonlinear Character of Temperature Distribution
Over the Thickness, are Obtained From the Generalized Principle of Virtual
Displacements.
Ключевые слова: дифференциальные уравнения, композитная оболочка,
термоупругое деформирование, связанная задача.
Key words: Differential Equations, Composite Shell, Thermoelastic Deformation,
Coupled Problem.
Рассматривается слоистая композитная оболочка постоянной толщины h,
собранная из m слоев также постоянной толщины, каждый из которых
1
2
армирован семейством однонаправленных волокон. Пусть x , x - внутренние
1
координаты отсчетной поверхности, в качестве которой примем «нижнюю»
лицевую поверхность тела оболочки, z - нормальная координата. Уравнения
поверхностей раздела j -го и ( j  1) -го слоев запишутся в виде
z  hj ,
j 1, 2, ..., m 1,
(1)
где h j  const, 0  h0  h1    hm  h. На поверхностях (1) поперечные
компоненты 
3
,  33 тензора напряжений и компоненты v , v3 вектора
перемещений должны удовлетворять условиям непрерывности [1], а
приращение  температуры – условиям идеального теплового контакта [2, 3].
Уравнения Дюамеля – Неймана анизотропного упругого тела запишем в виде
   A    c ,  3  p   3 ,  3  q 3  .

Здесь A
(2)
, p - компоненты тензоров эффективных тангенциальных и
поперечных сдвиговых жесткостей соответственно, c

- тензор
ij
температурных жесткостей. Обозначив через K эффективный тензор
коэффициентов теплопроводности композитного материала, запишем закон
i
Фурье для армированной среды ( q - компоненты вектора теплового потока):
q i   K ij  j .
(3)
Компоненты указанных тензоров определяются расчетным путем [1, 3] и
являются функциями механических и температурных характеристик элементов
композита и структурных параметров армирования.
Неклассические дифференциальные уравнения задачи термоупругого
деформирования многослойной анизотропной оболочки строим на основе
следующего допущения [1] о законе распределения поперечных компонент
тензора деформаций по толщине оболочки:
0( k )
 03  zh 1 03h   03   f ( z)  ,  330(k )  0.
 0(3k )  q
(4)
Здесь k  1, 2, ..., m - порядковый номер слоя,  0 ,  h - значения поперечных
3
3
сдвиговых напряжений на верхней и нижней лицевой поверхности
2
соответственно,
  ( x1 , x 2 , t ) - независимые кинематические характеристики,
учитывающие наличие поперечных сдвиговых деформаций, f (z ) функциональный параметр, задающий закон распределения этих деформаций
по толщине слоистой оболочки. Распределение компонент вектора
перемещений по толщине многослойного пакета, соответствующее закону (4),
определяется из соотношений деформации – перемещения [1] и имеет вид
v30 ( k )  w( x1 , x 2 , t ), v0 ( k )  0 ( k ) ( x1 , x 2 , z , t ) 
0( k )
 u ( x1 , x 2 , t )  z ( x1 , x 2 , t )  
( x1 , x 2 , z , t )  ( x1 , x 2 , t ).
(5)
Легко видеть, что представления (4), (5) удовлетворяют условиям межслоевого
контакта и условиям нагружения на поверхностях z  0, z  h оболочки.
Распределение тангенциальных деформаций по толщине пакета слоев,
отвечающее кинематическому закону (5), записывается в виде
(k )
2 
 (   zb )   ( k )  (   zb )  ( k )  (   zb )  u 
(6)
 (   zb ) u  z (   zb )   z (   zb ) 
 (   zb )  ( k )     (   zb )  ( k )    (   zb )  ( k )   
 (   zb ) ( k )    2(b  zb b ) w.
Предполагая, что на лицевых поверхностях z  0, z  h известны значения
температур:
при z  0   0 , при z  h   h ,
(7)
примем следующий закон распределения приращения температуры по толщине
h  0
s 1 K 33
 K133

1




s  g ( z ) ( x , x , t )  m 33
z

h

h

h

s 1
i
i 1   0 ,
 K 33
33
K1
K
i

1

i
 33 hi  hi1   s
i 1 K i
1
2
hs1  z  hs .
(8)
В (8) s  1, 2, ..., m - порядковый номер слоя, ( x , x , t )  независимая
1
2
характеристика, учитывающая отклонение поля температур от линейного
3
закона, g (z ) - заданная непрерывно дифференцируемая функция такая, что
g (0)  g (h)  0,
33


K 33
j g ( h j  0)  K j 1 g ( h j  0)
(9)
( j  1, 2,..., m  1),
K s33 - коэффициенты тензора теплопроводности s -го слоя. Легко убедиться,
что распределение температуры (8), удовлетворяет условиям нагрева (7) на
лицевых поверхностях оболочки и условиям идеального теплового контакта на
поверхностях раздела ее слоев. Функцию g (z ) можно взять в виде
z
g ( z )   t  h1 t  h2     t  hm 1 dt ,
(10)
0
если
h
g (h)   t  h1 t  h2  t  hm1 dt  0
(11)
0
(условие (11) выполнено, например, в случае двухслойной оболочки, геометрия
пакета слоев которой симметрична по толщине), и в виде
z
g ( z )   t  a t  h1 t  h2     t  hm1 dt,
(12)
0
если условие (11) нарушено. Параметр a в (12) однозначно определяется из
условия g (h)  0. Другой способ задания функции g (z ) предложен в [4].
Дифференциальные уравнения связанной задачи термоупругого
деформирования слоистой оболочки получим из обобщенного принципа
виртуальных перемещений [2]:
 W  P  D   L  H n dA,
(13)
A
где W 
c
1 ijkl
A


dV
,
P

 2 dV , D   ij H jH i dV , L   p ivi dA 


kl ij
2B
2T0 B
B
A
  F i  vi vi dV - потенциальная энергия деформации, температурный
B
потенциал, вариация функции диссипации D , виртуальная работа
поверхностных и массовых сил соответственно, q  To H   K  j i
4
i
ij
компоненты вектора теплового потока, c ,  - теплоемкость и плотность
композитного материала соответственно, B - объем, занятый телом оболочки,
A - поверхность, ограничивающая тело оболочки.
Используя эти соотношения, приводим уравнение (13) к виду



3
ij

       3  T c   i ( K  j )  T0c   dV 
B

0
g ( z)
 wdV    p v  p 3wdA 
  F   v v  F 3  w
B
(14)
A
1
  H n g ( z )dA   ni K ij  jg ( z )dA.
T0 A
A
Подставляя в вариационное уравнение (14) выражения (6) – (10), выполняя
интегрирование по поперечной координате z и приравнивая к нулю множители
при независимых вариациях обобщенных перемещений и множитель при
вариации поля температур отдельно в поверхностном и контурном интегралах,
приходим к системе дифференциальных уравнений связанной задачи
термоупругого деформирования слоистой оболочки
X   b Y   T  b  M   3   03  hb 3 ,
I    Y  b T     M   33   033   03    3     h  3  ,
Z    S  Q    3  (m) ( x1 , x 2 , h),

(15)

h
 
 
3
3
33
ˆ
ˆ
ˆ
c    K     K     K   K  T0  c  g ( z )dz W .
0
Здесь
T
m hk
Q   
k1 hk 1
X



( k)


,M ,S

.
.
   
m hk
k 1 hk 1


(k )
( k )
 zb  (
k ) 1, z ,    g a dz,

( k )
3
 3
( k )
 zb (
k )     ( k )  zb  ( k ) q f ( z )

g a dz ,
( k ) 
,Y , Z     ( k ) ( k )  u   z    
 1, z, 
( k )  g a dz,

m hk
k 1 hk 1
5
m hk
I  w   ( k )
k 1 hk 1
h
h
0
0
g a dz, Kˆ    g 2 ( z ) K  dz, Kˆ  3  Kˆ 3   g ( z ) g ( z ) K  3dz,
h
h
0
0
Kˆ 33   g ( z ) g ( z ) K 33dz, cˆ   c g 2 ( z )dz.
Замкнутая система (15) уравнений термоупругого деформирования слоистой
композитной оболочки, позволяет учесть сопряжение полей деформаций и
температур в теле оболочки, явление поперечных сдвигов в ее слоях,
нелинейный закон распределения температуры по поперечной координате.
Сопряжение полей деформаций и температур осуществляется подчеркнутым
членом в системе (15) и исчезает только в случае стационарного теплового
потока. Из вариационного уравнения (14) следуют также и краевые условия.
Кинематические и естественные граничные условия установлены в [1] и здесь
сохраняют свой вид. Тепловые краевые условия отражают взаимодействие
контура оболочки с окружающей средой. Этими условиями в точках контура
задается либо функция  , либо поток тепла через контур.
К достоинствам предлагаемой модели следует отнести независимость порядка
соответствующей системы дифференциальных уравнений и структуры этой
системы от количества слоев оболочки и строения пакета слоев.
Библиографический спсок
1. Андреев А.Н., Немировский Ю.В. Многослойные анизотропные оболочки
и пластины. Изгиб, устойчивость, колебания. Новосибирск: Наука, 2001.
288 с.
2. Новацкий В. Динамические задачи термоупругости. М.: Мир, 1970. 256 с.
3. Немировский Ю.В., Янковский А.П. Рациональное проектирование
армированных конструкций. Новосибирск: Наука, 2002. 487 с.
4. Немировский Ю.В., Бабин А.И. Теплопроводность многослойных
армированных оболочек// Нелинейные задачи расчета тонкостенных
конструкций. Матер. Всес. конф. Саратов, 7 – 9 июня 1988 г. - Саратов,
1989 г. - С. 126 – 130.
6
Андреев Александр Николаевич, доктор физико-математических наук,
профессор, заведующий кафедрой алгебры и геометрии Кемеровского
государственного университета; 650043, г. Кемерово, ул. Красная 6.
e-mail: algebra@kemsu.ru
7
Download