
- •механика
- •материалов
- •СТРУКТУРНАЯ МАТЕМАТИЧЕСКАЯ МОДЕЛЬ РЕОЛОГИЧЕСКОЙ ДЕФОРМАЦИИ ДВУХКОМПОНЕНТНОГО МАТЕРИАЛА
- •ПОДАТЛИВОСТЬ ОДНОНАПРАВЛЕННО АРМИРОВАННОГО НЕУПРУГОГО МАТЕРИАЛА*
- •ПРОЧНОСТИ
- •РАСЧЕТ ЧАСТОТ СОБСТВЕННЫХ КОЛЕБАНИЙ АНИЗОТРОПНЫХ ПЛАСТИН С ОТВЕРСТИЯМИ*
- •УСТОЙЧИВОСТЬ ПЛАСТИН ИЗ КОМПОЗИТОВ ПРИ СДВИГЕ В ПЛОСКОСТИ ДЛЯ СМЕШАННЫХ ГРАНИЧНЫХ УСЛОВИЙ
- •РАЗРЕШИМОСТЬ И ОЦЕНКА СОБСТВЕННЫХ ЧИСЕЛ СИСТЕМЫ НЕЛИНЕЙНЫХ УРАВНЕНИЙ ГРИГОЛЮКА—ЧУЛКОВА
- •МНОГОКРИТЕРИАЛЬНОЕ ПРОЕКТИРОВАНИЕ ТОНКОСТЕННЫХ ЭЛЕМЕНТОВ КОНСТРУКЦИЙ ИЗ КОМПОЗИТНЫХ МАТЕРИАЛОВ
- •ОБОБЩЕННАЯ МОДЕЛЬ МЕХАНИКИ НАМОТКИ ОБОЛОЧЕК ИЗ КОМПОЗИТНЫХ ПОЛИМЕРНЫХ МАТЕРИАЛОВ
- •РАВНОВЕСИЕ НИТИ С УЧЕТОМ ТРЕНИЯ ПРИ ХОРДОВОЙ НАМОТКЕ ДИСКОВ ИЗ КОМПОЗИТОВ
- •НЕСУЩАЯ СПОСОБНОСТЬ ПОЛЫХ СТЕРЖНЕЙ ИЗ КОМПОЗИТОВ ПРИ КРУЧЕНИИ
- •ВАРИАБЕЛЬНОСТЬ ДИНАМИЧЕСКОЙ ПРОЧНОСТИ ПОЗВОНОЧНИКА ЧЕЛОВЕКА
- •ВЛИЯНИЕ ФИЗИЧЕСКОГО СТАРЕНИЯ НА ДЕФОРМАТИВНЫЕ СВОЙСТВА ПОЛИКАРБОНАТА
- •J (**—T)miaai(mi+1>(x)c*T=l,
- •ФИЗИЧЕСКИЕ СВОЙСТВА НЕУПОРЯДОЧЕННЫХ СТРУКТУР
- •ПРОГНОЗИРОВАНИЕ ДЛИТЕЛЬНОГО СОПРОТИВЛЕНИЯ ПОЛИМЕРНЫХ МАТЕРИАЛОВ
- •ВЫСОКОМОЛЕКУЛЯРНЫЙ ПОЛИЭТИЛЕН
- •высокой плотности
- •ПОЛИИМИДЫ — НОВЫЙ КЛАСС ТЕРМОСТОЙКИХ ПОЛИМЕРОВ
У Д К 6 2 4 .0 7 3 .0 0 1 :6 7 8 .0 6 7
В. А. Поляков, В. В. Хитрое
УСТОЙЧИВОСТЬ ПЛАСТИН ИЗ КОМПОЗИТОВ ПРИ СДВИГЕ В ПЛОСКОСТИ ДЛЯ СМЕШАННЫХ ГРАНИЧНЫХ УСЛОВИЙ
Современные методы проектирования элементов конструкций типа тонкостенных панелей, изготовленных из слоистого композитного мате риала, основаны на всестороннем анализе решений задач прочности и устойчивости с учетом анизотропии свойств материала. Для последних принципиальным является вопрос о допустимом соответствии жесткости и размеров изделия. Использование высокомодульных волокон и техно логичность процесса изготовления тонкостенных пластин позволяют реа лизовать для них в широких пределах необходимые деформационные свойства.
Большая часть работ, относящихся к расчетной методике по устойчи вости тонкостенных анизотропных пластин, выполнена в предположении однородности свойств материала, причем преобладают исследования на случай нормальных нагрузок в плоскости пластины. Расчет при каса тельных усилиях, распределенных по контуру, проводился энергетиче ским методом в работе [1] и статическим методом в [2]. Обзор результа тов этих работ и анализ решения, полученного для полосы при сущест венной анизотропии материала, производится в работе [3]. Центральным вопросом в решении этих задач является выбор функции для прогиба при потере устойчивости пластины. Это обусловлено необходимостью рассмотрения ряда гармоник [1] или введения в аргумент тригономет рических функций угла наклона косых волн [2] с тем, чтобы обеспечить нетривиальную форму прогиба при изгибе от сдвига.
Смешанная задача, когда две противоположные стороны пластины защемлены, а две другие жестко заделаны, в известной авторам литера туре не рассматривалась. Такой случай, однако, может считаться ком промиссным при формулировке граничных условий по ребрам и торцам скручиваемого стержня замкнутого профиля, изготовленного из прямо угольных панелей (рис. 1—а).
В настоящей, работе рассматривается задача определения критиче ского значения касательных усилий для прямоугольной ортотропной пластины со смешанными граничными условиями. Исследуется зависи мость нагрузки от параметров анизотропии, характерных для продольно поперечной намотки и косоармированных композитных материалов, при различных геометрических размерах пластины.
Пусть прямоугольная ортотропная пластинка, толщиной h и разме рами в плане а и Ъ (рис. 1—б), деформируется касательными усилиями t = hxxy, распределенными равномерно по всем четырем сторонам. Глав ные оси упругости параллельны сторонам; стороны х = 0 и х = а жестко
Рис. 1. Кручение заделанного по торцу коробчатого стержня
(а) и схема касательных уси лий и размеры пластинки — грани стержня (б).
заделаны, а стороны у = 0, у=Ь свободно оперты. Критическое значение /1ф, при котором пластинка теряет устойчивость, будем определять энер гетическим методом. Согласно этому методу [3], следует приравнять по тенциальную энергию пластинки при изгибной форме равновесия П ра боте касательных усилий А при деформации сдвига, обусловленного из гибом пластины. При представлении прогиба в виде функционального ряда критическое значение нагрузки находится из условия вырожден
ное™ детерминанта следующей системы уравнений [3]: |
|
д с г = о; |
w |
где £/=П —A; Ci — коэффициент при i-й функции прогиба в его разло жении в ряд. Для смешанных граничных условий
№ = 0, |
dW |
=0 при х = 0,х = а; W=0 при у = 0, у=Ь (2) |
примем следующую двухчленную форму представления прогиба:
W'=Cn(p„(.v') sin^|- + Cn+icpn+1 (х) sin » (3)
в которой функции фп(х), л = 1,2 ,... являются собственными функциями уравнения заделанного стержня и, как отмечено в [4], были с этой целью использованы Ритцем. Эти функции удовлетворяют уравнению
dAфп |
—Яп4фп-, |
л — 1 , 2 , . . . |
(4) |
|
dx4 |
||||
|
|
|
||
и нормированы в промежутке [0, а]: |
|
|
||
а |
|
|
|
|
| (pn2dx=l; |
л = 1 , 2 , |
(5) |
||
о |
|
|
|
|
Запись их имеет следующий вид: |
|
|
||
(cos Хпх — ch Хпх) (sin pin—sh р„) —(sin Хпх — sh Хпх) X |
|
|||
Ф„ (x) |
X (cos р„ —ch pn) |
(6) |
||
] a (sin pn—sh pn) |
||||
|
|
|||
|
л = 1,2, 3,. |
|
||
где собственные числа Xn и p„ = aXn определяются из уравнения |
|
|||
|
cos pin ch p„= l. |
(7) |
Условия (2) при выборе прогиба в форме (3) с учетом (6) выполняются автоматически*.
При построении решения укажем [4], что первый корень уравнения 3
(7) имеет значение pi=^ л+0,0176, а для последующих его значений
справедлива с точностью более 0,1% следующая приближенная фор мула:
2л+1 |
(о) |
р « = — ^— л; п=2>3’ |
Работа касательных усилий, соответствующая малым прогибам, под
считывается с учетом (3) по формуле |
|
|
|
а Ь |
4 |
а |
|
Г f d\V dW |
Г |
(9) |
|
A = t J J |
dxdy = — tCnCn+i J (cpncp^+i —fp/n(pn+i) dx. |
* Функции фп имеют размерность [Z.-1/:], так что коэффициенты С„, Сп-и имеют со ответственно размерность [UU\.
Прежде чем вычислить интеграл в правой части (9), установим, опи раясь на условия (2), (4), (7), некоторые соотношения между интегра лами от произведений собственных функций ф„(х) и их производных. Используя уравнение (4) и условие фп|*=о,а = ф/п|*=о,а=0, которое вы полняется в силу граничных условий (2), получим после трехкратного интегрирования по частям следующее интегральное соотношение:
а |
|
а |
|
V J фпф'гЖ^ + Атж4 | |
ф»+1ф'п<*Х= [ф,,яф/'п+1]® . |
(Ю) |
|
0 |
0 |
|
|
Составим также очевидное в силу условий (2) соотношение
аа
J ф п ф ' п + l ^ + J |
ф п + 1 ф ' п ^ Л ; = 0 . |
( 1 1 ) |
0 0
Решая совместно (10), (11) относительно интегралов, получим
Г |
, |
, |
- |
Г , |
. |
[ф//пф//п+1]а |
(12) |
||
J фпф n + \ d x = |
J |
ф n * f n + l d x = |
■ |
------ т— . |
|||||
« |
|
|
|
п |
|
|
Ап |
+1 |
|
Расчет граничных значений вторых производных собственных функ ций (6) с учетом (7), (8) сводится к соотношениям
// |
2JW2 |
„ . . |
2А.П2 cos |
sh p „ -sin |
р„ ch рп |
|
Ф п (0) = ------ |
— |
; ф „ (а) = |
---у ---------------------a |
;-------------------------------- |
|
|
|
Уа |
|
sin |.iri—sn цп |
|
||
- |
2Ъп2 ( - l) n sin p „ - tg р„ |
= ( - I )t .n2 K |
2 |
(13) |
||
|
уa |
sin p„ —(—l) n tg p„ |
Уа, |
|
При записи последнего выражения в (13) было использовано следующее из анализа (7) соотношение sh р„= ( —1)" tg р„; л = 1 ,2 ,... Таким обра зом, окончательная запись выражения (12) имеет вид
о
и работа касательных усилий согласно (9) выразится формулой
64 Яп2Ап+12
с пс 71+1 t.
За Кп2—Я.п+12
Потенциальная энергия изгиба при потере устойчивости пластины имеет запись
аb
П = - М |
J Ф (x,y)dxdy, |
(14) |
z о |
о |
|
в которой для подынтегральной функции в случае ортотропии материала справедливо выражение
Ф= £>1 |
d2W дПУ |
+ 4Dh |
(15) |
+ 2DiVi2 дх2 ду2 + 7)2 |
где цилиндрические жесткости определяются через упругие константы материала в главных осях упругой симметрии:
' |
/г3 |
Ег о |
/23 |
(16) |
|
D\,2= —г* ----- :— ; |
7)/1=— -G 12; D^=D\v\2 + 2Dk. |
||
|
12 |
1—V12V21 |
12 |
|
Здесь h — толщина пластины. В целях упрощения последующих резуль татов последним выражением в (16) приведена комбинированная жест
кость, которая в случае изотропии материала вырождается, так что
^3 = ^1,2- Подставив (3) в (15), при интегрировании по (14) используем соот
ношение (5) и следующие из него и (2) интегральные зависимости для собственных функций:
|
|
|
|
а |
|
|
|
|
|
а |
<pn2dx=Xn4; |
п=1 ,2 ,... |
|
|
|
|
|||
|
|
|
J |
(p"2ndx= Xn4 j |
|
|
|
(17) |
|||||||||||
|
|
|
|
о |
|
|
|
|
|
о |
|
|
|
|
|
|
|
|
|
|
а |
ф'2ndx= - |
а |
ф„ф,/„^=J |
|
|
|
|
|
|
|
|
|||||||
|
J |
|
Л „2Хп2- |
; |
п= 1 ,2 ,... |
|
|
(18) |
|||||||||||
|
о |
|
о |
|
|
|
|
|
|
|
|
|
а |
|
|
|
|
|
|
г, |
о |
|
моч |
параметр |
|
|
COS(X„-ch(X„ |
|
sin (An |
|
|
от- |
|||||||
Введенный в (18) |
хп= —--------- ;----- = |
—;— гг-------------- |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
sin pn -sh р.„ |
( —1)“ —cos |
|
|
|
|||
личается от единицы менее чем на 0,1% при |
2. |
|
|
|
|
||||||||||||||
В результате, учитывая соотношения (17), (18), получим интегриро |
|||||||||||||||||||
ванием по формуле (14) следующее выражение для потенциальной энер |
|||||||||||||||||||
гии пластины: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
п |
|
I f ^ |
|
D i b |
|
|
It D 2n 4 |
, / ) 3 ( р п 2Х п 2 - 2 р , п > С п ) п 2 |
|
|
|
||||||||
п = |
2 VС4 |
1 |
? |
•*" + 1 5 Г |
+ -----------М,------------ |
|
|
|
|||||||||||
|
|
|
|
Цп+14 + |
16Z)2jt4 |
|
4£>з(|хп-н2Хп+г |
2|in+i^n+i) л- |
|
|
|
||||||||
|
|
|
|
~ 2 Ь ~ + |
|
аЧ> |
|
|
|
|
|||||||||
Приравнивая детерминант системы двух уравнений (1) |
при i=n, |
л+1 |
|||||||||||||||||
нулю, получим окончательное выражение для критического значения |
|||||||||||||||||||
сдвигающего усилия при потере устойчивости пластины: |
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
tкр — |
|
3 |
Цп+14 — Цп4 l / D l D j |
|
|
|
|
|
|||||
|
|
|
|
|
|
128 |
|
р„2р„+12 |
Ь'2 |
|
|
|
|
|
|||||
X VivD2 С3 |
|
|
|
|
|
|
|
|
|
||||||||||
П |
|
' Di |
|
|
С |
|
yDiD2 -I |
х |
|
||||||||||
|
|
|
D |
1 Цп4 |
| |
|
|
у |
Д 2 | |
2 я 2(р п2Хп2- 2 р п К п ) |
D i |
1 |
|
||||||
X |
ГDl ] |
/ |
Ц |
п1+ ‘41 |
6 |
л |
' |
4с° 21 |
/ | |
8 |
n 2 |
( P n + |
l2K n + l2 |
- |
2 |
p n + l X |
|||
' |
D2 |
С3 |
’ |
Л |
|
Di |
|
|
С |
T/nI/D1.DnJ2 ■’ |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(19) |
где для относительной длины принят параметр С—а/Ь.
Следует иметь в виду, что в зависимости от порядка собственных функций п критическая нагрузка, найденная по формуле (19), будет со ответствовать различающимся между собой изгибным формам потери устойчивости пластины. При я=1 прогиб W изменяется вдоль длины пластины согласно двум функциям — ф] и ф2, первая из которых описы вает одну полуволну, вторая — две с нулем посередине длины. С увели чением п число полуволн пропорционально возрастает.
Для одной и той же нагрузки ^кр при соответствующем соотношении размеров пластины в плане могут существовать две смежные формы по верхности при потере устойчивости. Значение параметра С=Сппр, при котором происходит изменение формы поверхности, соответствующей функциям ф„, фп+ь к форме поверхности, соответствующей функциям
фп+ь фп+2, определяется из уравнения |
|
||
^ t |
^ |
4)2 Рп (С) Qn+l (С) = - - n+2V n+;4)2 Рп+i (С) Q„+2(С), |
(20) |
Р п |
JIn+1* |
(ln + lV n + 24 |
|
где полиномы имеют вид
Р п { С ) = п * ~ \ 1 — - С * + |
2 л 2( ц п2х п2 — 2 ц пХп ) D 3 |
|
д а т |
||
1 Dt |
Qn(C) = 16я4 ] / - ~ с ^ + 8л2 (цп2хп2 —2Ц п Х п )D ifD\D2
Между двумя значениями С„,_I"P и Сппр находится такое значение С*„, при котором критическая нагрузка, приводящая к потере устойчивости пластины по форме, описываемой двумя собственными функциями срп и фп+1, будет минимальной*. Это значение относительно параметра опреде
ляется как С*п= Ух, где х — наименьший вещественный корень урав нения
4 / D \ 1 6j.ln^ “1“ P'n+l^ Z?32 (Ц?1^%712 2(Ll7lX?l) ( Ц7? + 1“%7? + 1 2 (ln + lX77-f 1) \ 2
* _ \ ~D 2 |
+ z v |
^ |
/ * " |
D3DI |
|
|
|
4 n W [ ( |
2|XnX?i) [In-f-14 ”T 4 (|l7i + l“%?i+ l“ |
2(.ln + iX n + l ) |Я ^ ]Х |
|
|
3 |
ZV |
|
|
16я8 |
Z)22Цп4Р'?г+14 —0. |
|
Сравнение расчетных значений критической нагрузки с результатами, полученными при использовании приближенных формул в [1], показы вает, что выражение (19) дает несущественное ее отличие (порядка 3—5%) в случае поперечного армирования пластинки, т. е. при укладке всех волокон в направлении оси у. Этот случай соответствует предпо сылке о существенной анизотропии материала, принятой в [1]. Вырож дение (19) на случай изотропии материала дает завышенные оценки для ^кр по сравнению с приведенными в [5] (отличие порядка 25%). Для по следнего случая, очевидно, недостаточно двухчленного представления прогиба в форме (3). Таким образом, выражение (19) можно рекомендо вать для качественного анализа изменения критической нагрузки при потере устойчивости ортотропной пластинки от сдвига.
На рис. 2 представлены кривые изменения приведенных критических напряжений пластин из углепластика в зависимости от относительной длины С= а/Ь для различной анизотропии свойств материала. Построе ние расчетных кривых приводилось при минимизации тКр по порядку соб ственных функций п. В случае продольно-поперечной укладки волокон в углепластике (см. рис. 2—а) изменение анизотропии упругих свойств характеризуется параметром й2, равным относительному содержанию волокон, уложенных в направлении оси у. В этом случае монотонное из менение цилиндрических жесткостей таково, что между их граничными значениями имеется следующая связь: D\(й2 = 0) =£>2(й2 = 1) и £>3(й2) =
= const. В случае косоугольной |
укладки волокон под углами ±0 (см. |
рис. 2—б) для цилиндрических |
жесткостей справедливо соотношение |
Z)I (0 = O)=D2 |0 = ту"| = £ I № = 0); н о D3(0) ^const. С увеличением доли
волокон в поперечном направлении у или с увеличением угла 0 критиче ская нагрузка возрастает в несколько раз. Так, изменение Ткр при С < 3 от смены укладок с продольной на поперечную составляет ~400% . Ана логичная тенденция имеет место и для случая свободного опирания по всему периметру пластины (см. [2]). Отметим, что в обоих случаях (см. рис. 2) минимизация нагрузки при переборе порядка функций п обнару-
* В случае n= 1 левая граница С0пр условно принимается равной нулю, а не нахо дится из уравнения (20).
Рис. 2. Влияние относительной длины с=а)Ь пластин и схемы укладки однонаправлен ного углепластика (£ц= 13,2 • 105 кгс/см2, £ ± =1,22*105 кгс/см2, “ 0,25 • 105 кгс/см2,
v llJL = 0,33) на приведенные критические напряжения: а — ортогонально армированный материал с различной долей й2 поперечных волокон; б — косоармироваиный материал с
различным углом армирования ±0 по отношению к оси х. |
1—7 — значения п. |
Л2 = 0 (I); |
0,2 (II); 0,4 (III); 0,8 (IV); 1,0 (V). 0 = 15 (I); 30 (И); |
45 (III); 60 (IV) и |
75° (V). |
живает немонотонное изменение Ткр по степени анизотропии упругих
свойств материала |
(см. кривые на рис. 2 при йг= 0,8 и |
0= ±6О°). |
|
Сопоставление |
критических напряжений |
типовой |
пластины (h/b = |
= 0,01) из косоармированного углепластика |
с приведенными в [5] дан |
ными по прочности при сдвиге подобного материала в плоскости арми рования 700—1200 кгс/см2 показывает, что в широком диапазоне 0 не сущая способность пластин будет определяться их жесткостью.
Таким образом, анализ влияния анизотропии позволяет оценить опас ность разрушения от потери устойчивости исследуемых пластин при сдвиге.
СП И С О К Л И Т Е Р А Т У Р Ы
1.Балабух Л. И. Устойчивость фанерных пластинок. — Техника воздуш. флота, 1937, Ко 9.
2.Секерж-Зснькович Я. И. К расчету на устойчивость листа фанеры, как анизотроп ной пластинки. — В кп.: Тр. ЦАГИ, 1931, № 76.
3.Лехницкий С. Г. Анизотропные пластинки. М., 1957. 463 с.
4.Канторович Л. В., Крылов В. И. Приближенные методы высшего анализа. М.; Л. 1962. 708 с.
5. Работное Ю. Н., Данилова Н. Н., Полилов А. Н., Соколова Г. В., Карпсйкин И. С., Вайнбсрг М. В. Исследование прочности намоточных эпоксидных угле- и стек лопластиков при кручении, растяжении и поперечном изгибе. — Механика полимеров, 1978, № 2, с. 219—225.
Институт механики полимеров |
Поступило в редакцию 22.01.82 |
АН Латвийской ССР, Рига |
|
М Е Х А Н И К А К О М П О З И Т Н Ы Х М А Т Е Р И А Л О В , 1982, Лг° 5, с. 839— 843
У ДК 624.074.001:678.067
В. П. Трошин
ВЛИЯНИЕ ПРОДОЛЬНОГО РАССЛОЕНИЯ в СЛОИСТОЙ ЦИЛИНДРИЧЕСКОЙ ОБОЛОЧКЕ НА ВЕЛИЧИНУ КРИТИЧЕСКОГО ВНЕШНЕГО ДАВЛЕНИЯ
Одним из основных дефектов, присущих слоистым композитам, яв ляется расслоение. Возникнув в результате несовершенства технологии изготовления или воздействия случайных факторов при эксплуатации, расслоение может существенным образом снизить несущую способность конструкции. Влияние нескольких видов расслоений на устойчивость стержней и оболочек было рассмотрено в работах [1, 2]. Эксперимен тально-теоретическое исследование устойчивости при внешнем давлении цилиндрических оболочек с кольцевыми и прямоугольными расслоени ями проведено в работе [3]. Теоретические результаты работы [3] по лучены на основе приближенной схемы, использующей на участке рас слоения эквивалентную жесткость.
Одним из неисследованных видов расслоений является продольное расслоение в цилиндрической оболочке, нагруженной внешним давле нием. Рассмотрим продольное расслоение вдоль всей длины оболочки (регулярное расслоение). Считаем, что при докритических нагрузках зона дефекта не распространяется. Межслоевая трещина (расслоение) разделяет всю оболочку на несколько независимых цилиндрических па нелей, определенным образом сопряженных. В силу симметрии рассмот
рим половину поперечного сечения оболочки |
(рис. 1). |
|
||
Уравнения |
нейтрального равновесия цилиндрических панелей I, II, |
|||
III имеют вид |
[4] |
|
|
|
|
TUA + T2U2 + N ,R -' = Q- |
1 |
+ |
|
|
МиЛ+М21,2- ^ = 0; |
/ |
чГ |
U |
N \11+ N2,2 “ Т22^?2_I + *7n —0.
Здесь (.. ,) { = д/д1й gb g2 — продольная и окружная координаты; R 2 =
= /?, R1= 00; R — радиус цилиндрической оболочки. Символ 1,2 обозна
чает, что невыписанные формулы получаются из приведенных переста новкой индексов 1 и 2. Подставляя значение нормальной составляющей внешней поверхностной нагрузки qn— q^2R (где к2 — изменение кри визны нейтральной поверхности в окружном направлении) в последнее уравнение системы (1), получим уравнения устойчивости оболочки при
действии внешнего давления q. |
|
|
|
|
В смежной форме равновесия связь между усилиями и моментами |
||
Tij, |
N u Mij, с одной стороны, и перемещениями щ, |
w, с другой, задается |
|
с помощью соотношений теории тонких ортотропных оболочек [4] |
|
||
^ 1 1 = ^ П [ ^ 1 , 1 -\~ w R i 1+ V2 1 (^2 , 2 + w R 2 !) ] i |
|
|
|
Т\2=*В\2[u2,\ + ^1,2+ "7Г h2Rz~l (-W, 12+^2,1^?2_1 + ^1,2^1_I) ] i |
(2) |
||
|
О |
|
|
11 |
~ ^11 [^,11 —^l.l^l-1 + V21 (w,22~ U2>2R 2~ l) ] ’> |
|
|
12=^ i2 (ДО,12~~ U2,\R2~X— U\i2R l_1) •
Здесь Bi^EihKl-VijVji); B ij=Gijh; Dij=— B ijh2; Eh Gijt v0- — модули

Интегрирование системы дифференциальных уравнений (7) с гра ничными условиями (8) проводим методом Кутта—Мерсоиа с промежу точным ортонормированием векторов решения [5, 6]. Значение крити ческой нагрузки находится из системы алгебраических уравнений (5). Алгоритм, реализующий программу поиска критической нагрузки, раз работан для БЭСМ-6.
Параметрический анализ проведем на примере круговой цилиндри ческой оболочки (jR/Am = 100) из композитного материала GI2/£' = 0,25 (Е\ = Е2=Е , vi2 = V2i =v = 0,25). Величину критического внешнего давле
ния для оболочки длиной I определим по формуле [7] |
|
|
#о= |
^ п^223(1—V12V21). |
(9) |
Исследуем влияние величины и места расположения межслоевой тре щины по толщине оболочки на значение критического внешнего дав ления.
Вначале рассмотрим расслоение, расположенное на нейтральной по верхности оболочки. Параметр h* = hI/hm i определяющий место распо ложения трещины по толщине оболочки, в этом случае равен А* = 0,5. На рис. 2 сплошными линиями обозначены зависимости критического параметра q*= qlq<n от параметра ширины трещины l*2 = h ! ^ R при раз личных значениях параметра /*= ///?, определяющего длину оболочки и расслоения. Здесь q01— величина критического внешнего давления обо лочки длиной /*= 1, определяемая по формуле (9); /2, 1\ — ширина и длина расслоения = Штриховыми линиями на рис. 2 обозначены критические нагрузки нереализующихся форм потери устойчивости.
Как видно из рис. 2, изменение длины оболочки существенно влияет на абсолютные значения критического внешнего давления. Однако об щий характер зависимостей q*= q*{l*2) при варьировании параметра длины Г сохраняется. Все эти зависимости при изменении /*2 от нуля до /*20 = 0,015—0,02 имеют горизонтальный участок, в пределах которого критический параметр q*= const. Следовательно, значение 1*2о опреде ляет предельную ширину допустимого расслоения. Меньшие по разме рам расслоения не снижают первоначального уровня критической на грузки q. Большие по размерам расслоения приводят к снижению критического внешнего давления и к локализации общей формы волно образования.
Общая потеря устойчивости происходит по одной из двух локальных форм: симметричной (рис. 3—а) и кососимметричной (рис. 3—б) (име ется в виду симметрия относительно центра расслоения). Изменение от носительных прогибов w* = w/w0 на рис. 3 показано для четверти попе речного сечения оболочки (w0 — наибольший прогиб). Исследование форм потери устойчивости показывает, что при симметричной форме мо-
Рис. 3. Общие формы потери устойчивости: |
а — симметричная (/*2= 0,12; /* = 2); |
б — кососиммстрнчная |
(/*2 = 0,05; /* = 2). |
Рис. 4. Зависимость критического параметра от размеров трещины, различно распо ложенной по толщине оболочки: /г*= 0,05 (/); 0,1 (2); 0,15 (5). /*=1.

ниже наименьшего значения критической нагрузки местной формы обо лочки с расслоением. В этом случае пересечение названных линий от сутствует, т. е. местная форма потери устойчивости не реализуется при любых размерах расслоений.
Для рассмотренных ранее оболочек длиной /* = 1 на рис. 6 представ лена зависимость параметра предельно допустимых расслоений /*20 = = l2o/2 nR от параметра Л*, характеризующего место расположения тре щины по толщине оболочки. Сплошная линия соответствует расслое ниям, приводящим к общей кососимметричной форме. Можно считать, что эта линия симметрична относительно вертикальной оси, проходящей через значение /г* = 0,5. Следовательно, расслоения, расположенные сим метрично относительно нейтральной поверхности оболочки в интервале 0,35—0,65, равноопасны. В общем случае на критический параметр и размер предельных расслоений влияет кривизна поверхности. По этой причине расслоение, расположенное ближе к внутренней поверхности оболочки, в большей степени снижает величину критического внешнего давления, чем расслоение, находящееся ближе к наружной поверхности. В соответствии с этим различаются и предельные размеры расслоений, обозначенные на рис. 6 штриховыми линиями. Эти размеры определены на основе смешанных форм потери устойчивости. Штрихпунктирные ли нии на рис. 6 соответствуют расслоениям, расположенным вблизи на ружной или внутренней поверхностей оболочки. Предельный размер та ких расслоений установлен исходя из возможной потери устойчивости по местной форме. Так как при местной форме следует ожидать неустой чивости только одного тонкого слоя, то потери несущей способности всей оболочки в этом случае не происходит (при условии нераспространения трещины). Следовательно, найденные на основе местной формы потери устойчивости предельные размеры расслоений (штрихпунктирные ли нии) не являются ограничениями по несущей способности.
Напомним, что все результаты получены при значении параметра толщины оболочки R/hJ1I= 100. Как показывают дополнительные рас четы, изменение этого параметра в широких пределах (R/hm = 30-^300) не меняет общего характера представленных зависимостей.
СП И С О К Л И Т Е Р А Т У Р Ы
1.Болотин В. В., Зебельян 3. X., Курзин А. А. Устойчивость сжатых элементов с дефектами типа расслоений. — Проблемы прочности, 1980, № 7, с. 3—8.
2.Трошин В. П. К устойчивости цилиндрических оболочек с расслоениями. — Ме
ханика композит, материалов, 1981, № 4, с. 729—731.
3. Андреев JI. В., Ободан Н. И. Задачи устойчивости цилиндрической оболочки с переменной жесткостью при внешнем давлении. — Прикл. механика, 1968, т. 4, вып. 12, с. 82—88.
4.Новожилов В. В. Теория тонких оболочек. Л., 1962. 430 с.
5.Годунов С. К. О численном решении краевых задач для систем обыкновенных
линейных дифференциальных уравнений. — Успехи мат. наук, 1961, т. 16, вып. 3(39),
с.171— 174.
6.Кармиишн А. В., Лясковец В. А., Мяченков В. И., Фролов А. Н. Статика и ди
намика тонкостенных оболочечных конструкций. М., 1975. 376 с.
7. Королев В. И. Слоистые анизотропные пластинки и оболочки из армированных пластмасс. М., 1965. 272 с.
Поступило в редакцию 28.09.81