求助：充满液体的液罐流固耦合分析

wjiangq

还是这样贴出来比较直观，高手帮忙看看有什么问题阿，比较急，多谢了
*KEYWORD
*TITLE
                                                                                
$
*DATABASE_FORMAT
         0
$
*NODE
       1-3.250459894E-01 8.000000000E-01 6.849626586E-01       0       0
       2 6.732808012E-02 8.000000000E-01 1.019890375E+00       0       0
       3-2.524707700E-01 8.000000000E-01 9.674021270E-01       0       0
                       ...............................
    2866-2.812430068E-01 8.000000000E-01 1.058341047E+00       0       0
    2867-3.191978444E-01 8.000000000E-01 1.015554003E+00       0       0
$
*SECTION_SOLID_ALE
         1         6
         
*SECTION_SHELL
         2         2    1.0000       2.0       0.0       0.0         0
0.600E-02 0.600E-02 0.600E-02 0.600E-02  0.00   
$
*CONTROL_ALE
3,1,2,-1

*CONSTRAINED_LAGRANGE_IN_SOLID
         2         1         1         1         3         4         2         1
         0         2       0.1                             

$
*ALE_REFERENCE_SYSTEM_GROUP
1,1,4,0,0,0,0,
0,0,0,1.3
*MAT_NULL
         1  998.     -10.0     0.868E-03  0.00      0.00      0.00      0.00   
*EOS_GRUNEISEN
         1  0.00      0.00      0.00      0.00      0.00      0.00      0.00   
  0.00   
*MAT_ELASTIC
         2 0.783E+04 0.205E+12  0.300000       0.0       0.0       0.0
$
*PART
Part          1 for Mat         1 and Elem Type         1
         1         1         1         1         0         0         0
$
*PART
Part          2 for Mat         2 and Elem Type         2
         2         2         2         0         0         0         0
*ELEMENT_SOLID
       1       1     452     451     238     237    1052    1051     838     837
       2       1    1052    1051     838     837    1288    1287    1074    1073
                            ................................
    1565       1    1665     630     629    1664    1901     605     604    1900
    1566       1    1901     605     604    1900     116     782     781     117
*ELEMENT_SHELL
    1567       2    2017    2019    2081    2076
                ....................   
    2414       2    2750    2238    2239    2764
    2415       2    2749    2237    2238    2750
*SET_SEGMENT
         1     0.000     0.000     0.000     0.000
      2017      2019      2081      2076     0.000     0.000     0.000     0.000
                         ..........................
      2750      2238      2239      2764     0.000     0.000     0.000     0.000
      2749      2237      2238      2750     0.000     0.000     0.000     0.000
*DEFINE_CURVE
         1         0     1.000     1.000     0.000     0.000
  0.000000000000E+00  0.000000000000E+00
  1.000000000000E-01  1.960000000000E+01
*SET_NODE_LIST
         1     0.000     0.000     0.000     0.000
      2017      2018      2019      2020      2021      2022      2023      2024
      2025      2026      2027      2028      2029      2030      2031      2032
                          ..............................
      2857      2858      2859      2860      2861      2862      2863      2864
      2865      2866      2867
*BOUNDARY_PRESCRIBED_MOTION_SET
         1         2         1         1     1.000         0 0.000     0.000   
*SET_NODE_LIST
         2     0.000     0.000     0.000     0.000
      2041      2073      2191      2192      2193      2194      2195      2228
      2233      2266      2267      2268      2269      2270
*BOUNDARY_SPC_SET
         2         0         1         0         1         0         0         0
*CONTROL_ENERGY
         2         2         2         2
*CONTROL_SHELL
  20.0             1        -1         1         2         2         1
*CONTROL_BULK_VISCOSITY
  1.50     0.600E-01
*CONTROL_TIMESTEP
    0.0000    0.9000         0  0.00      0.00   
*CONTROL_TERMINATION
   10.             0   0.00000   0.00000   0.00000
*DATABASE_BINARY_D3PLOT
100.0   
*DATABASE_BINARY_D3THDT
10.00   
*DATABASE_EXTENT_BINARY
         0         0         3         1         0         0         0         0
         0         0         4         0         0         0
*SET_NODE_LIST
         3     0.000     0.000     0.000     0.000
         1         2         3         4         5         6         7         8
         9        10        11        12        13        14        15        16
                           ...........................
      2001      2002      2003      2004      2005      2006      2007      2008
      2009      2010      2011      2012      2013      2014      2015      2016
*INITIAL_VELOCITY
         3         0         0
  0.00      16.7      0.00      0.00      0.00      0.00   
*SET_NODE_LIST
         4     0.000     0.000     0.000     0.000
      2017      2018      2019      2020      2021      2022      2023      2024
   
   
      2857      2858      2859      2860      2861      2862      2863      2864
      2865      2866      2867
*INITIAL_VELOCITY
         4         0         0
  0.00      16.7      0.00      0.00      0.00      0.00   
*END

wjiangq

果然没人回答
不过我自己已经搞出来了

johnie2004

有个疑问：

*EOS_GRUNEISEN
1 0.00 0.00 0.00 0.00 0.00 0.00 0.00
0.00

1后面的0.00恐怕不是吧，应该是100吧！
如果高速流体，应该是要1500。

如果说错了，请大家指教！

wjiangq

谢谢你的提醒阿，不过我将数据更改之后计算发现，结果似乎好像没什么变化。
我再看看到底是怎么回事，关键这几个参数是和水在冲击波作用下的某些性质相联系的，我手头又没有这方面的资料。

johnie2004

流体的状态方程最难对付了，一直没有时间看相关的资料！（每次都用数字试试，缺乏严谨性）一起学习进步！以前有大侠提到要看相关的书籍！哎，太懒！

wjiangq

我没有看到相关的贴子阿
如果你还记得，在这里回复一下好么。

johnie2004

看看这个帖子

http://www.simwe.com/forum/viewthread.php?tid=313094

wjiangq

多谢了