SPH接触爆炸,粒子飞散讨论

浏览:3192 评论:6 收藏:2
炸药和薄板均是SPH建模,在爆炸过程中,薄板应该向下弯曲,但是在模拟时出现如图所示的破裂,粒子快速飞散,炸药粒子甚至直接穿透薄板飞出。炸药和薄板参数均选自AUTODYN模型库。后来我发现,建模后,将薄板向下移动一个粒子尺寸的距离,这种现象就消失了。但这不是真实的原因,因为在autodyn里接触建模是没有任何问题的。请教造成这种现象的真实原因是什么?应该如何改进?附上我的K文件,单位cm-g-us
模型是二维的,Shell填充粒子生成的。文件太大,上传不了,我就把关键字段附上,删去了NODE 和ELEMENT SPH两部分的信息。
*CONTROL_BULK_VISCOSITY
$# q1 q2 type btype
1.5 0.06 1 0
*CONTROL_ENERGY
$# hgen rwen slnten rylen
2 2 2 2
*CONTROL_SPH
$# ncbs boxid dt idim memory form start maxv
0 01.00000E20 2 150 5 0.00.69999999
$# cont deriv ini ishow ierod icont iavis
0 0 0 0 0 0 0
*CONTROL_TERMINATION
$# endtim endcyc dtmin endeng endmas
20.0 0 0.0 0.01.000000E8
*CONTROL_TIMESTEP
$# dtinit tssfac isdo tslimt dt2ms lctm erode ms1st
0.00.67000002 0 0.0 0.0 0 0 0
$# dt2msf dt2mslc imscl unused unused rmscl
0.0 0 0 0.0
*DATABASE_GLSTAT
$# dt binary lcur ioopt
0.2 0 0 1
*DATABASE_MATSUM
$# dt binary lcur ioopt
0.2 0 0 1
*DATABASE_SPHOUT
$# dt binary lcur ioopt
0.2 0 0 1
*DATABASE_BINARY_D3PLOT
$# dt lcdt beam npltc psetid
0.02 0 0 0 0
$# ioopt
0
*DATABASE_EXTENT_BINARY
$# neiph neips maxint strflg sigflg epsflg rltflg engflg
0 0 3 1 1 1 1 1
$# cmpflg ieverp beamip dcomp shge stssz n3thdt ialemat
0 0 4 1 1 1 2 1
$# nintsld pkp_sen sclp hydro msscl therm intout nodout
0 0 1.0 0 0 0STRESS
$# dtdt resplt neipb
0 0 0
*PART
$# title
SphNode
$# pid secid mid eosid hgid grav adpopt tmid
3 2 3 3 0 0 0 0
*SECTION_SPH
$# secid cslh hmin hmax sphini death start
2 1.2 0.2 2.0 0.01.00000E20 0.0
*MAT_HIGH_EXPLOSIVE_BURN
$# mid ro d pcj beta k g sigy
30.930999990.41600001 0.0515 0.0 0.0 0.0 0.0
*EOS_JWL
$# eosid a b r1 r2 omeg e0 vo
3 0.4946 0.01891 3.9070001 1.1180.33333334 0.02484 1.0
*PART
$# title
SphNode
$# pid secid mid eosid hgid grav adpopt tmid
4 2 4 4 0 0 0 0
*MAT_JOHNSON_COOK
$# mid ro g e pr dtf vp rateop
4 8.93999960.46000001 0.0 0.0 0.0 0.0 0.0
$# a b n c m tm tr epso
9.00000E-4 0.00292 0.31 0.0 1.09 1356.0 293.0 0.0
$# cp pc spall it d1 d2 d3 d4
4.60000E-6 -9.0 3.0 0.00.54000002 4.8899999 -3.03 0.014
$# d5 c2/p erod efmin
1.12 0.0 01.00000E-6
*EOS_GRUNEISEN
$# eosid c s1 s2 s3 gamao a e0
4 0.391 1.51 0.0 0.0 1.97 0.0 0.0
$# v0
0.11206
*INITIAL_DETONATION
$# pid x y z lt
3 0.0 1.1 0.0 0.0
*ELEMENT_SPH
$# nid pid mass
18313 3 3.723993e-004
18314 3 3.723999e-004
35806 4 8.939996e-004
35807 4 8.939996e-004
35808 4 8.939999e-004
35809 4 8.939999e-004
35810 4 8.939999e-004
35811 4 8.939999e-004
35812 4 8.939999e-004

技术邻APP
工程师必备
工程师必备
- 项目客服
- 培训客服
- 平台客服
TOP

6
2