Jump to content

Построение тела, ограниченного поверхностями, средствами Ansys и APDL


Дмитрий Ермоченко

Recommended Posts

Дмитрий Ермоченко

Здравствуйте! Вопрос жизни и смерти - необходимо построить тело, ограниченное шестью поверхностями:

 

1) x^2+z^2=6*y
2) x^2+z^2=1,2
3) y=0,5
4) y=4,5
5) x=0
6) z=0
 
 Полученное тело необходимо разбить сеткой КЭ (преимущественно упорядоченной) и вычислить объем полученного тела и площади ограничивающих его поверхностей для геометрической и КЭ модели. Так же необходимо представить алгоритм решения задачи в виде кода.
Если с поверхностями-плоскостями все понятно (можно просто создать точку с нужными координатами и отрезать рабочей плоскостью), то с первыми двумя поверхностями я справиться не могу.
Заранее спасибо!
Link to post
Share on other sites


UnPinned posts
piden
5 hours ago, Дмитрий Ермоченко said:

Вопрос жизни и смерти

Наверно, проще всего будет построить линии пересечения этих функций с плоскостью XY, создать Area, замешить ее квадами, потом прокрутить на pi/2 вокруг Y.

image.png

Link to post
Share on other sites
Борман
14 минут назад, piden сказал:

Наверно, проще всего будет построить линии пересечения этих функций с плоскостью XY, создать Area, замешить ее квадами, потом прокрутить на pi/2 вокруг Y.

Вторую поверхность неправильно построил.

  • Нравится 1
Link to post
Share on other sites
piden
19 minutes ago, Борман said:

Вторую поверхность неправильно построил.

Да, там цилиндр.

Link to post
Share on other sites
Борман

Попробовал построить первую (параболу) по двум точкам и наклонам в них. Всегда думал что такой сплайн должен получится параболой... но увы. Немного не лег на параболу. Остается только коническое сечение делать, если строгая парабола нужна :)

Link to post
Share on other sites
karachun

@Борман А в APDL нельзя сразу построить узлы по уравнению параболы а потом по ним уже делать плоскую сетку и проворачивать ее?

Link to post
Share on other sites
Ветерок

А использовать для этого CAD - религия не позволяет?

30 минут назад, Борман сказал:

Попробовал построить первую (параболу) по двум точкам и наклонам в них.

Для задания параболы нужно ПЯТЬ условий.

Link to post
Share on other sites
karachun
5 минут назад, Ветерок сказал:

Для задания параболы нужно ПЯТЬ условий.

У меня в Femap почему-то их только три.

Image 5.png

Link to post
Share on other sites
Дмитрий Ермоченко

Все таки мне удалось разобраться с геометрией.

1 час назад, Graf Kim сказал:

Типа такого?

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

i9bLGRg1kNc.jpg?size=907x855&quality=96&

 

Далее, чтобы получить данные по модели, воспользовался командой GET и снял данные по объему:

VSUM,FINE
*GET,PAR10,VOLU,3,VOLU,

И после по поверхностям:

ASUM,FINE
*GET,PAR90,AREA,0,AREA

 

Сетку делал не упорядоченной, но достаточно мелкой:

MnTKRgwaGQY.jpg?size=997x761&quality=96&

 

По итогу дело осталось за малым:

Мне теперь через команду GET необходимо снять суммарную объем и площадь поверхности элементов. Если я правильно понимаю, то мне необходимо снять номер общего количества элементов, номер минимального и максимального элемента, а затем просто через цикл это как то просуммировать. Я разобрался как снять номера элементов и их параметры, но не до конца понимаю как их просуммировать. Т.е:

*GET,PAR0,ELEM,0,COUNT,,,
*GET,PAR1,ELEM,0,NUM,MIN,,
*GET,PAR2,ELEM,0,NUM,MAX,,

 

А дальше цикл от PAR1 до PAR2 с шагом один. Только вот как присвоить через i уникальное имя каждому снимаемому параметру с элемента(а у меня их почти 48 тысяч) и потом их просуммировать я не до конца понимаю. Спрашиваю, наверное, банальщину, но на экзамене все приходится схватывать на лету.

Код построений прилагаю (в том числе и с целью помочь какому-нибудь несчастному студенту, которому это задание попадется на экзамене как и мне):

/PREP7
!Кривая 1
!Параметры
x0=0
z0=0
!Начало цикла
*do,i,0,25
z=z0+0.15*i
x=x0+0.15*i
!Уравнение поверхности
y=(x**2+z**2)/6
!Создание точек
K,,x,y,Z
!Завершаем цикл
*enddo
!Блок команд, отвечающих за построение модели
FLST,3,i+1,3
!Цикл создания линий
*do,j,1,i+1
!Создаем массив
FITEM,3,j
*enddo
!Завершаем цикл
!Создаем сплайн по точкам
BSPLIN,,P51X
!Создаем ось вращения
k,,0,0.5,0
!Создаем поверхность вращения
FLST,2,1,4,ORDE,1   
FITEM,2,1   
FLST,8,2,3  
FITEM,8,27  
FITEM,8,1   
AROTAT,P51X, , , , , ,P51X, ,360, ,

!Поверхность y=4.5
!Создаем точку с необходимыми координатами
k,,0,4.5,0
!Ориентируем рабочую плоскость
KWPLAN,-1,      28,      29,      30
!Переносим раб. плоскость в созданную точку
KWPAVE,      31
!Обрезаем созданные грани по верхней границе
FLST,2,4,5,ORDE,2   
FITEM,2,1   
FITEM,2,-4  
ASBW,P51X
!Удаляем лишние поверхности, точки и линии
!Удаление лишних поверхностей
FLST,2,4,5,ORDE,2   
FITEM,2,9   
FITEM,2,-12
ADELE,P51X
!Удаление оставшихся точек и линий
FLST,2,8,4,ORDE,4   
FITEM,2,5   
FITEM,2,-8  
FITEM,2,17  
FITEM,2,-20
LDELE,P51X, , ,1
!Создаем поверхность по линиям
FLST,2,4,4  
FITEM,2,11  
FITEM,2,13  
FITEM,2,15  
FITEM,2,16  
AL,P51X

!Создаем объем по готовым поверхностям
FLST,2,5,5,ORDE,3   
FITEM,2,1   
FITEM,2,5   
FITEM,2,-8  
VA,P51X

!Поверхность x**2+z**2=1.2 является цилиндром с радиуосом
r=sqrt(1.2)
!Тогда
!Возвращаем рабочую плоскость в изначальную системы координат
WPCSYS,-1
WPAVE,0,0,0
!Создаем круг с заданным радиусом
PCIRC,1.095, ,0,360,
!Создаем направляющую линию для вытягивания
k,,0,0,5
k,,0,0,-5
LSTR,1,36
LSTR,1,37
!Вытягиваем поверхность круга вдоль направляющей
VDRAG,       2, , , , , ,       5   
VDRAG,       2, , , , , ,       6   
!Вычитаем объем цилиндра из основного тела
FLST,3,2,6,ORDE,2   
FITEM,3,2   
FITEM,3,-3  
VSBV,       1,P51X
!Удаляем все лишнии линии и точки направляющей
FLST,2,2,4,ORDE,2   
FITEM,2,5   
FITEM,2,-6  
LDELE,P51X, , ,

!Поверхность y=0.5
!Ориентируем раб. плоскость
KWPLAN,-1,      33,      34,      35
!Переносим раб. плоскость в нужную точку
KWPAVE,      27
!Обрезаем объем тела раб. плоскостью
VSBW,       4  
!Удаляем лишний объем и объекты
!Два лишних объема
FLST,2,2,6,ORDE,2   
FITEM,2,1   
FITEM,2,-2  
VDELE,P51X, , ,1
!Оставшиеся от направляющей точки
FLST,2,2,3,ORDE,2   
FITEM,2,36  
FITEM,2,-37
KDELE,P51X

!Плоскости x=0 & z=0
!Возвращаем рабочую плоскость в изначальную системы координат
WPCSYS,-1
WPAVE,0,0,0
!Делим объем тела пополам
VSBW,       3
!Ориентируем раб. плоскость
KWPLAN,-1,      51,      47,      31
!Делим изначальный объем на 4 части
VSBW,       1  
!Удаляем лишние объемы
FLST,2,2,6,ORDE,2   
FITEM,2,2   
FITEM,2,4   
VDELE,P51X, , ,1
!Возвращаем рабочую плоскость в изначальную системы координат
WPCSYS,-1
WPAVE,0,0,0

!Объем тела
VSUM,FINE
*GET,PAR01,VOLU,3,VOLU,
!TOTAL VOLUME =   46.277

!Площади поверхностей
ASUM,FINE
*GET,PAR02,AREA,0,AREA
!S=82.62044852

!Задаем тип конечных элементов
ET,1,PLANE183   
ET,2,SOLID186     
!Разбиваем свободной сеткой
ESIZE,0,20,
MSHAPE,1,3D
MSHKEY,0
!*  
CM,_Y,VOLU  
VSEL, , , ,       3
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  

!Анализ геометрии с сеткой КЭ
*GET,PAR03,ELEM,0,COUNT,,,
*GET,PAR04,ELEM,0,NUM,MIN,,
*GET,PAR05,ELEM,0,NUM,MAX,,

*do,i,PAR04,PAR05
???????
*enddo

 

 

 

 

 

 

 

Link to post
Share on other sites
Graf Kim
5 часов назад, Дмитрий Ермоченко сказал:

Почти, только вы немного неправильно ориентировали цилиндр

Вы уверены? Для обоих поверхностей ось вращения — OY, если смотреть на уравнения. 

Link to post
Share on other sites

Join the conversation

You can post now and register later. If you have an account, sign in now to post with your account.
Note: Your post will require moderator approval before it will be visible.

Guest
Reply to this topic...

×   Pasted as rich text.   Paste as plain text instead

  Only 75 emoji are allowed.

×   Your link has been automatically embedded.   Display as a link instead

×   Your previous content has been restored.   Clear editor

×   You cannot paste images directly. Upload or insert images from URL.

  • Recently Browsing   0 members

    No registered users viewing this page.



  • Сообщения

    • Vdmitriev
      Да не, в некотором роде пыталс обмануть математику вот этот парень. Конструктор от бога, явно инженер в душе. https://www.yaplakal.com/forum2/topic1252620.html и два https://www.yaplakal.com/forum46/topic655903.html Конструктор и экономный парень, но...   Сейчас общаюсь с таким же (по возрасту кстати такой же) - руки и голова есть, внешность хорошая, а с женщинами, как он говорит, не сложилось.
    • stanki74
      Компания ООО "АРЕАЛ" Станочное оборудование" г.Челябинск, предлагает рассмотреть перечень металлообрабатывающего оборудования с коротким сроком поставки. Стоимость оборудования включает: Таможенную очистку, НДС   Доставку до склада Покупателя   ШМР и ПНР на холостом ходу   Инструктаж по работе с оборудованием   Гарантийные обязательства 12 месяцев   Выгодные предложения от лизинговых компаний     Вертикальный обрабатывающий центр Extron M1100 Система ЧПУ Fanuc 0iMF Размер стола 1260х600 Направляющие скольжения Максимальная допустимая нагрузка на стол 1000кг Мощность привода шпинделя, постоянная 7,5/11 кВт Ременной привод 8 000 об/мин Конус шпинделя BT-40/50 Количество инструмента 20 (24*/30*) Вес  6500 кг Стоимость оборудования:  9 900 000 руб. с НДС     Вертикальный обрабатывающий центр SINO VMC850P Система ЧПУ Fanuc 0iMF Размер стола 1000х500 Направляющие качения Максимальная допустимая нагрузка на стол 750кг Мощность привода шпинделя, постоянная 11/15 кВт Прямой привод 8 000 об/мин Конус шпинделя BT-40 Количество инструмента 24 Вес  5500 кг Стоимость оборудования: 6 500 000 руб. с НДС     Вертикальный обрабатывающий центр SINO SVD-500 Система ЧПУ Fanuc 0iMF Размер стола 700х400 Направляющие линейные Максимальная допустимая нагрузка на стол 300кг Мощность привода шпинделя, постоянная 5,5/7,5 кВт Частота вращения 12 000 об/мин Конус шпинделя BT-40 Количество инструмента 16 Вес  3900 кг Стоимость оборудования: 4 820 000 руб. с НДС     Токарный обрабатывающий центр FOCUS FBL-300 Система ЧПУ Fanuc 0i-TF Диаметр обработки над станиной 510 мм Длина обработки 594 мм Диаметр прутка 77 мм  Торец шпинделя А2-8 Мощность  15 / 18,5 кВт Частота вращения  2 500 об/мин Револьверная головка 10 инструментов Направляющие скольжения  Маслоотделитель  Вес 5 800 кг Стоимость оборудования: 8 490 000 руб. с НДС     Токарный станок  BHC-2680 Система ЧПУ SIEMENS 828D Basic Диаметр обработки над станиной 800 мм Диаметр обработки над суппортом 430 мм Расстояние между центрами 2000 мм Отверстие шпинделя 107 мм  Торец шпинделя А2-11 Мощность 11 кВт Частота вращения  шпинделя 40-2350 об/мин Вес  5800 кг Стоимость оборудования:  6 900 000 руб. с НДС     Токарный обрабатывающий центр FOCUS FBL-230MC Система ЧПУ Fanuc 0i-TF Диаметр обработки над станиной 550 мм Длина обработки 322 мм Диаметр прутка 52 мм  Торец шпинделя А2-6 Мощность  11 / 15 кВт Частота вращения  3 500 об/мин 12-ти поз. приводная револьверная головка BMT55 Направляющие скольжения  Маслоотделитель  Вес  3 500 кг Стоимость оборудования: 10 000 000 руб. с НДС     Токарный станок 16Р25 (УЦИ – опция) Диаметр обработки над станиной 520 мм Диаметр обработки над суппортом 340 мм Расстояние между центрами 1500 мм Мощность  9 кВт (опция 11кВт) Частота вращения  10-1800 об/мин Вес  3590 кг Стоимость оборудования:  2 500 000,00 руб. с НДС     Вертикальный консольно-фрезерный станок 6Т13 Размер стола 1700х400 Максимальная допустимая нагрузка на стол 800кг Мощность привода шпинделя, постоянная 11 кВт Диапазон ступеней вращения шпинделя 30-1500 об/мин Устройство цифровой индикации (3-х координатное) Вес  4300 кг Стоимость оборудования:  2 900 000 руб. с НДС     Токарный станок  CMAPT Оперативная система управления, сенсорная HMI Диаметр обработки над станиной 500 мм Диаметр обработки над суппортом 300 мм Расстояние между центрами 1000 мм Класс точности установленных ШВП  P3 Трехступенчатый механический редуктор Отверстие шпинделя 107 мм  Торец шпинделя A2-8/11 Мощность 7,5 кВт Частота вращения  шпинделя 20-2000 об/мин Станина изготовлена из высококачественного чугуна Вес  2300 кг Стоимость оборудования:  2 800 000 руб. с НДС    
    • Golem
      g54; t0; t0808; g0 x48; z10; m21; g28 c0; g28 h0; g97 s2000 m3; G94 M8 G12.1 G01 z-3 f5000 G01 x41.9 c25 f2000 G01 x41.9 c-25 f150 G01 z10 f 2000 G13.1 G0 x48 M21 G0 c180 G97 s2000 m3 G94 M8 G12.1 G01 z-3 f5000 G01 x41.9 c25 f2000 G01 x41.9 c-25 f150 G01 z10 f 2000 G13.1 G0 x48 M21 G0 c90 G97 s2000 m3 G94 M8 G12.1 G01 z-3 f5000 G01 x35.9 c25 f2000 G01 x35.9 c-25 f150 G01 z10 f 2000 G13.1 G0 x48 M21 G0 c270 G97 s2000 m3 G94 M8 G12.1 G01 z-3 f5000 G01 x35.9 c25 f2000 G01 x35.9 c-25 f150 G01 z10 f 2000 G13.1 M22 G0 x100 z100 M30   Все равно пришлось с телефона писать) По поводу энкодера, но ведь в режиме оси C он позиционируется нормально.  
    • SergeyS
    • kkk
      Судя по всему кому-то ЧСВ не дает правильно оценивать происходящее, отсюда "много букафф" про личные обиды и засорение темы. Ну да ладно, уже выше писал, что мои сообщения лишь опосредованно к Вам относятся, кому будет полезно, возьмет на вооружение.   Кстати про полезности - скрины то мои смотрели, пробовали так сделать, какие результаты?
    • Docuk1
      Скорректировать кусок в коде спекэдитор.свп,  в этом окне   чтоб было так        
    • Anat2015
      Можно ссылку вставить на какой-нибудь ресурс, в который вы закачаете свою инфу.
    • Claus
      Если о простом и банальном то: например делаю сначала прямоугольник из осей, потом расставляю по углам отверстия, это как я лично считаю быстрее чем накидать отверстия и расставить по взаимосвязям. И и.п и т.д. В данном случае наверно понятнее обсуждать на примерах, вероятно нет общих механизмов, хотя есть конечно наука ОПУП и в частности логистика, можно решить это сетевыми графиками но это не наш случай.
    • Андрей2165
      Поправил, но все равно есть зарез. Где ошибка? СБОРКА2.prz
    • REDGARS
      Спасибо про напоминание о горячих клавишах. Но вопрос немного в другом. Интересуют приемы по сокращению количества определений, вспомогательных поверхностей и т.п.
×
×
  • Create New...