*HEADING
Cylindrical element testing (CCL18)
*PREPRINT, ECHO=YES, MODEL=YES, HISTORY=YES
**
**********************************************************
** Parameter definition
**
*PARAMETER
INRAD=1.5
EXRAD=2.0
HEIGHT=1.0
**
MIDHEIGHT=HEIGHT/2
MID=(INRAD+EXRAD)/2
**
***********************************************************
** Model definition
**
*NODE, NSET=SECTION
1, <INRAD>,0,0
2, <EXRAD>,0,0
3, <EXRAD>,0,<HEIGHT>
4, <INRAD>,0,<HEIGHT>
5,<MID>,0,0
6,<EXRAD>,0,<MIDHEIGHT>
7,<MID>,0,<HEIGHT>
8,<INRAD>,0,<MIDHEIGHT>
9,<MID>,0,<MIDHEIGHT>
*NCOPY, CHANGE NUMBER=100, OLD SET=SECTION, NEW SET=REST, SHIFT, MULTIPLE=7
0,0,0
0,0,0,0,0,1,45
*NSET, NSET=ALL
REST,SECTION
**
*NSET, NSET=INT, GENERATE
1,701,100
4,704,100
8,708,100
*NSET, NSET=EX, GENERATE
2,702,100
3,703,100
6,706,100
*NSET, NSET=TOP, GENERATE
3,703,100
4,704,100
7,707,100
*NSET, NSET=BOTTOM, GENERATE
1,701,100
2,702,100
5,705,100
*NSET, NSET=AR1, GENERATE
1,701,100
*NSET, NSET=AR2, GENERATE
2,702,100
*NSET, NSET=AR3, GENERATE
3,703,100
*NSET, NSET=AR4, GENERATE
4,704,100
**
*ELEMENT, TYPE=CCL18H, ELSET=LAST
7, 1,2,3,601,602,603,701,702,703,5,6,9,605,606,609,705,706,709
8, 1,3,4,601,603,604,701,703,704,9,7,8,609,607,608,709,707,708
*ELEMENT, TYPE=CCL18H, ELSET=FIRST
1, 201,202,203,1,2,3,101,102,103,205,206,209,5,6,9,105,106,109
2, 201,203,204,1,3,4,101,103,104,209,207,208,9,7,8,109,107,108
*ELGEN, ELSET=FIRST3
1,3,200,2
2,3,200,2
*ELSET, ELSET=ALL
FIRST3,LAST
*SOLID SECTION, ELSET=ALL, MATERIAL=MYMAT
*MATERIAL, NAME=MYMAT
*HYPERELASTIC,N=2
1.E5,0.5E5,0.5E5,0.8E5,0.75E5,1.E-7,1.E-7
*SURFACE, NAME=PSURF
2,S4
4,S4
6,S4
8,S4
************************************************************
** Boundary Condition
*TRANSFORM, TYPE=C, NSET=ALL
0,0,0, 0,0,1
*BOUNDARY
BOTTOM,1,3,0
ALL,2,2,0
************************************************************
** History data
**
*STEP, nlgeom
*STATIC
*DSLOAD
PSURF,P,10e3
*CONTROLS,PAR=FIELD,FIELD=DISP
 1.0e-5
*OUTPUT,HIST,FREQ=0
*OUTPUT, FIELD
*NODE OUTPUT,NSET=ALL
 U
*NODE OUTPUT,NSET=BOTTOM
 RF
*ELEMENT OUTPUT
 E
*END STEP