静水压力调试经验贴

蓝牙

有段时间比较关心水的初始静水压力分布，搜索了本版，也找到了一些线索，如下面的帖子

http://forum.simwe.com/forum.php?mod=viewthread&tid=979281&highlight=%E9%9D%99%E6%B0%B4%E5%8E%8B%E5%8A%9B
http://forum.simwe.com/forum.php?mod=viewthread&tid=885262&highlight=%2BLOAD%5C_DENSITY%5C_DEPTH
http://forum.simwe.com/forum.php?mod=viewthread&tid=737072&highlight=%2BLOAD%5C_DENSITY%5C_DEPTH
http://forum.simwe.com/forum.php?mod=viewthread&tid=846470&highlight=%E9%9D%99%E6%B0%B4%E5%8E%8B

现在把我关于静水压分布的调试结果分享给大家
主要的调试内容包括如下几个方面：
1.材料参数
2.质量阻尼
3.采用动态松弛施加初始静压
4.采用*LOAD_DENSITY_DEPTH施加初始静压

模型介绍：（采用了cnetuwmj提供的静水压ppt的模型参数
，再此向cnetuwm表示感谢，在此基础上修订了材料参数，）

材料参数：
    水和空气采用的的相对大气压（通过调整c0 和E0可以进行调整，本人参考了很多资料，觉得这套参数可行并且有效）
*MAT_NULL
$      MID        RO        PC        MU     TEROD     CEROD    YMBEAM    PRBEAM
         2 0.998E+03  -10.00e3  8.684E-4      0.00      0.00      0.00      0.00   
*EOS_LINEAR_POLYNOMIAL
$    EOSID        C0        C1        C2        C3        C4        C5        C6
         2       0E5   2.25E+9      0.00     0.000      0.00      0.00      0.00
$       E0        V0  
      0.00       0.0
*MAT_NULL
$      MID        RO        PC        MU     TEROD     CEROD    YMBEAM    PRBEAM
         1    1.1845     -10.0 1.8444E-5       0.0       0.0       0.0       0.0         
*EOS_LINEAR_POLYNOMIAL
$    EOSID        C0        C1        C2        C3        C4        C5        C6
         1      -1E5      0.00      0.00      0.00      0.40      0.40      0.00   
$       E0        V0  
     2.5E5         1   
$
调试结果：
1.重力+无初始化+无质量阻尼

                                压力云图

                        压力时程曲线
波形震荡非常厉害

2.重力+无初始化+质量阻尼（采用4*pi/基频周期）

质量阻尼采用*DAMPING_PART_MASS关键字

                        压力时程曲线
波形有明显改善，附k文件
3.重力+动力松弛静水压力初始化+质量阻尼（采用4*pi/基频周期）
*DEFINE_CURVE开启动态松弛，见k文件，底部压力正好是pgh

                     压力时程曲线
4.重力+LOAD_DENSITY_DEPTH+质量阻尼（采用4*pi/基频周期)

                              压力云图

                                压力时程
由云图和压力时程可以看出，底部压力为2.8E5Pa，刚好是pgh的两倍，说明在用*LOAD_DENSITY_DEPTH施加初始静水压力时，不能和重力加速度一起使用！
k文件

总结：
    1.如果要涉及重力加速度，采用*LOAD_DENSITY_DEPTH施加静水初压力是不对的！这个我觉得应该有什么设置可以控制的，但是现在还没有找到问题的解决方法，大家也帮帮忙！
    2.采用动力松弛施加重力加速度是可行的，但是也有一个缺陷，规模太大不容易收敛，在新版本中方法已经有了改进
    3.施加质量阻尼有利于消减波形震荡！


结贴

感谢长脚版主的加分鼓励，呵呵

tudouqd

*LOAD_DENSITY_DEPTH中的CR项是不是应该是-9.81？

songerking

我看斑竹和那位美女都是对水施加的静压力，实际上静水压力就是个球张量，除了水很多地方都存在这种无偏应力等梯度的环境，比如岩层深处，不知道这种办法对岩体是否具有相同功效，敬盼佳音。
另：计算静态结果，通过initial_stress_solid相关的命令实现的初始化静压是否效果一样。

蓝牙

songerking 发表于 2012-3-24 21:51
我看斑竹和那位美女都是对水施加的静压力，实际上静水压力就是个球张量，除了水很多地方都存在这种无偏应力 ...

岩石不可以只考虑静水压力，偏应力不可忽略，*LOAD_DENSITY_DEPTH应该不能用于你的问题
动力松弛应该可以，是要你材料考虑了偏应力，动态松弛是不管是水还是岩石的。
动力松弛实际就是求解重力载荷下的静平衡问题，实际上也可以采用ansys的隐式功能，求解出初始的应力场，代入到dyna中。

wephy51

下载学习了，好东西。

xth07

请教楼主，k文件里哪段文字表示“采用4*pi/基频周期”？ 谢谢。

xth07

楼主的文件float11.k在我的机子上运行了两个小时，动态松弛convergence factor徘徊在0.1左右，甚至有上升的趋势，无法收敛。请问楼主运行时，动态松弛收敛用了多长时间？

蓝牙

xth07 发表于 2013-2-27 11:04
楼主的文件float11.k在我的机子上运行了两个小时，动态松弛convergence factor徘徊在0.1左右，甚至有上升的 ...

动力松弛所用的时间很多，采用R5.5以上版本是有改进方法的

你参考21-hydroP （我1#的第一个附件）那个ppt中的后继方法可以极大改善收敛速度

simulin_k

cnetuwmj提供的ppt中说理想气体的C0不要定义，请教版主是参考什么定义的C0=1E+5?

simulin_k

手册上说这种做法好像是错误的

风速

请问楼主用过你提供的ppt中*INITIAL_ALE_HYDROSTATIC和*ALE_AMBIENT_HYDROSTATIC关键字联合设定初始静压力吗？我尝试用了一下，做出来的结果明显不对，最终下方的水流到上方的空气域中了

wjwme

不错不错

90357

真的是好东西，谢谢老师

90357

tudouqd 发表于 2013-11-25 08:16
*LOAD_DENSITY_DEPTH中的CR项是不是应该是-9.81？

朋友，这个重力加速的问题你理解透没？到底是正还是负啊

90357

风速 发表于 2014-8-20 21:10
请问楼主用过你提供的ppt中*INITIAL_ALE_HYDROSTATIC和*ALE_AMBIENT_HYDROSTATIC关键字联合设定初始静压力 ...

朋友，我最近也是卡在这儿了。如果我的流体边界定义无反射边界条件，则流体会流出模型；但是用*ALE_AMBIENT_HYDROSTATIC来模拟远场边界的话，水四周的part该怎么整？（不都是水材料的哇，那怎么可能划分为不同的part），请朋友指点迷津，谢谢啦！

qysim

几年之后再看这个帖子依然很有意义，想问一下楼主第一个附件的pdf是从哪个dyna的网站上找到的，我再官网并没有找到，想下载几篇这种官方教程学习。

90357

qysim 发表于 2016-5-22 15:08
几年之后再看这个帖子依然很有意义，想问一下楼主第一个附件的pdf是从哪个dyna的网站上找到的，我再官网并 ...

朋友，斑竹PDF里面的*initial_hydrostatic_ale和*ale_ambient_hydrostatic理解透了？（*aah里面的part即流体周围的环境边界的part怎么生产啊），请朋友指点一下

韦小勺

按楼主和ppt中的方法，发现只有在四周加载‘pool like’边界条件时压力才会稳定，假如设置了这种边界条件那么边界条件不能成为无反射边界条件。现在的问题是怎么在设置无反射边界条件的同时设置静水压力。谢谢大家。

振动狂魔

同问楼上的问题，应该可以探究下其他形式的加载方式。在模拟自由面气泡时需要设置边界条件为无反射，但是在楼主这种加载方式中，本人还未找到合理的加压方式。

41439937

谢谢分享