3 LSDYNA 爆破模拟出现能量比问题，提示 termimation due to energy，应该如何处理
现在注册，即可获得 4 信用和 4 瓦能量注册奖励。
爆炸模拟 K 运行出现了能量比的错误（如下图），得如何处理呢？
模型如下，炸药外面一层是空气，再外面是岩石，用流固耦合的方式模拟
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
$
- 2020.05.16 17:04
- 5622
- 0
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}\) 上，出现该问题的原因在于过度的沙漏，过度接触等造成的人为吸收了太多的能量。
