Firefighter

Сравнение сеток

46 сообщений в этой теме

Уважаемые профессионалы в области КЭ анализа!

Давно интересует Ваше мнение по поводу -><noindex>Сравнение сеток</noindex>.

Больше всего поражает разница в напряжениях (~30% - и это только от выбора формы элемента!).

Объясните, пожалуйста, может ли такое быть, и чего тогда стоят все эти автоматические генераторы сеток?

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах


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

В области концентрации напряжений всего пара элементов и много на однородных по структуре областях. Хороший пример неудачного разбиения.

Всегда лучше использовать квадратичные элементы. Кубики в принципе имеют большие возможности для приспособления к полю перемещений чем пирамидки. Так как интерполирующие полиномы имеют больше членов (бедный Василий Иванович).

Такой разброс по напряжениям вполне может быть в приведенной ситуации.

Для преодоления этого придуманы сгущения в интересующих областях. :rolleyes:

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Насчет сгущений - верно замечено. Вообще корень зла во всех автоматических мешерах в том, что алгоритмы, применяемые для разбиения на сетку, понятия не имеют, ЧТО именно Вас интересует в детальке, вдобавок разбиение достаточно "хаотично". Т.е. большинство мешеров будет одинаково бить как концентраторы напряжений, так и относительно "спокойные" области модели. Попробуйте побить в ANSYS элементарный объем coarse free mesh'ем, а потом "нормальным" 4-sided, разница очевидна...

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

Поэтому FEDORовские полиномы, интегрирования, кубики, пирамидки, интерполяции и члены - это несколько из другой оперы.

2fedor:

>[..]Так как интерполирующие полиномы имеют больше членов (бедный Василий Иванович)[..]

- Василий Иванович - это Ваше настоящее имя?

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

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

Да все из той же о конечных элементах. Если Вы имеете функции в объеме способные хорошо приспосабливаться, то в принципе безразлично, как разбиваете. Пирамидка обычно имеет одну точку интегрирования для линейных функций, поэтому осредняет грубо. Квадратичный - уже заметно лучше. Посмотрите обложку книжки Сегерлинда.

В пределе при методе Галеркина вообще достаточно одного элемента, только он будет иметь сложную геометрию, и много координатных функций. Обычноя математика, однако. Василий Иванович - это народный герой, который завалился на математике, так как не мог представить квадратный трехчлен. www.pinega.da.ru :rolleyes:

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

по представленным фрагментам можно сказать следующее:

1. Очень похоже, что работали в Патране

2. Тетраэдры следовало использовать 10 узловые, т.е. со срединным узлом. Тогда можно говорить о сравнении результатов.

3. "FEDORовские полиномы, интегрирования, кубики, пирамидки, интерполяции и члены" как раз из оперы, и слова из песни не выкинешь

4. Данный пример наглядно демонстрирует, от чего существует предубеждение (а равно заблуждение) по тетраэдным элементам.

5. Цитата с сайта: "данная TETRA-модель не позволяет определить истинный уровень концентрации напряжений в локальной зоне". Ну, что тут сказать - правильный вывод :smile: . И автомат можно зарядить холостыми :doh: ....

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Тетраэдры следовало использовать 10 узловые, т.е. со срединным узлом. Тогда можно говорить о сравнении результатов.

Говорить-то можно, а сравнивать как?

Какой результат точнее ?

Какова погрешность решения ?

Существует ли метод численной оценки погрешности ?

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

Ответов на эти вопросы руководства к программным продуктам практически не содержат.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

lks

Вроде бы эта тема затрагивалась на форуме. И на сайте Fedora статья выложена...

А так - количественный способ оценки погрешности, обусловленный сеткой, состоит в сопоставлении усредненных узловых результатов и значений в точках интегрирования. В ANSYS'e это команды PLNS и PLES, в продуктах MSC.Patran/Nastran, /Mark есть подобные инструменты. Есть еще субъективный подход - свой опыт и свое понимание физического явления (Типа: необычные "картинки" результатов - чаще всего неверные результаты)

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Начинаем с одного, переходим на другое. Господа, вопрос был о генерации сетки, а не об элементах (в конце вопроса была фраза- "чего тогда стоят все эти автоматические генераторы сеток?". Вот если бы там было - "чего тогда стоят кубики с пирамидками?" - или - "что лучше пирамидки или кубики?" - тогда - другое дело). Если ошибки в расчетах появляются только из-за тех или иных элементов, то какой смысл каждый раз советовать друг другу бить модель вручную и сетовать на тупость мешера в той или иной программе.

Так что, когда у Вас спрашивают, как пройти на улицу Горького, не советуйте брать такси до Тверской, - покажите человеку куда идти. А советы почитать Сегерлинда... Кто там еще на книжной полке у Fedor'a? Зенкевич, Галлагер? У которых, между прочим, информации об алгоритмах генерации сетки - ноль...

Т.е. я не против советов вроде "почитайте про методы Delaunay, Point insertion, Advancing Front, Structured и Unstructured Hex Meshing", - они, на мой взгляд, были бы уместны.

>"В пределе при методе Галеркина вообще достаточно одного элемента, только он будет иметь сложную геометрию, и много координатных функций. Обычноя математика, однако."

- Уважаемый Fedor, покажите мне хоть одну программу, которая таким способом РЕАЛЬНО рассчитывает РЕАЛЬНЫЕ модели и дает РЕАЛЬНЫЕ результаты. Вот это вот Ваше упоминание "обыкновенной математики" применительно к данной ситуации - не может не настораживать. Вы ведь вроде написали некий код, связанный с МКЭ. Так почему же Вы так примитизируете? Возьмите Зенкевича (Галлагера, Сегерлинда) и найдите там что-нибудь, связанное, к примеру, с т.н. hourglassing modes. Нашли? Нет? Все так же советуете всем читать их дальше?

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

еще один м-а-а-а-ленький вопрос: то ли с движком форума что-то не то, то ли просто совпадение...Но вот странные цифры:

Fedor Дата Jul 7 2004, 7:38

asds Дата Jul 7 2004, 7:58

lks Дата Jul 7 2004, 8:14

asds Дата Jul 7 2004, 8:36

Это, как вы понимаете, время, когда поступали реплики от перечисленных участников форума. Странное ощущение, что люди собираются пораньше с утра обсудить мой ответ вместо того, чтобы чистить зубы, идти на работу, и пр. и пр. Либо у г. Пинежанинова наблюдается размножение личности. Примечательно, что все три участника довольно единодушны в высказываниях, а последний (как это обожает делать Fedor) дает ссылочку на страничку опять таки на федоровскую страничку. Стоит также отметить, что обычно люди подписываются более-менее осмысленными никами, а не инициалами типа lks, asds. Странно, что с легкой руки г.Пинежанинова не появились персонажи qwerty, zxc, etc.

Кстати, Fedor, на кадфемовском форуме Вы спрашивали об отключении вывода графики на экран во время геометрических построений. Ответ написан. Это к вопросу о чтении книг и прочей документации.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

PS: я могу ошибаться. В этом случае приношу свои ГЛУБОЧАЙШИЕ извинения господам asds и lks.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

"Больше всего поражает разница в напряжениях (~30% - и это только от выбора формы элемента!)" - вроде это был основной вопрос. :rolleyes:

Ага, заговор устроили. :doh:

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Galitsky

1. По поводу времени - обратитесь к админу. Сейчас московское время 9:16. Работа начинается у нас в 08:30 (Нахожусь я в Питере). И удобнее мне с утречка посмотреть свежую информацию и новости (нравиться мне этот форум :smile: )

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

3. Что касается "совпадения во взглядах" - бывает... Опять, чтоль поворять - были бы это КЭ второго порядка смотрел бы результаты дальше, а так... Ну получили 30% отклонений - и чего Вы этим хотите сказать...

4. По поводу ника - не придал я ему значения, забасал удобные для себя клавиши.

Чистый флейм у нас дальше пошел. :poster_offtopic:

P.S. Если КАЖЕТСЯ - креститься надо :smile: (Была у меня соседка по комуналке - все ей странности казались. эта фраза оттуда)

P.P.S Извинения приняты.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

"РЕАЛЬНО рассчитывает РЕАЛЬНЫЕ модели" по сути вся технология подконструкций (суперэлементов) и т.д. :rolleyes:

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

"hourglassing modes" , а это что за зверь, по русски нельзя обозвать? :g:

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

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

По-моему форум у нас, в силу географического положения, сидит не на московском времени, а на Екатеринбургском, если есть такое :smile: Одним слово на час-два от Москвы сдвинут

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

"Так почему же Вы так примитизируете? " Помните как у Штирлица: "потому и без оружия" :rolleyes:

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Начинали с того, что вопрос о генерации сетки, а не о форме элементов, а теперь уже обсуждаем в какое время пришли сообщения, да у кого какой ник. Классический русский разговор на кухне.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Кстати, время сдвинуто ровно на час.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Да, случайности случаются :)... И не заговор, а простая шизофрения мне почудилась у одного из участников форума. Тьфу-тьфу, вроде все в порядке оказалось.

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

>В пределе при методе Галеркина вообще достаточно одного элемента, только он будет иметь сложную геометрию, и много координатных функций.

>"РЕАЛЬНО рассчитывает РЕАЛЬНЫЕ модели" по сути вся технология подконструкций (суперэлементов) и т.д. технология подконструкций (суперэлементов)

- Возвращаю Вас, уважаемый Fedor, к теме генерации сетки. Т.е. по-Вашему выходит, что мешер в данном случае не нужен? Задал элемент и ладно?

Substructuring is a procedure that condenses a group of finite elements into one element represented as a matrix. А узлы этой самой матрицы-суперэлемента-подконструкции как, руками в файл набивать?

Я к тому, что мало-мало инструкцию к программе читал, про суперэлементы слышал вообще-то. Только вот подменять понятия НЕ НАДО, Fedor.

Равно как и цитировать с умным видом обложку (!) книги Сегерлинда.

Про hourglassing modes. Открываете все тот же мануал ANSYS (ох уж этот RTFM, Fedor!). Набиваете в строке поиска HOURGLASSING. Нажимаете на клавишу Enter. Выбираете в появившемся списке строку Chapter 9. Hourglassing. Справа, если Ваш Windows еще не завис, появится соответствующий текст. Аккуратно выделяете текст мышкой (для этого надо установить курсор в начало первой строки текста, нажать левую кнопку мышки и, не отжимая кнопку, перевести курсор в правый нижний угол. После чего нажимаете сочетание клавиш Ctrl-C, запускаете программу-переводчик (например, Promt), выбираете направление перевода (ОБЯЗАТЕЛЬНО выберите "англо-русский"), нажимаете клавиши Ctrl-V и выбираете в главном меню программы-переводчика пункт "перевести весь текст". В результате работы программы, должен появится перевод интересующего Вас текста.

Если у Вас, Fedor, по каким-либо причинам нет документа "Ansys x.x Documentation" или не установлена программа-переводчик или Вы просто не владеете иностранными языками, сообщаю, что в русскоязычной переводной литературе термин hourglassing по моим данным отсутствует. А речь идет о возникновении ошибок при расчете, вследствие применения математических методов (читайте учебные пособия по высшей математике, книги по методу конечных элементов, вообще, материалы касающиеся повышения скорости вычислений на ЭВМ, - Вы ведь, я думаю, не HP-машину используете...), которые приводят к состояниям, лишенным физического смысла (читайте статьи авторов Belytschko, Liu, Flanagan и других).

Для Вас термин hourglassing - новый? А вот теоретические мануалы к ЛС-Дайне, Абакусу и пр. его нередко поминают. В Ansys Documentation он по некоторым причинам мало встречается, но в процессе расчетов часто используется. Самому термину более 20 лет. Делайте выводы об актуальности своих знаний.

RTFM

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

"возникновении ошибок при расчете" - вопрос вообще-то не новый, смотришь на число обусловленности и так все ясно. Скучно даже обсуждать. И других подходов тьма, если говорить об ошибках вычислений. При 64 разрядах обычно все в порядке.

Нас интересует ошибка аппроксимаций (лемма Сеа у Сьярле вроде). Их обычно оценивают по сравнению внутренней энергии и работы внешних сил, или локально по перепаду напряжений на границах элементов. У Ansys вроде тоже что-то в этом духе.

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

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

За объяснения как работают клавиши, спасибо.

Вы наверное из педагогов, они обычно любят тривиальные вещи обьяснять.

Ну да об этом лучше Федора Сологуба не напишешь. :rolleyes:

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

"ЛС-Дайне", так вот это откуда. Года 2-3 назад смотрел, но после того как увидел как варварски осуществляется интегрирование, потерял интерес. Попытки использовать "скользящий полином" для исправления плохих решений или измерений или "цифровая фильтрация" давно используются.

Еще один очень грамотный специалист говорил, что теорема Вариньена из теор. меха у них вроде нарушается. Это не добавило интереса.

В последней книжке Зенкевича хорошо показано к чему приводит дефицит точек интегрирования. Впрочем об этом уже писалось на сайте.

О "одном элементе" полистайте, например, Михлина. Раньше вообще функционалы минимизировались с помощью функций определенных сразу по всей области, это уж потом Курант предложил работать на подмножествах для технологической простоты вычислений. Еще раз повторю "~30% - и это только от выбора формы элемента" -это представляет интерес. Все остальное мелочи. Сеточки, автоматизация разбиения и пр. :rolleyes:

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Насчет шизо полностью согласен. Часто встречается в переходные периоды. Раскол он и в сознании раскол.

Многие, например, считают, что обычные прямые методы вариационного исчисления чем-то принципиальным отличаются от мкэ, кроме как специфической формой координатных| базисных|формы функций.

Что решение систем, как-то зависит от того как их сформировали.

Что программы не просто проекция математических идей на "железо".

В такое время живем. :rolleyes:

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Ну что, Fedor? Опять в слова играемся?

>Вы наверное из педагогов, они обычно любят тривиальные вещи обьяснять.

-По поводу объяснения работы клавиш: это был недошедший до Ленинграда (Вы ведь вроде ленинградец?) по линиям связи сарказм. Жаль, что Вы этого не поняли. Честное слово, первый раз встречаю человека, который упорно не желает использовать примитивнейшие возможности компьютера. Как же Вы ищете информацию в сети? Ах да, заходите на форум (лучше сразу на несколько) и ждете когда Вам ответят... Вы многого можете не узнать, если не научитесь "тривиальным вещам". И касается это не только сети.

>И других подходов тьма, если говорить об ошибках вычислений. При 64 разрядах обычно все в порядке.

(Насчет других подходов - могли бы назвать пару, чтобы я знал, что имеется здесь в виду). А Вы используете 64 разрядную систему? Ого! И дома у Вас 64-разрядка? Ух ты! И у всех Ваших знакомых тоже? Ну и ну! Так почему же Вы, Fedor, обсуждаете такие скучные вещи, как работа с wintel'ом, пишите какие-то программулины под нее? А потом, 64 разряда - ведь тоже не так много. Там, где вчера хватало нескольких миллионов элементов, завтра не будет хватать миллиардов. Что попишешь, прогресс.

>после того как увидел как варварски осуществляется интегрирование, потерял интерес.[про LS Dyna]

Fedor, - это не варварство. Это прмерно то, что принято называть оптимизацией (подхода, метода расчета, машинного кода и т.д.). One-point integration - (это варварство?) - используется не только в дайне (вернее сказать, используется и в ней тоже). Стыдно, Fedor. Если бы Вы сделали хоть одно телодвижение в сторону перечисленной мной литературы, материалов, а не вспоминали каждый раз выученных наизусть Сегерлинда и Норри, Вы бы вряд ли себе позволили прохаживаться по поводу дайны (почему абакус забыли? ведь я его упомянул. Смелее, Fedor, там ведь тоже варварский подход!). Вы серьезно считаете себя специалистом уровня Flanagan'a, Belytschko (которые посвятивили исследованию этих проблем известное количество усилий и времени), что говорите о варварском подходе?

>...один очень грамотный специалист говорил.... у них вроде нарушается...

...один... вроде... Вы к точным наукам вообще отношения не должны иметь. Выясните поточнее, потом говорите вволю. Кстати, LS Dyna тоже не совсем неграмотными писалась.

>Еще раз повторю "~30% - и это только от выбора формы элемента" -это представляет интерес.

Сравнивать "кубики" с "пирамидками" - Вам ЭТО интересно? Это даже не вчера, это, - я так думаю, даже не прошлый век. Это - из справочника пользователя Ansys, который по своему определению на полноту и фундаментальность не претендует.

>Все остальное мелочи. Сеточки, автоматизация разбиения и пр.

На этих мелочах люди миллионы делают. Если не знаете, смотрите словосочетание ICEM CFD, ParaMesh etc. Кстати, по приведенному в самом начале разговора адресу я увидел прежде всего неграмотную сетку. То, что там говорят про недостатки элементов, исходя из некорректной постановки вопроса,- так на рекламных щитах тоже много как бы правдивой информации, помеченной звездочками и мелкими-мелкими сносками.

>Часто встречается в переходные периоды.

Оптимист. Уже дети успели вырасти за это время. А Вы все думаете о "переходном периоде". У нас в науке "переходный период" продлится еще лет... - не знаю, Fedor! Сколько Вы дадите отечественной фундаментальной науке, чтобы она вылезла оттуда, где сейчас находится? и насколько это вообще возможно? и насколько хватит советского запаса?

>Многие, например, считают,

Я дал Вам повод так думать о моей скромной персоне? Укажите где.

>Что программы не просто проекция математических идей на "железо".

Именно потому, Fedor, и появились программисты. Вам, верно, и невдомек, что Ваше "просто" - СОВСЕМ непросто (извините за плохой каламбур... и не вспоминайте, пожалуйста, поручика Ржевского из своей подборки бородатых анекдотов...). Математическая идея, как правило, претерпевает огромные видоизменения для успешной реализации в коде. Отсюда и one-point integration, и вытекающий из него hourglass control и многое, поверьте мне, многое другое.

Хотя... Если рассчитывать модели кубиков и пирамидок, то применяйте full integration и программируйте свои элементы хоть в Excel'e. Такой подход свойствен людям, не понимающим, что метод расчета - ничто без самого расчета - "Как же, ведь все по книжке Сегерлинда!", "Самые качественные полиномы, самые гладкие интерполяции!". ЧТО ДЕЛАТЬ С ПРОГРАММОЙ, В КОТОРОЙ МОЖНО СЧИТАТЬ ТОЛЬКО ОДНОЭЛЕМЕНТНЫЕ МОДЕЛИ? Я ответ знаю. А Вы?

Аббревиатура RTFM - появилась неслучайно. Спрашивать про давно известные всем вещи, которые к тому же находятся в свободном доступе - руку протяни и читай? К тому же еще и пытаться полемизировать, не обладая определенными знаниями в данной области?

Зачем Вы это делаете, Fedor? Даете пример, как делать НЕ НАДО? Или создаете себе ауру непонятого гения? (как Вы пустили реплику на одном из форумов, - что-то вроде - "ну вот, пришел большой знаток, и никому сказать против нечего", - цитирую не дословно, но смысл фразы хорошо запомнился).

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Уважаемые специалисты!

С большим вниманием и легкой грустью наблюдаю за развернувшейся дискуссией. Бесспорно, что в спорах рождается Истина, но нельзя же так... :bash:

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

Ребята! . . . ДАВАЙТЕ ЖИТЬ ДРУЖНО!

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Согласен с Firefighter полностью.

Цепляться к словам - умеют многие; лучше заниматься делом, а не словесными перепалками. Свое многословие объясняю прежде всего тем, что Fedor частенько позволяет себе советовать почитать классиков, причем часто не к месту, причем часто уже читавшим их людям. А вопросы, поднимаемые им на форумах - нередко действительно интересны и охватывают практически все области, где только можно применить МКЭ (иногда диву даешься КАК это человек успевает одновременно заниматься вопросами грунтов, моделями горения, прочностью, гидрогазодинамики и пр.).

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Готов ответить на поставленные Вами вопросы по существу, а не "расплываясь мыслью по древу", но где-нибудь во флейме. Заведите топик если хотите и почетче сформулируйте вопросы. :rolleyes:

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

" люди миллионы делают" на спичках тоже.

На многие вопросы Вы найдете ответы почитав предыдущие обсуждения. :rolleyes:

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Добрый день коллеги,

мне кажется вы немного ушли от темы, но на вопрос

по моему так и не ответили. Постараюсь конкретизировать

вопрос не опираясь на действительно безполезные картинки.

1. Какое количество элементов => узлов => DOF обычно

требуется для моделирования чуть более сложной модели для вашей области.

(вопрос расплывчатый, но для фирм где МКЭ практикуют есть вполне реальный ответ для 32 и 64 платформ. Сам встречал. Для h элементов, советую приводить примеры для элементов со средним узлом. Без среднего узла рекомендуется только для контактных элементов и то не для всех решателей, напр. ABAQUS)

- для h ?????

- для p ?????

элементов соответственно.

2. Сколько времени пройдет счет для этого количества элементов?

И сколько времени нужно для построения сетки?

Будет очень хорошо если приведете название кода в котором это делается.

Можно также привести примеры для явного кода для сравнения, с коротким описанием задачи: кол элементов, DOF, шаг по времени, кол. шагов и время счета.

3. По поводу тетраэдров (со средним узлом). На данный момент кроме слов что: "это плохо" меня ни кто не убедил. Я вижу только следующие проблемы:

-Нужно большее количество элементов для качественного разбиения и кубики лучше моделируют пов. контакта.

НО!!!!!!

Если задача влезает в память или почти влезает.

ТО!!!!!!

зачастую строить хорошую сетку несколько дней и решать ее 0,5 дня хуже, чем. Построить сетку за 0,5 дня и решать ее несколько дней. Иногда да, нужно отлаживать решение и время счета более критично, но если вам после решения придеться перестроить сетку, то это мягко говоря неприятно.

4. Ошибку решения можно оценить и по разности в элементном решении для 1 узла и по градиенту внутри элемента. Если кто то готов поделиться цифрами очень буду признателен.

5. Более интересно описание влияния неидеальности формы на точность решения. К стати это касается H-разбиения на наших рисунках. Если кто то готов поделиться цифрами очень буду признателен.

PavelN

PS

Galitsky - спасибо за приток свежей воды в этот форум. Много мутных направлений я здесь слышал по поводу целостности функционала, а мне приходиться решать реальные задачи. И зачастую мы все работаем до первой поломки того что мы считаем, т.к. эксперименты и испытания давно не проводятся. Есть множество вопросов на которые прсто не найти ответов.

Но все равно, будте осторожны, не всем наравится революция, даже во взглядах. И очень может быть, что те кто пишет ответ в двух сообщениях с интервалом в 5 мин. может быть более любим, чем борцы за правду.

Удачи Всем.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Свою точку зрения я изложил в статье об интегрировании. Добавить нечего, взгляды не изменились.

Придерживаюсь мнения Зенкевича, что и нескольких хороших хватает. :rolleyes:

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

to FEDOR:

по существу и коротко:

Вопрос не я задавал, а Firefighter. Вопрос неплохой. Причина возникновения вопроса - "специалисты". Просто "специалистов" много в последнее время.

Разговаривать с Вами или, тем более, задавать Вам вопросы (если мне нужна будет Ваша помощь или Ваши советы... - Нет, такой вариант я лучше не буду рассматривать) во ФЛЕЙМЕ? Тем более не "расплываясь мыслью по древу"?! я-то, грешный раньше думал, что флейм, в отличие от дискуссии, подразумевает беспочвенные обвинения и словопрения типа "кто выше пИсает". Вам, наверно, виднее. Боюсь догадываться даже, какие аргументы Вы приводите "во флеймах", если в "нормальных" топиках прибегаете к подмене понятий и облыжным обвинениям.

С Зенкевичем и Ко - спорить никогда и в мыслях не было. Просто уповать на то, что в Зенкевиче есть ВСЕ - достаточно нелепо (я, например, не атеист, но не буду искать в Библии ответ на вопрос, какими элементами лучше бить ту или иную модель, хотя и существует мнение, что в Библии есть ответы на ВСЕ вопросы ).

Прошу заметить, что поправлять Вас цитатами из мануалов (инструкций по эксплуатации) к ANSYS, LS Dyna и Abaqus - могу не только я (времени жалко, да и мало его), а любой желающий, знакомый с английским языком.

Считаю, что запудривание мозговых извилин как начинающим свое знакомство с МКЭ, так и уже вполне опытным людям, цитированием классических трудов к месту и не к месту, - довольно вредное занятие. А постоянные ссылки на www.pinega.da.ru - считаю прямой саморекламой в неясных мне целях.

to Caepavel:

- Точной статистики по количеству не вел. Обычно это были модели порядка 20-50 тыс. элементов, причем ограничения накладывал, естественно, хард. Начинал (как и многие) - с плоских задач, - 1-10 тыс. элементов.

Сейчас, по всей видимости, предстоит работа над более "солидными" 3D моделями - от 100 тыс. до млн. элементов (в зависимости от выбранного подхода к моделированию). Основная неприятность (помимо поиска все того же харда) - введение своих моделей поведения материала (т.е. очередные траблы с UPF). Сколько времени будет занимать расчет таких моделей мне по понятным причинам сказать пока затруднительно. "Неидеальность формы" или лучше - "своеобразие" сетки, - связана в данном случае непосредственно с самим методом моделирования и разбиения (штатный мешер ansys не будет использоваться).

По поводу применения тех или иных элементов - крайне простой ответ - у всех есть свои плюсы и свои недостатки. В противном случае вместо пары сотен типов мы бы имели в ANSYS список из нескольких элементов, что-то вроде "... coupled membrane, coupled shell, coupled solid.." и т.д. Смотрю на быстродействие, на сходимость, на правдивость результатов на части модели с разными элементами. Честно говоря, достаточно эмпирический подход.

Спасибо Вам за теплые слова.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

"мутных направлений я здесь слышал по поводу целостности функционала"

Мне все одно не поверят, поэтому несколько цитат.

"В действительности, как может заметить читатель, рассмотренные в предыдущей главе процессы являются частным случаем мкэ, когда используется только один элемент..."

" Если базисные функции определены кусочным образом, то выгодно поставить им с соответствие некоторый малый носитель, полагая их равными нулю всюду, кроме рассматриваемого элемента"

c.94 Зенкевич, Морган. Конечные элементы и аппроксимация.

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

Тогда можно перейти к матрицам жесткости элемента и легко строиться каждая по всей области, " эта неэффективность объясняется тем, что матрица жесткости отдельного элемента ... имеет такое же число строк и столбцов, что и глобальная матрица. Как видно из формул большинство коэф. в матрице элемента равно нулю " с 105 Сегерлинд

Т. е. матрица жесткости элементов строится каждая по всей области и складываются. Далее уже описываются приемы позволяющие не хранить нули.

"... Поэтому методы подобного типа также наз. вариационно-разностными, проекционно-разностными, методом Галеркина. Проекционно-сеточные методы допускают и подпространства в которых отдельные базисные функции могут и не иметь локального носителя " с 283 Математический энкиклопедический словарь.

Думаю в комментариях не нуждается.

По поводу интегрирования, стараюсь быть законопослушным гражданином страны "математика" и не собираюсь бунтовать против ее короля. :rolleyes:

Участникам дискуссии большое спасибо, давно так не веселился. :rolleyes:

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

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

Саморекламы нет, просто хотел показать затурканным "ученым", что можно не боятся подсчетчиков времени и свободно излагать свои взгляды не только на "коммунальной кухне". Что у нас свободная страна. Впрочем это всегда личное, я себя скованным и раньше не ощущал, просто "мух и жучков" побольше было и позловреднее чем сейчас. Парткомы всякие "облико-морале".

Кафедральные посилелки, нишки. Грызня за деньги. В общем ничего оригинального.

Спасибо за разыгранную комедию.

Удачи

:smile:

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

"флейм, в отличие от дискуссии, подразумевает беспочвенные обвинения и словопрения типа "кто выше пИсает" ну, я это понимаю иначе.

Просто свободное обсуждение, как в бане "ни святых ни генералов". Ну и как классик учил "поэтом можешь ты не быть, а гражданином быть обязан".

Обсуждаем же не личности, а идеи.

Элементы коммунальной кухни иногда возникают, но всегда найдутся люди, которые поправят. Впрочем не хотите, так не хотите. У нас свободная страна. :rolleyes:

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Вот Ваш замечательный полив, Федор Пинежанинов:

>Саморекламы нет, просто хотел показать затурканным "ученым", что можно не боятся подсчетчиков времени и свободно излагать свои взгляды не только на "коммунальной кухне". Что у нас свободная страна. Впрочем это всегда личное, я себя скованным и раньше не ощущал, просто "мух и жучков" побольше было и позловреднее чем сейчас. Парткомы всякие "облико-морале". Кафедральные посилелки, нишки. Грызня за деньги. В общем ничего оригинального.Спасибо за разыгранную комедию.Удачи

А вот мой ответ:

Большего маразма мне не доводилось слышать ни от одного выжившего из ума инженера. Были конечно экземпляры, но Вы, Федор, превзошли ВСЕХ. Вы пытались обвинить меня, - мол, я люблю оскорблять инженеров, ненавижу фортран и классиков МКЭ, прочие мелочи (?) опустим. И вообще, - иди, мол, читай Зенкевича, Сегерлинда и математическую энциклопедию, пацан.

Спасибо, Fedor. Читал. Только не вижу необходимости читать Зенкевича в Вашем изложении. Или в Вашем (довольно посредственном, кстати) пересказе на PINEGA.DA.RU. Очередные Ваши цитаты - Вы думаете - к месту?

Могу Вам посоветовать не увлекаться Зенкевичем, а начать с более усвояемого произведения, например, "Курочки-Рябы", в котором Вам бы было все понятно.

Свободная страна? Вы в Швейцарии живете, Fedor? Или СПб уже отделился от РФ? А если она такая свободная, что ж Вы так шугаетесь "подсчетчиков времени" (лучше уж подсчитать количество воды, которое Вы вылили на форумах). И почему до сих пор поминаете жучков, парткомы? Больше вспомнить нечего? Или Вы у нас великий диссидент?

Грызня за деньги? - А вот это уже довольно глупое высказывание, г.Пинежанинов. С кем я тут грызусь за деньги, скажите, не стесняйтесь. С Caepavel? может быть с FireFighter? Или с lks? asds? Точно также я могу сказать, что Ваши речи оплачены ЦРУ, МИ6 или Моссадом, в целях разваливания российской науки изнутри.

Потом, прозвучал конкретный вопрос, - СКОЛЬКО ВЫ ЭЛЕМЕНТОВ ИСПОЛЬЗУЕТЕ? Вы ответили? - Нет. Проще в очередной раз сослаться на свою страничку PINEGA.DA.RU. Там есть указанная информация? Я не нашел.

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

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Galitsky

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

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

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Galitsky , Fedor

Коллеги! Вы явлетесь вполне квалифицированными людьми и, если бы вам пришлось совместно решать конкретную проблему, то, я думаю, вы легко бы договорились. Зачем же изливать друг на друга в таких количествах потоки желчи, уходя все дальше и дальше от обсуждаемой темы. Если вам приятно ругаться, вы иожете сделать это с помощью e-mail. Данный форум -это место, где люди делятся опытом и знаниями, а не отрицательными эмоциями.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

"легко бы договорились"

"это вряд ли" :rolleyes:

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

to lks:

-Оценить до расчета сходимость - для простых, известных задач - можно сравнивать, анализировать. Для сложных - невозможно. Меня не устроил бы ответ - "да, наверно, сходимость будет ничего", а методов точной оценки - не существует по определению, - здесь мы возвращаемся к вопросу о "плохих" и "хороших" элементах. Почему "плохие" элементы используются? Да потому что в первую очередь это выгодно по объему вычислений. Почему "хорошие" используются? - Ответы очевидны.

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

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

______________________________________________________________

2 Fedor:

Коллегой я Вас никогда не считал, рассуждать с Вами далее - не вижу никакого смысла, ввиду полного отсутствия у Вас своих идей (до сих пор с трудом заставляю себя оставить сообщение, когда "попадаю" на автоответчик; с Вами - ситуация еще хуже). Как программист Вы ничего из себя не представляете, - набить пару матричных операций в автоматически генерируемом MFC коде - может и ребенок. Не давайте глупых и вредных советов начинающим, - что посеешь, то и пожнешь. Поищите в интернете другие способы добавлять посетителей на свой сайт, т.к. то, чем Вы сейчас занимаетесь - дурной вкус.

До свидания, Fedor.

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

"никогда не считал" и правильно делали :rolleyes:

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

NO COMMENTS

Поделиться сообщением


Ссылка на сообщение
Поделиться на других сайтах

Создайте аккаунт или войдите для комментирования

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

Создать аккаунт

Зарегистрируйтесь для получения аккаунта. Это просто!


Зарегистрировать аккаунт

Войти

Уже зарегистрированы? Войдите здесь.


Войти сейчас

  • Сейчас на странице   0 пользователей

    Нет пользователей, просматривающих эту страницу