Метод сплайнов в расчетах стержней на устойчивость Текст научной статьи по специальности «Физика»

ISSN 1992-6502 (P ri nt)_

2016. Т. 20, № 4 (74). С. 45-53

Ъьомт, QjrAQnQj

ISSN 2225-2789 (Online) http://journal.ugatu.ac.ru

УДК 539.374

Метод сплайнов в расчетах стержней на устойчивость

в. п.Павлов

[email protected] ФГБОУ ВО «Уфимский государственный авиационный технический университет» (УГАТУ)

Поступила в редакцию 09.10.2016

Аннотация. Рассматриваются возможности метода сплайн-функций степени 5 дефекта 1 для расчета критических сил при потери устойчивости прямого стержня, нагруженного продольной силой, при различных способах закрепления его концов. На примере решения ряда тестовых задач, имеющих точное аналитическое решение, анализируется точность предлагаемого метода. Показано, что реализованный алгоритм применения метода сплайн-функций степени 5 позволяет определять значения критических сил с весьма высокой точностью, ограничением которой являются только возможности современных компьютеров. Результаты расчетов по методу сплайнов сопоставляются с результатами расчетов критических сил в пакете ДМБУБ. Наблюдается ряд преимуществ по точности предлагаемого метода сплайнов в сравнении с методами, реализованными в ДМБУБ.

Ключевые слова: стержни, потеря устойчивости, критические силы, сплайны, ДМБУБ, математическое моделирование

ВВЕДЕНИЕ

В настоящее время широко распространено мнение, что приобретя современный вычислительный комплекс типа ЛК8У8, можно, не утруждая себя углубленным изучением механики деформируемого твердого тела, получать решения задач о прочности, жесткости, потери устойчивости и колебаниях разнообразных конструкций с очень высокой точностью.

При этом, естественно, на первое место выходит задача поиска необходимых финансовых ресурсов, так как все пакеты типа ЛК8У8 весьма дороги.

Конечно, ведущие конструкторские и исследовательские организации такие средства находят, и работники таких организаций начинают считать, что теперь у них не будет проблем с расчетами на прочность, жесткость и устойчивость конструкций любого вида. При этом они уверены, что необходимая точность расчетов уже гарантирована наличием сертификатов на применяемые вычислительные комплексы. В итоге большинство пользователей вычислительных комплексов им верят и не пытаются провести самостоятельную проверку точности применяемых программ на тестах, отличающихся от тестов создателей программ.

В нашей стране широкое распространение получил вычислительный комплекс ЛК8У8, который считается одним из самых мощных.

В данной статье, отказавшись от безусловного преклонения перед пакетом ЛК8У8, было проведено объективное исследования его точности при решении задачи об устойчивости прямого стержня в зависимости от числа применяемых конечных элементов.

Также для решения задачи устойчивости был реализован разрабатываемый автором метод сплайнов и оценена его точность в сравнении с точностью пакета ЛК8У8. Результаты исследований представлены в данной статье.

1. РЕЗУЛЬТАТЫ ОПРЕДЕЛЕНИЯ КРИТИЧЕСКИХ СИЛ ДЛЯ ПРЯМОГО СТЕРЖНЯ В ПАКЕТЕ А^УБ

В документации к Л№У8 приводится текст программы УМ 127, на основе которой дается оценка точности пакета ЛК8У8 при решении задачи об устойчивости прямого стержня постоянного поперечного сечения при защемлении нижнего конца стержня и приложении продольной критической силы к верхнему свободному концу стержня.

Исходные данные в задаче УМ127 заданы в американской системе единиц: длина стержня I = 100 т, поперечное сечение квадратное

b х h при b = h = 0,5 in , модуль упругости материала E =30-106 psi. В документации по ANSYS численное решение сопоставляется с известным точным решением, для которого при рассматриваемых условиях закрепления стержня

С1 (m)

значения критических сил г определяются формулой [1, 2]:

F(m)_

кр

22

221

(2m -1)2, m = 1,2,...:

(1)

где т - номер формы потери устойчивости.

Чтобы объективно подойти к вопросу оценки точности, выполним по программе VM127 расчеты критических сил для 1-ой, 2-ой, 3-ей,

5-ой и 10-ой форм потери устойчивости на сетках различной размерности при числе конечных элементов М = 10, 20, 50, 100, 200, 500,

1-103, 2-103, 5-103,1-104, 2-104, 5-104,1 105.

Точность расчетов оценим десятичным логарифмом абсолютного значения относительной

ошибки ^ | 5^кр |:

(Т) _ С Л / Т7(Т) I (2)

lg|ôFJ= lg | (FkP ) - Fkp )/FKPT )|

кР

кр / кр

где

F_ - расчетное значение критической

кр

силы при потере устойчивости стержня; ^^) -

точное значение критической силы при потере устойчивости стержня,

Результаты численных расчетов в пакете AN-SYS представлены для различных форм потери устойчивости на рис. 1 графиками зависимостей

логарифма относительной ошибки ^ | 5^кр | от

десятичного логарифма ^ М числа конечных

элементов М по длине стержня.

iglsFj

ПОГРЕШНОСТЬ ВЫЧИСЛЕНИЯ КРИТИЧЕСКИХ СИЛ В ПАКЕТЕ ANSYS

Форма потери устойчивости

О - m = 1

А. - m = 2

О - m = 3

■ - m = 5

v - m = 10

Квадратное сечение Защемление снизу, свободный сверху

3 4

Рис. 1.

igM

Анализируя результаты расчетов в пакете ANSYS (см. рис. 1), видим, что при числе конечных элементов в пределах от 10 до 100 наблюдается достаточно высокая точность расчета значений критических сил (в пределах 2-5 значащих цифр). А затем при росте числа конечных элементов погрешность расчетов начинает резко снижаться и при числе элементов 1-104 точность расчетов определяется всего 1-3 значащими цифрами.

2. ОСНОВНЫЕ СООТНОШЕНИЯ ПРИ РАСЧЕТЕ СТЕРЖНЕЙ НА УСТОЙЧИВОСТЬ МЕТОДОМ СПЛАЙНОВ

2.1. Дифференциальные соотношения, описывающие деформирование оси стержня при потере устойчивости

В работах [1, 2] изогнутая ось стержня при потере устойчивости описывается дифференциальным уравнением четвертого порядка:

d

2 С

dx2

2

EI,

d_w

dx2

+ F„

d2 w

dx2

= 0,

(3)

где w = w(x) - функции перемещения точек осевой линии вдоль оси Z ; E - модуль упругости материала; Iy - осевой момент инерции поперечного сечения стержня относительно оси Y ; FКр - сжимающая критическая сила, направленная вдоль оси стержня, совпадающей в недеформи-рованном состоянии с направлением оси X .

В случае постоянных E = const и

I = const уравнение (3) приводится к виду

d^W + F d^w = 0

y dx4 ^ dx2 ■

(4)

Уравнение (4) содержит производные до четвертого порядка включительно, поэтому при его решении целесообразно применять сплайны степени 5 дефекта 1, обеспечивающие непрерывность производных до четвертого порядка включительно. Методика построения и применения таких сплайнов для решения задач о статическом деформировании стержней изложена в работах [3-8].

Учитывая, что до настоящего времени сплайны пятой степени не применялись для анализа устойчивости стержней, в данной статье была решена задача оценки возможностей метода сплайнов степени 5 дефекта 1 для расчета критических сил сжатого прямого стержня. При этом в качестве эталона при оценке точности решения выбрана задача Эйлера об устойчивости стержней [1, 2].

2

5

6

2.2. Основные положения метода сплайнов степени 5 дефекта 1

При построении сплайна пятой степени [3] на отрезке [а, Ь] формируется сетка А :

А: а = <х2 <... <х^ = Ь, (5)

имеющая N узлов. На данной сетке строится сплайн-функция х) степени 5 дефекта 1,

имеющая N = N + 4 степеней свободы [3].

В пределах каждого отрезка [х1, хм], 1 = 1, ..., N -1 сплайн-функция

1 (х) является многочленом пятой степени

^5Д( х) = £ а«( х - х Г,

(6)

а=0

х е [х, х,+11 1 = 1, ..., N -1.

(10)

Сформируем из них векторы-столбцы:

V/ = (/ , 1 = 1,..., ^ , V/ = (/!, 1 = 1,..., N) т, V,/ = / , 1 = 1,..., Л)Т, V*/ = (Л11, 1 = 1,..., N)т,

V,4/ = (/Г, 1 = 1,..., N)т.

Векторы узловых значений сплайн-функции У^51( х) и ее производных определим, согласно [3], матричными выражениями:

V,

V

й/

V

й 2/

V

й 3/

V

й 4/

М / о,

-- Мй/ о

М й 2 / о,

М й 3 / о,

М й 4 / о,

(11)

Согласно [3], параметры, определяющие сплайн, сведены в вектор-столбец о из

N = N+4 параметров сплайна:

о = (дк, к = 1,2,..., N + 4)т,

(7)

где

Чх = ЩЛ хl), Чг =

Чз =

{Ч1+4 =

й^ЗдСх)

йх2

й 5,1 (х )

йх4

Ч4 =

й^5Д( х1) йх

й 31^д( х,)

йх3

(8)

/ = 1,2,..., N.

Число компонент вектора о равно числу степеней свободы сплайна.

Для узловых значений сплайн-функции

1/У5Д (х^) , а также производных до четвертого порядка включительно в узлах выбранной сетки А введем обозначения:

/г = х, ),

_й 2^5Л( х,) 1 йх2

йх 4

// =

щл хг)

йх

г=

йХ 1(хг) (9)

/г =

йх3 1 = 1, ..., N.

где М/, М/, Мй2/, Мйз/, Мй4/ - прямоугольные матрицы размера N х (N + 4) , формируемые по методике, изложенной в [3].

На основе соотношений (11) построим дискретный аналог дифференциального уравнения

(4) для вычисления критической силы —^ :

Е1уМй 4 / о + - Мй 2 / о = 0. (12)

М„,,, М ,

3. МОДЕЛЬНЫЕ ЗАДАЧИ, ИМЕЮЩИЕ ТОЧНОЕ РЕШЕНИЕ

3.1. Общее решение уравнения потери устойчивости стержня

Для оценки точности метода сплайнов при расчете критических сил рассматривается прямой стержень длиной I = 2 м , имеющий прямоугольное поперечное сечение шириной Ь = 0,03 м и высотой к = 0,02 м . Материал стержня - сталь с модулем упругости

Е = 2-1011 Па. Минимальный осевой момент инерции поперечного сечения стержня

1у = Ьк3/12 = 2 -10-8 м4.

Рассматриваются стержни с различными закреплениями по концам.

При традиционном обозначении [1]

-

а2 = —^ ЕТ,

<

<

уравнение (4) приводится к виду d4 w

■ + a

2 d2 w=0

(14)

dx4 dx2 Общий интеграл уравнения (14) имеет вид

w = C sin ax + C2 cos ax + C¿x + C4, (15)

где C1, C2, C3, C4 - произвольные постоянные, определяемые из краевых условий.

3.2. Стержень, шарнирно закрепленной по концам

Шарнирно закрепленный по концам стержень имеет краевые условия:

d2 w

w-

dx

■ 0 при x = 0 и x = l. (16)

Для краевых условий (16) на основе (15) получено:

с2 = 0, Q = 0, Q = 0, sin al = 0. (17)

Из последнего выражения (17) следует, что

(18)

тк

a =-, т = 1,2, ....

т i , , ,

При подстановки (18) в (13) получается выражение для определения точных значений критических сил, соответствующих потере устойчивости:

„ , т 2к2 EI v FT = a 2 EJ, = y

кр

m y

l2

m = 1,2,....

(19)

3.3. Стержень, защемленный на левом конце и шарнирно опертый на правом

Рассматривается стержень, защемленный на левом конце при X = 0 и шарнирно опертый на правом при X = I. Для такого стержня краевые условия имеют вид:

dw

w =-= 0 при x = 0,

dx

w = ■

d2 w dx2

(20)

= 0 при x = l.

При условиях (20) значения ат параметра а являются [1] решением трансцендентного урав-

нения:

al cos al — sin al = 0.

(21)

По найденным значениям ат рассчитываются точные значения критических сил:

^ = атЫу, т = 1,2,.... (22)

3.4. Стержень, защемленный по концам

Рассматривается стержень, защемленный на левом при X = 0 и правом при X = I концах. Для такого стержня краевых условия имеют вид:

dw

w =-= 0 при x = 0,

dx

dw

w =-= 0 при x = l.

dx

(23)

При условиях (23) значения ат параметра а являются [1] решением трансцендентного урав-

нения:

sin-

al 2

íal al — cos —

V 2 2

. al Л

sin

2

= 0.

(24)

При найденных значениях ат по формуле

(22) рассчитываются точные значения критических сил:

3.3. Стержень с жестким закреплением слева и свободным правым концом

Рассматриваемый стержень жестко закреплен на левом конце при X = 0 и имеет свободный правый конец при X = I. Для такого стержня краевые условия имеют вид [1, 2]:

dw w = — = 0 dy

d2 w d3 w 9 dw dy

dy2 dy3

■ + a

при y = 0,

= 0 при y = l.

(25)

При условиях (25) значения ат параметра а являются [1] решением уравнения

cos al = 0 и принимают значения:

к

ат =27 (2m — 1), т = 1,2,...

(26)

(27)

Подставляя (27) в (13), получаем выражения для вычисления критических сил:

F (m) = К 2 EJx

кр

2 ;2

221

(2m — 1)2, т = 1,2,

(28)

<

4. ДИСКРЕТНЫЙ АНАЛОГ УРАВНЕНИЯ ДЕФОРМИРОВАНИЯ ШАРНИРНО ЗАКРЕПЛЕННОГО ПО КОНЦАМ СТЕРЖНЯ

ПРИ ПОТЕРЕ УСТОЙЧИВОСТИ, И РЕЗУЛЬТАТЫ ЧИСЛЕННЫХ РАСЧЕТОВ КРИТИЧЕСКИХ СИЛ

Введем обозначение

х =

= —а =--

ЕТ,

(29)

на основе которого уравнение (12), являющееся дискретным аналогом дифференциального уравнения (4), принимает вид:

м, 4 / о -щ 2 / о=0.

(30)

Рассматривая стержень, шарнирно закрепленный по концам, запишем систему из четырех уравнений краевых условий

АО = 0. (31)

где А - прямоугольная матрица размера 4 х (N + 4) , заполняемая компонентами соответствующих строк матриц Му и М

й 2 / •

А, = К,

А, = Ч,, А, = К!,

(32)

, = 1,..., N + 4.

На концах шарнирно закрепленного стержня, согласно краевым условиям (16), равны нулю вторые производные

^ = 0 при X = 0 и X = I. (33)

Подставив (33) в (14), получим равные нулю четвертые производные

й4 w

dx

0 при X = 0 и X = I.

(34)

Запишем матричный аналог дифференциальных уравнений (34) в виде системы двух алгебраических уравнений

ТО = 0, (35)

где Т - прямоугольная матрица размера 2 х (N + 4) :

Т = (Т, I = 1,2, , = 1,2,...,N+4).

(36)

заполняемая компонентами соответствующих строк матрицы М,4у :

= М , >

\Т2, = МN4!, , = 1,..., N + 4.

(37)

Объединим системы уравнений (31), (35) в единую систему:

ио = 0 (38)

где и - прямоугольная матрица размера 6 х (N + 4)

и = (ии, 7 = 1,2,..., 6, , = 1,2, ..., N + 4) , (39) заполняемая компонентами матриц А и Т :

'К, = А,, 7 = 1, . ., 4, ] = 1,..., N + 4, К, = Т1,,, (40)

К, = Т2,,, , = 1,..., N + 4.

Из первых четырех столбцов и еще двух столбцов с номерами 5 и N + 4 матрицы и сформируем квадратную матрицу В размера 6 х 6:

В = (В,, 7 = 1,...,6 , = 1,...,6), (41)

компоненты которой определим соотношениями:

[К, = ии, 7 = 1,...,6, , = 1, ..., 4, В, = 7 = 1,...,6, (42)

{Вг,6 = и7,N+4 , 7 = 1,...,6.

Затем определим матрицу С, обратную к матрице В :

С = В-1 . (43)

Умножим слева матрицу и на матрицу С и получим матрицу и~

и = Си, (44)

особенностью которой является то, что в первых четырех столбцах, а также в столбцах с номерами 5 и N + 4 расположены только нули, за исключением того, что в каждом из этих столбцов расположено еще только по одной единице:

и1Д = 1, и 2,2 = 1, ихз = 1,

и 44 = 1, и 5,5 = 1, и 6, N +4 = 1.

<

Введем в рассмотрение вектор номеров ограничений, накладываемых на применяемый сплайн посредством матрицы и :

V, = V, к = 1,...,6),

(46)

$р1 V к'

где в рассматриваемом конкретном случае:

V = 1, У = 2, ^3 = 3, (47)

V = 4, V = 5, V = N + 4.

Используя матрицу и~ , обнулим в матрицах Мй2/, Мй4/, применяемых в матричном уравнении (30), первые четыре столбца, а также столбцы под номерами 5 и N + 4 . При этом используем формулы:

мйу = ы*у -Ум*2гик,,

ь 1 1 1ук к, 1'

к=1 6

мй= мй-умй4/ик,,

^ 1 ¡, 1 1ук к, 1'

(48)

к=1

1 = 1,..., N, 1 = 1,..., N + 4.

В преобразованных матрицах Мй2/, Мй4/ столбцы с номерами 1, 2, 3, 4, 5, N + 4 заполнены только нулями.

Сформируем на основе преобразованных

матриц Мй2/, Мй4/ квадратные матрицы

м

й 2/

, 4у. размера (N - 2) х (N - 2) :

Мй2/ = (м*2/, 1,1 = 1,...,N - 2),

= (А/?/

(49)

■й4/ = (м,, , 1,1 = 1,...,N - 2),

исключая из матриц Мй2/, Мй4/ нулевые столбцы с номерами 1, 2, 3, 4, 5, N + 4 и строки с номерами 1, N, воспользовавшись соотношениями:

М1 = +5, М1 = +5,

(50)

1 = 1, ..., N - 2, 1 = 1, ..., N - 2.

Введем в рассмотрение вектор-столбец параметров сплайна

о = к = 1,2,..., N - 2)т, (51)

компоненты которого определяются по компонентам вектора о :

& = Чк+5, к = 1,2,..., N - 2.

(52)

На основе Мй2/, Мй4/ и вектора о уравнение (30) принимает вид:

М й4/ о -Шй2 / о = 0 . (53)

Умножив уравнение (53) на матрицу ,

обратную к матрице , получим:

М-1 / й 4 / о ХЕо = 0, (54)

где Е - единичная матрица. Введем обозначения:

О = / М~й4 /, (55)

где О - квадратная матрица размера (N - 2) х (N - 2) .

При обозначении (55) уравнение (54) принимает вид

(О -ХЕ )о = 0. (56)

Система (56) имеет решения только при условии, что Х являются собственными числами матрицы О .

Определив значения Х, на основе (29) запишем формулу для вычисления критической силы потери устойчивости стержня:

-р =-ХЕТх. (57)

Точность расчетов оценим десятичным логарифмом абсолютного значения относительной

ошибки 1§ | 5—кр |, определяемым формулой (2).

Результаты численных расчетов методом сплайнов представлены на рис. 2 графиком зависимости 1§ | 5—кр | от десятичного логарифма 1§(м) числа отрезков сетки узлов сплайна

м = N-1.

-8

-2

ПОГРЕШНОСТЬ ВЫЧИСЛЕНИЯ КРИТИЧЕСКИХ СИЛ МЕТОДОМ КОЛЛОКАЦИЙ МСС®

4

Прямоугольное сечение Шарниры по концам стержня

Форма потери устойчивости

О - т = 1

▲ - т = 2

О - т = 3

■ - т = 4

V - т = 5

• - т = 6 А - т = 7

♦ - т = 8 □ - т = 9 ▼ - т = 10

34 Рис. 2.

|дм

2

5

Видно (см. рис. 2), что при шарнирном закреплении стержня по концам реализованный метод сплайнов характеризуется вторым порядком сходимости [2], и вплоть до числа отрезков между узлами сплайнов М = 16384 наблюдается непрерывное снижение погрешности расчетов с ростом М . При этом для М = 16384 при первой форме потери устойчивости (т = 1) модуль относительной погрешность расчета критической силы

| 8^кр |< 1 -10 8. Наблюдается (см. рис. 2), что с

ростом номера формы потери устойчивости погрешности расчетов возрастает, но и при десятой форме потери устойчивости (т = 10) модуль относительной погрешность расчета критической

силы также весьма мала | 5^ < 1-10 6.

5. РЕЗУЛЬТАТЫ ЧИСЛЕННЫХ РАСЧЕТОВ КРИТИЧЕСКИХ СИЛ ДЛЯ СТЕРЖНЕЙ С РАЗЛИЧНЫМИ СПОСОБАМИ ЗАКРЕПЛЕНИЯ ПО КОНЦАМ

5.1. Стержень, защемленный по концам

На рис. 3 в координатах ^ | 5Е |~ ^М

представлены зависимости ошибок численного расчета критических сил при первых десяти формах .для стержня, защемленного по концам.

|д|5^Кр|

ПОГРЕШНОСТЬ ВЫЧИСЛЕНИЯ

КРИТИЧЕСКИХ СИЛ МЕТОДОМ КОЛЛОКАЦИЙ МСК5

Прямоугольное сечение Защемление по концам стержня

Форма потери устойчивости

О - т = 1

▲ - т = 2

ф - т = 3

■ - т = 4

V - т = 5

• - т = 6

А - т = 7

« - т = 8

□ - т = 9

▼ - т = 10

34 Рис. 3.

|дм

Видно (см. рис. 3), что и для стержня, защемленного по концам, метод сплайнов характеризуется вторым порядком сходимости [2], и вплоть до числа отрезков между узлами сплайна М = 16384 наблюдается непрерывное снижение погрешности расчетов с ростом М . При этом при числе отрезков М = 16384 для первой формы потери устойчивости (т = 1) модуль относительной погрешность расчета критической

силы | 5^кр |» 1-10 8. Также видно (см. рис. 3), что с ростом номера формы потери устойчивости

погрешности расчетов возрастает, но и при десятой форме потери устойчивости (т = 10) погрешность расчета критической силы также весьма мала (| 5^кр |< 1 -10 6).

При этом отмечаем, что не было достигнуто пороговое значение величины М , при превышении которого начинается рост погрешности расчетов из-за величины погрешностей арифметических вычислений.

5.2. Стержень, защемленный на левом конце и шарнирно опертый на правом

На рис. 4 в координатах ^ | 5^кр |~ ^М

представлены зависимости ошибок численного расчета критических сил при первых десяти формах .для стержня, защемленного на левом конце и шарнирно опертого на правом.

|д|5^ф1

-4

ПОГРЕШНОСТЬ ВЫЧИСЛЕНИЯ Форма

КРИТИЧЕСКИХ СИЛ потери

МЕТОДОМ КОЛЛОКАЦИЙ МСК5 устойчивости О - т = 1

ж - т = 2

0 - т = 3

■ - т = 4

V - т = 5

• - т = 6

А - т = 7

4 - т = 8

□ - т = 9

Прямоугольное сечение ▼ - т = 10

Защемление слева, шарнир справа

3 4 Рис. 4.

|дм

Видно (см. рис. 4), что для стержня, защемленного слева и шарнирно закрепленного справа, имеет место второй порядок сходимости [2], и вплоть до М = 16384 наблюдается непрерывное снижение погрешности при росте М . При этом при числе отрезков М = 16384 для первой формы потери устойчивости (т = 1) модуль относительной погрешность расчета критической

силы | 5^кр |» 1-10 8. А при десятой форме потери устойчивости (т = 10) модуль относительной погрешность расчета критической силы

|5^|< 1-10"6.

5.3. Стержень, с жестким закреплением слева и свободным правым концом

На рис. 5 в координатах ^ | 5^кр |~ ^М

представлены зависимости ошибок численного расчета критических сил при первых десяти формах для стержня, защемленного на левом конце и сводного на правом на правом.

о

2

5

0

1

2

5

|д15^кр1

ПОГРЕШНОСТЬ ВЫЧИСЛЕНИЯ КРИТИЧЕСКИХ СИЛ МЕТОДОМ КОЛЛОКАЦИЙ МСК5

Прямоугольное сечение т -

Защемление слева, свободный справа

Форма потери устойчивости

о - т = 1

А - т = 2

О - т = 3

■ - т = 4

т = 5

т = 6

т = 7

т = 8

т = 9

т = 10

V -• -

А -♦ -

□ -

34

Рис. 5.

|дм

Видно (см. рис. 5), что для стержня, защемленного на левом конце и свободного на правом конце, имеет место второй порядок сходимости вплоть до момента, когда начинает проявляться влияние ошибки арифметических вычислений.

Процесс снижения точности расчетов из-за ошибки арифметических вычислений особенно заметен при вычислении критической силы при первой форме потери устойчивости для т = 1. Видно (см. рис. 5), что при м > 103 для т = 1 вначале снижается порядок сходимости, а затем при м > 4096 погрешности расчетов начинают расти.

В целом же (см. рис. 5) при м = 16384 для всех форм потери устойчивости погрешность

расчета критической силы | 5—кр |< 1-10 6 .

6. СОПОСТАВЛЕНИЕ РЕЗУЛЬТАТОВ РАСЧЕТА КРИТИЧЕСКИХ СИЛ В ПАКЕТЕ А^УБ И ПО МЕТОДУ СПЛАЙНОВ

Для стержня защемленного снизу при свободном верхнем конце погрешности расчетов в пакете Л№У8 представлены на рис. 1, а по методу сплайнов степени 5 дефекта 1 для такого же способа закрепления - на рис. 5.

Видно, что метод сплайнов демонстрирует явную сходимость к точному решению вплоть до м = 1 -103 а далее начинает сказываться влияние погрешностей арифметических вычислений. Реально достигнутая точность расчетов при м = 1 -104 лежит в пределах 6-8 значащих цифр. При этом следует отметить, что современные компьютеры выполняют расчеты с шестнадцатью значащими цифрами. Следовательно, при появлении компьютеров, обеспечивающих расчеты с большим числом значащих цифр, метод сплайнов позволит получать более точные результаты.

Рассматривая результаты расчетов в пакете ЛК8У8, видим, что при числе конечных элементов не выше 100 наблюдается достаточно высокая точность расчета значений критических сил (в пределах 2-5 значащих цифр). А затем при росте числа конечных элементов погрешность расчетов начинает резко снижаться и при числе элементов 1-104 точность расчетов определяется всего 1-3 значащими цифрами. А это очень низкая точность. Кроме этого, результаты, получаемые в пакете ЛК8У8 вообще не имеют каких-либо традиционных признака сходимости при увеличении числа конечных элементов

ЗАКЛЮЧЕНИЕ

В данной статье представлены результаты исследования возможностей метода сплайнов степени 5 дефекта 1 для расчета критических сил при потере устойчивости прямого стержня.

Проведена оценка точности предложенного варианта метода сплайнов при решений тестовых задач о потере устойчивости продольно сжатых прямых стержней с различными способами их закрепления по концам. Показано, что разработанная автором методика применения сплайнов обеспечивает второй порядок сходимости и точность расчета критических сил при различных формах потери устойчивости с 6.. .8 значащими цифрами на сетке с числом отрезком сплайна в пределах 104 .

Для сравнения возможностей метода сплайна с возможностями пакета ЛК8У8 выполнены расчеты критических сил в пакете ЛК8У8 и выяснено, что у программы ЛК8У8 есть определенный недостаток, заключающийся в том, что сростом числа конечных элементов резко возрастает погрешность расчета и уже при числе конечных элементов порядка 104 результаты расчета имеют точность лишь в пределах 1 %.

Из этого можно сделать вывод, что пакет ЛК8У8 возможно устроит пользователей при выполнении инженерных расчетов с точностью до 1 %, но может не устроить исследователей-механиков, имеющих намерение получать результаты с более высокой точностью.

В этом случае им предлагается обратить внимание на метод сплайнов, изложенный в данной статье.

По мнению автора, преимущества метода сплайнов по сравнению с традиционным методом конечных элементов заключается в том, что на его основе непосредственно решается дифференциальное уравнение, описывающее изучаемое явление потери устойчивости стержня.

0

2

5

СПИСОК ЛИТЕРАТУРЫ

1. Биргер И. А., Мавлютов Р. Р. Сопротивление материалов: Учебное пособие. М.: Наука, 1986. 560 с. [I. A. Birger., R. R. Mavlyutov. Mechanics of materials. Tutorial. M.: Science, 1986. 560 p. ]

2. Прочность, устойчивость, колебания. Справочник в 3 Т. / под ред. д-ра техн. наук И. А. Биргера и чл.-корр. АН Латвийской ССР Я. Г. Пановко. М.: Машиностроение, 1968. 568 с. [ Strength, stability, oscilations. Volume 3. Pod red. I. A. Birgera iA. G. Panovko. M.: Mashinostroenie, 1968. 568 p.]

3. Павлов В. П. Метод сплайнов и другие численные методы решения одномерных задач механики деформируемых твердых тел. Уфа: Уфимск. Гос. Авиац. техн. ун-т Уфа. 2003. 197 с. [V. P. Pavlov. Spline-methods and other numerical methods for solving one-dimensional problems of deformable solids mechanics. Ufa: USATU. 2003. 197 p. ]

4. Абдрахманова А. А., Павлов В. П. Математическое моделирование напряженно-деформированного состояния стеклопластикового стержня при различных жесткостях опор // Вестник УГАТУ, Уфа, 2007, Т.9, №5(23). С. 87-92. [ A. A. Abdrakhmanova, V. P. Pavlov. Mathematical modeling of the stress-strain state of glass-plastic beam with different stiffnesses of supports // Vestnik UGATU, vol. 9, No. 5(23), 2007, pp. 87-92. ]

5. Павлов В. П. Оценка возможностей метода сплайн-функций при решении уравнений, описывающих изгиб стержней // межвузовский научн. сб. «Механика деформируемых тел и конструкций». 1998. Уфа. С. 15-20. [V. P. Pavlov. Assessment the possibilities of method of spline-functions for solving equations describing the bending of beams. // Mezhvu-zovskij nauchn. sb. «Mekhanika deformiruemyh tel i kon-strukcij». Ufa, 1998, pp. 15-20.]/

6. Павлов В. П., Абдрахманова А. А., Абдрахманова Р. П. Задача расчета стержней одномерным сплайном пятой степени дефекта два // математические заметки ЯГУ. 2013. Т.20, вып.1, с. 50-59. [V. P. Pavlov, A. A. Abdrahmanova, R. P. Abdrahmanova. The problem of calculating beams by one-dimensional spline fifth degree of the defect two // Matematiches-kie zametki JaGU. vol.20, No.1, Jakutsk, 2013. pp.50-59. ]

7. Павлов В. П., Абдрахманова А. А. Расчет стержней на устойчивость методом сплайнов // материалы российской науч.-техн. конф. «Мавлютовские чтения». Уфа. Уфимск. гос. авиац. техн. ун-т, 2016. Т. 1, С. 123-127. [V. P. Pavlov, A. A. Abdrahmanova. Calculation of beams stability by sing spline-method // Mavlyutovskie chteniya: materialy Rossijskoj auchno-tekhnicheskoj konferencii. Ufa: UGATU, 2016, Vol. 3, pp. 23-127.].

8. Павлов В. П., Абдрахманова А. А. Расчет частот собственных колебаний стержней методом сплайнов // Материалы российской науч.-техн. конф. «Мавлютовские чтения». Уфа. Уфимск. гос. авиац. техн. ун-т, 2016. Т. 1, с. 128-132. [V. P. Pavlov, A. A. Abdrahmanova. The calculation of natural frequencies of the beam by splines // Mavlyutovskie chteniya: materialy Rossijskoj nauchno-tekhnicheskoj konferencii. Ufa: UGATU, 2016, Vol. 3, pp. 128-132. ]

ОБ АВТОРЕ

ПАВЛОВ Виктор Павлович проф. каф. сопр. мат. Дипл. инж. по авиац. двигателям (УАИ, 1973). Д-р техн. наук по динамике и прочности (УГАТУ, 2005). Иссл. в обл. механики композиционных материалов.

METADATA

Title: The Method Of Splines In The Calculation Of The Studs

For Stability Author: V. P. Pavlov Affiliation:

Ufa State Aviation Technical University (UGATU), Russia. Email: [email protected] Language: Russian.

Source: Vestnik UGATU (scientific journal of Ufa State Aviation Technical University), vol. 20, no. 4 (74), pp. 45-53, 2016. ISSN 2225-2789 (Online), ISSN 1992-6502 (Print). Abstract: Discusses the possibilities of the method of spline functions of degree 5 defect 1 for the calculation of critical forces during buckling of a straight rod loaded by a longitudinal force at a time, personal ways of fastening its ends. On the solution of several test problems with exact analytical solution, analyze the accuracy of the proposed method. It is shown that the implemented algorithm is the application of the method of spline functions of degree 5 allows to determine the values of critical forces with a very high accuracy, restriction where you are only the possibilities of modern computers. The results of calculations by the method of splines compared with the results of calculations of the critical-ical forces in ANSYS. There are a number of advantages in the accuracy of the proposed spline method in comparison with methods implemented in ANSYS. Key words: studs, buckling, critical force, splines, ANSYS, mathematical modeling

About author:

PAVLOV, Victor Pavlovich, Prof., Dept. of Strength of Materials. Dipl. Engineer for Aircraft Engines (Ufa Aviation Inst., 1973). Dr. of Tech. Sci. (UGATU, 2005).