CC BY
88
Таблица 5
Возможные магнитные структуры для упорядоченного №5Се04(В03)2 (пространственная группа Рп21т (№ 31))
Т Направления магнитных моментов на ионах
1 Є1х + Є2Х - Є3Х - Є4Х + Є5х - e6x - e13x - e14x + e15x + e16x - eWx + e1^ e1y + e2y + e3y + e4y + e5y + e6y - e13y - e14y - e15y - e16y - e17y - e18y; e1;, + e2z - e3z - e4z + e5z - e6z + e7z + e8z + e9z - e12z + e13z + e14z - e15z - e16z + e17z - e18z - e19z - e20z + e22z + e23z
2 e1x + e2x - e3x - e4x + e5x - e6x + e'x + e8x + e9x - e12x + e13x + e14x - e15x - e16x + e17x - e18x - e19x - e20x + e22x + e23x; e1y + e2y + e3y + e4y + e5y + e6y + e7y + e8y + e9y + e12y + e13y + e14y + e15y + e16y + e17y + e18y + e19y + e20y + e22y + e23y; e1. + e2z - e3z - e4z + e5z + e6z - e13z - el4l + e15z + e16z - eiyz + e^
3 e1x + e2x + e3x + e4x + e5x + e6x + e7x + e8x + e9x + e12x + e13x + e14x + e15x + e16x + e17x + e18x + e19x + e20x + e22x + e23x є1. + e2y - e3y - e4y + e5y + e6y + e7y + e8y + e9y - e12y + e13y + e14y - e15y - e16y + e17y - e18y - e19y - e20y + e22y + e23y; e1;, + e2z + e3z + e4z + e5z + e6z - e13z - e14z - e'5z - e‘6z - ez - ez
4 e1 x + e2 x + e3 x + e4 x + e5 x + e6 x - e13 x - e14 x - e15 x - e16 x - e1' x - e18 xi є1. + e2y + e3y + e4y + e5y + e6y - e13y - e14y - e15y - e16y - e17y - e18y e1z + e2z + e3z + e4z + e5z + e6z + e7z + e8z + e9z + e12z + e13z + e14z + e15z + e16z + e17z + e18z + e19z + e20 + e22z + e23z
Таким образом, в ходе проделанной работы с помощью симметрийного анализа были построены возможные магнитные структуры для полностью неупорядоченного и различных случаев упорядочения магнитных ионов кристалла №5ве04(В03)2. Как показывают результаты рентгеновских исследований кристаллической структуры №5ве04(В03)2, упорядочения N1 и ве в позиции 4g даже после длительного отжига не наблюдается, магнитные ионы в позиции 4g располагаются случайным образом. Несмотря на то, что то N1 и ве распределены равновероятно, возможно образование локального (ближнего) порядка в небольших областях кристалла. В этом случае различные области могут иметь несколько различную магнитную структуру. Как показали результаты расчета, если упорядочение происходит в рамках одной кристаллографической ячейки, то направления магнитных моментов либо лежат в плоскости аЬ, либо направлены вдоль оси с. В случае удвоенной ячейки различные типы упорядочения магнитных ионов дают более сложную картину: для некоторых представлений в этом случае магнитный момент для некоторых ионов имеет все три компоненты. Кроме этого, в ряде случаев при единственном фазовом переходе магнитоупорядоченная структу-
ра может быть только ферромагнитной или наоборот только антиферромагнитной. Вероятнее всего реальная магнитная структура исследуемого соединения довольно сложная и может зависеть от условий роста кристалла (поскольку может изменяться распределение магнитных ионов в позиции 4g). Полученные в данной работе результаты помогут при анализе и интерпретации экспериментальных данных.
Библиографические ссылки
1. Ковалев О. В. Неприводимые представления пространственных групп : науч. изд. Киев : Изд-во АН УССР,1961.
2. Изюмов Ю. А., Найш В. Е., Озеров Р. П. Нейтронография магнетиков : науч. изд. М. : Атомиздат, 1981.
3. Толедано Ж.-К.; Толедано П. Теория Ландау фазовых переходов : науч. изд. М. : Мир. 1994.
4. International Tables for Crystallography (2006). Vol. A. : sci. ed. / edit. Theo Hahn - Dordrecht (Netherlands): Springer, 2005.
© Назаренко И. И., Софронова С. Н., 2013
УДК 519.62
МОДАЛЬНЫЙ КОНЕЧНО-ЭЛЕМЕНТНЫЙ АНАЛИЗ БАЛОК,
ПОДАТЛИВЫХ ПРИ ТРАНСВЕРСАЛЬНОМ СДВИГЕ
В. А. Нестеров
Сибирский государственный аэрокосмический университет имени академика М. Ф. Решетнева Россия, 660014, г. Красноярск, просп. им. газ. «Красноярский рабочий», 31. Е-шаіі; gorkolya@mail.ru
Рассматривается задача о собственных колебаниях балки с низкой трансверсальной сдвиговой жесткостью. На основе теории пластин Рейсснера-Мидлина разработана конечно-элементная модель, в которой в качестве узловых неизвестных фигурируют углы трансверсального сдвига. Представлены результаты модальных расчетов изотропных, трехслойных и композитных балок при учете неклассических граничных условий.
Ключевые слова: балка, собственные колебания, трансверсальный сдвиг, метод конечных элементов.
MODAL FINITE ELEMENT ANALYSIS OF SHEAR FLEXIBLE BEAMS
V. A. Nesterov
Siberian State Aerospace University named after academician M. F. Reshetnev 31“Krasnoyarskiy Rabochiy” prospect, Krasnoyarsk, 660014, Russia. E-mail: gorkolya@mail.ru
The author considers a problem of free oscillations of a beam with low value of transverse shear stiffness. On the base of Reissner-Midlin plate theory, the finite element of a shear flexible beam is developed. There are transverse shear strains among the basic kinematic nodal variables. The urgency of the model is shown by test modal analysis of beams. Vibrations of sandwich beams and thick layered composite beams under non-classical boundary conditions are investigated.
Keywords: beam, free vibrations, transverse shear strains, finite element method.
Композиционные материалы часто используются в производстве ракетно-космической техники, поскольку они, обладая высокими удельными механическими свойствами, позволяют изготавливать конструкции с большой степенью весового совершенства. Для композитов характерна низкая трансверсальная сдвиговая жесткость. Эту особенность поведения композитов необходимо учитывать при проектировании и расчете. Математическая модель при этом усложнится по сравнению с той, что используется при расчете конструкций из традиционных материалов.
Рассматривается задача о собственных колебаниях балок, податливых при трансверсальном сдвиге. Решение выполняется с помощью метода конечных элементов (МКЭ). Математическая модель базируется на теории пластин Рейсснера-Мидлина [1; 2]. Для вывода матриц жесткости и инерции используется вариационный подход. Для этого выражение полной энергии конечного элемента колеблющейся балки варьируется по узловым компонентам, среди которых присутствуют осредненные по высоте сечения углы трансверсального сдвига. Разработанная конечноэлементная модель протестирована на примерах решения классических задач. Ее актуальность показана расчетами, выполненными для монолитных изотропных балок с низкой трансверсальной сдвиговой жесткостью, для трехслойных балок с податливым заполнителем, а также для балок, изготовленных из однонаправленного композита (ортотропная модель). Частоты и формы собственных колебаний подтверждаются конечно-элементными решениями в среде пакета С08М08/М, где эти задачи рассматривались в более строгой (двумерной) постановке.
Исследуемая конечно-элементная модель допускает описание граничных условий неклассического вида, когда на торцах балки регламентируется поведение углов трансверсального сдвига. Показано, что для некоторых типов граничных условий задача о собственных колебаниях балки не может быть решена при использовании классических расчетных моделей.
Получим разрешающие уравнения для задачи о свободных поперечных колебаниях балки с низкой трансверсальной сдвиговой жесткостью. Задачу будем решать с помощью МКЭ, поэтому сначала необходимо записать выражение энергетического функционала. Полная энергия колеблющейся балки скла-
дывается из потенциальной энергии деформации и кинетической энергии движения:
E = О + Г . (1)
Здесь
'А + ) dxdydz ,
0 = 21 1 1(х + &*
(V)
Т=-2111р(х, у, Х) V2 йхйуйх , (2)
2 (V)
где р - плотность материала в точке с координатами (х, у, Х); V - скорость движения этой точки тела; V - объем тела.
Последовательно применяя гипотезы о поведении балки с низкой трансверсальной сдвиговой жесткостью [3], вместо (1) получим
1 \ Глг Ои О9 ^
^ I Ы~Г + М~Т + Оч I йх -
20 ^ ах ах )
1 к(«2 + »2)+(*-°Х)+ор[ч-°х) }
х dx = 0 .
(3)
где Ы, М, О - внутренние балочные усилия (Ы - продольная сила, М - изгибающий момент, О - перерезывающая сила), вычисляемые по формулам
N = B- + С
dx
С ,1 ,12 \
d у d w
dx dx2
M = С^ + D
dx
2
d у d w
dx dx2
Q = K у,
(4)
где Bp, Ср, Dp - параметры, характеризующие инерционные свойства
ь_
2 h-s
b_
2 h-s
Bp = J J pdydz, Ср = J J pzdydz,
_b -s 2
b 2 h-s
Dp = J J P z2 dydz.
(5)
_ b -s 2
Здесь И и Ь - высота и ширина сечения балки соответственно.
Перейдем теперь непосредственно к формированию матрицы жесткости балочного конечного элемента. Будем полагать, что прогиб изменяется по кубическому закону вдоль длины элемента, а угол сдвига у и продольное перемещение и - по линейному закону:
^ (х) = а1 +а2х + а3х2 +а4х3, у (х) = а5 + а6х ,
и (х) = а7 +а8 х. (6)
Введем в рассмотрение вектор основных кинематических переменных
о I ,
о = { w — у и ^ . ёх
Сх
ёх
26
^36
77
В В В В
=12 7 , к12 = 6 “2 , !2 к15 =-12 у, к16=6 і2
В В С В
= 4— і , к23 = , 23 і 24 = і ’ к25 = -6 Ц2
В В СВ кі 1 С
= 2— і , к27 = "У > к28 = —, к33 = —+ ■ і 33 і , к34 — 3 34 і
= В В кі С В
= і ’ л^37 — + 37 і 6 , к38 =і Л к44 = і
С С В В
=7 ’ к47 = - 7 ’ к48 =- 7 ’ к55 = 12 ~3
В В В С
= -6 ? ^ к66 =4 7, к67 = - 7 ’ к68 =- Т
= В = і 13Вр кі С 3 , к78 = і > к88 і2 + 42Вр = В = і ’ 11Вр і2 Вр Вр
р р 1 о • _ р , р _ р
35і
210 10
Ср
Р 3(3^р 12 -28Др) _ 13В>р і2 Вр
514 = 2 , ^15 = , ^ = + "
43
2
Вр
2
70і
Вр 5 = Ср
2 , 18 2
46
420 10
517 = — , 517 =—р , 518 = —- , 522 = —^Вр 1 +—Вр і
17 „ > 17 „ > 18 „ > 22 Ю5 р !5 р
5*23 =
524 =
Ее!
12
1 3 1
526 =----------Вр і3----------Вр і
26 140 р 30 р
СВ
12
527 =
Вр 1
533 =
537 =
Вр 1
534 =
Ср 1
12
Вр
13Вр і2 420
528 =
535 = "
536 =
(7)
Вектор узловых кинематических параметров произвольного элемента балки представим в следующем виде:
8 = {wl Г ~г} У1 и\ W2 í У 2 и21 . (8)
66
577 =
.Вві 6 ’ Ср і 538 =~Г, Вр 1 544 = ~^ 44 3 , 545
Ср і Ср і Вр і 13Вр і2
—— , 547 1^ 47 = , , 548 6 = 6 , 555 3:
11Вр і2 - Вр - = 7 5
210 10 , 558 :
Вр і3 105 р +—Вр і , 15 р = Вр1 , 568 =
567 12
_Вр
10
Ср!
12
Вр 1
Ср
Ср
С1
12
Вр!
3
578 =
Ср!
3
588 =
Вр!
3
Следуя вариационному подходу к МКЭ, получим систему уравнений для произвольного балочного элемента, которая в матричном представлении имеет следующий вид:
К, 6, - ш2 Б, 6, = 0, (9)
где К* и 8е - матрица жесткости и матрица инерци-альных параметров балочного элемента балочного элемента со следующими ненулевыми компонентами:
С помощью матриц К* и Бе формируются глобальные матрицы жесткости и инерции Ке и Бе, фигурирующие в глобальных уравнениях состояния
Ке А - ш2 А = 0, (10)
где глобальный вектор узловых неизвестных имеет вид
А = {8х 62 ... 6; ... 6я+х}7, б; - вектор узловых неизвестных /-го узла
(11)
8,- = {w¡
dw
Сх
Уі и
Описанный алгоритм вывода матрицы жесткости и матрицы инерции реализован в специальной программе формирования фундаментальных матриц и векторов теории МКЭ, предназначенных для расчета балок, податливых при трансверсальном сдвиге [4].
Выполним тестовое решение для стальной (Е = 210 ГПа, д = 0,3, р = 7 800 кг/м3) балки длиной 1 м, с прямоугольным поперечным сечением (И = 1 см, Ь = 1 мм). Будем полагать, что координатная ось X проходит через центры тяжести сечений. В этом случае смешанные жесткости (С), фигурирующие в матрице жесткости элемента К,, ,и параметр инерции (Ср), присутствующий в матрице Б,, равны нулю.
Сначала в нашем КЭ расчете зададим граничные условия свободного опирания, причем углы сдвига у на торцах не фиксируются:
при х = 0: w = 0, и = 0, при х = Ь: w = 0. (12)
Результат сравним с классическим решением, определяемым по формуле
(13)
где 3 - момент инерции сечения.
С помощью программы [5] решена обобщенная задача на собственные значения для 50-узловой моде-
3
3
2
2
2
5лс =
46
56
2
г
22
46
12
2
ли стальной балки длиной 1 м с прямоугольным поперечным сечением (h = 1 см, b = 1 мм) и получено значение первой круговой частоты, равное 147,811 (рад/с). Вычисленное по формуле (13) значение равно 147,833 (рад/с). Форма колебаний при этой частоте показана на рисунке.
Представлены первые пять частот собственных колебаний стальной балки, определенные по нашей модели, вычисленные по формуле (13) и найденные в результате решения в пакете COSMOS/M с использованием элементов PLANE2D (табл. 1). Почти идеальное совпадение по частотам подтверждает правильность рассматриваемого здесь подхода. Совпадение имеет место и по формам колебаний. Отметим, что результаты модального анализа в пакете COSMOS/M выдаются в виде циклической частоты f измеряемой в Гц (с-1) и связанной с круговой частотой формулой ю = 2 п f
Первая форма колебаний балки
Таблица 1
Частоты собственных колебаний для первых пяти мод
Важно отметить, что в данном случае частоты первых мод колебаний являются первыми в списке собственных значений задачи, т. е. наименьшими.
Выполним модальный анализ, уменьшив в 100 раз (по сравнению с предыдущим расчетом) жесткость балки при трансверсальном сдвиге. Все остальные параметры жесткости, размеры и граничные условия оставим без изменения. Результаты решения по нашей модели представлены списком первых пяти собственных значений задачи, расположенных во втором столбце (табл. 2). Здесь же для сравнения приводится
решение (табл. 2, третий столбец), выполненное в пакете COSMOS/M, в котором для учета уменьшенного параметра трансверсальной жесткости задействуется опция ортотропных свойств материала (для элементов PLANE2D). При этом фактически уменьшается соответствующий модуль сдвига в списке механических свойств материала. Вновь отмечаем хорошее совпадение результатов. Подчеркнем необходимость учета трансверсального сдвига в модальном анализе. Сравнение соответствующих данных (см. табл. 1 и 2) показывает заметное расхождение частот собственных колебаний, определенных с учетом трансверсального сдвига, и теми же величинами, вычисленными по классической формуле, где этот эффект не принимается во внимание.
Выполнено решение задачи о собственных колебаниях балки с шарнирно-опертыми концами, на которых зафиксирован угол трансверсального сдвига. В этом случае рассматривается разрешающая система алгебраических уравнений со следующими граничными условиями:
при х = 0: w = 0, и = 0, у = 0,
при х = L: w = 0, у = 0. (14)
Результаты решения (первые пять частот) приведены в четвертом столбце (табл. 2). Они показывают, что для рассматриваемой здесь шарнирно-опертой балки (с данным набором геометрических и механических параметров) расхождение в собственных значениях, определенных по двум моделям (с фиксацией угла у на торцах и без такового), невелико. Тем не менее качественное расхождение - верное: фиксация угла у на торцах делает балку более жесткой, что отражается в повышении величин собственных частот.
Таблица 2
Частоты колебания шарнирно-опертой балки
Номер моды Торцы со свободным сдвигом Решение в пакете COSMOS/M Торцы с фиксированным сдвигом
1 146,271 146,058 146,308
2 567,477 564,397 568,005
3 1 218,091 1 204,820 1 220,416
4 2 040,959 2 006,992 2 047,117
5 2 982,133 2 915,593 2 994,430
Параметры инерции балки в данном случае для прямоугольного сечения с шириной Ь и высотой И определяются следующим образом:
Ь И3
ВР=РЬИ, ср = ^ £р=р —.
Рассмотрим трехслойную балку с податливым заполнителем. Решим задачу о ее собственных колебаниях. Назначим следующие геометрические параметры: длина - 1 м, толщина несущих слоев - 1 мм, толщина слоя заполнителя - 3 см. Материал несущих слоев - сталь (Е = 210 ГПа, д = 0,28, р = 7 700 кг/м3). Изотропный материал заполнителя обладает модулем сдвига Озап = 2,9 МПа. С помощью программы [6] получены результаты решения задачи на собственные
Номер моды Тестируемое решение Аналитическое решение Решение в пакете COSMOS/M
1 147,811 147,833 147,809
2 590,982 591,332 590,956
3 1 328,730 1 330,497 1 328,605
4 2 359,758 2 365,328 2 359,361
5 3 682,269 3 695,826 3 681,293
значения при ширине сечения в 1 мм (табл. 3). Параметры инерции трехслойной балки вычисляются по формулам
Вр =(н К +Рз И) , Ср =
вр=(рн (н3 - ИЗ )+р з н\)),
где рн и рз - плотности материалов несущих слоев и заполнителя; Ин и Из - толщины несущих слоев и слоя заполнителя; Н - общая толщина пакета, Н = 2 Ин + Из.
Таблица 3
Частоты колебания трехслойной шарнирно-опертой балки
Для сравнения приводится результат (табл. 3, третий столбец) решения задачи в пакете COSMOS/M, где трехслойная конструкция моделировалась элементами PLANE2D (заполнитель) и BEAM2D (несущие слои).
Хорошее совпадение частот собственных колебаний (табл. 3), вычисленных с помощью двух разных моделей, подтверждает правильность разрабатываемого нами подхода. Здесь же приводятся результаты расчета для модели с фиксированным сдвигом на торцах (см. табл. 3, четвертый столбец). Небольшое отличие частот для балок с фиксацией сдвига и без такового характерно для граничных условий шарнирного опирания.
Выполнен модальный расчет трехслойной конструкции в пакете COSMOS/M с моделированием балки элементами слоистой пластины SHELL4L. На торцах -условия шарнирного опирания: слева - неподвижный шарнир, справа - подвижный в направлении продольной оси. Ширина сечения балки - 50 мм. Остальные размеры и механические свойства материалов такие же, как в балке предыдущего расчета.
Результаты расчетов (табл. 4) сравниваются с данными, определенными по нашей модели. Как видим, имеется заметное расхождение, что свидетельствует о том, что в модели пакета COSMOS/M не учитываются рассматриваемые нами эффекты. В данном случае большего доверия заслуживают наши результаты, поскольку они почти точно совпадают с решением, полученным в том же пакете COSMOS/M при строгом моделировании поведения заполнителя (см. табл. 3).
Рассмотрим трехслойную балку с граничными условиями жесткого защемления торцов.
при х = 0 и при х = L: w = 0,
и = 0, у = 0, dw/dx = 0. (15)
Сохраним геометрические параметры предыдущей задачи: длина - 1 м, толщина несущих слоев - 1 мм,
толщина слоя заполнителя - 3 см, ширина сечения -1 мм. Материал несущих слоев - сталь (Е = 210 ГПа, д = 0,28, р = 7 700 кг/м3). Модуль сдвига изотропного материала заполнителя примем Озап = 290 МПа.
Таблица 4
Частоты колебания трехслойной полосы
Номер моды Торцы со свободным сдвигом Решение в пакете COSMOS/M (SHELL4L)
1 20,964 22,227
2 43,292 46,028
3 65,338 69,487
4 87,306 92,829
5 109,243 116,097
Рассмотрим результаты решения задачи (табл. 5). Для сравнения приводятся частоты (см. табл. 5, третий столбец), вычисленные с помощью пакета COSMOS/M, где трехслойная конструкция моделировалась элементами PLANE2D (заполнитель) и BEAM2D (несущие слои). Защемление торцов обеспечено фиксацией угловых узлов на краях. Вновь следует отметить весьма хорошее совпадение и в количественном (табл. 5) и в качественном отношении результатов, полученных по нашей модели и в пакете COSMOS/M.
Выполнен модальный расчет для балки с защемленными торцами, на которых, однако, освобожден угол сдвига. Значения частот собственных колебаний показаны в предпоследнем столбце табл. 5. Они существенно расходятся с соответствующими величинами, вычисленными для жестко защемленной балки (табл. 5, второй столбец). Следовательно, в модальном расчете балки со вторым видом защемления для точного определения частот и форм собственных колебаний требуется рассматриваемая в нашей работе модель (постановка).
Для сравнения выполнено решение задачи для балки с шарнирным опиранием торцов, где, кроме того, допущен сдвиг (см. табл. 5, последний столбец). Как видим, результаты отличаются от тех, что получены для защемленной балки со свободным сдвигом на торцах.
Последнее решение задачи об определении частот собственных колебаний защемленной (со свободным сдвигом на торцах) трехслойной балки выполнено для весьма податливого заполнителя Озап = 2,9 МПа. Получены следующие величины пяти первых частот: f = 21,262 Гц, f2 = 43,833 Гц, f3 = 66,136 Гц, f4 = 88,365 Гц, f = 110,567 Гц. Сравнивая эти значения с частотами колебания шарнирно-опертой балки со свободным сдвигом на торцах (см. табл. 4), можно заключить, что в случае достаточно малой жесткости балки при трансверсальном сдвиге защемление торцов со свободным сдвигом можно моделировать шарнирным опиранием (тоже со свободным сдвигом).
В заключение выполним модальный анализ монолитной балки, изготовленной из однонаправленного композита, например углепластика, у которого следующие актуальные для данной задачи механические
Номер моды Торцы со свободным сдвигом Решение в пакете COSMOS/M Торцы с фиксированным сдвигом
1 20,99 20,388 21,220
2 43,304 42,278 43,805
3 65,346 64,473 66,119
4 87,312 87,22 88,341
5 109,248 111,049 110,555
свойства: Е = 180 ГПа, 0 = 5 ГПа и плотность таты этих вычислений в виде первых пяти частот собр = 1500 кг/м3. Выполним серию расчетов для различ- ственных колебаний (табл. 6). ных вариантов граничных условий. Получены резуль-
Таблица 5
Частоты колебания трехслойной балки, защемленной по торцам
Номер моды Жесткое защемление торцов Решение в пакете COSMOS/M Защемление со свободным сдвигом на торцах Шарнирно-опертые торцы со свободным сдвигом
1 134,105 130,53 77,563 70,122
2 309,411 298,818 254,025 245,059
3 518,337 497,572 477,927 467,133
4 742,287 710,749 716,150 703,270
5 972,990 931,689 956,270 941,120
Таблица 6
Частоты колебания многослойной балки при различных вариантах закрепления торцов
Номер моды Шарнирно-опертые торцы со свободным сдвигом Шарнирно-опертые торцы с фиксированным сдвигом Защемление торцов
Авторская модель COSMOS/M фиксированный сдвиг свободный сдвиг
1 49,598 49,587 49,599 111,787 96,267
2 197,494 197,335 197,521 305,515 269,699
3 441,062 440,279 441,197 592,381 533,552
4 776,122 773,732 776,538 966,163 884,950
5 1 197,208 1 191,62 1 198,192 1 420,876 1 319,846
Как и в предыдущих примерах, здесь рассматривались следующие варианты граничных условий: классическое шарнирное опирание со свободным транс-версальным сдвигом на торцах, шарнирное опирание с фиксацией сдвига на торцах, жесткое защемление торцов, защемление со свободным сдвигом на торцах.
Здесь же приведены данные расчета, выполненного для балки с классическим шарнирным опиранием в пакете С08М08/М, где монолитная композитная конструкция, колебания которой происходят в одной плоскости, моделируется элементами РЬЛКБ2Б.
Как следует из представленных значений (см. табл. 6), имеется хорошее совпадение данных нашего расчета с результатами вычислений, выполненных в пакете. Напомним, что наше решение соответствует одномерной балочной модели, а расчеты в пакете С08М08/М проводились с помощью двумерной модели, чтобы адекватно учесть поведение конструкции с низкой трансверсальной сдвиговой жесткостью.
Анализ представленных данных (см. табл. 6) подтверждает ранее сделанный вывод о том, что в общем случае модальный расчет композитных балок с неклассическими граничными условиями (шарнирно-опертые торцы с фиксированным сдвигом и защемление со свободным сдвигом на торцах) приводит к результатам, отличным от тех, которые получаются при классических вариантах этих граничных условий. Это обстоятельство подчеркивает актуальность разработанной модели.
В заключение сформулируем выводы.
1. Получен энергетический функционал для решения задачи о собственных колебаниях балок с низкой трансверсальной жесткостью.
2. Разработана конечно-элементная модель податливой при сдвиге балки, вектор узловых неизвестных которой включает углы трансверсального сдвига. Сформированы матрица жесткости и матрица инерции для соответствующего балочного элемента.
3. В результате проведенного численного исследования на примере расчета изотропных, ортотроп-ных композитных и трехслойных балок показана актуальность разработанной модели при проведении модального анализа конструкций с низкой трансвер-сальной сдвиговой жесткостью, а также при учете граничных условий неклассического вида.
Библиографические ссылки
1. Reissner E. The Effect of Transverse Shear Deformation on the Bending of Elastic Plates // Trans. ASME. J. of Applied Mechanics. 1945. Vol. 12 (2). P. 69-77.
2. Midlin R. D. Influence of Rotary Inertia and Shear on Flexural Motions of Elastic Plates // Trans. ASME. J. of Applied Mechanics. 1951. Vol. 18. P. 31-38.
3. Васильев В. В. Механика конструкций из композиционных материалов. М. : Машиностроение, 1988.
4. Комплекс программ для получения основных матриц и векторов теории МКЭ для конечных элементов балки, пластины и оболочки, податливых при трансверсальном сдвиге: свидетельство о гос. регистрации программы для ЭВМ № 2010611594 / 26.02.2010 / Нестеров В. А. Заяв. № 2009617627, 31.12.2009.
5. Конечно-элементное исследование НДС, устойчивости и собственных колебаний балки с учетом трансверсальной податливости: свидетельство о гос.
регистрации программы для ЭВМ № 2010611595 / 26.02.2010 / Нестеров В. А. Заяв. № 2009617628, 31.12.2009.
6. Исследование НДС, устойчивости и собственных колебаний трехслойной балки с податливым заполнителем с помощью МКЭ при учете трансвер-
сального сдвига в качестве узлового неизвестного: свидетельство о гос. регистрации программы для ЭВМ № 2010611781 / 09.03.2010 / Нестеров В. А. Заяв. № 2010610178, 18.01.2010.
© Нестеров В. А., 2013
УДК 512.54
СТРОЕНИЕ БЕСКОНЕЧНОЙ СИЛОВСКОЙ ПОДГРУППЫ В НЕКОТОРЫХ ГРУППАХ ШУНКОВА*
В. И. Сенашов
Сибирский государственный аэрокосмический университет имени академика М. Ф. Решетнева Россия, 660014, г. Красноярск, просп. им. газ. «Красноярский рабочий», 31 Институт вычислительного моделирования СО РАН 660036, Красноярск, Академгородок, 50. Е-mail: senll 12home@mail.ru
Изучаются группы, введенные В. П. Шунковым в 1975 г., и названные в его честь группами Шункова в 2000 г. Используется методика исследования бесконечных групп, разработанная в красноярской школе по теории групп. Цель работы - установить строение бесконечной силовской 2-подгруппы в группе Шункова, не обладающей почти слойно конечной периодической частью, когда в группе нормализатор любой конечной нетривиальной подгруппы обладает почти слойно конечной периодической частью. Доказано, что если некоторая си-ловская 2-подгруппа такой группы бесконечна, то она является расширением квазициклической 2-группы при помощи обращающего автоморфизма. Этот результат найдет применение при изучении бесконечных групп с условиями конечности. Строение искомой бесконечной силовской 2-погруппы полностью установлено.
Ключевые слова: группа, силовская подгруппа, инволюция, слойная конечность.
STRUCTURE OF THE INFINITE SYLOV SUBGROUP IN SOME SHUNKOV GROUPS
V. I. Senashov
Siberian State Aerospace University named after academician M. F. Reshetnev 31 “Krasnoyarskiy Rabochiy” prospect, Krasnoyarsk, 660014, Russia Institute of Computational Modeling of Siberian Branch of RAS ICM SB RAS Akademgorodok, 660036, Krasnoyarsk, Russia
The author studies groups entered by V. P. Shunkov in 1975 and named in his honor Shunkov groups in 2000, for the study the author uses a technique of infinite groups, developed at Krasnoyarsk School on group theory. The aim of the work was to establish the structure of an infinite Sylov 2-subgroup in Shunkov that does not have an almost layer-finite periodic part, when the normalizer of any finite non-trivial subgroup has an almost layer-finite periodic part. The author proves that if a Sylov 2-subgroup of the group is infinite, then it is an extension of a quasi-cyclic 2-group by reversing automorphism. This result will be used in the study of infinite groups with finiteness conditions. The structure of the unknown infinite Sylov 2-subgroup was completely revealed.
Keywords: group, Sylov subgroup, involution, layer-finiteness.
Слойно конечные группы впервые появились без названия в статье С. Н. Черникова [1], а затем в его последующих публикациях за ними закрепилось название слойно конечных групп. Группа называется слойно конечной, если множество ее элементов любого данного порядка конечно. Почти слойно конечные группы - это расширения слойно конечных групп с помощью конечных групп. Класс почти слойно конечных групп значительно шире класса слойно конечных групп. В то время как только некоторые чер-
никовские группы слойно конечны, все черниковские группы являются почти слойно конечными.
Здесь мы изучаем группы с условием: нормализатор любой конечной нетривиальной подгруппы обладает почти слойно конечной периодической частью. Класс групп, удовлетворяющий этому условию, довольно широк. В нем содержатся свободные бернсай-довские группы нечетного периода > 665 [2] и группы, построенные А. Ю. Ольшанским [3].
Определение. Группа называется группой Шунко-
Работа выполнена при поддержке гранта РФФИ (проект 10-01-00509) и гранта Министерства образования России (проект - алгебро-логические структуры и комплексный анализ с приложениями к передаче и защите информации).
