Пружинки (COMBIN)

[Avksent 10]

Доброго времени года!
Тестирую работу с пружинами (в частности combin14).

Суть: создаю две точки, на одной масса, вторая закреплена. Между точками указанный combin14. Даю ускорение, всё супер, а вот как вытащить output по пружинке - не пойму.

 

Пробовал зайти через *GET:

*GET,My_force, ELEM, 1, SMISC, 1, , 

и через ETATABLE:

ETABLE, My_force, SMISC, 1,

 

Выдаёт, что нихт в database:  Item SMIS has not been stored in the database.  The *GET command is ignored.

 

Прошу помочь и объяснить где я не докурил HELP.

 

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

 

Листинг APDL:

FINISH  
/PREP7

ET,1,MASS21 
R,1,25,25,25, , , , 
 

ET,2,COMBIN14   
KEYOPT,2,1,0
KEYOPT,2,2,3
KEYOPT,2,3,0
R,2,25000, , , , , ,

 

n,1,0,0,0    
n,2,0,0,0.250  

real,1  
type,1 
e,2

 

real,2  
type,2 
e,1,2

d,1,all,0

 

acel,0,0,9.81

 

FINISH  
/SOL
solve

 

FINISH  
/POST1
plnsol,u,sum

[Борман 2 397]

@Avksent

*GET,My_force, ELEM,, SMISC, 1, , 

 

Изменено 19 августа 2022 пользователем Борман

[Avksent 10]

13 часов назад, Борман сказал:

@Avksent

*GET,My_force, ELEM,, SMISC, 1, , 

 

*GET,My_force, ELEM, 2, SMISC, 1, , - выдаёт ту же ошибку.

 

Возможна ли ошибка из-за того, что я запрашиваю данные не входя в узловую систему координат?

"The model contains an element type ( COMBIN14 ) that operates entirely in the nodal coordinate system.  Viewing nodal displacements or forces in other than the nodal coordinate system may be invalid.  See the elements manual description for more information."

Изменено 19 августа 2022 пользователем Avksent

[Борман 2 397]

55 минут назад, Avksent сказал:

*GET,My_force, ELEM, 2, SMISC, 1, , - выдаёт ту же ошибку.

Фигасе, на 17.1 работает, на 22 не работает.

 

 

Этой падле надо явно ткнуть сюда.

 

[Avksent 10]

19.08.2022 в 20:44, Борман сказал:

Фигасе, на 17.1 работает, на 22 не работает.

 

 

Этой падле надо явно ткнуть сюда.

 

 

Неожиданно.

Благодарю!

 

Если что, итоговый код (мало ли):

 

FINISH  
/PREP7

ET,1,MASS21 
R,1,25,25,25, , , , 
 

ET,2,COMBIN14   
  
KEYOPT,2,1,0
KEYOPT,2,2,3
KEYOPT,2,3,0
R,2,25000, , , , , ,

DOF,UX,UY,UZ,ROTX,ROTY,ROTZ 

n,1,0,0,0    
n,2,0,0,0.250  

real,1  
type,1 
e,2

real,2  
type,2 
e,1,2

d,1,all,0

/eshape,1,1
/eshape,0,0  

d,2,UZ,0.05

FINISH  
/SOL
outres,misc,all,
solve

FINISH  
/POST1
plnsol,u,sum

SET,FIRST

*GET,My_force, ELEM, 2, SMISC, 1, , 

[Avksent 10]

Коллеги, дабы не засорять темы.

Начал смотреть элемент COMBIN37.

В APDL работает нормально (немного другая задача).

В WB реализовать не могу.

Суть: есть два тела (кубика), температуру на одном контролируем узлом K COMBIN37, узлы I и J заделываем на поверхностях кубиков. Хочу добиться того, чтобы COMBIN37 включался при определённой температуре.

ContrBody - NS на котором контролируем температуру.

Температуру на теле задаём Thermal Condition.

Файл проекта прилагаю.

 

 

Проблема: на узел K не передаётся температура от тела (равна 0 при любой температуре тела).

Если температуру задавать в узле K через D,...,TEMP,моя_температура - всё работает.

При этом перемещения в узле J совпадает с перемещением тела, к которому его привязываю.

Забираю командной вставкой

/POST26

/SHOW,png
/TITLE, UK-Temperature
ESOL,2,numElem,,NMISC,8,UK-Temperature
/AXLAB,Y,UK
/YRANGE, -100, 100, 
/XRANGE, 1, 5,
PLVAR,2
/SHOW,close

Текстовка командной вставки:

/PREP7

 

DOF,UX,UY,UZ,ROTX,ROTY,ROTZ,TEMP

 

oon1 = 10                                            !назначаем минимальную температуру работы элемента COMBIN37
oof1 = 100                                          !назначаем максимальную температуру работы элемента COMBIN37

 

!Определение номеров узлов для I, J, K
*GET,lastNumberNode,node,,num,max
NODEI=lastNumberNode+1
NODEJ=lastNumberNode+2
NODEK=lastNumberNode+3

 

!Присвоение номеров уравнений связи узлам I, J, K
LastCP=0
*GET,lastCP,CP,0,max
cpI=lastCP+1
cpJ=lastCP+2
cpK=lastCP+3

 

!Присвоение номера КЭ COMBIN37
*GET,lastType,etyp,0,num,max
numType=lastType+1

 

!Присвоение RealConstant COMBIN37
numReal=numType

 

!Соединение узла I с целевым узлом (можно с поверхностью) через CP - общие степени свободы
NSEL,s,node,,Selection_I
N,NODEI,0,0,10
CP,cpI,UZ,all
ALLSEL,all

 

!Соединение узла J с целевым узлом (можно с поверхностью) через CP - общие степени свободы
NSEL,s,node,,Selection_J
N,NODEJ,0,0,100
CP,cpJ,UZ,all
ALLSEL,all

 

!Соединение узла K с целевым объёмом второго тела (похоже тут проблема)
NSEL,s,node,,ContrBody
N,NODEK,0,0,0
CP,cpK,TEMP,all
ALLSEL,all

 

!Назначение аргументов и Real Constant COMBIN37
ET,numType,COMBIN37
KEYOPT,numType,2,8              !Degree of freedom for control nodes (K and L): TEMP
KEYOPT,numType,3,3              !Degree of freedom for active nodes (I and J): UZ (along nodal Z)
KEYOPT,numType,4,1              !ON-OFF range behavior: Unique ranges
KEYOPT,numType,5,0              !OFF-either-ON (or OFF-ON-OFF if unique)

R,numReal,10, , ,oon1 ,oof1 , 0  
RMORE, , 0, , , , ,  
RMORE, ,

 

!Создание элемента COMBIN37 между I и J с управлением по узлу K (без L, но всё тоже самое)
real,numReal
type,numType
e,NODEI,NODEJ,NODEK

 

*GET,numElem,elem,0,num,max     !Сохраняем номер элемента COMBIN37 для построцессора

 

/SOLU

 

 

02.wbpz

Изменено 23 августа 2022 пользователем Avksent

[Avksent 10]

Немного дополню:

В элементах куба температура соответствует заданным значениям через Thermal Condition

 

А вот в узле куба - температуры нет...

[soklakov 1 485]

3 часа назад, Avksent сказал:

заданным значениям через Thermal Condition

Гляньте инпутник. Это BFE по идее. Не температура. Вот и не передаёт.

[Avksent 10]

43 минуты назад, soklakov сказал:

Гляньте инпутник. Это BFE по идее. Не температура. Вот и не передаёт.

Termal Condition записан как:

/com,*********** Define Body Force Temperature ***********
esel,none
esel,a,type,,2
cm,_CM37,elem
esel,all
*DIM,_loadvari37,TABLE,6,1,1,TIME,
! Time values
_loadvari37(1,0,1) = 0.
_loadvari37(2,0,1) = 1.
_loadvari37(3,0,1) = 2.
_loadvari37(4,0,1) = 3.
_loadvari37(5,0,1) = 4.
_loadvari37(6,0,1) = 5.
! Load values
_loadvari37(1,1,1) = -100.
_loadvari37(2,1,1) = -100.
_loadvari37(3,1,1) = 50.
_loadvari37(4,1,1) = 150.
_loadvari37(5,1,1) = 50.
_loadvari37(6,1,1) = 150.

..................................................

cmsel,s,_CM37
nsle
bf,all,temp,%_loadvari37%
nsel,all
esel,all

 

Но у меня вопрос: а в чём разница между заданием BF (Defines a nodal body force load) и заданием "Температуры"?

[Борман 2 397]

37 минут назад, Avksent сказал:

и заданием "Температуры"?

Если она задается как d,x,temp,x то в этом случае вы явно задаете степень свободы, а в случае bf - объемную нагрузку (источником которой, среди прочего, может быть и температура)

[Avksent 10]

18 часов назад, Борман сказал:

Если она задается как d,x,temp,x то в этом случае вы явно задаете степень свободы, а в случае bf - объемную нагрузку (источником которой, среди прочего, может быть и температура)

Так, вроде начал понимать.

Тогда вопрос немного с другой стороны - допустим есть связанный расчёт: термический -> МДТТ.

Существует ли возможность передать данные из термического расчёта в механический так, чтобы можно было использовать COMBIN37 как контрольный?

 

Пока попробовал сделать этот расчёт напрямую и через external data - в обоих случаях данные на тела передаются через функцию BF.

Технически ведь external data берёт температуру с узлов - наверняка есть возможность передать температуру в узлы в механическом расчёте как степень свободы, а не объёмной нагрузкой.

[Борман 2 397]

46 минут назад, Avksent сказал:

Существует ли возможность передать данные из термического расчёта в механический так, чтобы можно было использовать COMBIN37 как контрольный?

Ну наверное если контрольный узел прицепить к термической сетке. ХЗ.

[soklakov 1 485]

3 часа назад, Avksent сказал:

Технически ведь external data берёт температуру с узлов - наверняка есть возможность передать температуру в узлы в механическом расчёте как степень свободы, а не объёмной нагрузкой.

В механическом расчете нет такой степени свободы - температура. Возможно, вы захотите термомеханический full coupled. И передавать ничего не придется.

[Avksent 10]

6 часов назад, soklakov сказал:

В механическом расчете нет такой степени свободы - температура. Возможно, вы захотите термомеханический full coupled. И передавать ничего не придется.

Да, уже посмотрел в эту сторону, но в итоговой модели предполагается использование балочных элементов, которые не завезли в coupled анализ))) Может и придётся, но сделаю ещё пару попыток.

Пока думаю делать так: Делать термический расчёт, переносить их external dat'ой в structural, выгружать результаты по узлам в файл, открывать файл и искать температуру по ближайшему узлу к контрольной точке, приравнивать температуру в K в этой точке. Костыль конечно жуткий, но как есть.

 

Была ещё идея воровать переданные данные с помощью вставки

NSEL,s,node,,выбранное_тело
*GET,узлы_в_теле,node,,num,max
*GET,температура_в_элементе,NODE,узлы_в_теле,BFE,TEMP
D,NODEK,TEMP,температура_в_элементе

 

Но не взлетел.

Изменено 25 августа 2022 пользователем Avksent

[soklakov 1 485]

15 часов назад, Avksent сказал:

в итоговой модели предполагается использование балочных элементов, которые не завезли в coupled анализ)

В ВБ действительно не подвезли. Но по факту то есть?

Хм... А может и нет) pipe разве что)