2# jiandan32
/PREP7
!*
ET,1,SOLID117
!*
ET,2,SOLID117
!*
ET,3,SOLID117
!*
ET,4,SOLID117
!*
ET,5,SOLID117
!*
ET,6,SOLID117
!*
ET,7,SOLID117
!*
ET,8,SOLID117
!*
ET,9,SOLID117
!*
KEYOPT,1,1,0
KEYOPT,1,2,0
KEYOPT,1,5,0
!*
!*
EMUNIT,MKS
!*
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,MURX,1,,1
*CSET,1,2, 1,2,
MPCOPY, ,1,2
TBCOPY,ALL,1,2
!*
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDE,MURX,2
MPDATA,MURX,2,,1
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,RSVX,2,,3.0e-8
*CSET,1,2, 1,3,
MPCOPY, ,1,3
TBCOPY,ALL,1,3
!*
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDE,MURX,3
MPDATA,MURX,3,,5000
TB,BH,4,1,25,
TBTEMP,0
TBPT,,100,0.46512
TBPT,,200,0.72993
TBPT,,300,0.9009
TBPT,,400,1.0204
TBPT,,500,1.1086
TBPT,,600,1.1765
TBPT,,700,1.2302
TBPT,,800,1.2739
TBPT,,900,1.31
TBPT,,1000,1.3405
TBPT,,1400,1.4257
TBPT,,1800,1.4778
TBPT,,2200,1.5131
TBPT,,2600,1.5385
TBPT,,3000,1.5576
TBPT,,3400,1.5726
TBPT,,3800,1.5847
TBPT,,4200,1.5945
TBPT,,4600,1.6028
TBPT,,5000,1.6098
TBPT,,7000,1.6332
TBPT,,9000,1.6465
TBPT,,11000,1.6551
TBPT,,13000,1.6611
TBPT,,15000,1.6656
!*
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,RSVX,4,,2.0e-7 /// 前处理之前的部分
/PREP7
FLST,2,18,6,ORDE,2
FITEM,2,1
FITEM,2,-18
WPSTYLE,,,,,,,,1
CYLIND,65, ,0,50,0,120,
CYLIND,75,65,0,50,0,120,
CYLIND,150,75,0,50,0,120,
CYLIND,190,150,0,50,0,120,
BLOCK,-20,20,0,200,0,50,
wprot,-60
BLOCK,-20,20,0,200,0,50,
FLST,2,2,6,ORDE,2
FITEM,2,3
FITEM,2,-4
FLST,3,4,5,ORDE,4
FITEM,3,28
FITEM,3,-29
FITEM,3,34
FITEM,3,-35
VSBA,P51X,P51X
FLST,2,4,6,ORDE,4
FITEM,2,11
FITEM,2,13
FITEM,2,-14
FITEM,2,16
VDELE,P51X, , ,1
FLST,2,2,6,ORDE,2
FITEM,2,5
FITEM,2,-6
VSBA,P51X, 8
FLST,2,2,6,ORDE,2
FITEM,2,3
FITEM,2,11
VDELE,P51X, , ,1
FLST,2,10,6,ORDE,8
FITEM,2,1
FITEM,2,-2
FITEM,2,4
FITEM,2,7
FITEM,2,-10
FITEM,2,12
FITEM,2,-13
FITEM,2,15
VGLUE,P51X
BLOCK,-35,-20,120,140,0,50,
BLOCK,20,35,120,140,0,50,
wprot,60
BLOCK,20,35,120,140,0,50,
BLOCK,-35,-20,120,140,0,50,
FLST,2,7,6,ORDE,5
FITEM,2,2
FITEM,2,-4
FITEM,2,6
FITEM,2,-8
FITEM,2,14
VOVLAP,P51X
BLOCK,-35,35,120,140,50,65,
BLOCK,-35,35,120,140,-15,0,
wprot,-60
BLOCK,-35,35,120,140,-15,0,
BLOCK,-35,35,120,140,50,65,
FLST,2,18,6,ORDE,4
FITEM,2,1
FITEM,2,-14
FITEM,2,16
FITEM,2,-19
VGLUE,P51X
NUMCMP,ALL///建模完毕
/PREP7
FLST,5,2,6,ORDE,2
FITEM,5,1
FITEM,5,10
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 4, , 1, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
FLST,5,3,6,ORDE,3
FITEM,5,5
FITEM,5,14
FITEM,5,16
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 1, , 1, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
FLST,5,5,6,ORDE,5
FITEM,5,4
FITEM,5,6
FITEM,5,8
FITEM,5,17
FITEM,5,-18
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 3, , 1, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
CM,_Y,VOLU
VSEL, , , , 13
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 2, , 2, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
CM,_Y,VOLU
VSEL, , , , 7
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 2, , 3, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
CM,_Y,VOLU
VSEL, , , , 15
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 2, , 4, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
CM,_Y,VOLU
VSEL, , , , 9
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 2, , 5, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
CM,_Y,VOLU
VSEL, , , , 12
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 2, , 6, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
CM,_Y,VOLU
VSEL, , , , 2
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 2, , 7, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
CM,_Y,VOLU
VSEL, , , , 11
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 2, , 8, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
CM,_Y,VOLU
VSEL, , , , 3
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 2, , 9, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
ALLSEL,ALL
FLST,5,189,4,ORDE,2
FITEM,5,1
FITEM,5,-189
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1,10, , , , , , ,1
!*
/REPLO
MSHAPE,1,3D
MSHKEY,0
!*
FLST,5,18,6,ORDE,2
FITEM,5,1
FITEM,5,-18
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
ALLSEL,ALL ///分完网格之后的部分
/SOL
FLST,2,9,5,ORDE,7
FITEM,2,24
FITEM,2,27
FITEM,2,32
FITEM,2,53
FITEM,2,58
FITEM,2,64
FITEM,2,-67
DA,P51X,AZ,0,0
ALLSEL,ALL
!*
ANTYPE,4
!*
TRNOPT,FULL
!*
EQSLV,SPAR,1E-9
*SET,T,1E-8
*SET,C,0
*SET,N,60
*SET,CON,0.002/N
*SET,PI,2*ASIN(1)
*SET,I,1760000*SIN(T)
ESEL,S,TYPE,,2
FLST,2,220,2,ORDE,2
FITEM,2,8941
FITEM,2,-9160
BFE,P51X,JS,1,I,0,0, ,
ESEL,S,TYPE,,3
FLST,2,166,2,ORDE,2
FITEM,2,6531
FITEM,2,-6696
BFE,P51X,JS,1,0,0,-I, ,
ESEL,S,TYPE,,4
FLST,2,232,2,ORDE,2
FITEM,2,10248
FITEM,2,-10479
BFE,P51X,JS,1,-I,0,0, ,
ESEL,S,TYPE,,5
FLST,2,165,2,ORDE,2
FITEM,2,7770
FITEM,2,-7934
BFE,P51X,JS,1,0,0,I, ,
ESEL,S,TYPE,,6
FLST,2,221,2,ORDE,2
FITEM,2,8720
FITEM,2,-8940
BFE,P51X,JS,1,-I*SIN(PI/6),I*COS(PI/6),0, ,
ESEL,S,TYPE,,7
FLST,2,154,2,ORDE,2
FITEM,2,1937
FITEM,2,-2090
BFE,P51X,JS,1,0,0,-I, ,
ESEL,S,TYPE,,8
FLST,2,202,2,ORDE,2
FITEM,2,8518
FITEM,2,-8719
BFE,P51X,JS,1,I*SIN(PI/6),-I*COS(PI/6),0, ,
ESEL,S,TYPE,,9
FLST,2,160,2,ORDE,2
FITEM,2,2091
FITEM,2,-2250
BFE,P51X,JS,1,0,0,I, ,
ALLSEL,ALL
NEQIT,30
*CREATE,LOAD
TIME,T
T=T+CON
C=C+1
OUTRES,ALL,1
*IF,C,EQ,((N*.75)+1),THEN
OUTPR, ,1
*ELSEIF,C,EQ,(N+1),THEN
OUTPR, ,1
*ELSE
OUTPR, ,0
*ENDIF
SOLVE
*END
*DO,I,1,61
*USE,LOAD
*ENDDO
FINISH |