! 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