/prep7 ET,1,45 csys,1 pi=4*atan(1) rad1=20 rad2=25 delta=rad2-rad1 dicke=1 amplitude=2 anz=10 teiler=360 elsize=360/teiler *do,i,1,teiler,1 ys=i zs=amplitude*sin(anz*pi*i/180) k,i,rad1,ys,zs !construction keypoints k,i+teiler,rad2,ys,zs !construction keypoints *enddo *do,i,1,teiler,1 *if,i,eq,teiler,then CSKP,100+i,0,i,1,i+teiler *else CSKP,100+i,0,i,i+1,i+teiler *endif CSYS,100+i k,2*teiler+i,0,-delta/2,-dicke/2 ! 5-----------4 k,3*teiler+i,0, delta/2,-dicke/2 ! | | k,4*teiler+i,0, delta/2, dicke/2 ! | {X(i)} | k,5*teiler+i,0,-delta/2, dicke/2 ! | | *enddo ! 2-----------3 ! /NERR,0,,,,0 *do,i,1,teiler,1 *if,i,eq,teiler,then V,2*teiler+i,3*teiler+i,4*teiler+i,5*teiler+i,2*teiler+1,3*teiler+1, 4*teiler+1,5*teiler+1 *else V,2*teiler+i,3*teiler+i,4*teiler+i,5*teiler+i,2*teiler+1+i,3*teiler+1+i,4*teiler+1+i,5*teiler+1+i *endif *enddo CSDELE,all KDELE,1,2*teiler,1 !delete construction keypoints !esize,elsize !MSHAPE,0,3D !MSHKEY,1 !vmesh,all