求助！ansys apdl命令流编写的三个小程序怎么整合到一起？

第一个：
!定义骨料半径范围
rmin=1
rmax=2
!定义该组骨料个数
num=100
!定义两点间最大角度和最小角度
da_max=90
da_min=10



cum=0
*do,i,1,100000000
*if,cum,eq,num,then
*exit
*endif
x=rand(xmin,xmax)
y=rand(ymin,ymax)
r=rand(rmin,rmax)
*if,x-1.05*r,gt,xmin,and,x+1.05*r,lt,xmax,then
*if,y-1.05*r,gt,ymin,and,y+1.05*r,lt,ymax,then
*if,cum,eq,0,and,index,eq,0,then
cum=cum+1
cir(cum,1)=x
cir(cum,2)=y
cir(cum,3)=r
*else
sum=0
*do,j,1,cum+index
dist=sqrt((cir(j,1)-x)**2+(cir(j,2)-y)**2)
*if,dist,lt,1.05*(cir(j,3)+r),then
*exit
*else
sum=sum+1
*endif
*enddo
*if,sum,eq,cum+index,then
cum=cum+1
cir(cum+index,1)=x
cir(cum+index,2)=y
cir(cum+index,3)=r
*endif
*endif
*endif
*endif
*enddo

index=index+num


第二个：

*afun,deg
/prep7
*do,i,1,index,1
n=nint(3+rand(0,1)*(ne-3))
*dim,ang,array,n
*dim,delta,array,n
*do,j,1,10000000,1
*do,k,1,n,1
ang(k)=rand(0,1)*360
*enddo
*vscfun,amax,max,ang
*vscfun,amin,min,ang
dmax=amax-amin
dmin=360
*do,jj,1,n,1
a=ang(jj)
*voper,delta,ang,sub,a
delta(jj)=360
*do,kk,1,n,1
delta(kk)=abs(delta(kk))
*enddo
*vscfun,dmin0,min,delta
*if,dmin0,lt,dmin,then
dmin=dmin0
*endif
*enddo
*if,dmax,ge,da_max,and,dmax,le,360-da_min,then
*if,dmin,ge,da_min,then
*exit
*endif
*endif
*enddo
*do,kk,1,n-1,1
*do,p,kk+1,n,1
*if,ang(p),lt, a=ang(kk)
ang(kk)=ang(p)
ang(p)=a
*endif
*enddo
*enddo
*do,kkk,1,n,1
xx=cir(i,3)*cos(ang(kkk))+cir(i,1)
yy=cir(i,3)*sin(ang(kkk))+cir(i,2)
k,,xx,yy
k%kkk%=_return
*enddo
*if,n,eq,3,then
a,k1,k2,k3
*elseif,n,eq,4,then
a,k1,k2,k3,k4
*elseif,n,eq,5,then
a,k1,k2,k3,k4,k5
*elseif,n,eq,6,then
a,k1,k2,k3,k4,k5,k6
*elseif,n,eq,7,then
a,k1,k2,k3,k4,k5,k6,k7
*elseif,n,eq,8,then
a,k1,k2,k3,k4,k5,k6,k7,k8
*endif
*del,,4
*del,,7
*enddo

第三个：

!定义坐标范围
xmin=0
xmax=100
ymin=0
ymax=100
!定义骨料总个数
total=10000
!存放骨料位置及半径的数组
*dim,cir,array,total,3
index=0
!多边形的最大边数
ne=8




希望大神能帮忙解决！！！