gimbal_al_slv_fix
assemble=1
gimbal_slv_rough
csys,CSgim0
nsel,s,loc,x,0
d,all,ux
nsel,s,node,,nm-67417
d,all,all
init_entity_num,1
cpl_nodes,'srcaset','trgaset',CSslv1,0.100,'ux'
nsel,s,node,,nm-18878
d,all,uz
init_entity_num,1
finish
/solu
analtemp=-60
rdelta=0.002
t_ref=70
parm_val,'gimb_slv','Dgb01','d'
tref,t_ref
antype,0
eqslv,pcg,1e-6,2
mat_extract,'alpx','ala356'
gimalpx=alpx
mat_extract,'alpx','c416'
slvalpx=alpx
rtemp=(gimalpx-slvalpx)*(analtemp-t_ref)*(d/2)
rdelta=-rdelta
tfinal=(rdelta+rtemp)/((d/2)*(gimalpx-slvalpx))+t_ref
tunif,tfinal
*go,:skip
solve
:skip