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

книги / Прямые и обратные задачи механики упругих композитных пластин и оболочек вращения

..pdf
Скачиваний:
9
Добавлен:
12.11.2023
Размер:
17.75 Mб
Скачать

10.2. Влияние выбора структурных моделей КМ

223

уменьшаются в 5 раз, в продольной арматуре в 3 раза, в окружной арматуре и прогибы в 2,5 раза.

Влияние механических характеристик. Исследуем влияние выбора модели КМ на НДС жестко защемленного комбини­ рованного сосуда давления для различных материалов армату­ ры и связующего. На рис. 10.10 представлены максимальные

безразмерные прогибы W

и интенсивности напряжений

в эле­

ментах

КМ в зависимости от соотношения между механиче­

скими

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

арматуры и связующего О =

Е {а/ Е хс.

Результаты получены по уточненной теории [9] при структуре (0,90,60,-60). Сплошным линиям соответствуют h\ = /13 = 0,1/г, штриховым — h\ = /13 = 0,4h.

Приведенные данные показывают незначительное влияние выбора модели КМ на величины максимальных интенсивностей напряжений и прогибов при hi = /13 = 0,4h и П > 30. Однако при П = 5 величины, полученные по МОВ и МДВ, различаются для связующего на 25%, для продольной арматуры на 35%. При hi = /13 = 0,lh различие более существенно и составляет при П = 30 для связующего 40%, для про­ дольной арматуры 20%, для спиральной арматуры 60%, для прогибов

224 Гл. 10. Сопряженные сосуды давления

30%. При увеличении отличия модулей Юнга арматуры и связующего влияние выбора модели КМ для связующего и прогибов уменьшается, а для арматуры увеличивается.

Уменьшение толщины слоя с окружной арматурой при О = 100 приводит к уменьшению интенсивностей напряжений в связующем в пять раз, в продольной арматуре в два раза, в спиральной арматуре в четыре раза, прогибов — в два раза.

Расчет нагрузок начального разрушения. Рассмотрим, как вли­ яет выбор моделей КМ на уровень нагрузок начального разрушения сопряженного сосуда давления.

На рис. 10.11 приведены зависимости нагрузки начального разру­ шения от угла укладки спирального семейства арматуры для углепла­ стикового сосуда давления, левый край которого подкреплен жесткой крышкой, а правый жестко защемлен. Геометрические и структурные параметры соответствуют рис. 10.9. Результаты получены по уточ­ ненной теории [9]. Кривым 2 соответствуют результаты, полученные по МОВ, кривым 3 — по уточненной модели КМ с одномерными волокнами и кривым 4 — по МДВ.

66

44

22

0

0

30

60

Ф

 

Рис. 10.11

 

 

Из рис. 10.11 видно, что при h\ = hz = 0,1 h для нагрузок начального разрушения существует максимум. При ф = 45° оболочка может выдер­ живать нагрузку в 8 раз большую, чем при ф > 70°. Разница величин, полученных по МОВ и МДВ, не превышает 15%. При h\ = hz = 0,4h зависимость нагрузки начального разрушения от угла спирального армирования становится монотонной. При ф = 10° нагрузка начального разрушения в три раза больше, чем для h\ = /13 = 0,1 h.

Влияние порядка расположения армированных слоев. Иссле­ дуем НДС комбинированно конструкции в зависимости от порядка расположения армированных слоев.

На рис. 10.12 приведены зависимости максимальных безразмерных прогибов, нагрузок начального разрушения, приведенных интенсивно-

10.2. Влияние выбора структурных моделей КМ

225

стей напряжений в связующем и продольной арматуре от угла спираль­ ного армирования жестко защемленной углепластиковой оболочки под действием постоянного внутреннего давления при hi = h$ = h/Ъ. Ре­ зультаты получены по уточненной теории [9] при использовании МДВ. Кривым 1 соответствует порядок расположения армированных слоев

(0,90, ф, —ф), кривым 2 — (0, ф, —ф, 90), кривым 3 — (90,

0).

0,6

ч

0,4

' 1

/ 2 / /

0.2

0 J- Л- L L

80

60

•10

20

На рис. 10.12 видно, что при ф « 40° расположение окружной арматуры во внутреннем слое (кривые 3) позволяет за счет пони­ жения уровня интенсивностей напряжений в продольной арматуре увеличить нагрузку начального разрушения в 1,7 раза по сравнению с другими структурами. Расположение слоя с окружной арматурой во внешнем слое (кривые 2 ) по сравнению с расположением в среднем слое (кривые У) несущественно влияет на НДС оболочки. Порядок рас­ положения армированных слоев на величины максимальных прогибов практически не влияет.

8 С. К. Голушко, Ю. В. Немировский

226 Гл. 10. Сопряженные сосуды давления

10.3. Анализ достоверности численных решений

В табл. 10.2-10.4 приведены максимальные относительные раз­ ности в равномерной метрике для величин, характеризующих НДС углепластикового жестко защемленного сопряженного сосуда давления и рассчитанных с помощью методов сплайн-коллокации и дискретной ортогонализации при параметрах: h \ = h z = 0,4/г, (0,10,-10,90).

Параметры

T O L

J

ю - 4

300

10"6

600

10"8

1200

T O L

J

10"4

300

10"6

600

10"8

1000

 

 

 

 

Т а б л и ц а 10.2

Относительная разность 8 по компонентам

W '

W

 

и

M l

Классическая линейная тео зия

 

1,21 • 10“ 2

6,95 • 10“ 3

6,9310“ 3

8,59- 10-3

8,7810"4

4,3910"4

5,26 • 10-4

6,02 • 10-4

3,9910"8

6,5410"8

4,61 • IQ"8

8,64 • 10-9

Классическая нелинейная теория

 

1,11 • 10"2

6,55 • 10_3

о

1

8,78

• 10-3

о

8,8310"4

4,57 - 10-4

сл со

1

6,32

• 10-4

о

4,02 • 10-8

6,86 • 10-8

4,87 - 10-8

8,75

• 10-9

Параметры

T O L J

10“4

300

10_6

600

10“8

1000

T O L

J

10"4

300

10“6

600

10"8

1000

Т а б л и ц а 10.3

Относительная разность 8 по компонентам

W '

 

W

 

и

 

Mi

 

Линейная теория Тимошенко

 

2,41 • 10“ 2

Vi сл

1

7,43

• 10“ 3

9,62

■10“ 3

о

9,76

• 10-4

5,21 • 10"4

6,87

10-4

7,54

■10"4

ОО СЛ

1

9,77

• 10“ 9

5,56

10-8

7,82

• 10“ 9

О ОО

 

Нелинейная теория Тимошенко

 

2,65

• 10~2

7,41 • 10"3

7,62

• 10"3

9,87

• 10"3

9,89

• 10“ 4

5,56

• 10-4

6,97

• 10“ 4

7,87

• 10“ 4

8,73

• 10-8

9,45

• 10“ 9

6,01 • 10“ 8

7,32

• 10"9

Результаты

получены при использовании МДВ. Результаты

в табл. 10.2

получены

по

линейной

и

нелинейной

классиче­

ской теориям,

в табл.10.3 по линейной и нелинейной

теориям

Тимошенко, в

табл. 10.4

по

линейной

и

нелинейной

теориям

Андреева-Немировского [9].

 

 

10.3. Анализ достоверности численных решений

 

227

 

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а 10.4

Параметры

Относительная разность 5 по компонентам

T O L

J

W '

 

 

W

 

П

 

Mi

Линейная теория Анд эеева-Немировского [9]

 

 

 

 

1 0

 

" 4

250

1,64-

1 0

" 2

8 , 1 1

1 0 " 3

4,1710" 2

1,28

1 0 " 2

1 0

6

500

9,74-

10

4

5,05

• 10“ 4

2,24

10~ 3

6,27

10- 4

1 0

8

1 0 0 0

8,31 •

10-

6

5,3710~ 8

2,32

10_ 6

3,81

■10“ 8

T O L

J

Нелинейная теория Андреева-Немировского [9]

1 0 ~ 4

250

1,82-

1 0

" 2

8,65- 1 0 " 3

4,51 •

10- 2

1,33-

10“ 2

1 0 ~ 6

500

9,87-

10"

4

5,44

■10" 4

2,56

10" 3

6,38

10- 4

1 0

- 8

1 0 0 0

3,43 •

10_ 6

6,38

■10- 8

1,56

-

10“ 6

2,14

10“ 8

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

8:

Г л а в а 11

АНАЛИЗ НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ МНОГОСЛОЙНЫХ ПОЛИАРМИРОВАННЫХ ЭКСЦЕНТРИЧЕСКИХ КОЛЕЦ

11.1. Получение разрешающей системы уравнений

Рассматривается тонкое многослойное армированное эксцентриче­

ское кольцо (рис. 11.1)

внешнего

радиуса г\

с отверстием радиуса го

и эксцентриситетом d. В кольце К

= 2d + 1 слоев, которые ограничены

поверхностями раздела

—h/ 2 = h - d -1 < h-d

< ■■■< hd = h / 2 .

Пластина нагружена равномерно распределенным поперечным внешним давлением интенсивности Р. Используются два вида краевых условий: жестко защемленный внутренний край, на внешнем контуре приложена равномерная растягивающая нагрузка То; оба контура жестко защемлены.

Эксцентрическое кольцо рассматривается в ортогональной биполяр­ ной системе координат, которая вводится так, чтобы две координатные линии системы ап = а\, ап = а." совпадали с контурными линиями пластины. Связь с декартовыми координатами имеет вид

 

mshai

,

m sin ап

, 1Ч

х = - ----------

у = ----------- -— .

(11.1)

 

cnai + coscc2

 

chai + совссг

 

Координатные

линии ап

Q!Y = const > О представляют

собой

 

 

Г

с центром на горизонтальной оси

окружности радиусом г = т / sha^

в точке (m cotlm p

0); линии 0:2 = а® = const — окружности с центром

на вертикальной оси, взаимно ортогональные с окружностями ап = оР{.

11.1. Получение разрешающей системы уравнений

229

Для определения параметра т получаем систему трех уравнений

m^thc*! —ctha") = d,

m /sh ai = n ,

m / she*'/ = ro,

решив которую находим

Выражения для параметров Ламе при биполярной системе коорди­

нат имеют вид

А\ = А 2 =

chai + cos «2

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

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

При решении задачи определения НДС композитного эксцентри­ ческого кольца использовались теория Кирхгофа-Лява и уточненная теории Андреева-Немировского.

Для пластины симметричного относительно срединной поверхности строения при выборе в качестве отсчетной этой поверхности задача определения НДС распадется на две независимые задачи растяжения и изгиба. Для теорий Кирхгофа-Лява и Андреева-Немировского урав­

нения, описывающие задачу растяжения, совпадут.

 

При выводе разрешающие

систем

уравнений

используются

безразмерные параметры

т,

А,

А и А 2,

Г03, Г23, 9 ^ ,

h]3mtp, /<£\ f lpm

( i , j , l , m = 1,2):

 

 

 

 

 

п

 

(chai + cosc*2) ’

 

230 Гл. 11. Анализ НДС многослойных эксцентрических колец d

£ ( ! ' ( z ) f q lmdz = f ^ h s/ E i ' \

fc= —dL

 

 

 

Пк- 1

 

 

 

f ( z ) q lmdz = fj,mh3 /Ei'K

В следующих соотношениях (j

= 0 ,.... К,

к = l , .... К, п = I , ... ,7)

P i = h j / h ,

= £ « / £ < » ,

=

b f = a f / E ^ \

ё ^ = д Р .

4 k ) = g \* \

Разрешающая система уравнений задачи^астяжения. Задача рас­ сматривалась в безразмерных переменных Tij,

T i j = P r \ T i j ,

щ =

(i , j = 1 , 2 )

Представим решение в виде суммы конечного числа членов тригоно­ метрического ряда:

~

=

N

[zn cos(na2) + z4N+4+n sin(no!2)],

Тц

zo + ^

 

 

7 1 = 1

 

_

 

 

N

 

T\ 2

=

ZN+\ + ^ 2

lz N + \+ n cos(na2) + Z5N+4+n sin(nQ!2)],

 

 

 

7 1 = 1

 

 

 

 

N

 

u \

=

Z2N +2 +

^ 2

[z 2N + 2+ n C O s (n a 2) + 2бДГ+4+п sin (n O :2 )] ,

 

 

 

71= 1

 

 

 

N

 

U2

=

Z^N + 3 + ^ 2

[2злг+з+п cos(na2) + Z7 N+4 +n sin(na2)] •

n= 1

В этом случае разрешающая система уравнения будет иметь вид

(/ = 0,..., 3, ц —0,... ,N, i2 = l, ... ,N)

3

rfZj(AT+l)+ii

da i

У " X i lJo Z j ( N + \ ) +

J = 0

N 3

+ E E K n 2 J ( N + l ) + n + ^ n ^ i V + S + J A T + n ] ,

71— 1 J = Q

11.1. Получение разрешающей системы уравнений

231

dzA ( N + l ) + I N + i 2

= y ^ , >Ch.OZJ(N+l) +

 

dct\

 

J = О

 

 

 

 

3 N

 

 

Wi2nzJ{N+\)+n + >Ci2l JZ4N+z+jN+n]

(1 1.2)

J=0n=l

Порядок полученной системы обыкновенных дифференциальных урав­ нений равен SN + 4. Коэффициенты системы (11.2) определяются со­ отношениями (если г = 0, то t = 1, иначе t = 2)

2 7 Г

Гх?°

х ?о

х 20

30 ]

=

_L

[р6,77>

9l,ni 9Ъ,т p4,n] COS(ZQ!2)C?Q!2I

|/* 7 Т 7 >

'S n '

 

^ i n j

 

^

 

 

 

 

 

 

0

 

 

 

 

 

 

 

2 7 7

 

Гх 91

х>*

X21

X?1!

=

[^6,т

k7 .n1 k3 .n1 &4.n] COS(ZQ!2)^Q:2I (11 -3)

[Л т '

■ 'Чп»

> г 7Т7>

Л гп J

 

^

 

 

 

 

 

 

0

 

 

 

 

 

 

 

2 7 7

 

Гк02

к 12

к 22

Ъ,321

[б7,П) ^8, m 64 n , e5,n] COs(iQ!2)rfQ!2>

 

 

^гп»

^ in J

 

^

 

 

 

 

 

 

0

 

 

 

 

 

 

 

2 7 Г

 

X 03

X 13 x 2 3

x ? 3

_1_

_

in>

^гп»

^ г п

tir

1

 

 

 

0

[f 7 . n 1 f 3 . n 1 f 4 . n 1 fb .n \c o s (ic x 2 ')d c x 2 ,

 

 

 

 

 

 

 

 

 

 

 

2 7 Г

 

 

 

 

 

 

 

 

^ 7 7 7

^ 7 7 7

 

 

 

-3 0

]

=

s

[^—6,n>

9—7.ni 9 ~3 .ni 9 4,n] sin(ia2)rfa2,

 

 

 

 

^ 7 7 7

 

 

 

X - 10

X - 20

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

]

=

2TT

 

 

 

 

 

: ^

01.

^ г п 11- ^ 7 7 1

 

>

х

гп

5

[^ —6,71)

k —7,ri)

k —2,.ni

k — 4 TJJ sin (iQ !2 )c^Q!2 j

 

 

 

 

х “ 21

 

- 31

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 7 7

 

 

 

 

 

x ~

m

’1

Х Г 12

- 2 2

 

 

- 3 2

 

 

 

[ б —7,Ti)

б _ 8 , П )

6 _ 4 n ,

e _ 5 , n ] s i l l ( 2012) б?0!2 >

 

 

^ 7 7 7

] = S

, 1 7 1

 

гп » ^ 7 7 7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

277

 

 

 

 

 

' x Z 0 3 ,

X - 13

х -

2 3

х - 3 3

 

 

 

[ / —7,Ti)

/ —8 ,Ti)

/ —4,Ti)

/ —5,TI] sin(2Q ;2) 6?C^2■

гп

■^777

^ 7 7 7

 

 

гп

 

 

27Г

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Задача изгиба. Теория Knpxj^£a--JlHBa. Задача рассматривалась

в безразмерных переменных N if М ^,

w (i ,j

= 1,2):

 

 

 

1Vi = P n N i,

 

 

М у = Р г ? М у ,

 

=

 

 

;U).

 

 

 

 

 

 

 

 

7

 

 

 

 

 

 

 

£ ch3

232 Гл. 11. Анализ НДС многослойных эксцентрических колец

Представим решение в виде суммы конечного числа членов триго­ нометрического ряда:

N 1

 

JVI +

N

 

 

zo + ^2 zn cos(na2),

 

*

 

 

 

 

 

 

A V& 2

 

 

 

 

П — 1

 

____ _

 

 

N

 

М и

=

z N+[ +

^ 2 z N + l+ n cos(nct2),

 

 

 

 

П=1

 

_ _

 

 

N

 

$ 1

=

z 2 N + 2 +

z 2 N + 2 + n C O s ( n a 2 ) ,

 

 

 

 

7 1 = 1

 

 

 

 

N

 

w

=

Z 3 N + 3 +

y ; гздг+3+n C O s ( n a 2 ) .

 

 

 

 

71= 1

 

В этом случае разрешающая система уравнения будет иметь вид

(7 = 0,

= 0,... ,N )

 

 

 

 

da i

N 3

(11.4)

 

 

ЕЕ ^ i n z J { N + \ ) + n -

 

dzi(N+l)+it

 

 

n=0 J=0

Порядок полученной системы обыкновенных дифференциальных уравнений равен 8N + 4. Коэффициенты системы (11.4) определяются соотношениями (если %= 0, то t = 1, иначе t = 2 )

lx?0

х ш

х 20

х 301 =

 

 

 

 

 

 

L zri’

^гп»

 

^ъп \

 

 

 

 

 

 

 

 

 

 

2 7 Г

 

Ai

 

~

 

 

 

 

 

 

 

 

 

 

cos(m2)^Q!2.

 

 

 

tir

—— cos(na2). ^4cos(no:2). 0* 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 n

 

 

 

 

Г „01

„11

21

31 1 _

1

 

C4,n

C4p,7l>

 

 

\ ^ i n '

^ i n

 

J

^

 

 

 

 

 

=-(сз,п - cos(na2)) - 2 ^ 0 2 ,n -

2c2p,n,

Ac\,n, 0

cos(ia2 )da2,

 

 

 

 

 

 

 

2TT

 

 

 

Г„02

„12

„22

„32 1

_

1

[

 

 

 

\_^in>^in'

^in> ^in \

 

^

\[^7.711^8,n> 64ri> e5n] cos(ia2 )dai2,

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

2 7 Г

 

 

 

[x93

x 13

x 23

x 331

=

L/V.n» Ув.ти / 4,TII

/ 5 ,71] COs(iO!2)dQ:2•

L

г п ’

-^гп *

 

J

 

^

0