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

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

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

12.1. Разрешающая система уравнений

253

Ненулевой компонент вектора свободных членов di = - t / 8 2.

Разрешающая система уравнений теории Григолюка-Чулкова

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

 

 

 

 

§ = Ш

 

 

 

 

 

 

(12.3)

где в качестве разрешающих выбраны функции

 

 

 

 

У0 = Т П,

Vl = Ni,

Ук+1 = Ф п

,

U N + 2 = Ti2,

V N + k + 2 =

^ 1 2 .

У2АГ+3 = U,

V 2 N + 4 =

W ,

V 2 N + k + 4

=

0 \

,

УЗАГ+5 =

V,

У ш + к + 5 = 0 2

Выпишем вспомогательные соотношения

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

к= 1

 

 

 

 

 

 

 

 

 

 

 

= 4

4

.

= a « 7 i,

c<*> = ag>~3

 

 

 

 

 

 

 

 

 

 

 

 

г]

'Ук >

 

 

f

= 4 j ] -

hk-<ag} + Y 1

Xlt5ij ’

 

 

 

 

 

 

 

 

 

 

1=1

 

 

 

 

 

f l - m ) = €

( ^ m ) - Ь п - ^ Z b W

-

 

 

 

 

+

 

 

 

 

+ Xkmlbg' - h t - i a g g + Хтк(Ц™' ~ Л т -1Oy"1) +

.

 

 

 

 

 

 

 

 

 

 

 

 

*=i

 

= 36(7fe —7fe(/i + /ifc + /ife-i) +7fc(^fc^ + ^fc-i^ + ^fc-i^fc) ~ Ikhk-ihkh)

°fc

 

 

 

/i3(/ifc- / i fc_,)3

 

 

 

 

Afc = I t (^27fc - 2^7fc + 7fc).

£fc =

4

(^7fc -

7fc).

 

 

/iD

 

 

 

 

 

 

 

 

 

 

 

 

 

^fcfc =

6____

 

T. = A _

*2*

 

 

 

 

5(/ifc —lifc-i) ’

 

fc

 

 

Afcfc ’

 

 

 

 

 

_Afc

A0fc

(fc)

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

Efc=i

Afcfc

 

 

 

^fcm _ _

TfcTm

.

cm, 1

Qij ~ Tij

------------2T +

°k

 

(ТЦТ22 -

r,22)

Afcfc

N

^

o(fc)

( g W g W - (gfj0)2)

N

e .. =

v ^ fc)

m= 1

fc=l

Вектор-функции правой части системы уравнений (12.3) имеют вид

9 о = - ( у о - Т2 2 )/г,

д \ = - у \ / г - Р ,

9N+2 = -2 у я + 2 /г .

9 2 N + 4 = -Ч?1.

254 Гл. 12. Анализ НДС многослойных круглых пластин и колец N

9k+1 =

~{Ук+1 —^22) /Г —922^1 "I"

(922У2ЛГ+г+4 —^{гУЗЛГ+1+б) >

 

 

 

 

 

 

 

г= 1

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

9N+k+2 =

-2yAT+fc+2/r + $?201 -

X ] №2/2ЛГ+г+4 + УпУЗЛГ+г+б) ,

 

 

 

 

 

 

 

г= 1

 

 

 

 

92N+3

&00У0 + ^01УИ+2 — S0E 22 - 0.5t?f +

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

+

^ 2

(^£'0(г+1)2/г+1 +

^0(ЛГ+г+1)УЛГ+г+2

S I QK 22) ,

 

 

 

 

г=1

 

 

 

 

 

 

У 2^+*:+4 =

&(к+1)0У0 +

^(fc+ l)iyjV + 2

~

S fc+ 1 ^ 2 2 +

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

+ ^ 2

(^(*+1)(*+1)УН-1 +

k(k+l)(N+i+l)VN+i+2 -

8^к+\)К г22) ,

 

 

г=:1

 

 

 

 

 

 

 

 

УЗЛГ+5

=

&10У0 + ^п У л Г+2 — Sl-^22 + У Ш + ь / г +

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

+

^ 2

(^1(Н -1)Уг+1

+

^1(Л^+г+1)УЛГ+г+2 ~

,

 

 

 

 

г = 1

 

 

 

 

 

 

УЗАГ+fc+5 — k(N+fc+1)0Уо + k(N+k+l)iyN+ 2

~ sN+k+\E22 + y3jv+fc+5/r +

N

 

 

 

 

 

 

 

 

 

 

 

^ у (^(^+&+1)(г+1)Уг+1 ”1" ^ ( N + k + l ) ( N + i + l ) y N + i + 2

~ s i ( N + k + 1) ^ 22) >

г= 1

 

 

 

 

 

 

 

 

 

 

 

где:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

=

2/1 / (е 22 + Уо) + У^(<?22У2ЛГ+г+4 —912УЗАГ+г+5)/(е22 + Уо).

 

 

 

 

 

г= 1

 

 

 

 

 

 

 

т 22

=

< 2(у2лт+з + 0.5t9f) + а\2 Е 22

+ а*2ьЕ 12 +

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

+ Y 1 <%292N+l+4 +

К 22 + 4 б К \2,

 

 

 

 

 

 

г=1

 

 

 

 

 

 

 

« 4

=

 

а й > ( * 2JV+3 + 0.5я?1) +

^ 2 2 ^ 2 2

+ ^ 2 6 ^ 1 2

+

 

 

 

 

+ Е

ЛгН W

i + 4 +

ЙгЫ>^ 2

+ Ш

]К пк ,

 

 

 

 

г= 1

 

 

 

 

 

 

 

 

 

^22 = У2ЛГ+з/г,

i^22 = У2ЛГ+г+4/ г ,

 

 

^ 1 2 = УЗЛГ+5 - УЗАГ+б/г.

ЛГ}2 = УЗАГ+г+5 ~ УзЛГ+г+ 5 / г >

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

255

где \kij \ = \ki:j\

коэффициенты матрицы kij

 

^оо = a ii.

^oi =

= a*6,

fcn =

ag6, fyi-HKj-H) =

 

k ( i + \ ) ( N + j + \ )

= ^(г+1)(ЛГ+.;'+1)

— 7 j6 ^ *

k ( N + i + l ) ( N + j + l ) ~

/ 5 6

12.2. Влияние выбора структурной модели КМ на расчетные характеристики НДС и уровень нагрузок начального разрушения

В рамках теории ломаной линии Григолюка-Чулкова [129] рас­ сматривается задача определения НДС армированной однослойной кольцевой боралюминиевой пластины (Е с = 70,0 ГПа, Е а = 363 ГПа, а* = 0,55 ГПа, сг* = 3,1 ГПа) при использовании различных моделей композиционного материала: МОВ, МОВУ, МДВ и МБГ.

Рис. 12.2

На рис. 12.2 показана зависимость компонент матрицы жесткости ||ajj|| и матрицы поперечной сдвиговой податливости \\ра/з\\ = ||<7а/з ||-1 от угла укладки арматуры. При ш = 0,59, a>z — 0,92, а для МБГ

256 Гл. 12. Анализ НДС многослойных круглых пластин и колец

 

 

 

 

 

 

 

 

Т а б л и ц а 12.1

Модели

 

Коэффициенты моделей КМ (-10

9)

 

КМ

ан

Oi2

0 2 2

013

0.23

Озз

ри

Pi2

Р22

 

 

 

Угол а эмирования 0

°

 

 

 

МОВ

2 0 2 , 2

1,9

6 , 2

0

0

2 , 2

336,5

0

336,5

МОВУ

228,8

1,9

6 , 2

0

0

2 , 2

336,5

0

336,5

МДВ

240,9

42,2

140,6

0

0

49,2

77,8

0

47,7

МБГ

240,5

40,9

136,4

0

0

57,9

57,9

0

47,7

 

 

 

Угол армирования 30°

 

 

 

МОВ

116,4

38,6

18,4

63,7

2 1 , 2

38,9

336,5

0

336,5

МОВУ

131,4

43,6

2 0 , 1

72,3

24,1

43,9

336,5

0

336,5

МДВ

197

61

146,9

32,6

10,9

6 8

70,3

13

55,2

МБГ

2 0 2 , 6

52,8

150,5

29,4

15,7

69,8

55,4

4,4

50,3

 

 

 

Угол армирования 45°

 

 

 

МОВ

55,2

50,9

55,2

49

49

51,2

336,5

0

336,5

МОВУ

61,8

57,5

61,8

55,7

55,7

57,8

336,5

0

336,5

МДВ

165,7

67,3

165,7

25,1

25,1

74,3

62,8

15

62,8

МБГ

172,6

56,7

172,6

26

26

73,8

52,8

5,1

52,8

 

 

 

Угол армирования 60°

 

 

 

МОВ

18,4

38,6

116,4

2 1 , 2

63,7

38,9

336,5

0

336,5

МОВУ

2 0 , 1

43,6

131,4

24,1

72,3

43,9

336,5

0

336,5

МДВ

146,9

61

197

10,9

32,6

6 8

55,2

13

70,3

МБГ

150,5

52,8

2 0 2 , 6

15,7

29,4

69,8

50,3

4,4

55,4

Угол армирования 90°

МОВ

6

, 2

1,9

2 0 2 , 2

0

0

2

, 2

336,5

0

МОВУ

6

, 2

1,9

228,8

0

0

2

, 2

336,5

0

МДВ

140,6

42,2

240,9

0

0

49,2

47,7

0

МБГ

136,4

40,9

240,5

0

0

57,9

47,7

0

336,5

336,5

77,8

57,9

объемная интенсивность укладки равна 0,54. Здесь и далее сплош­ ные линии 1 соответствуют МОВ, пунктирные линии 2 — МОВУ, штрихпунктирные линии 3 — МДВ, штрихштрихпнуктирные ли­ нии 4 — МБГ.

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

257

В табл. 12.1 приведены числовые значения Ца^Ц и ||ра/з||- Из табл. 12.1 и рис. 12.2 видно, что значения коэффициентов, рас­

считанных по МДВ и МБГ,

близки,

как и коэффициентов по МОВ

и МОВУ.

 

 

Для коэффициентов ап,

а 12, а22,

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

МОВ и МОВУ, меньше, чем рассчитанные по двумерным моделям. Особенно эта разница заметна для ап при углах близких к 90° и для а22 при углах близких к 0°. Это объясняется тем, что МОВУ не учитывает слагаемые, отвечающие за работу волокна в поперечном направлении, а МОВ, кроме того, не учитывает работу связующего между волокнами.

В табл. 1.1 приведены экспериментальные данные для эффективных модулей Ei, Е 2 для заданных параметров. Видно, что различие между экспериментальными данными и полученными оценками Е\ составляет около 17% для МОВ, 7% для МОВУ, а для МДВ и МБГ не более 2%; в случае Е 2 97% для МОВ и МОВУ, а для МДВ и МБГ не более

2 % .

Значения компонент матрицы поперечной сдвиговой податливости, рассчитанные по одномерным моделям КМ, значительно больше, чем рассчитанные по МДВ и МБГ.

Жестко защемленная на контурах пластина, нагруженная рас­ пределенным внешним давлением. На рис. 12.3-12.9 приведены за­ висимости характеристик НДС от радиуса жестко защемленной на контурах пластины для различных углов укладки арматуры при на­ гружении равномерно распределенным внешним давлением. Параметры расчета: h = 0,01 м, г0 = 0,03 м, г\ = 0,3 м, р = 106 Па.

Для углов армирования близких к 0° разница между напряжениями а *1 = сгц/р, рассчитанными по различным моделям, мала (рис. 12.3), так как при совпадении направлений армирования и нагружения учет работы волокна и связующего происходит похоже для всех моделей. Существенное различие наблюдается для углов близких к 90°, так как в данном направлении для одномерных моделей не учитывается работа

волокна.

 

Для углов укладки близких к 0° значения напряжений

= <т22/р

больше при использовании двумерных моделей (рис. 12.4),

так как

в окружном направлении для данных углов МОВ и МОВУ не учитыва­ ют работу волокна. Для углов укладки близких к 90° значения напря­ жений по одномерным моделям значительно превышают напряжения по двумерным. Это связано с тем, что для данных углов компоненты матрицы деформации для одномерных моделей существенно больше компонент, рассчитанных по двумерным моделям (рис. 12.5).

На рис. 12.6 показана зависимость напряжений <7J*2 = <т12/р. Видно, что для данного напряжения значения, рассчитанные с использованием одномерных моделей, превышают рассчитанные по двумерным моде­ лям.

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

258 Гл. 12. Анализ НДС многослойных круглых пластин и колец

Аналогичная ситуация складывается и для поперечных напряжений (рис. 12.7-12.8).

На рис.12.9 приведена зависимость прогибов в пластине от ради­ уса. Минимальные прогибы соответствуют укладке арматуры вдоль радиуса. Наибольшая разница между прогибами для различных углов наблюдается между значениями по МДВ и одномерным моделям, эта разница значительно меньше для МБГ и МДВ. При использовании МОВ и МОВУ наиболее сильно выражено влияние углов укладки арматуры.

Важной характеристикой НДС пластины является интенсивность напряжений в компонентах КМ. В рамках МБГ можно получить

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

259

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

Видно существенное влияние выбора модели КМ при углах укладки близких к 90°. Уровень напряжений в связующем, рассчитанный по м д в , значительно ниже рассчитанного по одномерным моделям, так как в МДВ нагрузка в связующем дополнительно воспринимается этим компонентом между волокнами в армированном слое.

На рис. 12.11 показана зависимость интенсивности напряжений в арматуре для различных углов армирования. Для углов ~ 0° все модели дают близкую оценку, но чем больше угол армирования, тем больше различие между результатами полученными по различным

9:

260 Гл. 12. Анализ НДС многослойных круглых пластин и колец

Рис. 12.5

моделям. Начиная с 60° виден разный характер распределения напря­ жений по МДВ и МБГ. Это связано с тем, что при МБГ уровень напряжений в волокне восстанавливается в результате предположения о том, что волокно работает как нить (аналогично предположениям, принятым при моделировании волокна в МОВ).

Жестко защемленная на внутреннем контуре пластина, нагру­ женная распределенным внешним давлением и растягивающим усилием по внешнему контуру. Равномерно распределенное давле­

ние р = 5

• 104 Па, растягивающая нагрузка по внешнему контуру

Т = 5 - 104

Па/м.

Величины напряжений и прогибов, рассчитанные по одномерным моделям, близки, как и результаты, полученные для двумерных мо­ делей КМ. Относительная разность напряжений не превышает для одномерных моделей 1-3%, для двумерных моделей 3-5% . Различие прогибов для одномерных моделей составляет 1012%, для двумерных не более 1 %. Значения кинематических характеристик, рассчитанных по одномерным моделям, для данного КМ и вида нагружения превы­ шают полученные в рамках МДВ и МБГ. Показано, что влияние угла укладки арматуры гораздо более выражено для одномерных моделей.

Для данного варианта нагружения приведены зависимости интен­ сивностей напряжений в матрице (рис. 12.12) и арматуре(рис. 12.13). Хорошо видно влияние выбора моделей КМ: уровень напряжений

12.3. Влияние выбора теории пластин

261

в связующем, рассчитанный по МДВ, значительно ниже рассчитан­ ного по одномерным моделям, так как в МДВ нагрузка в связующем дополнительно воспринимается этим компонентом между волокнами

вармированном слое.

Вслучае сложного нагружения различие между результатами, по­

лученными по одномерным моделям и МДВ с МБГ, значительно для всех углов армирования (см. рис. 12.13). При увеличении угла армирования возрастает разница между напряжениями, полученными по МДВ и МБГ. Это связано с методом определения напряжений в армирующих волокнах.

Таким образом, для данных видов нагрузок и данного компози­ ционного материала расчет по моделям с одномерными волокнами дает результаты оценки НДС «с запасом» При оценке необходимости использования более сложных моделей следует учитывать направление наибольшего нагружения. Если основная нагрузка будет приложена вдоль волокна, то одномерные модели не дадут величин, сильно от­ личающихся от полученных по двумерным моделям, в иных случаях различие может доходить до 80%.

12.3. Влияние выбора теории пластин на расчетные характеристики НДС и уровень нагрузок начального разрушения

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

Влияние выбора теорий на распределение компонент тензора напряжений по толщине пластины. Рассматривается трехслойная уг­ лепластиковая пластина со следующими механическими параметрами: Е с = 3 • 109 Па, Е а = 300 • 109 Па, vc = va = 0,3, а* = 0,09 • 109 Па, а* = 3- 109 Па.

На рис. 12.14-12.17 приведены графики распределения напряже­ ний по толщине для трехслойного кольца в сечении г = 0,25г! при го/г 1 = 0,2. Кривые 1 представляют значения, рассчитанные по теории Кирхгофа-Лява, 2 — Тимошенко, 3 — Андреева-Немировского, 4 — Григолюка-Чулкова. Результаты, полученные послойно для структуры армирования (±15°, 90°, ±15°), отображены сплошными линиями, для структуры (90°, ±15°, 90°) штриховыми.

На рис. 12.14 представлено распределение напряжений для жестко защемленной на внутреннем контуре пластине, нагруженной по внеш­ нему контуру растягивающим усилием То = 3 • 105 Н/м.

Различие напряжений сгц, 022, полученных по различным теориям, не превышает 1 %. Для касательных напряжений п з в рамках теории

Кирхгофа-Лява

не

удается получить

оценку, в

теории Тимошенко

не выполняется

их

равенство нулю на

свободной

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

262 Гл. 12. Анализ НДС многослойных круглых пластин и колец

удовлетворяются условия межслоевого контакта. Уровень напряжений Ti3 по теории Тимошенко близок к нулю, что составляет 0,001 % от максимального уровня напряжений п з в теории Григолюка-Чулкова.

Оценка уровня касательных напряжений при использовании гипо­ тезы ломаной линии в крайних слоях существенно выше, чем при гипо­ тезе Андреева-Немировского. Для среднего слоя оценки, полученные по обеим теориям, близки.

На рис. 12.15 представлено распределение напряжений для жестко защемленной на внутреннем контуре пластины, нагруженной равно­ мерно распределенным давлением интенсивности р = 2,5 • 104Н/м2.

Различие

напряжений сгп, 022

полученных по различным тео­

риям, снова

не превышает 1 %.

Теории Андреева-Немировского

и Григолюка-Чулкова дают достаточно близкую оценку касательных напряжений, однако зависимость, полученная при использовании теории ломаной линии, имеет несколько локальных максимумов. Изменение структуры армирования приводит к исчезновению с графика Ti3 локальных максимумов во внешних слоях.

На рис. 12.16 приведено распределение напряжений для жестко за­ щемленной на внутреннем контуре пластины, нагруженной равномерно распределенным давлением интенсивности р = 2,5 • 104 Н/м2 и нагру­ женной по внешнему контуру растягивающим усилием Т0 3 • 105 Н/м.