Подскажите по Ansys'у...

[Chingachguk 4]

Всем, Добрый вечер!

Решаю задачу устойчивости подкрепленной пластины в линейной постановке. Дана пластина (моделирую shell181) подкрепленная двумя ребрами (моделирую beam188). Составил код, сейчас как есть - пластина и ребра работают раздельно, а нужно совместно. Каким образом их можно связать в единую конструкцию или beam188 и shell181 не связываются?

Спасибо.

/PREP7 !Входим в препроцессор
!        Геометрические параметры
t=0.02  !Толщина пластины
a=0.4   !Расстояние между ребрами
b=3.6   !Длина пластины
h1=0.06 !Высота ребра
b1=0.02 !Ширина ребра
q=1     !Единичная сила
!        Задание ключевых точек
K,1,0,0,0,  
K,2,b,0,0,  
K,3,b,0,a,  
K,4,b,0,a*2,
K,5,b,0,a*3,
K,6,0,0,a*3,
K,7,0,0,a*2,
K,8,0,0,a,
!        Построение линии
LSTR, 1,2
LSTR, 3,8
LSTR, 4,7
LSTR, 5,6
LSTR, 2,3
LSTR, 3,4
LSTR, 4,5
LSTR, 6,7
LSTR, 7,8
LSTR, 8,1
!
!        Построение поверхности
AL,1,5,2,10
AL,2,6,3,9
AL,3,7,4,8
!
!        Задание свойств материала
!        Для плиты
UIMP,1,EX,,,2e11, !Модуль Юнга
UIMP,1,NUXY,,,0.3, !Коэффициент Пуассона
!        Для ребра
UIMP,2,EX,,,2e11, !Модуль Юнга
UIMP,2,NUXY,,,0.3, !Коэффициент Пуассона
!
!        Задание типа элемента
ET,1,SHELL181
ET,2,beam188
!
!        Блок задания констант элемента
!Для пластины
sect,1,shell,,shell181
secdata, t,1,0.0,3
secoffset,MID
seccontrol,,,, , , ,
!
!Для ребра
SECTYPE,2,BEAM,RECT,Beam_188,0
SECOFFSET,CENT 
SECDATA,h1,b1,0,0,0,0,0,0,0,0,0,0
!
!Разбиваем на элементы
LSEL,S,LINE,,1,4,
LESIZE,All,,,45,1 !Присвоить 45 делений на элементы
LSEL,S,LINE,,5,10
LESIZE,All,,,5,1 ! Присвоить 5 делений на элементы
!
!        Присвоим атрибуты поверхности
ASEL,S,AREA,,All !Выделение поверхности
AATT,1,,1,0,1    !Присвоить материал 1
LSEL,S,LINE,,2,3 !Выделение линии
LATT,2,,2,,,,2   !Присвоить материал 2
!
AMESH,All !Построить сетку на всех поверхностях
!
LMESH,All !Построить сетку на всех линиях
NUMMRG,NODE, , , ,LOW !Сшить совпадающие узлы
!
FINISH !Выход из препроцессора
!
/SOLU  !Вход в процессор решения
!
!      Задание условий закрепления
LSEL,S,LINE,,All
LSEL,U,LINE,,2,3
DL,All,,UY, !Закрепление всех выбранных линий
LSEL,S,LINE,,8,10
DL,All,,UX, !Закрепление всех выбранных линий
DL,All,,UZ, !Закрепление всех выбранных линий
!
!      Задание нагрузок
LSEL,S,LINE,,5,7
SFL,All,PRES,q
ALLSEL,ALL 
FINISH !Выход из препроцессора решения
!
/SOLU !Вход в препроцессор решения
PSTRES, ON !Установить вычисление напряжённого состояния
SOLVE !Запуск на решение
FINISH !Выход из препроцессора решения
/SOLU  
ANTYPE, 1 !Анализ на устойчивость в линейной постановке
BUCOPT, SUBSP,10 !Определить 10 форм потери устойчивости
MXPAND, 10,,,0
SOLVE !Запуск на решение
FINISH  
/POST 1 !Вход в постпроцессор
SET, FIRST !Прочитать первый ряд расчётных значений
PLDISP, 0 !Вывести деформированную форму графически
SET,LIST 

[AlexKaz 1 822]

Не заданы контакты, см. элементы Conta###, Target###

[soklakov 1 475]

Каким образом их можно связать в единую конструкцию или beam188 и shell181 не связываются?

Можно без контактов , сливайте узлы. Склеивайте геометрию в один part.

[Chingachguk 4]

Даа, не так все просто оказывается...)))

Вроде имеются узлы (nodes) определяющие конечные элементы оболочек и стержневых конечных элементов, узлы принадлежат как Оболочкам так и Балкам. Но почуму-то работают раздельно??? Пробовал такую задачу решить в другом комплексе (типа Лира) нет проблем...

to AlexKaz вреде это не контактная задача? У Вас есть пример подобный, не пойму как их применить то...

to soklakov Спасибо за примерчик, но есть две проблемы 1) не смог просчитать выкидывает ошибку типа:

cannot open message catalog "1049\libiomp5ui.dll"

2) А каким образом можно СЛИТЬ узлы в Mechanical APDL? есть там такая функция в APDL или через интерфейс...? 

 

PS: Ansys 14.0

PPS: нашел один пример правда не понятно в какой версии Ansys (думаю очень ранней), но там без контактов и слития узлов, просто объединением узлов через NUMMRG,NODE, , , ,LOW , используются элементы shell181 и beam 4. Попробовал реализовать через beam4 но не решает, видимо много изменений было...

Изменено 13 апреля 2015 пользователем Chingachguk

[kristeen 355]

не смог просчитать выкидывает ошибку типа: cannot open message catalog "1049\libiomp5ui.dll"

http://www.edr.no/blogg/ansys_bloggen/ansys_v14_bug
 

[Chingachguk 4]

Спасибо kristeen уже нашел))) первый вариант с игнором проще всего.

[Борман 2 375]

http://forum.dwg.ru/showthread.php?t=119564

[AlexKaz 1 822]

Не могу разобраться с командами сохранения массива значений в файл из-под LS-Dyna pre/post. Значения на графике закинул в массив, не понятно что делать дальше. Метод http://www.cae-club.ru/content/vyvod-dannykh-v-protsesse-resheniya у меня не работает, "cfopen is not recognized POST1 command".

[AlexKaz 1 822]

А, всё норм, через диалог Write to File всё работает

[AlexKaz 1 822]

Rigid Body возможно промоделировать, повысив E на пару порядков?

[soklakov 1 475]

Rigid Body возможно промоделировать, повысив E на пару порядков?

В статике и вообще в неявных схемах - да.

В явной динамике - нет. 

[Hulget 0]

Подскажите, пожалуйста, как пользоваться функцией Multiply? У меня есть решение задачи в КЭ постановке. Как, собственно, вытащить значения температуры и размера(в моем случае) для каждого КЭ?

[AlexKaz 1 822]

В постпроцессоре выделить элементы, внизу списка пост-комманд есть пункт 'закинуть в таблицу'. Потом из меню датабейз есть возможность сохранить таблицу в файл (см. выше два моих сообщения). Ну или все можно сделать быстрее и проще, но то есть в бездонной справке...

[Борман 2 375]

Подскажите, пожалуйста, как пользоваться функцией Multiply? У меня есть решение задачи в КЭ постановке. Как, собственно, вытащить значения температуры и размера(в моем случае) для каждого КЭ?

 

Тут вы найдете ответы на оба вопроса

help/ans_apdl/Hlp_P_APDLget.html 

[AlexKaz 1 822]

Явный солвер LS-Dyna решает задачу с временным шагом, пропорциональным размеру элемента. Допустим (как в статике в хэлпе ANSYS) ради ускорения расчёта детальку увеличиваем на N порядков (размер элемента растёт), свойства материалов изменяем соответственно. Интересует, возможно ли так поступить с явной задачей?

[soklakov 1 475]

Интересует, возможно ли так поступить с явной задачей?

Нет.

Потому как

свойства материалов изменяем соответственно.

 

Это если точно такой же результат нужен.

Но можно сильно завышать плотность материала для сведения задачи к квазистатической. Но если уж Вы взялись решать явную динамику, то, видимо, статика Вас не интересует.

Явный солвер LS-Dyna решает задачу с временным шагом, пропорциональным размеру элемента.

Не только.

[AlexKaz 1 822]

А квазистатической задача является до приблизительно каких скоростей, не постесняюсь спросить?  И т.е. один из самых важных параметров в динамике - плотность? У меня разлетающихся за горизонт ошмёток металла не предполагается, кинетическая энергия почти никакая - индентор елозит по поверхности пластины со скоростью 0.13 м/с. Хорошо бы пронаблюдать волну деформирования, но пока что не вижу связи деформации с плотностью - статика с задачками Герца считается без плотности, что в теории, что в МКЭ.

[Борман 2 375]

явную динамику

Если говорить о LS-Dyna, то почему после всего джентельментского минимума мероприятий, описанного в хелпе, эффект Hourglassing все равно остается ?

 

В частности, вот эта "фича" вообще не работает.

KEYOPT(1)

Element formulation:

 

0, 1 -- Constant stress solid element (default)

 

2 -- Fully integrated selectively-reduced solid 

 

[soklakov 1 475]

Хорошо бы пронаблюдать волну деформирования, но пока что не вижу связи деформации с плотностью

 

А вот для этого плотность нужна, и нужна правильная. Про скорость звука в твердых телах почитайте.

 

почему после всего джентельментского минимума мероприятий, описанного в хелпе, эффект Hourglassing все равно остается ?

Отправьте на портал  Вопрос интересный, но с ним явно возиться надо. И не мне. 

[AlexKaz 1 822]

"Пока что не вижу связи деформации с плотностью"

Увидел у Мейза в "Теории и задачах МСС", правда что-то мне кажется при такой низкой скорость уравнения не сильно будут отличаться от статики...

 

Что-то сплайн по узлам совсем не похож на более-менее гладкий, узлы на месте, а функция скачет. Задавал в path разное количество точек и делений отрезков, не помогает.