Численное моделирование распространения упругих волн в слоистых средах с податливыми прослойками Текст научной статьи по специальности «Физика»

тенны для выбора разрешённых значений параметров n и r, что подробно рассмотрено в [4].

Последовательные формы изогнутого гибкого элемента обода при r = 2, 3, 4 для случая п= 18, построенные по аналитическим выражениям (8), приведены на рисунке 2.

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

Литература

1. Попов Е.П. Теория и расчёт гибких и упругих стержней / Попов Е.П. - М.: Наука, 1986. - 296 с.

2. Захаров Ю.В. Нелинейный изгиб тонких упругих стержней / Захаров Ю.В., Охоткин К.Г. // ПМТФ. 2002. Т. 43. №5. С. 122-131.

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

4. Геометрическинелинейная модель трансформируемого обода большой космической антенны с гибким и композитными лентами / Лопатин А.В., Захаров Ю.В., Охоткин К.Г. // ВестникСибГАУ. 2012. Вып.5(45). С. 75-80.

Analitical model of flexible transformable hoop of spec aerial

Valentina Valentinovna Viliyanin, master

Siberia State Aerospace University after academician M.F. Reshetnev

Alexandr Konstantinovich Nikulin, graduate student Siberia State Technological University

In this article was used model of non-linear pendulum^ for^ finding optimal shape of a curved strip. The shapes with bend point presents via elliptic functions. Result of the article allows to improve construction of an aerial of satellite.

Keywords: geometric nonlinear, elliptic functions, column's Eyler.

УДК 539.374

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

Мария Александровна Похабова,аспирант Тел.: 8923 3498959, е-та11:роккаЪоуа-таг1уа@та11.ги Институт вычислительного моделирования СО РАН Ы1р://1ст.кта8п.гп

гибкого элемента с защемлением при r = 2, 3, 4 для случая n = 18

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

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

1. Введение

Характерным для горных пород является их блочное строение. При этом горный массив рассматривают как систему вложенных блоков разных масштабных уровней, от размеров кристаллических зёрен до блоков горного массива, выделяемых крупными разломами, соединённых прослойками, состоящими из более слабых трещиноватых пород [1].

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

Цель настоящей статьи состоит в исследовании особенностей формирования упругих волн в блочных средах на основе математических моделей, учитывающих сложные реологические свойства материала прослоек, и при их импульсном возбуждении проверить справедливость гипотезы о возможности существования волн маятникового типа экспериментальному и теоретическому исследованию которых, посвящены работы [2-5].

Следует отметить, что теория неоднородных сред со слоистой и блочной микроструктурой - область механики, интенсивно развивающаяся в течение более чем полувека. К настоящему времени в этой области построены адекватные математические модели, развиты эффективные аналитические методы исследования, разработаны численные алгоритмы для решения квазистатических и динамических задач, основанные на конечно-элементной аппроксимации непрерывных моделей [6-9]. Однако в задачах о распространении в слоистых средах высокочастотных волн, длины которых сравнимы с размерами блоков и прослоек, такие методы оказываются практически непригодными из-за методических погрешностей, обусловленных влиянием аппроксимационной вязкости, нередко превосходящей вязкость физическую, подлежащую учету в рамках используемой математической модели.

2. Уравнения одномерных движений слоистой среды

На рис. 1 изображена условная схема иерархического строения горной породы, которая в идеальном случае представляет собой вложенную слоистую структуру с инвариантным отношением характерных масштабов блоков и прослоек. Рассмотрим сначала отдельный фрагмент этой структуры в одномерном приближении - чередующуюся систему из п упругих слоев толщины И и уп-

л /п п 5

И

Рис. 1. Иерархическая слоистая структура горной породы

ругих прослоек толщины 8. Пусть р и р о, с и с о, а =

,2

и ао =

рс

Р 0 с0

- плотности, скорости звука и упругие

податливости материалов слоя и прослойки соответственно. Внутри слоя с номером к выполняются одномерные уравнения теории упругости

Р-

дvк до

к

док д\>

к

дг

дх

а-

дг

дх к

(1)

записанные относительно продольной скорости V в направлении оси х ( х меняется от

0 до к в пределах каждого слоя) и нормального напряжения ок. Поведение материала прослойки описывается с

помощью системы обыкновенных дифференциальных уравнений

Ро

Ж

дvк+1 + Vк

к+1 к о —о

а0

Ж до к+1 + о к

V к+1 — V к

(2)

Ж 2 8 Ж 2 8

в которую входят граничные значения введённых выше скоростей и напряжений, левые для (( + 1)-го слоя и правые для к -го слоя. Такая система может быть получена методом осреднения уравнений упругой среды в приближении тонкой прослойки (8«к). Она учитывает инерционные свойства прослойки. Упрощённое описание без учёта инерционных свойств вытекает из (2) в пределе при стремлении удельной массы прослойки Р08 к нулю. Тогда, очевидно, выполняется условие непрерывности напряжения на границе раздела слоёв. Если модуль податливости прослойки а08 также устремить к нулю, то из системы (2) последуют условия непрерывности скорости. Таким образом, получится малоинтересная для целей данной статьи математическая модель однородной упругой среды, не описывающая блочную структуры.

К системе уравнений (1) и (2) добавляются начальные условия Vк =ок = 0, (( = 1,...,п) и граничные условия следующего вида:о1 (0,г)=—р(г), Vп(к,г) = 0, где р(г) -заданное внешнее давление.

Численное решение задачи строилось на основе схемы распада разрыва Годунова на равномерной сетке с выбором предельно допустимого по условию Куранта-

Фридрихса-Леви шага по времени т = . Схема в этом случае в пределах слоя не

обладает искусственной диссипацией энергии. При меньших значениях шага использовалась кусочно-линейная БКО-реконструкция (БввепРаПуКопОвсПЫогу) второго порядка точности, которая, как показали расчёты, существенно снижает эффект сглаживания пиков численного решения с соответствующим уточнением получаемых результатов.

Условия согласования на границах раздела слоев и прослоек в форме системы

уравнений (2) также рассчитывались по схеме Годунова. Для этого в каждой из искусственно введённых ячеек длины 8, имитирующих отдельную прослойку, была реализована схема распада разрыва с независимым шагом по времени

«т, предельным по условию

т 0 = 8

с 0

Рис. 2. Сеточно-характеристическая схема со сбалансированными шагами по времени

Куранта-Фридрихса-Леви для материала прослойки. В прослойках производился расчёт такого количества шагов по времени, которое необходимо для достиже-

1

1

ния очередного временного слоя t + т основной схемы.

Сеточно-характеристическая интерпретация применяемого метода схематически приведена на рис. 2.

На этапе решения системы (2) с шагом т 0 использовались уравнения распада разрыва на границах раздела сред ("предиктор" схемы в прослойке):

к+1 , _к+1

2оv+ - а+ = 2оV - а , zv + + а + = zv

+ а

кк

2оv_+а+ = 2оv + а , 2У-+а-= zv + а .

Здесь 2 = рс и 2 о = р о со - акустические импедансы материалов основного слоя и

прослойки, величины с верхними индексами относятся к приграничным ячейкам взаимодействующих слоёв, с нижним индексом "+" - к правой, а с индексом "-" - к левой границе прослойки. Последующий пересчёт решения ("корректор" схемы) производился по формулам:

V = v + (а+ _а_)

а = а + (+^_)

Ро5' 4 " «о5'

в которых величины с крышкой относятся к новому временному слою. Предик-торные значения сеточных величин основной схемы с шагом т на границах слоёв вычислялись путём осреднения значений, относящихся к границам ячеек на мелких шагах. «Корректор» основной схемы в слоях выполнялся обычным способом на основе интегральных аналогов дифференциальных уравнений (1).

3. Вязкоупругие прослойки

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

Вязкоупругое взаимодействие по модели Максвелла (рис. 3), согласно которой деформация прослойки складывается из упругой и вязкой составляющих, Рис.3. Реологические схемы таимо- описывается системой (2) после замены второго урав-действия блоков: модель Максвелла нения более общим уравнением

«о

л ^к+1 , к а а + а

Л

2

■ + ■

к+1 к а + а

2п

V к+1 -V к

5

(3)

где ^ — коэффициент вязкости материала прослойки. В этом случае при численном решении вместо (3) на шаге "корректор" использовалась аппроксимация у равнений вяз-коупругой среды по следующей схеме Кранка-Николсон:

Ро

V _ V

а+ _а_

а _ а а+а

«о-+ -

v+ _ V.

то 5 то 2п 5

Эта схема также реализуется с помощью явного вычислительного алгоритма.

4. Задача отражения и прохождения волны

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

Решение уравнений (1) при х<0 представляет собой суперпозицию падающей и отражённой волн:

V = / (г - х ) + В- ( г + С ), о = - г/ (г - С) + Вг/ \ г + х

где /(г)= /овтг - функция, определяющая форму падающей волны. Прих)0решение соответствует прошедшей волне: V = / \ г - Х |, ст = - г/ (г —) . В случае упругой

I с) с

прослойки при х = 0 выполняются уравнения (2):

р08(1 + А + В)) = 2г(1 - А = В), а08г(1 + А - В))- = 2(1 - А + В), йг йг

из которых следует, что

1 -аВ В-а /ш8 г^ ^ /ш8 г

А = 1-—ч, В = т—^-г-, где а =---, р =--.

(1 + а)(1 + Р) (1 + а)(1 + Р) 2с0 г 2с0 г0

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

В случае вязкоупругой прослойки решение может быть получено из приведённого выше решения заменой модуля упругой податливости а0 в выражении для коэффициента Р = — /ша0 г8 комплексным модулем, равным а0 +—— в модели Максвелла.

2 /шп

5. Уравнения блочной среды

Рассмотрим состояние плоской деформации блочного массива, образованного квадратными упругими блоками со сторонами длины И, параллельными координатным осям декартовой системы координат х1, х2 , и тонкими межблочными прослойками толщины 8. Блоки в массиве нумеруются парами индексов к1 и к2, пробегающих значения от1 до N и Ы2, соответственно. Внутри каждого блока выполняется система уравнений однородной изотропной упругой среды:

Р У1 =ст11,1 + ст12,2' Р У2 =ст12,1 + ст22,2 ' ст 12 = Р с2( У2,1 + у1,2) '

ст11 = Р с12 (у1Д + У2,2) - 2Р с2 У2,2 ' ст22 = Р с12 (у1Д + У2,2 ) - 2Р с2 У1Д •

Здесь с1 и с2 - скорости продольных и поперечных упругих волн, точка над символом и индексы после запятой означают производные по времени и по пространственным переменным, в остальном используются общепринятые обозначения.

Упругая прослойка между соседними блоками в горизонтальном направлении с номерами (к1, к2) и (к1 +1, к2) описывается обыкновенными дифференциальными уравнениями, учитывающими ее массу и продольную и поперечную жёсткости:

У+ + У'- =ст+1 -ст-1 ст+1 +ст-1 = ,,2 У1+- у-

Р -~-=-с-' -~-= Р с1

2 8 2 8

Р + У2 = ст+2 -ст-2 ст+2 + ст-2 = Р'с'2 У2 Р -^-=-с-' -^-= Р с2

2 8 2 8 Упругая прослойка между блоками в вертикальном направлении с номерами (к1, к2) и (к1, к2 +1) - уравнениями:

у2 + = ст+2-ст-2 ст+2 + ст-2 = Р'с'2 У2 -У2 Р -~-=-с-' -~-= Рс1

2 8 2 8

У+ + У'1 = а+2 _а12 а+2 + а 2 = р'с'2 Ъ+_ V

р -^-=-с-' -^-= р с2

2 5 2 5

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

Результаты расчётов показали, что под действием импульсной нагрузки образуются падающие волны (продольная и поперечная), конические поперечные волны, поверхностные волны, распространяющиеся со временем по всей области массива. Из-за наличия податливых прослоек возникает серия волн, отражённых от границ раздела, и серия поверхностных волн, образующихся при выходе отражённых волн на границу полуплоскости.

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

6. Заключение

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

Литература

1. Садовский М.А. Естественная кусковатость горной породы // Докл. АН СССР. 1979. 247. № 4. С. 829-831.

2. Курленя М.В., Опарин В.Н., Востриков В.И. О формировании упругих волновых пакетов при импульсном возбуждении блочных сред. Волны маятникового типа // Докл. АН СССР. 1993. 333. № 4. С. 3-13.

3. Александрова Н.И., Черников А.Г., Шер Е.Н. Экспериментальная проверка одномерной расчетной модели распространения волн в блочной среде // Физ.-техн. проблемы разработки полезных ископаемых. 2оо5. № 3. С. 46-55.

4. Александрова Н.И.,Шер Е.Н., Черников А.Г. Влияние вязкости прослоек на распространение низкочастотных маятниковых волн в блочных иерархических средах // Физ.-техн. проблемы разработки полезных ископаемых. 2оо8. № 3. С. 3-13.

5. Сарайкин В.А. Учет упругих свойств блоков в низкочастотной составляющей волны возмущений, распространяющейся в двумерной среде // Физ.-техн. проблемы разработки полезных ископаемых. 2оо9. № 3. С. 9-24.

6. Слепян Л.И. Нестационарные упругие волны. -Л.: Судостроение, 1972.

7. Бреховских Л.М. Волны в слоистых средах. -М.: Наука, 1973.

8. Кунин И.А. Теория упругих сред с микроструктурой. Нелокальная теория упругости. -М.: Наука, 1975.

9. Садовская О.В., Садовский В.М. Математическое моделирование в задачах механики сыпучих сред. -М.: Физматлит, 2668.

NumericalmodeHngofthepropagationofelasticwavesinstratifiedmediawithcompliantnterlayers

MariyaAlexandrovnaPokhabova, graduate student Institute of Computational Modeling SB RAS

Developed computer algorithms that simulate the dynamic interaction of elastic blocks through thin viscoelastic layer in structurally inhomogeneous media such as rocks. On the basis of rheological models are constructed for the deformation of the material layers of different levels of complexity. In the numerical solution of one-dimensional problems of applied monotonous grid- characteristic schemes with a balanced number of time steps in the blocks and interlayers . The problems of software implementation of algorithms on multiprocessor systems with graphics cards . Numerical results demonstrate the qualitative features of the propagation ofplane waves in materials with a layered microstructure.

Keywords: elastic waves, the block medium, microstructure, rheological scheme, the computational algorithm, mathematic modeling.

УДК 004.89

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

Галина Ахметовна Самигулина,д.т.н.

Тел.: 872722724617, e-mail:galinasamigulina@mail.ru Зарина ИльдусовнаСамигулина, PhD - докторант Тел.: 8 7272 2724617,e-mail:galinasamigulina@mail.ru лаборатория «Интеллектуальные системы управления и сети» Институт проблем информатики и управления МОН РК, http://wwwipic.kz

Работа посвящена проблеме построения оптимальной иммунносетевой модели для сложных инвестиционных проектов с целью прогнозирования рисков и возможности своевременной корректировки в процессе их реализации. Используется перспективный подход искусственных иммунных систем (ИИС), основанный на понятии формального пептида и процедуре молекулярного узнавания. Достоинствами применения ИИС при интеллектуальном анализе данных является возможность параллельной обработки большого количества мнений экспертов, прогнозирования рисков и своевременное управление ходом выполнения сложных проектов. Большое значение при реализации иммунносетевой технологии имеет процедура построения оптимальной модели на основе отбора наиболее информативных признаков. Для решения этой проблемы используется метод главных компонент.

Ключевые слова: интеллектуальные системы, прогнозирование рисков, искусственные иммунные системы, сложные инвестиционные проекты, метод главных компонент.

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

З.И. Самигулина