/batch
/filename,input

/prep7

TOP_TH=.0015
TOP_RAD=.022
TOP_WI=.002
TOP_WO=.0248 !.02485
INDRAFT=.2*3.14159/180.0
OUTDRAFT=.4*3.14159/180.0
BOT_TH=.0015

PSI2PA=6894.75729317
!  Contact modulus
KN= 1e8 !0.4E6*PSI2PA*1000 !17583015 

!  Element size of hub
ESIZE1=0.001

!  Element size of core and film 
ESIZE2=0.005


et,1,plane42,,,1
et,2,contac48,,,1


!  Material properties of hub
mp,ex,1,0.4e6*PSI2PA
mp,nuxy,1,0.40

!  Material properties of film
mp,ex,2,0.4e6*PSI2PA
mp,nuxy,2,0.40

!  Material properties of core
mp,ex,3,0.1e6*PSI2PA
mp,nuxy,3,0.40

!  Material properties of contact surface
r,1,KN,,,5.0/3.0
mp,mu,4,0.5

!  Define keypoints for the hub
k,1,TOP_RAD,0.06
k,2,TOP_WO,0.06
k,3,TOP_WO+TAN(OUTDRAFT)*TOP_TH,0.06-TOP_TH
k,4,TOP_RAD+TOP_WI,0.06-TOP_TH
k,5,TOP_RAD,0.06-TOP_TH
k,6,TOP_WO+TAN(OUTDRAFT)*(0.06-0.002-BOT_TH),0.002+BOT_TH
k,7,TOP_WO+TAN(OUTDRAFT)*(0.06-0.002-BOT_TH)+0.0045,0.002+BOT_TH
k,8,TOP_WO+TAN(OUTDRAFT)*(0.06-0.002-BOT_TH)+0.0045,0.001+BOT_TH
k,9,TOP_WO+TAN(OUTDRAFT)*(0.06-0.002-BOT_TH),0.001+BOT_TH
k,10,TOP_RAD+TOP_WI+(0.06-TOP_TH-0.001)*TAN(INDRAFT),0.001
k,11,TOP_WO+TAN(OUTDRAFT)*(0.06-0.002-BOT_TH),0.001
k,12,0.055,0.001
k,13,0.055,0.001+BOT_TH
k,14,0.06,0.001+BOT_TH
k,15,0.06,0
k,16,0.055,0

!  Define lines for the hub
l,1,2	!1
l,2,3	!2
l,3,4	!3
l,4,5	!4
l,5,1	!5
l,3,6	!6
l,6,9	!7
l,9,11	!8
l,11,10	!9
l,10,4	!10
l,6,7	!11
l,7,8	!12
l,8,9	!13
l,8,13	!14
l,13,12	!15
l,12,11	!16
l,13,14	!17
l,14,15	!18
l,15,16	!19
l,16,12	!20

!  Define areas for the hub
al,1,2,3,4,5		!1
al,3,6,7,8,9,10		!2
al,7,11,12,13		!3
al,8,13,14,15,16	!4
al,15,17,18,19,20	!5

! now calculate the cross-sectional area
allsel,all,all
asel,s,,,1,5,1
asum 
*get,ATOT,area,,area


!  Define keypoints for the core and film
k,101,0.025,0.002+BOT_TH+0.1525
k,102,0.025,0.06
k,103,0.025,0.002+BOT_TH
k,104,0.03,0.002+BOT_TH
k,105,0.03,0.002+BOT_TH+0.1525
k,106,0.0575,0.002+BOT_TH+0.1525
k,107,0.0575,0.002+BOT_TH


!  Define lines for the core and film
l,101,102	!21
l,102,103	!22
l,103,104	!23
l,104,105	!24
l,105,101	!25
l,105,106	!26
l,106,107	!27
l,107,104	!28


!  Define areas for the core and film
al,21,22,23,24,25	!6
al,26,27,28,24		!7




!  Set area attributes for the hub
asel,s,area,,1
asel,a,area,,2
asel,a,area,,3
asel,a,area,,4
asel,a,area,,5
aatt,1,0,1


!  Set area attributes for the film
asel,s,area,,7
aatt,2,0,1


!  Set area attributes for the core 
asel,s,area,,6
aatt,3,0,1


!  Set element attributes/shape and mesh hub
allsel
eshape,3
esize,ESIZE1
amesh,1,5,1


!  Set the element atrtributes/shape and mesh core and film
esize,ESIZE2
amesh,6,7,1


!  Assign attributes for contact elements
type,2
real,1
mat,4


!  Select and name nodes on contact surface lines
lsel,s,line,,6
nsll,s,1
cm,GAP_HUB,node
lsel,s,line,,22
nsll,s,1
cm,GAP_CORE,node
! nsel,all
lsel,s,line,,11
nsll,s,1
cm,BOT_GAPH,node
lsel,s,line,,23
nsll,s,1
cm,BOT_GAPC,node
nsel,all


!  Now generate the contact elements
gcgen,GAP_HUB,GAP_CORE,,ESIZE1*2
gcgen,GAP_CORE,GAP_HUB,,ESIZE1*2
gcgen,BOT_GAPH,BOT_GAPC,,ESIZE1*2
gcgen,BOT_GAPC,BOT_GAPH,,ESIZE1*2



!  Apply boundary conditions
! lsel,s,line,,19
! sfl,all,pres,24623.35,24623.35

ksel,s,kp,,10
nslk,r
d,all,uy

! constraints on ferrule and film/core top
nsel,s,loc,y,0
d,all,uy
nsel,s,loc,y,0.1525+.002+BOT_TH
d,all,uy


allsel,all,all
fini


! solve it
/solve
neqit,25
solve
fini

/post1
set
! check the forces along the hub/core interface
allsel,all,all
lsel,s,line,,6
cm,HUB,line
esel,s,mat,,1
cmsel,s,HUB
fsum
*get,nforce,fsum,,item,fx

allsel,all,all
nsort,s,eqv
*get,eqvmax,sort,,max
fini


! do the optimization
/opt
opvar,TOP_WO,dv,.024,.0255
opvar,TOP_WI,dv,.0015,.0025
opvar,TOP_RAD,dv,.0175,.022
opvar,TOP_TH,dv,.001,.002
opvar,OUTDRAFT,dv,.2*3.14159/180,.5*3.14159/180
opvar,INDRAFT,dv,.2*3.14159/180,.5*3.14159/180,
opvar,BOT_TH,dv,.001,.0025
opvar,eqvmax,sv,,.36e8
opvar,nforce,sv,200,5000
opvar,atot,obj

optype,first
opfrst,100,50,.1
!optype,subp
!opsubp,100
oploop,prep

opprnt,full
opexe
fini

!/exit,all




/eof