Определение основной частоты колебании ортотропнои пластины, закрепленной в углах

Библиографические ссылки

1. Large Deployable Reflector (LDR) / K. Shintate, K. Terada, M. Usui et al. // J. of the National Institute of Information and Communications Technology. 2003. Vol. 50, № 3/4. P. 33-39.

2. Thomson M. W. AstroMesh deployable reflectors for KU- and KA-band commercial satellites // Proc. of 20th AIAA Intern. Communication Satellite Systems Conf. and Exhibit. Montreal, Quebec, 2002.

3. Development of the large and high accuracy deployable antenna for the VSOP-2 mission [Electronic resource] / Y. Murata, H. Hirabayashi, M. C. Natori et al. URL:

http://www.ursi.org/Proceedings/ProcGA05/pdf/BP.11(01 277).pdf (date of visit: 12.12.2012).

4. Tibert G. Deployable Tensegrity Structures for

Space Applications : Ph.D. Thesis [Electronic resource] Stockholm, 2002. URL: http://www2.mech.kth.se/

~gunnart/TibertDocThesis.pdf (date of visit: 12.12.2012).

5. Lopatin A. V., Rutkovskaya M. A. Design of large space antenna composite rim // Composite Structures. 2006. № 76. P. 99-105.

6. Zakharov Yu. V., Okhotkin K. G. Nonlinear bending of thin elastic rods // J. of Applied Mechanics and Technical Physics. 2002. Vol. 43, № 5. P. 739-744.

A. V. Lopatin, Yu. V. Zakharov, K. G. Okhotkin, V. V Vilyanen, A. V. Pashkovskiy

GEOMETRICALY NONLINEAR MODEL OF AN TRANSFORMED RIM OF A BIG SPACE ANTENNA WITH FLEXIBLE COMPOSITE ELEMENTS

A new geometrical and nonlinear model of a transformable rim of a large space antenna constructed with the use of flexible composite strips, as basic elements, is considered. The decision found will allow to define the optimum form of the bent curvilinear strip and to estimate the stored energy.

Keywords:transformable space antenna, flexible bent composite element, geometrical nonlinearity.

© Лопатин А. В., Захаров Ю. В., Охоткин К. Г., Вильянен В. В., Пашковский А. В., 2012

УДК 517.972.5

И. А. Лопатин

ОПРЕДЕЛЕНИЕ ОСНОВНОЙ ЧАСТОТЫ КОЛЕБАНИИ ОРТОТРОПНОИ ПЛАСТИНЫ,

ЗАКРЕПЛЕННОЙ В УГЛАХ*

Обобщенным методом Галеркина решена задача определения основной частоты колебаний прямоугольной ортотропной пластины, в углах которой отсутствует прогиб.

Ключевые слова: ортотропная пластина, обобщенный метод Галеркина, основная частота колебаний.

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

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

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

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

* Исследование выполнено при поддержке Министерства образования и науки РФ, соглашение 14.В37.21.0405 «Моделирование композитных элементов крупногабаритных трансформируемых механических систем космических аппаратов связи и телекоммуникаций».

Рассмотрим прямоугольную ортотропную пластину, отнесенную к системе декартовых координат ху (рис. 1). Размеры пластины по осям х и у обозначим а и Ь соответственно. Пластина закреплена в угловых точках таким образом, что в них отсутствует прогиб.

2

S = J (T - U )dt,

D11

f я2 А2 d w

dx2

+ 2D1

d 2 w d 2 w 12 dx2 dy2

+ D2

f *2 А2 d w

dy2

+ 4D3

f я2 А2 d w

dxdy

dxdy,

T = Tmax cos at, U = Umax sln2 at,

где Ттах - максимальная кинетическая энергия движения пластины; итах - максимальная потенциальная энергия изгиба пластины. Величины Ттах и итах определяются следующим образом:

^ a b

Tmax = 2a2Bp JJw2dxdy,

' 0 0

D

f й2 А2 d w

3x‘

( ил у

+ 2D

д 2 w д 2 w 12 dx2 dy2

+ D.

f Л2 А2 d w

22

dy2

+ 4 D

f я2 А2 d w

33

dxdy

dxdy,

Рис. 1

Получим вариационное уравнение изгибных колебаний пластины. Воспользуемся для этого принципом Гамильтона. Рассмотрим интеграл действия Гамильтона [1]:

где Вр - инерциальный параметр.

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

t1 t2 =

2п

a

Подставляя (4) в (1) и учитывая соотношение (6), будем иметь

2я/ 2 п/

где t - время; /2 - Ь интервал времени, в течение которого происходит движение пластины; Т - кинетическая энергия движения пластины; и - потенциальная энергия изгиба пластины.

Кинетическая энергия и потенциальная энергия изгиба ортотропной пластины определяются следующим образом [2]:

I Ь / л. N 2

S = Tmax J cos2 atdt - Umax J si^ atdt =

= —

max max

a

Интегрируя по времени, найдем S = —

V max max / •

В соответствии с принципом Гамильтона в промежутке времени 2п/ю, интеграл действия £ для действительного движения ортотропной пластины имеет стационарное значение, т. е.

5S = 0,

где 5 - знак вариации.

С учетом равенства (8) из условия (9) будем иметь

8(Tmax - Umax ) = °.

где w = м>(х, у, t) - прогиб пластины; D11, D12, D22, D33 - изгибные жесткости ортотропной пластины.

Рассматривая свободные колебания пластины, представим ее прогиб в следующем виде:

w

где ю - круговая частота колебаний; w(x,y) - функция формы пластины при поперечных колебаниях. Подставляя (3) в (2), получим

Подставляя (5) в (10) и выполняя варьирование, получим

a b

Mx (^ y)Sl it

JJi( x, y)5wdx*dy + J

00 0

b

-J[Vx(x, y)8w]dy +J

0 0

- J\Vy (x,y)Sw]b dx + 2 [Mxy (x, y)5w]

My(x, y)s| dy

dy -

dx -

= 0.

Здесь

Д х, у) = £„ — + 2 (и + 2 £>33)—— +

дх дх ду

~,4

гл д V _ 2

+ £22 -"л-ВрЮ V,

22 ду4 р

д2 V д2V

Мх (х, у) = £11 — + £12 —,

^ 11 ах2 12 ду2

д3 V , . „ ч д3 V

V

дх

дхду

Му (х, у) = £12 2 + £22 2

дх ду

V- (^ у) = £22 + (°12 + 4£33 ) — 2 ду >

Мху (X, у) = 2£33

д2 V

дхоу

„ • пх, . . лу„„ . пх . пу

5^ = Б1П-----5А + Б1П — 5В + Б1П------Б1П — 50,

) п пхсл п пх . пу 5| — І = — соб—5А +— соб—Біп—50,

дх і а а а а Ь

Г дw ) п пу п . пх пу

51 ----- I =— 008--5В + Б1П 008 — 5С.

^ ду) Ь а Ь а Ь

Подставляя (14), (13) в (11), (12) и выполняя интегрирование, получим, учитывая произвольность вариаций 54, 5В, 5С однородную систему линейных алгебраических уравнений

.3 Г и2 \

С -

п4 п2 4п

2 Д1^- А + 16Д,^— В + -4-—

11 4 12 2.2 2.2

а а Ь а Ь

£11 2 + £12 V а I

В равенствах (12) функции Мх, Му соответствуют изгибающим моментам, Ух, Уу соответствуют обобщенным перерезывающим силам, Мху соответствует крутящему моменту.

Уравнение (11) представляет собой основное вариационное уравнение, которому должны удовлетворять собственные формы колебаний w(x,y) ортотроп-ной пластины.

Представим прогиб пластины (рис. 2).

- В„ю21 2А +16В + -0 | = 0, 1 п2 п

П п4 4п3 Г а2

16£12 2,2 А + 2£22Т4В + ТГ £22 72 + £12 а2Ь2 Ь а Ь . Ь2

2

0 -

Р I -п-2 '

\П

- В ю2 116 А + 2В+-0 |= 0,

4п

2 »2 а Ь

2

\

£11 2 + £12

к а2 ,

2

А +

4п

(

2,2 а Ь

2

а

£22 ,2 + £12 V Ь I

В +

2 »2 а Ь

2

£11 2 + 2 (12 + 2£33 —+ £22 , 2

а Ь

- Вр ю2 | - А + - В + 0 | = 0.

0 -

Преобразуем уравнения (15) таким образом, чтобы коэффициенты при неизвестных А, В, 0 стали безразмерными. Для этого разделим каждое из уравнений (15) на ^/£п£22. В результате преобразований полу-

Рис. 2

Выбранная аппроксимация

. . пх „ . пу . пх . пу

W = 4 Б1П----+ В Б1П — + С Б1П Б1П —,

а Ь а Ь

где А, В и С - неизвестные числа, не удовлетворяет статическим граничным условиям на свободных краях пластины. Вместе с тем представление прогиба в виде (13) позволяет вполне успешно определить основную частоту колебаний закрепленной в углах ор-тотропной пластины, если для решения этой задачи использовать обобщенный метод Галеркина. Исходным уравнением для обобщенного метода Галеркина является вариационное уравнение (11). Это уравнение автоматически обеспечивает приближенное выполнение статических граничных условий на свободных краях пластины.

Вариации функции прогиба (13) и ее первых производных определяются следующим образом:

чим

2п4аА + 16л2Р12В + 4п3 (а + Р12 ) 0 --П 2 А + 16 В + -01 = 0,

16л2р12 А + 2п4 — В + 4п31 — + Р12| 0 -а Vа і

-П^ -6 А + 2В +^^ 01|=0,

4п3 (а + Р12 ^^4 + 4п31— + Р12 В +

+ п4 а + 2 (Р12 + 2Р33) + —

_ а

Здесь

0-п| - А + - В + 0 | = 0.

а^£11 ^Г, Р12 =

Д-

£

•\/£11£2

Р33 = £33

•\/£11£22

+

а

Величины а, Р12, Рзз - безразмерные геометрические и жесткостные параметры ортотропной пластины.

Величина п в уравнениях (16) определяется равенством

П =

Вр ю2 а2Ь2

х/а'лГ

Это безразмерный частотный параметр. Таким образом, задача определения основной частоты изгиб-ных колебаний ортотропной пластины, закрепленной в углах, сведена к нахождению параметра п, при котором однородная система (16) будет иметь решение, отличное от нуля.

Вычисление частотного параметра п может быть выполнено несколькими способами. Представим уравнения (16) в следующем матричном виде:

Ж - ^Х = 0.

Здесь

2п4а 16п2р12 4п3

И = 16п2р12 2п4 — а 4пз (а+Р12 ^

4п3 (а + Р12 ) 4п3 (а+р12) 4 П а + 2 (Р12 + 2Рзз ) + _ а

G =

16

П2

2

16

п2

4 4 і

4')

П

4

П

X =

Г А')

В

С

Вр

где £ = д/л.

В качестве примера определим частотный параметр для изотропной пластины, толщина которой равна И. Упругими характеристиками изотропного

материала являются модуль упругости Е и коэффициент Пуассона V. Плотность материала обозначим как р. Изгибные жесткости и инерциальный параметр для изотропной пластины определяются следующим образом:

-И

£11 = £22 = Е

12

3

а2=Е——• 12 2 12

£»12 = Е V

Вр = РИ,

12

где Е = -

Е

1 - V

С учетом равенств (19) безразмерные параметры а, р12, Р33 примут значения

= 1 и = я _1 -V

а = 2 , Р12 =У, Рзз = ~ • с 2

где с = а/Ь - отношение сторон пластины.

Частота колебаний изотропной пластины определяется по формуле

При традиционном способе величина п определяется как наименьший корень кубического уравнения, которое получено из условия det(H - "^) = 0. В данной работе параметр п находится как наименьшее собственное число обобщенной задачи на собственные значения (18).

Из анализа коэффициентов (16) следует, что величина частотного параметра п зависит от безразмерных параметров а, р12, р33. Значения этих параметров определяются жесткостными и геометрическими характеристиками ортотропной пластины.

Когда частотный параметр п найден, основная частота колебаний может быть получена из формулы (17), т. е.

Из равенств (20) следует, что частотный параметр £ для изотропной пластины, при заданном коэффициенте Пуассона V, зависит только от отношения с. Полагая V = 0,3, исследуем влияние отношения сторон с, изменяющегося от 1 до 2,5, на величину частотного параметра £.

Приведем результаты вычисления частотного параметра с£ для некоторых значений с (см. таблицу). В справочнике [3] можно найти аналогичные результаты, полученные с помощью метода конечных разностей. Из сравнения данных следует, что величины с£, найденные двумя способами, отличаются не более чем на 1,5 %. Это дает основание заключить, что предложенный в статье подход дает вполне приемлемые значения частотного параметра.

Вычисление частотного параметра

с с^

1 7,22

1,5 9,05

1,0 9,10

2,5 9,49

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

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

Библиографические ссылки

1. Leissa A.W. Vibration of plates / Acoustical Society of America. Washington, D. C., 1993.

2. Vasiliev V. V. Mechanics of composite structures. Washington, D. C. : Taylor & Francis, 1993.

3. Blevins R. D. Formulas for natural frequency and mode shape. Malabar : Krieger Publishing Company, 1979.

I. A. Lopatin

DETERMINATION OF VIBRATIONS DOMINANT MODE OF AN ORTHOTROPIC PLATE FASTENED TO THE BEDDING ANGLES

The problem of determination of the fundamental frequency of vibration of a rectangular orthotropic plate, fastened to bedding angles, is solved by the generalized Galerkin method.

Keywords: orthotropic plate, generalized Galerkin method, fundamental oscillation frequency.

© Лопатин И. А., 2012

УДК 62.501

А. В. Медведев

И-МОДЕЛИ ДЛЯ БЕЗЫНЕРЦИОННЫХ СИСТЕМ С ЗАПАЗДЫВАНИЕМ

Рассматривается проблема моделирования нового класса процессов, имеющих «трубчатую» структуру в пространстве «входных-выходных» переменных. Модели процессов этого класса существенно отличаются от общепринятых параметрических моделей, представляющих собой поверхности в том же пространстве. Специально анализируется вопрос о моделировании многомерных систем при наличии малых объемов обучающих выборок. Для построения обучающихся параметрических моделей «трубчатых» процессов вводится соответствующий непараметрический индикатор. Приводится новый класс обучающихся параметрических моделей и некоторые результаты их численного исследования.

Ключевые слова: идентификация, непараметрические алгоритмы, непараметрическая модель, априорная информация, стохастические процессы, Н-модели.

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

Например, ряд переменных измеряется электрическим способом (в этом случае дискретность контроля At может быть достаточно мала), а другие переменные контролируются в ходе химического анализа или физико-механических испытаний (в этом случае дискретность контроля АТ велика, т. е. АТ >> At). Тогда исследуемый объект может быть представлен статическим с запаздыванием. Такой процесс целесообразно по соответствующему каналу представить в виде

где х^) - выходная переменная объекта; и^ -т) -входная переменная; т - запаздывание; ^У) - случайное возмущение, действующее на объект; t - непрерывное время.

Покажем каналы измерения с помехами Ии , Их и дискретностью измерения АТ >> At (рис. 1).

x

Рис. 1.