接触分析学习求助

cumtkw

我在本论坛接触板块看到好多命令流，自己下载下来进行了修改和并运行分析，得到了结果，但是不知道如何判别运行结果是不是正确，没有一个答案 ，不知道大家在论坛上找到实例进行学习时，如何判别自己运算的结果是否合理？下面的例子是我运行过的，两种运行下来结果不太一样，我就不知道，我没两个谁的错了？？？遇到这种情况怎么办？

论坛原始命令流：
eex=2e11
pprxy=0.3
ff=1000
fe=0.3  !摩擦系数
fk=5  !接触刚度
fd=3  !切向接触刚度
/prep7   
et,1,solid45
mp,ex,1,eex
mp,prxy,1,pprxy
MP,DENS,1,7.8e3  
MP,MU,1,fe  
block,0,1,0,1,0,1   
block,0,1,0,1,1,2   
lesize,all,0.1   
vmesh,all   
  
ET,2,170
ET,3,174
keyopt,3,10,2
R,2, , ,fk, , , ,
!RMORE, , , , , ,fd,   
  
ASEL,S,,,2   
  
NSLA,S,1
ESLN,S,0
r,1
TYPE,2  
ESURF,ALL   
allsel,all
  
ASEL,S,,,7  
NSLA,S,1
ESLN,S,0
r,2,,,fk,
RMORE, , , , , ,fd
type,3
ESURF,ALL   
allsel,all
save
fini
  
/SOLU   
dk,1,ALL,0
dk,2,all,0
dk,3,all,0
dk,4,all,0   
SFA,12,1,PRES,ff
sftran
ACEL,0,0,9.8, !施加重力加速度
nlgeom,on
auto,on
nsubst,100,500,50
outres,all,all
save
solve
fini


我自己修改之后的命令流：
eex=2e11
pprxy=0.3
ff=1000
fe=0.3  !摩擦系数
fk=1  !接触刚度
fd=3  !切向接触刚度
/prep7   
et,1,solid45
mp,ex,1,eex
mp,prxy,1,pprxy
MP,DENS,1,7.8e3  
MP,MU,1,fe  
block,0,1,0,1,0,1   
block,0,1,0,1,1,2   
lesize,all,0.1   
vmesh,all   
  
ET,2,170
ET,3,174
keyopt,3,10,2    !? 为什么加这一条命令
R,2, , ,fk,1 , , ,
RMORE, , , , , ,fd,   
  
ASEL,S,,,2  
NSLA,S,1
!ESLN,S,0
real,2
TYPE,2  
ESURF,ALL   
allsel,all
  
ASEL,S,,,7  
NSLA,S,1
!ESLN,S,0
real,2
type,3
ESURF,ALL   
allsel,all
save
fini
  
/SOLU   
dk,1,ALL,0
dk,2,all,0
dk,3,all,0
dk,4,all,0   
SFA,12,1,PRES,ff
sftran
ACEL,0,0,9.8, !施加重力加速度
nlgeom,on
auto,on
nsubst,100,500,50
outres,all,all
save
solve
fini

mayuew052

你这个程序里面修改了接触刚度的值，由5改成了1，这个没什么问题，1就是默认的值吧。挺好。然后你的接触实常数也修改成了共用实常数对，这个是你修改对的。keyopt，3,10,2这个是定义接触单元的关键字，意思是接触刚度是如何更新的，选择一个更新方式。
另外提一个建议，你实体单元选择的是solid45单元，这个单元是低阶单元，你的接触单元选择的是cont174单元，这个是不匹配的，你应该选择cont173单元。
别的没有细看。感觉这个做的没有什么工程价值。

cumtkw

我现在只是用这个来进行学习，只有小的做懂做好了，实际工程的才不会出错。呵呵，我不知道，我算出来的结果合理吗？我运行后，是收敛的，结果也能查看，但是就是那结果，怎么去辨别与真实情况的符合程度，难道要看经验吗？有限元做的结果怎么能知道与实际的符合程度？。。。。迷茫呀，做完一个例子不知道自己做的对不对，这样就不知道自己是不是学会了这个接触？？？