哪位大神能帮我把这三个程序整合一下，可以直接被ansys读入。跪谢！！

caojingyy的家

第一个：

!定义骨料半径范围
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,<ang(kk),then
                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



急！！自己学习了很多还是做不好。。希望大神能帮我！！