各位大神好,在下按书上例子用分离式建了一个钢筋混凝土简支梁可以收敛,但是计算结果是错的,命令流如下?
/CONFIG,NRES,2000
/PREP7
ET,1,SOLID65
KEYOPT,1,1,1
KEYOPT,1,3,2
KEYOPT,1,7,1
ET,2,LINK180
ET,3,solid45
MP,EX,1,21945E6
MP,PRXY,1,0.2
FC=23.1E6
FT=1.89E6
TB,CONCR,1
TBDATA,,0.5,0.95,FT,-1
TB,MISO,1,,10
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.0020,FC
MP,EX,2,2.0E11
MP,PRXY,2,0.3
TB,BISO,2
TBDATA,,500E6,0
MP,EX,3,2.0E11
MP,PRXY,3,0.25
TB,BISO,3
TBDATA,,200E6,0
垫板属性
mp,ex,4,2.1e11
mp,prxy,4,0.3
实常数
R,1
R,2,379.94E-6
R,3,50.24E-6
! 建立几何模型
BLC4,,,0.15,0.3,-1.5
wpoff,,0.3,-0.4
blc4,,,0.15,0.05,-0.1
wpcsys,-1
wpoff,,,-1.4
blc4,,,0.15,-0.05,-0.1
wpcsys,-1
wpoff,,,-0.4
vsbw,all
wpoff,,,-0.1
vsbw,all
wpcsys,-1
wpoff,,,-1.4
vsbw,all
wpoff,,,-0.1
vsbw,all
wpcsys,-1
*DO,I,1,10
WPOFF,,,-0.15
VSBW,ALL
*ENDDO
WPCSYS,-1
WPROTA,,90
WPOFF,,,-0.025
VSBW,ALL
WPOFF,,,-0.25
VSBW,ALL
WPROTA,,,90
WPOFF,,,0.025
VSBW,ALL
WPOFF,,,0.1
VSBW,ALL
ALLSEL,ALL
VGLUE,ALL
ELEMSIZ=0.05
LSEL,S,LOC,X,0.025
LSEL,R,LOC,Y,0.025
CM,ZJA,LINE
LATT,2,2,2
LESIZE,ALL,ELEMSIZ
LSEL,S,LOC,X,0.125
LSEL,R,LOC,Y,0.025
CM,ZJB,LINE
LATT,2,2,2
LESIZE,ALL,ELEMSIZ
LSEL,S,LOC,X,0.025
LSEL,R,LOC,Y,0.275
CM,ZJC,LINE
LATT,2,2,2
LESIZE,ALL,ELEMSIZ
LSEL,S,LOC,X,0.125
LSEL,R,LOC,Y,0.275
CM,ZJD,LINE
LATT,2,2,2
LESIZE,ALL,ELEMSIZ
LSEL,S,TAN1,Z
LSEL,R,LOC,Y,0.025,0.275
LSEL,R,LOC,X,0.025,0.125
LSEL,U,LOC,Y,0.35
LSEL,U,LOC,Z,-0.4
LSEL,U,LOC,Z,-0.5
LSEL,U,LOC,Y,-0.05
LSEL,U,LOC,Z,-1.4
CM,GJA,LINE
LATT,3,3,2
LESIZE,ALL,ELEMSIZ
CMSEL,S,ZJA
CMSEL,A,ZJB
CMSEL,A,ZJC
CMSEL,A,ZJD
CMSEL,A,GJA
LMESH,ALL
LSEL,ALL
/ESHAPE,1
EPLOT
VSEL,S,LOC,Y,0,0.3
VATT,1,1,1
MSHKEY,1
ESIZE,ELEMSIZ
VMESH,ALL
ALLSEL,ALL
VSEL,S,LOC,Y,0.3,0.35
VATT,4,1,3
MSHKEY,1
ESIZE,ELEMSIZ
VMESH,ALL
ALLSEL,ALL
VSEL,S,LOC,Y,-0.05,0
VATT,4,1,3
MSHKEY,1
ESIZE,ELEMSIZ
VMESH,ALL
ALLSEL,ALL
ASEL,S,LOC,Z,0
DA,ALL,SYMM
LSEL,S,LOC,Y,-0.05
LSEL,R,LOC,Z,-1.45
DL,ALL,,UY
FLST,2,1,1,ORDE,1
FITEM,2,1254
F,P51X,FY,-100000
ALLSEL,ALL
/solu
antype,0
NLGEOM,1
NROPT,FULL
outres,all,all
autots,on
nsubst,200
neqit,100
cnvtol,u,,0.02
solve