GAPCYL单元的应力问题

zbingxing

用gap单元做两根梁的接触分析，在一根梁上加垂直向上的力，提取mises应力以后在gap单元处应力为零，可是连续梁的应力应该是出处保持一致的啊。附上inp文件，求高手解释









*heading
test with gap contact
*****************************
*****************************
*node,nset=one
1,0,0,0
5,0,0,2
*ngen,nset=one
1,5
*element,type=pipe31h,elset=one
1,1,2
*elgen,elset=one
1,4
*beam general section,elset=one,section=pipe,density=7850,poisson=0.3
0.1619,0.0127
1,0,0,
2.07e11, 7.85e10
*Section Points
0.1619, 0., 0., 0.1619, -0.1619, 0., 0., -0.1619,0.1145,0.1145,-0.1145,0.1145,-0.1145,-0.1145,0.1145,-0.1145
*node,nset=two
1001,0,0,0
1005,0,0,2
*ngen,nset=two
1001,1005
*element,type=pipe31h,elset=two
1001,1001,1002
*elgen,elset=two
1001,4
*beam general section,elset=two,section=pipe,density=7850,poisson=0.3
0.11225,0.01383
1,0,0
2.07e11, 7.85e10
*Section Points
0.11225, 0., 0., 0.11225, -0.11225, 0., 0., -0.11225,0.07938,0.07938,-0.07938,0.07938,-0.07938,-0.07938,0.07938,-0.07938
*ELEMENT,TYPE=GAPCYL,ELSET=GAP1
2002,2,1002
2003,3,1003
2004,4,1004
*GAP,ELSET=GAP1
0.01,0.,0.,1,0.01
*STEP , name=Step-displace-static, NLGEOM=yes, INC=8000
   apply load
*STATIC
0.0005, 0.01, 5e-20, 0.01
*boundary
1   , ENCASTRE
*boundary
1001, ENCASTRE
*CLOAD
**5,2,1833000
5,3,1833000
*Restart, write, frequency=0
**
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, frequency=5
*Node output
CF,RF,U
*Element output,directions=YES
MISESMAX,S,SF
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT, frequency=1
*End Step
******