寻求流固耦合处理的方法

初级学徒

我最近在学习清华版时党勇、李裕春和张胜民编著《基于ANSYS/LS-DYNA进行显示动力分析》，练习炸药在土壤中的爆炸作用的例子时，建模后修改K文件确认和书上K文件一样后提交计算，结果发现混凝土板在爆炸过程中无变化且无应力，做了几遍都是一样的。实在想不出哪里出错。望做过此例子的前辈给点提示

hmily861012

如果k文件确保无误，应该不会是这种结果，建议你将k文件上传，这样别人才能替你查找原因

初级学徒

建议不错，希望有人能帮忙

初级学徒

K文件在此，麻烦哪位前辈帮下忙

初级学徒

有点大上传一些K文件关键字：
*KEYWORD
*TITLE
                                                                                
$
*DATABASE_FORMAT
         0
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                               NODE DEFINITIONS                               $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*NODE
       1 6.000000000E+00 6.000000000E+00 0.000000000E+00       0       0
       2 6.000000000E+00 6.000000000E+00 1.500000000E+00       0       0
       3 6.000000000E+00-6.000000000E+00 0.000000000E+00       0       0
       4 6.000000000E+00-4.500000000E+00 0.000000000E+00       0       0
       5 6.000000000E+00-3.000000000E+00 0.000000000E+00       0       0
       6 6.000000000E+00-1.500000000E+00 0.000000000E+00       0       0
       7 6.000000000E+00 0.000000000E+00 0.000000000E+00       0       0
       8 6.000000000E+00 1.500000000E+00 0.000000000E+00       0       0
       9 6.000000000E+00 3.000000000E+00 0.000000000E+00       0       0
。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。


  91366 2.480000000E+02 8.400000000E+01 1.500000000E+00       0       0
   91367 2.480000000E+02 8.600000000E+01 1.500000000E+00       0       0
   91368 2.480000000E+02 8.800000000E+01 1.500000000E+00       0       0
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             SECTION DEFINITIONS                              $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*SECTION_SOLID_ALE
         1         11
         
*SECTION_SOLID_ALE
         2         11
         
*SECTION_SOLID_ALE
         3         11
         
*SECTION_SOLID
         4         0   
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                               ALE DEFINITIONS                                $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*ALE_MULTI-MATERIAL_GROUP
        1          1
        2          1
        3          1
$
*CONSTRAINED_LAGRANGE_IN_SOLID
1,2,0,0,0,5,3,0
0,0,0.15

*SET_PART_LIST
         1
         4
*SET_PART_LIST
         2
         2        3
$
*CONTROL_ALE
2,1,2,-1.00,0.00,0.00,0.00
0.00,0.00,0.00   
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             MATERIAL DEFINITIONS                             $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*MAT_HIGH_EXPLOSIVE_BURN
1,1.63,0.693,0.27,0.000e+00   
*EOS_JWL
1,3.71,7.430e-02,4.15,0.950,0.300,7.00E-02,1.00     
*INITIAL_DETONATION
1,0.00,0.00,0.00,0.00
*INITIAL_DETONATION
1,0.00,0.00,1.50,0.0
*MAT_SOIL_AND_FOAM
2,1.80,1.601E-04,1.328E+02,3.300E-03,1.310E-07,0.1232,0.0
0.0,0.0
0.0,0.050,0.090,0.11,0.15,0.19,0.21,0.22   
0.25,0.30
0.0,3.420E-02,4.530E-02,6.760E-02,1.270E-01,2.080E-01,2.710E-01,3.920E-01
5.660E-01,1.230E+00   
*MAT_NULL
         3 0.125E-02  0.00      0.00      0.00      0.00      0.00      0.00   
*EOS_GRUNEISEN
3,0.344,0.00,0.00,0.00,1.40,0.00,0.00   
0.00   
*MAT_PLASTIC_KINEMATIC
         4  2.65     0.4      0.300000 0.100E-02 0.400E-02 0.500   
  0.00      0.00     0.800
$
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              PARTS DEFINITIONS                               $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$
*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
$
*PART
Part          3 for Mat         3 and Elem Type         3
         3         3         3         3         0         0         0
$
*PART
Part          4 for Mat         4 and Elem Type         4
         4         4         4         0         0         0         0
。。。。。。。。。。。。。。。。。。。。。。。。。。。。
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              CONTACT DEFINITIONS                             $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              CONTROL OPTIONS                                 $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*CONTROL_CPU
  0.00   
*CONTROL_ENERGY
         2         1         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.6000         0  0.00      0.00   
*CONTROL_TERMINATION
3000.             0   0.00000   0.00000   0.00000
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                                 TIME HISTORY                                                                                $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*DATABASE_BINARY_D3PLOT
20.00   
*DATABASE_BINARY_D3THDT
4.000
$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              DATABASE OPTIONS                                                                         $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*DATABASE_EXTENT_BINARY
         0         0         3         1         0         0         0         0
         0         0         4         0         0         0
*END

hmily861012

如果k文件中主要关键字无误差，我猜测你的流固耦合没起作用，混凝土板和空气又重合区域，不知道建模的时候忽略没

初级学徒

我是直接在空气和土壤接触处，建立混泥土的体在划分网格时给用select/entities把混泥土板单独显示进行网格划分的，不知道是不是这步有问题 6# hmily861012

初级学徒

流固耦合除了用关键字定义，还要需要在哪里控制，有些不清楚。我是初学者，恳请赐教 6# hmily861012