FEA? VB?

[alleonid]

Уважаемые участники форума!
Был ли у кого-нить опыт реализации метода конечного-элемента на VBasic?

[uhm]

Конечный-Элемент - это двойная фамилия?
Ты бы написал, что за метод, или, хотя бы, из какой области...

[Viper]

У меня был опыт реализации конечно-элементного анализа для решения многопорных балок. В принципе реализовать FEA на VB вполне реально

[alleonid]

Метод конечного-элемента (Finite Element) это такой универсальный численный метод для решения задач теории упругости и так далее.
Вопрос следовало бы формулировать так - насколько удобен VB для решения математических задач (линейные уравнения и т.д.). Есть ли к примеру какие то специальные математические библиотеки под VB и так далее...
Обычно такие вещи делают на Фортран или С++.
Но я инженер а не профессиональный программист, и вряд ли мне стоит коллекционировать языки - один бы нормально изучить.

Конечный-Элемент - это двойная фамилия?
Ты бы написал, что за метод, или, хотя бы, из какой области...

[Viper]

Математика на VB вполне реализуема. И линейные уравнения и FEA и многое многое другое. Если есть алгоритм, то реализация его на любом языке вполне возможна. Насчет библиотек - поиск. Но лично мне не попадались. Все делал своими руками.

[alleonid]

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

Математика на VB вполне реализуема. И линейные уравнения и FEA и многое многое другое.

[v000v]

Реализация на VB вполне возможна... (сам не пробовал )

Альтернативный вариант:
1. Приобрести Ansys (~45000$ - лицензионный, на рынке можно найти гораздо дешевле )
2. Выучить встроеный в него APDL, позволяющий в числе прочего добавлять свои КЭ к существующим (>100).

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

Удачи!

[alleonid]

ANSYS я использую и APDL тоже юзаю (единственный язык который ИМХО знаю более или менее, даже дилом на этот счет имею).
К сожалению знания APDL недостаточно чтобы писать свои КЭ.
На эту тему имеется хорошее пособие Александра Душина, продается в CAD-FEM. К сожалению не владею фортраном, поэтому мне даже пособие вряд ли поможет.
Понятно что качества ANSYS мне не добится, я просто собираюсь приобрести опыт написания собственных программ по КЭ.
Пишет же Федор Пинежанинов...
Хочу быть похожим на него.

Реализация на VB вполне возможна... (сам не пробовал )

Альтернативный вариант:
1. Приобрести Ansys (~45000$ - лицензионный, на рынке можно найти гораздо дешевле )
2. Выучить встроеный в него APDL, позволяющий в числе прочего добавлять свои КЭ к существующим (>100).

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

Удачи!

[v000v]

APDL недостаточно чтобы писать свои КЭ

Я и не сказал, что его будет достаточно...

Фортраном я тоже не владею, однако мне кажется, что он не намного сложнее VB.

я просто собираюсь приобрести опыт

Вряд ли этот опыт может где-то пригодиться. (Мое мнение)

Пишет же Федор Пинежанинов...

У него есть свой сайт www.pinega....

Время потраченное на написание собственной программы, гораздо полезнее будет потратить на общение с Ansys'om на более высоком уровне (задачи Вам прийдется решать довольно близкие), и этот опыт может не только пригодиться но и окупиться...

Все вышесказанное, мое субъективное мнение, решать Вам...

[alleonid]

Да я не буду увлекаться
Месяц-два, простая программка - мсключительно для самообразования.
Погуляю, побалуюсь и вернусь к любиой жене...

Все вышесказанное, мое субъективное мнение, решать Вам...

[alleonid]

Вообще приятно найти в таком месте юзера ANSYS.
Это моя любимая программа. I have a dream. Когда я стану Абрамовичем - куплю ANSYS
Кстати
я орыскал немного, но так и не нашел
форума
равного по мощи VBStreets
и посвященного
Fortran
Может вы знаете?

Фортраном я тоже не владею, однако мне кажется, что он не намного сложнее VB.

[Viper]

Не хватает стандартных конечных элементов ANSYS? Их же там немеряно, на все случаи жизни хватает.

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

Мне лично хватает MSC.NASTRAN, а сам метод конечных элементов я реализовывал на VB для решения многоопорных балок. Ничего сложного.

[alibek]

А ху ис Ansys?
Система такого плана, как MatLab?

[Viper]

ANSYS - это программа конечно-элементного анализа. Позволяет производить расчеты конструкций на прочность (и прочие подобные дела). К MatLab не имеет никакого отношения

[alleonid]

Понятно, что это не сложно, и не язык главное.
Но одно дело знать в теории, а другое дело самому попробовать...
В ANSYS много конечных элементов
но не все, которые нужны.
К примеру ни в ANSYS ни в Femap
нет вантового элемента в том приближении, как его используют для расчета мачт. LINK 10 - не учитывает провисание, и пригоден только при большом натяге.
Все остальные варианты - моделируют оттяжку точным методом, что не всегда нужно.
Combine 39 задает нелинейную диаграмму Натяжение-удлинение только от начала координат.
Кстати в книге Шимковича есть пример башни с канатом и моделированние оттяжки у него не вполне корректно.
То бищь необходимость в пользовательских элементах всегда имеется.
Кстит как и в специализированных программах - типа SUDM для тех же мачт.

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

[alleonid]

Да, типа MatLAB но только - один численный метод, а именно метод Конечного элемента.

А ху ис Ansys?
Система такого плана, как MatLab?

[alleonid]

Viper!
а функции оперирующие с матрицами вы не писали, а если да, не могли бы подбросить?
Или вы не в общем виде реализовывали МКЭ для балок?
(то есть плясали не от общих зависимостей МКЭ а конкретно прописывали свой частный случай?)

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

[Viper]

Конкретно для балок была написана программа, которая делит балку на необходимое количество элементов постоянной жесткости, после чего все решается через КЭМ. Все это писалось на VB.NET.

Что касается матричных дел, то у меня есть набор классов, которые позволяют делать с матрицами практически все что необходимо (в том числе определители, системы линейных уравнений и так далее). Где то в мае, июне я скидывал это в форум. Можно поискать. Тема была что-то вроде "Метод Гаусса". Там правда есть пара неоптимизированных мест, щас я это дело поправил. Могу сбросить. Но новая версия на VB.NET.

[alleonid]

Не смог найти по запросу "Метод Гаусса" по www.vbstreets.ru ничего...
Не могли бы Вы сбросить (на VB)

Можно поискать. Тема была что-то вроде "Метод Гаусса".

[Viper]

скину на днях

[Viper]

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

[alleonid]

Большое спасибо.

[alleonid]

Существуют ли программки для перевода матрицы жесткости из формата Harwell-Boeing
в обычный наглядный вид?

[Пользователь_ПЭВМ]

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

Крамер совсем плохо решает, однако. Но так и должно по идее быть.
Вообще есть одна программа (мне известная), МКЭ, реализованная на VB - это Lisa (lisa-fet дотком). Считает очень медленно. Судя по всему это угасающий проект.
А у меня вопрос к вам - можно ли перенумеровать узлы в матрице жесткости, чтобы считалось быстрее?

[Viper]

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

Крамер совсем плохо решает, однако. Но так и должно по идее быть.

Ну это само собой. Насколько помню там должны быть и методы по Гауссу. Если вдруг нет, то в "Кирпичах" лежит архив, где точно есть все.

Вообще есть одна программа (мне известная), МКЭ, реализованная на VB - это Lisa (lisa-fet дотком). Считает очень медленно. Судя по всему это угасающий проект.
А у меня вопрос к вам - можно ли перенумеровать узлы в матрице жесткости, чтобы считалось быстрее?

Это к кому вопрос то? Перенумеровать то можно, и будет быстрее.

[Пользователь_ПЭВМ]

. Насколько помню там должны быть и методы по Гауссу. Если вдруг нет, то в "Кирпичах" лежит архив, где точно есть все..

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

[Viper]

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