Here is my input file :
/FILNAME, bobine, 1
/PREP7
/TITLE, OPTIMISATION, Bobine Forme
ANTYPE, STATIC
ET,1,SOLID45,,1
MP,EX,1,2.1E5
MP,NUXY,1,0.3
/COM TABLE OF NONLINEARITY
TB,BKIN,1,1 ! DEFINE NON-LINEAR MATERIAL PROPERTY FOR STEEL
TBTEMP,0
TBDATA,1,86000,1.0E4
ALPHA = -5.0
RIN = 25.
EP = 10.0
LEN1 = 250.
REX = EP + RIN
BETA = 2.0/50./3.1415*360.
CSYS,1
LEN2 = 245.
CYL4,,,RIN,(90 - BETA),REX,(90 + BETA),LEN1
LESIZE,9,,,100
LESIZE,10,,,100
LESIZE,11,,,100
LESIZE,12,,,100
VMESH,1
/COM change the nodal coordinate systeme to cylindrical systeme
NROTAT,ALL
SAVE
SAVE,bobine,db,
FINISH
/SOLU
/com selectionner les points sur RIN
NSEL,S,LOC,X,RIN
D,ALL,UX
NSEL,S,LOC,Z,0.
D,ALL,ALL
NLGEOM,ON
PRED,ON
AUTOTS,ON
OUTRES,BASIC,10
OUTPR,,30
OUTPR,VENG,30
SAVE
TIME,1E-6
NSUBST,30
/COM imposed displacement
NSEL,S,LOC,Z,LEN1
D,ALL,UZ,(LEN2 - LEN1)
NSEL,R,LOC,Y,(90 + BETA)
D,ALL,UY,ALPHA
LSWRITE,1
NSEL,ALL
ESEL,ALL
SOLVE
TIME,1
NSUBST,30
/COM imposed displacement
NSEL,S,LOC,Z,LEN1
D,ALL,UZ,(LEN2 - LEN1)
NSEL,R,LOC,Y,(90 + BETA)
D,ALL,UY,ALPHA
LSWRITE,1
NSEL,ALL
ESEL,ALL
FINISH
/POST1
SET,1,LAST,1,