约束了啊,为什么会出现这种报错?以下是命令流文件:
!!!!!环境设置!!!!!!!!!!!!!
/filname,bolt
/title,flang_bolt
!!!!设置环境基本参数,单位N,mm,MPa,密度tons/mm3
*afun,deg
*set,bolt_rad,8/2 !螺栓半径
*set,bolt_len,12 !螺栓长度
*set,bolt_rad2,13/(2*sin(60)) !螺栓头半径
*set,bolt_h,5.5 !螺栓头高度
*set,nut_h,3 !螺母高度
*set,nut_rad, 6 !螺母半径
*set,flange_rad1,8 !法兰1半径
*set,flange_rad2,9 !法兰2半径
*set,flange_h1,4 !法兰1厚度
*set,flange_h2,4 !法兰2厚度
*set,forc_bolt,35000 !螺栓的预紧力
/prep7
et,1,solid45 !定义法兰和螺栓的单元为8节点6面体实体单元
et,2,targe170 !定义3d接触的目标面单元
et,3,conta174 !定义3d接触的接触面单元
keyopt,3,9,0 !接触包含初始化穿透
!!!!!!!!!!!!!!!定义材料
!!!!!定义法兰材料
mp,ex,1,2.16e5
mp,nuxy,1,.3
mp,dens,1,7.85e-9
!!!!定义螺栓螺母材料
mp,ex,2,2.11e5
mp,nuxy,2,.3
mp,dens,2,7.85e-9
!!!定义摩擦系数
mp,mu,3,.3
!!!!!!!创建下法兰
cyl4,,,bolt_rad,0,flange_rad1,360,flange_h1
vsel,all
cm,flange_down,volu
vsel,all
!!!!!!创建上法兰
wpoffs,,,flange_h1
cyl4,,,bolt_rad,0,flange_rad2,360,flange_h2
wpave,0,0,0
vsel,all
vsel,u,volu,,flange_down
cm,flange_up,volu
vsel,all
!!!!创建螺栓
cyl4,,,bolt_rad,0,,,bolt_len+1
vsel,all
vsel,u,volu,,flange_down
vsel,u,volu,,flange_up
cm,bolt1,volu
wpave,0,0,0
rpr4,6,,,bolt_rad2,,-bolt_h
vsel,all
vsel,u,volu,,flange_down
vsel,u,volu,,flange_up
vsel,u,volu,,bolt1
cm,bolt2,volu
vsel,all
vsel,s,volu,,bolt1
vsel,a,volu,,bolt2
cm,bolt,volu
vsel,all
!!!!!创建螺母
wpoffs,,,flange_h1+flange_h2
rpr4,6,,,nut_rad,,nut_h
vsel,all
vsel,u,volu,,flange_down
vsel,u,volu,,flange_up
vsel,u,volu,,bolt
cm,nut1,volu
vsel,all
vsbv,nut1,bolt,sepo,delete,keep
vsel,all
vsel,u,volu,,flange_down
vsel,u,volu,,flange_up
vsel,u,volu,,bolt
cm,nut,volu
!!!!!!!!!!!划分网格!!!!!!!!!!!
mat,1
esize,1
cmsel,s,flange_up
vsweep,all
esize,1
cmsel,s,flange_down
vsweep,all
esize,.5
cmsel,s,nut
vsweep,all
esize,.4
cmsel,s,bolt1
vsweep,all
esize,.5
cmsel,s,bolt2
vsweep,all
vsel,all
cmsel,s,bolt
nslv,s,1
nummrg,node
allsel
!!!!!网格模型划分完毕!!!!!!!!!
!!!!!!!!建立预紧单元!!!!!!!!!!!!
wpave,
!wpoffs,,,flange_h1
vsel,s,volu,,bolt1
psmesh,1,bolt_el,,all,,0,z,3,,,,pre_elem !建立预紧单元
allsel
wpave,
!!!!!!!!预紧单元建立完毕!!!!!!!!
!!!!!!!创建接触!!!!!!!!!!!!!!!
!!!创建螺栓和下法兰接触
cmsel,s,bolt2
aslv
asel,r,loc,z,0
cm,blo_ar,area
allsel
cmsel,s,flange_down
aslv
asel,r,loc,z,0
cm,fld_1,area
r,1
real,1
mat,3
cmsel,s,blo_ar
type,2
nsla,s,1
esln,s,0
esurf,all
allsel
cmsel,s,fld_1
type,3
nsla,s,1
esln,s,0
esurf,all
allsel
!!!!!!!!!!!!创建法兰和法兰接触
et,5,targe170
et,6,conta174
keyopt,6,9,0
r,2
real,2
mat,3
type,5
cmsel,s,flange_down
aslv
asel,r,loc,z,flange_h1
nsla,s,1
esln,s,0
esurf,all
type,6
cmsel,s,flange_up
aslv
asel,r,loc,z,flange_h1
nsla,s,1
esln,s,0
esurf,all
allsel
!!!!!!!!!!!创建上法兰和螺母接触
et,7,targe170
et,8,conta174
r,3
real,3
mat,3
type,8
cmsel,s,flange_up
aslv
asel,r,loc,z,flange_h1+flange_h2
nsla,s,1
esln,s,0
esurf,all
type,7
cmsel,s,nut
aslv
asel,r,loc,z,flange_h1+flange_h2
nsla,s,1
esln,s,0
esurf,all
allsel
finish
!!!!!!!!!接触创建完毕
/solu
allsel
lsclear,all
nropt,full
pred,on
lnsrch,on
antype,0
autots,1
nsubst,25,100,1
outres,all,all
eqslv,pcg,1e-3
sload,1,9,lock,forc,forc_bolt,1,2 !施加预紧力,第二步载荷锁住位移
finish
/prep7
csys,1
nrotat,all
cmsel,s,flange_down
aslv,s
nsla,s,1
nsel,r,loc,x,flange_rad1
d,all,all
allsel,all
csys,0
nrotat,all
allsel
finish
/solu
time,10
allsel
lswrite,1
time,20
lswrite,2
lssolve,1,2,1
finish
/post1
set,last
plnsol,s,eqv,0,1
cmsel,s,bolt1
eslv,s,1
nslv,s,1
plnsol,s,eqv,0,1
allsel
cmsel,s,nut
eslv,s,1
nslv,s,1
plnsol,s,eqv,0,1
[ 本帖最后由 cax 于 2008-12-26 01:05 编辑 ] |