在hekie等各位热心人的帮助下,我上次上传的关于异形板的问题后来改进后,采用了shell63单元进行分析,终于得出了结果,可是结果还是出现了一些问题:
1.网格划分后总是出现一块带斑的区域,而且总是警告单元形状有问题
2.由于我是在实体上加载约束的,然后粘结,划分网格,再求解,可是这样的结果却是57~60号点的约束神奇般的消失了,所以解出来的结果也当然有问题了,为什么这些约束会消失呢?
以下是我的命令流和x方向应力图:
/prep7
/pnum,kp,1
/pnum,line,1
/pnum,area,1
et,1,63
! c30
mp,ex,1,3.00e10
mp,nuxy,1,0.2
mp,dens,1,2.5e3
r,1,0.65
!
k,1,0,0
k,2,0,-36.3
k,3,45.9,-34.3
k,4,90.5,-35.1
k,5,112.3,-49.2
k,6,115.2,-49.2
k,7,117,-39
k,8,117,0
l,1,2
k,9,24,-40
k,10,67,-33
k,11,103,-44
larc,2,3,9,110
larc,3,4,10,150
larc,4,5,11,20
l,5,6
l,6,7
l,7,8
l,8,1
kdele,9
kdele,10
kdele,11
local,11,1,70.5,-30
csys,11
k,,1.1,-90
k,,1.1
k,,1.1,90
k,,1.1,180
l,9,10
l,10,11
l,11,12
l,12,9
local,12,1,83,-30
csys,12
k,,2.2,-90
k,,2.2
k,,2.2,90
k,,2.2,180
l,13,14
l,14,15
l,15,16
l,16,13
local,13,1,97.5,-29.4
csys,13
k,,3.25,-90
k,,3.25
k,,3.25,90
k,,3.25,180
l,17,18
l,18,19
l,19,20
l,20,17
csys,0
k,,13,-3
k,,26,-3
k,,39,-3
k,,52,-3
k,,65,-3
k,,78,-3
k,,91,-3
k,,104,-3
k,,13,-9
k,,26,-9
k,,39,-9
k,,52,-9
k,,65,-9
k,,78,-9
k,,91,-9
k,,104,-9
k,,13,-15
k,,26,-15
k,,39,-15
k,,52,-15
k,,65,-15
k,,78,-15
k,,91,-15
k,,104,-15
k,,13,-21
k,,26,-21
k,,39,-21
k,,52,-21
k,,65,-21
k,,78,-21
k,,91,-21
k,,104,-21
k,,13,-28.8
k,,26,-27.8
k,,39,-28.4
k,,52,-30.2
k,,65,-31.3
k,,78,-31.3
k,,91,-30.1
k,,104,-32.3
a,1,2,53,45,37,29,21
a,1,21,22,23,24,25,26,27,28,8
a,21,29,30,31,32,33,34,35,36,28
a,29,37,38,39,40,41,42,43,44,36
a,37,45,46,47,48,49,50,51,52,44
a,8,28,52,7
l,2,53
l,53,54
l,54,55
l,55,3
al,2,66,65,64,21
a,45,53,54,55,47
l,55,12
l,11,49
al,67,69,11,70,56,55
l,4,20
l,20,14
l,16,10
al,66,3,71,72,13,16,73,9,12,69
l,18,52
al,71,4,5,6,63,74,17,20
l,15,51
al,70,10,73,15,75,58,57
al,14,72,19,18,74,59,75
ksel,s,,,21,60,1
dk,all,all
acel,0,0,10
asel,all
aglue,all
mshkey,0
mshape,0,2d
esize,1
amesh,all
/solu
solve |