Перейти к публикации

Рекомендованные сообщения

Всем привет.

Кто хочет немножко поиграть, детская головоломка.

Подсказка: Используйте метод конечных разностей.

P.S. Да это не метод конечных элементов, но для начала пойдет, это намного проще.

Игра.png

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


UnPinned posts
46 минут назад, VICT0R_1945 сказал:

Если учесть простоту метода, то 75% точностью вполне можно довольствоваться в элементарных задачах (в задачах с несложной геометрией).

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

 

19 минут назад, VICT0R_1945 сказал:

вроде как есть задачи, которые решались еще в СССР с помощью ЭВМ, но к сожалению, утверждать ничего не буду

Может тогда это и было оправдано в условиях ограниченных компьютерных мощностей когда выбор был либо посчитать с погрешностью в 20% либо задача вообще не влезет в память. Но сейчас это интересно только в плане обучения.

 

Сейчас еще посмотрю какие температуры я получу у себя интерполирую результаты средствами Femap.

UPD. Попробовал - вышло не очень, точность плюс-минус лапоть, не буду результаты даже приводить.

Изменено пользователем karachun
Ссылка на сообщение
Поделиться на других сайтах

Матлаб вроде можно заставить решать такие уравнения на сетке. У него как то проще с этим дожн быть чем в ансисе.

Там вроде даже дэмка была встроенная.

Ссылка на сообщение
Поделиться на других сайтах
21 минуту назад, karachun сказал:

Может тогда это и было оправдано в условиях ограниченных компьютерных мощностей когда выбор был либо посчитать с погрешностью в 20% либо задача вообще не влезет в память.

Эта задача если не влезала) в оперативную память, то включалась виртуальная память. Но в этой задаче 1 степень свободы в узле. Проблем никаких нет.

Эта задача зависит от угловых температур. Проверьте.

Ссылка на сообщение
Поделиться на других сайтах
44 минуты назад, karachun сказал:

@VICT0R_1945 У меня получилось то же самое.

1.png

Но разница есть и даже на простом квадрате она доходит до 20%.

2.png

А теперь подумайте что вы там наинтерполировали на своем поршне?

Ну подумал уже)) Так там еще все хуже... Там распределение вообще по 1 оси, такое было бы там в разы лучше))

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

@karachun Кстати, если говорить о ПК-расчетах и ручном, то выходит значения у интерполяции (как вы ее называете, хотя по мне это название неверное, все же это аналог разностного метода Эйлера, по сути первооткрывателя вообще всех разностных методов) даже больше, что безусловно только +, придется упрочнять конструкцию, если говорить о точности это одно, а если о запасе прочности - то совершенно другое, выходит, что так может даже и надежнее, этот краевой камень все же нельзя просто отодвинуть и плюнуть на него без внимания.

Вообщем пока не доказано обратное (что при определенных границах или ГУ температуры ниже расчетных на ПК), как вы и сказАли это был бы неплохой материал для изучения, во всяком случае это имхо.

Изменено пользователем VICT0R_1945
Ссылка на сообщение
Поделиться на других сайтах
19 минут назад, VICT0R_1945 сказал:

(что при определенных границах или ГУ температуры ниже расчетных на ПК),

Смените знаки всех температур на границах :)

Ссылка на сообщение
Поделиться на других сайтах
1 минуту назад, Борман сказал:

Смените знаки всех температур на границах :)

Предположим, значения температур заданы по Кельвину - ничего не выйдет. Так как задача банально будет уже не физична.

Ссылка на сообщение
Поделиться на других сайтах
13 минут назад, VICT0R_1945 сказал:

Предположим, значения температур заданы по Кельвину - ничего не выйдет. Так как задача банально будет уже не физична.

Значит алгоритм будет зависеть от того, от чего зависеть не должен :) Мы же говорим об уравнении Лапласа в конечном итоге, а там вообще вторые производные.

Даже самый плохой алгоритм не может тебе это позволить.

Ссылка на сообщение
Поделиться на других сайтах
47 минут назад, VICT0R_1945 сказал:

Предположим, значения температур заданы по Кельвину - ничего не выйдет.

На входе задаете температуры по Кельвину и на выходе правильные  температуры по Кельвину.

Иначе и быть не может. :biggrin:

Ссылка на сообщение
Поделиться на других сайтах
6 часов назад, VICT0R_1945 сказал:

Вообщем пока не доказано обратное (что при определенных границах или ГУ температуры ниже расчетных на ПК),

Так давайте проверим. Я так понимаю что приведенный в учебнике метод работает только с узлами, лежащими на сетке с равным шагом но сама форма может быть не только квадратной. Давайте посчитаем что-то более сложное - например сектор окружности.

 

6 часов назад, VICT0R_1945 сказал:

это был бы неплохой материал для изучения

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

Ссылка на сообщение
Поделиться на других сайтах
3 часа назад, karachun сказал:

Так давайте проверим. Я так понимаю что приведенный в учебнике метод работает только с узлами, лежащими на сетке с равным шагом но сама форма может быть не только квадратной. Давайте посчитаем что-то более сложное - например сектор окружности.

 

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

@karachun Хорошо, вот вам задача по теме. Так сказать по заявкам)

P.S. Эта задача хорошо иллюстрирует, как нужно разбивать цилиндры и почему современные мещеры стремятся в эту сторону.

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

Задача 3.png

Изменено пользователем VICT0R_1945
Ссылка на сообщение
Поделиться на других сайтах
1 час назад, VICT0R_1945 сказал:

P.S. Эта задача хорошо иллюстрирует, как нужно разбивать цилиндры и почему современные мещеры стремятся в эту сторону.

О чем это вы ???

 

1 час назад, VICT0R_1945 сказал:

Задача на порядок сложнее,

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

Ссылка на сообщение
Поделиться на других сайтах
8 минут назад, Борман сказал:

О чем это вы ???

 

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

1) Если вы зацепляли мешингами цилиндрические модели, то должны были видеть, как плохо это делают многие* мешинги без человеческой помощи.

2) Задача решается внутри круга в этом примере для простоты вычислений, так-то можно уйти от равномерного шага, просто в примере решили таким не заниматься, однако в действительности фигура все ровно будет описываться прямоугольной сеткой, т.е. нет тут КЭ.

* ICEM CFD неплохо зацепляет, но даже ему нужна порой человеческая помощь.

P.S. По своей сути в этом и есть идея МКЭ, создать хороший генератор сетки, но за 50 лет ничего так и не создали более-менее путного... Сейчас поддерживают многогранные элементы, например во Fluint'е, но ранее такого вроде как не было, прорыв этого 10-летия.

Изменено пользователем VICT0R_1945
Ссылка на сообщение
Поделиться на других сайтах
1 минуту назад, VICT0R_1945 сказал:

1) Если вы зацепляли мешингами цилиндрические модели, то должны были видеть, как плохо это делают многие* мешинги без человеческой помощи.

Почему? O-grid придумали еще очень давно и его везде используют.
Ди а не O-grid ом единым. Например вот тут процентов 80 элементов близки к кубам.

123.png

Ссылка на сообщение
Поделиться на других сайтах
1 минуту назад, karachun сказал:

Почему? O-grid придумали еще очень давно и его везде используют.
Ди а не O-grid ом единым. Например вот тут процентов 80 элементов близки к кубам.

123.png

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

Ссылка на сообщение
Поделиться на других сайтах
Только что, VICT0R_1945 сказал:

может у вас продвинутый мешинг

Не, Фемап никогда не был сильным в 3Д сетках.

 

1 минуту назад, VICT0R_1945 сказал:

просто сразу все настраиваете.

только один размер на все тело.

 

но даже если у вас он делает сетку нерегулярную а как под номером 4 то это тоже хорошо.
регулярная сетка это конечно красиво но это не цель. главное чтобы качество самих элементов было хорошим.

 

Image 1.png

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

@karachun После того как начал делать регулярную, уже стремно смотреть на нерегулярную)) Ну и считается она вроде как быстрее и создается, кстати тоже.

P.S. Имею ввиду, править ее может и долго, но сам процесс зацепления быстрый. + Можно делать в центре большие элементы (объемы), а по границам - мелкие, хотя существенной экономии в этом на практике не выходит (20-30% по объему и соответственно времени зацеплении при изменении основного размера в 2-4 раза).

H-Grid - это вообще брак, а не сетка, на такой ничего не посчитаешь. 0-Type-grid - интересная, но не для цилиндров, а больше колец. 3 и 4 вариант норм, но 3 симпотный.

Изменено пользователем VICT0R_1945
Ссылка на сообщение
Поделиться на других сайтах
26 минут назад, VICT0R_1945 сказал:

нет тут КЭ.

Именно по этому уравнению Лапласа качество сетки до лампочки.

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

Наконец-таки смог добраться до этой задачи, посчитал сам, тоже в Ансис как и Борман, использовал Plane55 в постановке построения сетки от узлов к элементам.

В общем кому код нужен - выкину.

Вот результаты.

 

 

Головоломка X.png

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

Присоединяйтесь к обсуждению

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

Гость
Ответить в тему...

×   Вставлено в виде отформатированного текста.   Вставить в виде обычного текста

  Разрешено не более 75 эмодзи.

×   Ваша ссылка была автоматически встроена.   Отобразить как ссылку

×   Ваш предыдущий контент был восстановлен.   Очистить редактор

×   Вы не можете вставить изображения напрямую. Загрузите или вставьте изображения по ссылке.

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

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




  • Сообщения

    • mactepnew
      Вот я тоже так думал, пока не столкнулся с обратным) Как еще столкнусь, пришлю скриншот. Хаос еще в том, что порядок в дереве построения весьма условный. операция, которая стоит перед групой, имеет бОльший timestamp(пардон, не знаю как это в переводе), чем эта группа, и выполняется позже. При этом все операции в группе имеют timestamp на пару сотен меньше и выполняются раньше, чем две группы, стоящие перед ней. Логика?  Если в эскизе все элементы определены, то можно поставить еще один размер (хочешь переопределить эскиз - да пожалста!), изменить его - и NX самостоятельно выберет, какой размер или ограничение отпустить и весь эскиз летит к чертям. Это ли не хаос?   Назовете хоть один CAD, которым можно пользоваться так? Любой из встреченных мной имеет много неочевидных нюансов, которые конечно не проявляются при выполнении университетского или базового курса. В Creo я на протяжении 5 лет находил всё новые фишки, которые не прописаны в документации. Огромная благодарность Ruslan за пояснения тонкостей. Даже Компас и Solidworks, которые считаются самыми дружественными, оказалось непросто обуздать, когда проходил собеседования и делал тестовые. А ведь университетский курс знал на отлично!   Возвращаясь к NX - вот сегодняшний пример. Утром деталь была оранжевая, единая. Изменил пару уклонов и толщин, стала желтая с элементами оранжевого и частично твердотельной, частично листовой. Я даже не знаю, как сформулировать вопрос в поддержку. 
    • Mikhail_Podgornyj
      Добрый день, кто может поделиться SKF Bearing V1.6 (ANSYS App)? Заранее благодарен!
    • kkk
      Короче, как сделать автоматом не нашел. Раскрасил нужные грани нужными цветами. Затем на вкладке с цветом модели можно выбирать нужный цвет, тогда в модели выбираются грани, покрашенные этим цветом. Можно еще сохранить выбор этих граней, но смысла не вижу, т.к. все легко выбирается по цвету. Затем либо в инструменте Измерить смотрим площадь и как-то ее используем, либо можно макросом (что выкладывал Kir95 ) внести площадь в свойства (на вкладку Конфигурации). Можно еще внести выбранные грани в датчики, но как использовать их в свойствах не нашел.    Макрос прикрепляю.     Площадь выбранных граней в свойства.swp
    • IgnatvKudrin
      Для нашего производства требуются постоянно инструментальные пружины для штампов. Сейчас заказы обрабатываются довольно долго - до 2-х месяцев. А у нас сроки на штампы 1-2 месяца с испытанием, доработкой и т.п. Хотелось бы найти поставщика, у которого уже есть склад с данными пружинами или кто может доставить за 1 месяц.   Пружины стандартные по ISO, но на всякий случай приложил каталог. Вот список пружин, которые нужны очень срочно:   LR 32 x 44 - 1 шт. LR 25 x 51 - 2 шт. LB 63 x 89 - 7 шт. LB 32 x 102 - 10 шт. LB 32 x 89 - 20 шт. LB 32 x 76 - 12 шт. LB 32 x 51 - 6 шт. LB 25 x 127 - 9 шт.   Готов выслушать все предложения. Возможно даже есть те, кто сможет изготовить подобные пружины с такими же параметрами. 20231101.pdf
    • bubblyk
      Знатоки форума, подскажите, как смоделировать простую (на первый взгляд) задачу термомеханики в упругой постановке.   Есть брусок правильный формы ака кирпич, в котором генерируется НЕРАВНОМЕРНОЕ удельное тепловыделение (HGEN). Задняя и передняя поверхности конвективно охлаждаются с известными (но различными) коэффициентами теплоотдачи и температурами жидкости/газа. Свойства кирпича (теплофизика + механика) принимаются зависящими от температуры. Постановка - стационарная Собственно вопрос : как грамотно смоделировать температурное поле в кирпиче с учетом температурного расширения?   Исходный код приложил.     Надеюсь на коллективный разум и опыт! EPTH2D.mac p.s.   Если есть опыт решения этой задачи в Ansys WB, тоже можно сюда приложить командные вставки.
    • OutXL
      Это было давно использовал NC-Link  и данную инструкцию. The Maho.docx
    • Mikhail_Podgornyj
      Добрый день, кто может поделиться SKF Bearing V1.6 (ANSYS App)? Заранее благодарен!
    • alexqwq10
      В данных корпусах ничего лишнего быть не должно. Максимум шпильки заземления, в количестве предусмотренном для конкретного корпуса. Рассматривал очень много подобных корпусов разных производителей и данная проблема актуальна для подавляющего большинства.
    • MFS
      Коллеги, приветствую, нужно штамповать на горячую из алюминиевого сплава (АК6 или около того) заготовки размером 280*180*12. Подскажите, пожалуйста, параметры пресса для такой заготовки - размеры, усилия. А есл дадите пример такой марки и модели, более менее подходящей для этой работы (нашего или импортного), то будет и еще лучше. Заранее спасибо!
    • Elbaro
      Добрый день. Уже несколько дней пытаюсь найти на просторах инета макрос, позволяющий менять материал выделенных деталей в сборке. У меня есть макрос это делающий, но только в открытой детали, заставить его работать в сборке не получается. Может кто поделится ссылкой или часть кода? Спасибо.
×
×
  • Создать...