close

Вход

Забыли?

вход по аккаунту

?

Алгоритм расчета контактных границ с учетом эрозии конечных элементов при высокоскоростном взаимодействии тел.

код для вставкиСкачать
2016. Т. 21, вып. 3. Физика
УДК 539.3
DOI: 10.20310/1810-0198-2016-21-3-1275-1277
АЛГОРИТМ РАСЧЕТА КОНТАКТНЫХ ГРАНИЦ
С УЧЕТОМ ЭРОЗИИ КОНЕЧНЫХ ЭЛЕМЕНТОВ
ПРИ ВЫСОКОСКОРОСТНОМ ВЗАИМОДЕЙСТВИИ ТЕЛ
 П.А. Радченко
Томский государственный архитектурно-строительный университет,
e-mail: radchenko@live.ru
В работе предложен алгоритм расчета контактных границ взаимодействующих тел в полной трехмерной пост ановке. Данный алгоритм поддерживает высокую степень распараллеливания вычислений, что актуально для к ачественных конечно-элементных сеток. Эрозия элементов позволяет поддерживать высокую скорость численного счета, а также сохранять устойчивость решения и шаг по времени при больших скоростях нагружения.
Адекватность предложенного алгоритма подтверждена сопоставлением с рядом экспериментальных исследований.
Ключевые слова: математическое моделирование; разрушение; динамическое нагружение; метод конечных элементов; упругость; пластичность.
В настоящее время существует множество пакетов
программ для численного моделирования поведения
материалов и конструкций при динамическом нагружении (ANSYS, NASTRAN, SIMULIA и т. д.). Несмотря на отличия в интерфейсе пользователя и отдельных
пакетах, данные комплексы основываются на основных
уравнениях механики сплошных сред и алгоритмах
расчета поверхностей взаимодействующих тел. Как
правило, все эти программы используют алгоритмы
«элемент-узел» и «узел-узел» для определения возможного проникания одного тела в другое. В данной
работе авторами предлагается отличающийся скоростью работы и надежностью алгоритм типа «элементэлемент» в полной трехмерной постановке.
В начале работы авторского программного комплекса EFES [1] загружаются сетки конечных элементов из стороннего сеточного генератора NETGEN [2].
Каждое тело обладает нумерованными поверхностями,
причем нет ограничений на количество конечных элементов, а также сложность геометрии тел. Нумерация
поверхностей глобальная для всего ансамбля взаимодействующих тел, и пользователь программного комплекса EFES явно указывает пары поверхностей, которые будут взаимодействовать в ходе решения поставленной задачи. Наиболее емким, с точки зрения объема
вычислений модулем, в методике расчета динамического нагружения тел является расчет контактных границ, а именно, перебор всех возможных взаимодействий отдельных конечных элементов.
На каждом шаге интегрирования явной конечноразностной схемы подсчитывается количество треугольников, определяющих каждую из участвующих в
контактном взаимодействии поверхностей. Поскольку,
благодаря возможной эрозии расчетной сетки при разрушении, число элементов может уменьшаться, то и
число данных поверхностных треугольников может
варьироваться. Зная текущее количество треугольни-
ков, определяется размерность и выделяется память
под соответствующие динамические массивы для используемых в алгоритме контакта переменных. Далее
определяются коэффициенты A, B, C (направляющие
косинусы вектора нормали к поверхности), D уравнения Ax  By  Cz  D  0 , описывающего в пространстве плоскость каждой поверхностной грани конечного элемента, образуемую узлами n1, n2, n3 (рис. 1).
Определяются и заносятся в созданные массивы
координаты узлов поверхностных элементов, образующих ребра тетраэдров от поверхностных узлов к
внутренним: m1m4, m2m4, m3m4. Проверка проникания
поверхностного элемента одного тела в поверхностный
элемент другого сводится к задаче проверки пересечения трех ребер тетраэдра с треугольной поверхностной
гранью другого:
Рис. 1. Общая схема контактного взаимодействия «элементэлемент»
1275
ISSN 1810-0198. Вестник Тамбовского университета. Серия Естественные и технические науки
i 
( A  x4  B  y 4  C  z 4 )
,i  1,2,3
( A  ( x4  xi )  B  ( y4  yi )  C  ( z 4  zi ))
Здесь x, y, z с индексом 4 – координаты внутреннего
узла рассматриваемого тетраэдра, с индексом i – координаты поверхностных узлов m1, m2, m3, а A, B, C, D –
коэффициенты уравнения плоскости, образуемой узлами n1, n2, n3 соответственно. Если коэффициент
0 < i < 1, то соответствующее ребро проникло сквозь
поверхностную грань конечного элемента. Поскольку
точка пересечения ребра тетраэдра с плоскостью может
находится и вне грани n1n2n3, происходит проверка
принадлежности точки пересечения соответствующему
треугольнику. Для этого определяется сумма углов,
образуемых данной точкой с узлами n1, n2, n3. При ее
равенстве 2 точка лежит внутри грани данного тетраэдра. Так как через одну плоскость могут проникнуть
несколько ребер тетраэдра, но само ребро может пересекаться лишь с одним треугольником в пространстве,
то организовывая цикл перебора всех возможных
взаимодействий по узлам m1, m2, m3 каждого тела,
можно получить крайне высокое распараллеливание
данного алгоритма. По окончании данного цикла все
данные о полученных прониканиях элементов заносятся в соответствующий массив.
В соответствии с моделью скольжения поверхностей без трения [3] проникший узел m2 сносится на
граничную поверхность по нормали в положение m2'
(рис. 1), в дальнейшем обозначим эту точку как С. Изменение количества движения данного узла в нормальном направлении определяется соотношением:
Js 
рости узла m2 Vsn . Доли импульса Ri , R j , R k находятся из условий сохранения количества движения и
момента количества движения. Изменение скорости
узла m2 определяется из условия равенства нормальных
скоростей данного узла и поверхности элемента в точке С. Для определения скорости поверхности элемента
в положении m2' используется линейная аппроксимация. Тогда выражение для скорости поверхности элемента uc в направлении X имеет вид:
uc  qi ui  q j u j  qk uk ,
где qi  A jck Aijk . В данной формуле A jck является
площадью треугольника, образованного узлами с индексами j и k, а также точкой С, на которую сносится
узел m2. Данный треугольник является противолежащим по отношению к узлу с индексом i. Aijk – площадь всей грани, образованной узлами с индексами i, j,
k. Аналогично определяются константы qi и qk . Другие компоненты скорости ( vc и wc ) определяются из
выражения, аналогичного выражению для uc .
Приравнивая нормальные скорости узла m2 и поверхности элемента в положении C, находим:
Vsn 
Vcn  Vsn
,
1  qi Ri M s M i  q j R j M s M j  qk Rk M s M k
где Vcn и Vsn – нормальная скорость поверхности
элемента и скорость узла m2, нормальная к поверхности элемента:
M s Vsn ,
где M s – масса узла m2; Vsn – нормальное изменение скорости данного узла. Импульс J s передается
трем узлам элемента n1, n2, n3, обозначим их i, j, k соответственно. Доли импульса, передаваемые узлам треугольного элемента, обозначим
Ri  J i J s ,R j  J j J s ,Rk  J k J s ,
Vsn  Aus  Bv s  Cws ,
Vcn  Auc  Bv c  Cwc .
В данных уравнениях A, B и C – направляющие косинусы вектора нормали к поверхности элемента. Таким образом, сначала определяются Vsn , а затем
n
однозначно вычисляются Vi n , V j и Vkn .
где
J i  M i Vin ,J j  M j V jn ,J k  M k Vkn .
n
n
n
Здесь mi , m j , mk ,Vi , V j , Vk – масса и измене-
ние нормальных скоростей узлов i, j, k соответственно.
n
n
n
Тогда для определения Vi , V j , Vk имеют место
соотношения
Vi n  Ri M s Vsn mi ,
V jn  R j M s Vsn m j ,
Vkn  Rk M s Vsn mk .
n
n
n
Для определения Vi , V j , Vk необходимо оп-
ределить доли момента Ri , R j , Rk
1276
и изменение ско-
После того, как все проникшие узлы сносятся на
поверхности элементов, с которыми есть пересечение и
произошла корректировка скоростей, как этих узлов,
так и узлов, образующих грань, сквозь которую они
проникли, реализуется алгоритм эрозии для разрушенных ячеек расчетной сетки. Те из них, которые разрушились на текущем шаге интегрирования и находятся
на поверхности тела, поочередно удаляются из счета.
При этом элементы, которые имеют смежные с ними
грани ячеек, становятся поверхностными и включаются
в алгоритм корректировки контактных границ на следующих этапах расчета. При удалении из расчетной
сетки элемента масса оставшихся смежных узлов корректируется для сохранения массы всей расчетной области в целом.
Данный алгоритм позволяет адекватно описывать
поведение тел различной геометрии при численном
моделировании деформации и разрушения конструкций при импульсном нагружении [4]. Применение па-
2016. Т. 21, вып. 3. Физика
раллельных вычислений, многоуровневая оптимизация
данного алгоритма в рамках программного комплекса
EFES [1] позволяет существенно сократить время счета
и расширить масштаб решаемых с его помощью задач.
СПИСОК ЛИТЕРАТУРЫ
1.
2.
Радченко П.А., Батуев С.П., Радченко А.В. Трехмерное моделирование деформации и разрушения гетерогенных материалов и конструкций при динамических нагрузках (EFES 1.0). Свидетельство
о регистрации программы для ЭВМ № 2014614671.
Schöberl J. NETGEN – An advancing front 2D/3D-mesh generator
based on abstract rules // Computing and Visualization in Science.
1997. V. 1 (1). Р. 41-52.
3.
4.
Johnson, G. R. Tree-dimensional computer code for dynamic response
of solids to intense impulsive loads // Int. J. Numer. Methods Eng.
1979. V. 14. № 12. Р. 1865-1871.
Радченко П.А., Батуев С.П., Радченко А.В., Плевков В.С. Численное моделирование разрушения оболочки из бетона и фибробетона при импульсном воздействии // Омский научный вестник. 2015.
№ 3 (143). С. 345-349.
БЛАГОДАРНОСТИ: Работа выполнена при финансовой поддержке РФФИ (грант № 16-31-00125 мол_а,
грант № 16-38-00256 мол_а).
Поступила в редакцию 10 апреля 2016 г.
UDC 539.3
DOI: 10.20310/1810-0198-2016-21-3-1275-1277
ALGORITHM OF CALCULATION OF CONTACT BOUNDARIES
TAKING INTO ACCOUNT THE EROSION OF FINITE ELEMENTS
AT HIGH-VELOCITY INTERACTION OF BODIES
 P.A. Radchenko
Tomsk State University of Architecture and Building, Tomsk, Russian Federation,
e-mail: radchenko@live.ru
In this paper we propose an algorithm of calculating contact boundaries of the interacting bodies. This algorithm maintains a high degree of parallelization which is important for high-quality finite element meshes.
The erosion is to support high speed computations and also to maintain the stability of the solution and the
time step at high loading rates. The adequacy of the proposed algorithm is confirmed by comparison with
several experimental studies.
Key words: mathematical modeling; fracture; dynamic loading; finite element method; elasticity; plasticity.
REFERENCES
1.
2.
3.
4.
Radchenko P.A., Batuev S.P., Radchenko A.V. Trekhmernoe modelirovanie deformatsii i razrusheniya geterogennykh materialov i
konstruktsiy pri dinamicheskikh nagruzkakh (EFES 1.0). Svidetel'stvo o registratsii programmy dlya EVM no. 2014614671.
Schöberl J. NETGEN – An advancing front 2D/3D-mesh generator based on abstract rules. Computing and Visualization in Science,
1997, vol. 1 (1), pp. 41-52.
Johnson, G. R. Tree-dimensional computer code for dynamic response of solids to intense impulsive loads. Int. J. Numer. Methods Eng.,
1979, vol. 14, no. 12, pp. 1865-1871.
Radchenko P.A., Batuev S.P., Radchenko A.V., Plevkov V.S. Chislennoe modelirovanie razrusheniya obolochki iz betona i fibrobetona
pri impul'snom vozdeystvii. Omskiy nauchnyy vestnik – Omsk Scientific Bulletin, 2015, no. 3 (143), pp. 345-349.
GRATITUDE: The work is fulfilled under financial support of Russian Fund of Fundamental Research (grant no. 16-3100125 мол_а, grant no. 16-38-00256 мол_а).
Received 10 April 2016
Радченко Павел Андреевич, Томский государственный архитектурно-строительный университет, г. Томск, Российская Федерация, кандидат физико-математических наук, доцент кафедры прикладной математики, e-mail:
radchenko@live.ru
Radchenko Pavel Andreevich, Tomsk State University of Architecture and Building, Tomsk, Russian Federation, Candidate of Physics and Mathematics, Associate Professor of Applied Mathematics Department, e-mail: radchenko@live.ru
1277
Документ
Категория
Без категории
Просмотров
6
Размер файла
171 Кб
Теги
граница, конечный, высокоскоростная, алгоритм, элементов, взаимодействия, эрозии, расчет, тел, контактные, учетом
1/--страниц
Пожаловаться на содержимое документа