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

[kristeen 355]

http://cccp3d.ru/topic/51252-opredelenie-tcentra-mass-v-nds/#entry456837

[Борман 2 375]

Проект перенесен из WB в классику. Как теперь узнать МЦХ сетки в существующей СК ? Вроде бы можно выловить консольный листинг при запуске процесса решения, то есть не начав решать никак не узнать?

Численно взять интегральчик не боитесь ?

 

CENTRX(N), CENTRY(N), and CENTRZ(N) always retrieve the element centroid in global Cartesian coordinates

*GET, VOL_N, ELEM, N,VOLU

Ну и просуммировать в цикле.

[Борман 2 375]

Ну еще строго говоря надо выдернуть номер материала, а из номера метариала выдернуть плотность. :) В общем, это только глаза боятся...

[Fedor 1 597]

При вычислении матрицы жесткости вычисляются интегралы по объему элемента от функций. Вычислить просто объем это считать функции единицей, совсем просто. Для любого конечного элемента, не обязательно пирамидки. Вообще успех мкэ во многом связан с возможностью легко вычислять объемы по элементам. Поэтому в частности обычно используют и полиномы в качестве базисных функций. Хотя несложно построить и элементы из других наборов функций не обязательно полиномов. Можно и часть не полиномов задействовать, только зачем это вот в чем вопрос, вот в чем загвоздка, сказал бы принц датский :)

[Борман 2 375]

При вычислении матрицы жесткости вычисляются интегралы по объему элемента от функций. Вычислить просто объем это считать функции единицей, совсем просто. Для любого конечного элемента, не обязательно пирамидки. Вообще успех мкэ во многом связан с возможностью легко вычислять объемы по элементам. Поэтому в частности обычно используют и полиномы в качестве базисных функций. Хотя несложно построить и элементы из других наборов функций не обязательно полиномов. Можно и часть не полиномов задействовать, только зачем это вот в чем вопрос, вот в чем загвоздка, сказал бы принц датский :)

Мы по Гауссу-Остроградскому привыкши.:) 

[Fedor 1 597]

Да обычно численно по Гауссу через декартово произведение. К поверхностным не приводят. Хотя есть и варианты сразу объемы интегрировать. Немного поменьше надо точек интегрирования. Есть в Справочнике по специальным функциям под редакцией Абрамовица и Стигана. Том есть и гауссовы точки для точности 2 n-1 то есть наивысшей. Но это сейчас не особо актуально их легко и посчитать в Mathematica например :)

 

https://ru.wikipedia.org/wiki/%D0%A4%D0%BE%D1%80%D0%BC%D1%83%D0%BB%D0%B0_%D0%93%D0%B0%D1%83%D1%81%D1%81%D0%B0_%E2%80%94_%D0%9E%D1%81%D1%82%D1%80%D0%BE%D0%B3%D1%80%D0%B0%D0%B4%D1%81%D0%BA%D0%BE%D0%B3%D0%BE  это что-то не помню чтобы использовалось...

[Борман 2 375]

У нас что-то все в НИИ Механики при ННГУ очень любили мусолить тему про метод граничных элемнтов. Не знаю, чем у них все закончилось...

http://www.mech.unn.ru

[Fedor 1 597]

Есть несколько книжек. Банерджи П., Баттерфилд Р. Методы граничных элементов в прикладных науках. -М.: Мир, 1984 в ней кстати чудесно многие базисные функции описаны для мкэ.

Тема как-то не пошла. Матрицы получаются заполненными. Примерно то же самое можно получить и исключая внутренние узлы как в суперэлементах  делают. Или можно конструировать элементы исключая внутренние узлы. Например кубики из пирамидок и что-то подобное. Единственно интересное это бесконечные границы например для моделирования полупространства. Но обычно этим не заморачиваются. Просто берут побольше земли и все дела. Но в этой книжке есть пример сколько примерно надо брать земли в конечных элементах чтобы получать как в граничных ... Это многим у нас не нравится :)

Думаю что тихо увяло как везде :)

[AlexKaz 1 821]

Спасибо всем! Единственное, уточню в соответствии с http://ansys.net/tips_sheldon/STI08_PSOLVE.zip

 

! **** CENTER OF MASS, MASS, AND MASS MOMENTS OF INERTIA ****

! Быстрый подсчёт
/solu
irlf,0 ! отключить Inertia Relief, иначе не выводит МЦХ
/output,mass_output,txt
psolve,elform
/output
finish

!  *********** PRECISE MASS SUMMARY ***********
/solu
outpr,basic,all
irlf,-1
psolve,elform
psolve,elprep
irlist

[Fedor 1 597]

Кстати, Борман, с ансисом реализовать граничные элементы можно достаточно несложно. Сначала по модели из КЭ выделить только те грани элементов которые встречаются только один раз. Это и будут необходимые поверхностные. для модели. Дальше придется попрограммировать найти решение, уже для граничных элементов а потом посчитать перемещения и для внутренних узлов кэ и дальше можно рисовать уже средствами ансиса все :)    Можно посравнивать матрицы жесткости при редукции внутренних узлов КЭ и матрицы ГЭ . И вообще много чего любопытного можно поискать при желании :)

[soklakov 1 475]

Или можно конструировать элементы исключая внутренние узлы. Например кубики из пирамидок и что-то подобное. Единственно интересное это бесконечные границы например для моделирования полупространства. Но обычно этим не заморачиваются. Просто берут побольше земли и все дела. Но в этой книжке есть пример сколько примерно надо брать земли в конечных элементах чтобы получать как в граничных ... Это многим у нас не нравится :) Думаю что тихо увяло как везде :)

в HFSS при расчете высокочастотных электромагнитных полей помимо МКЭ может использоваться "интегральный метод моментов". Насколько удалось разобраться - все те же граничные элементы. Так что какое-никакое применение есть. Вот в задачи акустики примерить очень хочется, поскольку аналогия с антеннами на лицо.

Сначала по модели из КЭ выделить только те грани элементов которые встречаются только один раз. Это и будут необходимые поверхностные. для модели. Дальше придется попрограммировать найти решение, уже для граничных элементов а потом посчитать перемещения и для внутренних узлов кэ

А это больше суперэлементы напоминает. Не? Другое? 

[Fedor 1 597]

Суперэлементы или еще говорят редукция опирается на обычную технику конечных элементов, а граничные элементы на фундаментальные решения. Конечными элементами удобно моделировать конечные области, а граничными то, что вне их, то есть бесконечные. Для полей техника в самый раз. Где то у Зенкевича были описаны элементы с бесконечными границами вдоль части сторон, но в жизни такие не попадались для прочности. Какие-то были для полей в версии 5 или где-то рядом. Не знаю есть ли сейчас.   infinite - надо набрать в хелпах да поискать :)

[sapr3000 6]

а граничными то, что вне их, то есть бесконечные.

На самом деле полубесконечные.)

 

Конечными элементами удобно моделировать конечные области

Есть и полубесконечные конечные элементы. В научных работах Зенкевича встречались.)

 

Есть и аналитические конечные элементы и гибридные конечные элементы и т.д. и т.п..

Изменено 9 декабря 2014 пользователем sapr3000

[Fedor 1 597]

"На самом деле полубесконечные" нельзя быть полубеременной :)

Много всякой всячины народ придумывает, коль в моде тема. Гибриды всякие ежа и ужа. Как - то на конференции что-то говорили такое и даже табличку показали где энергия на графике почему-то прыгала, а не равномерно понижалась по мере разбиения.. Я сказал, что тогда должно существовать конечное разбиение при котором будет точное решение  представляемое бесконечными рядами для пластинки. Обещали прислать иностранные статьи где тоже такое есть, но до сих пор не прислали. :)

Сейчас как-то притихли искать новизну с актуальностью в этом направлении :)

[sapr3000 6]

Из 2-х или 100 или 1000 бесконечных элементов ничего не получится.

 

Бесконечность останется бесконечностью)

[Fedor 1 597]

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

Кроме того бывает потенциальной, а бывает актуальной :)

[sapr3000 6]

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

Энергия деформации растет по мере разбиения.)

[Fedor 1 597]

Растет, растет,  жесткость же понижается. Главное тут равномерность, монотонность. Там был график отклонения нормы от аналитического решения, вроде.  

[sapr3000 6]

Гибриды всякие ежа и ужа. Как - то на конференции что-то говорили такое и даже табличку показали где энергия на графике почему-то прыгала, а не равномерно понижалась по мере разбиения..

И в работах Зенкевича и в отечественной программе гибридные конечные элементы - совместные. Поэтому и сходятся монотонно.

[Fedor 1 597]

Что за гибриды. С производными что ли ? :)