analtemp=arg1
gimbal_al_slv,1
gimbal_slv_rough
/prep7
parm_val,'gimb_slv','Dgb01','D01',4.4265
parm_val,'gimb_slv','Dgb04','D08',2.6717
parm_val,'gimb_slv','Dgbz01','Dgbz01',0.711
parm_val,'gimb_slv','Dgbz02','Dgbz02',1.122
parm_val,'gimb_slv','Dgbz03','Dgbz03',2.815
parm_val,'roughslv','Dz07','Dz07',2.875
csys,CSslv1
ksel,s,loc,x,D08/2
ksel,r,loc,z,-Dz07
lslk,s,1
asll,s,0
asel,u,loc,z,-Dz07
cm,aset1,area
ksel,s,loc,x,D01/2
ksel,r,loc,z,0
lslk,s,1
asll,s,0
asel,u,loc,z,0
cmsel,a,aset1
cm,srcaset,area
csys,CSXdatm1
asel,s,loc,x,D08/2
asel,r,loc,z,-Dgbz03,-Dgbz02
cm,aset2,area
asel,s,loc,x,D01/2
asel,r,loc,z,-Dgbz01,0
cmsel,a,aset2
cmsel,u,srcaset
cm,trgaset,area
init_entity_num,1
cpl_nodes,'srcaset','trgaset',CSslv1,0.090/2,'ux'
finish
/solu
analtemp=70
tref,70
antype,0
eqslv,pcg,1e-6,2
mat_extract,'alpx','ala356'
gimalpx=alpx
mat_extract,'alpx','c416'
slvalpx=alpx
rdelta=0.001
rtemp=(gimalpx-slvalpx)*(analtemp-70)*(4.4263/2)
tfinal=-(rdelta+rtemp)/((4.4263/2)*(gimalpx-slvalpx))
tunif,tfinal
solve