О численном решении задач устойчивости плоских колец

Тип работы:
Реферат
Предмет:
Механика


Узнать стоимость

Детальная информация о работе

Выдержка из работы

УЧЕНЫЕ ЗАПИСКИ ЦАГИ Т о м XX 1989
№ 3
УДК 629.7. 015.4. 023
О ЧИСЛЕННОМ РЕШЕНИИ ЗАДАЧ УСТОЙЧИВОСТИ ПЛОСКИХ КОЛЕЦ
Р. А. Мирошник
Приводится метод численного решения задач устойчивости для замкнутых плоских колец, нагруженных в своей плоскости произвольными распределенными и сосредоточенными нагрузками, фиксированными, либо следящими. Рассматривается потеря устойчивости кольца как с выходом из его плоскости, так и без выхода из нее.
Рассматриваются различные числовые примеры решения задач устойчивости для колец на упругом основании, нагруженных постоянным давлением.
Полученные результаты согласуются для частных случаев с известными в литературе
1. В статье приводится метод численного решения задач устойчивости для замкнутых колец, осевая линия которых до нагружения есть плоская кривая.
Рассматривается потеря устойчивости кольца как с выходом из его плоскости, так и без выхода из нее, при этом кольцо в своей плоскости может быть нагружено произвольными распределенными и сосредоточенными нагрузками, а также иметь переменные по длине жесткости. Кольцо может иметь произвольное число промежуточных упругих, либо жестких опор. Силы и моменты, нагружающие кольцо, могут быть как фиксированными, т. е. сохраняющими после потери устойчивости свои направления относительно неподвижных осей, так и следящими, т. е. сохраняющими свои направления в связанных осях. Кроме этого, внешние силы могут зависеть от перемещения точек осевой линии кольца и от их производных.
Среди известных методов исследования устойчивости такого рода систем [1, 4, 6] наиболее распространены энергетические методы, однако эти методы эффективны при условии задания с достаточной точностью форм потери устойчивости, что не всегда возможно.
Изложенный в статье метод связан с численным интегрированием линейных дифференциальных уравнений равновесия [3] с помощью метода Ньютона [2, 7].
При выводе уравнений равновесия кольца после потери устойчивости предполагается, что форма кольца в критическом состоянии мало отличается от его формы в обычном состоянии. Критическое состояние кольца — это его неустойчивое напряженно-деформированное состояние, предшествующее потере устойчивости, когда осевая линия кольца еще остается плоской кривой.
Из общих уравнений равновесия пространственно-криволинейного стержня [5] для частного случая плоского кольца получены дифференциальные уравнения равновесия кольца после потери устойчивости в безразмерной форме. В связанной системе
координат (е1г е2, е3), где в качестве независимой переменной взят угол ф, эти уравнения имеют следующий вид:
+ & lt-?!,-да.
Q2 = Ql* М3 + Qi — APt,
М$ - Qz Ат3, (1. 1)
(r)з ~ Мъ «2 = - и1& gt- и = ы2″
& lt-?3 = - & lt-?$ Af, + & lt-??>- м2 — Др3,
Afj= Af2 — М3, Л1а — Дт1″
М2= (9Af3* - 1) Afj -f Q3 — Д/Иг» 2)
«- = ра#, + «2,
«3 = - *2-
Здесь Qj, Mj, ft у и Uj (/=1, 2, 3)-соответственно безразмерные внутренние силы, моменты, углы поворота связанных осей и перемещения точек осевой линии после потери устойчивости, Api, Am. j — безразмерные приращения приложенных к кольцу нагрузок при его переходе в новое равновесное состояние. Величины со звездочкой Q$, О2*'-1 ^1**' @2*'- ^3*- компоненты безразмерных векторов внутренних
сил и моментов, соответствующих критическому состоянию кольца. Компоненты со звездочкой в рассматриваемом случае предполагаются известными из решения соответствующей линейной задачи статики для кольца.
Безразмерные величины связаны с размерными (с индексом р) зависимостями:
Qi = QipM2x2, — Q2pjAg х2, Q3 = Q3plA3 *2& gt-
Мх = MlpjA2 х, М2 = М2 Р1А3 х, М3 — Мзр1А3 х,
Uj=Ujp-%, APl = Дptр! А2 Др2 = bp2pIA2 х2, Ар3 = Ар3рАз х2, Ату = Атхр/Л3 х Ат2 = Ат2р1А3х, Ат3 = Am3pjA2x, р = A2/At
Oil* = Qi*"/4& gt- х2, Q$ = ОздМз xs, ' '-
Qf*1 = QuplA4 *3& gt- Qf* = Q2*plA2 x2,
где x = 1/p, p — радиус кривизны кольца, Л2, Л3 — крутильная и изгибные жесткости стержня.
Рассмотрим случай, когда приращения внешних нагрузок, входящие в уравнения (1. 1), не зависят от неизвестных функций, входящих в уравнения (1. 2) и наоборот приращения нагрузок из уравнений (1. 2) не зависят от функций уравнения (1. 1). В этом случае система плоско-пространственных уравнений равновесия кольца распадается на две независимые системы уравнений (1. 1) и (1. 2), причем первая из них служит для определения критических нагрузок без выхода кольца из плоскости, а вторая — при нахождении критических нагрузок с выходом кольца из плоскости.
Критические нагрузки для указанных случаев находятся из решения соответствующих краевых задач для неоднородной системы уравнений путем движения по параметру нагрузки [3]. Если исходная система дифференциальных уравнений однородна, она искусственно превращается в неоднородную путем приложения малой допол-
нительной нагрузки в направлении потери устойчивости так, чтобы решение исследуемой системы уравнений всегда было отличным от тривиального.
Для замкнутого кольца условия периодичности решений имеют вид соответственно для систем уравнений (1. 1) и (1. 2) —
& lt-?, (0) = & lt-?, (271), & lt-?2 (0) = & lt-?2 (2*), Мй (0) = М3 (2те) 93 (0) = «3 (2 л), и, (0) = щ (2 те), ы2 (0) = и2 (2 л).
(1. 3)
(?з (0) = & lt-?3 (2х), Л*! (0) = М, (2т.), Л12 (0) = М2 (2я), (0) = «і (2те), й2 (0) = (2те), и3 (0) = а3 (2 те).
(1. 4)


Ф1 = Ф. 8
Фи
Ф. і
Фір
Векторы невязок Фі и Ф2 соответственно для систем (1. 1) и (1. 2) равны
& lt-?і(0)-<-гі(2^
Фз (0) — Ол (2 я)
М3{0)-М3(2к), ф2= Ф23 = Лї2 (0) — Л12(2 те) (1−5)
«з (0) —з (2 я) н, (0) — Ц] (2 л) і И2 (0) — И2 (2 те) /
/Ф2Д / Сїз (0)-(?3(2г.)
ф22 Мх (0) — Мі (2 те)
Ф23 Л*2 (0) — М2 (2 те)
Ф24 «і(0)-"і (2 те)
і®*! «2 (0) — & amp-2 (2 те)
Ф2в/ '- «з (0) — и3 (2 те) у
Векторы начальных параметров К1 и К2, удовлетворяющие условиям периодичности (1. 3) или (1. 4) для исследуемых систем уравнений
Г,=
находятся из зависимостей
& lt-2і (0) С?2(0) М3(0)
«з (0) щ (0)
1 и2(0)
Угі'- / 0з (0)
у22 / АГж (0)
У 23 М2(0)
У'-и (0)
У 25 / МО)
Уж/ «з (0),
(1. 6)
У^УЇ+КгЧГг)-*!^), Г2=У° + ^2−1(1^)-Ф2(К°),
(1. 7)
где Кр Кд-произвольные векторы начальных параметров, задаваемые в нача-ле численного решения задачи, №Г1(К®), ИМ^г) — якобианы, вычисляемые численно с помощью разностных формул и равные
имк?) =
'- дуп '- ¦ ' йуі6
ду1Х '- '- '- д-у]6
, (Г°) =
ЙФ21 дФ2і
Зу2і дУ2в
д$ 2 В дФ2б
дЪх дУж
(1. 8)
Для формирования матриц (1. 8) каждая из систем уравнений (1. 1) и (1. 2) должна быть проинтегрирована семь раз.
Критические нагрузки находятся путем увеличения внешних нагрузок, от
которых линейно зависят компоненты С?2*'-& gt- ^1*)& gt- Фг**» ^з*& gt- от нуля Д° зна*
чений, при которых решения уравнений (1. 1) или (1. 2) резко возрастают [3].
2. В качестве примера рассмотрим задачу об устойчивости плоской формы кругового шпангоута, радиуса показанного на рис. 1. В рассматриваемом случае критические нагрузки находятся из системы уравнений (1. 2) с условиями периодичности (1. 4). Шпангоут, подкрепляющий оболочку 3, на которую Действует внешнее давление Рр, состоит из нетонкостенного свободного пояса 2, тонкой стенки 1 и связанного с оболочкой пояса 4. Пунктиром на рисунке показана форма потери устойчивости (опрокидывания) шпангоута. Свободный пояс 2 рассмотрим как кольцо, лежащее на упругом основании с известными линейной жесткостью кир и угловой жесткостью Наличие упругого основания приводит к появлению слагаемых — кииг и — ?& lt->-81
А-А
Рис. 1
в
50
100
150. П.
Рис. 2
соответственно в величинах Ар3 и Ami в уравнениях (1. 2), где ku = kUpR3/A3, kfr = k$ R/A3. Компоненты внутренних сцловых факторов в критическом состоянии, р
входящие в уравнения, равны
= -Р, = М3* = 0 (р = Рр ЩА3).
На рис. 2 показаны зависимости критического безразмерного давления р от упругости основания k (для стенки 1 постоянной толщины принято k = ku, ku= *& amp-/3).
Кривая 1 построена для следящей нагрузки рр (для различных значений |3 все кривые совпадают с ней). Кривые 2 и 3 соответствуют {5=0 и оо для фиксированной нагрузки, когда
Ap3 = - ku иъ — р.
Всем промежуточным значениям 0& lt-р<-°о соответствуют кривые, располагаю щиеся между кривыми 2 и 3.
Из расчетов получено, что при k = 0 значения критического давления совпадают с известными зависимостями [4] для круговых колец: ркр = 3 для следящей нагрузки и р"р=9/(4 + р) для мертвой нагрузки.
При вычислении критического давления для кольца без упругого основания (k = 0) для того, чтобы якобиан W2 (был невырожденным, необходимо исключить возможность смещения и поворота кольца как жесткого целого. Для этого надо вычислять критическое давление при малой жесткости упругого основания (/г & gt- 0), что практически не меняет искомую величину критического давления.
Для превращения исходной системы (1,2) в неоднородную, необходимо загрузить кольцо самоуравновешенной нагрузкой
& lt-?з = .Д sin 2 ?, (2. 1)
где, А — малая величина.
Предлагаемый алгоритм численного решения задач устойчивости позволяет находить критические нагрузки для плоских колец при произвольном законе изменения внутренних силовых факторов в критическом состоянии.
ЛИТЕРАТУРА
1. А л футов Н. А. Основы расчета на устойчивость упругих систем.- М.: Наука, 1978.
2. Демидович Б. П., Марон И. А. Основы вычислительной математики.- М.: Наука, 1960.
3. Мирошник Р. А. О численном решении задач устойчивости прямых стержней с упругими опорами. — Ученые записки ЦАГИ, 1988, т. 19 № 1.
4. Прочность, устойчивость колебания. Т. 1. — /Под редакцией Биргера И. А. и Паповко Я. Г. — М: 1968.
5. СветлицКий В. А. Нелинейные задачи статической устойчивости стержней. Расчеты на прочность. — М.: 1986. вып. 27.
6. Тимошенко С. П. Устойчивость упругих систем.- М.- 1946.
7. Шаманский В. Е. Методы решения нелинейных краевых задач на ЭЦВМ.- Киев.: Наукова думка, 1966.
Рукопись поступила 10/VI 1987

ПоказатьСвернуть
Заполнить форму текущей работой