Ansys. Проблема разбивки

[sKoTs 0]

Имеется куб и распределенные в нем включения сферической формы. Несколько сфер соприкасаются с гранями куба.

Проблема:

Разбивка сфер, дело не хитрое. А вот куб, с вырезанными обьемами сфер разбить не получается:

*** ERROR *** CP = 21.672 TIME= 23:58:57

Volume 14 cannot be meshed. 98 location(s) found where non-adjacent

boundary triangles touch. Geometry configuration may not be valid or

smaller element size definition may be required.

*** ERROR *** CP = 21.672 TIME= 23:58:57

VMESH failure, perhaps due to:

(1) Poorly shaped triangular facets: Mesh areas with triangular

elements, modify mesh control or turn on SmartSizing to get good

triangular mesh, then VMESH again.

(2) Complex geometry: Subdivide volume & try again.

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

Тип КЭ - Solid95

Разбивка свободная на тетраэдры

Посоветуйте как разбить обьем куба?

Батничег для геометрии:

/PREP7

!* Строим обьем куба

BLOCK,0,400,0,400,0,400,

!* Строим сферы, перемещая центр координат на вектор

wpoff,70.604309,70.604309,185.845886

SPHERE,70.603310, ,0,360,

wpoff,259.627899,121.879761,144.386322

SPHERE,69.766777, ,0,360,

wpoff,-125.985855,137.094849,-0.653290

SPHERE,70.420090, ,0,360,

wpoff,10.558212,-259.267204,-0.026703

SPHERE,70.310715, ,0,360,

wpoff,-145.004074,226.477776,-95.623856

SPHERE,69.799591, ,0,360,

wpoff,61.954529,-121.734055,-164.557854

SPHERE,69.369537, ,0,360,

wpoff,81.043915,34.583969,142.289482

SPHERE,69.993210, ,0,360,

wpoff,-143.051315,120.612976,-141.912338

SPHERE,69.746620, ,0,360,

wpoff,166.785004,0.715118,-0.715157

SPHERE,69.031494, ,0,360,

wpoff,93.699707,-0.734619,149.124184

SPHERE,69.766670, ,0,360,

wpoff,-258.879112,-171.421005,111.629120

SPHERE,70.213211, ,0,360,

wpoff,258.382988,-88.548088,-128.544525

SPHERE,70.262787, ,0,360,

!Вырезаем их из куба с сохранением

VSBV,1,ALL,,DELETE,KEEP

[Форумный боец 1]

проблема известная.

[sKoTs 0]

проблема известная.

И имеет простое решение?

[_Fedor_ 16]

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

Конкретно -

Избавиться от касательной. Для этого или размер кубика немного увеличить - уменьшить, или шариков

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

Или через предельный переход, примерно как вводят обобщенные функции

[sKoTs 0]

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

[Форумный боец 1]

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

[_Fedor_ 16]

"случайный характер структуры ... конечно за точность"

Приличную случайность надо тщательно организовывать. У Кнута пол толстого тома этому посвящено

[a_schelyaev 367]

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

Делов то.

[_Fedor_ 16]

Что-то заклинило. Помогите.

Есть две линии в плоскости. Найти точку пересечения этих линий. Формально они не пересекаются. Не доходят друг до друга.

Но точка пересечения в принципе существует

[sKoTs 0]

Что-то заклинило. Помогите.

Есть две линии в плоскости. Найти точку пересечения этих линий. Формально они не пересекаются. Не доходят друг до друга.

Но точка пересечения в принципе существует

ну ващето функция прямой на плоскости y=kx+b ... берешь две точки линии, решаешь систему - получаешь функцию прямой, тоже самое для другой линии ... потом ищешь пересечение функций ...

[_Fedor_ 16]

Это то понятно, интересует конкретно в Ansys как делать. Например, точку на отрезке поставить не проблема, но между [0,1] параметрически. Надо бы экстраполировать, например, задать 2 или -3 но чертовы программеры поставили перехват и сбрасывают к 0 или 1. Ну выдавали бы предупреждение коль хочется, но мешать то зачем?

Обычно наши этим грешат, но и американцы оказывается too

"функция прямой на плоскости " - я вообще-то писал о линиях, а не прямых линиях

kbetw - help - бесконечный кладезь мудрости

[Islander 2]

LEXTND не подойдет (потом оверлап сделать)?

[_Fedor_ 16]

Спасибо. С Вашим предложением вроде удобнее работать