cdh840102 发表于 2007-12-9 21:27

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

我最近在学习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 发表于 2007-12-10 20:21

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

mututu_gui 发表于 2007-12-12 09:58

回复 #2 cdh840102 的帖子

我一直读下来,从/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 发表于 2011-9-7 21:57

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

uqhxie 发表于 2011-11-23 12:52

i ve been suffering on the same problem

uqhxie 发表于 2011-11-25 16:01

现在知道了 要新建一个不同名的文件
页: [1]
查看完整版本: 求助:ansys帮助中的流固偶合例子为什么会出错