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

III МОДЕЛЮВАННЯ ПРОЦЕС1В В МЕТАЛУРГП ТА

МАШИНОБУДУВАНН1

УДК 531; 539.3

Канд. техн. наук Ю. А. Лымаренко Запорожская государственная инженерная академия, г. Запорожье

ОПРЕДЕЛЕНИЕ БАЗОВЫХ ВАРИАНТОВ ДЕФОРМИРОВАНИЯ ПРОСТРАНСТВЕННОГО ДИСКРЕТНОГО ЭЛЕМЕНТА

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

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

Введение

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

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

Симметризация и антисимметризация вектора перемещений дискретного элемента

Рассмотрим дискретный элемент в форме куба (рис. 1). Перемещения вершин элемента в направлении координатных осей, сонаправленных ребрам параллелепипеда, обозначим через ц, У1, wi, 1= 1, 8. Тогда перемещение всего элемента будет задаваться вектором в двадцатичетырехмерном пространстве:

Элементы этого вектора представляют собой коэффициенты в разложении вектора и по базисным векторам:

X! =(1, 0, ..., 0),

X2 =(0, 1, ..., 0),

х24 =

© Ю. А. Лымаренко, 2016

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

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

Рассмотрим три плоскости симметрии: 1-2 - плоскость, проходящая перпендикулярно стороне 1-2 через ее середину; 1-4; 1-5 - аналогично.

Вначале выполняем симметризацию 8 и антисимметризацию А относительно плоскости 1-2 (рис. 2). Под симметризацией здесь и далее подразумевается операция

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

uS = (u1 + u2 )/2, ua = (u1 - u2 )/2 и т. д. При этом, чтобы не загромождать изображения, на рисунке 2 (и на всех последующих) верхние индексы «s» и «a», обозначающие соответственно результат выполнения операций симметризации и антисимметризации, опущены.

Затем для каждого из полученных случаев выполняем симметризацию S и антисимметризацию A относительно плоскости 1-4, получая уже четыре случая: SS, SA, AS, AA (рис. 3).

Далее для каждого из полученных случаев выполняем симметризацию S и антисимметризацию A относительно плоскости 1-5, получая восемь случаев: SSS, SSA, SAS, SAA, ASS, ASA, AAS, AAA (рис. 4).

Рис. 3. Второй этап симметризации перемещений: относительно плоскости 1-4

Рис. 4. Третий этап симметризации перемещений: относительно плоскости 1-5

В каждом из построенных восьми случаев выделяем варианты, содержащие только ненулевые значения и, только ненулевые значения у и только ненулевые зна-

чения те, получая в итоге 24 окончательных результата, представленных на рисунке 5.

Рис. 5. Геометрическая иллюстрация новых базисных

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

У1 =

У 4 = У 7 = У10 = У13 = У16 =

У19 =

У 22 =

0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0),

0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, -1, 0, 0, -1, 0, 0, -1, 0, 0, -1, 0, 0),

0, 0, 1, 0, 0, -1, 0, 0, -1, 0, 0, 1, 0, 0, 1, 0, 0, -1, 0, 0, -1, 0, 0),

0, 0, 1, 0, 0, -1, 0, 0, -1, 0, 0, -1, 0, 0, -1, 0, 0, 1, 0, 0, 1, 0, 0),

0, 0, -1, 0, 0, -1, 0, 0, 1, 0, 0, 1, 0, 0, -1, 0, 0, -1, 0, 0, 1, 0, 0),

0, 0, -1, 0, 0, -1, 0, 0, 1, 0, 0, -1, 0, 0, 1, 0, 0, 1, 0, 0, -1, 0, 0),

0, 0,-1, 0, 0, 1, 0, 0, -1, 0, 0, 1, 0, 0, -1, 0, 0, 1, 0, 0, -1, 0, 0),

0, 0, -1, 0, 0, 1, 0, 0, -1, 0, 0, -1, 0, 0, 1, 0, 0, -1, 0, 0, 1, 0, 0).

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

В отличие от базисных векторов (2) новые базисные вектора имеют наглядный физический смысл. Каждый

из этих векторов отвечает за один из базовых вариантов деформирования пространственного дискретного элемента:

У1, У2, У3 - поступательные перемещения; У6, У8, У13 - растяжения-сжатия;

УН , У 12 , У16 , У 18 , У 19 , У20 - изгибы;

УЮ , У 17, У21 - кручения;

У 22, У 2з, У 24 - самоуравновешенные деформации;

У4, У15, У5, У9, У7, У14 - попарно комбинация сдвига и поворота.

Из последней шестерки векторов можно с помощью линейной комбинации в явном виде выделить:

У4 + У15 , У5 + У9 , У7 + У14 - сдвиги;

У15 - У4 , У9 - У5 , У14 - У 7 - пов°р°тЫ.

Представим вектор (1) в виде разложения по новым базисным векторам:

24

U = Z sy .

i=1

Обратим внимание на то, что вектора уг являются взаимно ортогональными:

Уi • У] = 0 (, } = 1, ..., 24; г * ]).

Это позволяет находить коэффициенты разложения (з) следующим образом:

— • yi = siyi • yi.

Следовательно,

U • ys' = —(i = 1,..., 24).

y- •y- v ;

Выводы

Использование операции симметризации вектора перемещений позволило выделить базовые варианты

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

Список литературы

1. Шамровский А. Д. Решение плоских статических задач механики деформируемого твердого тела при помощи дискретных моделей, получаемых на основе экспериментальных данных / А. Д. Шамровский, Ю. А. Лыма-ренко, Д. Н. Колесник // Проблеми обчислювально! ме-ханжи i мщност конструкцш. Зб. наук. праць. - Дншро-петровськ : Лра. - 2011. - Вип. 17. - С. 274-288.

2. Шамровський О. Д. Дискретна модель плоского еле-менту сюнченних розмiрiв для ортотропного середови-ща / О. Д. Шамровський, Т. О. Мшяйло // Методи розв'я-зування прикладних задач мехашки деформiвного твердого тша. - Дншропетровськ : Лiра, 2012. - Вип. 13. -С. 428-436.

3. Шамровский А. Д. Напряжения в деформированном дискретном элементе / А. Д. Шамровский, Д. Н. Колесник // Восточно-Европейский журнал передовых технологий. Серия «Прикладная мехашка». - 2013. -Вып. 3/7 (63). - С. 12-15.

4. Шамровский А. Д. Усовершенствованный метод последовательных перемещений для расчета пространственных стержневых конструкций / А. Д. Шамровский, Т. А. Миняйло, Д. Н. Колесник // Новi матерiали i технологи в металургй та машинобудуванш. - 2012. - № 1. -С. 83-86.

Одержано 20.10.2016

Лимаренко Ю.О. Визначення базових itapiai гмв деформування просторового дискретного елемента

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

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

Lymarenko Yu. Defining a basic types of spatial discrete element deformations

Construction a discrete models of continuous elastic medium involves a problem of model elastic characteristics determination. Technique of finding a basis in the space of displacement vectors of spatial discrete element is proposed. Obtained basis will enable in future to evaluate the sought stiffnesses by comparing a results of solving an inverse problems in discrete and continues formulations.

Key words: continuous elastic medium, the spatial discrete element, displacement, basis, symmetrization and antisymmetrization vectors.