切换到宽版

 找回密码
 注册
Simdroid-非首页
Simwe仿真论坛(forum.simwe.com)，CAE/CAD/CAM/，FEA/FEM/有限元分析论坛---(手机验证注册)»论坛 › ..:: ANSYS专区 › C08：ANSYS--其他软件 › C型钢半刚性节点的静力非线性求解
返回列表 发新帖
查看: 115|回复: 2

[新手疑问区] C型钢半刚性节点的静力非线性求解

[复制链接]
army011
发表于 2006-10-23 15:07:54 | 显示全部楼层 |阅读模式 来自 内蒙古包头
本人做ansys半年余，求解总是不收敛，看了好长时间帖子，还是调不出来，命令流如下，请教各位大侠。

*SET,CA,150 !C型钢参数
*SET,CB,65
*SET,CC,20
*SET,CD,1.8
*SET,CR,6
*SET,BL,35  !螺栓端距
*SET,BR1,8  !螺栓半径
*SET,BR2,14 !垫片外半径
*SET,BR3,12 !螺帽半径
*SET,BR4,18 !接触区域半径
*SET,BR5,9  !螺栓孔半径
*SET,BR6,10 !垫片内半径
*SET,BH1,2  !垫片高度
*SET,BH2,6  !螺帽高度
*SET,BJ,1.2 !螺栓间距参数1.2;1.4;1.6
*SET,JPT,6  !节点板厚度
*SET,SD,10  !安装间隙
*SET,JJ,600 !支撑间距

!加载参数
*SET,PREFOR,70000!施加螺栓预应力
*SET,DISP,100    !施加位移


/PREP7  
ET,1,solid45

MPTEMP,,,,,,,,        !节点板
MPTEMP,1,0  
mp,ex,1,176000  
mp,prxy,1,0.3   
tb,kinh,1,1,4   
tbpt,,235/(176000),235  
tbpt,,0.002,235
tbpt,,0.02,235  
tbpt,,0.1,235/2

MPTEMP,,,,,,,,        !C型钢  
MPTEMP,1,0  
mp,ex,2,176000  
mp,prxy,2,0.3   
tb,kinh,2,1,4   
tbpt,,235/(176000),235  
tbpt,,0.002,235
tbpt,,0.02,235  
tbpt,,0.1,235/2

MPTEMP,,,,,,,,        !弯角  
MPTEMP,1,0  
mp,ex,3,176000  
mp,prxy,3,0.3   
tb,kinh,3,1,4   
tbpt,,235*1.2/(176000),235*1.2  
tbpt,,0.1,235*1.2   
tbpt,,0.02,235*1.2  
tbpt,,0.1,235*1.2/2

MPTEMP,,,,,,,,        !高强螺栓  
MPTEMP,1,0  
MPDATA,EX,4,,206000
MPDATA,PRXY,4,0.3,  
TB,BKIN,4,1,2,1
TBDATA,,640,20600,,,,   

MPTEMP,,,,,,,,      
MPTEMP,1,0  
mp,ex,5,176000  
mp,prxy,5,0.3   
tb,kinh,5,1,4   
tbpt,,235/(176000),235  
tbpt,,0.002,235
tbpt,,0.02,235  
tbpt,,0.1,235/2
MP,MU,5,0.25          !节点板与型钢  

MPTEMP,,,,,,,,        
MPTEMP,1,0  
mp,ex,6,176000  
mp,prxy,6,0.3   
tb,kinh,6,1,4   
tbpt,,235/(176000),235  
tbpt,,0.002,235
tbpt,,0.02,235  
tbpt,,0.1,235/2
MP,MU,6,0.15          !螺栓与螺孔

  
!柱子建模
CYL4,-CA/2+CR,-JPT/2-CR,CR,90,CR-CD,180,JJ-BL
CYL4,-CA/2+CR,-JPT/2-CB+CR,CR,180,CR-CD,270,JJ-BL
CYL4,CA/2-CR,-JPT/2-CR,CR,0,CR-CD,90,JJ-BL
CYL4,CA/2-CR,-JPT/2-CB+CR,CR,270,CR-CD,360,JJ-BL

BLOCK,-CA/2+CR,CA/2-CR,-JPT/2,-JPT/2-CD,0,JJ-BL
BLOCK,-CA/2+CR,CA/2-CR,-JPT/2,-JPT/2-CD,JJ-BL,JJ+BJ*CA+BL
BLOCK,-CA/2+CR,CA/2-CR,-JPT/2,-JPT/2-CD,JJ+BJ*CA+BL,JJ+BJ*CA+JJ   

BLOCK,-CA/2,-CA/2+CD,-JPT/2-CR,-JPT/2-CB+CR,0,JJ-BL
BLOCK,-CA/2,-CA/2+CD,-JPT/2-CR,-JPT/2-CB+CR,JJ-BL,JJ+BJ*CA+BL
BLOCK,-CA/2,-CA/2+CD,-JPT/2-CR,-JPT/2-CB+CR,JJ+BJ*CA+BL,JJ+BJ*CA+JJ

BLOCK,CA/2,CA/2-CD,-JPT/2-CR,-JPT/2-CB+CR,0,JJ-BL
BLOCK,CA/2,CA/2-CD,-JPT/2-CR,-JPT/2-CB+CR,JJ-BL,JJ+BJ*CA+BL
BLOCK,CA/2,CA/2-CD,-JPT/2-CR,-JPT/2-CB+CR,JJ+BJ*CA+BL,JJ+BJ*CA+JJ

BLOCK,-CA/2+CR,-CA/2+CC,-JPT/2-CB,-JPT/2-CB+CD,0,JJ-BL
BLOCK,-CA/2+CR,-CA/2+CC,-JPT/2-CB,-JPT/2-CB+CD,JJ-BL,JJ+BJ*CA+BL
BLOCK,-CA/2+CR,-CA/2+CC,-JPT/2-CB,-JPT/2-CB+CD,JJ+BJ*CA+BL,JJ+BJ*CA+JJ

BLOCK,CA/2-CR,CA/2-CC,-JPT/2-CB,-JPT/2-CB+CD,0,JJ-BL
BLOCK,CA/2-CR,CA/2-CC,-JPT/2-CB,-JPT/2-CB+CD,JJ-BL,JJ+BJ*CA+BL
BLOCK,CA/2-CR,CA/2-CC,-JPT/2-CB,-JPT/2-CB+CD,JJ+BJ*CA+BL,JJ+BJ*CA+JJ

wpoff,0,-3,JJ+BJ*CA/2   
wpro,,90.000000,
wpro,,,-90.000000   
wpro,90.000000,,
wpro,,90.000000,
CYL4,BJ*CA/2,CA/2-BL,BR4, , , ,CD
CYL4,BJ*CA/2,-(CA/2-BL),BR4, , , ,CD
CYL4,-BJ*CA/2,CA/2-BL,BR4, , , ,CD
CYL4,-BJ*CA/2,-(CA/2-BL),BR4, , , ,CD   
FLST,2,5,6,ORDE,3   
FITEM,2,6   
FITEM,2,20  
FITEM,2,-23
VOVLAP,P51X
CYL4,BJ*CA/2,CA/2-BL,BR5, , , ,6
CYL4,BJ*CA/2,-(CA/2-BL),BR5, , , ,6   
CYL4,-BJ*CA/2,CA/2-BL,BR5, , , ,6   
CYL4,-BJ*CA/2,-(CA/2-BL),BR5, , , ,6  
VSBV,      21,      25  
VSBV,      23,      27  
VSBV,      22,      26  
VSBV,      20,       6  
wpro,,-90.000000,
wpro,-90.000000,,
wpro,,,90.000000   
wpro,,-90.000000,
wpoff,0,JPT/2,-BJ*CA/2-BL   
CYL4,-CA/2+CR,-JPT/2-CR,CR,90,CR-CD,180,BJ*CA+BL*2
CYL4,-CA/2+CR,-JPT/2-CB+CR,CR,180,CR-CD,270,BJ*CA+BL*2
CYL4,CA/2-CR,-JPT/2-CR,CR,0,CR-CD,90,BJ*CA+BL*2
CYL4,CA/2-CR,-JPT/2-CB+CR,CR,270,CR-CD,360,BJ*CA+BL*2

wpoff,0,0,BJ*CA+BL*2
CYL4,-CA/2+CR,-JPT/2-CR,CR,90,CR-CD,180,JJ-BL
CYL4,-CA/2+CR,-JPT/2-CB+CR,CR,180,CR-CD,270,JJ-BL
CYL4,CA/2-CR,-JPT/2-CR,CR,0,CR-CD,90,JJ-BL
CYL4,CA/2-CR,-JPT/2-CB+CR,CR,270,CR-CD,360,JJ-BL
ALLSEL,ALL
VGLUE,all
CM,column,VOLU


!梁建模

wpoff,0,0,-CA/2
wpro,,90.000000,
wpoff,CA/2+SD,0,0  
wpro,,,-90.000000   
wpro,90.000000,,
CYL4,-CA/2+CR,-JPT/2-CR,CR,90,CR-CD,180,-(BJ*CA+BL*2)  
CYL4,-CA/2+CR,-JPT/2-CB+CR,CR,180,CR-CD,270,-(BJ*CA+BL*2)  
CYL4,CA/2-CR,-JPT/2-CR,CR,0,CR-CD,90,-(BJ*CA+BL*2)
CYL4,CA/2-CR,-JPT/2-CB+CR,CR,270,CR-CD,360,-(BJ*CA+BL*2)   

BLOCK,-CA/2+CR,CA/2-CR,-JPT/2,-JPT/2-CD,0,-(BL+BJ*CA+BL)   
BLOCK,-CA/2+CR,CA/2-CR,-JPT/2,-JPT/2-CD,-(BL+BJ*CA+BL),-(JJ+BJ*CA+BL)   
BLOCK,-CA/2,-CA/2+CD,-JPT/2-CR,-JPT/2-CB+CR,0,-(BL+BJ*CA+BL)   
BLOCK,-CA/2,-CA/2+CD,-JPT/2-CR,-JPT/2-CB+CR,-(BL+BJ*CA+BL),-(JJ+BJ*CA+BL)

BLOCK,CA/2,CA/2-CD,-JPT/2-CR,-JPT/2-CB+CR,0,-(BL+BJ*CA+BL)
BLOCK,CA/2,CA/2-CD,-JPT/2-CR,-JPT/2-CB+CR,-(BL+BJ*CA+BL),-(JJ+BJ*CA+BL)
BLOCK,-CA/2+CR,-CA/2+CC,-JPT/2-CB,-JPT/2-CB+CD,0,-(BL+BJ*CA+BL)
BLOCK,-CA/2+CR,-CA/2+CC,-JPT/2-CB,-JPT/2-CB+CD,-(BL+BJ*CA+BL),-(JJ+BJ*CA+BL)
BLOCK,CA/2-CR,CA/2-CC,-JPT/2-CB,-JPT/2-CB+CD,0,-(BL+BJ*CA+BL)
BLOCK,CA/2-CR,CA/2-CC,-JPT/2-CB,-JPT/2-CB+CD,-(BL+BJ*CA+BL),-(JJ+BJ*CA+BL)
wpoff,0,0,-(BJ*CA+BL*2)
CYL4,-CA/2+CR,-JPT/2-CR,CR,90,CR-CD,180,-(JJ-BL)  
CYL4,-CA/2+CR,-JPT/2-CB+CR,CR,180,CR-CD,270,-(JJ-BL)  
CYL4,CA/2-CR,-JPT/2-CR,CR,0,CR-CD,90,-(JJ-BL)
CYL4,CA/2-CR,-JPT/2-CB+CR,CR,270,CR-CD,360,-(JJ-BL)

wpoff,0,0,(BJ*CA+BL*2)
wpoff,0,-JPT/2,-(BL+BJ*CA/2)
wpro,,90.000000,
CYL4,CA/2-BL,BJ*CA/2,BR4, , , ,CD
CYL4,-(CA/2-BL),BJ*CA/2,BR4, , , ,CD
CYL4,CA/2-BL,-BJ*CA/2,BR4, , , ,CD
CYL4,-(CA/2-BL),-BJ*CA/2,BR4, , , ,CD   

FLST,2,5,6,ORDE,5   
FITEM,2,9   
FITEM,2,26  
FITEM,2,-27
FITEM,2,29  
FITEM,2,-30
VOVLAP,P51X

CYL4,CA/2-BL,BJ*CA/2,BR5, , , ,6
CYL4,-(CA/2-BL),BJ*CA/2,BR5, , , ,6   
CYL4,CA/2-BL,-BJ*CA/2,BR5, , , ,6   
CYL4,-(CA/2-BL),-BJ*CA/2,BR5, , , ,6  
VSBV,      26,       9  
VSBV,      27,      55  
VSBV,      30,      57  
VSBV,      29,      56   
ALLSEL,ALL
CMSEL,U,COLUMN
VGLUE,ALL
CM,beam,VOLU

!节点板建模

BLOCK,-CA/2,CA/2,-BL-BJ*CA/2,BL+BJ*CA/2+SD,0,-JPT/2
BLOCK,-BJ*CA+CA/2-2*BL,CA/2,BL+BJ*CA/2+CA+SD,BL+BJ*CA/2+SD,0,-JPT/2
BLOCK,-CA/2,-BJ*CA+CA/2-2*BL,-BL-BJ*CA/2,BL+BJ*CA/2+SD,0,-JPT/2
FLST,2,4,3  
FITEM,2,90  
FITEM,2,97  
FITEM,2,81  
FITEM,2,74  
A,P51X
VSBA,      12,      77  
VDELE,      13, , ,1
wpoff,0,0,-JPT/2
CYL4,CA/2-BL,BJ*CA/2,BR5, , , ,JPT/2
CYL4,-(CA/2-BL),BJ*CA/2,BR5, , , ,JPT/2
CYL4,CA/2-BL,-BJ*CA/2,BR5, , , ,JPT/2
CYL4,-(CA/2-BL),-BJ*CA/2,BR5, , , ,JPT/2
CYL4,CA/2-BL,BJ*CA/2+BL+SD+BL,BR5, , , ,JPT/2
CYL4,-BJ*CA+CA/2-BL,BJ*CA/2+BL+SD+BL,BR5, , , ,JPT/2
CYL4,CA/2-BL,BJ*CA/2+BL+SD+BL+CA-2*BL,BR5, , , ,JPT/2
CYL4,-BJ*CA+CA/2-BL,BJ*CA/2+BL+SD+BL+CA-2*BL,BR5, , , ,JPT/2
FLST,3,4,6,ORDE,2   
FITEM,3,17  
FITEM,3,-20
VSBV,      11,P51X
FLST,3,4,6,ORDE,4   
FITEM,3,12  
FITEM,3,-13
FITEM,3,15  
FITEM,3,-16
VSBV,      10,P51X
ALLSEL,ALL
CMSEL,U,COLUMN
CMSEL,U,beam
VGLUE,ALL   
CM,jointplate,VOLU


!螺栓建模

VSEL,none
CYL4,CA/2-BL,BJ*CA/2,BR1, , , ,CD+JPT/2+BH1+BH2  
wpoff,0,0,CD+JPT/2
CYL4,CA/2-BL,BJ*CA/2,BR2, , , ,BH1
CYL4,CA/2-BL,BJ*CA/2,BR3, , , ,BH1+BH2
VOVLAP,ALL
CM,Bbolt1,VOLU

VSEL,none
wpoff,0,0,-(JPT/2+CD)
CYL4,-(CA/2-BL),BJ*CA/2,BR1, , , ,CD+JPT/2+BH1+BH2  
wpoff,0,0,CD+JPT/2
CYL4,-(CA/2-BL),BJ*CA/2,BR2, , , ,BH1
CYL4,-(CA/2-BL),BJ*CA/2,BR3, , , ,BH1+BH2
VOVLAP,ALL
CM,Bbolt2,VOLU

VSEL,none
wpoff,0,0,-(JPT/2+CD)
CYL4,CA/2-BL,-BJ*CA/2,BR1, , , ,CD+JPT/2+BH1+BH2   
wpoff,0,0,CD+JPT/2
CYL4,CA/2-BL,-BJ*CA/2,BR2, , , ,BH1
CYL4,CA/2-BL,-BJ*CA/2,BR3, , , ,BH1+BH2
VOVLAP,ALL
CM,Bbolt3,VOLU

VSEL,none
wpoff,0,0,-(JPT/2+CD)
CYL4,-(CA/2-BL),-BJ*CA/2,BR1, , , ,CD+JPT/2+BH1+BH2  
wpoff,0,0,CD+JPT/2
CYL4,-(CA/2-BL),-BJ*CA/2,BR2, , , ,BH1
CYL4,-(CA/2-BL),-BJ*CA/2,BR3, , , ,BH1+BH2
VOVLAP,ALL
CM,Bbolt4,VOLU

VSEL,none
wpoff,0,0,-(JPT/2+CD)
CYL4,CA/2-BL,BJ*CA/2+BL+SD+BL,BR1, , , ,CD+JPT/2+BH1+BH2  
wpoff,0,0,CD+JPT/2
CYL4,CA/2-BL,BJ*CA/2+BL+SD+BL,BR2, , , ,BH1
CYL4,CA/2-BL,BJ*CA/2+BL+SD+BL,BR3, , , ,BH1+BH2
VOVLAP,ALL
CM,Cbolt1,VOLU

VSEL,none
wpoff,0,0,-(JPT/2+CD)
CYL4,-BJ*CA+CA/2-BL,BJ*CA/2+BL+SD+BL,BR1, , , ,CD+JPT/2+BH1+BH2
wpoff,0,0,CD+JPT/2
CYL4,-BJ*CA+CA/2-BL,BJ*CA/2+BL+SD+BL,BR2, , , ,BH1
CYL4,-BJ*CA+CA/2-BL,BJ*CA/2+BL+SD+BL,BR3, , , ,BH1+BH2
VOVLAP,ALL
CM,Cbolt2,VOLU

VSEL,none
wpoff,0,0,-(JPT/2+CD)
CYL4,CA/2-BL,BJ*CA/2+BL+SD+BL+CA-2*BL,BR1, , , ,CD+JPT/2+BH1+BH2   
wpoff,0,0,CD+JPT/2
CYL4,CA/2-BL,BJ*CA/2+BL+SD+BL+CA-2*BL,BR2, , , ,BH1
CYL4,CA/2-BL,BJ*CA/2+BL+SD+BL+CA-2*BL,BR3, , , ,BH1+BH2
VOVLAP,ALL
CM,Cbolt3,VOLU

VSEL,none
wpoff,0,0,-(JPT/2+CD)
CYL4,-BJ*CA+CA/2-BL,BJ*CA/2+BL+SD+BL+CA-2*BL,BR1, , , ,CD+JPT/2+BH1+BH2  
wpoff,0,0,CD+JPT/2
CYL4,-BJ*CA+CA/2-BL,BJ*CA/2+BL+SD+BL+CA-2*BL,BR2, , , ,BH1
CYL4,-BJ*CA+CA/2-BL,BJ*CA/2+BL+SD+BL+CA-2*BL,BR3, , , ,BH1+BH2
VOVLAP,ALL
CM,Cbolt4,VOLU
!!回原点
WPCSYS,-1,0
!创建选择组件

CMGRP,BOLT,BBOLT1,BBOLT2,BBOLT3,BBOLT4,CBOLT1,CBOLT2,CBOLT3,CBOLT4  

VSEL,none
CMSEL,S,BEAM
CMSEL,A,COLUMN
VSEL,R,LOC,Y,-JPT/2-CD,-(JPT/2-CD+CR)
CM,WANJIAO1,VOLU

VSEL,none
CMSEL,S,BEAM
CMSEL,A,COLUMN
VSEL,R,LOC,Y,-(JPT/2-CD+CB),-(JPT/2-CD-CR+CB)
CM,WANJIAO2,VOLU

CMGRP,WANJIAO,WANJIAO1,WANJIAO2,

ALLSEL,ALL
CMSEL,S,BEAM
ASLV,R  
ASEL,R,LOC,X,-CA/2,CA/2+2*BL+BJ*CA+SD
ASEL,R,LOC,Y,-JPT/2,-JPT/2
ASEL,R,LOC,Z,JJ-BL,JJ+BL+BJ*CA
CM,BJ,AREA

ALLSEL,ALL
CMSEL,S,COLUMN
ASLV,R  
ASEL,R,LOC,X,-CA/2,CA/2+2*BL+BJ*CA+SD
ASEL,R,LOC,Y,-JPT/2,-JPT/2
ASEL,R,LOC,Z,JJ-BL,JJ+BL+BJ*CA
CM,CJ,AREA

ALLSEL,ALL
CMSEL,S,jointplate
ASLV,R  
ASEL,R,LOC,Y,-JPT/2
ASEL,R, , ,      88  
CM,JC,AREA

ALLSEL,ALL
CMSEL,S,jointplate
ASLV,R  
ASEL,R,LOC,Y,-JPT/2
ASEL,R, , ,      63
CM,JB,AREA

CMGRP,CCJJ,BJ,CJ
CMGRP,JJCC,JB,JC

CMSEL,S,BBOLT1   
ASLV   
FLST,5,4,5,ORDE,4   
FITEM,5,133
FITEM,5,-134
FITEM,5,144
FITEM,5,153
ASEL,R, , ,P51X
CM,B1,AREA   

CMSEL,S,BBOLT2  
ASLV
FLST,5,4,5,ORDE,4   
FITEM,5,164
FITEM,5,-165
FITEM,5,185
FITEM,5,332
ASEL,R, , ,P51X
CM,B2,AREA  

CMSEL,S,BBOLT3  
ASLV
FLST,5,4,5,ORDE,4   
FITEM,5,343
FITEM,5,-344
FITEM,5,348
FITEM,5,353
ASEL,R, , ,P51X
CM,B3,AREA  

CMSEL,S,BBOLT4  
ASLV
FLST,5,4,5,ORDE,4   
FITEM,5,364
FITEM,5,-365
FITEM,5,369
FITEM,5,374
ASEL,R, , ,P51X
CM,B4,AREA

CMSEL,S,CBOLT1   
ASLV
FLST,5,4,5,ORDE,4   
FITEM,5,385
FITEM,5,-386
FITEM,5,390
FITEM,5,395
ASEL,R, , ,P51X
CM,B5,AREA  

CMSEL,S,CBOLT2  
ASLV
FLST,5,4,5,ORDE,4   
FITEM,5,406
FITEM,5,-407
FITEM,5,411
FITEM,5,416
ASEL,R, , ,P51X
CM,B6,AREA  

CMSEL,S,CBOLT3     
ASLV
FLST,5,4,5,ORDE,4   
FITEM,5,427
FITEM,5,-428
FITEM,5,432
FITEM,5,437
ASEL,R, , ,P51X
CM,B7,AREA  

CMSEL,S,CBOLT4   
ASLV
FLST,5,4,5,ORDE,4   
FITEM,5,448
FITEM,5,-449
FITEM,5,453
FITEM,5,458
ASEL,R, , ,P51X
CM,B8,AREA

CMSEL,S,BEAM
CMSEL,A,JOINTPLATE  
ASLV   
FLST,5,5,5,ORDE,4   
FITEM,5,71  
FITEM,5,75  
FITEM,5,285
FITEM,5,-287
ASEL,R, , ,P51X
CM,C1,AREA  

CMSEL,S,BEAM
CMSEL,A,JOINTPLATE  
ASLV  
FLST,5,5,5,ORDE,5   
FITEM,5,54  
FITEM,5,-55
FITEM,5,79  
FITEM,5,83  
FITEM,5,169
ASEL,R, , ,P51X
CM,C2,AREA  

CMSEL,S,BEAM
CMSEL,A,JOINTPLATE  
ASLV  
FLST,5,5,5,ORDE,5   
FITEM,5,89  
FITEM,5,-90
FITEM,5,181
FITEM,5,-182
FITEM,5,273
ASEL,R, , ,P51X
CM,C3,AREA  

CMSEL,S,BEAM
CMSEL,A,JOINTPLATE  
ASLV  
FLST,5,5,5,ORDE,5   
FITEM,5,94  
FITEM,5,-95
FITEM,5,173
FITEM,5,-174
FITEM,5,271
ASEL,R, , ,P51X
CM,C4,AREA  

CMSEL,S,COLUMN  
CMSEL,A,JOINTPLATE  
ASLV  
FLST,5,5,5,ORDE,5   
FITEM,5,98  
FITEM,5,102
FITEM,5,147
FITEM,5,-148
FITEM,5,149
ASEL,R, , ,P51X
CM,C5,AREA  

CMSEL,S,COLUMN  
CMSEL,A,JOINTPLATE  
ASLV
FLST,5,5,5,ORDE,5   
FITEM,5,106
FITEM,5,-107
FITEM,5,119
FITEM,5,-120
FITEM,5,135
ASEL,R, , ,P51X
CM,C6,AREA

CMSEL,S,COLUMN  
CMSEL,A,JOINTPLATE  
ASLV   
FLST,5,5,5,ORDE,5   
FITEM,5,111
FITEM,5,-112
FITEM,5,123
FITEM,5,-124
FITEM,5,139
ASEL,R, , ,P51X
CM,C7,AREA  

CMSEL,S,COLUMN  
CMSEL,A,JOINTPLATE  
ASLV  
FLST,5,5,5,ORDE,5   
FITEM,5,127
FITEM,5,-128
FITEM,5,131
FITEM,5,-132
FITEM,5,137
ASEL,R, , ,P51X
CM,C8,AREA  


ALLSEL,ALL  
SAVE

!!!!!!!划分网格!!!!!!!
LSEL,S,LOC,Y,0
LSEL,R,radius,,BR5
LESIZE,ALL, , ,4, , , , ,1

LSEL,S,LOC,Y,-JPT/2
LSEL,R,radius,,BR4
LESIZE,ALL, , ,4, , , , ,1
LSEL,S,LOC,Y,-JPT/2
LSEL,R,radius,,BR5
LESIZE,ALL, , ,4, , , , ,1

LSEL,S,LOC,Y,-JPT/2-CD
LSEL,R,radius,,BR4
LESIZE,ALL, , ,4, , , , ,1
LSEL,S,LOC,Y,-JPT/2-CD
LSEL,R,radius,,BR5
LESIZE,ALL, , ,4*2, , , , ,1
LSEL,S,LOC,Y,-JPT/2-CD
LSEL,R,radius,,BR2
LESIZE,ALL, , ,4, , , , ,1
LSEL,S,LOC,Y,-JPT/2-CD
LSEL,R,radius,,BR1
LESIZE,ALL, , ,4*2, , , , ,1
LSEL,S,LOC,Y,-JPT/2-CD
LSEL,R,radius,,BR3
LESIZE,ALL, , ,4, , , , ,1

LSEL,S,LOC,Y,-JPT/2-CD-BH1
LSEL,R,radius,,BR2
LESIZE,ALL, , ,4, , , , ,1
LSEL,S,LOC,Y,-JPT/2-CD-BH1
LSEL,R,radius,,BR3
LESIZE,ALL, , ,4*2, , , , ,1
LSEL,S,LOC,Y,-JPT/2-CD-BH1
LSEL,R,radius,,BR1
LESIZE,ALL, , ,4*2, , , , ,1

LSEL,S,LOC,Y,-JPT/2-CD-BH1-BH2
LSEL,R,radius,,BR3
LESIZE,ALL, , ,4, , , , ,1
LSEL,S,LOC,Y,-JPT/2-CD-BH1-BH2
LSEL,R,radius,,BR1
LESIZE,ALL, , ,4, , , , ,1

LSEL,S,radius,,CR
LSEL,A,radius,,CR-CD
LESIZE,ALL, , ,3, , , , ,1

LSEL,S,length,,CB-CR*2
LSEL,A,length,,CC-CR
LESIZE,ALL,15, , , , , , ,1

LSEL,S,length,,CA-CR*2
LESIZE,ALL,12, , , , , , ,1

LSEL,S,length,,BJ*CA+2*BL
LSEL,A,length,,BJ*CA+2*BL+SD
LSEL,A,length,,CA
LSEL,A,length,,BJ*CA+2*BL-CA
LESIZE,ALL,8, , , , , , ,1

LSEL,S,length,,JPT/2
LSEL,A,length,,CD
LSEL,A,length,,BH1
LSEL,A,length,,BH2
LESIZE,ALL,2, , , , , , ,1  

LSEL,S,length,,JJ-BL
LESIZE,ALL,30, , , , , , ,1

LSEL,S,LOC,Y,0,-JPT/2
LSEL,R,LOC,X,CA,CA+BJ*CA/2
LSEL,R,LOC,Z,JJ,JJ+BJ*CA-CA
LESIZE,ALL,15, , , , , , ,1

ALLSEL,ALL
MOPT,QMESH,MAIN

CMSEL,S,JOINTPLATE  
VATT,       1, ,   1,       0   
VSWEEP,ALL

VSEL,NONE
CMSEL,S,BEAM
CMSEL,A,COLUMN  
CMSEL,U,WANJIAO
VATT,       2, ,   1,       0   
VSWEEP,ALL

VSEL,NONE
CMSEL,S,WANJIAO
VATT,       3, ,   1,       0   
VSWEEP,ALL

VSEL,NONE
CMSEL,S,BOLT
VATT,       4, ,   1,       0   
VSWEEP,ALL


!!!!!!接触建立!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!节点板与C型钢!!!!!!!
ET,3,170
ET,4,174
KEYOPT,4,7,1   !时间步长二分
KEYOPT,4,9,2
KEYOPT,4,12,2
allsel,all
mat,5
R,5,,,-2000
REAL,5
!! 建立目标面
ASEL,S,,,JJCC
TYPE,3
NSLA,S,1
ESLN,S,0
ESURF,ALL   
!!建立接触面
ASEL,S,,,CCJJ
TYPE,4
NSLA,S,1
ESLN,S,0
ESURF,ALL   
ALLSEL
!!!!!!!螺栓杆与孔壁的接触!!!!!!!!!
*DO, n, 1, 8
ET,5,170
ET,6,174
KEYOPT,6,7,1  
KEYOPT,6,9,2
KEYOPT,6,12,1
allsel,all
mat,6
R,%n+5%,,,-2000
REAL,%n+5%
!! 建立目标面－螺栓杆
ASEL,S,,,B%n%
TYPE,3  
NSLA,S,1
ESLN,S,0
ESURF,all
!!建立接触面－孔壁
ASEL,S,,,C%n%
TYPE,4
NSLA,S,1
ESLN,S,0
ESURF,all
*ENDDO

!!!!建立预紧单元!!!!!
/PREP7   
ALLSEL,ALL
psmesh,BB1,pretens,,volu,18,0,y,0,,,,Bbolt-1  
psmesh,BB2,pretens,,volu,68,0,y,0,,,,Bbolt-2  
psmesh,BB3,pretens,,volu,74,0,y,0,,,,Bbolt-3  
psmesh,BB4,pretens,,volu,80,0,y,0,,,,Bbolt-4  
psmesh,CB1,pretens,,volu,86,0,y,0,,,,Cbolt-1  
psmesh,CB2,pretens,,volu,92,0,y,0,,,,Cbolt-2  
psmesh,CB3,pretens,,volu,98,0,y,0,,,,Cbolt-3  
psmesh,CB4,pretens,,volu,104,0,y,0,,,,Cbolt-4


!!!!边界约束条件!!!!!!!!!!!!!!!
NSEL,ALL
NSEL,s,LOC,Z,0   
D,ALL, , , , , ,UX,UY,UZ,

NSEL,ALL
NSEL,S,LOC,Z,2*JJ+BJ*CA
D,ALL, , , , , ,UX,UY,UZ,   

NSEL,S,LOC,Y,0
DSYM,SYMM,Y, ,

SAVE
ALLSEL,ALL

!!!施加预应力!!!
/SOLU   
ANTYPE,0
AUTOTS,ON
NEQIT,100
NLGEOM,ON
NSUBST,1,,,1
SLOAD,ALL,9,LOCK,FORC,PREFOR, 1,2
ALLSEL,ALL  
*do,n,1,2
solcontrol,on !优化非线性求解
time,n
EQSLV,PCG,1E-2,3
CNVTOL,F, ,0.05,2,0.04
LSWRITE,n
*enddo
LSSOLVE,1,2,1,
save

!!!单调荷载!!!
/solu
ANTYPE,,REST
ALLSEL,ALL
num=node(CA/2+SD+JJ+BJ*CA+BL,-JPT/2,JJ+BJ*CA-CA/2)
D,num,UZ,DISP
NSUBST,500,1000,20
allsel,all  
time,3  
NLGEOM,1
OUTRES,alL,1
allsel
KBC,0 !采用渐进式加载方式加载
lnsrch,off !关闭线性搜索
EQSLV,PCG,1E-2,3
CNVTOL,F, ,0.05,2,0.04
LSWRITE,3,  
FINISH  
/SOLU   
LSSOLVE,3,3,1,


上边是用solid45建的模型，由于对称性，模型只建了一半，用预紧单元给螺栓施加预应力，钢材间，螺杆螺孔间定义了接触对，柱子上下铰接，梁端施加位移。
还有一个是用shell181建的模型，螺栓使用beam4，螺栓和型钢使用combin39连接，柱子上下铰接，梁端施加位移。
两个模型都运行不出来，现求助高手帮忙。
不胜感激

*SET,CA,160!C型钢参数   
*SET,CB,65  
*SET,CC,20  
*SET,CD,1.8
*SET,CR,6   
*SET,EP,35!螺栓端距
*SET,BJ,1.25 !螺栓间距参数1.0;1.25;1.5
*SET,JPT,6  !节点板厚度  
*SET,SD,10  !安装间隙   
*SET,JJ,600 !支撑间距   

*SET,DISP,50     !施加位移

/PREP7  

ET,1,SHELL181   
R,1,CD, , , , , ,
R,2,JPT, , , , , ,

ET,2,COMBIN39,,,1   
R,3,0,0,0.1,5000,0.5,8000,  
RMORE,1,9000,1.5,12000,3,15000,

ET,3,COMBIN39,,,2
R,4,0,0,0.1,5000,0.5,8000,   
RMORE,1,9000,1.5,12000,3,15000,

ET,4,COMBIN39,,,3   
R,5,0,0,0.1,50000000,   

ET,5,BEAM4  
R,6,200.905,3210,3210,1000,1000, ,
RMORE, ,6421,10000,10000, , ,   

!!!!!!节点板!!!!!
mp,ex,1,176000  
mp,prxy,1,0.3   
tb,kinh,1,1,4   
tbpt,,235/(176000),235  
tbpt,,0.002,235
tbpt,,0.02,235  
tbpt,,0.1,235/2

!!!!!!C型钢!!!!!
mp,ex,2,176000  
mp,prxy,2,0.3   
tb,kinh,2,1,4   
tbpt,,235/(176000),235  
tbpt,,0.002,235
tbpt,,0.02,235  
tbpt,,0.1,235/2

!!!!!!螺栓!!!!!!
mp,ex,3,176000  
mp,prxy,3,0.3   
tb,kinh,3,1,4   
tbpt,,235*1.2/(176000),235*1.2  
tbpt,,0.1,235*1.2   
tbpt,,0.02,235*1.2  
tbpt,,0.1,235*1.2/2

!!!!!!截面属性!!!!!!
SECTYPE,1,SHELL,RECT
SECDATA,CD,2, ,5

SECTYPE,2,SHELL,RECT
SECDATA,JPT,1, ,5

SECTYPE,3,BEAM,CSOLID
SECDATA,BR,8,4

!!!!!开始建模!!!!!!!!

K,1,-(CA-CD)/2+CC,-JPT/2-CB+CD/2,0,
K,2,-(CA-CD)/2,-JPT/2-CB+CD/2,0,
K,3,-(CA-CD)/2,-(JPT+CD)/2,0,   
K,4,-(CA-CD)/2+EP,-(JPT+CD)/2,0,
K,5,(CA-CD)/2-EP,-(JPT+CD)/2,0,
k,6,(CA-CD)/2,-(JPT+CD)/2,0,
K,7,(CA-CD)/2,-JPT/2-CB+CD/2,0,
K,8,(CA-CD)/2-CC,-JPT/2-CB+CD/2,0,  
K,9,(CA-CD)/2-CC,-JPT/2-CB+CD/2,BJ*CA/2,
K,10,(CA-CD)/2-CC,-JPT/2-CB+CD/2,BJ*CA/2+EP,
K,11,(CA-CD)/2-CC,-JPT/2-CB+CD/2,BJ*CA/2+EP+JJ/12,  
LSTR,       1,       2  
LSTR,       2,       3  
LSTR,       3,       4  
LSTR,       4,       5  
LSTR,       5,       6  
LSTR,       6,       7  
LSTR,       7,       8  
LSTR,       8,       9  
LSTR,       9,       10
LSTR,       10,      11
LFILLT,1,2,CR-CD/2, ,   
LFILLT,2,3,CR-CD/2, ,   
LFILLT,5,6,CR-CD/2, ,   
LFILLT,6,7,CR-CD/2, ,   
LSEL,S,LOC,Z,0  
ADRAG,       ALL, , , , , ,      8  
LSEL,S,LOC,Z,BJ*CA/2
ADRAG,       ALL, , , , , ,      9  
LSEL,S,LOC,Z,BJ*CA/2+EP
ADRAG,       ALL, , , , , ,      10
K,1001,CA/2+SD,-JPT/2-CB+CD/2,-(CA-CD)/2+CC+(BJ*CA-CA+CD)/2+EP,
K,1002,CA/2+SD,-JPT/2-CB+CD/2,-(CA-CD)/2+(BJ*CA-CA+CD)/2+EP,
K,1003,CA/2+SD,-(JPT+CD)/2,-(CA-CD)/2+(BJ*CA-CA+CD)/2+EP,   
K,1004,CA/2+SD,-(JPT+CD)/2,-(CA-CD)/2+EP+(BJ*CA-CA+CD)/2+EP,
K,1005,CA/2+SD,-(JPT+CD)/2,(CA-CD)/2-EP+(BJ*CA-CA+CD)/2+EP,
k,1006,CA/2+SD,-(JPT+CD)/2,(CA-CD)/2+(BJ*CA-CA+CD)/2+EP,
K,1007,CA/2+SD,-JPT/2-CB+CD/2,(CA-CD)/2+(BJ*CA-CA+CD)/2+EP,
K,1008,CA/2+SD,-JPT/2-CB+CD/2,(CA-CD)/2-CC+(BJ*CA-CA+CD)/2+EP,  
K,1009,CA/2+SD+EP,-JPT/2-CB+CD/2,(CA-CD)/2-CC+(BJ*CA-CA+CD)/2+EP,   
K,1010,CA/2+SD+BJ*CA+EP,-JPT/2-CB+CD/2,(CA-CD)/2-CC+(BJ*CA-CA+CD)/2+EP,
K,1011,CA/2+SD+BJ*CA+EP+JJ/12,-JPT/2-CB+CD/2,(CA-CD)/2-CC+(BJ*CA-CA+CD)/2+EP,   
LSTR,       1001,       1002
LSTR,       1002,       1003
LSTR,       1003,       1004
LSTR,       1004,       1005
LSTR,       1005,       1006
LSTR,       1006,       1007
LSTR,       1007,       1008
LSTR,       1008,       1009
LSTR,       1009,       1010
LSTR,       1010,       1011
LFILLT,84,85,CR-CD/2, ,
LFILLT,85,86,CR-CD/2, ,
LFILLT,88,89,CR-CD/2, ,
LFILLT,89,90,CR-CD/2, ,
LSEL,S,LOC,X,CA/2+SD
ADRAG,       ALL, , , , , ,      91
LSEL,S,LOC,X,CA/2+SD+EP
ADRAG,       ALL, , , , , ,      92
LSEL,S,LOC,X,CA/2+SD+BJ*CA+EP   
ADRAG,       ALL, , , , , ,      93
ASEL,S,LOC,Z,BJ*CA/2+EP,BJ*CA/2+EP+JJ/12
ASEL,R,LOC,X,-CA/2,CA/2
AGEN,12,ALL, , , , ,JJ/12 , ,0  
ASEL,S,LOC,X,CA/2+SD+BJ*CA+EP,CA/2+SD+BJ*CA+EP+JJ/12
AGEN,12,ALL, , ,JJ/12, , , ,0   
ASEL,S,LOC,X,-CA/2,CA/2
ARSYM,Z,ALL,
ASEL,ALL
ARSYM,Y,ALL
/VIEW,1,1,1,1   
/ANG,1  
/REP,FAST   



K,2000,CA/2+SD+2*EP+BJ*CA,0,BJ*CA/2-CA+CD+EP,
wpro,,90.000000,

RECTNG,-(CA-CD)/2,(CA-CD)/2,BJ*CA/2+EP,-BJ*CA/2-EP,
RECTNG,(CA-CD)/2,CA/2+SD+2*EP+BJ*CA,BJ*CA/2+EP,BJ*CA/2-CA+CD+EP+EP,
A,1888,1891,1892,2000
RECTNG,-(CA-CD)/2+EP,(CA-CD)/2-EP,BJ*CA/2,-BJ*CA/2,
RECTNG,CA/2+SD+EP,CA/2+SD+EP+BJ*CA,BJ*CA/2,BJ*CA/2-CA+CD+EP+EP,

ASEL,S,LOC,Y,0  
AOVLAPE,ALL

alls
lsla
lsel,invert
ldele,all   
ALLS
KSLL
KSEL,INVERT
KDELE,ALL   


ALLSEL,ALL
NUMMRG,KP,0.1,0.1, ,LOW



!!!!!!!创建实体组件!!!!!!!!!!!!
ASEL,S,LOC,Y,0
CM,JP,AREA

ASEL,S,LOC,Y,-JPT/2,-JPT/2-CB
ASEL,R,LOC,X,-CA/2,CA/2
CM,COLUMN1,AREA

ASEL,S,LOC,Y,JPT/2,JPT/2+CB
ASEL,R,LOC,X,-CA/2,CA/2
CM,COLUMN2,AREA

ASEL,S,LOC,Y,-JPT/2,-JPT/2-CB
ASEL,R,LOC,X,CA/2+SD,CA/2+SD+BJ*CA+EP+JJ
CM,BEAM1,AREA

ASEL,S,LOC,Y,JPT/2,JPT/2+CB
ASEL,R,LOC,X,CA/2+SD,CA/2+SD+BJ*CA+EP+JJ
CM,BEAM2,AREA

ASEL,S,LOC,Y,JPT/2+CD,JPT/2+CD+CR
ASEL,A,LOC,Y,-(JPT/2+CD),-(JPT/2+CD+CR)
ASEL,A,LOC,Y,JPT/2-CD+CB,JPT/2-CD-CR+CB
ASEL,A,LOC,Y,-(JPT/2-CD+CB),-(JPT/2-CD-CR+CB)
CM,CORNER,AREA

!!!!!!!划分网格!!!!!!!  
LSEL,S,radius,,CR-CD/2  
LESIZE,ALL, , ,3, , , , ,1  
LSEL,S,length,,EP   
LSEL,A,length,,EP-CR+CD/2   
LESIZE,ALL, , ,2, , , , ,1  
LSEL,S,length,,CB-CD-2*(CR-CD/2)
LSEL,A,length,,CC-CR+CD/2   
LSEL,A,length,,CA-2*EP-CD   
LESIZE,ALL,15, , , , , , ,1
LSEL,S,length,,BJ*CA
LSEL,A,length,,BJ*CA/2  
LESIZE,ALL,15, , , , , , ,1
LSEL,S,length,,JJ/12
LESIZE,ALL, , ,1, , , , ,1  

LSEL,S,length,,BJ*CA+2*EP
LSEL,A,length,,CA-CD
LSEL,A,length,,BJ*CA
LSEL,A,length,,CD/2+SD+EP
LSEL,A,length,,CA-2*EP-CD
LSEL,A,length,,CD/2+SD+EP+BJ*CA+EP
LSEL,A,length,,EP
LSEL,A,length,,CA-EP-CD
LSEL,A,length,,BJ*CA+2*EP-(CA-EP-CD)
LESIZE,ALL,15, , , , , , ,1

LSEL,S,LOC,Y,0
LSEL,R,LOC,X,CA,CA+BJ*CA
LSEL,R,LOC,Z,-BJ*CA+EP/2,-BJ*CA+CA
LESIZE,ALL, , ,20 , , , , ,1

ASEL,NONE
CMSEL,S,JP  
AATT,1,2,1,,2
AMESH,ALL

ASEL,NONE
CMSEL,S,COLUMN1
CMSEL,A,COLUMN2
CMSEL,A,BEAM1
CMSEL,A,BEAM2
CMSEL,U,CORNER
AATT,2,1,1,,1
AMESH,ALL

ASEL,NONE
CMSEL,S,CORNER
AATT,3,1,1,,1
AMESH,ALL

ALLSEL,ALL
NUMMRG,NODE,0.1,0.1, ,LOW
NUMCMP,ALL


!!!!!!!!建立弹簧单元!!!!!!!!!
BN11=node(-(CA-CD)/2+EP,(JPT+CD)/2,BJ*CA/2)
BN21=node(-(CA-CD)/2+EP,(JPT+CD)/2,-BJ*CA/2)
BN31=node((CA-CD)/2-EP,(JPT+CD)/2,BJ*CA/2)
BN41=node((CA-CD)/2-EP,(JPT+CD)/2,-BJ*CA/2)
BN51=node(CA/2+SD+EP,(JPT+CD)/2,BJ*CA/2)
BN61=node(CA/2+SD+EP,(JPT+CD)/2,BJ*CA/2-CA+CD+EP+EP)
BN71=node(CA/2+SD+EP+BJ*CA,(JPT+CD)/2,BJ*CA/2)
BN81=node(CA/2+SD+EP+BJ*CA,(JPT+CD)/2,BJ*CA/2-CA+CD+EP+EP)

BN12=node(-(CA-CD)/2+EP,0,BJ*CA/2)
BN22=node(-(CA-CD)/2+EP,0,-BJ*CA/2)
BN32=node((CA-CD)/2-EP,0,BJ*CA/2)
BN42=node((CA-CD)/2-EP,0,-BJ*CA/2)
BN52=node(CA/2+SD+EP,0,BJ*CA/2)
BN62=node(CA/2+SD+EP,0,BJ*CA/2-CA+CD+EP+EP)
BN72=node(CA/2+SD+EP+BJ*CA,0,BJ*CA/2)
BN82=node(CA/2+SD+EP+BJ*CA,0,BJ*CA/2-CA+CD+EP+EP)

BN13=node(-(CA-CD)/2+EP,-(JPT+CD)/2,BJ*CA/2)
BN23=node(-(CA-CD)/2+EP,-(JPT+CD)/2,-BJ*CA/2)
BN33=node((CA-CD)/2-EP,-(JPT+CD)/2,BJ*CA/2)
BN43=node((CA-CD)/2-EP,-(JPT+CD)/2,-BJ*CA/2)
BN53=node(CA/2+SD+EP,-(JPT+CD)/2,BJ*CA/2)
BN63=node(CA/2+SD+EP,-(JPT+CD)/2,BJ*CA/2-CA+CD+EP+EP)
BN73=node(CA/2+SD+EP+BJ*CA,-(JPT+CD)/2,BJ*CA/2)
BN83=node(CA/2+SD+EP+BJ*CA,-(JPT+CD)/2,BJ*CA/2-CA+CD+EP+EP)


N,10011,-(CA-CD)/2+EP,(JPT+CD)/2,BJ*CA/2,,,,
N,10021,-(CA-CD)/2+EP,(JPT+CD)/2,-BJ*CA/2,,,,
N,10031,(CA-CD)/2-EP,(JPT+CD)/2,BJ*CA/2,,,,
N,10041,(CA-CD)/2-EP,(JPT+CD)/2,-BJ*CA/2,,,,
N,10051,CA/2+SD+EP,(JPT+CD)/2,BJ*CA/2,,,,
N,10061,CA/2+SD+EP,(JPT+CD)/2,BJ*CA/2-CA+CD+EP+EP,,,,
N,10071,CA/2+SD+EP+BJ*CA,(JPT+CD)/2,BJ*CA/2,,,,
N,10081,CA/2+SD+EP+BJ*CA,(JPT+CD)/2,BJ*CA/2-CA+CD+EP+EP,,,,

N,10012,-(CA-CD)/2+EP,0,BJ*CA/2,,,,
N,10022,-(CA-CD)/2+EP,0,-BJ*CA/2,,,,
N,10032,(CA-CD)/2-EP,0,BJ*CA/2,,,,
N,10042,(CA-CD)/2-EP,0,-BJ*CA/2,,,,
N,10052,CA/2+SD+EP,0,BJ*CA/2,,,,
N,10062,CA/2+SD+EP,0,BJ*CA/2-CA+CD+EP+EP,,,,
N,10072,CA/2+SD+EP+BJ*CA,0,BJ*CA/2,,,,
N,10082,CA/2+SD+EP+BJ*CA,0,BJ*CA/2-CA+CD+EP+EP,,,,

N,10013,-(CA-CD)/2+EP,-(JPT+CD)/2,BJ*CA/2,,,,
N,10023,-(CA-CD)/2+EP,-(JPT+CD)/2,-BJ*CA/2,,,,
N,10033,(CA-CD)/2-EP,-(JPT+CD)/2,BJ*CA/2,,,,
N,10043,(CA-CD)/2-EP,-(JPT+CD)/2,-BJ*CA/2,,,,
N,10053,CA/2+SD+EP,-(JPT+CD)/2,BJ*CA/2,,,,
N,10063,CA/2+SD+EP,-(JPT+CD)/2,BJ*CA/2-CA+CD+EP+EP,,,,
N,10073,CA/2+SD+EP+BJ*CA,-(JPT+CD)/2,BJ*CA/2,,,,
N,10083,CA/2+SD+EP+BJ*CA,-(JPT+CD)/2,BJ*CA/2-CA+CD+EP+EP,,,,



*DO, m, 1, 3
*DO, n, 1, 8
TYPE,2   
MAT,3
REAL,3     
SECNUM,3
TSHAP,LINE
E,BN%n%%m%,100%n%%m%
*ENDDO
*ENDDO


*DO, m, 1, 3
*DO, n, 1, 8
TYPE,3   
MAT,3
REAL,4     
SECNUM,3
TSHAP,LINE
E,BN%n%%m%,100%n%%m%
*ENDDO
*ENDDO


*DO, m, 1, 3
*DO, n, 1, 8
TYPE,4   
MAT,3
REAL,5     
SECNUM,3
TSHAP,LINE
E,BN%n%%m%,100%n%%m%
*ENDDO
*ENDDO



*DO, n, 1, 8
*DO, m, 1, 2
TYPE,5   
MAT,3
REAL,6   
SECNUM,3
TSHAP,LINE
E,100%n%%m%,100%n%%m+1%
*ENDDO
*ENDDO


!!!!边界约束条件!!!!!!!!!!!!!!!
NSEL,ALL
NSEL,s,LOC,Z,-(JJ+BJ*CA/2+EP)   
D,ALL,ALL,0,

NSEL,ALL
NSEL,S,LOC,Z,JJ+BJ*CA/2+EP
D,ALL,ALL,0,

NSEL,ALL
NSEL,S,LOC,X,CA/2+SD+EP+BJ*CA+JJ
D,ALL,UY,0,



NSEL,NONE
*DO, n, 1, 8
*DO, m, 1, 3
NSEL,A, , ,   100%n%%m%
*ENDDO
*ENDDO
D,ALL, ,0, , , ,ROTX,ROTY,ROTZ, , ,

NSEL,NONE
*DO, n, 1, 8
*DO, m, 1, 3
NSEL,A, , ,   BN%n%%m%
*ENDDO
*ENDDO
D,ALL, ,0, , , ,ROTX,ROTY,ROTZ, , ,

SAVE
FINISH
ALLSEL,ALL

/SOLU   
ANTYPE,0
AUTOTS,ON
NEQIT,200
NLGEOM,ON
NSUBST,500,1000,20
time,10
OUTRES,ALL,ALL
num1=node(CA/2+SD+EP+BJ*CA+JJ,-(JPT+CD)/2,0)
num2=node(CA/2+SD+EP+BJ*CA+JJ,(JPT+CD)/2,0)
D,num1,UZ,DISP
D,num2,UZ,DISP
KBC,0 !采用渐进式加载方式加载
lnsrch,ON !关闭线性搜索
EQSLV,PCG,1E-2,1
CNVTOL,F, ,0.05,2,0.04  

ALLSEL,ALL
SOLVE

其实实际的试验模型很简单，就是两个C型钢背靠背，中间夹一块钢板，用高强螺栓连起来，柱子两端固定后，两端加反复荷载。加载时，在平面内，型钢，钢板，螺栓都承受剪力。型钢和螺栓，型钢和钢板之间有摩擦。不知道ansys能不能解之类问题。
国外有类似的研究，用的是shell（壳）单元，beam（梁）单元，spring（弹簧）单元建模，将接触、摩擦都简化为弹簧的性质。

这种连接的接触状态比较复杂，在荷载较小的时候，剪力由预紧螺栓的摩擦力承受，力增大后，螺栓和型钢间出现滑移，待螺杆和孔壁接触后，螺栓既受剪力，又受孔壁的压力。共三种状态。再加上材料非线性，几何非线性，问题就更加复杂。
所以我现在倾向于使用shell单元，不使用接触，将接触状态的非线性直接定义为combin（弹簧）单元的属性。


这个shell模型估计还是自由度的问题，combin和beam能不能连接，怎么连接，需不需要偶合？
非线性, 刚性, 节点, 求解, 型钢

相关帖子
回复

使用道具 举报

zjj820917
发表于 2009-3-28 22:24:45 | 显示全部楼层 来自 湖南长沙
Simdroid开发平台
不错!学习!
要是能用简单的命令流搞出来就好了!
里面又太多的界面操作
回复 赞 不支持

使用道具 举报

ytdzqh
发表于 2011-12-14 16:45:03 | 显示全部楼层 来自 内蒙古包头
高手啊，佩服
回复 赞 不支持

使用道具 举报

返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 注册

本版积分规则

Simapps系列直播

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

GMT+8, 2024-9-29 17:28 , Processed in 0.055103 second(s), 13 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

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