/PREP7
K,1,-1.7,0,0,
K,2,-1.7,1.6,0,
K,3,0,7,0,
K,4,1.7,7,0,
K,5,1.7,5.4,0,
K,6,0,1.4,0,
K,7,0,4.4,0,
KWPAVE, 5
CSWPLA,11,1,1,1,
K,8,4.14,-8,0,
CSYS,0
KWPAVE, 8
CSYS,4
K,9,0,-1.8,0,
KWPAVE, 9
CSWPLA,12,1,1,1,
K,10,4.11138,195,0,
CSYS,0
K,11,-1.22,3.27,0,
K,12,-0.82,3.79,,
KWPAVE, 7
CSWPLA,12,1,1,1,
K,13,3.55,-15,,
KWPAVE, 12
CSWPLA,11,1,1,1,
K,14,4.2,-15,,
CSYS,0
K,15,3.54,3.09,,
k,16,0.66,2.34,,
k,17,0 ,0,0
LSTR, 17, 1
LSTR, 1, 2
LSTR, 17, 6
LSTR, 7, 3
LSTR, 3, 4
LSTR, 4, 5
LSTR, 5, 8
LSTR, 8, 9
LSTR, 9, 10
LSTR, 7, 13
LSTR, 12, 14
LARC, 2, 12, 11
LARC, 6, 10, 16
LARC, 14, 13, 15
FLST,2,14,4
FITEM,2,4
FITEM,2,5
FITEM,2,6
FITEM,2,7
FITEM,2,8
FITEM,2,9
FITEM,2,13
FITEM,2,3
FITEM,2,1
FITEM,2,2
FITEM,2,12
FITEM,2,10
FITEM,2,14
FITEM,2,11
AL,P51X
VOFFST,1,60, ,
ET,1,FLUID142
*SET,H,60
*SET,H1,20
*SET,H2,H-H1
k,100,(kx(7)+kx(12))/2-0.12,(ky(7)+ky(12))/2-0.05,kz(7)
bsplin,7,100,12
al,14,11,10,43
vext,17,17,,,,kz(29)
wpoffs,,,-H2
vsbw,2,,,1
vdele,4,,,1
vglue,all
LESIZE,5, , ,5, , , , ,1
LESIZE,14, , ,5, , , , ,1
LESIZE,1, , ,5, , , , ,1
LESIZE,8, , ,5, , , , ,1
LESIZE,7, , ,10, , , , ,1
LESIZE,10, , ,10, , , , ,1
LESIZE,30, , ,50, , , , ,1
FLST,5,2,6,ORDE,2
FITEM,5,2
FITEM,5,-3
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VSWEEP,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
ALLSEL,ALL
FLST,5,2,5,ORDE,2
FITEM,5,1
FITEM,5,17
ASEL,S, , ,P51X
NSLA,S,1
FINISH
/SOL
FLST,5,212,1,ORDE,4
FITEM,5,119
FITEM,5,-210
FITEM,5,4073
FITEM,5,-4192
CM,_Y,NODE
NSEL,R, , ,P51X
CM,_Y1,NODE
CMSEL,S,_Y
CMDELE,_Y
!*
/GO
!*
D,_Y1,PRES,0
!*
CMDELE,_Y1
ASEL,S, , , 2
NSLA,S,1
FLST,5,184,1,ORDE,32
FITEM,5,1
FITEM,5,-92
FITEM,5,4193
FITEM,5,-4197
FITEM,5,4247
FITEM,5,4297
FITEM,5,-4299
FITEM,5,4349
FITEM,5,-4358
FITEM,5,4408
FITEM,5,-4412
FITEM,5,4462
FITEM,5,-4471
FITEM,5,4521
FITEM,5,-4530
FITEM,5,4580
FITEM,5,-4583
FITEM,5,4633
FITEM,5,-4637
FITEM,5,4687
FITEM,5,-4689
FITEM,5,4739
FITEM,5,-4749
FITEM,5,4778
FITEM,5,4796
FITEM,5,-4800
FITEM,5,4822
FITEM,5,-4831
FITEM,5,4858
FITEM,5,-4862
FITEM,5,4895
FITEM,5,-4898
CM,_Y,NODE
NSEL,R, , ,P51X
CM,_Y1,NODE
CMSEL,S,_Y
CMDELE,_Y
/GO
D,_Y1,PRES,21000000
CMDELE,_Y1
!*
FLST,5,16,5,ORDE,6
FITEM,5,4
FITEM,5,-12
FITEM,5,18
FITEM,5,-22
FITEM,5,26
FITEM,5,-27
ASEL,S, , ,P51X
NSLA,S,1
FLST,5,3928,1,ORDE,10
FITEM,5,93
FITEM,5,-118
FITEM,5,211
FITEM,5,-3304
FITEM,5,4099
FITEM,5,-4162
FITEM,5,4167
FITEM,5,4180
FITEM,5,4191
FITEM,5,-4932
CM,_Y,NODE
NSEL,R, , ,P51X
CM,_Y1,NODE
CMSEL,S,_Y
CMDELE,_Y
D,_Y1,VX,0
D,_Y1,VY,0
D,_Y1,VZ,0
!*
CMDELE,_Y1
!*
ALLSEL,ALL
FLDATA12,PROP,DENS,0
FLDATA13,VARY,DENS,0
FLDATA12,PROP,VISC,21
FLDATA13,VARY,VISC,1
FLDATA12,PROP,COND,0
FLDATA13,VARY,COND,0
FLDATA12,PROP,SPHT,0
FLDATA13,VARY,SPHT,0
!*
FLDATA7,PROT,DENS,CONSTANT
FLDATA8,NOMI,DENS,1400,
FLDATA9,COF1,DENS,0
FLDATA10,COF2,DENS,0
FLDATA11,COF3,DENS,0
FLDATA7,PROT,VISC,POWL
FLDATA8,NOMI,VISC,40210,
FLDATA9,COF1,VISC,10,
FLDATA10,COF2,VISC,1,
FLDATA11,COF3,VISC,0.3,
FLDATA12,PROP,IVIS,40210,
FLDATA7,PROT,COND,CONSTANT
FLDATA8,NOMI,COND,-1,
FLDATA9,COF1,COND,0
FLDATA10,COF2,COND,0
FLDATA11,COF3,COND,0
FLDATA7,PROT,SPHT,CONSTANT
FLDATA8,NOMI,SPHT,-1,
FLDATA9,COF1,SPHT,0
FLDATA10,COF2,SPHT,0
FLDATA11,COF3,SPHT,0
FLDATA2,ITER,EXEC,100,
FLDATA2,ITER,OVER,0,
FLDATA2,ITER,APPE,0,
FLDATA3,TERM,VX,0.01,
FLDATA3,TERM,VY,0.01,
FLDATA3,TERM,VZ,0.01,
FLDATA3,TERM,PRES,1e-008,
FLDATA3,TERM,TEMP,1e-008,
FLDATA3,TERM,ENKE,0.01,
FLDATA3,TERM,ENDS,0.01,
FLDATA5,OUTP,SUMF,10,
SOLVE
FINISH
/POST1
SET,LAST
/post1
asel,s,area,,1,,,1
nsla,s,1
*get,nmin,node,,num,min
*get,nmax,node,,num,max
*SET,j,0
*SET,nvz,0
*dim,ennum,,nmax,1
*do,i,nmin,nmax,1
*if,nsel(i),gt,0,then
*SET,j,j+1
*get,ennum(j,1),node,i,v,z
*SET,nvz,nvz+ennum(j,1)
*endif
*enddo
asel,s,area,,17,,,1
nsla,s,1
*get,nmin,node,,num,min
*get,nmax,node,,num,max
*SET,x,0
*SET,nvz2,0
*dim,enum,,nmax,1
*do,i,nmin,nmax,1
*if,nsel(i),gt,0,then
*SET,x,x+1
*get,enum(x,1),node,i,v,z
*SET,nvz,nvz+enum(x,1)
*endif
*enddo
*SET,anvz,(nvz+nvz2)/(j+x)
*SET,Vobj,50-anvz |