!把钢管钢化
!计算结果无法收敛
!外面和里面为钢材,中间为橡胶
!荷载施加在节点上
!有接触存在
fini
/cle
/prep7
L1=80 !表示长度
d1=34 !表示钢管外径
d2=22 !表示橡胶内径
d3=40 !表示钢管外径
P0=500 !力N
p1=p0/(3.1412((d1/2)*(d1/2)-(d2/2)*(d2/2))) !表示相当于0.5KN时的初始压力
!p2=p0/(3.1412((d1/2)*(d1/2)-(d2/2)*(d2/2)))!表示注浆压力
p2=1000
cylind,d3/2,d1/2,10,-L1-10
cylind,d1/2,d2/2,0,-L1
cylind,d2/2,0,10,-L1-10
mu1=0.9 !设置摩擦系数
ET,1,SOLID92 !模拟钢管的弹性93模型
MP,EX,1,2.1E+005
MP,PRXY,1,0.3
ET,2,SOLID186 !使用超单元以便进行橡胶模拟
MP,EX,2,6e2 ! C0 and C1 terms for Young's modulus
MP,PRXY,2,0.4
mopt,split,on
mopt,split,err
mopt,split,warn
mopt,qmesh,main
asel,s,area,,1,16
asel,u,area,,7,8
AESIZE,ALL,4,
asel,s,area,,7,8
AESIZE,ALL,2,
VSEL,s, , , 1
vsel,a,,,3
VATT, 1, , 1, 0
MSHKEY,0
MSHAPE,1,3d
vmesh,1
vmesh,3
VSEL,s, , , 2
VATT, 2 , 0, 2
MSHKEY,0
MSHAPE,1,3d
vmesh,2
!*
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_LINECM,LINE
CM,_AREACM,AREA
/GSAV,cwz,gsav,,temp
MP,MU,2,mu1
MAT,2
R,3
REAL,3
ET,3,170
ET,4,174
RMODIF,3,1,,,1.0,0.1,,
RMODIF,3,7,,,1.0e20,0.0,1.0
KEYOPT,4,2,0
KEYOPT,4,6,0
KEYOPT,4,7,0
KEYOPT,4,8,
KEYOPT,4,9,0
KEYOPT,4,11,
KEYOPT,4,12,0
! Generate the target surface
ASEL,S,,,5
ASEL,A,,,6
ASEL,A,,,15
ASEL,A,,,16
CM,_TARGET,AREA
TYPE,3
NSLA,S,1
ESLN,S,0
ESURF,ALL
CMSEL,S,_ELEMCM
! Generate the contact surface
ASEL,S,,,9
ASEL,A,,,10
ASEL,A,,,11
ASEL,A,,,12
CM,_CONTACT,AREA
TYPE,4
NSLA,S,1
ESLN,S,0
ESURF,ALL
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,3
ESEL,A,TYPE,,4
ESEL,R,REAL,,3
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,3
ESEL,A,TYPE,,4
ESEL,R,REAL,,3
CMSEL,A,_NODECM
CMDEL,_NODECM
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM
CMSEL,S,_LINECM
CMDEL,_LINECM
CMSEL,S,_AREACM
CMDEL,_AREACM
/GRES,cwz,gsav
CMDEL,_TARGET
CMDEL,_CONTACT
fini
/solu
!分析类型
antype,static
NLGEOM,ON
!牛顿方程
nropt,full,,on
pred,on
autot,on
nsubst,10,100,5
!nsubst,10,100,5
outres,all,all
asel,s,area,,3,4
nsla,s
d,all,all,0
!F,P51X,FZ,0.1
!SFE,P51X,1,PRES, ,0.2, , ,
asel,s,area,,7
nsla,s
!*get,enum,elem,,count
*get,numna,node,,COUNT
f,all,fz,p2/numna
asel,s,area,,8
nsla,s
f,all,fz,-p2/numna-p1/numna
alls
save
solve
fini
/post1
esel,s,mat,,1
plnsol,s,1 |