ANSYS混凝土柱滞回为什么曲线是直线，是本构的问题还是模型的问题？

/CLEAR

/PREP7

AS0=615.44$AS1=78.5

ET,1,SOLID65

KEYOPT,1,1,0                          !考虑大变形

KEYOPT,1,5,1                          !线性解输出控制：给出每个积分点的解

KEYOPT,1,6,3                          !非线性解输出控制：给出每个积分点的解

KEYOPT,1,7,1                          !开裂后考虑应力松弛，有助于计算收敛

MP,EX,1,2.4E4$MP,PRXY,1,0.2$FC=25

TB,CONCR,1,1,9$TBDATA,,0.35,0.75,3.1125,-1

TB,MISO,1,,15

TBPT,,0.0002,4.8$ TBPT,,0.0004,9.375$ TBPT,,0.0006,13.51

TBPT,,0.0008,17.02$ TBPT,,0.001,19.83$ TBPT,,0.0012,21.95

TBPT,,0.0014,23.43$ TBPT,,0.0016,24.365$ TBPT,,0.0018,24.856

TBPT,,0.002,25$ TBPT,,0.0038,25$TBPLOT

MP,EX,2,2E5$MP,PRXY,2,0.25$TB,BISO,2$TBDATA,,360

MP,EX,3,2E5$MP,PRXY,3,0.25$TB,BISO,3$TBDATA,,210

 

 

 

!体积配筋率

R,1,3,AS1/(70*100),0,0,2, AS0/(70*70)$RMORE,90,0,3, AS1/(70*100),0,90

R,2,3, AS1/(70*100),0,90

R,3,3, AS1/(70*100),0,0

R,4

!实常数最多只能有六项，后六项需添加RMORE

BLC4,,,500,1500,500

WPOFFS,,,70$VSBW,ALL

WPOFFS,,,360$VSBW,ALL

WPROTA,,,90$WPOFFS,,,70$VSBW,ALL

WPOFFS,,,360$VSBW,ALL

 

VSEL,S,LOC,X,0,70$VSEL,R,LOC,Z,0,70

VATT,1,1,1$MSHAPE,0,3D$MSHKEY,1$LESIZE,ALL,70$VMESH,ALL$ALLSEL

VSEL,S,LOC,X,430,500$VSEL,R,LOC,Z,0,70

VATT,1,1,1$MSHAPE,0,3D$MSHKEY,1$LESIZE,ALL,70$VMESH,ALL$ALLSEL

VSEL,S,LOC,X,0,70$VSEL,R,LOC,Z,430,500

VATT,1,1,1$MSHAPE,0,3D$MSHKEY,1$LESIZE,ALL,70$VMESH,ALL$ALLSEL

VSEL,S,LOC,X,430,500$VSEL,R,LOC,Z, 430,500

VATT,1,1,1$MSHAPE,0,3D$MSHKEY,1$LESIZE,ALL,70$VMESH,ALL$ALLSEL

VSEL,S,LOC,X,0,70$VSEL,R,LOC,Z, 70,430

VATT,1,2,1$MSHAPE,0,3D$MSHKEY,1$LESIZE,ALL,120$VMESH,ALL$ALLSEL

VSEL,S,LOC,X,430,500$VSEL,R,LOC,Z, 70,430

VATT,1,2,1$MSHAPE,0,3D$MSHKEY,1$LESIZE,ALL,120$VMESH,ALL$ALLSEL

VSEL,S,LOC,X,70,430$VSEL,R,LOC,Z, 0,70

VATT,1,3,1$MSHAPE,0,3D$MSHKEY,1$LESIZE,ALL,120$VMESH,ALL$ALLSEL

VSEL,S,LOC,X,70,430$VSEL,R,LOC,Z, 430,500

VATT,1,3,1$MSHAPE,0,3D$MSHKEY,1$LESIZE,ALL,120$VMESH,ALL$ALLSEL

VSEL,S,LOC,X,70,430$VSEL,R,LOC,Z, 70,430

VATT,1,4,1$MSHAPE,0,3D$MSHKEY,1$LESIZE,ALL,120$VMESH,ALL$ALLSEL

 

!/DEVICE,VECTOR,1

!$/ESHAPE,1.0

!$/REPLOT

 

NSEL,S,LOC,Y,0

D,ALL,ALL

 

/SOLU

ANTYPE,0

NLGEOM,ON

OUTRES,ALL,ALL

AUTOTS,OFF

CNVTOL,F,,0.1,2,0.9

 

TIME,1$NSUBST,20

NSEL,S,LOC,Y,1500

D,ALL,UX,1

ALLSEL

SOLVE

 

TIME,2$NSUBST,40

NSEL,S,LOC,Y,1500

D,ALL,UX,-1

ALLSEL

SOLVE

 

TIME,3$NSUBST,40

NSEL,S,LOC,Y,1500

D,ALL,UX,2

ALLSEL

SOLVE

/POST26

NSOL,2,491,U,X

RFORCE,3,491,F,X

XVAR,2

PLVAR,3

PRVAR,2,3