CC BY
43
УДК 539.3:517.95:534.1
КРАЕВАЯ ЗАДАЧА ДЛЯ ГИБКИХ ОСЕСИММЕТРИЧНО НАГРУЖЕННЫХ СОСТАВНЫХ ОБОЛОЧЕК ВРАЩЕНИЯ И СТЕРЖНЕВЫХ СИСТЕМ
Э. Я. Еленицкий
Самарский государственный архитектурно-строительный университет,
443001, Россия, Самара, ул. Молодогвардейская, 194.
E-mail: elenit@list .ru
Предлагается эффективный алгоритм статического 'расчёта геометрически нелинейных составных тонкостенных конструкций. Используются линейные дифференциальные соотношения моментной теории. Нелинейность учитывается путём осуществления заранее не известного начального углового смещения каждого участка, сохраняющего форму образующей. Неизвестными разрешающей системы алгебраических уравнений являются произвольные постоянные общего решения и начальные углы поворота образующих. Линеаризация осуществляется итерационным методом Ньютона—Рафсона и обеспечивает высокую точность результатов.
Ключевые слова: геометрическая нелинейность, стержень, оболочка вращения, угловое смещение, дифференциальные уравнения, итерационный расчёт.
Введение. Большинство исследований конечных перемещений гибких конструкций основано на математических моделях, включающих системы нелинейных дифференциальных уравнений [1,2]. В настоящей статье предлагается альтернативный подход, позволяющий выполнять расчёты гибких плоских стержневых систем и осесимметрично нагруженных составных оболочек вращения без применения нелинейных дифференциальных зависимостей. Перемещения таких систем соизмеримы с размерами стержня (оболочки), а относительные деформации остаются малыми. Последнее обстоятельство указывает на то, что для каждого участка конструкции целесообразно выделить две группы перемещений, первая из которых соответствует его переносу без изменения формы образующей (нелинейная часть), а вторая — изгибу и продольному растяжению — сжатию (линейная часть). Первая группа перемещений определяется тремя степенями свободы: двумя линейными перемещениями и углом поворота начального сечения участка, что значительно упрощает последующий анализ напряжённо-деформированного состояния. Вторая группа перемещений и соответствующих им внутренних усилий описывается классическими линейными дифференциальными уравнениями и их интегралами [3]. Поскольку перемещения второй группы являются бесконечно малыми, возможна суперпозиция решений обоих этапов движения. Кинематические и статические условия сопряжения и внешнего закрепления границ участков обеспечивают формирование нелинейной системы алгебраических уравнений. При этом неизвестными задачи являются угловые смещения переносного движения, а также 6п (п — число участков) постоянных интегрирования общих решений дифференциальных уравнений. Каждый шаг построенного итерационного процесса включает решение линейной краевой задачи
Эдуард Яшевич Еленицкий (к.т.н., доц.), доцент, каф. сопротивления материалов и строительной механики.
для составной конструкции при угловых смещениях ее участков, полученных на предыдущей итерации, и уточнения величины угловых смещений следующего приближения. Сходимость нелинейного решения обусловлена тем, что процедуры линеаризации имеют простой геометрический смысл, и тем, что используются апробированные замкнутые линейные решения краевых задач. Ввиду того, что неизвестными задачи являются произвольные постоянные, отпадает необходимость построения матрицы жёсткости (податливости), как это предусмотрено в методе перемещений (сил) [4].
1. Построение решения. Рассмотрим элемент, вырезанный из системы, подверженной действию статических нагрузок. Будем относить понятие «элемент» к прямолинейному стержню на линейно упругом основании в условиях плоской деформации, а также к круговым коническим, цилиндрическим оболочкам или пластине при осесимметричном воздействии. Полагая перемещения большими, а деформации малыми, поставим задачу определения положения упругого элемента пк с учётом конечных перемещений, вызванных действием следящих внеузловых нагрузок р(ж), и условий закрепления на краях, допускающих действие силовых и кинематических факторов (рис. 1). При этом исходное положение элемента Поко длиной L в фиксированной глобальной системе отсчета XZ определяется декартовыми координатами и углом наклона а начального сечения. Прямая OZ является осью вращения, а отрезок щко — образующей пластины (a = 0), цилиндрической (a = тг/2) или конической (0 < a < тт/2) оболочки. Для стержня его исходное положение также совпадает с отрезком щко- Расчётное сечение элемента имеет координату ж, для которой
ж = хо + Ах, хо ^ х ^ L,
где для стержня и цилиндрической оболочки координата начального сечения Хо = 0, ДЛЯ круглой пластины И конической оболочки Хо = О\П0-
В рамках разработанного метода полагаем, что деформирование элемента происходит в два этапа: до некоторого положения noki и затем —до конечного положения пк. При этом угловое смещение на первом этапе не приводит к изменению формы образующей и определяется углом поворота <pi, получаемым в процессе расчёта таким образом, чтобы отрезки noh и пк оказались параллельными. Указанное смещение реализуется приложением к элементу трёхкомпонентного вектора дополнительных нагрузок рДж), имеющего такую же структуру, что и вектор р(ж). Здесь и в дальнейшем индекс г означает поворот — «rotation». Продольные и моментные нагрузки на рис. 1 условно не показаны.
Поскольку движение щко — noki сопровождается нелинейными перемещениями, используем систему отсчета XiZi, которая оказывается в конечном положении первой стадии деформирования. Введем подлежащие определению вектор-функции перемещений di и усилий fj:
di(x) = [щ(х) Wi(x) фг{х)]Т, fi(x) = [Ni(x) Qi(x) Mi{ж)]т, (1)
зависящие от выбранной системы координат следующим образом:
d0 = h(ipi)di, f0 = h(ipi) =
COS CPi — sin Lpi 0 sin Lpi COS СPi 0
0 0 1
(2)
Полагаем, что линейная часть соотношений (1) удовлетворяет решению дифференциальных уравнений равновесия моментной теории Кирхгофа—Лява:
Ьвіі(х) = р(ж), Їі(х) = (3)
В соотношениях (1)-(3) г/,і, иц, грі — продольные, поперечные и угловые перемещения расчётного сечения; Л^, Мі — соответствующие меридиональные усилия растяжения, поперечного сдвига и изгиба; Іі(ірі)—матрица преобразования векторов из подвижной в фиксированную XoZo (г =
= 0) систему координат; Ь, £/ —матрицы дифференциальных операторов соответственно для уравнений равновесия и физических соотношений между усилиями и перемещениями. Структура этих матриц имеет следующий вид [3]:
Ьп Ьі2 о" Ь/11 Ь/12 о’
ь = -^21 -^22 0 > Ь/ = 0 Ь/22 0
0 0 0 0 Ь/23 0
Движение элемента из исходного Поко в конечное положение пк характеризуется перемещением расчётного сечения по траектории гпогп^т^тл (рис. 2) Заменяем перемещение по дуге ШоШг движением вдоль отрезков гпогпг, Ш1Ш2, параллельных осям Zi, Xi. При этом обжатие участка погпо при его переходе в положение П0ГП1 компенсируется растяжением при последующем движении
Рис. 2. Траектория перемещения фиксированного сечения
до положения пот.2- В случае линейного решения перемещение ТОЧКИ то по траектории rno'mi не сопровождается обжатием участка погпо. Поэтому будем полагать, что последующее движение расчётного сечения до точки m2 происходит без увеличения продольных деформаций элемента, что компенсирует отмеченную выше потерю деформаций сжатия. Таким образом, перемещение расчётного сечения на участке riii'm2 характеризуется следующими зависимостями:
щ(х) = Дж(1 — cos <fi), Wi(x) = 0, Ni(x) = 0. (5)
Полные перемещения dj в координатах XjZj определяются векторами
dj = d ei + d ri + Adi, (6)
dri= [О Аж sin Lpi sin T , Adi= [Аж(1 - cos ipi) 0 (<^ - sint^)]T. (7)
Здесь сумма векторов d„ + Д<& — перемещения по траектории m^mim^
(этап 1), вектор dej — перемещения по траектории rn.2m.3m4 (этап 2). Первые два слагаемых в правой части равенства (6) представляют линейную часть решения, последнее учитывает конечность угла поворота ipi.
Перемещениям (6) соответствуют усилия, определяемые вторым соотношением (3), то есть
fj = Lf (dej + dri + Д<й) = £ei + fri + A/i) (8)
где с учётом зависимостей (4), (5), (7)
= Lf (Adi) = 0. (9)
Подстановка линейной части соотношений (6) в исходную систему уравнений равновесия (3) преобразует её к виду
Ldei = р - pri, pri = Ldri. (10)
Представим интегралы уравнений (10) следующим образом:
dei — dpi dri, fei — fl/jCj -j- fpi fri, (11)
где cidi — матрица размерностью 3x6 фундаментальных решений для элемента в положении noki; cifi — матрица той же размерности, получаемая с помощью второго равенства (3); Сi — шестикомпонентный вектор-столбец произвольных постоянных общего решения; dpi, dri, fpi, frj — частные интегралы, соответствующие нагрузкам р, pri.
Выражения для перемещений (6) и усилий (8) с учётом соотношений (9), (11) принимают окончательный вид:
di — 0-diCi dpi -\- di — (ifiCi fpi- (12)
Дальнейшие действия сводятся к формированию граничных условий и решению системы линейных алгебраических уравнений относительно неизвестных произвольных постоянных. Однако следует иметь в виду, что зависимости (2), (12) справедливы, когда линии noki и пк параллельны. Поскольку окончательное положение элемента заранее не известно, то угол ipi должен определяться на основе итерационного расчёта. Для очередной (г + 1)-й итерации имеем (см. рис. 1):
(fi+i = (fi + Aipi, tgAipi = (wki-wni)/L,
где wni, Wki — нормальные линейные перемещения концевых сечений элемента, полученные на г-том шаге расчёта.
Линеаризацию решения проводим методом Ньютона—Рафсона по схеме
APj = Рi — Ai Сi, Сг+1 = Ci + Ai 1APj, (13)
где Ai, Pi, APi — соответственно невырожденная матрица коэффициентов при неизвестных, столбец свободных членов и невязка разрешающей системы алгебраических уравнений для г-той итерации.
Процесс итераций по (13) начинается с линейного расчёта при ipi = 0 и производится до тех пор, пока величина Aipi не станет меньше заданной точности. Как известно, метод Ньютона—Рафсона обеспечивает быструю сходимость, а в случае решения задач с весьма сильной нелинейностью успешно сочетается с методом пошагового приращения нагрузки.
Окончательный результат для усилий после обеспечения сходимости итерационного процесса следует представить в виде
fi0 = h(tpi)~lfi,
где квадратная матрица в правой части равенства осуществляет преобразование полученных усилий fi на криволинейную ось элемента (см. рис. 1).
Поскольку приведённые процедуры основаны на замкнутых решениях дифференциальных уравнений, исключается влияние на результат дополнительных погрешностей, связанных с конечномерными аппроксимациями дифференциальных зависимостей. Если решение (2), (12) получено для элемента, длина которого не обеспечивает малость угла поворота dw/dx, следует произвести разбиение исходного интервала на некоторое количество участков и
сформировать систему алгебраических уравнений с учётом следующих условий сопряжения на внутренних границах участков:
Не—ldo;e—1 (Ье) — He&Q>e{XQe), Не— ifo;e— l(-t'e) — Не fo,e(^-Oe))
где do;e, fo,e — матрицы для элемента е, имеющие тот же смысл, что ив (2); Хое, Le — координаты крайних сечений элемента е; Не — матрица для перехода к глобальной системе координат XZ; е — 1, е — номера смежных элементов.
Замечаем, что наличие матриц Не позволяет также решать задачи для криволинейных стержней и оболочек вращения произвольной формы, образующие которых аппроксимируются ломаной линией с необходимым для заданной точности количеством участков. Кроме того, полученные соотношения допускают применение процедур формирования граничных условий не только для системы последовательно соединенных элементов, но и для составных конструкций ветвящегося типа [5].
2. Тестовые расчёты. Сравним результаты решений некоторых задач нелинейного деформирования гибких конструкций с аналогичными результатами, полученными на основе представленного метода углового смещения.
В качестве первого теста рассмотрим в рамках модели Бернулли консольный стержень постоянного сечения, изгибаемый моментом М (рис. 3, а). В линейной постановке кривизна определяется приближенно как вторая производная от функции прогиба, в нелинейной — точно, поскольку зависимость радиуса кривизны р от жёсткости стержня EJ приводит к следующему результату [6]:
1 / р = M/(EJ) = const. (14)
Из геометрических соображений с учётом равенства (14) имеем
u = x — psinip, w = р{1 — cos ф), Ф = х/р> (15)
где перемещения u, w, ф являются, как и ранее, компонентами вектора do.
Решение задачи на экстремум позволяет определить максимально возможный прогиб на конце стержня wmax и соответствующий этому перемещению момент Mq по формулам [6]:
Wmax = 0,7246L, М0 = 2,3311 EJ/L. (16)
При этом максимальный прогиб балки не зависит от её жёсткости, что является весьма удобным для последующего анализа.
На рис. 3, а приведены результаты решения по методике автора. Изображенные на нем перемещения показаны в масштабе размеров стержня, имеющего длину L = 10 м. Расчёты выполнялись при различном количестве итерационных циклов i с предварительным разбиением стержня на п = 5 равных участков. Перемещения стержня в линейной постановке (г = 0) представлены пунктирной кривой. Результат существенно изменяется после первой итерации (г = 1), а при г = 3 совпадает с точным решением (15), (16).
На рис. 3, б показаны графики зависимости относительной ошибки полученного решения е от количества участков разбиения п и количества итераций г. Следует отметить, что, несмотря на существенную нелинейность задачи, погрешность результатов не превышает 3 % уже при п = 3, г = 2.
Рис. 3. Изгиб консольного стержня
Рис. 4. Закритическое деформирование круглой пластины
Дальнейшее уточнение решения заданием п = 5, % = 2 снижает относительную ошибку до 1,5 %.
В качестве второго теста рассмотрена задача о закритическом поведении шарнирно опертой круглой пластины, находящейся под действием распределенных по контуру сжимающих радиальных усилий Р (рис. 4, а). В этом случае пластина способна после потери устойчивости нести возрастающую нагрузку, сохраняя при этом осесимметричную форму деформирования. Стратегия решения основана на том, что на первом шаге к пластине помимо сжимающих контурных усилий прикладывается небольшая по величине равномерно распределённая поперечная нагрузка. Последующий расчёт производится для конических элементов без участия этой нагрузки. На рис. 4, а, б представлены эпюры нормальных перемещений w и радиальных продольных усилий Nr, полученных для нагрузок Р различной интенсивности. В обозначениях использованы безразмерные параметры монографии [2]. Необходимо отметить, что независимо от метода расчёта наличие больших распорных усилий в сжато-изогнутой пластине существенно снижает сходимость итерационного процесса по сравнению со сходимостью «обычных» задач нелинейного деформирования. Приведенные графики получены при разбиении пластины по радиусу на п = 50 равных участков и количестве итераций i ^ 20. При этом погрешность для перемещений и усилий, сравниваемых с аналогичными величинами, полученными методом последовательных приближений по параметру [2], не превышает 1 %. На рис. 4, в представлены нормальные перемещений центра пластины в зависимости от различных значений п и Р. График, выделенный пунктирной линией, получен методом Бубнова—Га-лёркина с удержанием одного члена разложения [7,8]. Анализ результатов, полученных в центре пластины, показывает (рис. 4, г), что точность расчёта перемещений (кривая 1) и усилий (кривая 2) обеспечивается при различном количестве участков п.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Кармишин А. В., Лясковец В. А., Мяченков В. И., Фролов А.Н. Статика и динамика тонкостенных оболочечных конструкций. М.: Машиностроение, 1975. 376 с. [Karmishin А. V., Lyaskovets V.A., Myachenkov V.I., Frolov A.N. Statics and dynamics of thin-wall shell-type structures. Moscow: Mashinostroenie, 1975. 376 pp.]
2. Валишвили H. В. Методы расчёта оболочек вращения на ЭЦВМ. М.: Машиностроение, 1976. 279 с. [Valishvili N. V. Methods of analyzing shells of revolution on a digital computer: Mashinostroenie, 1976. 279 pp.]
3. Timoshenko S. P., Woinowsky-Krieger S. Theory of Plates and Shells. New York: McGraw-Hill, 1959. 595 pp.; русск. пер.: Тимошенко С. П., Войновский—Кригер С. Пластины и оболочки. М.: Наука, 1966. 636 с.
4. Смирнов А. Ф., Александров А. В., Лащеников Б. Я., Шапошников Н. Н. Строительная механика. Тонкостенные пространственные системы / ред. А. Ф. Смирнов. М.: Стройиз-дат, 1983. 488 с. [Smirnov A. F., Aleksandrov А. V., Lashchenikov В. Ya., Shaposhnikov N. N. Civil engineering mechanics. Thin-walled spatial systems / ed. A. F. Smirnov. Moscow: Stroyizdat, 1983. 488 pp.]
5. Еленицкий Э. Я. Расчёт свободных колебаний призматических систем с распределенными параметрами// Изв. вузов. Строительство, 1996. №7. С. 26-31. [Yelenitskiy Е. Ya. The calculation of free oscillations for prismatic systems with distributed parameters // Izv. Vuzov. Stroitel’stvo, 1996. no. 7. Pp. 26-31].
6. Лукаш П. А. Основы нелинейной строительной механики. М.: Стройиздат, 1978. 204 с.
[Lukash P. A. Basis of nonlinear civil engineering mechanics. Moscow: Stroyizdat, 1978. 204 pp.]
7. Вольмир А. С. Устойчивость деформируемых систем. М.: Наука, 1967. 984 с.
[Vol’mir A. S. Stability of Deformable Systems. Moscow: Nauka, 1967. 984 pp.]
8. Григолюк Э. И. К вопросу о поведении круглой пластины после потери устойчивости // Вестник инженеров и техников, 1949. №3. С. 103-106. [Grigoluk Е. I. On the question of circular plate behavior after buckling // Vestnik Inzhenerov i Tekhnikov, 1949. no. 3. Pp. ЮЗ-106].
Поступила в редакцию 02/V/2012; в окончательном варианте — 21 /XI/2012.
MSC: 35Q74; 74К20, 74К25
BOUNDARY VALUE PROBLEM FOR THE FLEXIBLE AXIALLY LOADED COMPOUND SHELLS OF REVOLUTION AND BEAMS SYSTEMS
E. Ya. Elenitsky
Samara State University of Architecture and Civil Engineering,
194, Molodogvardeyskaya St., Samara, 443001, Russia.
E-mail: elenit@list .ru
The effective algorithm of static calculation of geometrically nonlinear compound thin structures is offered. Linear differential equations of moment theory are used. Nonlinearity is considered by assigning unknown initial angular displacement of each segment retaining the form, of the generating line. Unknown values of algebraic equations resolving system are the arbitrary constants of the general solution and the initial angles of generating lines rotation. Linearization is realized by Newton-Raphson iterative method and provides the high precision of results.
Key words: geometrical nonlinearity, beam, shell of revolution, angular deflection, differential equations, iterative calculation.
Original article submitted 02/V/2012; revision submitted 21/XI/2012.
Eduard Ya. Elenitsky (Ph. D.(Techn.)), Associate Professor, Dept, of Materials Resistance and Civil Engineering Mechanics.
