型钢混凝土梁模拟，加载50kn梁就破坏了

wyslft

模拟锈蚀型钢混凝土梁时，加弹簧单元梁的承载力弱，超过50kn梁就破坏的很严重～
跟试验描述的中加载至200kn，相差太大了
试验描述：梁长2m，截面250×180，混凝土c30，型钢界面：100×8×6，型钢屈服强度210，纵筋4根12，HRB400,箍筋[email=150@6]150@6[/email]，HPB235
命令流如下，请各位高手指点迷津

Fcu=33   
Fc=22.1
Ft=2.926  
! 定义混凝土单元
ET,1,65
KEYOPT,1,1,1      !不考虑大变形
KEYOPT,1,7,1      !考虑考虑应力松弛，有助于计算收敛
! 定义混凝土的实常数
R,1,,,
! 定义混凝土的材料属性
MP,PRXY,1,0.2        
MP,DENS,1,2500E-9   
TB,MISO,1,1,18      
Eo=32000      
Io=0.002     
Es=Fcu*0.85/0.002      
*do,I,0.0002,0.005,0.0003
  D_equation=1+(Eo/Es-2)*(I/Io)+(I/Io)*(I/Io)  
*if,I,EQ,0.0002,THEN
  Ec=Eo/D_equation     
  MP,EX,1,Ec   
*elseif,I,EQ,0.0005,THEN
  TBPT,,I,Ec*I
*else
  TBPT,,I,Eo*I/D_equation     
*endif
*enddo     
TB,CONCR,1
TBDATA,,0.35,0.75,Ft,-1  

!定义型钢下翼缘单元
ET,2,SOLID45
R,2
!MPTEMP,1,0
MP,EX,2,1.97E5  
MP,PRXY,2,0.25
TB,BKIN,2,1,2      
TBDATA,,267.9,0     
                                 

!定义型钢上翼缘和腹板单元
ET,3,SOLID45
R,3
MP,EX,3,2.06E5
MP,PRXY,3,0.25
!MPDATA,DENS,3,,7.850E-6
TB,BKIN,3,1,2,1               
TBDATA,,270,0         
      
!定义箍筋单元
!箍筋直径和面积
PI=ACOS(-1)
DSTU=6
ASTU=0.25*PI*DSTU*DSTU
ET,4,LINK8
! 定义箍筋的实常数
R,4,ASTU
MP,EX,4,2.06E5
MP,PRXY,4,0.25
!MPDATA,DENS,4,,7.850E-6
TB,BKIN,4,1,2,1
TBDATA,,265,0         

!纵向受力钢筋直径和面积
DPULL_LW=12
APULL_LW=0.25*PI*DPULL_LW*DPULL_LW
ET,5,LINK8
R,5,APULL_LW
MP,EX,5,2.06E5
MP,PRXY,5,0.25
!MPDATA,DENS,5,,7850E-12
TB,BKIN,5,1,2,1
TBDATA,,552,0   

!定义垫板的单元
ET,6,SOLID45
R,6
MP,EX,6,2.06E5
MP,PRXY,6,0.25

!定义弹簧单元
!Z方向
ET,7,COMBIN39
KEYOPT,7,1,0
KEYOPT,7,2,0
KEYOPT,7,3,3

!上翼缘外侧   
R,7,-0.006,-375.54,-0.005,-376.20,-0.004,-376.86
rmore,-0.003,-377.53,-0.002,-378.19,-0.001,-378.85
rmore,0,0,0.001,378.85,0.002,378.19
rmore,0.003,377.53,0.004,376.86,0.005,376.20

R,8,-0.006,-525.75,-0.005,-526.68,-0.004,-527.61,
RMORE,-0.003,-528.54,-0.002,-529.46,-0.001,-530.39
RMORE,0,0,0.001,530.39,0.002,529.46
RMORE,0.003,528.54,0.004,527.61,0.005,526.68

R,9,-0.006,-347.37,-0.005,-347.99,-0.004,-348.60
RMORE,-0.003,-349.21,-0.002,-348.82,-0.001,-350.43
RMORE,0,0,0.001,350.43,0.002,349.82,
RMORE,0.003,349.21,0.004,348.60,0.005,347.99

R,10,-0.006,-558.61,-0.005,-559.60,-0.004,-560.58
RMORE,-0.003,-561.57,-0.002,-562.55,-0.001,-563.53
RMORE,0,0,0.001,563.53,0.002,562.55,
RMORE,0.003,561.57,0.004,560.58,0.005,559.60

R,11,-0.006,-591.46,-0.005,-592.51,-0.004,-593.55
RMORE,-0.003,-594.60,-0.002,-595.64,-0.001,-596.68
RMORE,0,0,0.001,596.68,0.002,595.64,
RMORE,0.003,594.60,0.004,593.55,0.005,592.51

!X方向
ET,8,COMBIN39
KEYOPT,8,3,1
R,100,-1E-6,-1E6,0,0,1E-6,1E6

!Y方向
ET,9,COMBIN39
KEYOPT,9,3,2
R,101,-1E-6,-1E6,0,0,1E-6,1E6

!型钢截面
L=2000
HS=100
BS=100
DS=6
TS=8
HB=250
BB=180         
K,1,-BS/2,HS/2,0
K,2,-BS/2,HS/2-TS,0
K,3,-DS/2,HS/2-TS,0
K,4,-DS/2,-HS/2+TS,
K,5,-BS/2,-HS/2+TS,0
K,6,-BS/2,-HS/2+0.1,0
K,7,BS/2,-HS/2+0.1,0
K,8,BS/2,-HS/2+TS,0
K,9,DS/2,-HS/2+TS,
K,10,DS/2,HS/2-TS,0
K,11,BS/2,HS/2-TS,0
K,12,BS/2,HS/2,0
A,1,2,3,4,5,6,7,8,9,10,11,12
K,13,-BB/2,-HB/2+0.1,0   
K,14,-BB/2,HB/2,0  
K,15,BB/2,HB/2,0
K,16,BB/2,-HB/2+0.1,0  
A,13,14,15,16
ASBA,2,1,SEPO,DELETE,KEEP         
NUMCMP,AREA
NUMCMP,KP
NUMCMP,LINE
                             
VEXT,1,,,,,1000
VEXT,2,,,,,1000

!支座垫板
KP_START_NUM=KPINQR(,14)                                
K,,BB/2,-HB/2+0.1,825   
K,,-BB/2,-HB/2+0.1,825
K,,-BB/2,-HB/2+0.1,1000
K,,BB/2,-HB/2+0.1,1000
A,KP_START_NUM+1,KP_START_NUM+2,KP_START_NUM+3,KP_START_NUM+4
AREA_START_NUM=ARINQR(,14)
VEXT,AREA_START_NUM,,,,-40

!集中荷载处垫板
KP_START_NUM=KPINQR(,14)
K,,-BB/2,HB/2,375
K,,BB/2,HB/2,375
K,,BB/2,HB/2,525
K,,-BB/2,HB/2,525
A,KP_START_NUM+1,KP_START_NUM+2,KP_START_NUM+3,KP_START_NUM+4   
AREA_START_NUM=ARINQR(,14)
VEXT,AREA_START_NUM,,,,40,

!垫板和混凝土共面
VGLUE,2,3
VGLUE,6,4

!分割出腹板、翼缘位置
WPSTYL,DEFA
WPOFFS,-BS/2
WPROTA,,,90
VSBW,ALL
WPOFFS,,,BS/2-DS/2
VSBW,ALL
WPOFFS,,,DS
VSBW,ALL
WPOFFS,,,BS/2-DS/2
VSBW,ALL
WPSTYL,DEFA
WPOFFS,,HS/2
WPROTA,,90,
VSBW,ALL
WPOFFS,,,TS
VSBW,ALL
WPOFFS,,,HS-2*TS
VSBW,ALL
WPOFFS,,,TS-0.1
VSBW,ALL

!分割出钢筋位置
C=25
WPSTYL,DEFA
WPOFFS,-BB/2+C,,
WPROTA,,,90
VSBW,ALL
WPSTYL,DEFA
WPOFFS,BB/2-C,,
WPROTA,,,90
VSBW,ALL
WPSTYL,DEFA
WPOFFS,,HB/2-C
WPROTA,,90
VSBW,ALL
WPSTYL,DEFA
WPOFFS,,-HB/2+0.1+C
WPROTA,,90
VSBW,ALL

!切割箍筋
GB=150
WPSTYL,DEFA  
WPOFFS,,,GB/2
VSBW,ALL
*do,I,1,6
WPOFFS,,,GB
VSBW,ALL
*enddo
WPSTYLE,,,,,,,,0
KSEL,S,LOC,Z,0
KSEL,A,LOC,Z,1000
NUMMRG,KP

!分配钢筋的单元特性
!纵向钢筋
LSEL,S,LOC,Z,0,L/2
LSEL,R,LOC,Y,HB/2-C
LSEL,R,LOC,X,-BB/2+C
LATT,5,5,5
ALLSEL,ALL
LSEL,S,LOC,Z,0,L/2
LSEL,R,LOC,Y,HB/2-C
LSEL,R,LOC,X,BB/2-C
LATT,5,5,5
ALLSEL,ALL
LSEL,S,LOC,Z,0,L/2
LSEL,R,LOC,Y,-HB/2+0.1+C
LSEL,R,LOC,X,BB/2-C
LATT,5,5,5
ALLSEL,ALL
LSEL,S,LOC,Z,0,L/2
LSEL,R,LOC,Y,-HB/2+0.1+C
LSEL,R,LOC,X,-BB/2+C
LATT,5,5,5
ALLSEL,ALL

!箍筋
!X=-BB/2+C
LSEL,S,LOC,X,-BB/2+C         
LSEL,R,LOC,Y,-HB/2+C,HB/2-C
LSEL,R,LOC,Z,GB/2
LATT,4,4,4
ALLSEL,ALL
*do,I,1,6
LSEL,S,LOC,X,-BB/2+C         
LSEL,R,LOC,Y,-HB/2+C,HB/2-C
LSEL,R,LOC,Z,GB/2+I*GB
LATT,4,4,4
ALLSEL,ALL
*enddo
!X=BB/2-C
LSEL,S,LOC,X,BB/2-C         
LSEL,R,LOC,Y,-HB/2+C,HB/2-C
LSEL,R,LOC,Z,GB/2
LATT,4,4,4
ALLSEL,ALL
*do,I,1,6
LSEL,S,LOC,X,BB/2-C         
LSEL,R,LOC,Y,-HB/2+C,HB/2-C
LSEL,R,LOC,Z,GB/2+I*GB
LATT,4,4,4
ALLSEL,ALL
*enddo
!Y=-HB/2+C
LSEL,S,LOC,Y,-HB/2+0.1+C      
LSEL,R,LOC,X,-BB/2+C,BB/2-C
LSEL,R,LOC,Z,GB/2
LATT,4,4,4
ALLSEL,ALL
*do,I,1,6
LSEL,S,LOC,Y,-HB/2+0.1+C      
LSEL,R,LOC,X,-BB/2+C,BB/2-C
LSEL,R,LOC,Z,GB/2+I*GB
LATT,4,4,4
ALLSEL,ALL
*enddo
!Y=HB/2-C
LSEL,S,LOC,Y,HB/2-C     
LSEL,R,LOC,X,-BB/2+C,BB/2-C
LSEL,R,LOC,Z,GB/2
LATT,4,4,4
ALLSEL,ALL
*do,I,1,6
LSEL,S,LOC,Y,HB/2-C      
LSEL,R,LOC,X,-BB/2+C,BB/2-C
LSEL,R,LOC,Z,GB/2+I*GB
LATT,4,4,4
ALLSEL,ALL
*enddo
!分配型钢的单元属性
VSEL,S,LOC,X,-BS/2,BS/2
VSEL,R,LOC,Y,HS/2-TS,HS/2
VATT,3,3,3
ALLSEL,ALL
VSEL,S,LOC,X,-TS/2,TS/2
VSEL,R,LOC,Y,-HS/2+TS,HS/2-TS
VATT,3,3,3
ALLSEL,ALL
VSEL,S,LOC,X,-BS/2,BS/2
VSEL,R,LOC,Y,-HS/2+TS,-HS/2+1
VATT,2,2,2
ALLSEL,ALL

!分配混凝土的单元属性
VSEL,S,LOC,X,-BS/2,-HB/2
VSEL,R,LOC,Y,-HB/2+0.1,HB/2
VATT,1,1,1
ALLSEL,ALL
VSEL,S,LOC,X,BS/2,HB/2
VSEL,R,LOC,Y,-HB/2+0.1,HB/2
VATT,1,1,1
ALLSEL,ALL
VSEL,S,LOC,X,-BS/2,-TS/2
VSEL,R,LOC,Y,-HS/2+TS,HS/2-TS
VATT,1,1,1
ALLSEL,ALL
VSEL,S,LOC,X,TS/2,BS/2
VSEL,R,LOC,Y,-HS/2+TS,HS/2-TS
VATT,1,1,1
ALLSEL,ALL
VSEL,S,LOC,X,-BS/2,BS/2
VSEL,R,LOC,Y,HS/2,HB/2
VATT,1,1,1
ALLSEL,ALL
VSEL,S,LOC,X,-BS/2,BS/2
VSEL,R,LOC,Y,-HB/2+0.1,-HS/2
VATT,1,1,1
ALLSEL,ALL

!钢垫板单元属性
VSEL,S,LOC,Y,HB/2,HB/2+40
VATT,6,6,6
ALLSEL,ALL
VSEL,S,LOC,Y,-HB/2-39.9,-HB/2+0.1
VATT,6,6,6
ALLSEL,ALL

!尺寸划分
LSEL,S,LENGTH,,TS     !=8            
LESIZE,ALL,,,1,,,
ALLSEL,ALL
LSEL,S,LENGTH,,TS-0.1   !=7.9
LESIZE,ALL,,,1,,,
ALLSEL,ALL
LSEL,S,LENGTH,,DS       !=6
LESIZE,ALL,,,1,,,
ALLSEL,ALL
!LSEL,S,LENGTH,,12.5    !集中荷载垫板
!LESIZE,ALL,,,1,,,
!ALLSEL,ALL
LSEL,S,LENGTH,,40    !垫板厚度
LESIZE,ALL,,,1,,,
ALLSEL,ALL
LSEL,S,LENGTH,,GB/2   !=75
LESIZE,ALL,,,2,,,
ALLSEL,ALL
LSEL,S,LENGTH,,GB       !箍筋间距=150
LESIZE,ALL,,,4,,,
ALLSEL,ALL
!LSEL,S,LENGTH,,125      !被垫板切割后的箍筋间距
!LESIZE,ALL,,,5,,,
!ALLSEL,ALL
LSEL,S,LENGTH,,C       !=25
LESIZE,ALL,,,1,,,
ALLSEL,ALL
LSEL,S,LENGTH,,HS-2*TS     !=84
LESIZE,ALL,,,4,,,
ALLSEL,ALL
LSEL,S,LENGTH,,(HB-HS-2*C)/2   !=50
LESIZE,ALL,,,2,,,
ALLSEL,ALL
LSEL,S,LENGTH,,(BS-DS)/2       !=47
LESIZE,ALL,,,2,,,
ALLSEL,ALL
LSEL,S,LENGTH,,(BB-BS-2*C)/2     !15
LESIZE,ALL,,,1,,,
ALLSEL,ALL

!网格划分
VSEL,S,TYPE,,6
MSHKEY,1
MSHAPE,0
VMESH,ALL
ALLSEL,ALL
LSEL,S,TYPE,,5
LMESH,ALL
ALLSEL,ALL
LSEL,S,TYPE,,4
LMESH,ALL
ALLSEL,ALL
VSEL,S,TYPE,,3
MSHKEY,1
MSHAPE,0
VMESH,ALL
ALLSEL,ALL
VSEL,S,TYPE,,2
MSHKEY,1
MSHAPE,0
VMESH,ALL
ALLSEL,ALL
VSEL,S,TYPE,,1
MSHKEY,1
MSHAPE,0
VMESH,ALL
ALLSEL,ALL


!创建弹簧单元
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!上翼缘外侧（Y=HS/2）
*DO,I,1,7
ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,Y,HS/2
NSEL,R,LOC,X,-BS/2
TYPE,7
REAL,7
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF

ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,Y,HS/2   
NSEL,R,LOC,X,-(BS-DS)/4-DS/2           
TYPE,7
REAL,8
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF

ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,Y,HS/2
NSEL,R,LOC,X,-DS/2
TYPE,7
REAL,9
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF

ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,Y,HS/2
NSEL,R,LOC,X,DS/2
TYPE,7
REAL,9
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF     
     
      
ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,Y,HS/2
NSEL,R,LOC,X,(BS-DS)/4+DS/2           
TYPE,7
REAL,8
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF

ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,Y,HS/2
NSEL,R,LOC,X,BS/2  
TYPE,7
REAL,7
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF


!上翼缘内侧 Y=HS/2-TS
ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,Y,HS/2-TS
NSEL,R,LOC,X,-BS/2  
TYPE,7
REAL,7
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF

ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,Y,HS/2-TS
NSEL,R,LOC,X,-(BS-DS)/4-DS/2           
TYPE,7
REAL,8
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF

ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,Y,HS/2-TS
NSEL,R,LOC,X,-DS/2
TYPE,7
REAL,10
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF

ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,Y,HS/2-TS
NSEL,R,LOC,X,DS/2
TYPE,7
REAL,10
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF

ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,Y,HS/2-TS
NSEL,R,LOC,X,(BS-DS)/4+DS/2           
TYPE,7
REAL,8
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF

ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,Y,HS/2-TS
NSEL,R,LOC,X,BS/2
TYPE,7
REAL,7
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF

!腹板 X=-DS/2            
ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,X,-DS/2
NSEL,R,LOC,Y,(HS-2*TS)/4
TYPE,7
REAL,11
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF

ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,X,DS/2
NSEL,R,LOC,Y,(HS-2*TS)/4
TYPE,7
REAL,11
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF

ALLSEL,ALL
NSEL,S,LOC,Z ,GB*I-GB/2
NSEL,R,LOC,X,-DS/2
NSEL,R,LOC,Y,0
TYPE,7
REAL,11
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF

ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,X,DS/2
NSEL,R,LOC,Y,0
TYPE,7
REAL,11
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF
ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,X,-DS/2
NSEL,R,LOC,Y,-(HS-2*TS)/4
TYPE,7
REAL,11
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF
ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,X,DS/2
NSEL,R,LOC,Y,-(HS-2*TS)/4
TYPE,7
REAL,11
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF

!下翼缘内侧 Y=-HS/2+0.1+TS
ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,Y,-HS/2+0.1+TS
NSEL,R,LOC,X,-BS/2  
TYPE,7
REAL,7
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF

ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,Y,-HS/2+0.1+TS
NSEL,R,LOC,X,-(BS-DS)/4-DS/2           
type,7
REAL,8
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF

ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,Y,-HS/2+0.1+TS
NSEL,R,LOC,X,-DS/2
TYPE,7
REAL,9
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF

ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,Y,-HS/2+0.1+TS
NSEL,R,LOC,X,DS/2
TYPE,7
REAL,9
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF

ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,Y,-HS/2+0.1+TS
NSEL,R,LOC,X,(BS-DS)/4+DS/2           !18.6667
TYPE,7
REAL,8
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF

ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,Y,-HS/2+0.1+TS
NSEL,R,LOC,X,BS/2
TYPE,7
REAL,7
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF


!下翼缘外侧
ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,Y,-HS/2+0.1
NSEL,R,LOC,X,-BS/2
TYPE,7
REAL,7
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF

ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,Y,-HS/2+0.1
NSEL,R,LOC,X,-(BS-DS)/4-DS/2           !-18.666667
TYPE,7
REAL,8
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF

ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,Y,-HS/2+0.1
NSEL,R,LOC,X,-DS/2
TYPE,7
REAL,9
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF

ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,Y,-HS/2+0.1
NSEL,R,LOC,X,DS/2
TYPE,7
REAL,9
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF

ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,Y,-HS/2+0.1
NSEL,R,LOC,X,(BS-DS)/4+DS/2           !18.666667
TYPE,7
REAL,8
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF

ALLSEL,ALL
NSEL,S,LOC,Z,GB*I-GB/2
NSEL,R,LOC,Y,-HS/2+0.1
NSEL,R,LOC,X,BS/2  
TYPE,7
REAL,7
EINTF
TYPE,8
REAL,100
EINTF
TYPE,9
REAL,101
EINTF
*ENDDO

!施加对称约束
NSEL,S,LOC,Z,0
DSYM,SYMM,Z, ,     !表面垂直与z方向

!施加支座约束
NSEL,S,LOC,Y,-HB/2-39.9
NSEL,R,LOC,Z,900
D,ALL, , , , , ,UX, , , , ,
D,ALL, , , , , ,UY, , , , ,
ALLSEL,ALL
SAVE
/SOL
KBC,1
P=50000         !30Kn
Q=P/HB/100
ASEL,R,LOC,Y,HB/2+40
SFA,ALL,1,PRES,Q
SFTRAN
/PSF,PRES,NORM,2,0,1
ALLSEL,ALL

        
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/SOL
ANTYPE,0           !选择分析类型
OUTPR,BASIC,ALL    !结果输出控制
OUTRES,ALL,ALL     !控制写入数据库和结果文件的结果数据
NLGEOM,ON          !是否计入大变形效应
AUTOTS,ON          !是否使用自动时间步长
PRED,ON       !在非线性分析中是否打开预测器
NROPT,FULL,     !指定牛顿拉夫逊法求解的选项。FULL:完全牛顿拉夫逊法
CNVTOL,F,,0.5,2,1 !力的收敛准则  
CNVTOL,U,,0.5,2,1 !位移的收敛准则
NEQIT,50           !在非线性分析中指定平衡迭代的最大次数
NCNV,2       !终止分析选项。2：如果求解不收敛，终止分析,但不终止程序
DELTIM,0.005,0.005,0.02
ALLSEL,ALL
SOLVE
FINISH

sdaulj

不知道楼主问题解决了没有？我也想做这方面的分析。可否交流一下？我的QQ13940210，期待指点！

LS_Terminator

学习一下，走了

3D主席

请教楼主combin39实常数怎么取的？

xingqiong1986

请问楼主的问题是否解决？我最近在做型钢梁的受弯模拟，也存在不收敛的问题，请教啊