WWW.DISUS.RU

БЕСПЛАТНАЯ НАУЧНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА - Авторефераты, диссертации, методички

 

Pages:     | 1 |   ...   | 2 | 3 || 5 | 6 |   ...   | 7 |

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ (ИнЭРТ-2012) Труды X Международного научно-технического форума Ростов-на-Дону 2012 ИнЭРТ-2012 УДК 621:502:005.591.6:658 И 66 И 66 Инновация, экология и ...»

-- [ Страница 4 ] --

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

Список используемой литературы 1. Власов В.Г. Новый энциклопедический словарь изобразительного искусства.

В 10 т. – СПб.: Азбука-классика, 2004-2009.

2. Гапоненко А.Л. Теория управления. – М.: Изд-во РАГС, 2003.

3. Райзберг Б.А., Лозовский Л.Ш., Стародубцева Е.Б. Современный экономический словарь. - 5-е изд., перераб. и доп. – М.: ИНФРА-М, 2006. – 495 с.

4. Демин Ю. Делопроизводство. Подготовка служебных документов. – М.: «Питер Пресс», 2007.

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

УДК 681. Синтез системы управления реактором периодического действия In the report a complex problem of the creation of new control methods over chemical batch reactors (BR) is viewed. Those reactors represent non-linear multilinked objects. An ideology of general synergetic control theory is a basis of these methods. The core of new control methods comes to the synthesis objective control laws that take account of the natural chemical behavior of BR & the aspirations of the technological control task in the best way. At that the aims of the synthesized control systems are the desired invariants – attractors that reflect the technological aspirations. Chemical reactor; technological invariant; synergetic-chemical control theory; control law synthesis.

В докладе рассматривается задача синтеза синергетических законов управления [1, 2] типовыми химическими реакторами периодического действия (рис. 1). В реактор поступают потоки растворов некоторых веществ – компоненты G A и GB с начальными концентрациями C A 0 и C B 0.

С помощью мешалки производится интенсивное перемешивание реагентов, что позволяет усреднить концентрации по объему и использовать модель идеального смешения. В реакторе происходит необратимая химическая реакция соединения второго порядка типа, провоAB D вращением.

В процессе химического синтеза необходимо управлять подачей потока G A для приготовления раствора вещества в течение желаемого времени.

Математическая модель реактора имеет следующий вид [3]:

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

концентрацией x10 C A 0 ; – площадь сечения реактора.

Иначе говоря, реактор имеет один канал G A подачи реагента ( A ), а другой канал после загрузки компонента закрыт. Можно показать [3, 4], что начальная стадия протекания проB исходит процесс накопления подаваемого компонента A, что может привести к неуправляемому возрастанию его концентрации на завершающей стадии. Такое явление обычно недопустимо, т. к. оно ведет к излишнему загрязнению продукта исходным компонентом A. В этой свяD зи необходимо изменять поток G A компонента таким образом, чтобы на финишной стадии его концентрация уменьшалась вместе с убыванием концентрации ранее загруженного компонента B. Иначе говоря, требуется обеспечить управляемое плавное достижение желаемых конечных концентраций компонентов и и, в частности, их нулевых значений. Из уравнений (1) также следует, что при U 0, т. е. при прекращении подачи компонента A, в реакторе происходит неуправляемое естественное завершение химического процесса [3, 4].

Объект (1) включает управление во все уравнения и при u 0 обладает статической характеристикой в виде прямой независящей от управления. Это означает, что синтезируемый закон управления u ( x1, x 2, x3 ) должен обеспечивать асимптотически устойчивое удержание объекта (1) в одной из точек указанной характеристики.

Поставим задачу синтеза закона управления u ( x1, x 2, x3 ), обеспечивающего желаемое значение концентраций x1 x1s ; x 2 x 2s, в том числе и x1s x 2s 0. В соответствии с методом АКАР введем следующую макропеременную:

где 1, A x10. Тогда, используя функциональное уравнение

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

в силу уравнений объекта (1) получаем закон управления Этот закон обеспечивает неизбежный перевод объекта (1) сначала на многообразие (2), а затем движение вдоль него вплоть до попадания в заданное состояние x1 x1s ; x 2 x 2s.

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

представлены графики переходных процессов замкнутой системы (1), (3) и график изменения управления (3), подтверждающий эффективность синтезированных законов управления. Начальные условия переменных состояния: x1 ( 0) 15 моль м 3, x 2 ( 0 ) 1000 моль м 3.

Таким образом, метод АКАР позволяет осуществить аналитический синтез нового класса синергетических законов управления нелинейными химическими объектами.

Список используемой литературы Колесников А.А. Синергетическая теория управления. М.: Энергоатомиздат, 1994.

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

2. Современная прикладная теория управления. Ч. II. Синергетический подход в теории управления / Под ред. А.А. Колесникова. ФЦП «Интеграция». Москва-Таганрог: Изд-во ТРТУ, 2000.

Нейдорф Р.А. Характерные нелинейные инварианты и математические модели объектов управления химических производств // Современная прикладная теория управления.

Ч.III. Новые классы регуляторов технических систем / Под ред. А.А. Колесникова. ФЦП «Интеграция». Москва-Таганрог: Изд-во ТРТУ, 2000.

Ковалев П.А. Анализ задач управления химико-технологическими системами // Проблемы автоматизации и управление производственными системами. Ростов-на-Дону: Изд.

Центр ДГТУ, 1999.

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

УДК 681. The report provides a synergistic synthesis control system of injection engine with the analytical design of aggregated regulators (ADAR) [1]. The main engine fuel management tasks: getting the maximum torque on the shaft of the motor and the compliance with the stringent environmental requirements. To meet these challenges, according to the method of dekompozirovanna injection model is ADAR engine, which defines the controls u1, u2, providing maximum torque on the shaft and maintaining optimum air-fuel in the cylinders of the engine.

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

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

В докладе представлен синергетический синтез системы управления инжекторным двигателем с применением метода аналитического конструирования агрегированных регуляторов (далее АКАР) [1]. Рассмотрены основные задачи управления инжекторным двигателем: получение максимального крутящего момента на валу двигателя и соблюдение жестких экологических требований. Для решения этих задач согласно методу АКАР строится декомпозированная модель инжекторного двигателя, из которой определяются управления u1, u2, обеспечивающих максимальный крутящий момент на валу и поддержание оптимального соотношения воздухтопливо в цилиндрах двигателя.

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

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

Рассмотрим модель синтеза, включающую в себя дифференциальные уравнения, изложенные ниже. Динамика инжекторного двигателя описывается следующими выражениями [2]:

где – скорость вращения двигателя; t – коэффициент эффективности заполнения впускного коллектора, c – коэффициент наполнения цилиндров, P – давление воздуха во впускном колa лекторе до дроссельной заслонки, – универсальная газовая постоянная, Ta – температура возP духа до дроссельной заслонки, – отношение теплоемкости при постоянном давлении к теплоемкости при постоянном объеме, – диаметр дроссельной заслонки, 0 – угол дроссельной заD слонки при которой она полностью перекрывает отверстие, 1 – функция, определяющая поток воздуха через единицу площади, 2 – площадь открываемого дроссельной заслонкой отверстия, – угол поворота дроссельной заслонки, fi – поток топлива, впрыскиваемого форсункой, ff – поток массы топлива в топливной пленке, ign – угол опережения зажигания, J – приведенный к коленчатому валу момент инерции двигателя, M L – приведенный к коленчатому валу момент внешней нагрузки, – коэффициент, учитывающий, какая часть распыляемого топлива отношение воздух-топливо, c f 0, c f 1, c f 2 – коэффициенты, которые находятся экспериментально.

Нелинейные функции 1 ( P), 2 ( ), 3 ( ), 4 ( ign, ) определяются следующими выражениями [2]

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

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

Поэтому для нахождения закона управления введём макропеременные: 1 z1, 2 P 2, удовлетворяющие решению 1 2 0 однородных дифференциальных уравнений При попадании изображающей точки системы в окрестность инвариантных многообразий 1 0, 2 0 поведение системы будет описываться редуцированной моделью

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

Введем макропеременную 3 z 2 и функциональное уравнение d 3 3 3 0, из совместноdt го решения которых с учетом модели (8) определим выражение для внутреннего управления Законы управления u1 fi, u2 определяются из решения функциональных уравнений с учетом математической модели синергетического синтеза (6), выражения для которых здесь не приводятся ввиду их громоздкости. Результаты моделирования замкнутой системы приведены на рис. 1.

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

Список используемой литературы 1. Колесников А.А. Синергетическая теория управления. – М.: Энергоатомиздат, 1994.

2. Герасимов Д.Н., Джавахериан Х., Ефимов Д.В., Никифоров В.О. Инжекторный двигатель как объект управления. I. Схема двигателя и синтез математической модели // Известия РАН. ТиСУ, 2010, № 5, с. 135-147.

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

УДК 539.3:534. The problem of reconstruction properties of an inhomogeneous skin based on the standard model of a viscoelastic body with the multi-layer is considered. It is assumed that the layer has three components, simulating subcutaneous fat, dermis and epidermis. Two-dimensional problem is reduced to one-dimensional inverse problem whose solution is implemented with the involvement of the iterative process, every step of which deals with Fredholm integral equations of the first and second kind. Computational experiments are presented.

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

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

Приведены результаты численных экспериментов.

Рассмотрим установившиеся сдвиговые колебания неоднородного по толщине вязкоупругого слоя, занимающего область S x1, x2 (, ), x3 [0, h ]. Нижняя грань слоя S ром pe i t, где р ( р1, 0,0). Также предполагается, что вязкоупругий слой в свою очередь состоит из трех слоев, толщины которых известны.

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

Соответствующие уравнение и граничные условия после отделения временного множителя принимают вид:

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

Осредняя введенные в задаче функции по переменной x1 и проводя несложные преобразования, запишем задачу (1) в виде:

В данной работе рассматривала обратная задача реконструкции функций 2 ( x3 ) и U ( x3, ), удовлетворяющих (2), по дополнительной информации

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

n, получим вспомогательную задачу, для которой и построим операторные соотh ношения.

Дополнительное условие примет вид:

Заметим, что подобная задача рассматривалась ранее в [3], где для ее решения строился итерационный процесс, использующийся и в этой работе, однако в вычислительных экспериментах предполагалось, что материальные характеристики моделируются непрерывными функциями. В настоящем исследовании эта же задача рассмотрена в предположении, что характеристики имеют точки разрыва первого рода в известных точках. Обратная задача заключается в определении функций h( x ), g ( x), u ( x ) ( r (x) полагается известной положительной функцией) из краевой задачи (4) с учетом дополнительного условия (5). Полученные определяющие соотношения в итерационном процесс аналогичны полученным в [3], поэтому основное внимание уделено вычислительным экспериментам.

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

Параметры для моделирования соответствуют участку кожного покрова в районе предплечья. Мгновенный и длительный модули предполагались кусочно-непрерывными квадратичными функциями на [0, 1]:

Начальные приближения выбирались методом минимизации на конечномерной сетке.

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

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

Далее приведены результаты восстановления (рис. 1, 2):

Для восстановления частотный диапазон выбирался в между первым и вторым экстремумами амплитудно-частотной характеристики. Погрешность реконструкции не превосходит % в окрестностях точек разрыва и 4 % остальных точках.

Работа выполнена при частичной поддержке Российского фонда фундаментальных исследований (грант №10-01-00194-а), ФЦП "Научные и научно-педагогические кадры инновационной России" на 2009 – 2013 годы (госконтракт П596).

Список используемой литературы 1. Балабанов Е. И Аналитический обзор. Кожа человека, механические свойства, теплопередача. М.: 2001.

2. Ватульян А.О. Обратные задачи в механике деформируемого твердого тела. М.:

Физматлит, 2007. 223 с.

3. Ватульян А. О., Явруян О. В., Богачев И. В. Идентификация упругих характеристик неоднородного по толщине слоя // Акустический журнал. Т. 57. № 6. 2011. С. 723-730.

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

УДК 621.391. Методика рационального проектирования конструкций погрузочного оборудования дорожных машин семейства «АМКОДОР» на основе конечно-элементного моделирования А. А. Боровик, В. В. Напрасников, Камран Каземпур Абдолреза Беларусь, БНТУ, [email protected], [email protected] A design procedure of the road-making machine loading equipment is considered on the base of finite-element modeling.

Создание новой техники в современном машиностроении проходит в несколько этапов:

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

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

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

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

Последовательность создания рациональной конструкции состоит из следующих этапов:

- разработка полной трехмерной твердотельной модели рабочего оборудования погрузчика. Результаты выполнения этого этапа представлены на рисунке 2;

- определение элементов конструкции, не оказывающих существенного влияния на распределение напряжений в материале и соответствующее упрощение модели. Такое упрощение позволяет существенно сократить размерность конечно-элементной модели, создаваемой в дальнейшем. В данной работе такое упрощение выполнено вручную. В настоящее время ведутТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

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

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

- выполнение тестового расчета для проверки адекватности созданной модели, на примере одного расчетного случая;

Рис. 2. Вид рабочего оборудования погрузчика без упрощений Рис. 3. Вид рабочего оборудования погрузчика после упрощения геометрии

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

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

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

- определение воздействия внешней по отношению к объекту среды на выбранную конструкцию;

- выполнение прочностных расчетов;

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

Рис. 4. Результаты расчетов двух вариантов конструкций

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

Рис. 5. Распределение напряжение по Мизесу в конструкции рабочего оборудования погрузчика В ходе реализации данной методики на примере дорожной машины «АМКОДОР – 320», в выбранной конструкции были выявлены места, где напряжения превышают допустимые. На основании полученных данных, была проведена работа по усовершенствованию существующей конструкции путем внесения следующих изменений:

- для увеличения прочности верхней части ковша изменены размеры верхней коробки;

- удалены внешние горизонтальные рёбра жёсткости;

- изменена геометрия проушин навески ковша;

- добавлены ребра жесткости внутри коробки ковша под проушины.

На рисунке 4 приведены результаты расчета (в виде диаграмм напряжений) первоначальной конструкции и конструкции после внесения изменений.

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

Рис. 6. Фотография эксплуатационных разрушений рабочего оборудования погрузчика В результате, внесенные изменения позволили снизить напряжения в опасных участках конструкции с 346МПА до 173МПА.

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

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

Данная методика внедрена в конструкторскую практику на предприятии «АМКОДОР» в Республике Беларусь. Отдельные ее этапы продолжают совершенствоваться.

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

Список используемой литературы 1. Боровик А. А. Особенности построения и подготовки конечно-элементных моделей сложных транспортных конструкций. 10-я Международная НТК «Наука – образованию, производству, экономике», Минск, 2012, 218 c.

2. Боровик А.А., Напрасников В.В., Камран Каземпур Абдолреза. Выбор рационального варианта системы погрузочного оборудования дорожной машины АМКОДОР - 320 на основе конечно-элементного анализа. 10-я Международная НТК «Наука – образованию, производству, экономике», Минск, 2012, 226 c.

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

УДК 539. A number of structural and composite materials, for which cracks are the characteristic defects, exhibits anisotropic and rheological properties. In this paper the problem of the identification of a single crack in the viscoelastic orthotropic layer by the known information about the measured at the part of the upper bound fields of displacements is considered.

The boundary value problem is reduced to a system of boundary integral equations, which is solved by the method of boundary elements. The viscoelastic properties of the material are considered at the frames of the correspondence principle.

Problem of the identification of the crack parameters is reduced to minimization the residual functional. The numerical implementation is effected by the genetic algorithm. Asymptotic analysis of the problem under the assumption that the characteristic crack size is small compared with the thickness of the layer is made.

Рассмотрим установившиеся колебания вязкоупругой ортотропной полосы S толщины h, нижняя грань которой x 2 0 лежит без трения на жестком основании. Верхняя граница свободна от напряжений за исключением области приложения нагрузки.

Вязкоупругие свойства учитываются на основе принципа соответствия [1]. В рамках концепции комплексных модулей физико-механические свойства вязкоупругого материала задаются комплексными функциями, зависящими от частоты колебаний.

Слой ослаблен внутренней трещиной. Трещина моделируется как математический разрез в области S с берегами l0, на которых компоненты вектора смещений терпят разрыв. Применяя теорию дислокаций, действие трещины заменяется действием фиктивных массовых сил f k ( k 1,2 ), представляющих собой конечные скачки компонент полей перемещений на трещине.

Считаем, что во время колебаний берега трещины не взаимодействуют друг с другом и свободны от напряжений.

После отделения временного множителя краевая задача принимает вид:

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

где kl – компоненты тензора напряжений, C klmn (i ) - компоненты тензора комплексных модулей, f k [C klmn (i ) nm n ( )], l дельта-функция Дирака, – координата, отсчитываемая фиктивные массовые силы, по нормали к поверхности l0, nl – компоненты единичных векторов внешних нормалей к поверхностям l0.

Задача идентификации трещины состоит в определении геометрии и местоположения дефекта по дополнительной информации («экспериментальным» данным) – полям перемещений, измеренным на части верхней границы полосы l2 :

На основе формул Сомильяны и представлений функций Грина в виде контурных интегралов получены интегральные представления полей смещений в виде [2]:

без дефекта вызванное действием заданной нагрузки, второе слагаемое (отраженное поле) обуm) словлено наличием трещины и определяется через функции раскрытия трещины, U i ( m ) ( x, ) – функции Грина и соответствующие им напряжения, вычисляемые по определяюij щим соотношениям с комплексными модулями типа обобщенного закона Гука и удовлетворяющие дополнительным граничным условиям на верхней и нижней границе слоя

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

Для определения компонент функции раскрытия трещины осуществлен предельный переход в (7) при y l0, далее, удовлетворяя условию отсутствия напряжений на трещине (5), получена система ГИУ относительно скачков смещений на трещине:

Главные части ядер интегральных операторов K ji ( x, y ), так же как и в упругом случае, являются гиперсингулярными, а соответствующие интегралы понимаются в смысле конечного значения по Адамару [3].

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

На основе решения прямой задачи решена задача идентификация трещины. Система нелинейных операторных уравнений относительно неизвестных i (x ) и l0, связывающих геометрию трещины с полем смещений на верхней грани, состоит из уравнений прямой задачи вида (10) и дополнительного интегрального соотношения вида:

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

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

В качестве входной информации заданы значения компонент полей перемещений g m ( i ) в M точках i l2 на частоте. Функционал невязки имеет вид:

где qk - параметры трещины простейшей конфигурации, в случае прямолинейной трещины Численный анализ проведен для ортотропного композита, состоящего из эпоксидной смолы с 60-процентным содержанием продольных волокон графита [4]. Минимизация функционала (12) осуществлена при помощи генетического алгоритма [5].

Относительные погрешности восстановления инвариантных параметров прямолинейной трещины не превышают 3-4 % при точных входных данных. При зашумлении входной информации с амплитудой в 2 % относительные погрешности восстановления редко превышают 5Для обеспечения точности, соответствующей порядку функционала невязки 0.001, количество вычислений функционала невязки, а, значит, и обращений к прямой задаче, составляет в среднем от 1000 до 2000 раз. Численный анализ показал, что для эффективной работы алгоритма требуется 8-12 точек съема данных, увеличение их числа не сказывается на точности восстановления параметров трещины.

Так как при использовании генетического алгоритма требуется многократное решение прямой задачи, которое достаточно сложно и связано с затратой машинного времени, то при малом характерном размере дефекта можно использовать асимптотические разложения по малому параметру. В этом случае удается значительно упростить вычисление полей перемещений, решение системы интегральных уравнений и свести задачу идентификации к поэтапному определению параметров трещины из трансцендентных уравнений [6]. Анализ вычислительных экспериментов с использованием асимптотического подхода показал, что три параметра, d, L при l 0.1h определяются с относительной погрешностью менее 5-7 %. Параметр l в силу того, что он определяется в последнюю очередь и на его определении оказывают влияние погрешности определения предыдущих параметров, определяется с погрешностью не превышающей 10-12 %. При этом время вычислений сокращается примерно на два порядка по сравнению с методом граничных интегральных уравнений.

Работа выполнена при поддержке ФЦП «Научные и научно-педагогические кадры инновационной России» на 2009-2013 годы (госконтракт № П596) и гранта РФФИ 10-01-00194-а.

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

Список используемой литературы 1. Кристенсен Р. Введение в теорию вязкоупругости. М.: Мир. 1974. 338 с.

2. Ватульян А.О., Азарова П.А., Явруян О.В. Идентификация параметров наклонной прямолинейной трещины. Механика композиционных материалов и конструкций // Механика композиционных материалов и конструкций. 2008. Т. 14, № 3 с.461-472.

3. Белоцерковский С.М., Лифанов И.К. Численные методы в сингулярных интегральных уравнениях. М.: Наука. 1985. 256 с.

4. Garnich M.R., Hansen A.C. A multicontinuum Approach to Structural Analysis of Linear Viscoelastic Composite Materials // J. of Applied Mechanics. December 1997. Vol. 64. P. 795Баранов И.В., Ватульян А.О., Соловьев А.Н. Об одном генетическом алгоритме и его применении в обратных задачах идентификации упругих сред // Вычислительные технологии. 2006. Т.11. №3. С.14-26.

6. Ватульян А. О. Обратные задачи в механике деформируемого твердого тела. М.: Физматлит. 2007. 223 с.

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

УДК 539.3:534. Об акустических методах определения неоднородных свойств упругих тел In this paper the inverse coefficient problem based on the findings obtained by acoustic sensing. A method of solving the inverse problem on the basis of difference approximations.

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

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

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

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

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

не вносят искажения в происходящие в материале процессы и в то же время несут информацию о состоянии исследуемого объекта.

1. Постановка задачи. Представлена классическая краевая смешанная задача для неоднородного эллиптического оператора типа Гельмгольца с переменными коэффициентами описывающая крутильные колебания стрежня с поперечным сечением S и переменным модулем сдвига [2].

Рассмотрим задачу решения дифференциального уравнения (1.1) с переменными коэффициентами в области S [0,1] [c, c] со следующими граничными условиями Теория эллиптических уравнения по сей день не является завершенной. К сожалению не все задачи этого класса допускают аналитическое решение, поэтому решение такого рода задач можно построить лишь приближенно на основе каких либо численных методов. В данной работе предложен метод разностных аппроксимаций.

2. Исследование прямой задачи. Прямая задача состоит в нахождении функции u x1, x2 при заданном законе изменения g x1, x 2 и удовлетворяющей краевым условиям (1.2) в некотором диапазоне изменения.

В методе разностных аппроксимаций применен пятиточечный шаблон для вторых производных и двухточечный для первых производных (левая разностная производная) [3]. После введения прямоугольной сетки с равным шагом по координатным осям x1, x2 была получена система линейных алгебраических уравнений с учетом граничных условий, которые позволяют исключить из этой системы граничные узловые точки. Таким образом полученная система имеет решение, причем единственное. Решая систему 2.1, находим узловые значения функции ui, j, которые при решении обратной задачи послужат как входная информация.

3. Исследование обратной задачи. Целью обратной задачи является восстановление неизвестной функции g x1, x2 по известным узловым значениям функции u x1, x2, полученным в результате решения прямой задачи методом разностных аппроксимаций, описанным в разделе 2. Основной проблемой на этом пути является некорректная задача вычисления производной от функции, заданной таблично.

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

Для оценки точности разностного подхода введем в рассмотрение относительную погрешность где gТ i, j — точное решение, g П i, j — полученное.

Для аналогии с реальными экспериментами функция u x1, x2, используемая как входная информация для восстановления g x1, x2, была зашумлена [4].

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

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

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

Список используемой литературы 1. Ватульян А. О., Денина О. В. Обратные коэффициентные задачи для стержней.

Методы определения неоднородных свойств упругих стержней на основе акустического зондирования. LAP LAMBERT Academic Publishing GmbH & Co. KG, 2011. 104 с.

2. Ватульян А.О. К теории обратных коэффициентных задач в линейной механике деформируемого тела. ПММ. – 2010. – №6. – С. 911-918.

3. Самарский А.А. Введение в теорию разностных схем. Москва: Наука, 1971. – 552 с.

4. Ватульян А.О., Бурьян А.Ю., Осипов А.В. Об идентификации переменной жесткости при анализе поперечных колебаний балки. Вестник Донского гос. техн. ун-та. – 2010. – Т. 10. – № 6. – С. 825-833.

5. Альберг Дж. Теория сплайнов и ее приложения. Москва: Мир. 1972. – 316 с.

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

УДК 621. In the article the mathematical model of the Ripper with battery energy. The analysis of the influence of elastic connection, geometrical parameters and the mass of the transmission coefficient. Compiled graphical relationship that establishes the relationship between the frequency of the kinematic excitation and transmission coefficient.

Из предложенного класса машин, объединенных принципом замыкания динамических нагрузок на рабочем оборудовании и грунте [4—20], рассмотрим возможность управления динамическим состоянием механической системы «базовая машина – аккумулятор энергии – рабочий орган – грунт» с целью перераспределения энергии колебаний в зону разрушения грунтового массива.

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

Уравнения вынужденных колебаний масс m1 и m2, если пренебречь демпфированием в элементах трансмиссии и гусеничного движителя, при кинематическом возбуждении системы имеют вид:

Если кинематическое возбуждение описывается гармонической функцией

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

то колебания масс m1 и m2 можно записать, используя гармонические законы Подставляя выражения (2) и (3) в (1) приведем последние к алгебраическому виду Используя полученные выражения, вычислим передаточные функции Для перехода к безразмерным величинам, поделим числитель и знаменатель подкоренного выражения уравнений (5) и (6) на c1c3. Полагая, что c 2 c 3 ; 2 1 K ; 2 2 с 2 2 ;

с 2 с1 K c ; m2 m1 K m, представим слагаемые уравнений в следующем виде:

В этом случае коэффициенты передачи 1 и 2 примут вид:

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

Сопоставляя графические зависимости коэффициентов передачи основного объекта 1 и колеблющейся массы 2 (рис. 2), можно сказать, что для различных частот кинематического возбуждения на графике существуют зоны соответствующие оптимальному значению K 1 2. В этих зонах коэффициент K стремится к своему минимальному значению. Например, при 2 0.82, K c 0,5 среднее значение K 0,038 0,833 0,043.

С увеличением 2 значение коэффициента K уменьшается. Например, в случае когда 2 0,182, Kc 0.1 K 0,098. Величина коэффициента свойств колеблющейся массы, а также от сил вязкого сопротивления, возникающих в зоне контакта колеблющейся массы с грунтовым массивом. В рассматриваемом случае величина диссипативных сил была принята постоянной и равной 2 20кН с / м.

Рис. 2. Зависимость коэффициента передачи от частоты кинематического Частота кинематического возбуждения связана со скоростным режимом движения основного объекта и физико-механическими свойствами среды в зоне взаимодействия с колеблющейся массой. Изменение указанных параметров влечет за собой рассогласование между коэффициентами передачи 1 и 2 в сторону ухудшения динамического состояния основного объекта. Адаптивное управление упруго-инерционными свойствами колеблющейся массы в ручном или автоматическом режиме позволит целенаправленно приводить механическую сисТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

тему к рациональному состоянию.

Список используемой литературы 1. Лойцянский Л.Г. Курс теоретической механики. Т. 2. Динамика / Л. Г. Лойцянский, А. И. Лурье. – М.: Наука, 1983. – 640 с.

2. Коренев Б.Г. Динамические гасители колебаний. Теория и технические приложения / Б.Г. Коренев, П.М. Резников. – М.: Наука, 1963. – 535 с.

3. Елисеев С.В. Динамический синтез в обобщенных задачах виброзащиты и виброизоляции технических объектов / С.В. Елисеев, Ю.Н. Резник, А.П. Хоменко, А.А. Засядько. – Иркутск: Изд-во Иркут. гос. ун-та, 2008. – 523 с.

4. Геллер Ю.А. Создание эффективной техники на примере машин для специальных земляных работ, действующих по принципу замыкания динамических нагрузок на рыхлительном оборудовании и грунте: монография. – Чита: ЗабГУ, 2011. – 217 с.

5. Геллер Ю.А. Способ механической разработки грунтов; Заявка на предполагаемое изобретение. RU 2008116379. Приоритет от 05.11. 2008. опубликована 27.10.2009. – Бюл. № 30.

6. А.с. 815169 СССР МКИ3 E 02 F 5/30. Рыхлитель / Н.П. Безручко, Ю.А. Геллер (СССР). – № 2727234/29-03; заявл. 22.02.79; опубл. 23.03.81, Бюл. № 11. – 2 с.: ил.

7. А.с. №889805 СССР, МКИ3 E02F 5/30. Рыхлитель / Н.П. Безручко, Ю.А. Геллер, А.А. Киричек, В.П. Козлов, А.П. Гаршин (СССР). –№2892666/29-03; заявл. 07.03.80; опубл.

15.12.81, Бюл. № 46. – 3 с.: ил.

8. А.с. 939672 СССР МКИ3 E 02 F 5/30. Рыхлитель / Н.П. Безручко, Ю.А. Геллер, А.А. Киричек (СССР). – № 3222893/29-03; заявл. 24.12.80; опубл. 30.06.82, Бюл. № 24. – 3 с.: ил.

9. А.с. 968558 СССР, МКИ3 E 02 F 5/30. Устройство для разработки прочных грунтов / Н.П. Безручко, Ю.А. Геллер, А.А. Киричек (СССР). – № 3272533/29-03; заявл. 03.04.81; опубл.

23.10.82, Бюл. № 39. – 4с.: ил.

10. А.с. 994650 СССР, МКИ3 E 02 F 5/30. Рыхлитель для разработки мерзлых и прочных грунтов / Н.П. Безручко, Ю.А. Геллер, А.А. Киричек (СССР). – № 2892665/29-03; заявл.

07.03.80; опубл. 07.02.83, Бюл. № 5. – 4 с.: ил.

11. А.с. 1016445 СССР, MКИ3 E 02 F 5/30. Рыхлитель / Ю.А. Геллер, А.А. Киричек, Н.П. Безручко, Г.Р. Круглов (СССР). – № 3399226/29-03; заявл. 24.02.82; опубл. 07.05.83, Бюл.

№ 17. – 4 с.: ил.

12. А.с. 1304465 (РФ), MКИ3 E 02 F 5/30. Рыхлитель / Ю.А. Геллер, А. А. Киричек, Н.Е. Курбатов, Е.П. Маккавеев (РФ). – № 3700504/29-03; заявл. 10.02.84; опубл. 15.12.86, Бюл.

№ 33. – 4 с.: ил.

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

13. Пат. 1176944 (РФ), MКИ3 E 02 F 5/30. Рыхлитель / Ю.А.Геллер (РФ); заявитель и патентообладатель Чит. гос. ун-т. – № 3709935/29-03; Заяв. 02.01.84; Опубл. 07.09.85. – Бюл.

№ 33. – 4 с: ил.

14. Пат. 2222669 (РФ), МКИ 7 Е 02 F 5/30. Вибрационный рыхлитель / Ю.А. Геллер (РФ); заявитель и патентообладатель Чит. гос. ун-т. – № 2001114130; Заяв. 23.05.01; Опубл.

27.01.04. – Бюл. № 3. – 5 с.: ил.

15. Пат. 2239689 (РФ), MКИ 7 E 02 F 3/00, G01M 19/00. Стенд для исследования рабочих органов землеройных машин / Ю.А. Геллер (РФ); заявитель и патентообладатель Чит. гос.

ун-т. – № 2002122136/03; Заяв. 13.08.02; Опубл. 10.11.04. - Бюл. № 31. – 9 с.: ил.

16. Пат. 2367747 (РФ), МКИ 7 Е 02 F 5/30. Вибрационный рыхлитель / Ю.А. Геллер (РФ); заявитель и патентообладатель Чит. гос. ун-т. - № 2008116382; Заяв. 24.04.08; Опубл.

20.09.09. - Бюл. № 26. – 7 с.: ил.

17. Пат. 2372447 (РФ), МКИ 7 Е 02 F 5/30. Рыхлитель ударного действия / Ю.А. Геллер (РФ); заявитель и патентообладатель Чит. гос. ун-т. – № 2008120282; Заяв. от 21.05.08; Опубл.

10.11.09. – Бюл. № 31. – 5 с.: ил.

18. Пат. 2380489 (РФ), МКИ 7 Е 02 F 5/30. Рыхлитель ударного действия / Ю.А. Геллер (РФ); заявитель и патентообладатель Чит. гос. ун-т. – № 2008116381; Заяв. 24.04.08; Опубл.

27.01.10. – Бюл. № 3. – 4 с.

19. Геллер Ю.А. Рыхлитель с пружинным аккумулятором энергии двухстороннего действия / Ю.А. Геллер; Положительное решение на выдачу патента по заявке RU 2010146238.

Приоритет от 12.11.2010.

20. Геллер Ю.А. Рыхлитель с жидкостным аккумулятором энергии двухстороннего действия / Ю.А. Геллер; Заявка на предполагаемое изобретение. RU 2012145318. Приоритет от 30.05.2012.

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

УДК 539. Компьютерный анализ прочности поврежденной костной ткани в модельной системе остеосинтеза с аппаратом Илизарова при статических и динамических воздействиях The finite element analysis of the damage bone tissue for model osteosynthesis system with Ilizarov external fixator was realized under static and dynamic loads.

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

Так, в [1] имеется обзор подобных исследований, в [2, 3] с использованием комплекса NASTRAN проведен линейный и нелинейный статический анализ имитатора кости с одним модулем и полной сборкой аппарата Илизарова. В [1, 4] МКЭ на основе программного кода MechanicsFE применялся для анализа упругих и пороупругих свойств поврежденной костной ткани при колебаниях кости с аппаратом Илизарова.

В настоящей работе объектом исследования являлась биомеханическая модель, состоящая из имитатора кости и аппарата Илизарова, а основная цель состояла в определении напряженно-деформированного состояния (НДС) костной ткани в месте перелома на различных этапах заживления и при варьировании внешних воздействий, связанных с осевой нагрузкой на кость и компрессией или дистракцией соединительных стрежней между репонирующими кольцами. При исследовании рассматриваемой системы использовалась техника метода конечных элементов, в программном комплексе ANSYS строились пространственные конечноТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

элементные модели и решались соответствующие статические и динамические задачи теории упругости.

Для упрощения анализа кость моделировалась в форме цилиндра, подразделенного на части для задания различных механических свойств в неповрежденной компактной и спонгиозной костных тканях и в месте перелома. Аппарат Илизарова включал наборы внешних опор в форме колец, соединительных стержней и cпиц, служащих связующими звеньями между костью и жесткой конструкцией аппарата. Построенная твердотельная модель представлена на рис. 1 в форме опорных линий и кривых. (При этом были построены дополнительные линии для удобства задания мест соединения различных объектов.) Для дальнейшего формирования экономичной конечно-элементной сетки в геометрической модели для имитатора кости с переломом использовались цилиндрические объемы, для базовых и репонирующих колец аппарата Илизарова – кольцевые поверхности, а для стержней и спиц – линии.

Рис. 1. Твердотельная модель (линии) Рис. 2. Конечно-элементная модель Модель имела следующие геометрические размеры: длина имитатора кости – 450 мм;

длина зоны перелома – 10 мм; внешний радиус имитатора кости – 15 мм; внутренний радиус для спонгиозной ткани – 7.5 мм; внешний и внутренний радиусы соединительных колец – 70 и 55 мм, соответственно; толщина соединительных колец – 5 мм; длина соединительных стержней между базовыми и репонирующими кольцами – 150 мм; длина соединительных стержней между репонирующими кольцами в зоне перелома – 100 мм; радиус соединительных стержней – 2.5 мм; длина спиц – 80 мм; радиус спиц – 0.75 мм.

Для достижения оптимального соотношения между точностью конечно-элементных расчетов и относительным минимумом вычислительных затрат указанные части конструкции моделировались в ANSYS соответствующими трехмерными структурными, оболочечными и баТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

лочными конечными элементами SOLID92, SHELL93 и BEAM189, соответственно. Такой подход представляется вполне оправданным с учетом большей жесткости аппарата в сравнении с костной тканью и тем, что здесь анализ НДС в зонах контакта костной ткани со спицами не являлся объектом исследования. Построенная конечно-элементная модель (рис. 2) позволяет задавать также ортотропные упругие свойства материалов имитатора кости (областей компактной и спонгиозной тканей), вводить неоднородные механические свойства регенерирующей ткани вблизи места перелома, изменять базовые геометрические и механические характеристики модели и задавать различные нагрузки на имитатор кости. Важным моментом исследования являлся учет посредством уравнений связей деформаций стержней с целью осуществления управляемого остеосинтеза с требуемой компрессией или дистракцией.

При статическом анализе нижняя торцевая часть имитатора кости считалась жестко закрепленной. Рассматривались следующие варианты внешних воздействий: 1) задана компрессия или дистракция соединительных стрежней между репонирующими кольцами в зоне перелома на 1 мм; 2) к верхнему торцу приложено давление 1 МПа; 3) одновременно задана компрессия на 1 мм и приложено давление 1 МПа; 4) одновременно задана дистракция на 1 мм и приложено давление 1 МПа.

Стержни, кольца и спицы аппарата Илизарова принимались изотропными материалами с модулем Юнга E 2.23 105 (МПа) и коэффициентом Пуассона 0.3 (сталь). В здоровой кости спонгиозная ткань считалась изотропным материалом с E 500 (МПа) и 0.25, а корковая (компактная) часть – ортотропным материалом со следующими параметрами: E x 8500, xz 0.12. В зоне перелома задавались различные варианты регенерирующей ткани: 1) мягкая гелеобразная ткань, как изотропный материал с E 0. как изотропный материал с E 3.78 (МПа) и 0.45 ; 3) спонгиозная ткань здоровой кости; 4) спонгиозная ткань здоровой кости во внутренней цилиндрической части и корковая ткань на этапе ее формирования во внешней цилиндрической части, как ортотропный материал с xz 0.12 ; 5) спонгиозная ткань здоровой кости во внутренней цилиндрической части и корковая ткань здоровой кости во внешней цилиндрической части.

Ниже представлены таблицы, в которых для каждой, заполняющей место перелома ткани, и для различных вариантов нагрузок приведены максимальные значения напряжений по Мизесу и соответствующие пределы прочности u, взятые из [5]. Данные таблицы позволяют оценить допустимые величины внешних нагрузок на кость и деформаций стержней аппарата

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

Илизарова в процессе регенерации кости, исходя из критериев допуска на максимальные характеристики напряжений в костной мозоли.

Напряжения по Мизесу во внешней части места перелома Характеристики Напряжения по Мизесу (MPa) при различных воздействиях u (MPa) места перелома Компрессия Давление Компрессия Дистракция и корковой части Напряжения по Мизесу во внутренней части места перелома Характеристики Напряжения по Мизесу (MPa) при различных воздействиях u (MPa) места перелома Компрессия Давление Компрессия Дистракция и корковой части В результате решения динамических задач для ряда наборов значений механических свойств костной мозоли в зависимости от периодов регенерации были определены также первые резонансные частоты системы, построены амплитудно-частотные характеристики в характерных точках кости и аппарата и проведен анализ переходных процессов при заданных динамических внешних воздействиях и условиях закрепления. Более полные данные по полученным результатам приводятся в докладе, а здесь опущены в силу ограниченности объема статьи.

Список используемой литературы 1. Математические модели и компьютерное моделирование в биомеханике / Под ред. А.В. Зиньковского и В.А. Пальмова. СПб.: Изд-во Политехн. ун-та, 2004. Гл. 10. Л.Б. Маслов, Д.В. Ликсонов. Биомеханические характеристики нижней конечности человека. С. 385-438.

2. Бушманов А.В., Барабаш С.А. Численное моделирование деформации спиц циркулярного аппарата // Медицинская информатика. 2007. № 1 (13). С. 33–40.

3. Бушманов А.В., Соловцова Л.А. Исследование жесткости аппарата Илизарова // Российский ж. биомеханики. 2008. Т. 12, № 3 (41). С. 97–102.

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

4. Маслов Л.Б. Математическое моделирование колебаний пороупругих систем.

Иваново: Изд-во ИГЭУ, 2010. 264 с.

5. Образцов И.Ф., Адамович И.С., Барер А.С. и др. Проблемы прочности в биомеханике. М., 1988. 311 с.

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

УДК 539. Контактная задача для трансверсально изотропного полупространства The spatial contact problem for a transversely isotropic half-space is investigated for an elliptic punch. The isotropy planes are perpendicular to the half-space boundary. The integral equation of the contact problem is derived. An exact solution is obtained for the case when the punch is an elliptic paraboloid and the contact zone is given. For an elliptic punch with a polynomial base the exact solution structure is established.

В декартовых координатах рассмотрим трансверсально изотропное упругое полупространство x0, граница которого перпендикулярна плоскостям изотропии zconst. Закон Гука имеет вид [1] Пусть в начале координат к границе полупространства приложена нормальная сосредоточенная сила Px. Граничные условия задачи Буссинеска для уравнений упругого равновесия можно записать в виде где (x) дельта-функция Дирака.

Подставив закон Гука (1) в уравнения равновесия в напряжениях, перейдем к уравнениям равновесия в перемещениях, решение которых будем искать в форме двойного преобразования Фурье по переменным y, z. В результате нормальное перемещение при x 0 найдем в виде [1]

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

Здесь 1, 2 являются положительными корнями уравнения (предполагается, что такие корни существуют) Рассмотрим контактную задачу о вдавливании эллиптического в плане штампа в трансверсально изотропное полупространство. Пусть основание штампа описывается функцией Штамп вдавливается без перекоса центрально приложенной силой P, испытывая осадку. Штамп имеет острые кромки, поэтому эллиптическая область контакта считается известной в виде При заданной функции (5), области контакта (6) и осадке требуется определить нормальное контактное давление q(y,z) в области и силу P. Рассмотрим случай, когда эллипс контакта вытянут вдоль оси z, т. е. a b (случай, когда эллипс контакта вытянут вдоль оси y, рассматривается аналогично). Вводя безразмерные обозначения используя решение задачи Буссинеска (3), интегральное уравнение контактной задачи можно записать в форме (штрихи далее опускаем) Решение уравнения (8), следуя идеям Л. А. Галина [2], будем искать в форме

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

Внесем (9) в (8). С учетом значений интегралов 3.715.10 и 6.554.2 [3] и перехода к полярным координатам интегральное уравнение контактной задачи запишется в виде ( ( y, z ) S ) Обозначим z cos k 1 y sin. Ясно, что внутри эллипса контакта | | 1. Используя интегралы 3.741.2, 3.741.3 и 3.784.7 [3], выполним квадрату по переменной в (11). Тогда придем к соотношению ( ( y, z ) S ) Заметим, что интеграл, содержащий zy sin 2, равен нулю. Приравнивая в (13) слева и справа члены при одинаковых степенях переменных y, z и свободные члены, придем к следующей системе трех линейных алгебраических уравнений для определения постоянных An (n0,1,2):

Коэффициенты системы (14) имеют вид

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

Решение системы (14), (15) не вызывает труда. Вдавливающая штамп сила находится затем по формуле Для эллиптического в плане штампа с полиномиальным основанием имеет место следующая теорема, обобщающая результат Л.А. Галина для изотропного случая [2].

Теорема. Решение интегрального уравнения (оно отличается от (15) лишь правой частью) где Rm(y,z) заданный полином степени m, имеет следующую структуру:

где Qm(y,z) полином степени m с неопределенными коэффициентами, которые можно найти после подстановки (18) в (17), взятия интегралов, тождественно удовлетворяя уравнению (17).

Исследование поддержано грантом РФФИ 12-01-00065.

Список используемой литературы 1. Fabrikant V.I. Non-traditional contact problem for transversely isotropic half-space // Quarterly Journal of Mechanics and Applied Mathematics. 2011. Vol. 64. No. 2. P. 151–170.

2. Галин Л.А. Контактные задачи теории упругости и вязкоупругости. — М.: Наука, 1980. — 303 с.

3. Градштейн И.С., Рыжик И.М. Таблицы интегралов, сумм, рядов и произведений. М.:

Наука, 1971. 1108 с.

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

УДК 539. Конечно-элементный дизайн ультразвуковых пьезопреобразователей из пористой Россия, ДГТУ, ЮНЦ РАН*, [email protected], [email protected] The methodology for modeling of ultrasound devices from porous piezoceramics was presented. As an example the hydroacoustic emitters made of solid ceramics with intermediate interjacent non-piezoelectric layers and also piezoelectric emitter from porous piezoceramics without such layers were considered.

В современное время пьезоустройства с композитными пьезоэлементами становятся все более привлекательными для разнообразных ультразвуковых применений. Основная выгода при использовании пьезокомпозитов заключается в возможности существенного улучшения основных рабочих параметров пьезоустройств, особенно, при их нагрузке на акустические среды. Так, экспериментальные исследования [1, 2] показывают, что пористая пьезокерамика обладает высокой объемной пьезочувствительностью в широкой полосе частот и более низким импедансом по сравнению с плотной пьезокерамикой. При этом указанные эффекты усиливаются при повышении степени пористости керамики, но важнейшие толщинные характеристики (пьезомодуль d 33, коэффициенты электромеханической связи k t, k33 ) для ряда пористых пьезокерамик практически не зависят от пористости, а соответствующие продольные величины d 31, k p, k31 ) быстро убывают с ростом пористости [1-3]. С целью дальнейшего анализа эффективности пористой пьезокерамики для гидроакустических применений в настоящей работе проведено компьютерное моделирование пьезоизлучателя из пористой пьезокерамики с согласующими упругими слоями, а также, пьезоизлучателя без переходных слоев.

На первом этапе исследования проводился расчет эффективных модулей пористых пьезокомпозитов различной связности с использованием метода эффективных модулей и конечноэлементного (КЭ) моделирования представительных объемов пористой пьезокерамики [3-5].

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

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

и электрической индукции однородного тела. Расчеты проводились с помощью КЭ пакета ANSYS и специально разработанных компьютерных программ, написанных на макроязыке APDL ANSYS. При этом использовались разработанные В.В. Ремизовым программы, реализующие следующие методы генерации представительных объемов двухфазного пьезокомпозита в виде куба с кубическими элементами: случайный, метод начальной концентрации, метод ОДА (ограниченная диффузией агрегация) Виттена-Сандера [6]. На примере пьезокерамики ПКР-8 в [5] было проведено сравнение вычисленных эффективных характеристик с экспериментальными данными, полученными в НИИ физики ЮФУ, и показано, что различные методы моделирования представительных объемов дают и различную степень точности в определении пьезомодулей d 31 и d 33. Также было установлено, что учет неоднородности поля поляризации существенно улучшает результаты КЭ моделирования, а наилучшим из рассмотренных методов можно считать метод ОДА Виттена-Сандера.

На втором этапе рассматривалась одномерная модель цилиндрического пьезоизлучателя, предназначенного для возбуждения в воде мощных ударных акустических импульсов короткой длительности. По методологии, представленной в [7,8] для плотной пьезокерамики, здесь были исследованы следующие типы пьезоизлучателей: 1) трехслойный пьезопреобразователь, состоящий из двух упругих согласующих слоев и пьезоэлемента, изготовленного из плотной керамики; 2) излучатель, состоящий из одного упругого согласующего слоя и пьезоэлемента, изготовленного из пористой керамики; 3) излучатель, состоящий из одного пьезоэлектрического слоя, изготовленного из высокопористой керамики.

В рабочем состоянии пьезоизлучатель генерировал импульсы в воду через наружный слой. Система возбуждалась разностью потенциалов, подаваемой на электродированные поверхности пьезоэлектрического диска. Во всех случаях излучатель состоял из дисков одинакового радиуса, причем для одномерной по пространственным координатам задачи существенны были только толщины дисков. В первом примере пьезоэлектрический диск был изготовлен из плотной керамики PZT-5H, поляризованной по толщине и характеризующейся плотностью, упругим модулем c c33, пьезомодулем e e33, диэлектрической проницаемостью, и добротностью Q1. Электроды были нанесены на торцы пьезоэлемента. Первый согласующий слой из стекла, приклеенный к пьезоэлементу, характеризовался толщиной h2, плотностью 2, скоростью распространения продольных волн v2 и добротностью Q2. Третий диск был выполнен из эбонита и характеризовался толщиной h3, плотностью 3, скоростью распространения продольных волн v3 и добротностью Q3. В случае двухслойного и однослойного пьезоэлектрических излучателей пьезоэлемент, выполненный из пористой керамики, характеризовался эффекТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

тивными модулями, полученными в первой части исследования для представительного объема, сгенерированного методом ОДА Виттена-Сандера, при 30 % и 60 % пористости.

Рис. 1. Давление в фиксированной точке, пористая керамика с одним согласующим слоем Рис. 2. Давление в фиксированной точке, пористая керамика без согласующих слоев В одномерной теории учитывались только толщинные движения пьезоизлучателя. Толщина пьезоэлектрического слоя выбирались равной половине длины продольной пьезожесткой волны на частоте ft и рассчитывалась по формуле h1 1 / 2 c /(2 f t ), а толщина j-го упругого слоя – по формуле h j j / 4 v j /( 4 f t ). Разность потенциалов на электродах задавалась в форме импульса, линейно возрастающего от нуля до 2 кВ за 1 мкс, постоянного значения 2 кВ в течение 1 мкс, и убывающего до нуля также за 1 мкс. В рассмотренных примерах пьезокерамический слой разбивался на 20 равных конечных элементов, упругие слои – на 10 элементов, и акустическая среда – на 60 элементов. Ниже приведены графики давлений в фиксированной точке акустической среды в зависимости от времени для случая пьезоизлучателя, имеющего один согласующий слой (рис. 1), а также излучателя без согласующих упругих слоев (Рис. 2), в сравнении с аналогичным результатом для плотной пьезокерамики с двумя согласующими слоями.

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

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

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

Список используемой литературы 1. Данцигер А.Я., Разумовская О.Н., Резниченко Л.А. и др. Многокомпонентные системы сегнетоэлектрических сложных оксидов: физика, кристаллохимия, технология. Аспекты дизайна пьезоэлектрических материалов. Ростов н/Д: Изд-во Рост. ун-та, 2002. Т. 2. 365 с.

2. Лопатин С.С., Лупейко Т.Г. Свойства пористой пьезоэлектрической керамики типа цирконата-титаната свинца // Изв. АН СССР. Сер. Неорг. Матер. 1991. Т. 27, № 9. С. 1948-1951.

3. Getman I., Lopatin S. Theoretical and experimental investigation of the porous PZT ceramics // Ferroelectrics. 1996. V. 186. P. 301-304.

4. Nasedkin A.V., Shevtsova M.S. Improved finite element approaches for modeling of porous piezocomposite materials with different connectivity. Nova Science Publishers, N.-Y., 2011.

Ch.7. P.231-254.

5. Domashenkina T.V., Nasedkin A.V., Remizov V.V., Shevtsova M.S. Finite element modeling of porous piezocomposite materials with different connectivity and applications for analysis of ultrasonic transducers // Proc. 7th GRACM Int. Congress on Comput. Mechanics, Athens, Greece, June 30 – July 2, 2011. CD. Paper 141. 10 p.

6. Witten T.A. Sander L.M. Diffusion-limited aggregation, a kinetic critical phenomenon // Phys. Rev. Lett. 1981. V. 47. P. 1400-1403.

7. Даниленко А.С., Наседкин А.В. Исследование импульсных характеристик многоcлойных пьезоизлучателей по МКЭ // Совр. пробл. мех. сплошной среды. Тр. V Межд.

конф., г. Ростов-на-Дону, 12-14 окт. 1999 г. Т. 2. Ростов-на-Дону: Изд-во СКНЦ ВШ, 2000.

С. 93-98.

8. Наседкин А.В., Даниленко А.С. Специальные формы одномерных конечных элементов для пьезоэлектрического анализа // Изв. вузов. Сев.-Кавк. регион. Естеств. науки. 2004.

Спецвыпуск. С. 78-82.

9. Рахимов В.Ф., Цихоцкий Е.С. Пьезоэлектрический преобразователь для силовой антенной решетки литотриптора // Пьезотехника-2005. Матер. Межд. научно-практич. конф.

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

Ростов-на-Дону, Азов, 23-26 августа 2005 г. Ростов-на-Дону: РГПУ, 2005. С. 154-156.

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

УДК 539.3:534. К идентификации неоднородных балочных конструкций при изгибно-крутильных The article describes the flexural-torsional vibrations of rods of variable stiffness. Develop a method for solving the inverse problem of reconstruction of heterogeneous characteristics: the dimensionless stiffness and dimensionless shear modulus.

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

Уравнения установившихся изгибно-крутильных колебаний консольно закрепленного стержня в безразмерных координатах имеют следующий вид:

где u (x) - смещение, (x ) - изгиб стержня, B (x) - безразмерная жесткость, C (x ) - безразмерный модуль сдвига, k, - безразмерные частоты колебаний, b1, b2, b3, b4 - коэффициенты пропорциональности, введенные при обезразмеривании и зависящие от геометрических характеристик стержня.

Будем считать, что балка на конце x 0 жестко закреплена, а на конце x l действует нагрузка в виде момента. Соответствующие граничные условия имеют вид:

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

Решение обратной задачи состоит в восстановлении функций B (x) и C (x ) по значениям функций u (x) и (x ), заданных в наборе точек. Функция B (x) восстанавливает из первого уравнения системы (1). Отметим, что данное дифференциальное уравнение относительно функции B (x) имеет второй порядок, и для ее определения необходимо решать задачу Коши.

Некорректность обратной задачи проявляется в том, что необходимо находить вторую производную от функции, заданной в наборе точек. Поэтому, на первом этапе аппроксимируем функции u (x) и (x ) сплайнами пятой степени, что позволяет находить вторую производную с достаточной степенью точности. На втором этапе построим операторное выражение, связывающее функцию B (x) с функциями u (x) и (x ). Дважды интегрируя первое уравнение системы (1) от 1 до x, получим операторное выражение для функции B (x) :

Проведена серия вычислительных экспериментов по восстановлению функции B (x) для различных видов неоднородностей: монотонно возрастающие функции, монотонно убывающие функции, немонотонные функции, кусочно-разрывные функции. На рисунках 1-2 приведены результаты восстановления некоторых функций.

Как видно на рисунках 1-2, погрешность восстановления функции B (x) не превышает двух процентов, что демонстрирует работоспособность предложенного метода.

Рис. 1. Решение обратной задачи для функции B x 1 x 2 (монотонно возрастающая функция)

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

Рис. 2. Решение обратной задачи для функции B x 1 e x (монотонно убывающая функция) Функцию С (x ) будем восстанавливать из второго уравнения системы (1) по заданной информации о функциях u (x) и (x ) и найденной функции B (x). Методом, описанным выше, получим операторное уравнение для функции С (x ) :

Из полученного операторного уравнения заметно, что предложенный способ не может быть применен для восстановления функции С (x ), так как в уравнении (4) в знаменателе стоит первая производная функции изгиба (x ), которая, согласном граничным условиям (2), в точке x 0 равна нулю. Следовательно, предложенный способ может быть применен только для восстановления безразмерной жесткости B (x).

Список используемой литературы 1. Филиппов А. П. Колебания деформируемых систем. Изд. 2-е, перераб. / А. П. Филиппов. – М.: Машиностроение, 1970. – 736 с 2. Ватульян А.О. Обратные задачи в механике деформируемого твердого тела. М.: Физматлит, 2007. 223с.

3. Березин И. С., Жидков Н. П. Методы вычислений./ И. С. Березин, Н. П. Жидков. - Государственное издательство физико-математической литературы. 1959. – 620 с.

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

УДК 539.3:534. Особенности конечно-элементного моделирования сейсмостойкости остова морской А. Н. Соловьев*, С. В. Красновская, В. В. Напрасников, М. А. Мирзаванд A frame of the offshore structure is represented by a piping system composed of pipes with various thickness, diameter and length. The proposed model which is as simple and efficient as a beam model can be used to obtain the seismic response of piping systems with an accuracy close to that obtained by a complicated finite element model. Seismic analysis of the piping system has been reviewed based on the applied horizontal and vertical accelerations to the ground nodes of the frame. The results of the experiments are presented.

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

В работе рассматривается методика моделирования напряженно-деформированного состояния остова морской платформы при различных сейсмических воздействиях на основе создания параметрических конечно-элементных моделей в среде ANSYS.

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

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

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

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

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

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

Для целей проектирования акселерограммы были нормированы по пиковому ускорению в соответствии с интенсивностью землетрясения. Оцифровки используемых акселерограмм нормировались по максимальной абсолютной величине 999. Оцифровка давалась с постоянным шагом времени 0,02 с. В ней также использовалась два масштабных коэффициента: первый множитель преобразует ее к исходной физической норме ускорений в долях ускорения свободного падения; второй — масштабирует акселерограмму к расчетному уровню.

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

Рассматривалась специфическая запись вертикального ускорения короткого сейсмического толчка с величиной 0,727g. Пиковое ускорение в акселерограмме приходится на момент

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

времени, равный 2,1 с. Максимальное перемещение при этом возникает на верхнем уровне. При этом превалирует составляющая суммарного перемещения в направлении аппликаты.

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

Рис. 2. Распределение эквивалентных напряжений по Мизесу Оцифровка акселерограммы длительностью 4 с представляет землетрясение магнитудой 6,8. Пиковое ускорение приходится на момент времени 2,1 с. Значения ускорений данной оцифровки, приведенные к расчетному уровню, использовались для моделирования землетрясений с вертикальным и с горизонтальным ускорениями. Максимальные напряжения и перемещения на момент приложения пикового ускорения в акселерограмме, равный 2,1 с, приведены в таблице 1.

Максимальное перемещение по компоМаксимальное напряженентам, м

UX UY UZ

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

При моделировании сейсмического воздействия с горизонтальным ускорением рассматривалась запись горизонтального ускорения при катастрофическом землетрясении магнитудой 8,0, произошедшем 16 сентября 1978 года в окрестностях иранского города Тебес. Акселерограмма представляла запись высокоинтенсивных сейсмических колебаний с амплитудой в 0.934g. Значения ускорений прикладывались к узлам, расположенным на нижнем уровне остова, по направлению оси ординат.

Оказалось, что для данного типа конструкции наиболее опасным является сейсмическое воздействие с горизонтальным ускорением. Результаты распределения напряжения по Мизесу в материале конструкции представлены на рисунке 2.

Как видно, максимальные напряжения составляют 310 МПа, что превышает предел текучести материала.

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

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

Выполнены расчеты для землетрясений магнитудой 6,8 и 8,0.

При воздействии землетрясения магнитудой 6,8 напряжения в материале конструкции не превышают предела текучести. При воздействии землетрясения магнитудой 8,0 напряжения в материале превышают предел текучести. Конструкция не выдержит моделируемого сейсмического воздействия. Землетрясение с такой магнитудой вызовет деформации в конструкции, которые приведут к ее разрушению.

Список используемой литературы 1. Якимуш И.С., Соловьев А.Н., Мирзаванд М.А. Особенности построения параметрических моделей для анализа прочности и сейсмостойкости морской нефтедобывающей платформы. 10-я Международная НТК «Наука – образованию, производству, экономике», Минск, 2012, с. 202.

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

УДК 539. Определение изменения формы поверхности непрерывно-неоднородного термоупругого This paper describes the axisymmetric quasi-static thermoelasticity problem for functionally-graded half-space which elastic modulus, Poisson’s ratio, heat conductivity and linear expansion coefficient continuously change in the subsurface layer. It is considered that circle area of radius is exposed to heat source and heat flow is constant in time. The surface outside this area is considered heat-isolated. The boundary problem of this problem is solved using Hankel integral transforms. The approximated analytical solution constructed by the asymptotical method is used to determine the heat flow and the half-space surface displacements. Numerical results are provided for different kinds of thermoelastic properties variation in the subsurface layer, and effect of such variation on surface displacements is thoroughly analyzed. The research is done with financial support from RFBR (projects 12-07-00639_a) and Federal Program on Scientific and Academic Personnel 2009-2013 (Agreement No. 14.740.21.1605).

Keywords: thermoelasticity, thermal deformation, nonhomogeneous materials.

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

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

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

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

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

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

Работа выполнена при финансовой поддержке грантов РФФИ (11-08-91168-ГФЕН_а, 12а), а также в рамках ФЦП «НиНПКИР 2009-2013» (соглашение № 14.B37.21.1632).

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

Список используемой литературы:

1. Коваленко А.Д. Введение в термоупругость. – Киев: «Наукова думка», 1965. – 204 с.

2. Айзикович С.М., Александров В.М., Васильев А.С., Кренев Л.И., Трубчик И.С Аналитические решения смешанных осесимметричных задач для функциональноградиентных сред. – М.: Физматлит, 2010. – 192 с.

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

УДК 539. Контактная задача о кручении упругой среды с неоднородным анизотропным покрытием The problem about torsion by a circular stamp of an elastic medium with an inhomogeneous coating is considered. The problem is reduced to dual integral equation using Hankel transform. Approximate analytical solution of high accuracy is constructed. We consider the case when the elastic properties of the coating vary in several times.

Distribution of contact stresses, strains and displacements are constructed not only on the surface, but also inside the material. Also we obtained simple approximate engineering formulas for calculation mechanical characteristics of the problem, which allows us to capture the qualitative differences of layered and functional-gradient materials. The accuracy of the approximate solutions is analyzed.

The research is done with financial support from State contract No. 11.519.11.3028 and Federal Program on Scientific and Academic Personnel 2009-2013 (Agreement No. 14.B37.21.1131).

Keywords: contact problems, coatings of complicated structure, torsion, functionally-graded coatings.

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

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

Приведен анализ точности полученных решений. На основе полученных результатов показано качественное отличие свойств слоистых и функционально-градиентных материалов.

Работа выполнена при финансовой поддержке ГК № 11.519.11.3028, а также в рамках ФЦП «НиНПКИР» (соглашение № 14.B37.21.1131).

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

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

Список используемой литературы:

1. Aizikovich S. Analytical solution of the spherical indentation problem for a half-space with gradients with the depth elastic properties/ S. Aizikovich, V. Alexandrov, J. Kalker, L. Krenev // International Journal of Solids and Structures 39 — 2002. — C. 2745–2772.

2. Aizikovich S. Evaluation of the elastic properties of a functionally-graded coating from the indentation measurements / S. Aizikovich, L. Krenev, I. Sevostianov, I. Trubchik //

ТРУДЫ X МЕЖДУНАРОДНОГО НАУЧНО-ТЕХНИЧЕСКОГО ФОРУМА

«ИННОВАЦИЯ, ЭКОЛОГИЯ И РЕСУРСОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ»

УДК 539. Математическое моделирование внедрения штампа в упругую среду с мягким покрытием This paper describes the mathematical model of contact between a punch and foundation, which corresponds to soft elastic layer lying on rigid strainless substrate. The problem is reduced to dual integral equation, which solution is constructed using dual-sided asymptotic method. The method described in this paper allows obtaining approximated solution of high accuracy. Numerical results are provided for cases when elastic properties of the coating and the substrate differ by 10, 100, 1000 times. Given results are compared to those obtained using other common methods.

The research is done with financial support from RFBR (projects 11-08-91168-GFEN_a, 12-07-00639_a) and Federal Program on Scientific and Academic Personnel 2009-2013 (Agreement No. 14.740.11.1605).

Keywords: contact problems, indentation, inhomogeneity, soft coatings.

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



Pages:     | 1 |   ...   | 2 | 3 || 5 | 6 |   ...   | 7 |


Похожие работы:

«A/66/661 Организация Объединенных Наций Генеральная Ассамблея Distr.: General 19 January 2012 Russian Original: English Шестьдесят шестая сессия Пункт 130 повестки дня Взаимодействие между Организацией Объединенных Наций, национальными парламентами и Межпарламентским союзом Вербальная нота Постоянного представительства Марокко при Организации Объединенных Наций от 12 января 2012 года на имя Генерального секретаря Постоянное представительство Королевства Марокко при Организации Объединенных...»

«в номере АкАдемик СкулАчев w w w.ek smo.ru АриАднА бориСовА издАтельСтво ЭкСмо — утверждАет: новый Автор в извеСтной победитель конкурСА Серии. душевное тепло продлить молодоСть — ревизор-2013 гАрАнтировАно возможно уже СейчАС 10 ок тябрь 2013 10 октябрь журна л распр ос траняется бесплатно АдреСА регионАльных СодержАние диСтрибуционных центров Новос ти изд ательс тва т орговый д ом ЭкСмо Фи ли А л ЭкСмо Ведущие проек ты изд ательс тва в роС т ове-н А-д он у 142701, Московская область, г....»

«корпорaтивный отдых остaвь корпорaтивный отдых остaновки корпорaтивный отдых остaновку корпорaтивный отдых остaться корпорaтивный отдых остром корпорaтивный отдых отвлеченные корпорaтивный отдых отвлеченный корпорaтивный отдых отдaть корпорaтивный отдых отделение корпорaтивный отдых отдых корпорaтивный отдых отзыв корпорaтивный отдых откровенный корпорaтивный отдых отличaющийся корпорaтивный отдых отменный корпорaтивный отдых отношение корпорaтивный отдых отношении корпорaтивный отдых отношения...»

«Список  доменов  Ru ­Center,  заблокированных  с  24  ноября  2010  г. ааааа.рф аааа.рф аанг.рф аарон ­авто.рф абажур.рф абакан ­автоматизация.рф абакана.рф абакан ­карта.рф абакан ­наутилус.рф абаков.рф абак.рф абактал ­инструкция.рф абактал.рф абап.рф абарис.рф аббревиатура.рф абб.рф абвгд.рф абвер.рф абдоминопластика.рф абд.рф абдулманов.рф абдулов ­александр.рф...»

«ФОРУМ: http://forum.babikov.com/ ТЕМА: Идеальное здоровье без лекарей и аптекарей  Don’t be stupid! (Не будь дураком!) Доктор ЗДОРОВЕНЬКИН ЧЕЛОВЕК – РАЗУМНЫЙ? ЛЕЧЕБНИК ДЛЯ МАЛОИМУЩИХ. 2005 г. ЦЕНА КНИГИ: Устанавливается Читателем  ОПЛАТА: Для [email protected] в системе Деньги@Mail.ru http://money.mail.ru/  Страница 1  ФОРУМ: http://forum.babikov.com/ ТЕМА: Идеальное здоровье без лекарей и аптекарей  ЛЕЧЕБНИК ДЛЯ МАЛОИМУЩИХ. ВВЕДЕНИЕ. Читатели, знакомые с моей книгой Терапия отчаяния,...»

«Главные новости дня 28 августа 2013 Мониторинг СМИ | 28 августа 2013 года Содержание СОДЕРЖАНИЕ ЭКСПОЦЕНТР 26.08.2013 Coffeetea.ru. Новости Московская кофейня на паяхъ выступила Генеральным спонсором United Coffee&Tea Industry Event и Командного Чемпионата Кофейных Энтузиастов United Coffee&Tea Industry Event (UCTIE) – главное индустриальной событие в России и других странах СНГ. Состоится 12-14 сентября в Москве, ЦВК Экспоцентр,. 7  27.08.2013 PublisherNews.ru. Новости предприятий и...»

«ПРАВА БЕЗДОМНЫХ КТО СЧИТАЕТСЯ БЕЗДОМНЫМ В Израиле, как и во все мире, не существует четкой и однозначной формулы, позволяющей определить, при каких условиях человек может считаться “бездомным”. Тем не менее, можно определить понятие “дом” как место жительства, отвечающее трем условиям: оно физически подходит для жизни человека; в нем у его обитателя есть возможность частной жизни; и, наконец, человек живет там на законых основаниях. Отсутствие одного или более из этих условий означает...»

«СТЕНОГРАММА круглого стола Комитета Государственной Думы по образованию на тему Вопросы здоровья в учреждениях профессионального образования: состояние и проблемы отрасли и законодательства Здание Государственной Думы. Зал 706. 7 июня 2012 года. 11 часов. Председательствует Гильмутдинов И.И. Дегтярёв А.Н. Добрый день, уважаемые коллеги, товарищи и друзья, соратники, все участники круглого стола! Комитет по образованию Государственной Думы Федерального Собрания Российской Федерации и подкомитет....»

«Хамрохон Зарифи Таджикистан в системе обеспечения региональной безопасности ОБСЕ Душанбе Ирфон 2011 ББК 66.4(2 тадж) + 66.5(2 тадж) + 66.5(0)6 3-34 Таджикистан в системе обеспечения региональной 3-34 Хамрохон Зарифи. безопасности ОБСЕ. - Душанбе, Ирфон, 2011 - 456 с. ISBN 978-99947-830-4-5 Книга содержит анализ деятельности ОБСЕ во всех трех измерениях, подход Организации в решении современных вызовов и угроз, отражает роль Республики Таджикистан в системе обеспечения региональной безопасности...»

«ПРОФЕССИОНАЛЬНЫЙ ФОРУМ AV FOCUS МОСКВА 5-6 СЕНТЯБРЯ Москва Генеральный спонсор: Спонсоры форума: УЧАСТНИКИ ФОРУМА ОРГАНИЗАТОР AV CLUB 127422, г. Москва, а/я 15 тел./факс: +7 495 780-0301 GSM: +7 962 935- E-mail: [email protected] www.avclub.ru AV FOCUS АВ Клуб - профессиональное сообщество на рынке AV индустрии - cредство коммуникации для профессионалов отрасли для обмена опытом и информацией об аудио-видео технологиях. Уникальное международное профессиональное сообщество, объединяющее...»

«ОРГАНИЗАТОРЫ ОРГАНИЗАТОР ФОРУМА Dealmakers forum ВЫСТАВКИ ИНВЕСТИЦИИ. СТРОИТЕЛЬСТВО. НЕДВИЖИМОСТЬ в РФ и странах СНГ 23 24 апреля 2008 года, Экспоцентр, Москва Аарон Голдштайн, Глава направления по развитию бизнеса на территории России и СНГ, NYSE Euronext Владимир Авдеев, партнер, генеральный директор, S.A.Ricci / King Sturge НОВЫЕ РЫНКИ, НОВЫЕ ВОЗМОЖНОСТИ КАЗАХСТАН ТАТАРСТАН Андрей Алешкин, директор департамента коммерческой недвижимости, Colliers International БЕЛАРУСЬ АЗЕРБАЙДЖАН в...»

«№ 17 198 А Н Т Р О П О Л О Г И Ч Е С К И Й ФОРУМ Мария Пироговская Ветлянская чума 1878–1879 гг.: санитарный дискурс, санитарные практики и (ре)формирование чувствительности Что чувства наши, или лучше сказать, что чувственность может быть изощреннее, то доказывали примеры чувств, из соразмерности своей болезнию выведенные [Радищев 1941: 139–140]. Воля к очищению требует противника своего масштаба. А для хорошо динамизированного материального воображения сильно загрязненная субстанция дает...»

«РеаСпоМед 2003 МАТЕРИАЛЫ 3 го Российского научного форума РеаСпоМед 2003 Москва, ЦДХ, 25 28 марта 2003 года Москва 2003 Материалы 3 го Российского научного форума РеаСпоМед 2003 М., Авиаиздат, 2003 216 с. Российская академия медицинских наук Мораг Экспо ISBN 5 94943 007 7 ©МОРАГ Экспо, 2003 ТЕЗИСЫ МИОТЕРАПИЯ ДЕТЕЙ С ПОСЛЕДСТВИЯМИ ПЕРИНАТАЛЬНОГО ПОРАЖЕНИЯ ЦЕНТРАЛЬНОЙ НЕРВНОЙ СИСТЕМЫ Аксенова А.М., Сереженко Н.П., Андреева В.В., Аксенова Н.И. Россия, г.Воронеж, государственная медицинская...»

«OFFSHORE MARINTEC RUSSIA _ СПРАВОЧНИК УЧАСТНИКА (Часть 2. Формы заявок) Offshore Marintec Russia 07 – 10 октября 2014 г. ЦВК ЭКСПОФОРУМ г. Санкт-Петербург Организатор мероприятия: Настоящий Справочник содержит всю информацию, необходимую для успешной подготовки к участию в выставке. Пожалуйста, найдите время прочесть его внимательно. Это позволит Вам избежать осложнений и дополнительных расходов на выставочной площадке. _ 1 07- 10 октября 2014 года, КВЦ Экспофорум г. Санкт-Петербург OFFSHORE...»

« МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РФ  Федеральное государственное бюджетное   образовательное учреждение   высшего профессионального образования   Пензенский государственный университет (ПГУ)    Пензенское региональное отделение  Общероссийской общественной организации  Ассоциация юристов России        Экстремизму – отпор!    Первый молодежный форум   Приволжского федерального округа    (г. Пенза, 2526 октября 2012 г.)    Материалы форума  Под редакцией  доктора юридических наук, доцента ...»

«Владимир Орлов Президент ПИР-Центра1 Терроризм как современная угроза глобальной безопасности: выводы для России и Индии и области для сотрудничества Доклад на V Дискуссионном форуме Россия и Индия: партнерство в глобальном формате 2 Москва 12 сентября 2011 г. Международное сообщество вступило в XXI век в сопровождении новых, нетрадиционных угроз глобальной безопасности. Не успев освободиться от страхов, которые в XX веке были вызваны гонкой вооружений и угрозой мировой войны с масштабным...»

«14-16 октября 2014 г. ТЕХНОПОЛИС МОСКВА www.forinnovations.ru Глобальная дискуссионная площадка в области инноваций Форум Открытые инновации 2014 Ключевая тема Форума 2014 Форум Открытые инновации посвящен новейшим технологиям и перспективам международной Созидательное разрушение: как сохранить конкурентоспособность в 21 веке кооперации в области инноваций. “ Скорость технологических изменений никогда еще не была настолько стремительной. Созидательно разрушая рынки с поСегодня ускоренный темп...»

«Непрерывное образование в сфере культуры №6/2009 Научно-практический форум ПРОЕКТ КОНЦЕПЦИЯ РАЗВИТИЯ ОБРАЗОВАНИЯ В СФЕРЕ КУЛЬТУРЫ И ИСКУССТВА В КРАСНОЯРСКОМ КРАЕ ДО 2020 ГОДА Введение Образование в сфере культуры и искусства является важнейшей составляющей культурно-образовательного пространства в современном российском обществе, национальным достоянием страны, обеспечивает сохранение высочайшего уровня профессионального искусства России. Оно социально востребовано как образование, органично...»

«АЗИАТСКО-ТИХООКЕАНСКИЕ ОРИЕНТИРЫ РОССИИ ПОСЛЕ САММИТА АТЭС ВО ВЛАДИВОСТОКЕ К ИТОГАМ ВТОРОГО АЗИАТСКО-ТИХООКЕАНСКОГО ФОРУМА №8 2013 г. Российский совет по международным делам Москва 2013 г. УДК 327(470:5) ББК 66.4(2Рос),9(59:94) А35 Российский совет по международным делам Редакционная коллегия Главный редактор: докт. ист. наук, член-корр. РАН И.С. Иванов Члены коллегии: докт. ист. наук, член-корр. РАН И.С. Иванов (председатель); докт. ист. наук, акад. РАН В.Г. Барановский; докт. ист. наук, акад....»

«Главные новости дня 10 июня 2013 Мониторинг СМИ | 10 июня 2013 года Содержание ЭКСПОЦЕНТР 10.06.2013 ТПП-Информ В Экспоцентре стартовал крупнейший смотр выставочной индустрии Через формат выставки для выставочников мы уже перешагнули. В 2012 году форум был удостоен знака Всемирной ассоциации выставочной индустрии UFI. И это говорит о том, что наш форум стал авторитетным, – сказал в своем приветственном слове генеральный директор Экспоцентра Сергей Беднов. 7  10.06.2013 Коммерсантъ Не салона...»










 
2014 www.av.disus.ru - «Бесплатная электронная библиотека - Авторефераты, Диссертации, Монографии, Программы»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.