关于ALE分析老出错，提示没找到load set~~
本帖最后由 wyc412721 于 2011-7-23 22:14 编辑关于ALE Analysis分析出错
大家好：
我在学FLUID中的ALE分析，在运行帮助中的例子
7.7. ALE Analysis of a Moving Cylinder时，程序见下面，网格划分没问题，自己再原例子上加了个后处理：
/post1
set,4,2
/dscal,,1.0
pldi,2
finish
但是在后处理/post1时，总是提示我：”load set not found on result file”,请大家帮我看看，多谢啦！！！！
附源程序：
FINISH
/CLEAR
/batch,list
/verify,dv-1101c
/title, dv-1101c, Flotran Remeshing: Cylinder pass through a channel
/COM,****************************************************
/COM,* Cylinder pass through a channel (pure fluid flow)
/COM,* Verify partial domain re-meshing
/COM,****************************************************
r1 = 2.5
l1 = 5.0
l2= 20.0
l3=60.0
h = 10.0
rx0= l1+(l2-l1)/2.0
ry0 = h/2.0
nl1 = 5
nl3 = 10
nr=8
nc=8
/PREP7
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNC_C1
*SET,_FNCNAME,'DIS'
*DIM,_FNC_C1,,1
*SET,_FNC_C1(1),1.0
! /INPUT,force.func
*DIM,%_FNCNAME%,TABLE,6,3,1
!
! Begin of equation: a*{TIME}
*SET,%_FNCNAME%(0,0,1), 0.0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 1, 17, 3, 1
*SET,%_FNCNAME%(0,2,1), 0.0, 99, 0, 1, -1, 0, 0
*SET,%_FNCNAME%(0,3,1),
0
! End of equation: a*{TIME}
RECTNG,0,l1,0,h,
RECTNG,l1,l2,0,h,
RECTNG,l2,l3,0,h,
aglue,all
CYl4,rx0,ry0,0,45,r1,135
CYl4,rx0,ry0,0,135,r1,225
CYl4,rx0,ry0,0,225,r1,315
CYl4,rx0,ry0,0,315,r1,405
asel,s,,,2
asel,a,,,3
asel,a,,,6
asel,a,,,7
allsel,below,area
aglue,all
allsel
FLST,3,4,5,ORDE,3
!Specifies data required for a picking operation (GUI).
FITEM,3,2
!Identifies items chosen by a picking operation (GUI).
FITEM,3,8
!Subtracts areas from areas.
FITEM,3,-10
ASBA,
4,P51X
aplot
et,1,141
KEYOPT,1,4,1
type,1
!Sets the element type attribute pointer.
mat,1
asel,s,,,1
esize,.5
mshape,1,2d
mshkey,0
amesh,all
et,2,141
KEYOPT,2,4,1
type,2
mat,1
asel,s,,,3
esize,.5
mshape,1,2d
mshkey,0
amesh,all
et,3,141
KEYOPT,3,4,1
type,3
mat,1
asel,s,,,5
esize,2.0
mshape,1,2d
mshkey,0
amesh,all
nsel,s,loc,x,0
D,ALL,VX,0.0
D,ALL,VY, 0.0
D,ALL,UX,0.0
D,ALL,UY, 0.0
nsel,s,loc,x,l3
D,ALL,pres,0.0
D,ALL,UX,0.0
D,ALL,UY, 0.0
nsel,s,loc,y,0
D,ALL,vy,0.0
D,ALL,UX,0.0
D,ALL,UY, 0.0
nsel,s,loc,y,h
D,ALL,vy,0.0
D,ALL,UX,0.0
D,ALL,UY,0.0
lsel,s,,,5
lsel,a,,,23
lsel,a,,,25
lsel,a,,,27
nsll,s,1
!select the nodes on lines,margins of the round hole
!nplot
d,all,UX,%DIS%
d,all,UY,0.0
d,all,VX,1.0
d,all,VY,0.0
allsel
cdwrite,db,fluid,cdb,
!Writes geometry and load database items to a file.
fini
/solu
FLDATA30,QUAD,MOMD,2,
!Controls the quadrature orders.控制面积积分的阶次，亦即积分点的个数
FLDATA30,QUAD,MOMS,2,
FLDATA30,QUAD,PRSD,2,
FLDATA30,QUAD,PRES,2,
FLDATA30,QUAD,THRD,0,
FLDATA30,QUAD,THRS,0,
FLDATA30,QUAD,TRBD,0,
FLDATA30,QUAD,TRBS,2,
/solu
FLDATA4,TIME,NUMB,100000
FLDATA4,TIME,TEND,44.0
FLDA,SOLU,ALE,T
! ALE solution
Controls which features of the solution algorithm are activated.
FLDATA1,SOLU,FLOW,1
FLDATA1,SOLU,TRAN,1
!FLDATA1,SOLU,TURB,1
!Activates the turbulence model
FLDATA2,TIME,GLOB,5
!Sets iteration and output controls for steady state analyses.
FLDATA7,PROT,DENS,CONSTANT
!Specifies the type of fluid property.
FLDATA8,NOMI,DENS,1.0
!Specifies the NOMI coefficient of the fluid property equation.
FLDATA7,PROT,VISC,CONSTANT
FLDATA8,NOMI,VISC,1.0
!FLDA,BULK,BETA,1.0e5
FLDA,TIME,STEP, 2.0
FLDATA,OUTP,TAUW,T
!?
FLDATA,TIME,METH,NEWM
FLDATA,TIME,DELT,0.5
!
!! Set ANSYS-STRUCTURE commands
SAVE
/COM
/COM
Re-meshing Commands
/COM
FLDATA,REMESH,ELEM,ALL
! all defined element re-meshing
FLDATA,REMESH,ARMA,10.0
! maximum aspect ratio
FLDATA,REMESH,VOCH,5.0
! maximum volume change
FLDATA,REMESH,ARCH,5.0
! maximum aspect ratio change
SOLVE
SAVE
finish
/post1
!set,STEP,SUBSTEP
! sets the set number for processing
set,4,2
/dscal,,1.0
! sets the scaling for display
pldi,2
! display the deformed mesh overlayed with undeformed
Displays the displaced structure.
finish
! edges
代码编辑 BY wyc412721
页:
[1]