книги / Прямые и обратные задачи механики упругих композитных пластин и оболочек вращения
..pdf10.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