有个帖子oscillating plate with two way fluid structure interaction 写得很清楚
按步骤来 可以出来 我忘记帖子是从哪里来的 只有纸质的 照着做了一下 例子倒是看懂了 为什么要设置声速啊? 有没有做大坝-水的耦合作用的呢? 网上找到的
水坝空库
/BATCH
KEYW,PR_SET,1
KEYW,PR_STRUC,1
KEYW,PR_THERM,0
KEYW,PR_FLUID,0
KEYW,PR_ELMAG,0
KEYW,MAGNOD,0
KEYW,MAGEDG,0
KEYW,MAGHFE,0
KEYW,MAGELC,0
KEYW,PR_MULTI,0
KEYW,PR_CFD,0
/GO
!*
/prep7
!* define material proterties
mp,dens,1,2650. !mat 1 for dam
mp,ex,1,3.15e10
mp,prxy,1,.167
!* define element type
et,1,PLANE42,,,2
!* define geometry
lwater=618.
k, 1,0.0,0.0,0.0
k, 2,70.2,0.,0.
k, 3,0.,66.5,0.
k, 4,21.9875,66.5,0.
k, 5,0.,103.,0.
k, 6,14.8,103.,0.
a, 1, 2, 4, 3
a, 3, 4, 6, 5
asel,s,loc,x,0.,100.
aatt,1,,1
cm,adam,area
allsel,all
!* mesh geometry
ESIZE,0,10
lsel,s,loc,y,0.1,66. !坝下部剖分分数
lesize,all,,,15
mshape,0,2D
mshkey,1
allsel,all
amesh,all
finish
/solu
antype,modal
MODOPT,LANB,30
MXPAND,30, , ,0
esel,s,mat,,1 !坝体约束
nsle,s
nsel,r,loc,y,-1.0,1.0
d,all,ux,0.
d,all,uy,0.
/pbc,all,,1
/pnum,type,1
/number,1
gplot
allsel,all
save
solve 满库
/BATCH
KEYW,PR_SET,1
KEYW,PR_STRUC,1
KEYW,PR_THERM,0
KEYW,PR_FLUID,0
KEYW,PR_ELMAG,0
KEYW,MAGNOD,0
KEYW,MAGEDG,0
KEYW,MAGHFE,0
KEYW,MAGELC,0
KEYW,PR_MULTI,0
KEYW,PR_CFD,0
/GO
!*
/prep7
!* define material proterties
mp,dens,1,2650. !mat 1 for dam
mp,ex,1,3.15e10
mp,prxy,1,.167
mp,dens,2,1000. !mat 2 for water
mp,sonc,2,1440
!* define element type
et,1,PLANE42,,,2
et,2,29
et,3,29,,1
!* define geometry
lwater=618.
k, 1,0.0,0.0,0.0
k, 2,70.2,0.,0.
k, 3,0.,66.5,0.
k, 4,21.9875,66.5,0.
k, 5,0.,103.,0.
k, 6,14.8,103.,0.
k, 7,-1.*lwater,103.,0.
k, 8,-1.*lwater,66.5,0.
k, 9,-1.*lwater,0.,0.
a, 1, 2, 4, 3
a, 3, 4, 6, 5
a, 8, 3, 5, 7
a, 9,1, 3, 8
asel,s,loc,x,0.,100.
aatt,1,,1
cm,adam,area
asel,s,loc,x,-1*lwater,0.
aatt,2,,3
cm,awater,area
allsel,all
!* mesh geometry
ESIZE,0,10
lsel,s,loc,y,0.1,66. !坝下部剖分分数
lesize,all,,,15
lsel,s,loc,x,-0.1,-1*lwater-1. !水体长度方向剖分分数
lsel,r,loc,y,-0.1,67.
lesize,all,,,40,0.5
lsel,s,loc,x,-0.1,-1*lwater-1. !水体长度方向剖分分数
lsel,r,loc,y,100.,104.
lesize,all,,,40,2.0
mshape,0,2D
mshkey,1
allsel,all
amesh,all
!更改与水体接触的单元类型
esel,s,type,,3
nsle,s
nsel,r,loc,x,-1.,1.
esln,r
emodif,all,type,2
allsel,all
finish
/solu
antype,modal
MODOPT,UNSYM,30
MXPAND,30, , ,0
/pnum,type,1
/number,1
nsel,s,loc,x,-1.,1. !库水坝体FSI面
esln,s
esel,r,mat,,2
sf,all,fsi
esel,s,mat,,1 !坝体约束
nsle,s
nsel,r,loc,y,-1.0,1.0
d,all,ux,0.
d,all,uy,0.
esel,s,mat,,2 !水体表面约束
nsle,s
nsel,r,loc,y,102.,104.
d,all,pres,0.
/pbc,all,,1
/psf,fsi,,2
allsel,all
gplot
save
solve 这是自己做的一个作业
ansys做的,半浸没梁振动算例
以半浸没与水中的桥墩模态问题为背景,并假设:
1. 桥墩为实心等截面的实体,实际桥墩模型应该是空心壳体,截面尺寸也非常复杂,因而需要分块划分单元。
2.不考虑桥墩中的钢筋,简化为由单一的混凝土构成桥墩,密度1.265e4 kg/m3,弹性模量3.15e10 N/m2,泊松比0.16
3.不考虑水面重力影响,忽略自由表面波,水的密度1e3 kg/m3,音速 1.44e3 m/s
4.坝体和水的交界面上,固体节点和液体节点位移耦合在一起,在单元中液体和固体粘在一起。
附件在
http://forum.vibunion.com/forum/thread-51033-1-1.html
请教渡槽水体与渡槽的边界条件怎么施加?
我也正在做渡槽流固耦合方面东西,请教一下流固耦合边界、滑移壁面边界和自由表面边界具体怎么用命令流在ANSYS中实现,请各位指教,不胜感激! 回复 怡人 的帖子昨晚我也遇到这个问题了。后来将shell63改为solid45就可以了。
我暂时知道的,声分析里面匹配单元是二维plane42+fluid29,三维的是solid45+fluid30;ansys+cfx分析里面,shell63+fluid142。
08的贴,10年回,呵呵。送给也像我这样在摸索的同学...... 我现在也在学习流固耦合的东西,谢谢各位的讨论。受益匪浅,学习了。 现在开始做一些流固耦合相关的东西,这个帖子的东西真的很好,好好学习~ 我碰到一个很奇怪的问题,我使用了FLUID130单元后,就只能采用Damped提取模态
错误提示为:
Whenusing the infinite acoustic elements (FLUID129,FLUID130),the only applicable modal analysis method is the Damped ansys workbench中怎么实现双向流固耦合,其中固体区域的橡胶的模型怎么设置? 先看看 。。。。