!  3D geometry
! Rotation of letter R
! Not an input file


/PREP7
!  Input data in SI system

! Geometrical sizes of the letter "R"

HL=1.4 	 ! height 
WL=0.8   ! width 
WWL=0.2 ! width of the walls of "R"
HAL=0.7 ! additional parameter

PEL=(HAL/2)/(WL-WWL) ! Ellipticity parameter for the curved part of the figure
LOCAL,11,1,WWL,HL-HAL/2,,,,,PEL ! Elliptical coordinate system 11
CSYS,0	! Change active coordinate system to Cartesian 

! Construct a 2D-domain

! Define main keypoints for external boundary
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


 !   Define lines by end keypoints
L,1,2 $ L,2,3  $ L,3,4 $ L,4,5 $ L,5,6
CSYS,11    ! Change active coordinate system to coordinate system 11
L,6,7 $ L,7,8
CSYS,0 !  Change active coordinate system to Cartesian 
L,8,9 $ L,9,1 
AL,1,2,3,4,5,6,7,8,9 ! Define area 1 by lines

!Define the hole - area 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 ! Substract area 2 from area 1

! Define keypoints for an axis of rotation (parallel to y-axis) 
K,1001,-WWL,0
K,1002,-WWL,HL

! Create volume by rotating area 3 around the axis defined by points 1001 and 1002 for an angle of 45 degrees
VROTAT,3,,,,,,1001,1002,45 

! parameter for element size
dmesh=WWL/4

! set finite element size
esize,dmesh


ET,1,SOLID70  ! 3-D 8 node thermal solid, degree of freedom: TEMP (temperature)
! ET,1,SOLID45 ! 3-D 8 node structural solid, degrees of freedom: ux,uy,uz (displacements)

ET,2,MESH200,6 ! FE of type 2: auxiliary mesh-only element, no degrees of freedom
TYPE,2 ! set element type 2 (to mesh 2D domain)
! Select existing 2D-domain
ASEL,S,LOC,Z,0
! Mesh 2D-domain with auxiliary elements to define mesh pattern
AMESH,ALL
TYPE,1 ! set element type 1 (to mesh a volume), main FEf
! Fill an existing unmeshed volume with elements by sweeping the mesh from an adjacent area through the volume
VSWEEP,ALL
Asel,all ! select all areas