CC BY
33
Если этот алгоритм выдает t, то функция f будет иметь более высокий порядок КИ, нежели f Известно, что аффинная добавка никак не влияет на порядок КР.
Итак, изложены результаты реализации и исследования двух генераторов частично бент-функций, описаны возможности по улучшению порядков КИ и КР таких функций. Работоспособность алгоритмов проверена на практике, изучены свойства генерируемых функций.
Автор выражает благодарность И. А. Панкратовой и К. В. Сафонову за постоянное внимание к работе и ценные замечания.
Библиографические ссылки
1. Carlet C. Boolean Functions for Cryptography and Error Correcting Codes [Electronic resource] // INRIA. 2010. URL: http://www-rocq.inria.fr/codes/Claude.Carlet/ chap-fcts-Bool-corr.pdf. (дата обращения: 08.06.2010).
2. Токарева Н. Н. Бент-функции: результаты и приложения. Обзор работ // Прикл. дискретная математика. 2009. № 1 (3). С. 15-36.
3. Bent-function [Electronic resource] // Wikipedia. 2010. URL: http://en.wikipedia.org/wiki/Bent_function (дата обращения: 08.06.2010).
4. Zheng Y., Zhang X. M. Plateaued functions // ICICS’99. Lecture Notes in Computer Science. 1999. Vol. 1726. P. 284-300.
5. Carlet С. Partially-bent functions // Design, Codes and Cryptography. 1993. Vol. 3. №. 2. P. 135-145.
6. Токарева Н. Н. Обобщения бент-функций. Обзор работ // Дискрет. анализ и исследование операций. 2010. Т. 17. № 1. С. 34-64.
7. Агафонова И. В. Криптографические свойства нелинейных булевых функций [Electronic resource] // Семинар по дискретному гармоническому анализу и геометрическому моделированию DHA & CAGD. 2007. URL: http://dha.spb.ru/PDF/cryptoBOOLEAN.pdf (дата обращения: 08.06.2010).
M. V. Naumov
THE GENERATION OF PARTIALLY-BENT FUNCTIONS
Most important characteristics of cryptographic functions are balancedness, nonlinearity, propagation criterion, correlation immunity, degree and non-existence of nonzero linear structure. Partially-bent functions form super-class of the class of bent functions. These functions may achieve desirable characteristics.
Two algorithms for generation of partially-bent functions were supposed and studied. Second algorithm may improve cryptographic characteristics of generated functions.
Keywords: Boolean functions, partially-bent functions, correlation immunity.
© Наумов М. В., 2010
УДК 519.62
В. А. Нестеров
ПОПЕРЕЧНЫЕ КОЛЕБАНИЯ ПЛАСТИНЫ, ПОДАТЛИВОЙ ПРИ ТРАНСВЕРСАЛЬНОМ СДВИГЕ
Рассматривается конечно-элементный модальный расчет пластины с низкой трансверсальной сдвиговой жесткостью. В каждом из четырех узлов прямоугольного конечного элемента пластины в качестве основных кинематических параметров присутствуют углы трансверсального сдвига. На примере анализа собственных колебаний композитной пластины показана актуальность разработанной конечно-элементной модели. Представлены результаты численного исследования влияния граничных условий неклассического вида на величины частот собственных колебаний.
Ключевые слова: пластина, метод конечных элементов, трансверсальный сдвиг, модальный анализ.
Композитные конструкции, обладающие высокой степенью весового совершенства, часто используются в производстве космической техники. Композитные пластины отличаются низкой сдвиговой жесткостью по отношению к трансверсальным напряжениям. Учет указанной особенности при реализации численного расчета приводит к повышению порядка разрешаю-
щих уравнений за счет введения в рассмотрение углов трансверсального сдвига.
Разрешающие уравнения теории метода конечных элементов (МКЭ) для задачи о собственных колебаниях пластины с низкой трансверсальной сдвиговой жесткостью получим вариационным способом. Для этого запишем выражение полной энергии колеблющейся пластины
E = U + T ,
(1)
где и - потенциальная энергия деформации; Т - кинематическая энергия движения.
Выражение для потенциальной энергии деформации пластины как трехмерного тела имеет вид
U = 2 í í í (Öe +ÖЛ +Özez +
(V )
+1 e + 1 e + 1 ё I dxdydz,
xy xy xz xz yz yz f у ’
düx
dx
düy
dy
düz
dz
du*.+dü±
dy dx
üx = üx (x, y, z, t); üy = üy (x, y, z, t); üz = üz (x, y, z, t) .
Подставляя (5) в (4), получим: ё = e sin rot,
E„
E..
az °y ax
z■; e = — — u x
yy
E
E
E
®z ®x °y Txy
e =—^ — u — u^~ ; e =——
z e e zy E G
1 yz
eyz = G-
Gyz
G,_
(2)
где напряжения и деформации являются функциями координат и времени
ст = ст( х, у, z, t) ; ё = ё( х, у, z, t), (3)
гДе ст - любой параметр из ах, сту, а,, Тху, Ху,, ^ ;
ё - любой параметр из ёх, ёу, ё,, ёХу, ё^, ёа .
Г еометрические соотношения имеют вид
где Ех(у,г) - модуль упругости соответствующего направления; Gxy(yг,xг) - модуль сдвига в соответствующей плоскости; Ду, Дух, Дг, Ду, Дг Дх - коэффициенты Пуассона.
Имеет место свойство симметрии упругих постоянных:
ЕхРху = Еу Vух ; ЕУ Руг = Е Ргу '; ЕхРхг = Ег Ргх . (12)
Используя принятые в [1] допущения, выражение потенциальной энергии деформации (2) можно привести к виду
дйу дйъ дй дй
е =—- + —-; е =—- + —-, (4)
дъ ду дъ дх
где йх, йу, й2 - проекции перемещения произвольной
точки на соответствующие оси координат. Они также являются функциями координат и времени:
U = U sin2 rot, где U определяется выражением
U=і íllL fe+ z ^ )+„ y fe+z ^
2 J J -М I dx dx ¡ ^dv dy
(13)
(V) +1.
dx
dü dv í dO dO
- + — + z
dy dx ^ dy dx
+ Txz'Vx +1 yzVy } dxdydz.
(14)
(5)
При гармонических колебаниях закон изменения перемещений по времени можно представить в следующем виде:
üx = их sinrot; üy = иу sinrot; üz = uz sinrot, (6)
где ю - круговая частота колебаний, а перемещения их, uy, uz зависят только от координат:
üx = üx (х,y, z); иу = uy (x, y, z); uz = uz (x, y, z). (7)
(8)
где e - любой параметр из ex, ey, ez, exy, eyz, exz, для которых справедливы геометрические соотношения
ди диу ди ди диу
e = —-; e = —— ; e = —- ; e = —- H------------- ;
x ^ y z xy *
дх ду dz ду dx
диу ди2 дих ди2
ez =—^ + —-; eXI =—- + —-. (9)
yz ^ ~ xz ^ ’
dz ду dz дx
В силу линейности физических соотношений
можно записать
д = а sin го/, (10)
где с - любой из параметров Ох, Оу, Oz, Txy, Tyz, Txz, для
которых справедливы физические соотношения
В^іражение кинетической энергии колебания пластины имеет вид
T = \ i i ¡(vX + vy + )р dxdydz ’ (15)
2 (V)
где р - плотность материала (в общем случае является функцией координат); vx, vy, vz - проекции вектора
скорости на соответствующие оси координат, связанные с перемещениями дифференциальными зависимостями
_ дйх дйу _ düz
vx =—-; vy =—-; vz =—-. (16)
x dt y dt z dt
Подставляя (5) в (16), а затем результат - в (15), получим
T = T cos2 rot, (17)
где
T =1 í í í ü + ü2y + üI^ pro2 dxdydz .
2 (V)
(18)
Приравнивая максимальные значения потенциальной и кинетической энергии, которые определяются выражениями (14) и (18), получим
y ííí
du cQ,
ax I — + z—- I + & x
xI dx dx 1 x
+ir
E
(V)
dü dv í dQx dQy — + — + z I —x + —y dy dx ^ dy dx
- + z-
dy dy
+ 1xzVx +1 yzVy }>
e
ex =
ez =
y
+
y
+
+
+
a
a
a
y
x
z
:dxdydz-—] ] |р®2 [(и + гвх)2 +
+ (у + гв у )
(V) 2
+ м
dxdydz = 0.
При выводе (19) учтено, что
их = и + гвх ; иу = V + гв ; и2 = м,
(19)
(20)
дм дм
в = у---------------; в = у-------------------.
х х дх у у ду
(21)
Проинтегрируем (19) по толщине пластины, т. е. по координате г в пределах от -5 до h - s. В результате преобразований получим
1 а Ь
2 И ( Ы 8 х + N 8 у + Нху 8 ху + Мх г х +
2 0 0
+Му г у + Мху 1ху + х + у ) ^у -
а Ь
т ю2 И{вр(и2+v 2+м2)+
—ю 2
(22)
+2С
дм) ( дм
и\ Ух -~ М уу - —
дх
D„
дм ) ( дм
Ух ) +(уу -ау
■ dxdy = 0.
Здесь Ы, М, Q - внутренние погонные усилия, определяемые следующими соотношениями:
Ых = В118х + В128у + СИ%х + С12гу ;
Ыу = В218х + В228у + С21гх + С22гу ;
Ыху = Взз8ху + сзИху; & = кх у х; & = ку у у;
Мху = С338ху + D33гxy ;
Мх = С118х + С128у + Мдх + Dl2гу ; Му = С218х + С228у + D2lгх + D22гу ,
(23)
ди
ду
ди ду
дв
х ^ ; у ^ ; ху ^ ^ ; ^^х ^
дх ду ду дх дх
св..
50 дв у
г у ; г ху --у+^. (25)
В итоге с учетом (21) получим
] ] к (ди ) + 2см ди (ду, -^)+^1 +
\дх) дх \ дх дх ) \ дх дх )
где и и V - перемещения точек начальной плоскости вдоль осей X и 7 соответственно; м - прогиб точек начальной плоскости; 0х, 0у - углы наклона сечения, определяемые соотношениями
ди ду
' дх ду
ди ( ду д2м
дх ( ду ду
ду ( дух 52м
ду ( дх дх
+ 2Д,^1(5Уу-52м )+В. (| ) +
дх дх И ду ду'
+2С22 —
ду ( ду д2м
ди ду
ду ( ду ду 2
2 | + D22
ду у д- м
-— | + (26)
+В33\ду + дх1 + 2С33 [ду + дх
ди ду
ду ду2
5ух , дуу - 2 д
дх ду дхду
+ D,.
-2
м
дх ду дхду
2 ] ] \Вр(и2 +у м) + 2Ср
+ Кх + Ку у2у } ^у -
дм
дм
и\ ух -дх |+у 1уу -^у
+D„
дм )2 ( дм^
ух-- 1+1уу ~~ду
дх
>сЬсф = 0.
где В, С, D и К - параметры мембранной, смешанной, изгибной и трансверсальной жесткости соответственно; Вр, Ср, Dp - параметры инерции пластины, вычисляемые по формулам
И-5 И-5 И-5
Вр = ] р dz ; Ср = ] р zdz ; Dр = ] р г2 dz. (24)
- 5 - 5 - 5
Сделаем последние преобразования в функционале (22). Приведем его к виду, удобному для реализации КЭ процедуры. Для этого с помощью физических соотношений (23) выразим усилия через деформации, которые также выразим через перемещения и углы трансверсального сдвига с помощью геометрических соотношений
Функционал (26) позволяет получить разрешающие уравнения для задачи о свободных колебаниях пластины.
Выполним модальный расчет пластины с помощью МКЭ. Будем рассматривать четырехузловой конечный элемент пластины, для которого вектор узловых кинематических параметров имеет вид
бе ={ 61 62 6з Ь4}Т, (27)
где бг- (I = 1, 2, 3, 4) - вектор неизвестных г-го узла
6 ' Ь (I), [| 1 ухг уи ■ (28>
Кинематические переменные внутри элемента представим следующим вектором
I дм дм V
6 ' г * ¿у у- уу и уI. (29)
Его компоненты определяются через узловые значения (28) с помощью соотношения
6 = Р 6е, (30)
где Р - матрица размерностью 7 х 28. Это так называемая матрица функций формы.
Подставим выражения для компонент вектора 6 (29) в функционал (26), выполним интегрирование по площади элемента, в результате получим функцию полной энергии колебания конечного элемента пластины. Эта функция зависит от компонент вектора узловых кинематических параметров (28), минимизация по которым приводит к системе однородных алгебраических уравнений:
+
2
м
+
+
+
К 6е - ш2 Se 6е = 0, (31)
где Ке - матрица жесткости; Se - матрица инерциаль-ных параметров конечного элемента пластины (обе -симметричные, размером 28х28).
В выражениях для компонентов матрицы Ке фигурируют параметры жесткости, а в выражениях для компонентов матрицы Se - параметры инерции. Например:
к 2 14Dзз а2Ь2 + 5D1УЬ2 а2 + 10D11Ь4 + 10Dууa4
11 = 5 аь ’
s, =
1 727Bpb a2 + 5 520Dpb + 5 520Dpa2
12 600 аЬ
Если структура пакета пластины симметрична относительно срединной плоскости и система координат связана со срединной плоскостью, то в этом случае смешанные жесткости пластины С и инерционный параметр Ср становятся равным нулю, и структура матрицы инерции упрощается за счет появления в ней множества нулевых компонентов.
Представим матрицу жесткости и матрицу инерции конечного элемента пластины в блочном виде:
Ki-i Ki-ii KI-III Ki-iv
Kn-I Kii-ii Kii-iii KII-IV
Kiii-i KIII-II Kiii-iii KII-IV
Kiv-i KIV-II KIV-III KIV-IV
Se.
Si-i Si-ii Si-iii SI-IV
Sii-i Sii-ii Sii-iii SII-IV
Siii-i Siii-ii Siii-iii SIII-IV
SIV-I SIV-II SIV-III SIV-IV
(32)
где К- и Sч (V = I, II, III, IV) - подматрицы размером 7 х 7. Ненулевые компоненты этих матриц определяются теми же выражениями, что и соответствующие компоненты матриц жесткости и инерции, фигурирующих в уравнениях системы (31).
Глобальная система разрешающих уравнений в задаче модального расчета пластины имеет вид
К А - ш2 SE А = 0, (33)
где А - глобальный вектор узловых неизвестных,
А ' (б! б2 ... 6 ... 6*}Т, (34)
где N - общее число узлов в системе; К - глобальная матрица жесткости; S■z - глобальная матрица инерции.
Система (33) в математическом смысле является обобщенной задачей на собственные значения, решая которую (после учета в ней граничных условий), определим частоты и формы собственных колебаний пластины.
Протестируем полученный алгоритм решения на примере модального расчета изотропной квадратной пластины, жестко защемленной по контуру. Зададим следующие размеры пластины и механические свойства материала: длина стороны пластины а = Ь = 1 м,
толщина пластины h = 3 мм, модуль упругости E = 210 ГПа, коэффициент Пуассона ц = 0,3, плотность р = 7 800 кг/м3.
В результате расчета, выполненного по нашей модели с двадцатью элементами вдоль каждой из сторон пластины, определены первые пять значений собственных частот (табл. 1) и соответствующие им формы колебаний. Здесь в качестве формы колебаний будем рассматривать распределение функции прогибов. Для оценки полученных результатов выполнен расчет в пакете COSMOS/M, где для моделирования использованы элементы тонкой оболочки SHELL4.
Таблица 1
Значения частот собственных колебаний
Номер моды Частоты собственных колебаний, Г ц
Тестируемое решение Решение COSMOS/M
1 26,908 26,976
2 54,831 55,017
3 54,831 55,017
4 80,483 81,112
5 98,287 98,630
Результаты модального расчета, выполненного с помощью нашей модели (рис. 1), имеют хорошее совпадение с решением в пакете COSMOS/M (рис. 2). Это касается и частот собственных колебаний (см. табл. 1) и форм мод (рис. 3-6).
Прогиб
Рис. 1. Форма колебаний, соответствующая первой собственной частоте
Рис. 2. Форма колебаний, соответствующая первой собственной частоте (COSMOS/M)
Ke =
Рис. 3. Форма колебаний, соответствующая второй собственной частоте
Выполним численное исследование влияния вида граничных условий на величины частот собственных колебаний квадратной композитной пластины. Будем сравнивать два вида защемления: классическое (I) и со свободным сдвигом (II). Примем длины сторон пластины а = Ь = 1 м, а ее толщину И = 10 мм.
Будем полагать пластину изготовленной из однонаправленного композита с углами укладки ф = ±45о. Зададим следующие механические свойства (углепластик):
- модуль упругости вдоль волокон Е— = 180 ГПа;
- модуль упругости поперек волокон Е2 ' 6,2 ГПа;
- модуль сдвига 012 = 5 ГПа;
- коэффициент Пуассона = 0,007;
- плотность р = 1 500 кг/м3.
Рис. 5. Форма колебаний, соответствующая четвертой собственной частоте
Рис. 4. Форма колебаний, соответствующая третьей собственной частоте
Рис. 6. Форма колебаний, соответствующая пятой собственной частоте
Упругие параметры пластины, которую будем считать однослойной и условно однородной, вычислим по формулам
A11 = E1 cos4 ф + E2 sin4 ф + 2 (Е1ц12 + 2G12j sin2 ф cos2 ф; A22 = E1 sin4 ф + E2 cos4 ф + 2 ( Е1ц12 + 2G12 j sin2 ф cos2 ф; A2 = A = E1^12 +|E1 + E2 - 2 ( E^ + 2G12 jj sin2 ф cos2 ф; A33 =^E1 + E2 - 2 Е1ц12 j sin2 ф cos2 ф + G12 cos2 2ф,
где Е— и Е2 - приведенные модули упругости, вычисляемые по формуле
1(2)
Е = -
1(2) , -
1 -М-12^21
а коэффициент Пуассона ^ определяется из условия симметрии упругих постоянных:
,, =^12 Е1 г^21 р .
Е2
Модули трансверсального сдвига примем равными
Коэффициенты жесткости пластины вычислим по формулам
В1— = А1— И ; В22 = А22И ; В21 = А21И ;
„ , А1—И3 А22И3
В33 = А33И ; D11 = ——; Д_ = - 22 •
12
12
12
12
^ = G12Ä:
12
СЦ = 0 (Ц = 11, 22, 12, 21, 33),
а параметры инерции - по формулам
Вр = ph ; Ср = 0 ; Dp = р^/12
Результаты вычислений в виде значений первых пяти частот собственных колебаний поместим в таблицу (табл. 2). Здесь же приведены частоты для пластин с толщинами 30 и 50 мм.
Формы колебаний соответствуют тем, что приведены на рис. 2-6.
Анализ представленных результатов говорит о том, что в задаче модального расчета результаты, полученные при защемлении со свободным сдвигом краев, отличаются от тех, что определены при классическом защемлении. Это различие нарастает по мере увеличения толщины пластины, достигая максимума (по первой частоте) в рассматриваемой модели в 10,6 %.
Проведенная работа позволяет сформулировать следующие выводы.
1. Получен энергетический функционал для решения задачи о собственных колебаниях пластин с низкой трансверсальной жесткостью.
2. Разработана конечно-элементная модель податливой при сдвиге пластины, вектор узловых неизвестных которой включает углы трансверсального сдвига. Сформированы матрица жесткости и матрица инерции для соответствующего элемента пластины.
3. В результате проведенного численного исследования на примере расчета композитной пластины показана актуальность разработанной модели при проведении модального анализа конструкций с низкой трансверсальной сдвиговой жесткостью, а также при учете граничных условий неклассического вида.
Библиографическая ссылка
1. Васильев В. В. Механика конструкций из композиционных материалов. М. : Машиностроение, 1988.
Таблица 2
Значения частот собственных колебаний композитной пластины
Номер моды Значения собственных частот, Гц
h = 10 мм h = 30 мм h = 50 мм
Полное защемление Защемление со свободным сдвигом Полное защемление Защемление со свободным сдвигом Полное защемление Защемление со свободным сдвигом
1 111,24 110,42 323 306 509 455
2 222,09 220,72 634 608 975 900
3 222,09 220,72 634 608 975 900
4 349,33 347,40 977 943 1 457 1 370
5 380,66 378,74 1 066 1 033 1 592 1 512
V. A. Nesterov
VIBRATION OF A PLATE WITH LOW TRANSVERSE SHEAR STIFFNESS
Modal finite element analysis of a plate with low transverse shear stiffness is considered. There are two transverse shear strains at each offour nodes. In modal analysis of a composite plate the urgency of the developed finite element model is shown. Results of numerical research of influence of nonclassic boundary conditions are presented.
Keywords: a plate, finite element method, transverse shear strains, modal analysis.
© Нестеров В. А., 2010
