Контактная задача

[intruder777 3]

Как определить площади контакта при давлении одного тела на другое, в ansys версии 18?https://fex.net/294333308941?fileId=433521091

Изменено 23 апреля 2018 пользователем intruder777

[piden 2 874]

8 hours ago, intruder777 said:

в ansys версии 18?

 

!contact area macro 
!This macro extracts the contact area from a group of elements
!The resulting value is listed in the command object in the solution branch 
!
!User must have a named selection called cnt_press_edges, scoped to one or more edges or faces in this case
!These faces must have contact elements otherwise the results will return zero
!If user is using program controlled for contact behavior, ANSYS may internally flip the contact and target (auto symmetric). If this happens, the results will be returned as zero
!In this case the user has to make sure that edges picked for the named selection have contact elements.

!Here are some options;
!If behavior is set to program controlled and results come back as zero, then go to the named selection and pick the opposing face
!If behavior is set to asymmetric, then just make sure the named selection points to the same faces as contact faces in the contact definition
!If behavior is set to symmetric, then run and record the contact area, save this number, then in the named selection, select the opposing face instead,
!rerun, extract the contact area, add to the first contact area from the prior step, and divide by two. This will give an average of the contact area
!
!Sean Harvey - ANSYS Inc.
!Rev 1.0 -  Original release  11/4/2013
!Rev 1.1 -  Modified by piden 22/4/2018

!user assumes all responsbility from use of the macro

resume
set,last                                !set the last results, change as needed!!!  Especially if multi results sets

allsel

*get,num_el,elem,0,count                !get number of contact elements selected
*get,start_id,elem,0,num,min            !get id of first contact element
*dim,cnt_area,array,num_el              !dimension the contact area array

cmsel,s,cnt_press_edges                 !select nodes of interest. Need face name selection on area of interest
esln,s,1                                !select elements attached to these nodes            
esel,r,ename,,174                       !reselect the contact elements. Change element type if needed

*vget,elem_sel,elem,start_id,esel       !create mask vector of selected 1=selected -1=not selected

allsel

*vabs,1
*vmask,elem_sel                         !turn on mask
*vmask,elem_sel                         !turn on mask
*vscfun,elem_tot_area,sum,elem_area     !sum the area of the contact elements

    
*vmask,elem_sel                         !turn on mask
etable,carea,nmisc,58                   !get the contact area. Refer to HELP to find proper nmisc value for CAREA for chosen element
*vmask,elem_sel                         !turn on mask
*vget,cnt_area,elem,start_id,etab,carea !store in array
*vmask,elem_sel                         !turn on mask                   
*vscfun,tot_area,sum,cnt_area           !sum the contact area


my_contact_area = tot_area              !output variable tot_area

 

Ansys 19.0

 

 

[intruder777 3]

@piden Спасибо, а можно понять что такое 174 и 58? почему здесь 

"esel,r,ename,,174 " две запятых? 
Реально крутая штука спасибо! Буду пробовать разобраться! 
А что с сеткой? У меня правильно была?

Как насчет настроек контакта? Они правильные?

Изменено 23 апреля 2018 пользователем intruder777

[piden 2 874]

2 minutes ago, intruder777 said:

Как насчет настроек контакта? Они правильные?

Ну, контакт с ними считается. А площадь - я не пробовал, свои поставил. Изменил behavior = unsymmetric.

3 minutes ago, intruder777 said:

Спасибо, а можно понять что такое 174 и 58?

Да, можно. Открой хелп и прочти про элемент CONTA174. И найди, что для этого элемента выводит nmisc 58.

Желтым было выделено то, что я изменил в исходном скрипте.

4 minutes ago, intruder777 said:

А что с сеткой? У меня правильно была?

Приходилось видеть и красивее... Но как по мне - покатит.

9 minutes ago, intruder777 said:

две запятых? 

в APDL если параметр пропускаешь, то используется значение по-умолчанию.

Например

k,1,0,0,0 - создается точка с номером 1 и коорд. (0,0,0)

а

k,,0,0,0 создаст еще одну точку в этом же месте, но номер ей присвоит автоматом. В данном случае 2.

 

[intruder777 3]

@piden

cmsel,s,cnt_press_edges                !select nodes of interest. Need face name selection on area of interest
esln,s,1     

Нужно ли изменить здесь cnt_press_edges на 

cnt_area?

[piden 2 874]

Нет! Просто скопируй скрипт себе и замени тот, что был, на этот. И больше в скрипте ничего не трогай)

[Борман 2 375]

1 час назад, piden сказал:

k,,0,0,0 создаст еще одну точку в этом же месте, но номер ей присвоит автоматом. В данном случае 2.

А что делает команда k,,,, ?

[piden 2 874]

20 minutes ago, Борман said:

А что делает команда k,,,, ?

создается точка с коорд. (0,0,0). Номер присваивается автоматически.

 

1 hour ago, piden said:

в APDL если параметр пропускаешь, то используется значение по-умолчанию.

А команда

k,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,

делает за тебя расчет, пишет отчет и получает зп. Главное, пару лишних зпт не добавить. А то пойдет и спустит всё на женщин и вино...

Изменено 23 апреля 2018 пользователем piden

[Борман 2 375]

12 минуты назад, piden сказал:

А команда

k,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,

делает за тебя расчет

по умолчанию..
 

12 минуты назад, piden сказал:

пишет отчет

по умолчанию...
 

12 минуты назад, piden сказал:

получает зп

по умолчанию.

Номер присваивается автоматически.

[Mrt23 173]

13 минуты назад, piden сказал:

создается точка с коорд. (0,0,0). Номер присваивается автоматически.

 

А команда

k,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,

делает за тебя расчет, пишет отчет и получает зп. Только не всегда срабатывает(

А если я з/п сам хочу получить, какие запятые не нужно ставить?

[piden 2 874]

Just now, Mrt23 said:

А если я з/п сам хочу получить, какие запятые не нужно ставить?

Вот ради этих сакральных знаний люди и изучают apdl годами!

[Mrt23 173]

Провал...

Не тем софтом занимался:(

[Борман 2 375]

Только что, piden сказал:

Вот ради этих сакральных знаний люди и изучают apdl годами!

 

2 минуты назад, Mrt23 сказал:

А если я з/п сам хочу получить, какие запятые не нужно ставить?

Нигде. Надо в параметрах запуска вместо "-$" использовать ключ "+$"

Например "+$1000"

[piden 2 874]

3 minutes ago, Борман said:

Надо в параметрах запуска вместо "-$" использовать ключ "+$"

Нда. А говорят рыбаки суеверный народ 

 

4 minutes ago, Mrt23 said:

Не тем софтом занимался:(

Вишь, а в бабакусе-то вашем все ключи через "-"!  Тота и оно!

6 minutes ago, Борман said:

Надо в параметрах запуска вместо "-$" использовать ключ "+$"

Например "+$1000"

Борман походу к БД бухгалтерии доступ получил....

[Борман 2 375]

6 минут назад, piden сказал:

 БД бухгалтерии

bd.db

[Mrt23 173]

Ну да...

Не зря @Борман так любит apdl... :)

[piden 2 874]

2 minutes ago, Борман said:

bd.db

fini

/clear

/exit

- следы замети 

[hr4d 87]

@pidenА заделись ресурсами, которые читаешь/мониторишь по поводу CAE, а? 

Если не сложно конечно. 

[piden 2 874]

5 minutes ago, hr4d said:

заделись ресурсами, которые читаешь/мониторишь по поводу CAE

Да я никакие не мониторю) Когда что-то нужно - в гугл вбиваю и ищу.

Мониторю только официальные порталы ансиса / симулии / альтаира на предмет выхода новых туторов. Но я их, увы, не выполняю, а складываю в сундучок)

[hr4d 87]

2 минуты назад, piden сказал:

Да я никакие не мониторю) Когда что-то нужно - в гугл вбиваю и ищу.

Мониторю только официальные порталы ансиса / симулии / альтаира на предмет выхода новых туторов. Но я их, увы, не выполняю, а складываю в сундучок)

Блин,  такая же фигня. Папка customer portal уже неприлично разрослась((