各位大哥好~
本人在用非线性接触模拟,两个叠放的承压空心板(上板厚40紧压在下板厚40上,上板(接触面)把压力传递给下板,下板(目标面)底部有环形竖向约束),结果提示错误,上板的边界处位移过大,1/4部分的命令流如下,请了解的大神们指教~~ 谢谢了。
finish
/clear
/filname,chengyaban-case1,1
/TITLE, case1
/PREP7
blc5,12.5,12.5,25,25
blc5,17.5,17.5,35,35
blc5,22.5,22.5,45,45
blc5,30,30,60,60
wpoff,,,40
blc5,12.5,12.5,25,25
blc5,22.5,22.5,45,45
asbl,all,all
adele,1,5,4
nummrg,all
numcmp,all
adele,1,2,1
nummrg,all
numcmp,all
asel,s,area,,1,3,1
cm,n1,area
asel,s,area,,4
cm,n2,area
asel,s,area,,1,2,1
cm,n3,area
allsel
!!!定义实常数!!!!
ET,1,SHELL181
ET,2,SHELL181
ET,3,TARGE170 !!目标单元.
ET,4,CONTA173
KEYOPT, 4, 2, 0
KEYOPT, 4, 9, 1
KEYOPT, 4, 11, 1
KEYOPT, 4, 12, 1
R,1,40,,,, , , !N1.
R,2,40,,,, , , !N2.
R,3
R,3,,,0.1,0.1,,
!定义钢材材料.
MP,EX,1,2.06e5
MP,PRXY,1,0.3
MP,DENS,1,7.85E-6
/VIEW, 1, 1, 1, 1
allsel
nummrg,kp
numcmp,kp
numcmp,line
numcmp,area
!加载荷载.!!约束.
lsel,s,line,,5,6,1
dl,all,,uz
dl,all,,uy
dl,all,,ux
!dl,all,,rotx.
!dl,all,,roty.
!dl,all,,rotz.
!!!划分网格.
allsel
lsla,s !相关线的选择.
esize,5 !线的单元尺寸定义.
vsel,none
asel,none
esel,none
cmsel,s,n1
aatt,1,1,1
amesh,all
cm,e1,elem
esel,none
cmsel,s,n2
aatt,1,2,2
amesh,all
cm,e2,elem
esel,none
allsel
save
!定义目标面.
ASEL,S,area,,1,2,1
NSLA,S,1
ESLN,S,0
TYPE, 3
REAL, 3
ESURF,all
ALLS
!定义接触面.
ASEL,S,area,,4
NSLA,S,1
ESLN,S,0
TYPE, 4
REAL, 3
ESURF,all
!!荷载.
allsel
cmsel,s,n2
time,1 !定义载荷步1.
sfa,all,,pres,-17.857
allsel
/PBC,ALL, ,1
/PSF,PRES,NORM,2,0,1
!!显示面上线力方向.
allsel
/SOLU !进入求解器.
Dl,17,,SYMM !定义面的对称位移边条.
Dl,18,,SYMM
Dl,19,,SYMM
Dl,20,,SYMM
Dl,21,,SYMM
Dl,22,,SYMM
Dl,23,,SYMM
Dl,24,,SYMM
solve |