24小时热门版块排行榜    

查看: 1035  |  回复: 0

杨扬弃_y

新虫 (初入文坛)

[求助] 关于钢筋混凝土梁粘结滑移的计算,请大牛帮我看一段命令流并指出问题和解决方法。

finish
/clear
/prep7  

et,1,solid65
keyopt,1,7,1
mp,ex,1,28500
mp,prxy,1,0.2
fck=30
ftk=3
tb,concr,1
tbdata,,0.5,0.95,ftk,-1
tb,miso,1,,11
tbpt,,0.0002,fck*0.19
tbpt,,0.0004,fck*0.36
tbpt,,0.0006,fck*0.51
tbpt,,0.0008,fck*0.64
tbpt,,0.001,fck*0.75
tbpt,,0.0012,fck*0.84
tbpt,,0.0014,fck*0.91
tbpt,,0.0016,fck*0.96
tbpt,,0.0018,fck*0.99
tbpt,,0.002,fck
tbpt,,0.0033,fck

et,2,link8
mp,ex,2,2e5
mp,prxy,2,0.3
tb,biso,2
tbdata,,350,0

mp,ex,3,2e5
mp,prxy,3,0.25
tb,biso,3
tbdata,,200,0

et,3,combin39
keyopt,3,3,3

et,4,solid45
mp,ex,4,2e6
mp,prxy,4,0.3

pi=acos(-1)
r,1
r,2,0.25*pi*14*14
r,3,0.25*pi*8*8
r,4,0.25*pi*6.5*6.5
r,5,0.25*pi*6.5*6.5/2
r,6,0,0,0.5E-5,4497.60,1.0E-5,7126.36
rmore,1.5E-5,8567.28,2.0E-5,9324.77,2.5E-5,9726.64
rmore,3.0E-5,9924.14,3.5E-5,9891.92,4.0E-5,9428.07
rmore,4.5E-5,8154.07,5.0E-5,5514.84
r,7


blc4,,,82.5,415,2075
/view,1,1,1,1   
/replot
wpoffs,,,2075   
*do,i,1,12,1
wpoffs,,,-170   
vsbw,all
*enddo  
wpcsys,-1   
wpoffs,,25  
wprota,,90  
vsbw,all
wpoffs,,,-365   
vsbw,all
wpcsys,-1   
wpoffs,25   
wprota,,,90
vsbw,all
wpcsys,-1   
wpoffs,,,35
blc4,,,82.5,-25,170
wpcsys,-1
/replot


vatt,1,1,1  
lsel,s,loc,x,0
lsel,r,loc,y,0
lesize,all,34
lsel,all
lsel,s,loc,z,0
lsle,r,loc,x,0
lesize,all,25
lsel,all
vsel,s,loc,y,0,415
mshkey,1
vmesh,all   
allsel,all  

vatt,4,7,4  
lsel,s,loc,x,0
lsel,r,loc,y,-25,0
lsel,r,loc,z,35
lesize,all,25
lsel,all
vsel,s,loc,y,-25,0
mshkey,1
vmesh,all   
nummrg,all  



lsel,s,tan1,z   
lsel,r,loc,x,25,82.5
lsel,r,loc,y,25,390
lsel,u,loc,x,82.5
lsel,u,loc,z,0  
lsel,u,loc,z,125
lsel,u,loc,z,1425
lsel,u,loc,z,1525   
lsel,u,loc,z,2075
latt,3,4,2
lesize,all,25
type,2  
real,4  
mat,3   
lmesh,all  
lsel,all


lsel,s,loc,z,2075   
lsel,r,loc,x,25,82.5
lsel,r,loc,y,25,390
lsel,u,loc,x,82.5  
latt,3,5,2  
lesize,all,25
type,2  
real,5  
mat,3   
lmesh,all
lsel,all

lsel,s,loc,x,25
lsel,r,loc,y,390
lsel,r,loc,z,35,2075
latt,3,3,2
lesize,all,25
type,2  
real,3  
mat,3   
lmesh,all   
lsel,all

esel,s,type,,2
/eshape,1   
eplot  
allsel,all
eplot

lsel,none
k,10000,25,25,35   
k,10001,25,25,2075  
l,10000,10001
lesize,all,34  
latt,2,2,2  
lmesh,all

cpintf,ux !创建节点的耦合自由度UX
*get,cpmax,cp,,max !获取藕合组总数
type,3
real,6
*do,i,1,cpmax
*get,n1,cp,i,term,1,node !获取第i组的第一个节点号
*get,n2,cp,i,term,2,node !获取第i组的第二个节点号
e,n1,n2
*enddo
cpintf,uy

asel,s,loc,y,-25
asel,r,loc,z,35,205
da,all,uy
allsel,all

asel,s,loc,x,82.5   
da,all,symm
asel,s,loc,z,2075   
da,all,symm
allsel,all

*SET,p0,60e3
*SET,q0,p0/82.5/170
asel,s,loc,z,1395,1565  
asel,r,loc,y,415
sfa,all,1,pres,q0   
allsel,all

/solu
antype,0
outres,all,all  
autots,on   
neqit,60
nsubst,800,,500
cnvtol,u,,0.05  
cnvtol,f,,0.05
solve
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 杨扬弃_y 的主题更新
信息提示
请填处理意见