【原创】读接触问题的接触力。

ansyssimulation

通过cp，轻松读接触力。
  
/clear
/prep7
*AFUN,DEG
ch=25.4/1000
  
ET,1,43
mp,ex,1,19.3e10
mp,prxy,1,0.3
  
et,5,45
mp,ex,5,19.3e10
mp,prxy,5,0.3
TB,BKIN,5
TBDATA,1,1030E6,0    ! 301 stainless steel
  
!-----------contact-----------------
ET,9,170
ET,10,173   
R,9,,,0.1,0.1,,  
MP,MU,9,0.0
  
*AFUN,DEG
  
NUMSTR,KP,1
NUMSTR,LINE,1
NUMSTR,AREA,1
k,2,0.193559*ch,-0.328883*ch  
k,19,0.145883*ch,-0.396972*ch  
l,2,19
k,3,0.199293*ch,-0.332898*ch  
k,20,0.193559*ch,-0.328883*ch  
l,3,20
k,4,0.151617*ch,-0.400987*ch  
k,21,0.199293*ch,-0.332898*ch  
l,4,21
k,5,0.115358*ch,-0.399642*ch  
k,22,0.020282*ch,-0.304567*ch  
l,5,22
k,6,0.015333*ch,-0.309517*ch  
k,23,0.110408*ch,-0.404592*ch  
l,6,23
k,7,0.007000*ch,-0.267036*ch  
k,24,0.000000*ch,-0.269869*ch  
l,7,24
k,8,0.000000*ch,-0.077500*ch  
k,25,0.007000*ch,-0.077500*ch  
l,8,25
k,9,0.000000*ch,-0.022500*ch  
k,26,0.007000*ch,-0.022500*ch  
l,9,26
k,10,0.007000*ch,-0.267036*ch  
k,27,0.007000*ch,-0.272500*ch  
l,10,27
k,11,0.000000*ch,-0.269869*ch  
k,28,0.000000*ch,-0.272500*ch  
l,11,28
k,12,0.007000*ch,-0.267036*ch  
k,29,0.007000*ch,-0.077500*ch  
l,12,29
k,13,0.000000*ch,0.000000*ch  
k,30,0.007000*ch,0.000000*ch  
l,13,30
k,14,0.000000*ch,-0.077500*ch  
k,31,0.000000*ch,-0.269869*ch  
l,14,31
k,15,0.007000*ch,-0.077500*ch  
k,32,0.007000*ch,-0.022500*ch  
l,15,32
k,16,0.007000*ch,-0.022500*ch  
k,33,0.007000*ch,0.000000*ch  
l,16,33
k,17,0.000000*ch,0.000000*ch  
k,34,0.000000*ch,-0.022500*ch  
l,17,34
k,18,0.000000*ch,-0.022500*ch  
k,35,0.000000*ch,-0.077500*ch  
l,18,35
k,36,0.052349*ch ,-0.272500*ch  
k,60,0.052349*ch ,-0.272500*ch ,1
k,64,0.007000*ch ,-0.272500*ch  
Circle,36,0.045349*ch ,60,64,45,1
k,42,0.129500*ch ,-0.385500*ch  
k,61,0.129500*ch ,-0.385500*ch ,1
k,65,0.115358*ch ,-0.399642*ch  
Circle,42,0.020000*ch ,61,65,100,1
k,48,0.129500*ch ,-0.385500*ch  
k,62,0.129500*ch ,-0.385500*ch ,1
k,66,0.110408*ch ,-0.404592*ch  
Circle,48,0.027000*ch ,62,66,100,1
k,54,0.052349*ch ,-0.272500*ch  
k,63,0.052349*ch ,-0.272500*ch ,1
k,67,0.000000*ch ,-0.272500*ch  
Circle,54,0.052349*ch ,63,67,45,1
KSEL,S,KP,,1,10000
NUMMRG,KP,0.001*ch,0.001*ch, ,LOW
  
FLST,2,4,4   
FITEM,2,12   
FITEM,2,15   
FITEM,2,16   
FITEM,2,8   
AL,P51X  
FLST,2,4,4   
FITEM,2,14   
FITEM,2,17   
FITEM,2,8   
FITEM,2,7   
AL,P51X  
FLST,2,4,4   
FITEM,2,11   
FITEM,2,13   
FITEM,2,7   
FITEM,2,6   
AL,P51X  
FLST,2,12,4  
FITEM,2,4   
FITEM,2,19   
FITEM,2,1   
FITEM,2,2   
FITEM,2,3   
FITEM,2,20   
FITEM,2,5   
FITEM,2,21   
FITEM,2,10   
FITEM,2,18   
FITEM,2,9   
FITEM,2,6   
AL,P51X  
  
ratio=0.6
  
lesize,20,0.0015*ch*ratio
ESIZE,0.003*CH*ratio
  
TYPE,1
REAL,1
MAT,1
  
AMESH,ALL
  
ALLS
  
TYPE,5
REAL,5
MAT,5
  
ESIZE,,4
VEXT,ALL,,,,,0.024*CH
  
ASEL,S,LOC,Z,0
ACLEAR,ALL
  
!-------------------------------TARGET AREA-------
  
local,11,0,,-0.4125*ch,,,-90
csys,11
wpcsys,,11
  
rectng,0.1*ch,0.2*ch,0.01*CH,-0.03*ch
  
csys,0
wpcsys,,0
  
!---------mesh target area--------
  
TYPE,9
REAL,9
MAT,9
  
esize,0.004*ch*ratio
ASEL,S,,,30
amesh,all
  
!--------mesh contact area---------
  
TYPE,10
REAL,9
MAT,9
  
ASEL,S, , ,25
NSLA,S,1
ESLN,S,0
ESURF,ALL   
  
!-----------boundary condition-----------
  
vsel,s,,,1,3
nslv,,1
D,ALL,UX
D,ALL,UY
D,ALL,UZ
  
ASEL,S,,,30
NSLA,,1
  
CP,11,UY,ALL
cplgen,11,ux,uz
*GET,nmin,NODE,,NUM,MIN, , , ,   
  
d,nmin,uz
d,nmin,ux
  
fini
  
!-------------solution -------------
  
/solu
  
antype,static
  
nlgeom,on
autots,off        !   automatic time step  
neqit,2000       !   max # of equi iter in each substep
!cnvtol,F,,0.01    !   set convergence criteria = 0.02
  
ALLS
/eof
  
save
  
y_dis=0.185-0.117
x_con=0.13
y_con=-0.41
  
!----------  first load  -------------
  
time,1
  
D,Nmin,Uy,y_dis*ch
  
nsubst,20
outres,all,all
solve
  
fini
  
/post1
  
*DIM,LABEL,CHAR,6
LABEL(1) ='RF_2(gf)'
/OUTPUT,RES_SIM,TXT,,APPEND
*VWRITE,LABEL(1)
(1X,A8)
/output
  
*do,i,1,20
set,1,i
nsel,s,,,nmin
*get,rf_2,node,NMIN,rf,fy
RF_2=RF_2/9.8*1000
/output,RES_SIM,TXT,,APPEND
*VWRITE,RF_2,i
(1X,F10.4,F10.4)
/OUTPUT
*enddo
  
PLNS,S,EQV   
ANDATA,0.5, ,2,0,0,1,1,1
  
fini
  
/eof

henlee99

请楼主解释下

sw2004

谢谢分享

逍遥游01

既然是分享,有些地方做一下注释还是必须的!

kevin210099

thanks you a lot!

lhw_heavy

感谢！！学习中~