找回密码
 注册
Simdroid-非首页
查看: 485|回复: 12

[接触分析] 一个螺栓连接的命令流

[复制链接]
发表于 2007-9-24 17:04:49 | 显示全部楼层 |阅读模式 来自 上海松江区

这是我看到的一个命令流,希望对大家有点用,谢谢。
!ANSYS6.1计算test3(对称模型,只建一半模型,试验前分析)
!/BATCH  
!/input,menust,tmp     ,'',,,,,,,,,,,,,,,,1  
/PREP7
!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!
参数定义!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!   
hc=400 !
柱截面高度
bc=200 !
柱截面宽度
tcf=10 !
柱翼缘厚度
tcw=8  !
柱腹板厚度
lc=1100 !
柱构件伸出长度
tep=20 !
端板厚度
bep=bc+20  !
端板宽度
HEP1= hc+200   !
端板高度
tst=10    !
端板外伸部分加劲肋厚度
hst=80    !
端板外伸部分加劲肋高度
bst=bc/2-5   !
端板外伸部分加劲肋高度
lbt=2*tep     !
螺栓杆长度
dbt=20      !
螺栓杆直径或有效直径
dbth=31.4     !
螺栓头和螺母直径
lbth=12.5    !
螺栓头厚度
preten=155000   !
螺栓施工预拉力  
miu=.4   !
端板间抗滑移系数
hb=400    !
梁截面高度   
bb=200    !
梁截面宽度   
tbf=8      !
梁翼缘厚度  
tbw=6     !
梁腹板厚度   
lb1=870    !
梁构件伸出长度  
lb2=200     !
梁构件伸出的水平加载断长度
hb1=362       !
梁最左端高度
dh0=dbt+2     !
螺栓孔直径
randa=0.05      !
梁的坡度
aa=50      !
螺栓中心到梁翼缘边缘(非受力方向)的距离
aa1=50      !
螺栓中心到梁翼缘边缘(受力方向的距离
ab=120      !
一二排螺栓间距
displa=-50    !
施加的位移荷载大小   

                                                                        
!!!!!!!!!!!
建模开始!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!
梁模型!!!!!!!!!!
block,o,tep, -hep1/2,hep1/2,0,bep/2
block,-tep,0,-hep1/2-100+4*tcf,hep1/2,0,bep/2
k,20,tep,hb/2
k,21,tep,hb/2-tbf
k,22,tep+lb1,hb/2-tbf+lb1*randa
k,23,tep+lb1,hb/2+lb1*randa
k,24,tep,hb/2,bb/2
k,25,tep,hb/2-tbf,bb/2
k,26,tep+lb1,hb/2-tbf+lb1*randa,bb/2
k,27,tep+lb1,hb/2+lb1*randa,bb/2
v,20,21,22,23,24,25,26,27
k,28,tep,hb/2-tbf
k,29,tep,-hb/2+tbf
k,30,tep+lb1,-hb/2+tbf+lb1*randa
k,31,tep+lb1,hb/2-tbf+lb1*randa
k,32,tep,hb/2-tbf,tbw/2                  
k,33,tep,-hb/2+tbf ,tbw/2                 
k,34,tep+lb1,-hb/2+tbf+lb1*randa ,tbw/2   
k,35,tep+lb1,hb/2-tbf+lb1*randa,tbw/2   
v,28,29,30,31,32,33,34,35

block,tep+lb1,tep+lb1+lb2,hb/2-tbf+lb1*randa,hb/2+lb1*randa,0,bb/2
block,tep+lb1,tep+lb1+lb2,-hb/2+tbf+lb1*randa,hb/2-tbf+lb1*randa,0,tbw/2
vgen,2,3,5,2,,-hb+tbf,, ,0

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!
颜色区分volumn!!!!!!!
/PNUM,KP,0  
/PNUM,LINE,0
/PNUM,AREA,0
/PNUM,VOLU,1
/PNUM,NODE,0
/PNUM,TABN,0
/PNUM,SVAL,0
/NUMBER,1   
/PNUM,ELEM,0
/REPLOT

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!
柱模型!!!!!!!!!!
!!!!!!!
端板对接斜坡!!!!!!!!!!!!!!!!!!!!
!!!!!!
端板对接斜坡!!!!!!!!!!!!!!!!!!!!
k,80,-tep,-hep1/2-100+4*tcf
k,81,-tep,-hep1/2-100
k,82,-tep+tcf, -hep1/2-100
k,83,0,-hep1/2-100+4*tcf  
k,84,-tep,-hep1/2-100+4*tcf,bep/2      
k,85,-tep,-hep1/2-100,bep/2              
k,86,-tep+tcf, -hep1/2-100,bep/2         
k,87,0,-hep1/2-100+4*tcf,bep/2   
v,80,81,82,83,84,85,86,87

!!!!!!
柱翼缘!!!!!!!!!!!!!!!!!!!!
block,-tep,-tep+tcf,-hep1/2-100,-hb/2-lc,0,bc/2
block,-tep-hc+tcf,-tep-hc+2*tcf,-hb/2-lc,hb/2-(hc-2*tcf)*randa,0,bc/2         
!!!!!!
柱腹板!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!
k,89,-tep-hc+2*tcf,-hb/2-lc
k,90,-tep,-hb/2-lc
k,92,-tep-hc+2*tcf,hb/2-(hc-2*tcf)*randa-tbf   
k,91,-tep,hb/2-tbf
k,93,-tep-hc+2*tcf,-hb/2-lc,tcw/2                                      
k,94,-tep,-hb/2-lc,tcw/2                             
k,96,-tep-hc+2*tcf,hb/2-(hc-2*tcf)*randa-tbf,tcw/2   
k,95,-tep,hb/2-tbf,tcw/2
v,89,90,91,92,93,94,95,96                        
k,98, -tep,hb/2
k,97, -tep-hc+2*tcf,hb/2-(hc-2*tcf)*randa
k,100, -tep,hb/2,bb/2                              
k,99, -tep-hc+2*tcf,hb/2-(hc-2*tcf)*randa,bb/2  
v,91,92,97,98,95,96,99,100
block,-tep-hc+2*tcf,-tep,-hb/2,-hb/2+tst,tcw/2,tcw/2+bst  
     
!!!!!!
端板外伸加劲肋!!!!!!!!!!!!!!!!!!!!
k,109,tep,hb/2
k,110,tep+hst,hb/2+hst*randa
k,111,tep,hb/2+hst
k,112,tep,hb/2,tst/2               
k,113,tep+hst,hb/2+hst*randa,tst/2  
k,114,tep,hb/2+hst,tst/2            
V, 109,110,111,112,113,114
k,115,tep,-hb/2                       
k,116,tep+hst,-hb/2+hst*randa         
k,117,tep,-hb/2-hst                  
k,118,tep,-hb/2,tst/2                     
k,119,tep+hst,-hb/2+hst*randa ,tst/2      
k,120,tep,-hb/2-hst,tst/2                  
V, 115,116,117,118,119,120
k,121,-tep,hb/2                          
k,122,-tep-hst,hb/2-hst*randa                     
k,123,-tep,hb/2+hst                     
k,124,-tep,hb/2,tst/2                  
k,125,-tep-hst,hb/2-hst*randa ,tst/2   
k,126,-tep,hb/2+hst ,tst/2              
v, 121,122,123,124,125,126            

!!!!!!
梁加劲肋!!!!!!!!!!!!!!!!!!!!
wpoff,tep+lb1+lb2/2,-hb/2+tbf+lb1*randa,0
block,-tst/2,tst/2,0,hb-2*tbf,tbw/2,tbw/2+bst
block,75-tst/2,tst/2+75,0,hb-2*tbf,tbw/2,tbw/2+bst
block,-75-tst/2,-75+tst/2,0,hb-2*tbf,tbw/2,tbw/2+bst
wplane,,0,hb/2+100-aa,bep/2-aa1,0,hb/2+aa+10,bep/2-aa1,0,hb/2-aa+100+10,bep/2-aa1+10
!!!!!!
端板打孔!!!!!!!!!!!!!!!!!!!
cylind,dh0/2,,-tep-10,tep+10,0,360
vsel,S,LOC,Y,hb/2+100-aa  
vgen,2,all, , ,,-ab , , ,0
vgen,2,all, , ,,-(hep1-2*aa-ab) , , ,0
cm,bolt,volu


!!!!
布尔运算!!!!!!!
vsel,s,loc,x,-tep+5,tep   
VSEL,U, , ,bolt
cm,endplate,volu
cmsel,all
VSBV,endplate,bolt,keep

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!
建螺栓!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

cylind,dbt/2,,-tep,tep,0,360
cylind,dbth/2,,-tep-lbth,-tep,0,360
cylind,dbth/2,,tep,tep+lbth,0,360
wpoff,-ab,0
cylind,dbt/2,,-tep,tep,0,360
cylind,dbth/2,,-tep-lbth,-tep,0,360
cylind,dbth/2,,tep,tep+lbth,0,360
wpoff,-(hep1-2*aa-ab),0
cylind,dbt/2,,-tep,tep,0,360
cylind,dbth/2,,-tep-lbth,-tep,0,360
cylind,dbth/2,,tep,tep+lbth,0,360
wpoff,ab
cylind,dbt/2,,-tep,tep,0,360
cylind,dbth/2,,-tep-lbth,-tep,0,360
cylind,dbth/2,,tep,tep+lbth,0,360
allsel,all
NUMMRG,ALL, , , ,LOW

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!
粘贴柱腹板翼缘端板和相应的加劲肋!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
vsel,s,loc,x,-tep-hc,0  
vsel,u,loc,x,0,tep
VGLUE,all  
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!
粘贴梁腹板翼缘端板和相应的加劲肋!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
vsel,s,loc,x,0,lb1+lb2+tep
vsel,u,loc,x,-tep,0
VGLUE,all

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!
粘贴螺栓!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
vsel,s,loc,z,bep/2-aa1
VGLUE,all
cm,bolt,volu
ALLSEL,ALL  
NUMMRG,ALL, , , ,LOW

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!
定义材料和单元属性!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!
下面定义材料1Q345钢材用于梁柱端板!!!!!!!!!
!!!!!!!!!!!!!!!!
et,1,solid92
mp,ex,1,206e3
mp,prxy,1,0.3
tb,biso,1,1,2
tbtemp,0
tbdata,,345,0,,,,

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!
下面定义材料2,用于高强螺栓!!!!!!!!!
!!!!!!!!!!!!!!!!
mp,ex,2,206e3
mp,prxy,2,0.3
tb,biso,2,1,2
tbtemp,0
tbdata,,940,0,,,,
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!
给单元附材料号和属性!!!!!!!!!
!!!!!!!!!!!!!!!!
!!!!
高强螺栓附属性!!!!!!!!
cmsel,s,bolt
VATT, 2,,1,0   
!!!!!!!!!!!!!!!!
!!!!
梁柱端板附属性!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!
ALLSEL,ALL
cmsel,u,bolt
VATT, 1,,1,0  

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!
网格划分!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!mesh
螺栓杆!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!
cmsel,s,bolt
vsel,r,loc,x,0
LSEL,s,radius,,dbt/2
LESIZE,all, , ,4, , , , ,1  
LSEL,s,length,,2*tep
LESIZE,all, , ,10, , , , ,1  
MSHAPE,1,3D
MSHKEY,0
VMESH,all

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!mesh
螺栓头和螺母!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!
cmsel,s,bolt
vsel,u,loc,x,0
LSEL,s,radius,,dbth/2
LESIZE,all, , ,4, , , , ,1   
LSEL,s,radius,,dh0/2
LESIZE,all, , ,8, , , , ,1
LSEL,s,length,,lbth
LESIZE,all, , ,2, , , , ,1
VMESH,all
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!mesh
端板!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!
先建endplatecomponent!!!!!!!!!!!!!!!!!!!
allsel,all
vsel,s,loc,x,-tep+tcf/2,tep
vsel,u,loc,x,0
cm,endplate,volu
ESIZE,bep/10
VMESH,all
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!mesh
!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
vsel,s,loc,x,-tep-hc,0
vsel,r,loc,z,0,tcf
vsel,u,loc,y,hb/2,hb/2+hst
cm,col_web,volu
eSIZE,hc/8
VMESH,all
vsel,s,loc,x,-tep-hc,0
cmsel,u,bolt
cmsel,u,endplate
cmsel,u,col_web
eSIZE,bc/6
VMESH,all


!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!mesh
!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
vsel,s,loc,x,tep,tep+lb1+lb2
cmsel,u,bolt
lsel,s,loc,x,tep+lb1+lb2
LSEL,r,length,,bb/2
LESIZE,all,bb/6
lsel,s,loc,x,tep+lb1+lb2
LSEL,r,length,,hb-2*tbf
LESIZE,all,hb/8
VMESH,all

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!
接触建立!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
allsel,all
cmsel,s,endplate
ASLV,S
asel,r,loc,x
asel,u,loc,y,-hep/2,-hep/2-10
cm,target1,area
allsel,all
cmsel,s,endplate
ASLV,S
asel,r,loc,x
asel,r,loc,y,-hep/2,-hep/2-10
cm,contact1,area
cmsel,s,target1
cmsel,a,contact1

MP,MU,1,0.4
MAT,1  
R,3
REAL,3  
ET,3,170
ET,4,174
KEYOPT,4,9,0

ASEL,S,,,target1
TYPE,3  
NSLA,S,1
ESLN,S,0
ESURF,ALL   

ASEL,S,,,contact1
TYPE,4  
NSLA,S,1
ESLN,S,0
ESURF,ALL   
ALLSEL  
!

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!
对高强螺栓psmesh!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
cmsel,s,bolt
vsel,r,loc,x,
cm,boltshank1,volu
PSMESH,1,bolt1, ,v,1, 0,x,
PSMESH,2,bolt2, ,V,22, 0,x,
PSMESH,3,bolt3, ,V,27, 0,x,
PSMESH,4,bolt4, ,V,30, 0,x,

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!
下面施加荷载和位移条件进行求解
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!
下面施加loadstep 1,边界约束条件、柱轴向压力和螺栓预拉力
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
ALLSEL,ALL  
FINISH  
/SOLU   
ANTYPE,0
NLGEOM,1
!!NSUBST,10,0,0   
NSUBST,12,0,0   



!!NEQIT,200

NEQIT,400  
CNVTOL,F, ,0.0 1,2, ,


ASEL,S,LOC,Z,   
DA,all,SYMM

nsel,s,loc,y,-(hb/2+lc+1)
D,all, ALL

ALLSEL,ALL
SLOAD,ALL,9,LOCK,FORC,155000, 1,2

FINISH  
/SOLU   
LSWRITE,1,
solve
save
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!
下面施加loadstep 2,粱端荷载
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
ALLSEL,ALL

FINISH  
/SOLU   
ANTYPE,,REST,,,0  
!*  
NLGEOM,1
NSUBST,20,5000,1
OUTRES,alL,1

ALLSEL,ALL  
num=node(tep+lb1+lb2/2,hb/2+(lb1+lb2/2)*randa,0)   
D,num, uy,displa  
LSWRITE,2,  
ALLSEL,ALL
EPLOT
SAVE
solve
save
发表于 2007-11-12 14:40:58 | 显示全部楼层 来自 江苏无锡
Simdroid开发平台
你不能把别人的命令流放在这个地方,别人还有问题呢?你解决了吗?容易误导。
回复 不支持

使用道具 举报

发表于 2008-3-7 09:47:30 | 显示全部楼层 来自 湖南湘潭
有作螺栓连接分析,螺纹分析的就加这个群,我也是研究螺栓的
群号是42723503
回复 不支持

使用道具 举报

发表于 2009-5-7 22:22:01 | 显示全部楼层 来自 江苏无锡
好东西~多谢~
回复 不支持

使用道具 举报

发表于 2009-5-21 13:07:30 | 显示全部楼层 来自 浙江杭州
好东西,楼主辛苦
回复 不支持

使用道具 举报

发表于 2009-10-30 10:49:56 | 显示全部楼层 来自 陕西西安
谢谢啦!辛苦啦!
回复 不支持

使用道具 举报

发表于 2009-11-2 02:31:43 | 显示全部楼层 来自 广东广州
真诚的微笑,献给好帖子
否则就要责怪
回复 不支持

使用道具 举报

发表于 2012-10-14 15:04:34 | 显示全部楼层 来自 湖南长沙
感谢分享@
回复 不支持

使用道具 举报

发表于 2012-11-7 22:40:18 | 显示全部楼层 来自 陕西西安
好例子,值得分享
回复 不支持

使用道具 举报

发表于 2013-1-14 15:33:29 | 显示全部楼层 来自 美国
XUE XI LE           
回复 不支持

使用道具 举报

发表于 2013-1-14 21:56:54 | 显示全部楼层 来自 天津
还没试,明天具体试了再来提问题~  辛苦!我正在学习关于螺栓连接的建模!
回复 不支持

使用道具 举报

发表于 2013-1-20 10:30:00 | 显示全部楼层 来自 湖北武汉
要早看到这贴就好了   前天刚刚花了好长时间把这个例子从pdf文件打成了txt格式
回复 不支持

使用道具 举报

发表于 2013-2-3 14:28:03 | 显示全部楼层 来自 广西柳州
ok!!好啊!
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-5-23 19:41 , Processed in 0.079587 second(s), 13 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表