球体侵彻多层靶,接触问题,网格自身重叠?
各位专家、牛人大咖您们好,我最近在做一个侵彻例子。球体和两层靶板之间的接触为侵蚀接触,两层靶板之间的接触为自动面面接触,并且各part间都留有间隙。调试了很久,依旧无法解决网格自身重叠的问题。球体与第一层靶板接触部分区域网格并没有侵蚀掉,导致最后靶板网格交叉重叠在一起,而且计算也没有报错停止,大家看图吧。模型也重新建了,时间步长我也调小过(曾经调到0.1),材料参数也没有问题。我估摸着肯定是接触定义有问题,但我实在不知道从哪入手调试了,附上我的控制文件,接触定义。希望各位能够给予帮助。part1为球体,part2为第一层靶板,part3为第二层靶板。
*CONTACT_ERODING_SURFACE_TO_SURFACE
$# cid title
$# ssid msid sstyp mstyp sboxid mboxid spr mpr
2 1 3 3 0 0 1 1
$# fs fd dc vc vdc penchk bt dt
0.150000 0.100000 0.000 0.000 20.000000 2 0.000 1.0000E+7
$# sfs sfm sst mst sfst sfmt fsf vsf
1.000000 1.000000 0.000 0.000 1.000000 1.000000 1.000000 1.000000
$# isym erosop iadj
1 1 0
$# soft sofscl lcidab maxpar sbopt depth bsort frcfrq
2 0.100000 0 1.025000 0.000 2 0 0
$# penmax thkopt shlthk snlog isym i2d3d sldthk sldstf
0.000 0 0 0 0 0 0.000 0.000
$# igap ignore dprfac dtstif unused unused flangl cid_rcf
2 0 0.000 0.000 0 0 0.000 0
*CONTACT_ERODING_SURFACE_TO_SURFACE_ID
$# cid title
2
$# ssid msid sstyp mstyp sboxid mboxid spr mpr
3 1 3 3 0 0 1 1
$# fs fd dc vc vdc penchk bt dt
0.150000 0.100000 0.000 0.000 20.000000 2 0.000 1.0000E+7
$# sfs sfm sst mst sfst sfmt fsf vsf
1.000000 1.000000 0.000 0.000 1.000000 1.000000 1.000000 1.000000
$# isym erosop iadj
1 1 0
$# soft sofscl lcidab maxpar sbopt depth bsort frcfrq
2 0.100000 0 1.025000 0.000 2 0 0
$# penmax thkopt shlthk snlog isym i2d3d sldthk sldstf
0.000 0 0 0 0 0 0.000 0.000
$# igap ignore dprfac dtstif unused unused flangl cid_rcf
2 0 0.000 0.000 0 0 0.000 0
*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE
$# cid title
$# ssid msid sstyp mstyp sboxid mboxid spr mpr
3 2 3 3 0 0 1 1
$# fs fd dc vc vdc penchk bt dt
0.150000 0.100000 0.000 0.000 20.000000 2 0.000 1.0000E+7
$# sfs sfm sst mst sfst sfmt fsf vsf
1.000000 1.000000 0.000 0.000 1.000000 1.000000 1.000000 1.000000
$# soft sofscl lcidab maxpar sbopt depth bsort frcfrq
2 0.100000 0 1.025000 0.000 2 0 1
$# penmax thkopt shlthk snlog isym i2d3d sldthk sldstf
0.000 0 0 0 0 0 0.000 0.000
$# igap ignore dprfac dtstif unused unused flangl cid_rcf
2 0 0.000 0.000 0 0 0.000 0
*CONTROL_CONTACT
$ MID RO G SIGY EH PC FS
$# slsfac rwpnal islchk shlthk penopt thkchg orien enmass
2.000000 0.000 2 0 1 0 1 0
$# usrstr usrfrc nsbcs interm xpene ssthk ecdt tiedprj
0 0 10 0 4.000000 0 0 0
$# sfric dfric edc vfc th th_sf pen_sf
0.000 0.000 0.000 0.000 0.000 0.000 0.000
$# ignore frceng skiprwg outseg spotstp spotdel spothin
0 0 0 0 0 0 0.000
$# isym nserod rwgaps rwgdth rwksf icov swradf ithoff
0 0 1 0.000 1.000000 0 0.000 0
$# shledg pstiff ithcnt tdcnof ftall unused shltrw
0 0 0 0 0 0 0.000
有什么问题, 单元刚度退化叠在一起有问题么? 充塞破坏也是一种破坏形式,没有谁说弹体前方的单元就一定要消失? 我觉得这个结果挺好