Оптимизация плоской задачи

[Dandy 3]

Всем доброго времени суток !

Решаю двумерную задачу оптимизации и столкнулся с проблемой. При просмотре множества примеров, в качестве State variable используют объем (т.е. при решении 2D задач используют опцию элемента Plane stress with thickness с заданием толщины), а в качестве минимизируемой функции Objective напряжения. Мне также необходимо минимизировать напряжения, только в моем случае я рассматриваю чисто плоскую задачу и в качестве State variable я хочу взять площадь поверхности, но нигде не могу найти как ее посчитать и забить в параметры :(

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

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

Заранее всем спасибо.

Изменено 3 февраля 2010 пользователем Dandy

[Dandy 3]

Небольшая поправка:

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

Подскажите, как быть в этом случае случае?

Изменено 3 февраля 2010 пользователем Dandy

[Romich_Lukin 0]

Я так понял сила прикладывается не нормально к поверхности?

Если так то смотри команду SELTOL, <размер погрешности выборки> (это только для команд селект).

И задай ее чуть больше размера элемента, сетку при этом можно сделать поменьше.

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

По поводу площади

Main Menu >> parametrs > Get scalar data > Model data + Areas > Имя параметра, номер ареаса, и выбрать данные "площадь" (Area)

или макросом

*GET,[имя параметра],AREA,[номер поверхности],AREA

Небольшая зарисовка на тему приложения силы

seltol,1.5*[размер элемента]
NSEL,S,LOC,X,0  
NSEL,R,LOC,Y,0  
*GET,nodecount,NODE,,COUNT, , , ,   ! Получить кол-во элементов
F,ALL,[направление силы],[величина силы]/nodecount
ALLSEL,ALL  

[Dandy 3]

Romich_Lukin - вы мой спаситель !!! Не знал про команду SELTOL. Все вроде бы работает.

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