王思敏教材中钢混梁无法计算，提示本构关系曲线斜率小于0，求大神解释。
浏览：92715 回答：3
!EXAMPLE 钢混简支梁2
FINISH$/CLEAR$/CONFIG,NRES,2000$/PREP7
!定义单元与材料性质
ET,1,SOLID65,,,,,,,1
ET,2,LINK180
MP,EX,1,13585
MP,PRXY,1,0.2
FC=14.3
FT=1.43
TB,CONC,1
TBDATA,1,0.5,0.95,FT,-1
TB,MISO,1,,11
TBPT,,0.0002,FC*0.19
TBPT,,0.0004,FC*0.36
TBPT,,0.0006,FC*0.51
TBPT,,0.0008,FC*0.64
TBPT,,0.0010,FC*0.75
TBPT,,0.0012,FC*0.84
TBPT,,0.0014,FC*0.91
TBPT,,0.0016,FC*0.96
TBPT,,0.0018,FC*0.99
TBPT,,0.002,FC
TBPT,,0.0033,FC*0.85
MP,EX,2,2.0E5
MP,PRXY,2,0.3
TB,BISO,2
TBDATA,,300,0
PI=ACOS(-1)
R,1,0.25*PI*22*22
R,2,0.25*PI*22*22/2
R,3,0.25*PI*10*10
R,4,0.25*PI*10*10/2
!创建几何模型并切分
BLC4,,,150/2,300,2000/2
*DO,I,1,9
WPOFF,,,100
VSBW,ALL
*ENDDO
WPCSYS,-1
WPOFF,,,50
VSBW,ALL
WPCSYS,-1
WPROTA,,-90
WPOFF,,,30
VSBW,ALL
WPOFF,,,240
VSBW,ALL
WPCSYS,-1
WPOFF,30
WPROTA,,,90
VSBW,ALL
WPCSYS,-1
!划分钢筋网格
ELEMSIZ=50
LSEL,S,LOC,X,30
LSEL,R,LOC,Y,30
CM,ZJ,LINE
LATT,2,1,2
LESIZE,ALL,ELEMSIZ
LSEL,S,LOC,X,75
LSEL,R,LOC,Y,30
CM,ZJB,LINE
LATT,2,2,2
LESIZE,ALL,ELEMSIZ
LSEL,S,LOC,X,30
LSEL,R,LOC,Y,270
CM,JLJ,LINE
LATT,2,3,2
LESIZE,ALL,ELEMSIZ
LSEL,S,TAN1,Z
LSEL,R,LOC,Y,30,270
LSEL,R,LOC,X,30,70
LSEL,U,LOC,Z,50
CM,GJ,LINE
LATT,2,3,2
LESIZE,ALL,ELEMSIZ
LSEL,S,LOC,Z,0
LSEL,R,LOC,Y,30,270
LSEL,R,LOC,X,30,70
CM,GJB,LINE
LATT,2,4,2
LESIZE,ALL,ELEMSIZ
LSEL,ALL
CMSEL,S,ZJ
CMSEL,A,ZJB
CMSEL,A,JLJ
CMSEL,A,GJ
CMSEL,A,GJB
CM,GJ,LINE
LMESH,ALL
LSEL,ALL
!划分混凝土单元网格
VATT,1,,1
MSHKEY,1
ESIZE,ELEMSIZ
VMESH,ALL
ALLSEL,ALL
!施加荷载和约束
LSEL,S,LOC,Y,0
LSEL,R,LOC,Z,900
DL,ALL,,UY
ASEL,S,LOC,Z,0
DA,ALL,,SYMM
ASEL,S,LOC,X,75
DA,ALL,,SYMM
P0=180000
Q0=P0/150/100
ASEL,S,LOC,Z,0,50
ASEL,R,LOC,Y,300
SFA,ALL,1,PRES,Q0
ALLSEL,ALL
!求解控制设置
/SOLU$ANTYPE,0$NSUBST,100$OUTRES,ALL,ALL$AUTOS,ON
NEQIT,50
CNVTOL,U,,0.015
SOLVE
!进入POST1查看结果
/POST1$SET,LAST$PLDISP,1
ESEL,S,TYPE,,2
ETABLE,SAXL,LS,1$PLLS,SAXL,SAXL
ESEL,S,TYPE,,1
/DEVICE,VECTOR,ON$PLCRACK
!进入时程后处理
/POST26
NSOL,2,205,U,Y$PROD,3,2,,,,,,-1
PROD,4,1,,,,,,P0/1000
XVAR,3$PLVAR,4
FINISH$/CLEAR$/CONFIG,NRES,2000$/PREP7
!定义单元与材料性质
ET,1,SOLID65,,,,,,,1
ET,2,LINK180
MP,EX,1,13585
MP,PRXY,1,0.2
FC=14.3
FT=1.43
TB,CONC,1
TBDATA,1,0.5,0.95,FT,-1
TB,MISO,1,,11
TBPT,,0.0002,FC*0.19
TBPT,,0.0004,FC*0.36
TBPT,,0.0006,FC*0.51
TBPT,,0.0008,FC*0.64
TBPT,,0.0010,FC*0.75
TBPT,,0.0012,FC*0.84
TBPT,,0.0014,FC*0.91
TBPT,,0.0016,FC*0.96
TBPT,,0.0018,FC*0.99
TBPT,,0.002,FC
TBPT,,0.0033,FC*0.85
MP,EX,2,2.0E5
MP,PRXY,2,0.3
TB,BISO,2
TBDATA,,300,0
PI=ACOS(-1)
R,1,0.25*PI*22*22
R,2,0.25*PI*22*22/2
R,3,0.25*PI*10*10
R,4,0.25*PI*10*10/2
!创建几何模型并切分
BLC4,,,150/2,300,2000/2
*DO,I,1,9
WPOFF,,,100
VSBW,ALL
*ENDDO
WPCSYS,-1
WPOFF,,,50
VSBW,ALL
WPCSYS,-1
WPROTA,,-90
WPOFF,,,30
VSBW,ALL
WPOFF,,,240
VSBW,ALL
WPCSYS,-1
WPOFF,30
WPROTA,,,90
VSBW,ALL
WPCSYS,-1
!划分钢筋网格
ELEMSIZ=50
LSEL,S,LOC,X,30
LSEL,R,LOC,Y,30
CM,ZJ,LINE
LATT,2,1,2
LESIZE,ALL,ELEMSIZ
LSEL,S,LOC,X,75
LSEL,R,LOC,Y,30
CM,ZJB,LINE
LATT,2,2,2
LESIZE,ALL,ELEMSIZ
LSEL,S,LOC,X,30
LSEL,R,LOC,Y,270
CM,JLJ,LINE
LATT,2,3,2
LESIZE,ALL,ELEMSIZ
LSEL,S,TAN1,Z
LSEL,R,LOC,Y,30,270
LSEL,R,LOC,X,30,70
LSEL,U,LOC,Z,50
CM,GJ,LINE
LATT,2,3,2
LESIZE,ALL,ELEMSIZ
LSEL,S,LOC,Z,0
LSEL,R,LOC,Y,30,270
LSEL,R,LOC,X,30,70
CM,GJB,LINE
LATT,2,4,2
LESIZE,ALL,ELEMSIZ
LSEL,ALL
CMSEL,S,ZJ
CMSEL,A,ZJB
CMSEL,A,JLJ
CMSEL,A,GJ
CMSEL,A,GJB
CM,GJ,LINE
LMESH,ALL
LSEL,ALL
!划分混凝土单元网格
VATT,1,,1
MSHKEY,1
ESIZE,ELEMSIZ
VMESH,ALL
ALLSEL,ALL
!施加荷载和约束
LSEL,S,LOC,Y,0
LSEL,R,LOC,Z,900
DL,ALL,,UY
ASEL,S,LOC,Z,0
DA,ALL,,SYMM
ASEL,S,LOC,X,75
DA,ALL,,SYMM
P0=180000
Q0=P0/150/100
ASEL,S,LOC,Z,0,50
ASEL,R,LOC,Y,300
SFA,ALL,1,PRES,Q0
ALLSEL,ALL
!求解控制设置
/SOLU$ANTYPE,0$NSUBST,100$OUTRES,ALL,ALL$AUTOS,ON
NEQIT,50
CNVTOL,U,,0.015
SOLVE
!进入POST1查看结果
/POST1$SET,LAST$PLDISP,1
ESEL,S,TYPE,,2
ETABLE,SAXL,LS,1$PLLS,SAXL,SAXL
ESEL,S,TYPE,,1
/DEVICE,VECTOR,ON$PLCRACK
!进入时程后处理
/POST26
NSOL,2,205,U,Y$PROD,3,2,,,,,,-1
PROD,4,1,,,,,,P0/1000
XVAR,3$PLVAR,4
邹正刚 
