流固耦合模态分析的典型事例

xiaowujun111

请教欧阳大师，我做了一个模拟重力坝的流固耦合例子，算出的基频怎么跟没有水一样呢！？我弄不明白，请大师指点。命令流如下：
/prep7
et,1,plane42
et,2,solid65
et,3,fluid30
mp,ex,1,3.0e10
mp,nuxy,1,0.167
mp,dense,1,2400
mp,dense,2,1000
mp,sonc,2,1517
mp,mu,2,0
rectng,0,80,0,120
rectng,-240,0,0,110
lsel,s,,,1,3,2
lesize,all,,,8
lsel,s,,,2,4,2
lesize,all,,,12
amesh,1
lsel,s,,,5,7,2
lesize,all,,,24
lsel,s,,,6,8,2
lesize,all,,,11
amesh,2
extopt,esize,5,0,
extopt,aclear,1
type,2,
mat,1
vext,1,,,,,50
eplot
extopt,esize,5,0,
extopt,aclear,1
type,3,
mat,2
vext,2,,,,,50
eplot
vglue,all
nummrg,all
numcmp,all
allsel
/solu
antype,modal
modopt,unsym,10
eqslv,front
mxpand,10,,,yes,
outp,all
outres,all,all
/output,zhijieouhe1,txt
nsel,s,loc,z,0
nsel,a,loc,z,50
d,all,uz
nsel,s,loc,y,0
nsel,r,loc,x,0,80
d,all,all
nsel,s,loc,y,0
nsel,r,loc,x,-240,0
d,all,uy
nsel,s,loc,x,0
nsel,r,loc,y,0,110
nplot
sf,all,fsi
allsel
eplot

da,11,pres,0
sfgrad,pres,,y,,-9810
da,12,pres,9810*120
!nsel,s,loc,x,-240
!d,all,all


allsel
eplot
solv
fini

/post1
set,next
plnsol,u,sum,2,1


先谢过欧阳大师。

xiaowujun111

无水时：
SET   TIME/FREQ    LOAD STEP   SUBSTEP  CUMULATIVE
     1  2.5473             1         1         1
     2  7.4958             1         2         2
     3  8.5593             1         3         3
     4  17.150             1         4         4
     5  20.112             1         5         5
     6  20.220             1         6         6
     7  20.235             1         7         7
     8  21.140             1         8         8
     9  21.348             1         9         9
    10  23.148             1        10        10

直接耦合
   SET   TIME/FREQ    LOAD STEP   SUBSTEP  CUMULATIVE
     1  2.5473             1         1         1
     2  0.0000             1         1         1
     3  4.1892             1         2         2
     4  0.0000             1         2         2
     5  6.2158             1         3         3
     6  0.0000             1         3         3
     7  7.4958             1         4         4
     8  0.0000             1         4         4
     9  8.5593             1         5         5
    10  0.0000             1         5         5
    11  8.8748             1         6         6
    12  0.0000             1         6         6
    13  10.697             1         7         7
    14  0.0000             1         7         7
    15  11.541             1         8         8
    16  0.0000             1         8         8
    17  11.991             1         9         9
    18  0.0000             1         9         9
    19  13.307             1        10        10
    20  0.0000             1        10        10

欧阳中华

.
    不好意思，我做得仅仅是加筋板结构流固耦合振动问题，也是自己用FORTRAN编程，ANASYS和土木中的重力坝都不了解，希望有这方面经验的网友出来帮助讨论...

xiakeqian

你好，请问装水后计算出的频率为什么有20阶？模态提取不是设置为10阶吗，

1991

欧阳中华 发表于 2006-4-19 17:35
很高兴能看到这样的例题。问题是在：

1、单元的选择；

那请问一下这里为什么要定义声速呢？

mxlzhenzhu

1991 发表于 2014-5-10 16:48
那请问一下这里为什么要定义声速呢？

流体属性定义要有声速吧；

你参考Virtual Lab模块软件下的流固耦合问题，算起来很简单。

欧阳中华

1991 发表于 2014-5-10 16:48
那请问一下这里为什么要定义声速呢？

。
   流体方程建立时需要涉及声在流体里的速度.... ...

欧阳中华

1991 发表于 2014-5-10 16:48
那请问一下这里为什么要定义声速呢？

。
   流体方程建立时需要涉及声在流体里的速度.... ...

1991

欧阳中华 发表于 2014-5-13 17:53
。
   流体方程建立时需要涉及声在流体里的速度.... ...

谢谢！明白了，那会儿刚刚开始理论知识不足，谢谢您的解答！希望以后碰到类似问题还能向您请教！

1991

mxlzhenzhu 发表于 2014-5-10 18:22
流体属性定义要有声速吧；

你参考Virtual Lab模块软件下的流固耦合问题，算起来很简单。

嗯是的，这个问题现在看有点白痴啊！谢谢您的热心解答！谢谢！

1991

xiakeqian 发表于 2012-10-31 10:53
你好，请问装水后计算出的频率为什么有20阶？模态提取不是设置为10阶吗，

解答：选用的特征值求解器不是常规的求解结构的求解器，特征值表示有差别，第一个为一阶频率的虚部，第二个为一阶的实部，……