Решение задач динамики и модальный анализ

[Борман 2 376]

Борман не томите.

Слету не получилось. Вот частоты без химии и с полной матр. масс.

1 0.14759

2 0.15714

3 0.30034

4 0.33721

5 0.48540

6 0.51476

7 0.59621

8 0.72508

9 0.81578

10 0.89806

11 0.99229

12 1.4587

Кстати, мои элементы, это вырожденные гексы. А ваши ?

[Борман 2 376]

С полной матрицей совпадает.

Ну тогда вот вам с диаональной

1 0.11435

2 0.12485

3 0.18656

4 0.23443

5 0.24179

6 0.30511

7 0.32084

8 0.38218

9 0.38590

10 0.43053

11 0.57218

12 0.69056

Командный фпйл приложил, и комменты написал.

Описание всех команд здесь <noindex>http://www.kxcad.net/ansys/ANSYS/ansyshelp/Hlp_C_CmdTOC.html</noindex>

FILE.txt

[Борман 2 376]

Кстати, на APDL можно написать простенький экспортёр в формат ИСПА. Только вот не уверен, поддерживает ли APDL русский язык. Найти бы еще тот док, который вы мне летом скидывали...

[Гость ISPA]

Борман, посмотрите на правильном ли я пути.

3.rar

[Борман 2 376]

Борман, посмотрите на правильном ли я пути.

По синтаксису должно работать. Дома проверю правильно ли расставлены узлы в элементе.

ET,х,45

KEYOPT,х,1,1 - этим выключается химия для типа х

KEYOPT,х,1,0 - а вот так включается (так по умолчанию)

TYPE,х

[Борман 2 376]

Борман, посмотрите на правильном ли я пути.

Все работает.

[Борман 2 376]

Вроде наладил экспорт в ваш формат. Сетку приложил.

Маленькие глюки с русским языком. Надеюсь вы не обидитесь.

Возможно у вас пойдут глюки... Эту сетку я целмком заполнил 8-узл. элементами, и перебирал узлы так, как они перебираются в АНСИС. Возможно там есть 8-узл. тетраедры и треугольные призмы. В общем смотрите, пишите комменты. Кстати скрипт для экспорта оказался очень даже простенький..

EXPORT.rar

[Борман 2 376]

У них не 6 узлов, а 8.

Не понял. У меня и так 8. Или вам нужен элемент у которого 6 узлов ? Может возникнуть мелкая проблемка... из картинки все понятно. У меня сетки нет перед глазами, посмотрите, совпадет с чем узел элемента-гексы (?), с узлом соседнего гекса "1", или с узлом "2" ? Проверьте так каждый узел пирамиды...

Короче говоря одинаковые ли номера соответсвующих узлов у гекс вокруг пирамиды ? Думаю в АНСИСе они разные...

С утра не проснулся еще.. мыслю сформулировать не могу..

[Борман 2 376]

Думаю в АНСИСе они разные...

Хотя сетка из ГМ. В общем камень преткновения...

[Fedor 1 597]

"Только вот не уверен, поддерживает ли APDL русский язык" - вывод то можно организовать, но проще на cpp. Когда-то давно я написал приятелю такую примочку-утилитку с конвертацией туда-сюда элементов :)

[Борман 2 376]

Тьфу, понял. Вопросов нет. Я думал у моей призмы все 8 узлов разные...

Короче, экспорт очень простой...

[Борман 2 376]

Понятно.

Только вот у элемента нет признака формы, по-этому придется делать проверку узлов на совпадение. Это замедлит и без того долгий экспорт..

Господа, кто-нибудь в курсе, какая команда быстрее отрабатывает ? NENEXT, или аналогичный *GET ?

[Борман 2 376]

Вот новый экпорт. Добавил, точнее попытался, добавить нужные блоки. Только вот "КОНЕЦ" никак не пишется. Не пойму почему. Русский до кучи глючит... Правда это мое первое обращение к задаче вывода в файл через ADPL, подумаю еще, может тонкостей каких не знаю.

EXPORT.rar

[Борман 2 376]

Запишите тетра модель для проверки.

EXPORT3.rar

[Борман 2 376]

ISPA

а у вас как дела с экспортом в АНСИС ?

[Fedor 1 597]

" "КОНЕЦ"" - а Вы попропуйте 3.141592653589793-ц :) Неужели трудно по человечески {... } блоки оформлять?

[Борман 2 376]

Вернее как правильно задать нагрузку собственным весом. Откуда АНСИС знает величину g.

Не знаю как правильно. Я это делаю так...

Сначала задается плотностть в свойствах материала...

MP,DENS,x,7800

Затем задается ускорение всех системы (если Y вверх)

ACEL,0,9.81,0 - это, собственно, и моделирует ускорение (св. падения).

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

Понятно, что 7800 и 9.81 - в соответствующих единицах должны быть.

To simulate gravity (by using inertial effects), accelerate the structure in the direction opposite to gravity. For example, apply a positive ACEL_Y to simulate gravity acting in the negative Y direction. Units are length/time2.

[Борман 2 376]

Тогда как в размерности кГ, см при определении собственных частот получить Гц. Что-то подсказывает мне что без g не обойтись.

Не понял слово "тогда". В лобовом частотном анализе нагрузки не учитываются. Для их учета надо решать доп. задачу.

Как я выше писал... так решаются статические задачи с учетом веса.

[Гость ISPA]

Я по другому спрошу. Я определяю собственные частоты в размерности кГ, см. Задаю 2000000, 0.00785. Частоты в АНСИС я получу в Гц?

[Vova 419]

Я по другому спрошу. Я определяю собственные частоты в размерности кГ, см. Задаю 2000000, 0.00785. Частоты в АНСИС я получу в Гц?

Вроде вначале этой темы, этот вопрос уже поднимался.

Модуль Юнга 2.0E+6

Плотность 7.85e-6

Тогда получаются частоты. Так в Настране, в Ансисе наверняка тоже так.