Как изменить тип КЭ?

[Rigid 0]

Подскажите, можно ли изменить тип КЭ-элемента, выбранного через ESEL?

Проблема вот в чем (коротко). Некий объем (цилиндр+усеченный конус+отверстия в нем) представляет собой жидкость, налитую в бак. Объем разбивается на гексы и используется элемент fluid80. При запуске программа ругается на degenerate bricks, видимо потому, что fluid не может быть пирамидальным или клиновидным. Пирамид там нет, а вот клиновидные элементы встречаются. Сделать сетку из чистых гексов не представляется возможным (мне), ввиду особенностей объема (отверстия и пр.) Видится такой путь: те элементы, на которые ansys ругается, указать как solid185. Для всего объема заданы свойства жидкости, поэтому может такая замена не сильно исказит общую картину, учитывая, что проблемных элементов около сотни а всего их свыше 20000.

[Борман 2 391]

Вот так

<noindex>http://kxcad.net/ansys/ANSYS/ansyshelp/Hlp_C_EMODIF.html</noindex>

[Rigid 0]

Вот так

<noindex>http://kxcad.net/ansys/ANSYS/ansyshelp/Hlp_C_EMODIF.html</noindex>

Спасибо, я подозревал этот путь.

[Rigid 0]

а вот такой вопрос: можно ли программно определить форму 3D-элемента... то есть какой он получился, pyramid, prism или hex? Тогда можно будет прогнать все элементы и те, которые призмы, заменить.

[Влад. 6]

У хексы первого порядка 8 узлов - i, j, k ... p или 1, 2, 3 ... 8

Командой *GET, Par, ELEM, N, NODE можно в цикле по всем элементам считывать номера узлов.

Далее сравниваешь номера. Если у узлов o и p, а также узлов k и l номера одинаковые - элемент призматический.

Аналогично ищешь тетры.

[Rigid 0]

У хексы первого порядка 8 узлов - i, j, k ... p или 1, 2, 3 ... 8

Командой *GET, Par, ELEM, N, NODE можно в цикле по всем элементам считывать номера узлов.

Далее сравниваешь номера. Если у узлов o и p, а также узлов k и l номера одинаковые - элемент призматический.

Аналогично ищешь тетры.

понял, а по общему кол-ву нельзя?... если узлов меньше 8 то призма.... или совпадающие узлы все равно физически существуют и учитываются?

[Влад. 6]

Хелп нужно порыть, может и есть что проще.

В принципе за пол часа можно тот макрос набрать.

Кстати можно! Выделяешь в цикле элемент, принадлежащие ему узлы и считываешь их количество той же командой *GET.

*do

esel

nsle

*get,n,node,count

*if

*enddo

[Rigid 0]

Хелп нужно порыть, может и есть что проще.

В принципе за пол часа можно тот макрос набрать.

Кстати можно! Выделяешь в цикле элемент, принадлежащие ему узлы и считываешь их количество той же командой *GET.

*do

esel

nsle

*get,n,node,count

*if

*enddo

ага... я уже попробовал... если 6 узлов выводит, то призма.. спасибо !

[Борман 2 391]

Вот общий случай... разбирался когда делал экспорт (по линейным КЭ) в формат ИСПА.

*IF,NELEM(ne,5),EQ,NELEM(ne,6),THEN !ТЕТРАЭДР

...

*ELSEIF,NELEM(ne,3),EQ,NELEM(ne,4),THEN !ПРИЗМА

...

*ELSE ! ГЕКСАЭДР

...

*ENDIF

Вроде работало.