книги / Теория наведенной неоднородности и ее приложения к проблеме устойчивости пластин и оболочек
..pdfБифуркация равновесных состояний при разупрочнении 171
Al(F>u S , r l \ r ^ 4 ^ u S ) + q + p ( § ^ , J + p ^ +k2] =
A j(F .u § .rW )= ° ;
Аз(р,и?.т:,г5)= °;
А^ ^ , г м И ьМ |
=0- |
|
(8.2.1) |
||
|
|
||||
Рассмотрим в качестве примера материал, |
обладающий свой |
||||
ствами трансверсальной изотропии: |
|
|
|||
Al(F.U-iY“7°2) = k |
, g +kjg |
+i^ l(> 4 c (x 3 A)(x? -^ )d xJ- |
|||
У гГ , |
Н |
, Н |
|
|
|
5 х ? |
d x id x l |
5 х ? 5 х 2 |
5 х з |
,a4uj| |
|
|
|
,Г^и?,|2 |
|||
|
|
L dxi |
d x id x i |
5x25X2 |
j |
. . к м ) . !
A3(u-l,YpY2)= ^ ( Ь ^ х з . т ^ х з - ^ ! dX3 *
hS
172 У стойчивость конструкций с учетом разупрочнения м атериал а
dxi |
ах?5x1 5x2J |
|
5х? 5x2 |
5x2 5xi |
5xi 5x25xi 5x2 |
Здесь полагаем также, что воздействие жидкой среды является |
||
равномерным по плану и симметричным |
относительно срединной |
поверхности оболочки, тогда оператор Ц и з ,и э ) получается заменой функции F на и? в выражении для L(uS»F); Е0модуль Юнга; Gu, G23модули сдвига материала в поперечном направлении.
Величина концентрации жидкости в материале с(х:>,т), содер жащаяся в интегралах по толщине, входящих в (8.2.2), определяется, как и ранее, по (8.1.2), то есть значением аргумента интеграла оши бок pz.
Исследование устойчивости оболочек состоит из двух этапов: определение исходного равновесного состояния и поиск точки раз ветвления исходного решения. Исходное состояние может быть оп ределено из решения общих нелинейных уравнений или каких-либо упрощенных вариантов, в частности, из линейных уравнений, полубезмоментных и других.
Упрощения в определении исходного состояния могут быть сделаны только при анализе конкретной решаемой задачи. Эти во просы для некоторых осесимметричных задач обсуждались в /58/.
Для поиска точек бифуркации решения используем линеари зованные уравнения устойчивости в вариациях, полученные путем дифференцирования по Фреше нелинейных операторов краевой за дачи, описывающей напряженно-деформированное состояние обо лочки под действием приложенных нагрузок и внешней среды.
Для получения уравнений устойчивости оболочек с низкой сдвиговой жесткостью, уравнения состояния которых построены на основе уточненной кинематической гипотезы, построим уравнения в вариациях:
Бифуркация равновесных состояний при разупрочнении |
173 |
/ a 4 j3 5 * 8 F ^ g W S F |
: ^SF эЧЛ ] |
Ч Зх ? 3x1 3x5 Эх? |
3xi 3X2ЭХ| 3X2/ |
,8F) - ( |
|
|
Л и ? 3*U°) 0 |
|
|
л Л и ? |
|
Л |
Эх? Эх? |
Эх? Эх? |
53x13Xi9х25xi 5xJ |
|
|
|
(8.2.3) |
В уравнения |
(8.2.3) не входят в явном виде у]1 и у°, но они |
||
входят неявно в и? |
и F через решение докритических уравнений. |
Заметим также, что, в отличие от (8.2.1), в (8.2.3) входят вариа ции искомых величин.
В дальнейшем будем считать, что все нагрузки, действующие на оболочку, изменяются пропорционально некоторому числовому
параметру, обозначаемому |
Р, то есть |
q (x i,x 2)=P*qs(xi,X2), |
Pi=P*Px(x2), P2=P*Py(xi), где |
q s, Рх, Pf заданные функции аргументов |
|
хь х2. |
|
|
Разрешая уравнения исходного состояния, определяем F и и® как функции от параметров Р, £> сф г): F = F(P,Cc(Pz)),
U§(P,Cc(Pz)).
Подставляя это решение в уравнения устойчивости (8.7), имеем нулевое решение. Ставится задача-найти такие значения параметра Р, при которых нарушается единственность, то есть, наряду с реше нием {U3°=0, F=0, у[‘=0, у2=0}, имеется хотя бы еще одно, отличное
от нулевого.
Рассмотрим шарнирно закрепленную, сжатую торцевыми уси лиями по образующей цилиндрическую панель из полимерного ма териала, подвергающегося разупрочнению под действием жидкой среды. Исходное состояние в этом случае является однородным и задача (8.2.3) превращается в линейную краевую задачу на собствен ные значения с переменными коэффициентами.
1 7 4 У стойчивость конструкций с учетом разуп рочнен ия м атериала
Переменность коэффициентов обусловлена их зависимостью от параметра концентрации жидкой среды в материале оболочки.
Для решения задачи используем приближенный метод сведе ния краевой задачи на собственные значения к алгебраической. В качестве такого метода примем метод Бубнова-Галеркина, представ ляя искомые функции в виде
F = XByФн(xi)0 2j(x2^ |
F =XAy03i(xi)04j(x2); |
|||
y? = S Bjj<D 5i(xi)a>6j(x2) |
Уз = £C ij<I>7i(xi)<I>8j(x2), |
|||
где {Фы}, (k=l,2,...,8) - |
системы координатных |
функций, удовле |
||
творяющих известным требованиям /59, 187/. |
|
|
||
В частном случае, |
приняв |
разложение |
в |
ряд по синусам и |
удерживая два первых (нечетный |
и четный) |
члена, удовлетворим |
граничным условиям, соответствующим шарнирному опиранию на гибкие несжимаемые (нерастяжимые) в касательной плоскости ребра, записываемым для полных функций.
Рассматривая только симметричное в плане оболочки дефор мирование, построим бифуркационную диаграмму деформирования при t=0 и t->°o. На рис.72 показан характер изменения кривых закритического деформирования Р(и§) и график критического време ни потери устойчивости в результате разупрочнения материала обо лочки под действием жидкой среды. Ось нагрузки на оболочку можно разбить на три характерных уровня нагружения (Рф - фикси
рованная во времени нагрузка): р* >pjp - исходная форма равнове
сия оболочки неустойчива при t=0; р*р<Рз<р"р - оболочка устойчи ва на интервале 0 < t < tkp; pj<p"p- оболочка устойчива на бесконеч
ном интервале времени. |
|
|
Форма потери устойчивости оболочки определяется |
значени |
|
ем безразмерного параметра |
кривизны кб (кб = ka2/(2h), |
а - размер |
оболочки в плане). |
_______ |
|
При кб=ккр (ккр=10я2л/О(с)В(с) , где D(c) и В(с) - жесткости |
оболочки на изгиб и растяжение-сжатие) имеет место случай крат ных значений критических нагрузок, соответствующих симмет ричной и несимметричной формам потери устойчивости.
Бифуркация равновесных состояний при разупрочнении |
175 |
При кб > ккр несимметричной форме потери устойчивости соот ветствует меньшее значение критической нагрузки. На рис.73 кри вые 1 и 2 соответствуют симметричной и несимметричной формам закритического деформирования оболочки. Значения Рс и Р„ (рис.73.) соответствуют нагрузкам потери устойчивости по симметричной и несимметричной формам, соответственно.
р
|
1 |
|
|
|
=2kkp *=0.3 1 |
||
|
л |
|
1 |
|
-} |
|
|
|
V - |
|
|
|
5 Чо |
л |
15 |
|
\ |
|
|
|
|
а___ |
|
Рис.72 |
Рис.73 |
|
|
Интерес представляет случай, когда кб<ккр. При этом мини |
|||
мальная критическая |
нагрузка соответствует симметричной форме |
||
потери устойчивости. |
Однако, как показывают расчеты, с течением |
времени в результате разупрочнения материала вследствие диффузии жидкой среды "опасной" формой потери устойчивости может стать несимметричная форма. Критическая нагрузка, ей соответствующая, становится минимальной. На рис.74,75 приведены графики критиче ских нагрузок, соответствующих симметричной форме потери ус тойчивости (сплошная кривая) и несимметричной форме потери устойчивости (штрихпунктирная кривая) в зависимости от парамет ра Pz (8.1.2).
На рис.75 графики Ркрфг) получены с использованием кинема тической гипотезы Кирхгофа-Лява. При значении кб = 0,8к«р в ре зультате процесса разупрочнения происходит смена "опасной" фор мы потери устойчивости с симметричной на несимметричную.
Рассмотрим влияние учета низкой сдвиговой жесткости мате риала оболочки на форму потери устойчивости. Сдвиговую жест кость оболочки будем характеризовать безразмерным параметром
176 У стойчивость конструкций с учетом разуп роч нен ия м атериала
a=0,5(l+v)a2Gi3/Eoh2, (Gi3 = G23). При а - х » деформации попереч ного сдвига отсутствуют, при а=30 деформация поперечного сдвига оказывает существенное влияние на устойчивость оболочки во
времени. На рис.76 приведены графики Ркр(Рг) для |
значения кб5 |
=0.5ккр. Как видно из приведенных результатов, учет |
деформаций |
поперечного сдвига позволяет существенно уточнить величину ми нимальной критической нагрузки и соответственно критическое время потери устойчивости оболочки с низкой сдвиговой жестко стью в условиях разупрочнения материала.
8.3. Выпучивание при взаимодействии вязко-упругого материала с жидкой средой
При рассмотрении вопроса об устойчивости конструкций из вязко-упругого материала выделим два этапа: этап нагружения мо нотонно изменяющимися силами и этап поведения конструкций во времени под действием постоянных сил.
Для тонкостенных конструкций с развивающейся неодно родностью материала на втором этапе проявляется влияние ухуд шения механических свойств материала, выражающееся в медлен ном движении конструкции (типа ползучести), которое называют выпучиванием.
Выпучивание при взаимодействии вязко-упругого материала |
177 |
Изучение выпучивания во времени, как возмущенного процес са, позволяет судить об устойчивости исходного состояния равнове сия. Процесс выпучивания, рассматриваемый как невозмущенное, основное состояние неидеальной системы, сам может являться объ ектом исследования на устойчивость.
Процесс выпучивания описывается уравнением:
f=0(f,N ,X i), |
(8.3.1) |
|
где f - характерное перемещение; |
N - параметр нагрузки; h - пара |
|
метры возмущения. |
|
|
При изучении устойчивости, |
как известно, определяется, какие |
состояния равновесия на мгновенной бифуркационной диаграмме, описывающей квазистатический процесс нагружения, устойчивы, а какие - неустойчивы в условиях выпучивания.
Рассмотрим приложение теории бифуркаций Пуанкаре к зада чам устойчивости полимерных вязко-упругих тонкостенных конст рукций, взаимодействующих с жидкой средой, приводящей к разу прочнению материала, на примере выпучивания нелинейной беско нечно длинной пластины.
Особенностью полимерных материалов является наличие вяз ких свойств, а также упругости, возникающей за счет раскручива ния цепных молекул /270/.
Для развития данного вида упругости требуется некоторый промежуток времени.
Описать поведение полимерного материала под нагрузкой можно с помощью схем-моделей, состоящих из системы упругих и вязкого элементов. В частности, это позволяет сделать модель стан
дартного вязко-упругого тела: |
|
cr+ Tp<T=Ee+TpHe, |
(8.3.2) |
где Тр - время релаксации; Н,Е - мгновенный и длительный модуль
упругости.
Уравнение в приращениях, описывающее поведение пластины во времени, имеет вид:
178 У стойчивость конструкций с учетом разупрочнения м атериала
(8.3.3)
здесь Uo - начальная погибь пластины; U, AU - прогиб пластины и его приращение; N, AN - продольная нагрузка и ее приращение; AM - приращение изгибающего момента:
Здесь, как и |
ранее, D*,D - изгибные жесткости в возмущенном |
и невозмущенном |
состоянии, a R*1, AR*1- кривизна и приращение |
кривизны, определяемые выражениями (8.1.4).
При изучении выпучивания с постоянной нагрузкой полагаем AN = 0; U = fsin(rcs/l); Uo = fosin(7ts/l), где 1- ширина пластины, и ис пользованы процедуры Бубнова-Галеркина и метода Рунге-Кутта.
Для решения вопроса об устойчивости состояний равновесия на мгновенной бифуркационной диаграмме рассмотрим поведение изо бражающей точки Р на фазовой прямой f вблизи положений равно весия 0(f,N,Xj) = 0.
На рис.77 приведена вспомогательная плоскость \|/=fl>(f,N,Xj) и показано поведение изображающей точки Р на фазовой прямой f.
Результаты получены при следующих значениях безразмер ных параметров: Н/Е = 2; f0/h = 10“3; (1/h)2 = 400; N6 = N /pJp, где h -
-толщина пластины, pjp - мгновенная Эйлерова критическая нагруз
ка.
Для значения Ne, равного половине от мгновенной критической нагрузки, существует одно положение равновесия, вблизи которого функция у меняет знак с плюса на минус, что говорит о его устойчи вости. При незначительном возрастании параметра нагрузки N6 = =0,505 на бифуркационной диаграмме длительного нагружения суще ствуют три положения равновесия, что отражается нулевыми значе ниями функции у на фазовой прямой.
Выпучивание при взаимодействии вязко-упругого материала 179
Рис. 77 Исследование процесса выпучивания пластины во времени
при Ne = 0,505 показывает, что состояние равновесия, характери зующееся малыми приращениями, является неустойчивым. Удален ные закритические состояния равновесия устойчивы, по Ляпунову, на бесконечном интервале времени.
Полагаем, как и ранее, что воздействие жидкой среды приво дит к изменению длительных упругих свойств полимерного мате риала, в частности Е(с). Если разупрочнение материала носит огра ниченный характер, то возможно решать задачу устойчивости на не ограниченном интервале времени. При двухстороннем воздействии среды расчетная схема сводится к трехслойной пластине со слоями переменной во времени жесткости.
Равновесие такой трехслойной пластины определяется не толь ко величиной нагрузки, но и уровнем разупрочнения материала, поэтому в качестве фазовых прямых примем характерное перемеще ние f и T *=xD z/h 2 - безразмерный параметр времени взаимодейст вия, где Dz - коэффициент диффузии жидкой среды в материал пла стины; h - толщина пластины.
В фазовом пространстве (\|/,£т*) построим бифуркационную поверхность длительного нагружения с учетом разупрочнения мате риала (рис.78). Степень разупрочнения материала, зависящая от времени взаимодействия т*, существенно влияет на форму бифурка ционной поверхности.
180 У стойчивость конструкций с учетом разуп роч нен ия м атериала
Исследование выпучивания пластины во времени t как возму щенного процесса позволяет определить устойчивые состояния рав новесия и критическое время взаимодействия материала и среды Ткр, по истечении которого происходит "обмен устойчивостью" и появление сжато-изогнутых состояний равновесия. При отсутствии агрессивного воздействия (т*=0) состояние равновесия (fsO , Ne =0,405) является асимптотически устойчивым. Исходное состояние равновесия остается устойчивым и в условиях воздействия среды для времени взаимодействия т*<т*к-р. При достижении временем взаи модействия т* значения т*кр происходит бифуркация состояния рав новесия. Устойчивыми становятся сжато-изогнутые состояния рав новесия, исходная равновесная форма за точкой бифуркации (т*=т*кр) является неустойчивой.