/out,g_31,out fini /prep7 *get,_etyp,etyp,,num,max _etyp=_etyp+1 esel,,mat,,5 *get,elmx,elem,,num,max *get,r_31,elem,elmx,attr,real nsle *get,rdmn,node,,mnloc,x *get,_ymx,node,,mxloc,y *get,_ymn,node,,mnloc,y a_face=3.14*rdmn*(_ymx-_ymn) nsel,,ext nsel,r,cp,,1,10000 cm,cn1,node cm,cn2,node *get,_nct,node,,count a_node=a_face/_nct et,_etyp,31 type,_etyp r,r_31,a_node,1,1,5.67e-8 real,r_31 mat,5 *do,_i11,1,_nct *get,ndmn,node,,num,min nsel,all nsel,u,,,ndmn _nj=node(nx(ndmn),ny(ndmn),0) cmsel,,cn1 nsel,a,,,_nj e,_nj,ndmn cmsel,,cn2 nsel,u,,,ndmn cm,cn2,node *enddo /out