!  3D geometry
! Вращение буквы R
! Часть входного файла
!	Геометрия области - тело в форме буквы "R"

/PREP7
!  Входные данные в системе СИ

! Геометрические размеры тела в форме буквы "R"
HL=1.4 	 ! Высота 
WL=0.8   ! Ширина 
WWL=0.2 ! Ширина стенок буквы "R"
HAL=0.7 ! Дополнительный размер по высоте

PEL=(HAL/2)/(WL-WWL) ! Параметр эллиптичности для криволинейной части фигуры
LOCAL,11,1,WWL,HL-HAL/2,,,,,PEL ! Эллиптическая система координат 11
CSYS,0	! Переход в основную декартову систему координат

! Создание двумерной области

! Определение основных опорных точек для внешней границы
K,1,0,0
K,2,WWL,0
K,3,WWL,HL-HAL-WWL
K,4,WL-WWL,0
K,5,WL,0
K,6,WWL,HL-HAL
K,7,WL,HL-HAL/2
K,8,WWL,HL
K,9,0,HL

!   создание линий по точкам
L,1,2 $ L,2,3  $ L,3,4 $ L,4,5 $ L,5,6
CSYS,11    ! Переход в систему координат с номером 11
L,6,7 $ L,7,8
CSYS,0 ! Переход в основную декартову систему координат
L,8,9 $ L,9,1 
AL,1,2,3,4,5,6,7,8,9 ! Определение площади 1 по линиям

! Определение отверстия - площади 2
K,10,WWL,HL-HAL/2-(WL-2*WWL)*PEL
K,11,WL-WWL,HL-HAL/2
K,12,WWL,HL-HAL/2+(WL-2*WWL)*PEL
CSYS,11
L,10,11 $ L,11,12
CSYS,0
L,12,10
AL,10,11,12
ASBA,1,2 ! Вырезать из области 1 область 2

! Определение точек оси вращения, отстоящей от плоской области слева на WWL (параллельна оси y) 
K,1001,-WWL,0
K,1002,-WWL,HL
! Создать объем вращением области 3 вокруг оси с точками 1001 и 1002 на 45 градусов
VROTAT,3,,,,,,1001,1002,45  

! Параметры для триангуляции
DMESH=WWL/4
! Установка общего размера конечного элемента 
esize,dmesh

ET,1,SOLID70  ! 3-D 8-узловой элемент для температурного анализа, степени свободы: TEMP (температура)
! ET,1,SOLID45 ! ! 3-D 8-узловой элемент для структурного анализа, степени свободы: ux,uy,uz (перемещения)

ET,2,MESH200,6 ! КЭ типа 2: вспомогательный геометрический КЭ для построения сетки, без степеней свободы

TYPE,2 ! Выбрать тип КЭ для разбиения плоской области (вспомогательный КЭ)
!Выбрать исходную двумерную область
ASEL,S,LOC,Z,0
!Разбить эту область вспомогательным элементом
AMESH,ALL
TYPE,1 ! Выбрать тип КЭ для разбиения объемов (основной КЭ)
! Протянуть разбиение плоской области вдоль объема
VSWEEP,ALL
Asel,all ! Выбрать все области