求助:ansys帮助中的流固偶合例子为什么会出错

cdh840102

我最近在学习FSI流固耦合分析方法，我把HELP里面的例子的命令流复制下来，用文本格式导入ansys中，可老是出错，但由于是新手，不知道错误在哪边，希望高手能指点一下。

This example problem considers the fluid-solid interaction of a simple piston with a dissimilar mesh interface. The piston displacements are prescribed by a sinusoidal distribution in time. Figure 3.2: "Simple Piston" shows the mesh and boundary conditions for the problem. You determine the displacement and velocity distributions on the fluid side at the fluid-solid interface.


3.6.2. Problem Specifications
The cylinder and piston dimensions are:

Cylinder length = 1.0
Cylinder diameter = 0.25
Piston height = 0.10

The solid material properties are:

Young's modulus = 1.0e11
Density = 10
Poisson's ratio = 0.3

The fluid properties are:

Bulk modulus = 10000
Density = 1.0
Viscosity = 0.001
Thermal conductivity = 1.0
Specific heat = 1.0

The piston motion is sinusoidal with a maximum displacement of 0.85

/BATCH
MAXDISP = 0.85
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNC_C1
*DEL,_FNC_C2
*SET,_FNCNAME,'uxdsp'
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*SET,_FNC_C1(1),MAXDISP
*SET,_FNC_C2(1),2.5
! /INPUT,time.func
*DIM,%_FNCNAME%,TABLE,6,4,1
! Begin of equation: a*sin(omega*{TIME})
%_FNCNAME%(0,0,1)= 0.0, -999
%_FNCNAME%(2,0,1)= 0.0
%_FNCNAME%(3,0,1)= %_FNC_C1(1)%
%_FNCNAME%(4,0,1)= %_FNC_C2(1)%
%_FNCNAME%(5,0,1)= 0.0
%_FNCNAME%(6,0,1)= 0.0
%_FNCNAME%(0,1,1)= 1.0, -1, 0, 1, 18, 3, 1
%_FNCNAME%(0,2,1)= 0.0, -1, 9, 1, -1, 0, 0
%_FNCNAME%(0,3,1)=   0, -2, 0, 1, 17, 3, -1
%_FNCNAME%(0,4,1)= 0.0, 99, 0, 1, -2, 0, 0
! End of equation: a*sin(omega*{TIME})
!-->
/PREP7  
ET,1,141
KEYOPT,1,4,1
ET,2,42
RECT,0.1,1.0,,0.25  
TYPE,1
LESIZE,4,,,4
LESIZE,2,,,4
LESIZE,1,,,8
LESIZE,3,,,8
ALLSEL  
AMESH,ALL
ALLSEL,BELOW,AREA   
NSEL,R,LOC,X,0.1
SF,ALL,FSIN,1   
ALLSEL,BELOW,AREA  
NSEL,R,LOC,X,1  
D,ALL,UX,0  
D, ALL, UY,0
D, ALL, PRES,0
ALLSEL, BELOW, AREA  
NSEL,R,LOC,Y,0  
D, ALL, UY,0  
D, ALL, VY,0  
ALLSEL, BELOW, AREA  
NSEL,R,LOC,Y,0.25   
D,ALL,UY,0  
D,ALL,VY,0  
ALLSEL  
RECT,,0.1,,0.25
LESIZE,5,,,2
LESIZE,7,,,2
LESIZE,8,,,4
LESIZE,6,,,4
TYPE,2  
MAT,1
ALLSEL  
ASEL,S,,,2  
ALLSEL,BELOW,AREA   
AMESH,ALL   
ALLSEL,BELOW,AREA   
NSEL,R,LOC,X,0.1
SF,ALL,FSIN,1   
ALLSEL,BELOW,AREA   
NSEL,R,EXT  
D,ALL,UX,%uxdsp%
D,ALL,UY,0  
ALLSEL
! Solid Material properties
MP,EX,1,2.0E11
MP,DENS,1,7800
MP,PRXY,1,0.3
FLDATA1,SOLU,FLOW,1
FLDATA1,SOLU,TEMP,0
FLDATA1,SOLU,TRAN,1
FLDATA1,SOLU,ALE,1
FLDATA2,TIME,GLOB,10,
! Fluid Material properties
FLDATA7,PROT,DENS,AIR-SI
FLDATA7,PROT,VISC,AIR-SI
FINISH
/SOLU
FSAN,ON                            !FSI analysis on
FSCO,ALL,1.0e-2               !Convergence for all FSI interface quantities
FSRE,ALL,1.0                 !Relaxation for all FSI interface quantities
FSOU,1                                        !Output frequency for results file
DELTIM,0.0125               ! Time increment for SOLID analysis
FLDA, TIME, STEP, 0.0125         !Time increment for FLUID analysis
FSIT, 10                                   ! Number of stagger iterations per time step
FSTI, 1.0                            ! End time for FSI analysis
FSDT,0.0125                   ! Time increment for FSI analysis
FSOR,SOLID                           ! Order for FSI analysis (SOLID first)
FSTR,SOLID,                   ! Transient analysis for SOLID analysis
FSTR,FLUID,            ! Transient analysis for FLUID analysis
SAVE
SOLVE
SAVE
FINISH
/SOLU
FSRS,TIME,1.0                        ! Restart time for FSI analysis
FSTI,1.25                            ! End time for FSI analysis
FSDT,0.0125                          ! Time increment for FSI analysis
SOLVE
SAVE
FINISH
/POST26
FILE,,rfl
NSOL,2,23,UX
NSOL,3,23,VX
PRVAR,2
PLVAR,2
PRVAR,3
PLVAR,3
FINISH

cdh840102

不知道怎么把图片弄上来.就把错误消息打出来一下.
error attempting tp back up unUseable results file. check read/write permission . FILE.RFO.说结果文件不可用,但是不知道问题出在哪边.

mututu_gui

我一直读下来，从/batch到finish,没有出错，在Main Menu>TimeHist Postpro>Settings>File里面读一下以你相应的RFL文件，再进行
NSOL,2,23,UX
NSOL,3,23,VX
PRVAR,2
PLVAR,2
PRVAR,3
PLVAR,3
FINISH
就没问题了。

tanzouqing

楼主正解,可原因还是不清楚?

uqhxie

i ve been suffering on the same problem

uqhxie

现在知道了 要新建一个不同名的文件