嵌入问题!

xhill

我用拉格朗日方法计算厚靶的侵彻问题!用侵蚀接触!接触刚度改为1.0了!可是计算后还是发生了嵌入问题!请各位指点一下,是怎么回事?多谢了!
下面是文件:
*KEYWORD
*TITLE
balls7.                                                                              
$
*DATABASE_FORMAT
         0
$
*Include
1_Node.k
*Include
2_Element.k
*Include
3_Bound.k
*Include
4_Contact.k
*Include
5_velocity.k
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             SECTION DEFINITIONS                              $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*SECTION_SOLID
         1         1
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             MATERIAL DEFINITIONS                             $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*MAT_RIGID
         1  7.80      2.40      0.289000       0.0       0.0       0.0
  1.00      0.00      0.00     
  
*MAT_JOHNSON_COOK
$      MID        RO     G         E        PR        DTF       VP  
         2   2.20000 0.125     0.296      0.180000  0.100000   0.0
$        A         B     N         C         M         TM       TR      EPSO
  0.238E-03 0.174E-01  1.00     0.474      1.00     0.300E+04  20.0    1.0e-6
$       CP        PC  SPALL       IT        D1        D2        D3        D4
  0.690E-05  -4.4e-4   3.00      1.00      0.04      1.00      0.00      0.00     
$   D5  
  0.00     
$
*EOS_GRUNEISEN
         2  0.23      0.83      0.00      0.00      0.66     -0.66      0.00     
  0.00   
*MAT_ADD_EROSION
$      MID   EXCL
         2  0.00
$ PFAIL     SIGP1    SIGVM      EPSP1    EPSSH     SIGTH     IMPULSE
  0.00      0.00      0.00      0.31      0.00     
*MAT_RIGID
         3  3.80      3.30      0.240000       0.0       0.0       0.0
  1.00      0.00      0.00     
  
$
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              PARTS DEFINITIONS                               $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$
*PART
Part          1 for Mat         2 and Elem Type         1
         1         1         2         2         0         0
$
*PART
Part          2 for Mat         1 and Elem Type         1
         2         1         1         0         0         0
$
*PART
Part          3 for Mat         3 and Elem Type         1
         3         1         3         0         0         0
$
*PART
Part          4 for Mat         4 and Elem Type         1
         4         1         3         0         0         0
$
*PART
Part          5 for Mat         5 and Elem Type         1
         5         1         3         0         0         0
$
*PART
Part          6 for Mat         6 and Elem Type         1
         6         1         3         0         0         0
$
*PART
Part          7 for Mat         7 and Elem Type         1
         7         1         3         0         0         0
$
*PART
Part          8 for Mat         8 and Elem Type         1
         8         1         3         0         0         0
$
*PART
Part          9 for Mat         9 and Elem Type         1
         9         1         3         0         0         0
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              CONTROL OPTIONS                                 $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*Control_Contact
$   slsfac    rwpnal    islchk    shlthk    penopt    thkchg     orien
       1.0       1.0         1         1         2
$   usrstr    usrfac     nsbcs    interm    xpenen     SSTHK      ECDT   TIEDPRJ
                             1                  4.
*CONTROL_BULK_VISCOSITY
       1.0      0.06         1
*CONTROL_HOURGLASS
         1      0.10
*CONTROL_ENERGY
         2         2         2         2
*CONTROL_SHELL
  20.0             1        -1         1         2         2         1
*CONTROL_TIMESTEP
    0.0000    0.9000         0  0.00      0.00     
*CONTROL_TERMINATION
  0.150E+04         0   0.00000   0.00000   0.00000
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              DATABASE OPTIONS                                $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*DATABASE_BINARY_D3PLOT
  150.0     
*DATABASE_EXTENT_BINARY
         0         0         3         1         0         0         0         0
         0         0         0         0         0         0
*END
接触控制:
*keyword
$弹、球和靶体的接触
*CONTACT_ERODING_SURFACE_TO_SURFACE
         1         2         3         3         0         0         1         1
    0.0000    0.0000    0.0000    0.0000    0.0000         0    0.00000.1000E+08
    1.0000    1.0000    0.0000    0.0000    1.0000    1.0000    1.0000    1.0000
         1         1         1
*CONTACT_ERODING_SURFACE_TO_SURFACE
         1         3         3         3         0         0         1         1
    0.0000    0.0000    0.0000    0.0000    0.0000         0    0.00000.1000E+08
    1.0000    1.0000    0.0000    0.0000    1.0000    1.0000    1.0000    1.0000
         1         1         1
*CONTACT_ERODING_SURFACE_TO_SURFACE
         1         4         3         3         0         0         1         1
    0.0000    0.0000    0.0000    0.0000    0.0000         0    0.00000.1000E+08
    1.0000    1.0000    0.0000    0.0000    1.0000    1.0000    1.0000    1.0000
         1         1         1
*CONTACT_ERODING_SURFACE_TO_SURFACE
         1         5         3         3         0         0         1         1
    0.0000    0.0000    0.0000    0.0000    0.0000         0    0.00000.1000E+08
    1.0000    1.0000    0.0000    0.0000    1.0000    1.0000    1.0000    1.0000
         1         1         1
*CONTACT_ERODING_SURFACE_TO_SURFACE
         1         6         3         3         0         0         1         1
    0.0000    0.0000    0.0000    0.0000    0.0000         0    0.00000.1000E+08
    1.0000    1.0000    0.0000    0.0000    1.0000    1.0000    1.0000    1.0000
         1         1         1
*CONTACT_ERODING_SURFACE_TO_SURFACE
         1         7         3         3         0         0         1         1
    0.0000    0.0000    0.0000    0.0000    0.0000         0    0.00000.1000E+08
    1.0000    1.0000    0.0000    0.0000    1.0000    1.0000    1.0000    1.0000
         1         1         1
*CONTACT_ERODING_SURFACE_TO_SURFACE
         1         8         3         3         0         0         1         1
    0.0000    0.0000    0.0000    0.0000    0.0000         0    0.00000.1000E+08
    1.0000    1.0000    0.0000    0.0000    1.0000    1.0000    1.0000    1.0000
         1         1         1
*CONTACT_ERODING_SURFACE_TO_SURFACE
         1         9         3         3         0         0         1         1
    0.0000    0.0000    0.0000    0.0000    0.0000         0    0.00000.1000E+08
    1.0000    1.0000    0.0000    0.0000    1.0000    1.0000    1.0000    1.0000
         1         1         1
$弹和球的接触
*CONTACT_ERODING_SURFACE_TO_SURFACE
         5         2         3         3         0         0         1         1
    0.0000    0.0000    0.0000    0.0000    0.0000         0    0.00000.1000E+08
    1.0000    1.0000    0.0000    0.0000    1.0000    1.0000    1.0000    1.0000
         1         1         1
*CONTACT_ERODING_SURFACE_TO_SURFACE
         6         2         3         3         0         0         1         1
    0.0000    0.0000    0.0000    0.0000    0.0000         0    0.00000.1000E+08
    1.0000    1.0000    0.0000    0.0000    1.0000    1.0000    1.0000    1.0000
         1         1         1
*CONTACT_ERODING_SURFACE_TO_SURFACE
         7         2         3         3         0         0         1         1
    0.0000    0.0000    0.0000    0.0000    0.0000         0    0.00000.1000E+08
    1.0000    1.0000    0.0000    0.0000    1.0000    1.0000    1.0000    1.0000
         1         1         1
$球之间的接触
*CONTACT_ERODING_SURFACE_TO_SURFACE
         3         4         3         3         0         0         1         1
    0.0000    0.0000    0.0000    0.0000    0.0000         0    0.00000.1000E+08
    1.0000    1.0000    0.0000    0.0000    1.0000    1.0000    1.0000    1.0000
         1         1         1
*CONTACT_ERODING_SURFACE_TO_SURFACE
         4         5         3         3         0         0         1         1
    0.0000    0.0000    0.0000    0.0000    0.0000         0    0.00000.1000E+08
    1.0000    1.0000    0.0000    0.0000    1.0000    1.0000    1.0000    1.0000
         1         1         1
*CONTACT_ERODING_SURFACE_TO_SURFACE
         5         6         3         3         0         0         1         1
    0.0000    0.0000    0.0000    0.0000    0.0000         0    0.00000.1000E+08
    1.0000    1.0000    0.0000    0.0000    1.0000    1.0000    1.0000    1.0000
         1         1         1
*CONTACT_ERODING_SURFACE_TO_SURFACE
         7         6         3         3         0         0         1         1
    0.0000    0.0000    0.0000    0.0000    0.0000         0    0.00000.1000E+08
    1.0000    1.0000    0.0000    0.0000    1.0000    1.0000    1.0000    1.0000
         1         1         1
*CONTACT_ERODING_SURFACE_TO_SURFACE
         8         7         3         3         0         0         1         1
    0.0000    0.0000    0.0000    0.0000    0.0000         0    0.00000.1000E+08
    1.0000    1.0000    0.0000    0.0000    1.0000    1.0000    1.0000    1.0000
         1         1         1
*CONTACT_ERODING_SURFACE_TO_SURFACE
         9         8         3         3         0         0         1         1
    0.0000    0.0000    0.0000    0.0000    0.0000         0    0.00000.1000E+08
    1.0000    1.0000    0.0000    0.0000    1.0000    1.0000    1.0000    1.0000
         1         1         1
$
*end

xhill

初始图形

xhill

嵌入图形局部放大

xhill

用segment接触可以避免吗?

xhill

那位大虾帮帮忙啊!

belldingding

最好把K文件以附件的形式给出，这样别人就可以方便给你调试了