【精华——原创】一个复杂零件激光快速堆积成形的APDL优化设计的例子

sunclark

本APDL程序对复杂形状零件的快速（焊接成形）成形进行了路径优化，并运用生死单元技术，高级编程技术如数组及矩阵在ANSYS中实现了复杂建模，并成功进行了温度场的优化分析。是我发表在《焊接学报》上的一篇文章的程序，奉献给大家，共同进步！！并欢迎讨论！本程序很长，已成功运行，大家放心试验。
FINISH
/CLEAR
/title,cap deposition manufacturing
/PREP7
/units,si
*AFUN,DEG
PRECISION,0

WPOF,,-0.006,
WPOF,-0.008,,
BLOCK,0,0.016,0,0.012,0,0.009,
WPOF,,0.006,
WPOF,0.008,,
CYL4,0,0,0.003, , , ,0.009
CYL4,0,0,0.016, , , ,0.009
WPOF,,,0.0085
/TRIAD,OFF
WPRO,,-8.8,
WPOF,,0.017,
WPRO,,90.00,
WPOF,,0.0911,
CYL4,0,0,0.0949, ,0.094, ,0.034
VINV,3,4
NUMCMP,VOLU
VOVLAP,1,2
VOVLAP,4,3
/PNUM,VOLU,1
VDELE, 2, , ,1
VDELE, 8, , ,1
VDELE, 1, , ,1
NUMCMP,VOLU
NUMCMP,KP   
NUMCMP,AREA
NUMCMP,LINE
/TRIAD,ORIG

CSYS,0  
WPAVE,0,0,0
wpro,,-90.000000,
WPRO,,8.8
CYL4,0,0,0.010, ,0.011, ,0.010
VOVLAP,1,4  
VDELE,      10, , ,1
VDELE,      11, , ,1

VSEL,S, , ,9
CYL4,0,0,0.011, ,0.012, ,0.010
VOVLAP,1,9
VDELE,      11, , ,1
VDELE,      13, , ,1

VSEL,S, , ,4
CYL4,0,0,0.012, ,0.013, ,0.010
VOVLAP,1,4
VDELE,      13, , ,1
VDELE,      14, , ,1

VSEL,S, , ,9
CYL4,0,0,0.013, ,0.014, ,0.010
VOVLAP,1,9
VDELE,      14, , ,1
VDELE,      15, , ,1

VSEL,S, , ,4
CYL4,0,0,0.014, ,0.015, ,0.010
VOVLAP,1,4
VDELE,      15, , ,1
VDELE,      16, , ,1
ALLSEL,ALL
VADD,2,3
WPOF,,,0.001
VSBW,1
WPOF,,,0.001
VSBW,3
WPOF,,,0.001
VSBW,4
WPOF,,,-0.003
save
CSYS,1  
WPRO,,,90
VSEL,S, , ,       9
VSEL,A, , ,       10
VSEL,A, , ,       11
VSEL,A, , ,       12
VSEL,A, , ,       13
VSEL,A, , ,       14
CM,V0,VOLU
VSBW,  V0

WPRO,,60,
VSEL,S, , ,       4
VSEL,A, , ,       16
VSEL,A, , ,       17
VSEL,A, , ,       18
VSEL,A, , ,       19
VSEL,A, , ,       20
VSEL,A, , ,       21
VSEL,A, , ,       22
VSEL,A, , ,       23
VSEL,A, , ,       24
VSEL,A, , ,       25
VSEL,A, , ,       26
CM,V1,VOLU
VSBW,  V1

WPRO,,60,
VSEL,S, , ,       9
VSEL,A, , ,       10
VSEL,A, , ,       11
VSEL,A, , ,       12
VSEL,A, , ,       13
VSEL,A, , ,       14
VSEL,A, , ,       28
VSEL,A, , ,       29
VSEL,A, , ,       30
VSEL,A, , ,       31
VSEL,A, , ,       32
VSEL,A, , ,       33
VSEL,A, , ,       34
VSEL,A, , ,       35
VSEL,A, , ,       36
VSEL,A, , ,       37
VSEL,A, , ,       38
VSEL,A, , ,       39
VSEL,A, , ,       40
VSEL,A, , ,       41
VSEL,A, , ,       42
VSEL,A, , ,       43
VSEL,A, , ,       44
VSEL,A, , ,       27
CM,V2,VOLU
VSBW,  V2
WPRO,,-120,
WPRO,,,-90
CSYS,0
NUMCMP,VOLU
NUMCMP,KP   
NUMCMP,AREA
NUMCMP,LINE
VPLOT
/REPLOT
SAVE

ET,1,SOLID70
MPTEMP,,,,,,,,,,,,,  
MPTEMP,1,30
MPTEMP,2,200
MPTEMP,3,300
MPTEMP,4,400
MPTEMP,5,500
MPTEMP,6,600
MPTEMP,7,700
MPTEMP,8,800
MPTEMP,9,900
MPTEMP,10,1000
MPTEMP,11,1300
MPTEMP,12,1400
MPDATA,c,1,,461   
MPDATA,c,1,,476   
MPDATA,c,1,,490   
MPDATA,c,1,,506   
MPDATA,c,1,,528  
MPDATA,c,1,,555   
MPDATA,c,1,,595   
MPDATA,c,1,,654   
MPDATA,c,1,,699  
MPDATA,c,1,,760   
MPDATA,c,1,,1400   
MPDATA,c,1,,900   
MPTEMP,,,,,,,,,,,,,  
MPTEMP,1,20
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,600
MPTEMP,8,700
MPTEMP,9,800
MPTEMP,10,900
MPTEMP,11,1000
MPTEMP,12,1100
MPDATA,KXX,1,,11.72   
MPDATA,KXX,1,,12.56   
MPDATA,KXX,1,,14.65   
MPDATA,KXX,1,,16.74   
MPDATA,KXX,1,,19.26  
MPDATA,KXX,1,,21.35   
MPDATA,KXX,1,,23.44   
MPDATA,KXX,1,,25.53   
MPDATA,KXX,1,,27.63  
MPDATA,KXX,1,,30.14   
MPDATA,KXX,1,,32.43  
MPDATA,KXX,1,,35.68  
MPTEMP,,,
MPTEMP,1,20
MPTEMP,2,1100
MPDATA,DENS,1,,8360
MPDATA,DENS,1,,8310
MPTEMP,,,
MPTEMP,1,0
MPDATA,HF,1,,33   
WPSTYLE,,,,,,,,0
ALLSEL,ALL
VSEL,S, , ,       1
VSEL,A, , ,       2
VSEL,A, , ,       3
CM,V3,VOLU
ESIZE,0.0006,0,
VSWEEP,V3
ALLSEL,ALL
VPLOT  
VSEL,U, , ,       1
VSEL,U, , ,       2
VSEL,U, , ,       3
ESIZE,0.000603,0,
VSWEEP,ALL
NUMMRG,NODE, , , ,LOW   
NUMMRG,KP, , , ,LOW
NUMMRG,KP, , , ,LOW
NUMCMP,NODE
NUMCMP,KP   
NUMCMP,LINE
NUMCMP,AREA
NUMCMP,VOLU
/PNUM,LINE,1
/PNUM,AREA,1
/PNUM,VOLU,1
save
!The element number is 9428
!The nodes number is 12773

FINISH
/SOLU
ANTYPE,4
TUNIF,20
!IC,ALL,ALL,112,
NSEL,ALL,EXT
SF,ALL,CONV,33,20
r=0.006
len1=0.007
len2=0.012
len3=0.016
ang1=0.5
v=0.006
kcoe=0.68
n=len3/(r/7.5)
X1=0
Y1=0

ALLSEL,ALL
VPLOT
*DIM,VNAME,CHAR,44
VNAME(1)='a1','b1','c1','d1','e1','f1','g1','h1','i1','j1','k1','l1','m1','n1'
VNAME(15)='a2','b2','c2','d2','e2','f2','g2','h2','i2','j2','k2','l2','m2','n2'
VNAME(29)='a3','b3','c3','d3','e3','f3','g3','h3','i3','j3','k3','l3','m3','n3'
VNAME(43)='a4','b4'
*DO,I,1,44
VSEL,S,,,I
VPLOT
ESLV,R
NSLE,R
CM,VNAME(I),NODE
CMSEL,,VNAME(I)
ESLN,S,1,ALL
EKILL,ALL
ESEL,S,LIVE
NSLE,S
*ENDDO

ESEL,S,LIVE
EPLOT
/VIEW,1,1,1,1   
/ANG,1  
/REP,FAST   


ESEL,ALL
CMSEL,S,VNAME(2)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
EPLOT
!The first layer
ALLSEL,ALL
    *DO,b,1,3
         *DO,c,1,10
         step1=2*(len3/n)
         x1=x1+step1
         LOCAL,11,0,0.008-x1,-0.006+(b-1)*0.004,0,,,,,,
         NSEL,S,LOC,X,0,0.0016
         NSEL,R,LOC,Y,0,0.002
         NSEL,R,LOC,Z,0,0.001
         NPLOT
         BF,ALL,HGEN,kcoe*6.479899014E10
         ESEL,ALL
         NSLE,ALL
         ALLSEL,ALL

         TIME,(b-1)*len3/v+c*step1/v
         ESTIF,1.0E-50
         NLGEOM,ON
         NROPT,FULL,,ON
         AUTOTS,1
         KBC,1
         lnsrch,1
         TSRES,ERASE
         OUTPR,NSOL,LAST,
         TINTP,,,,1,,,,,
         SOLVE

         ESEL,S,LIVE
         NSLE,S         
         NPLOT
         BF,ALL,HGEN,0
         ESEL,S,LIVE
         NSLE,S
         SAVE
         *ENDDO
    x1=0
    *ENDDO
save
x1=0
ALLSEL,ALL
    *DO,b,1,3
         *DO,c,1,10
         step1=2*(len3/n)
         x1=x1+step1
         LOCAL,11,0,-0.008+x1,0.004-(b-1)*0.004,0,,,,,,
         NSEL,S,LOC,X,-0.0016,0
         NSEL,R,LOC,Y,0,0.002
         NSEL,R,LOC,Z,0,0.001
         NPLOT
         BF,ALL,HGEN,kcoe*0.50*6.479899014E10
         ESEL,ALL
         NSLE,ALL
         ALLSEL,ALL

         TIME,(b-1)*len3/v+c*step1/v+8
         ESTIF,1.0E-50
         NLGEOM,ON
         NROPT,FULL,,ON
         AUTOTS,1
         KBC,1
         lnsrch,1
         TSRES,ERASE
         OUTPR,NSOL,LAST,
         TINTP,,,,1,,,,,
         SOLVE

         ESEL,S,LIVE
         NSLE,S         
         NPLOT
         BF,ALL,HGEN,0
         ESEL,S,LIVE
         NSLE,S
         SAVE
         *ENDDO
    x1=0
    *ENDDO
save
x1=0
ESEL,ALL
CMSEL,S,VNAME(1)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
EPLOT
!The second layer
ALLSEL,ALL
    *DO,b,1,4
         *DO,c,1,10
         step1=2*(len2/n)
         Y1=Y1+step1
         LOCAL,11,0,0.008-(b-1)*0.004,-0.006+Y1,0.001,,,,,,
         NSEL,S,LOC,X,-0.002,0
         NSEL,R,LOC,Y,-0.0012,0
         NSEL,R,LOC,Z,0,0.001
         NPLOT
         BF,ALL,HGEN,kcoe*0.86*4.859924261E10
         ESEL,ALL
         NSLE,ALL
         ALLSEL,ALL

         TIME,(b-1)*len2/v+c*step1/v+16
         ESTIF,1.0E-50
         NLGEOM,ON
         NROPT,FULL,,ON
         AUTOTS,1
         KBC,1
         lnsrch,1
         TSRES,ERASE
         OUTPR,NSOL,LAST,
         TINTP,,,,1,,,,,
         SOLVE

         ESEL,S,LIVE
         NSLE,S         
         NPLOT
         BF,ALL,HGEN,0
         ESEL,S,LIVE
         NSLE,S
         SAVE
         *ENDDO
    Y1=0
    *ENDDO
save
Y1=0
ALLSEL,ALL
    *DO,b,1,4
         *DO,c,1,10
         step1=2*(len2/n)
         Y1=Y1+step1
         LOCAL,11,0,-0.008+(b-1)*0.004,0.006-Y1,0.001,,,,,,
         NSEL,S,LOC,X,0,0.002
         NSEL,R,LOC,Y,0,0.0012
         NSEL,R,LOC,Z,0,0.001
         NPLOT
         BF,ALL,HGEN,kcoe*0.74*4.859924261E10
         ESEL,ALL
         NSLE,ALL
         ALLSEL,ALL

         TIME,(b-1)*len2/v+c*step1/v+24
         ESTIF,1.0E-50
         NLGEOM,ON
         NROPT,FULL,,ON
         AUTOTS,1
         KBC,1
         lnsrch,1
         TSRES,ERASE
         OUTPR,NSOL,LAST,
         TINTP,,,,1,,,,,
         SOLVE

         ESEL,S,LIVE
         NSLE,S         
         NPLOT
         BF,ALL,HGEN,0
         ESEL,S,LIVE
         NSLE,S
         SAVE
         *ENDDO
    Y1=0
    *ENDDO
save
Y1=0
ESEL,ALL
CMSEL,S,VNAME(3)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
EPLOT
!The third layer
ALLSEL,ALL
    *DO,b,1,3
         *DO,c,1,10
         step1=2*(len3/n)
         x1=x1+step1
         LOCAL,11,0,0.008-x1,-0.006+(b-1)*0.004,0.002,,,,,,
         NSEL,S,LOC,X,0,0.0016
         NSEL,R,LOC,Y,0,0.002
         NSEL,R,LOC,Z,0,0.001
         NPLOT
         BF,ALL,HGEN,kcoe*0.68*6.479899014E10
         ESEL,ALL
         NSLE,ALL
         ALLSEL,ALL

         TIME,(b-1)*len3/v+c*step1/v+32
         ESTIF,1.0E-50
         NLGEOM,ON
         NROPT,FULL,,ON
         AUTOTS,1
         KBC,1
         lnsrch,1
         TSRES,ERASE
         OUTPR,NSOL,LAST,
         TINTP,,,,1,,,,,
         SOLVE

         ESEL,S,LIVE
         NSLE,S         
         NPLOT
         BF,ALL,HGEN,0
         ESEL,S,LIVE
         NSLE,S
         SAVE
         *ENDDO
    x1=0
    *ENDDO
save
x1=0
ALLSEL,ALL
    *DO,b,1,3
         *DO,c,1,10
         step1=2*(len3/n)
         x1=x1+step1
         LOCAL,11,0,-0.008+x1,0.004-(b-1)*0.004,0.002,,,,,,
         NSEL,S,LOC,X,-0.0016,0
         NSEL,R,LOC,Y,0,0.002
         NSEL,R,LOC,Z,0,0.001
         NPLOT
         BF,ALL,HGEN,kcoe*0.56*6.479899014E10
         ESEL,ALL
         NSLE,ALL
         ALLSEL,ALL

         TIME,(b-1)*len3/v+c*step1/v+40
         ESTIF,1.0E-50
         NLGEOM,ON
         NROPT,FULL,,ON
         AUTOTS,1
         KBC,1
         lnsrch,1
         TSRES,ERASE
         OUTPR,NSOL,LAST,
         TINTP,,,,1,,,,,
         SOLVE

         ESEL,S,LIVE
         NSLE,S         
         NPLOT
         BF,ALL,HGEN,0
         ESEL,S,LIVE
         NSLE,S
         SAVE
         *ENDDO
    x1=0
    *ENDDO
save
x1=0
ESEL,ALL
CMSEL,S,VNAME(15)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
EPLOT
!The fourth layer
ALLSEL,ALL
    *DO,b,1,4
         *DO,c,1,10
         step1=2*(len2/n)
         Y1=Y1+step1
         LOCAL,11,0,0.008-(b-1)*0.004,-0.006+Y1,0.003,,,,,,
         NSEL,S,LOC,X,-0.002,0
         NSEL,R,LOC,Y,-0.0012,0
         NSEL,R,LOC,Z,0,0.001
         NPLOT
         BF,ALL,HGEN,kcoe*0.94*4.859924261E10
         ESEL,ALL
         NSLE,ALL
         ALLSEL,ALL

         TIME,(b-1)*len2/v+c*step1/v+48
         ESTIF,1.0E-50
         NLGEOM,ON
         NROPT,FULL,,ON
         AUTOTS,1
         KBC,1
         lnsrch,1
         TSRES,ERASE
         OUTPR,NSOL,LAST,
         TINTP,,,,1,,,,,
         SOLVE

         ESEL,S,LIVE
         NSLE,S         
         NPLOT
         BF,ALL,HGEN,0
         ESEL,S,LIVE
         NSLE,S
         SAVE
         *ENDDO
    Y1=0
    *ENDDO
save
Y1=0
ALLSEL,ALL
    *DO,b,1,4
         *DO,c,1,10
         step1=2*(len2/n)
         Y1=Y1+step1
         LOCAL,11,0,-0.008+(b-1)*0.004,0.006-Y1,0.003,,,,,,
         NSEL,S,LOC,X,0,0.002
         NSEL,R,LOC,Y,0,0.0012
         NSEL,R,LOC,Z,0,0.001
         NPLOT
         BF,ALL,HGEN,kcoe*0.91*4.859924261E10
         ESEL,ALL
         NSLE,ALL
         ALLSEL,ALL

         TIME,(b-1)*len2/v+c*step1/v+56
         ESTIF,1.0E-50
         NLGEOM,ON
         NROPT,FULL,,ON
         AUTOTS,1
         KBC,1
         lnsrch,1
         TSRES,ERASE
         OUTPR,NSOL,LAST,
         TINTP,,,,1,,,,,
         SOLVE

         ESEL,S,LIVE
         NSLE,S         
         NPLOT
         BF,ALL,HGEN,0
         ESEL,S,LIVE
         NSLE,S
         SAVE
         *ENDDO
    Y1=0
    *ENDDO
save
Y1=0

ALLSEL,ALL
*DO,r,1,50
TIME,0.4*r+64
ESTIF,1.0E-50
AUTOTS,1
KBC,1
lnsrch,1
TSRES,ERASE
OUTPR,NSOL,LAST,
TINTP,,,,1,,,,,
SOLVE
*ENDDO

ALLSEL,ALL
*DO,t,1,16
TIME,1*t+84
ESTIF,1.0E-50
AUTOTS,1
KBC,1
lnsrch,1
TSRES,ERASE
OUTPR,NSOL,LAST,
TINTP,,,,1,,,,,
SOLVE
*ENDDO
x1=0

!The fifth layer
ALLSEL,ALL
    *DO,b,1,3
         *DO,c,1,10
         step1=2*(len3/n)
         x1=x1+step1
         LOCAL,11,0,0.008-x1,-0.006+(b-1)*0.004,0.004,,,,,,
         NSEL,S,LOC,X,0,0.0016
         NSEL,R,LOC,Y,0,0.002
         NSEL,R,LOC,Z,0,0.001
         NPLOT
         BF,ALL,HGEN,kcoe*0.9*6.479899014E10
         ESEL,ALL
         NSLE,ALL
         ALLSEL,ALL

         TIME,(b-1)*len3/v+c*step1/v+100
         ESTIF,1.0E-50
         NLGEOM,ON
         NROPT,FULL,,ON
         AUTOTS,1
         KBC,1
         lnsrch,1
         TSRES,ERASE
         OUTPR,NSOL,LAST,
         TINTP,,,,1,,,,,
         SOLVE

         ESEL,S,LIVE
         NSLE,S         
         NPLOT
         BF,ALL,HGEN,0
         ESEL,S,LIVE
         NSLE,S
         SAVE
         *ENDDO
    x1=0
    *ENDDO
save
x1=0
ALLSEL,ALL
    *DO,b,1,3
         *DO,c,1,10
         step1=2*(len3/n)
         x1=x1+step1
         LOCAL,11,0,-0.008+x1,0.004-(b-1)*0.004,0.004,,,,,,
         NSEL,S,LOC,X,-0.0016,0
         NSEL,R,LOC,Y,0,0.002
         NSEL,R,LOC,Z,0,0.001
         NPLOT
         BF,ALL,HGEN,kcoe*0.75*6.479899014E10
         ESEL,ALL
         NSLE,ALL
         ALLSEL,ALL

         TIME,(b-1)*len3/v+c*step1/v+108
         ESTIF,1.0E-50
         NLGEOM,ON
         NROPT,FULL,,ON
         AUTOTS,1
         KBC,1
         lnsrch,1
         TSRES,ERASE
         OUTPR,NSOL,LAST,
         TINTP,,,,1,,,,,
         SOLVE

         ESEL,S,LIVE
         NSLE,S         
         NPLOT
         BF,ALL,HGEN,0
         ESEL,S,LIVE
         NSLE,S
         SAVE
         *ENDDO
    x1=0
    *ENDDO
save
x1=0

!The sixth layer
ALLSEL,ALL
    *DO,b,1,4
         *DO,c,1,10
         step1=2*(len2/n)
         Y1=Y1+step1
         LOCAL,11,0,0.008-(b-1)*0.004,-0.006+Y1,0.005,,,,,,
         NSEL,S,LOC,X,-0.002,0
         NSEL,R,LOC,Y,-0.0012,0
         NSEL,R,LOC,Z,0,0.001
         NPLOT
         BF,ALL,HGEN,kcoe*0.85*4.859924261E10
         ESEL,ALL
         NSLE,ALL
         ALLSEL,ALL

         TIME,(b-1)*len2/v+c*step1/v+116
         ESTIF,1.0E-50
         NLGEOM,ON
         NROPT,FULL,,ON
         AUTOTS,1
         KBC,1
         lnsrch,1
         TSRES,ERASE
         OUTPR,NSOL,LAST,
         TINTP,,,,1,,,,,
         SOLVE

         ESEL,S,LIVE
         NSLE,S         
         NPLOT
         BF,ALL,HGEN,0
         ESEL,S,LIVE
         NSLE,S
         SAVE
         *ENDDO
    Y1=0
    *ENDDO
save
Y1=0
ALLSEL,ALL
    *DO,b,1,4
         *DO,c,1,10
         step1=2*(len2/n)
         Y1=Y1+step1
         LOCAL,11,0,-0.008+(b-1)*0.004,0.006-Y1,0.005,,,,,,
         NSEL,S,LOC,X,0,0.002
         NSEL,R,LOC,Y,0,0.0012
         NSEL,R,LOC,Z,0,0.001
         NPLOT
         BF,ALL,HGEN,kcoe*0.9*4.859924261E10
         ESEL,ALL
         NSLE,ALL
         ALLSEL,ALL

         TIME,(b-1)*len2/v+c*step1/v+124
         ESTIF,1.0E-50
         NLGEOM,ON
         NROPT,FULL,,ON
         AUTOTS,1
         KBC,1
         lnsrch,1
         TSRES,ERASE
         OUTPR,NSOL,LAST,
         TINTP,,,,1,,,,,
         SOLVE

         ESEL,S,LIVE
         NSLE,S         
         NPLOT
         BF,ALL,HGEN,0
         ESEL,S,LIVE
         NSLE,S
         SAVE
         *ENDDO
    Y1=0
    *ENDDO
save
Y1=0

*DO,p,1,50
ALLSEL,ALL
TIME,0.6*p+132
ESTIF,1.0E-50
AUTOTS,1
KBC,1
lnsrch,1
TSRES,ERASE
OUTPR,NSOL,LAST,
TINTP,,,,1,,,,,
SOLVE
*ENDDO

!The seventh layer
ALLSEL,ALL
    *DO,b,1,3
         *DO,c,1,10
         step1=2*(len3/n)
         x1=x1+step1
         LOCAL,11,0,0.008-x1,-0.006+(b-1)*0.004,0.006,,,,,,
         NSEL,S,LOC,X,0,0.0016
         NSEL,R,LOC,Y,0,0.002
         NSEL,R,LOC,Z,0,0.001
         NPLOT
         BF,ALL,HGEN,kcoe*1.04*6.479899014E10
         ESEL,ALL
         NSLE,ALL
         ALLSEL,ALL

         TIME,(b-1)*len3/v+c*step1/v+162
         ESTIF,1.0E-50
         NLGEOM,ON
         NROPT,FULL,,ON
         AUTOTS,1
         KBC,1
         lnsrch,1
         TSRES,ERASE
         OUTPR,NSOL,LAST,
         TINTP,,,,1,,,,,
         SOLVE

         ESEL,S,LIVE
         NSLE,S         
         NPLOT
         BF,ALL,HGEN,0
         ESEL,S,LIVE
         NSLE,S
         SAVE
         *ENDDO
    x1=0
    *ENDDO
save
x1=0
ALLSEL,ALL
    *DO,b,1,3
         *DO,c,1,10
         step1=2*(len3/n)
         x1=x1+step1
         LOCAL,11,0,-0.008+x1,0.004-(b-1)*0.004,0.006,,,,,,
         NSEL,S,LOC,X,-0.0016,0
         NSEL,R,LOC,Y,0,0.002
         NSEL,R,LOC,Z,0,0.001
         NPLOT
         BF,ALL,HGEN,kcoe*1.04*6.479899014E10
         ESEL,ALL
         NSLE,ALL
         ALLSEL,ALL

         TIME,(b-1)*len3/v+c*step1/v+170
         ESTIF,1.0E-50
         NLGEOM,ON
         NROPT,FULL,,ON
         AUTOTS,1
         KBC,1
         lnsrch,1
         TSRES,ERASE
         OUTPR,NSOL,LAST,
         TINTP,,,,1,,,,,
         SOLVE

         ESEL,S,LIVE
         NSLE,S         
         NPLOT
         BF,ALL,HGEN,0
         ESEL,S,LIVE
         NSLE,S
         SAVE
         *ENDDO
    x1=0
    *ENDDO
save


ESEL,ALL
CMSEL,S,VNAME(8)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
EPLOT
*DO,d,1,2
   *DO,e,1,10
      LOCAL,11,1,0,0,0,,,,,1,1
      NSEL,S,LOC,X,0.008+(d-1)*0.001,0.008+d*0.001
      NSEL,R,LOC,Y,-ASIN(0.6)+7.2*(e-1)-ang1,-ASIN(0.6)+7.2*e+ang1
      NSEL,R,LOC,Z,0.002,0.011
      NPLOT
      BF,ALL,HGEN,kcoe*0.78*5.091624385E10
      ESEL,ALL
      NSLE,ALL
      ALLSEL,ALL

      TIME,0.011309733*(d-1)/v+0.0011309733*e/v+178
      ESTIF,1.0E-50
      NLGEOM,ON
      NROPT,FULL,,ON
      AUTOTS,1

      KBC,1
      lnsrch,1
      TSRES,ERASE
         OUTPR,NSOL,LAST,
      TINTP,,,,1,,,,,
      SOLVE

      ESEL,S,LIVE
      NSLE,S         
      NPLOT
      BF,ALL,HGEN,0
      ESEL,S,LIVE
      NSLE,S
   *ENDDO
   SAVE
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(6)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
EPLOT
*DO,f,1,4
   *DO,g,1,15
      LOCAL,11,1,0,0,0,,,,,1,1,
      NSEL,R,LOC,X,0.006+(f-1)*0.001,0.006+f*0.001
      NSEL,R,LOC,Y,ASIN(0.6)+7.2*(g-1)-ang1,ASIN(0.6)+7.2*g+ang1
      NSEL,R,LOC,Z,0.0020,0.011
      NPLOT
      BF,ALL,HGEN,kcoe*0.8*5.527836649E10
      ESEL,ALL
      NSLE,ALL
      ALLSEL,ALL

      TIME,0.0169646*(f-1)/v+0.0011309733*g/v+181.769911
      ESTIF,1.0E-50
      NLGEOM,ON
      NROPT,FULL,,ON
      AUTOTS,1

      KBC,1
      lnsrch,1
      TSRES,ERASE
         OUTPR,NSOL,LAST,
      TINTP,,,,1,,,,,
      SOLVE

      ESEL,S,LIVE
      NSLE,S         
      NPLOT
      BF,ALL,HGEN,0
      ESEL,S,LIVE
      NSLE,S
   *ENDDO
   SAVE
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(7)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
EPLOT

*DO,h,1,2
   *DO,i,1,10
      LOCAL,11,1,0,0,0,,,,,1,1,
      NSEL,R,LOC,X,0.008+(h-1)*0.001,0.008+h*0.001
      NSEL,R,LOC,Y,180-ASIN(0.6)+7.2*(i-1)-ang1,180-ASIN(0.6)+7.2*i+ang1
      NSEL,R,LOC,Z,0.0020,0.011
      NPLOT
      BF,ALL,HGEN,kcoe*0.78*5.091624385E10
      ESEL,ALL
      NSLE,ALL
      ALLSEL,ALL

      TIME,0.011309733*(h-1)/v+0.0011309733*i/v+193.0796443
      ESTIF,1.0E-50
      NLGEOM,ON
      NROPT,FULL,,ON
      AUTOTS,1

      KBC,1
      lnsrch,1
      TSRES,ERASE
         OUTPR,NSOL,LAST,
      TINTP,,,,1,,,,,
      SOLVE

      ESEL,S,LIVE
      NSLE,S         
      NPLOT
      BF,ALL,HGEN,0
      ESEL,S,LIVE
      NSLE,S
   *ENDDO
   SAVE
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(5)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
EPLOT

*DO,j,1,4
   *DO,k,1,15
      LOCAL,11,1,0,0,0,,,,,1,1,
      NSEL,R,LOC,X,0.006+(j-1)*0.001,0.006+j*0.001
      NSEL,R,LOC,Y,180+ASIN(0.6)+7.2*(k-1)-ang1,180+ASIN(0.6)+7.2*k+ang1
      NSEL,R,LOC,Z,0.0020,0.011
      NPLOT
      BF,ALL,HGEN,kcoe*0.76*5.527836649E10
      ESEL,ALL
      NSLE,ALL
      ALLSEL,ALL

      TIME,0.0169646*(j-1)/v+0.0011309733*k/v+196.8495553
      ESTIF,1.0E-50
      NLGEOM,ON
      NROPT,FULL,,ON
      AUTOTS,1

      KBC,1
      lnsrch,1
      TSRES,ERASE
         OUTPR,NSOL,LAST,
      TINTP,,,,1,,,,,
      SOLVE

      ESEL,S,LIVE
      NSLE,S         
      NPLOT
      BF,ALL,HGEN,0
      ESEL,S,LIVE
      NSLE,S
   *ENDDO
   SAVE
*ENDDO
SAVE




!the first map
ESEL,ALL
CMSEL,S,VNAME(26)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.010,0.011
   NSEL,R,LOC,Y,-6*(m-1)-30+ang1,-6*m-30-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.74*4.25303639E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.0013823*5/6*m/v+208.1592886
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(43)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.010,0.011
   NSEL,R,LOC,Y,-6*(m-1)-150+ang1,-6*m-150-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.78*4.25303639E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.0013823*5/6*m/v+208.1592886+1.919861111
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(32)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.010,0.011
   NSEL,R,LOC,Y,-6*(m-1)+90+ang1,-6*m+90-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.78*4.25303639E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.0013823*5/6*m/v+208.1592886+1.919861111*2
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(44)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.010,0.011
   NSEL,R,LOC,Y,6*(m-1)-30+ang1,6*m-30-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.78*4.25303639E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.0013823*5/6*m/v+208.1592886+1.919861111*3
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(25)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.010,0.011
   NSEL,R,LOC,Y,6*(m-1)+90+ang1,6*m+90-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.78*4.25303639E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.0013823*5/6*m/v+208.1592886+1.919861111*4
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(31)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.010,0.011
   NSEL,R,LOC,Y,6*(m-1)-150+ang1,6*m-150-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.78*4.25303639E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.0013823*5/6*m/v+208.1592886+1.919861111*5
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE


!the second map
ESEL,ALL
CMSEL,S,VNAME(18)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.011,0.012
   NSEL,R,LOC,Y,-6*(m-1)-30+ang1,-6*m-30-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.79*4.316831936E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.001507964*5/6*m/v+219.6784553
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(36)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.011,0.012
   NSEL,R,LOC,Y,-6*(m-1)-150+ang1,-6*m-150-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.79*4.316831936E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.001507964*5/6*m/v+219.6784553+2.094394444
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(11)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.011,0.012
   NSEL,R,LOC,Y,-6*(m-1)+90+ang1,-6*m+90-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.79*4.316831936E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.001507964*5/6*m/v+219.6784553+2.094394444*2
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(35)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.011,0.012
   NSEL,R,LOC,Y,6*(m-1)-30+ang1,6*m-30-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.79*4.316831936E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.001507964*5/6*m/v+219.6784553+2.094394444*3
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(17)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.011,0.012
   NSEL,R,LOC,Y,6*(m-1)+90+ang1,6*m+90-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.79*4.316831936E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.0013823*5/6*m/v+219.6784553+2.094394444*4
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(12)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.011,0.012
   NSEL,R,LOC,Y,6*(m-1)-150+ang1,6*m-150-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.79*4.316831936E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.001507964*5/6*m/v+219.6784553+2.094394444*5
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

!the third map
ESEL,ALL
CMSEL,S,VNAME(20)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.012,0.013
   NSEL,R,LOC,Y,-6*(m-1)-30+ang1,-6*m-30-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.78*4.381584415E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.001633628*5/6*m/v+232.244822
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(38)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.012,0.013
   NSEL,R,LOC,Y,-6*(m-1)-150+ang1,-6*m-150-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.78*4.381584415E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.001633628*5/6*m/v+232.244822+2.268927778
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(13)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.012,0.013
   NSEL,R,LOC,Y,-6*(m-1)+90+ang1,-6*m+90-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.78*4.381584415E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.001633628*5/6*m/v+232.244822+2.268927778*2
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(37)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.012,0.013
   NSEL,R,LOC,Y,6*(m-1)-30+ang1,6*m-30-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.78*4.381584415E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.001633628*5/6*m/v+232.244822+2.268927778*3
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(19)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.012,0.013
   NSEL,R,LOC,Y,6*(m-1)+90+ang1,6*m+90-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.78*4.381584415E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.001633628*5/6*m/v+232.244822+2.268927778*4
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(14)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.012,0.013
   NSEL,R,LOC,Y,6*(m-1)-150+ang1,6*m-150-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.78*4.381584415E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.0013823*5/6*m/v+232.244822+2.268927778*5
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE


!the fourth map
ESEL,ALL
CMSEL,S,VNAME(22)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.013,0.014
   NSEL,R,LOC,Y,-6*(m-1)-30+ang1,-6*m-30-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.8*4.447308181E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.001759291*5/6*m/v+245.8583887
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(40)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.013,0.014
   NSEL,R,LOC,Y,-6*(m-1)-150+ang1,-6*m-150-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.8*4.447308181E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.001759291*5/6*m/v+245.8583887+2.443459722
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(27)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.013,0.014
   NSEL,R,LOC,Y,-6*(m-1)+90+ang1,-6*m+90-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.8*4.447308181E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.001759291*5/6*m/v+245.8583887+2.443459722*2
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(39)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.013,0.014
   NSEL,R,LOC,Y,6*(m-1)-30+ang1,6*m-30-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.8*4.447308181E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.001759291*5/6*m/v+245.8583887+2.443459722*3
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(21)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.013,0.014
   NSEL,R,LOC,Y,6*(m-1)+90+ang1,6*m+90-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.8*4.447308181E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.001759291*5/6*m/v+245.8583887+2.443459722*4
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(28)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.013,0.014
   NSEL,R,LOC,Y,6*(m-1)-150+ang1,6*m-150-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.8*4.447308181E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.0013823*5/6*m/v+245.8583887+2.443459722*5
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE


!the fifth map
ESEL,ALL
CMSEL,S,VNAME(24)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.014,0.015
   NSEL,R,LOC,Y,-6*(m-1)-30+ang1,-6*m-30-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.8*4.514017804E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.001884955*5/6*m/v+260.519147
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(42)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.014,0.015
   NSEL,R,LOC,Y,-6*(m-1)-150+ang1,-6*m-150-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.8*4.514017804E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.001884955*5/6*m/v+260.519147+2.617993056
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(29)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.014,0.015
   NSEL,R,LOC,Y,-6*(m-1)+90+ang1,-6*m+90-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.8*4.514017804E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.001884955*5/6*m/v+260.519147+2.617993056*2
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(41)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.014,0.015
   NSEL,R,LOC,Y,6*(m-1)-30+ang1,6*m-30-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.8*4.514017804E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.001884955*5/6*m/v+260.519147+2.617993056*3
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(23)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.014,0.015
   NSEL,R,LOC,Y,6*(m-1)+90+ang1,6*m+90-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.8*4.514017804E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.001884955*5/6*m/v+260.519147+2.617993056*4
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(30)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.014,0.015
   NSEL,R,LOC,Y,6*(m-1)-150+ang1,6*m-150-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.8*4.514017804E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.001884955*5/6*m/v+260.519147+2.617993056*5
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE


!the sixth map
ESEL,ALL
CMSEL,S,VNAME(16)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.015,0.016
   NSEL,R,LOC,Y,-6*(m-1)-30+ang1,-6*m-30-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.8*4.739718694E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.002010619*5/6*m/v+276.2271053
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(34)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.015,0.016
   NSEL,R,LOC,Y,-6*(m-1)-150+ang1,-6*m-150-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.8*4.739718694E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.002010619*5/6*m/v+276.2271053+2.792526389
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(9)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.015,0.016
   NSEL,R,LOC,Y,-6*(m-1)+90+ang1,-6*m+90-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.8*4.739718694E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.002010619*5/6*m/v+276.2271053+2.792526389*2
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(33)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.015,0.016
   NSEL,R,LOC,Y,6*(m-1)-30+ang1,6*m-30-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.8*4.739718694E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.002010619*5/6*m/v+276.2271053+2.792526389*3
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(4)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.015,0.016
   NSEL,R,LOC,Y,6*(m-1)+90+ang1,6*m+90-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.8*4.739718694E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.002010619*5/6*m/v+276.2271053+2.792526389*4
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE

ESEL,ALL
CMSEL,S,VNAME(10)
ESLN,S,1,ALL
EALIVE,ALL
ESEL,S,LIVE
NPLOT
   *DO,m,1,10
   LOCAL,11,1,0,0,0,,,,,1,1
   NSEL,S,LOC,X,0.015,0.016
   NSEL,R,LOC,Y,6*(m-1)-150+ang1,6*m-150-ang1
   NSEL,R,LOC,Z,0.0020,0.011
   NPLOT
   BF,ALL,HGEN,kcoe*0.8*4.739718694E10
   ESEL,ALL
   NSLE,ALL
   ALLSEL,ALL

   TIME,0.002010619*5/6*m/v+276.2271053+2.792526389*5
   ESTIF,1.0E-50
   NLGEOM,ON
   NROPT,FULL,,ON
   AUTOTS,1
   KBC,1
   lnsrch,1
   TSRES,ERASE
   OUTPR,NSOL,LAST,
   TINTP,,,,1,,,,,
   SOLVE

   ESEL,S,LIVE
   NSLE,S         
   NPLOT
   BF,ALL,HGEN,0
   ESEL,S,LIVE
   NSLE,S
   *ENDDO
   SAVE   
*ENDDO
SAVE


ESEL,S,LIVE
NSLE,S         
NPLOT
BF,ALL,HGEN,0
ESEL,S,LIVE
NSLE,S
*DO,z,1,40
ALLSEL,ALL
TIME,0.5*z+292.9822636
ESTIF,1.0E-50
AUTOTS,1
KBC,1
lnsrch,1
TSRES,ERASE
OUTPR,NSOL,LAST,
TINTP,,,,1,,,,,
SOLVE
*ENDDO
save

*DO,w,1,40
ALLSEL,ALL
TIME,1*w+312.9822636
ESTIF,1.0E-50
AUTOTS,1
KBC,1
lnsrch,1
TSRES,ERASE
OUTPR,NSOL,LAST,
TINTP,,,,1,,,,,
SOLVE
*ENDDO
save

*DO,t,1,40
ALLSEL,ALL
TIME,1.5*t+352.9822636
ESTIF,1.0E-50
AUTOTS,1
KBC,1
lnsrch,1
TSRES,ERASE
OUTPR,NSOL,LAST,
TINTP,,,,1,,,,,
SOLVE
*ENDDO
save

*DO,x,1,40
ALLSEL,ALL
TIME,2*x+412.9822636
ESTIF,1.0E-50
AUTOTS,1
KBC,1
lnsrch,1
TSRES,ERASE
OUTPR,NSOL,LAST,
TINTP,,,,1,,,,,
SOLVE
*ENDDO
save
FINISH
!temperature of the fourth layer of the substrate will increased dramaticaly , !so we cool it for 40 second
!after three layers of substrate

[ 本帖最后由 sunclark 于 2007-3-7 14:01 编辑 ]

sunclark

金属零件的无模成形是快速成型技术 (Rapid Prototyping, RP) 的重要发展方向。等离子熔积成形中，因温度分布不均产生热应力，会导致成形件翘曲变形和裂纹等缺陷。预测金属零件成形过程的温度场分布及其变化是应力应变分析的前提和成形性控制的基础。由温度场进而预报应力场的分布和变化，及时采取措施以减小和降低某些应力集中部位的应力累积，从而避免成形件的冷裂和热裂，最终改善和提高零件的成形性与成品率。等离子熔积无模成形属高能束急速加热快速凝固过程，温度场复杂且实测困难；用有限元法可方便快捷模拟任意工艺条件下成形件的温度场，从而揭示三维瞬态温度场与应力应变场的时空演变历程，对成形工艺的优化设计具有指导意义。

[ 本帖最后由 sunclark 于 2007-2-9 22:43 编辑 ]

sunclark

现有直接制造高密度零件的快速成形技术如激光近终成形 (Laser Engineered Net Shaping, LENS) 及等离子熔积过程的温度场模拟，多采用简单的平板或对称有限元模型、恒定热输入，未考虑熔积路径及冷却条件的优化，且以二维模拟居多，而复杂形状零件的高能束直接快速成形过程三维瞬态温度场的有限元分析还少见报道。论文基于有限元软件ANSYS开发材料逐步叠加和填充的 ( ANSYS Parametric Design Language, APDL) 子程序，优化设计零件的拓扑结构及网格划分，并用“生死单元”和自适应网格技术实现了复杂形状零件三维“虚拟增材制造”过程的温度场模拟。

sunclark

这篇文章定于07年7月刊发，请大家多多引用，非常感谢
计算中对零件基座预热温度从理论上进行了推导，并特别处理了复杂的热对流和辐射边界条件。自认为是目前模拟计算高能束快速成形零件做的工作量最大的一个！但由于保密原因，对其中网格自适应没有提及，欢迎交流！
wushengchuan@gmail.com
因上传文件大小限制，无法上传该篇文章全文，抱歉！
王桂兰，吴圣川，张海鸥. 复杂零件等离子熔积无模成形的温度场模拟[J]. 焊接学报，2007

[ 本帖最后由 sunclark 于 2007-2-9 22:54 编辑 ]

zjjxsjy

好文章，谢谢了。能否发给我？
zjjxsjy@163.com

jinqinghai

楼主的分析技术不错，奉献精神可嘉！
顶一个!

sunclark

不知3M可以不可以传上来！
谢谢斑竹的鼓励

[ 本帖最后由 sunclark 于 2007-3-1 10:02 编辑 ]

kitty827

谢谢楼主的分享，我的课题正好与这个相关，是做激光烧结快速成型的，很感谢楼主

xbzhang

先看看，好像很厉害的命令流哦。
不知能不能正常运行？

xbzhang

把命令流拷贝成TXT文件读入ANSYS10，观察输出窗口，发现很多小主元问题的提示。好像没有给约束。
楼主在计算时有没有这些小主元的提示，是否影响结果？下面给出这些提示的截图？
我现在还没有算完，想看看最后的结果如何？

xbzhang

ANSYS算到下午13：46，突然退出，不给算了，不知何故。

sunclark

没有问题的,可以正常计算.贴出的命令流是我最后的版本,应该没有问题的.
请xbzhang把error文件中的错误信息拷贝一下,我看看,是什么原因?!

xbzhang

我把错误文件的（.err）开头拷贝过来了，以后就全是小主元的警告，后来就退出了。楼主可以看一下计算时间，我想可能版上的命令流有遗漏。

/COM,ANSYS RELEASE 10.0    UP20050718       08:26:35    03/03/2007

*** WARNING ***                         CP =      12.797   TIME= 08:35:36
Assembly V0 is deleted.                                                

*** WARNING ***                         CP =      14.344   TIME= 08:35:38
Assembly V1 is deleted.                                                

*** WARNING ***                         CP =      27.906   TIME= 08:35:52
Shape testing revealed that 8 of the 6734 new or modified elements      
violate shape warning limits.  To review test results, please see the   
output file or issue the CHECK command.                                 

*** WARNING ***                         CP =      28.109   TIME= 08:35:53
After NUMMRG,NODE, node 6041 (and possibly others) is associated with   
more than one solid model entity.  Future commands which depend on the  
node to solid model connectivity (meshing, mesh clearing, solid BC      
transfer, etc.) may not operate properly.  NUMMRG,KP may correct this   
problem.                                                               

*** WARNING ***                         CP =      28.188   TIME= 08:35:53
After NUMMRG,KP, node 6041 (and possibly others) is associated with     
more than one solid model entity.  Future commands which depend on the  
node to solid model connectivity (meshing, mesh clearing, solid BC      
transfer, etc.) may not operate properly.                              

*** WARNING ***                         CP =      28.234   TIME= 08:35:53
After NUMMRG,KP, node 6041 (and possibly others) is associated with     
more than one solid model entity.  Future commands which depend on the  
node to solid model connectivity (meshing, mesh clearing, solid BC      
transfer, etc.) may not operate properly.                              

*** WARNING ***                         CP =      29.781   TIME= 08:35:55
The requested EPLOT command cannot be executed.                        
Either elements do not exist or the selected elements                  
are not supported by PowerGraphics.                                    
To correct this situation, try selecting different                     
elements,etc, or turning off PowerGraphics.                             

*** WARNING ***                         CP =      29.969   TIME= 08:35:56
Large deflection effects not included for element SOLID70.              

*** WARNING ***                         CP =      30.000   TIME= 08:35:56
Previous testing revealed that 8 of the 9428 selected elements violate  
shape warning limits.  To review warning messages, please see the      
output or error file, or issue the CHECK command.                       

*** WARNING ***                         CP =      30.062   TIME= 08:35:56
The program chosen initial timestep/load-factor is arbitrary.  It is   
necessary for the user to supply a suitable initial                     
timestep/load-factor through the NSUB or DELTIM command for            
convergence and overall efficiency.                                    

*** WARNING ***                         CP =      30.500   TIME= 08:35:58
Material property C of material 1 of element 3611 is evaluated at a     
temperature of 20, which is below the supplied temperature range.      
Temperature range checking terminates.                                 

*** WARNING ***                         CP =      49.312   TIME= 08:36:18
There are 11237 small equation solver pivot terms.  May be a            
transitional effect caused by the Newton-Raphson logic.                 

*** WARNING ***   SUPPRESSED MESSAGE    CP =      68.422   TIME= 08:36:37
There are 11237 small equation solver pivot terms.  May be a            
transitional effect caused by the Newton-Raphson logic.                 

最后的警告信息是：

*** WARNING ***                         CP =   17690.422   TIME= 13:43:03
There are 4896 small equation solver pivot terms.  May be a            
transitional effect caused by the Newton-Raphson logic.                 

*** WARNING ***                         CP =   17693.734   TIME= 13:43:06
There are 4896 small equation solver pivot terms.  May be a            
transitional effect caused by the Newton-Raphson logic.  

总共的错误提示为8637行，在UltraEdit中显示的。

sunclark

如果再有问题，请贴上来，也可能是您的计算机有些于设定没有做。
这个结果文件有几个G，所以，你先要用命令/config,fsplit,300  
/config,nres,2000   
/config,nbuf,16
为提高求解速度在计算前再手动输入EQSLV,PCG,1E-4，使用PCG求解器

xbzhang

再请问一下你用的是ANSYS的什么产品，多学科吗？就是multiphysics

sunclark

对，是multiphysics
选择不同学科，可能运行会有问题

menyueqi

真是很好的帖子，多跟楼主学习了呀。分析技术很好呀。

liny0901

华科大材料系的牛人!!!!!!

wangdswangds

:) :) :) :)

xuyong415

太强了，哎看不东