3 LSDYNA 爆破模拟出现能量比问题,提示 termimation due to energy,应该如何处理

爆炸模拟 K 运行出现了能量比的错误(如下图),得如何处理呢?attachments-2020-05-jqwAVhdW5ebfab8da4a47.png

模型如下,炸药外面一层是空气,再外面是岩石,用流固耦合的方式模拟

attachments-2020-05-Y2BwiEY55ebfac896b5e3.png

K 文件如下

*SECTION_SOLID
$   1SECID   2ELFORM      3AET
         1        13          
*SECTION_SOLID
$   1SECID   2ELFORM      3AET
         2        11    
               
*SECTION_SOLID
$   1SECID   2ELFORM      3AET
         3        11          
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             MATERIAL DEFINITIONS                             $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*MAT_PLASTIC_KINEMATIC
1,1768,2E+9,0.3,2E+5,0,0
0,0,0,0                          
*MAT_NULL
2,1.25,-1,1.7456e-5
*EOS_LINEAR_POLYNOMIAL
1,0,0,0,0,0.4,0.4,0
2.5E+5,1                            
*MAT_HIGH_EXPLOSIVE_BURN
$      MID        RO         D       PCJ      BETA         K         G      SIGY
3,1000,4500,5.06E+9
*EOS_JWL
$    EOSID         A         B        R1        R2      OMEG        E0        V0
2,178.85E+9,0.311E+9,4.75,1.05,0.18,2.25E+9   
*INITIAL_DETONATION
$      PID         X         Y         Z        LT
         3         0         0         0         0   
*ALE_MULTI-MATERIAL_GROUP_PART                                                   
2                                                                               
*ALE_MULTI-MATERIAL_GROUP_PART                                                   
3      
*CONTROL_ALE
$$     DCT      NADV      METH      AFAC      BFAC      CFAC      DFAC      EFAC
         2         1         2      -1.0       0.0       0.0       0.0          
$$   START       END     AAFAC     VFACT    VLIMIT       EBC      PREF   NSIDEBC
       0.0       0.0       0.0    
*CONSTRAINED_Lagrange_IN_SOLID
$#   slave    master     sstyp     mstyp     nquad     ctype     direc     mcoup
         1         2         1         1         3         4         3         0
$#   start       end      pfac      fric    frcmin      norm   normtyp      damp
     0.000       2.0  0.100000     0.000  0.300000         0
$#      cq      hmin      hmax     ileak     pleak   lcidpor     nvent  blockage
     0.000     0.000     0.000         0  0.010000   
*CONSTRAINED_Lagrange_IN_SOLID
$#   slave    master     sstyp     mstyp     nquad     ctype     direc     mcoup
         1         3         1         1         2         4         2         0
$#   start       end      pfac      fric    frcmin      norm   normtyp      damp
     0.000       2.0  0.100000     0.000  0.300000         1
$#      cq      hmin      hmax     ileak     pleak   lcidpor     nvent  blockage
     0.000     0.000     0.000         0  0.010000                               
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              PARTS DEFINITIONS                               $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*PART
$  HEADING
SYS\__
$     1PID    2SECID      3MID    4EOSID     5HGID     6GRAV   7ADPORT     8TMID
         1         1         1         0         0                              
*PART
$  HEADING
SYS\__
$     1PID    2SECID      3MID    4EOSID     5HGID     6GRAV   7ADPORT     8TMID
         2         2         2         1         0                              
*PART
$  HEADING
SYS\__
$     1PID    2SECID      3MID    4EOSID     5HGID     6GRAV   7ADPORT     8TMID
         3         3         3         2         0                              
$      
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             COORDINATE SYSTEMS                               $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*DEFINE_COORDINATE_SYSTEM
$     1CID       2XO       3YO       4ZO       5XL       6YL       7ZL
         1     0.187      0.45     0.809     1.187      0.45     0.809
$      1XP       2YP       3ZP
     0.187      1.45     0.809
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                               LOAD DEFINITIONS                               $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*DEFINE_CURVE
$    1LCID     2SIDR      3SFA      4SFO     5OFFA     6OFFO   7DATTYP
         1                                                            
$                 1A                  2O
                   0              1e+008
                0.01              1e+008
                 0.1              1e+008
*DEFINE_CURVE
$    1LCID     2SIDR      3SFA      4SFO     5OFFA     6OFFO   7DATTYP
         2                                                            
$                 1A                  2O
                   0             9.80665
                0.01             9.80665
                 0.1             9.80665
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              CONTACT DEFINITIONS                             $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*CONTACT_AUTOMATIC_SINGLE_SURFACE
$    1SSID     2MSID    3SSTYP    4MSTYP   5SBOXID   6MBOXID      7SPR      8MPR
         0         0         5         0                             1         1
$      1FS       2FD       3DC       4VC      5VDC   6PENCHK       7BT       8DT
         0         0         0         0        10                              
$     1SFS      2SFM      3SST      4MST     5SFST     6SFMT      7FSF      8VSF
                                                                                
$    1SOFT   2SOFSCL   3LCIDAB   4MAXPAR    5SBOPT    6DEPTH    7BSORT   8FRCFRQ
         2                                       3         5                    
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                                CONTROL OPTIONS                               $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*CONTROL_TERMINATION
$  1ENDTIM   2ENDCYC    3DTMIN   4ENDENG   5ENDMAS
     0.001  10000000      0.01        10         0
*CONTROL_TIMESTEP
$  1DTINIT   2TSSFAC     3ISDO   4TSLIMT    5DT2MS     6LCTM    7ERODE    8MS1ST
         0       0.6         0         0         0         1         1         0
*CONTROL_Hourglass
$     1IHQ       2QH
         1       0.1
*CONTROL_BULK_VISCOSITY
$      1Q1       2Q2     3TYPE
       1.5      0.06        -2
*CONTROL_CONTACT
$  1SLSFAC   2RWPNAL   3ISLCHK   4SHLTHK   5PENOPT   6THKCHG    7ORIEN   8ENMASS
         0         0         1         1         1         0         2         0
$ 1USRSTRC   2USRFRC    3NSBCS   4INTERM    5XPENE    6SSTHK     7ECDT  8TIEDPRJ
        
$   1SFRIC    2DFRIC      3EDC      4VFC       5TH    6TH_SF   7PEN_SF
       
$  1IGNORE   2FRCENG  3SKIPRWG   4OUTSEG  5SPOTSTP  6SPOTDEL 7SPOTTHIN
         2         0         0         1         0         1       0.5
*CONTROL_SOLID
$   1ESORT   2FMATRX  3NIPTETS   4SWLOCL
         1   
*Damping_GLOBAL
$    1LCID   2VALDMP      3STX      4STY      5STZ      6SRX      7SRY      8SRZ
         0         0      
*CONTROL_ENERGY
$    1HGEN     2RWEN   3SLNTEN    4RYLEN
         2         1         2         2
*CONTROL_ACCURACY
$     1OSU      2INN   3PIDOSU
         1         4 
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                                 TIME HISTORY                                 $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*DATABASE_GLSTAT
$      1DT   2BINARY     3LCUR    4IOOPT    5DTHFF    6BINHF
      8E-6
*DATABASE_MATSUM
$      1DT   2BINARY     3LCUR    4IOOPT    5DTHFF    6BINHF
      8E-6
*DATABASE_NODOUT
$      1DT   2BINARY     3LCUR    4IOOPT    5DTHFF    6BINHF
      8E-6
*DATABASE_ELOUT
$      1DT   2BINARY     3LCUR    4IOOPT    5DTHFF    6BINHF
      8E-6
*DATABASE_BINARY_d3plot
$      1DT     2LCDT     3Beam    4NPLTC
      8E-6
*DATABASE_BINARY_RUNRSF
$      1DT       2NR
      5000
$
  • 163
请先 登录 后评论

2 个回答

首先了解一下 energy ratio 是怎么回事,数值仿真中,energy ratio 的计算方法如下:

\[eratio=\frac{E_{total}}{E_{total}+W_{ext}}\]

这其中,\(E_{total}\) 可按如下公式计算:

\[E_{total}=E_{kin}^0+E_{int}^0+W_{ext}\]

那么,energy ratio 的表达式,即可写成如下的形式:

\[eratio=\frac{E_{kin}^0+E_{int}^0+W_{ext}}{E_{kin}^0+E_{int}^0+2W_{ext}}\]

表达式中,有三个变量,分别为,\(E_{kin}^0\)、\(E_{int}^0\)、\(W_{ext}\),前两个变量分别为弹性应变能、变形功,第三个变量为外部功,包括力、压力的做功,以及速度、位移、加速度等于边界条件上的作用功。

从报错信息看,错误点应该在能量功 \(W_{ext}\) 上,出现该问题的原因在于过度的沙漏,过度接触等造成的人为吸收了太多的能量。

请先 登录 后评论
柳捷 - 自由

attachments-2020-05-1QW9byj45ec0c974ecf6b.png

1. 接触关键字 修改试试

2. 试试质量缩放系数

3. 六面体网格

请先 登录 后评论