命令流如下：

FINISH

/CLEAR

/PREP7

!PARAMETER OF CONCRETE

*SET,EC,31218.5

*SET,FC,26.8

*SET,FT,2.39

!PARAMETER OF FIBER

*SET,EF,72000

*SET,FFT,1600

!DEFINE ELEMENT TYPE

!CONCRETE

ET,1,SOLID65

KEYOPT,1,1,1

KEYOPT,1,3,2

KEYOPT,1,7,1

!FIBER

ET,2,LINK8

!CONNECTION BETWEEN CONCRETE AND FIBER

ET,3,COMBIN39

KEYOPT,3,1,0

KEYOPT,3,3,0

!DEFINE MATERIAL MODELS

!CONCRETE

MP,EX,1,EC

MP,PRXY,1,0.2

TB,CONCR,1                        !定义混凝土材料参数，关闭压碎

TBDATA,,0.4,0.95,FT,-1

TB,MISO,1,,11                     !定义混凝土应力-应变曲线

TBPT,,0.1*0.0015888,4.96          !应力应变曲线点

TBPT,,0.2*0.0015888,9.26

TBPT,,0.3*0.0015888,13.5

TBPT,,0.4*0.0015888,17.3

TBPT,,0.5*0.0015888,20.5

TBPT,,0.6*0.0015888,22.9

TBPT,,0.7*0.0015888,24.8

TBPT,,0.8*0.0015888,26

TBPT,,0.9*0.0015888,26.6

TBPT,,0.0015888,26.8

TBPT,,2.492*0.0015888,26.8

!FIBER

MP,EX,2,EF

MP,PRXY,2,0.3

!DEFINE REAL CONSTANTS

R,1

R,2,0.975

R,3,0,0,0.46,11.69*7,1.84,13.15*7

RMORE,4.35,7.22*7,5,7.22*7

!PARAMETERS OF SPECIMENS

*SET,A,20    !FIBER LENGTH IN CONCRETE

*SET,A1,20

*SET,D,5

*SET,H,15    !SPECIMEN THCKNESS

*SET,L,200    !SPECIMEN TOTAL LENGTH

*SET,B,80    !SPECIMEN WIDTH

!MODEL CREATING

BLOCK,0,B,0,A1+A,0,H

BLOCK,0,B,A1+A+D,A1+A+D+L,0,H

WPOFF,B/2,0,0

WPROTA,0,0,90

VSBW,ALL

WPCSYS

WPOFF,0,A1,0

WPROTA,0,90,0

VSBW,ALL

WPCSYS

WPOFF,0,0,H/2

VSBW,ALL

WPCSYS

!MESH THE CONCRETE

ALLSEL

LESIZE,ALL,2

LSEL,ALL

MSHKEY,1

VMESH,ALL

!MESH THE FIBER

LSEL,NONE

K,50,B/2,A1,H/2

K,51,B/2,A1+A,H/2

K,52,B/2,A1+A+D,H/2

K,53,B/2,A1+A+D+L,H/2

L,50,51

L,51,52

L,52,53

LESIZE,ALL,2

LATT,2,2,2

LMESH,ALL

ALLSEL

!DEFINE THE COMBINE ELEMENT

CPINTF,UX

*GET,CPMAX,CP,,MAX

TYPE,3

REAL,3

*DO,I,1,CPMAX

*GET,N1,CP,I,TERM,1,NODE

*GET,N2,CP,I,TERM,2,NODE

E,N1,N2

*ENDDO

CPINTF,UY

*GET,CPMAX,CP,,MAX

TYPE,3

REAL,3

*DO,I,1,CPMAX

*GET,N1,CP,I,TERM,1,NODE

*GET,N2,CP,I,TERM,2,NODE

E,N1,N2

*ENDDO

CPINTF,UZ

*GET,CPMAX,CP,,MAX

TYPE,3

REAL,3

*DO,I,1,CPMAX

*GET,N1,CP,I,TERM,1,NODE

*GET,N2,CP,I,TERM,2,NODE

E,N1,N2

*ENDDO

!DEFINE THE DISPLACEMENT

ASEL,S,LOC,Y,A1+A+D+L                 !选取对称面

DA,ALL,SYMM                     !施加对称约束

!DEFINE THE LOAD

ALLSEL

ASEL,S,LOC,Y,0

SFA,ALL,,PRES,750/(H*B)

FINISH

/SOL

ANTYPE,0                       !分析类型为静力分析

NSUBST,50                      ! 荷载步50

OUTRES,ALL,ALL

AUTOS,ON

NEQIT,50                       ! 迭代次数50

CNVTOL,U,,0.05                  !收敛误差为5%

SOLVE                           !开始求解

FINISH

/POST26

NSOL,2,2206,U,X                ! 选择变量

PROD,3,2,,,,,,-1

PROD,4,1,,,,,,750

XVAR,3                           !横轴变量

PLVAR,4                          !绘制变量

FINISH

处理过程中，出现错误！错误命令见图片。求大神指导！为何 UZ方向是非受力方向，加载后却位移超过限制呢？

QQ图片20140704091034.jpg