主要的调试内容包括如下几个方面:
1.材料参数
2.质量阻尼
3.采用动态松弛施加初始静压
4.采用*LOAD_DENSITY_DEPTH施加初始静压
模型介绍:(采用了静水压ppt的模型参数
在此基础上修订了材料参数)
材料参数:
水和空气采用的的相对大气压(通过调整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
$ 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/基频周期)
压力云图
压力时程