王新敏《ANSYS工程结构数值分析》p483求解错误。

sdsnjxll

finish
/clear
/config,nres,2000
/prep7
/view,1,1,2,3
et,1,65,,,,,,,1
et,2,8
mp,ex,1,13585
mp,prxy,1,0.2
fc=14.3
ft=1.43
tb,concr,1
tbdata,,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.98
tbpt,,0.002,fc
tbpt,,0.0033,fc*0.85

mp,ex,2,2e5
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
/eshape,1
eplot

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/15000
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
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

sdsnjxll

lidenf

谁有这本书的电子版啊，求~~~

simulin_k

miso不能有下降段

civilman狼人

我找到了源代码，运行没问题，我使用的是ANSYS11.0。下面两个图是计算结果，和书中的结果一致，你试试我附件的命令流，看看情况如何。
  

shifang157300

sdsnjxll 发表于 2012-10-17 19:58

这是由在定义MISO材料曲线时定义了下降段造成的。将MISO曲线丁一部分最后一项改成：
tb,concr,1
tbdata,,0.5,0.95,ft,-1
tb,miso,1,,11
tbpt,,0.0002,fc*0.19
tbpt,,0.0004,fc*0.363
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.911
tbpt,,0.0016,fc*0.961
tbpt,,0.0018,fc*0.98
tbpt,,0.002,fc
tbpt,,0.0033,fc*1.01

n3ch

shifang157300 发表于 2012-11-6 19:31
这是由在定义MISO材料曲线时定义了下降段造成的。将MISO曲线丁一部分最后一项改成：
tb,concr,1
tbdata,,0 ...

遇到了同样的问题，但实际分析时，需要考虑下降段该怎么处理？如果没有下降段，混凝土的本构就不完整了