2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




Начать новую тему Ответить на тему На страницу Пред.  1 ... 5, 6, 7, 8, 9
 
 Re: COMSOL, ANSYS etc.
Сообщение08.05.2019, 12:49 
Аватара пользователя


31/10/08
1244
Alex-Yu в сообщении #1391641 писал(а):
Ну это как раз просто: для каждой тройки координат (только на поверхности) определяем два узла сетки.

А разве мы тогда не получим деление на 0? Длина между узлами нулевая. Удваивать стоит количество коэффициентов.

Alex-Yu в сообщении #1391641 писал(а):
. И то, что на одной стороне, может быть довольно сложным образом связано с тем,что на другой.

Почему сложным?
https://ru.wikipedia.org/wiki/Произвольный_разрыв

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

У вас есть условие связи двух границ. В момент контакта оно накладывается на общие узлы. Для линейных систем так вообще всё просто и дописывается с права. Для нелинейных берём известные формулы, для неизвестных оно записывается пользователем. Далее линеаризуем и получаем численную схему.

Alex-Yu в сообщении #1391630 писал(а):
При этом построение базисных функций (а следовательно и их матриц перекрытия и т.п.) запутаются с граничными условиями.

Так как это уже не граница они просто выключаются из работы. Так как за них работает формула "разрыва" или "контакта".

 Профиль  
                  
 
 Re: COMSOL, ANSYS etc.
Сообщение08.05.2019, 12:57 
Заслуженный участник


21/08/10
2462
Pavia в сообщении #1391652 писал(а):
разве мы тогда не получим деление на 0? Длина между узлами нулевая.


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

-- Ср май 08, 2019 16:58:42 --

Pavia в сообщении #1391652 писал(а):
Почему сложным?
https://ru.wikipedia.org/wiki/Произвольный_разрыв


Во-первых это дико частный случай: гидродинамика. Во-вторых это разрыв не в том смысле (в смысле разрыв течения, а не граница двух сред).

-- Ср май 08, 2019 17:01:09 --

Pavia в сообщении #1391652 писал(а):
Вы говорите что на границе у вас должно быть равенство частных производных что соответствует непрерывности потока.



Может, кто-то это и говорил, но точно не я. И какой еще поток? В общем случае понятия потока может вообще не быть.

-- Ср май 08, 2019 17:02:36 --

Pavia в сообщении #1391652 писал(а):
Так как это уже не граница они просто выключаются из работы. Так как за них работает формула "разрыва" или "контакта".



amon говорил нечто противоположное. Во всяком случае я его так понял.

 Профиль  
                  
 
 Re: COMSOL, ANSYS etc.
Сообщение08.05.2019, 15:48 
Аватара пользователя


31/10/08
1244
Alex-Yu в сообщении #1391656 писал(а):
Ребра между этими узлами устраивать не надо. Все ребра только с одной стороны поверхности (в смысле что нет ребер, начинающихся с одной стороны и заканчивающихся с другой).

Раз рёбер нет. Значит нет связи. Нет формулы которая их объединяет так?
Тогда зачем они нам нужны, зачем мы их удваеваем? Что мы такого делаем с новыми узлами того, что не делаем со старыми?

Допустим 1D случай. Две лодочки едут на встречу сталкиваются. Лодочки изготовлены из метаматерьяла у одной на границе функция $f()$ у другой $g()$.

Alex-Yu в сообщении #1391656 писал(а):
И какой еще поток? В общем случае понятия потока может вообще не быть.

Если потока нет, то это не FEM метод. У вас есть энергия и у неё есть поток выраженный дифуром. Не хотите поток энергию пусть будет перенос информации. Нету потока, следовательно нет дифура. Что решать собрались?

Alex-Yu в сообщении #1391656 писал(а):
Во-первых это дико частный случай: гидродинамика.

Не вижу проблем для обобщения.

Цитата:
amon говорил нечто противоположное.

Пояснения от amon тоже не помешали бы.

 Профиль  
                  
 
 Re: COMSOL, ANSYS etc.
Сообщение08.05.2019, 15:50 
Заслуженный участник
Аватара пользователя


01/08/06
3131
Уфа
Я COMSOL не видел. Но по опыту работы с другими системами (в т.ч. с другой стороны баррикад :-) ) могу сказать, что для настолько больших и сложных систем не бывает (и не может быть) руководств, исчерпывающе отвечающих на все вопросы эксплуатации. Для платных продуктов есть поддержка со стороны разработчиков (включающая в себя исправления ошибок и доработки под заказчика), для бесплатных — форумы, сообщества энтузиастов (тут тоже иногда можно добиться изменений). Без этого такими сложными системами на практике невозможно пользоваться. Alex-Yu, я не исключаю, что некоторые вещи, о которых вы пишете, ещё не реализованы "из коробки", и их можно добиться только "грязными хаками". Потому что, как бы ни были круты разработчики и как бы ни были велики их возможности и опыт, нельзя объять необъятного, как говорится.

 Профиль  
                  
 
 Re: COMSOL, ANSYS etc.
Сообщение08.05.2019, 16:39 
Заслуженный участник


21/08/10
2462
worm2 в сообщении #1391704 писал(а):
о которых вы пишете, ещё не реализованы "из коробки",


Вроде утверждается, что в Comsol можно решать свои уравнения, какие хочешь. И с какими хочешь гранусловиями.

-- Ср май 08, 2019 20:40:33 --

Pavia в сообщении #1391702 писал(а):
У вас есть энергия и у неё есть поток выраженный дифуром.



Очень странное утверждение. Как будто дифуравнения могут выражать поток и ничего более... Вот вам уравнение равновесия тонкой пластинки из кубического материала, нагруженной плотностью силы $F$:

$$
\label{eq:zetadiffeq}
A_1\zeta_{,\alpha,\alpha,\beta,\beta}  + (A_2+A_3)\zeta_{,x,x,y,y} =
F (x,y) \, ,
$$

А вот гранусловия к этому уравнению:

$$
\left\{
\begin{array}{l}
\displaystyle
(A_2-2A_1)\zeta_{,\alpha ,\beta,\beta}n_{\alpha} -
A_2\zeta_{,\alpha,\beta,\gamma}n_{\alpha}n_{\beta}n_{\gamma}
  \\
\displaystyle
-  2(A_2+A_3) (\zeta_{,x,x,y}n_y^3 + \zeta_{,y,y,x}n_x^3 )   \\
\displaystyle
+ \frac{1}{R} % R^{-1}
\left[ A_2 (\zeta_{,\alpha,\alpha} - 2\zeta_{,\alpha,\beta}n_{\alpha}n_{\beta} )
+4(A_2+A_3) \zeta_{,x,y}n_x n_y \right]  \\
\displaystyle
 = 0 \\ \\
\displaystyle
(2A_1+A_2) \zeta_{,\alpha,\alpha} -
A_2\zeta_{,\alpha,\beta}n_{\alpha}n_{\beta}   \\
\displaystyle
+2(A_2+A_3)  \zeta_{,x,y}n_x n_y = 0
\, .
\end{array}
\right.
$$

Где здесь поток???

Обозначения типа принятых в ОТО: индекс через запятую --- производная, по повторяющимся индексам суммирование. $\zeta$ -- прогиб пластинки, $n$ -- вектор нормали к границе пластинки, $F$ -- нагрузка, остальное -- константы.

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

-- Ср май 08, 2019 21:03:09 --

Pavia в сообщении #1391702 писал(а):
Не вижу проблем для обобщения.



Да? Не видите проблем? Ну вот выше довольно простой пример (бывает намного хуже). Если сумеете, напишите статью. То, что при соответствующем качестве изложения пройдет "на ура" например в Int. J. Solids and Struct. (откуда списаны эти уравнения, с небольшими упрощениями) почти гарантирую :-)

 Профиль  
                  
 
 Re: COMSOL, ANSYS etc.
Сообщение08.05.2019, 19:14 


16/04/19
161
worm2
спасибо за объяснение
Pavia в сообщении #1391702 писал(а):
Если потока нет, то это не FEM метод. У вас есть энергия и у неё есть поток выраженный дифуром. Не хотите поток энергию пусть будет перенос информации. Нету потока, следовательно нет дифура. Что решать собрались?

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

 Профиль  
                  
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  [ Сообщений: 126 ]  На страницу Пред.  1 ... 5, 6, 7, 8, 9

Модераторы: photon, whiterussian, profrotter, Jnrty, Aer, Парджеттер, Eule_A, Супермодераторы



Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group