流固耦合问题讨论

pengweicai

不要一次性要命令流 拷进去 运行。

建议分段 导入。比如 前处理的， 求解的， 后处理的。...

cdh840102

分段导入也不行...
我自己把整个命令走了一遍,还是同样的错误提示.实在不知道问题出在哪边,在几个论坛上面问了,没人回答.
pengweicai  不知道你有没有看过流固耦合中那圆环在水中固有频率的例子,我想能不能参考这个例子,仿真出泵膜一面在有限水中,一面在空气中的固有频率呢,另外,在这个例子,它只设置了外圆环与水接触,内圆环是否能理解成在空气中,能否把无限吸收边界去掉,换成有限吸收边界,或者干脆不设呢.由于是新手,异想天开一下,不知道能不能行.下面是圆环例子的命令流:
finish
/clear
/PREP7
!定义单元类型
ET,1,PLANE42 ! structural element
ET,2,FLUID29 ! acoustic fluid element with ux & uy
ET,3,129 ! acoustic infinite line element
r,3,0.31242,0,0
ET,4,FLUID29,,1,0 ! acoustic fluid element without ux & uy
!材料属性
MP,EX,1,2.068e11
MP,DENS,1,7929
MP,NUXY,1,0
MP,DENS,2,1030
MP,SONC,2,1460
! 创建四分之一模型
CYL4,0,0,0.254,0,0.26035,90
CYL4,0,0,0.26035,0,0.31242,90
! 选择属性,网格划分
ASEL,S,AREA,,1
AATT,1,1,1,0
LESIZE,1,,,16,1
LESIZE,3,,,16,1
LESIZE,2,,,1,1
LESIZE,4,,,1,1
MSHKEY,1
MSHAPE,0,2D ! mapped quad mesh
AMESH,1
ASEL,S,AREA,,2
AATT,2,1,2,0
LESIZE,5,,,16,1
LESIZE,7,,,16,1
LESIZE,6,,,5
LESIZE,8,,,5
MSHKEY,0
MSHAPE,0,2D ! mapped quad mesh
AMESH,2
! 关于Y轴镜像
nsym,x,1000,all ! offset node number by 1000
esym,,1000,all
! 关于y轴镜像
nsym,y,2000,all ! offset node number by 2000
esym,,2000,all
NUMMRG,ALL ! merge all quantities
esel,s,type,,1
nsle,s
esln,s,0
nsle,s
esel,inve
nsle,s
emodif,all,type,4
esel,all
nsel,all
! 指定无限吸收边界
csys,1
nsel,s,loc,x,0.31242
type,3
real,3
mat,2
esurf
esel,all
nsel,all
! 标识流固交接面
nsel,s,loc,x,0.26035
esel,s,type,,2
sf,all,fsi,1
nsel,all
esel,all
FINISH
/solu
antype,modal
modopt,damp,10
mxpand,10,,,yes
solve
finish
为了便于对比，也对圆环在空气中做了模态分析
finish
/clear
/PREP7
!定义单元类型
ET,1,PLANE42 ! structural element
!材料属性
MP,EX,1,2.068e11
MP,DENS,1,7929
MP,NUXY,1,0
! 创建四分之一模型
CYL4,0,0,0.254,0,0.26035,90
! 选择属性,网格划分
ASEL,S,AREA,,1
AATT,1,1,1,0
LESIZE,1,,,16,1
LESIZE,3,,,16,1
LESIZE,2,,,1,1
LESIZE,4,,,1,1
MSHKEY,1
MSHAPE,0,2D ! mapped quad mesh
AMESH,1
! 关于Y轴镜像
nsym,x,1000,all ! offset node number by 1000
esym,,1000,all
! 关于y轴镜像
nsym,y,2000,all ! offset node number by 2000
esym,,2000,all
NUMMRG,ALL
/solu
antype,modal
modopt,lanb,10
mxpand,10,,,yes
solve
finish