本帖最后由 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
|