Как заменить матрицу жёсткости в ANSYS

[Igoshechka 1]

Необходимо провести расчёт для одномерной задачи со скачком. Есть стержень с областями  V1 и V2, между ними граница. На границе узлы имеют двойную нумерацию. Есть два способа, первый - прописать в ANSYS условие контакта на границе, но я не знаю как это сделать. (Буду рад, если кто-нибудь подскажет.) Знаю, что условие контакта имеет вид:

u_i (на левой границе) = u_i (на правой границе) - U_i, где U_i - скачок в перемещениях, который задан по условию задачи.

 

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

 

Я получил из ANSYS Mechanical APDL глобальную матрицу жёсткости с помощью следующих двух команд:
*DMAT, MatKD, D, IMPORT, FULL, file.full, STIFF
*PRINT, MatKD, Kdense.matrix

 

После чего отредактировал её и хочу загрузить обратно. Какие команды для этого использовать?

Также есть ещё один вопрос: как получить из ANSYS вектор нагрузки F? (Который используется вместе с матрицей жёсткости A*U=F.)

Заранее спасибо.

[Fedor 1 601]

"А контакты - нелинейная задача)" смотря какие. Например уместная здесь задача о вдавливании штампа вполне себе линейная. Задача Герца нелинейная, только из-за переменной области контакта... :)

[soklakov 1 475]

7 часов назад, Fedor сказал:

Например уместная здесь задача о вдавливании штампа вполне себе линейная. Задача Герца нелинейная, только из-за переменной области контакта... :)

Так и у штампа область контакта переменная. И поверхности скользят друг по другу.

[Fedor 1 601]

Обычно под этим понимается плоский штамп, да и трением пренебрегают или считают склеенным если численно моделируют  ... 

[Fedor 1 601]

https://new.chitai-gorod.ru/product/osnovy-teorii-uprugogo-diskretnogo-kontakta-2535795   Надо почитать ... 

[soklakov 1 475]

4 часа назад, Fedor сказал:

Обычно под этим понимается плоский штамп

Когда кастрюлю из плоского диска штампуют - это плоский штамп?

[Fedor 1 601]

Тут обычный для теории упругости смысл вдавливания в полупространство. Надо чувствовать контекст. Смекать .  http://www.pinega3.narod.ru/mg/mg.htm  :)