重力坝—库水—地基流固耦合问题

vankor

小弟毕设用ADINA做了一个重力坝—库水—地基的简单模型，用的是势流体，分析系统在水平地震作用下的动力响应。坝体材料参数为：弹性模量E=24GPa，泊松比 ，密度 ，坝体阻尼比为5%；地基材料参数为：弹性模量E=14.34GPa，泊松比 ，密度 ；库水体积模量参数为 ，密度 。对模型施加水平向地震加速度波。

模型采用传统的固定边界无质量地基模型
模型采用平面四节点四边形等参单元进行离散，每条边界分为20等份。离散后共剖分1891个结点，1814个单元，其中坝体单元214个，地基单元1600个，水体单元1200个。


在ADINA中，各个边界条件设置如下：地基两侧边及底边为固定边界，不考虑地基质量，但考虑地基刚度；库尾水边界设置为Fluid—Infinite Region，即无穷远尾水；水面不考虑表面重力波影响，设置为自由水面。整个系统重力荷载不考虑，只计算纯动力荷载。
模型中坝体阻尼采用瑞利比例阻尼，求得阻尼系数α=1.022 β=0.002， 。
在模型中考虑坝水相互作用，坝—基相互作用和水—基相互作用，忽略水体的可压缩性。计算时间步长为0.01秒，地震持续10秒，共1000步。模型解法采用Newmark隐式解法。
我是参照一个老师的博士论文，可是算出来的结果相差很大。总体来说我的结果偏大，比如坝踵动水压力（我提取的是FE—pressure）打概大到30%左右；坝踵处竖向应力大到三倍左右；而且坝面动水压力和加速度分布也不规律。我找了两个多周的问题也没找到原因，导师说可能是势流体与坝体流固耦合面以及势流体与库底耦合面设置有问题，我看了一些资料都是说自动生成的……
所以现在很纠结，那个大侠小弟救一命……跪谢！
我的in文件由于改了无数次，因此很长，但是模型很简单~~大家慢慢看，多谢！！
我的in文件在所发连接的的附件中“dam[url]http://www.simwe.com

vankor

in文件连接：[url]http://www.simwe.com

vankor

in文件在附件中

fengmo2194184

请问问题最终解决了吗？