（求助）如何处理多场耦合求解器MFS算出来得流固耦合结果？

zhoulzu

我用MFS求解器算了一个流固耦合问题, 问题描述见"Load Transfer Coupled Physics Analysis"这一章的算例部分。

在ansys帮助里，这个算例是用“基于PHYSICS Enviroment”的方法求解耦合问题的，我现在重新用MFS求解器对这个问题进行了求解。

问题的基本描述是这样的：流体流过一个弹性障碍物（称为Gasket的东西, 见附件1的紫色部分），求稳态解

问题1：
粘性系数小的时候，例如1.0e-2时，问题不收敛；原因是，障碍物的自由端位移达到了1e6量级，显然不对了，流体速度则超过1e9量级，也不对；只有较大的粘性系数才能收敛；但是，ansys的原算例给的粘性系数很小，为4.6e-4，而且按照这个算例地原方法是能够收敛的
具体错误提示为：
  *** ERROR ***                           CP =      27.640   TIME= 14:29:48
The value of UY at node 898 is 111811807.  It is greater than the      
current limit of 1000000.  This generally indicates rigid body motion   
as a result of an unconstrained model.  Verify that your model is      
properly constrained.

问题2：后处理
画出来的流场压力、速度、流线contour图仍然是未变形的流场区域(见附件)，怎么才能画出障碍物变形之后的流场压力、速度、流线图？


下面是命令流：

/prep7
a=10e-2
b=30e-2
Vin=3.5 !m/s
!!!!!!!流体常数：
rho=1e3 !kg/m^3
nu=1e-1 !kg*s/m , 黏性
t=2e-2 ! Gasket 厚度
w=8e-2 !
h=14e-2 !
!!!!!!!!! Gasket弹性常数：
mu=0.3
E=2.0e11 ! Pa
!!!!!!!!!!!!! Create Geometry
k,1,0,0
k,2,a,0
k,3,a,b
k,4,0,b
k,5,a,h
k,6,a,h+t
k,7,a-w,h+t
k,8,a-w,h

!!!!!!!!Fluid Model
ET,1,Fluid141
allsel,all
flda,solu,flow,1
flda,advm,mome,supg ! 平流方程的离散方法,  
flda,iter,exec,100 !
flda,outp,sumf,10
FLDATA5,OUTP,TAUW,1
FLDATA34,MIR,MOME,0.1  ! 松弛系数，

fldata7,prot,dens,constant
fldata7,prot,visc,constant
fldata8,nomi,dens,rho   !kg/m3
fldata8,nomi,visc,nu !kg*s/m  
flda,conv,pres,1.E-9

!!!FLUID area
a,1,2,5,8,7,6,3,4
!!!!!meshing control
lsel,s,loc,y,h+t
lsel,a,loc,x,a-w
lsel,a,loc,y,h
lesize,all,t/10,,,1.0,,-1,-2
!!!!!!!AATT, MAT, REAL, TYPE, ESYS, SECN
aatt,1,,1
smrtsize,3
amesh,all
!!!B.C.
!!!!LSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP
lsla
lsel,s,loc,x,0
lsel,a,loc,x,a
dl,all,,vx,0.,1
dl,all,,vy,0.,1
lsel,s,loc,y,h !!!!!!Gasket壁
lsel,a,loc,y,h+t
lsel,a,loc,x,a-w
dl,all,,vx,0.,1
dl,all,,vy,0.,1
lsel,s,loc,y,b !!!!!!!出口压力:0
dl,all,,pres,0.,1
lsel,s,loc,y,0 !!!!!进口速度
dl,all,,vx,0,1
dl,all,,vy,Vin,1
allsel,all
!!!!!!!载荷传递设置
lsel,s,loc,y,h
sfl,all,fsin,1
lsel,s,loc,y,h+t
sfl,all,fsin,2
lsel,s,loc,x,a-w
sfl,all,fsin,3

!!!!!!!!!!!!!!!!!! Structure Model
ET,2,PLANE42
mp,EX,2,E
mp,nuxy,2,mu
!mp,prxy,2,0.2
!!!!!!!!Gasket area
!a,5,6,7,8
blc5,a-w/2,h+t/2,w,t
asel,s,area,,2
lsla
smrtsize,4
!!!!!!AATT, MAT, REAL, TYPE, ESYS, SECN
aatt,2,,2,,2
lsel,s,loc,y,h+t
lsel,a,loc,y,h
lesize,all,t/5!,,,1.0,-1,-2
lsel,s,loc,x,a-w
lesize,all,t/5!,,,1.0,-1,-2
lsel,s,loc,x,a
lesize,all,t/5!,,,1.0,-1,-2

amesh,all
!!!!!!!B.C.
lsla
lsel,r,loc,x,a
dl,all,,ux,0.,1
dl,all,,uy,0.,1
!!!!!!!!!!!载荷传递
lsla
lsel,r,loc,y,h
sfl,all,fsin,1
lsla
lsel,r,loc,y,h+t
sfl,all,fsin,2
lsla
lsel,r,loc,x,a-w
sfl,all,fsin,3
allsel,all
finish
/solu
!!!!!!!!!!!!!!!!!!!!!!!多场耦合设置
!!  1. 激活MFS
mfanalysis,on
!!  2. 物理场编号
mfelem,1,1
mfelem,2,2
!!  3. 定义求解顺序
mforder,1,2
!!  4. 收敛条件
mfconv,all,1e-5
!!   5. 设置每个物理场的文件
mffname,1,fluid
mffname,2,struct
!!   6. 为每一个物理场设置求解选项
antype,static !!Fluid Field
morph,on
mfcmmand,1,fluid

antype,static !!! Struct field
deltim,1
morph,off
kbc,1
mfcmmand,2,struct


!!!!!  9. 设置载荷传递内容
mfsurface,1,1,forc,2
mfsurface,1,2,disp,1
mfsurface,2,1,forc,2
mfsurface,2,2,disp,1
mfsurface,3,1,forc,2
mfsurface,3,2,disp,1
!!!!!! 10. 全局设置：如载荷传递方式
mfinter,cons ! globally conservative load transfer
!!!!! 11. Stagger setting
mfiter,20
mfrelax,all,0.5 !!!!!松弛因子
mfdtime,1

!!! 12. 时间历程设置：最外层循环
mftime,1
!!!! 13. 输出间隔设置
mfcalc,1,1
mfcalc,2,1
mfoutput,1

!!!! 14. 求解
solve
save
finish

[ 本帖最后由 zhoulzu 于 2008-2-20 16:56 编辑 ]

hulifofo

顶一个 呵呵

liangmengbuct

学习了:lol