Анализ способа решения плоской контактной задачи методом конечных элементов. Разработка алгоритмов определения контактирующих пар "узел–ответная поверхность". Особенности формирования глобальной матрицы жесткости для программного комплекса Mathcad.
Аннотация к работе
Алгоритм решения плоской контактной задачи Периодически в расчетной практике возникает ситуация, когда необходимо более точно определить зону приложения нагрузки - решить контактную задачу [1,2]. Известные программные комплексы имеющие возможность анализа контактного взаимодействия не всегда применимы ввиду имеющихся каких-либо особенностей. Перемещения по нормали к поверхности контакта совпадают из условия непроникновения одного из контактирующих тел (далее домена) в другое [5,9], по касательной к поверхности зоны контакта - из условия равновесия действия сил трения между контактирующими поверхностями, тогда для контактирующих поверхностей: (1) где, ucont , vcont - горизонтальное и вертикальное перемещение штампа соответственно, utarg , vtarg - соответственно горизонтальное и вертикальное перемещение тела основной конструкции. Т.е. при выполнении расчета с применением МКЭ необходимо связать перемещения соответствующих узлов обоих доменов, а так как расчет проводится итерационным способом, то, соответственно, необходимо выполнять обновление связей между узлами на каждой итерации в соответствии с перемещениями, полученными на предыдущей итерации.