24小时热门版块排行榜

返回列表

feng135yun

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 0.5
帖子: 2
在线: 22分钟
虫号: 3729999
注册: 2015-03-11

[交流] ansys求帮助已有4人参与

谁有CFRP钢管混凝土命令流或者帮忙检查下我的命令流哪里错了
finish
/clear
/PREP7

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义材料属性
et,1,181
KEYOPT,1,1,2
r,1,0.168
mp,ex,1,2.3e5
MP,PRXY,1,0.3
tb,bkin,1,1,2,1
tbdata,,3550

et,2,65
mp,ex,2,3.0e4       !混凝土弹模
mp,prxy,2,0.2          !混凝土泊松比
tb,conc,2             !定义非线性材料特性
tbdata,,0.5,0.9,2.55,-1,,  !张开系数，闭合，单轴抗拉强度，单轴抗压
tb,miso,2,1,12          !多线性等向强化
tbpt,,0.0001,3.0    !应变应力关系
tbpt,,0.0003,7.076
tbpt,,0.0005,11.156
tbpt,,0.0007,14.726
tbpt,,0.0009,17.786
tbpt,,0.0011,20.336
tbpt,,0.0013,22.376
tbpt,,0.0015,23.906
tbpt,,0.0017,24.926
tbpt,,0.002,25.5
tbpt,,0.0022,25.5
tbpt,,0.0023,25.5

et,3,45
MP,EX,3,2.06E5
MP,PRXY,3,0.3
TB,bkin,3,1,2, 1       !双线性随动强化模型
TBDATA,,235

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!设置参数区

R1=73
R2=75
bb=4/3
zzz=300

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!建立模型

csys,1
k,1, R1/bb,0
k,2, R1/bb,90
k,3, R1/bb,180
k,4, R1/bb,270
k,5, R1,0
k,6, R1,90
k,7, R1,180
k,8, R1,270
k,9, R2,0
k,10, R2,90
k,11, R2,180
k,12, R2,270
csys,0
l,1,2
l,2,3
l,3,4
l,4,1
csys,1
l,5,6
l,6,7
l,7,8
l,8,5

l,9,10
l,10,11
l,11,12
l,12,9

csys,0
l,1,5
l,2,6
l,3,7
l,4,8
l,5,9
l,6,10
l,7,11
l,8,12

al,1,2,3,4
al,1,13,5,14
al,2,14,6,15
al,3,15,7,16
al,4,16,8,13

al,17,9,18,5
al,18,10,19,6
al,19,11,20,7
al,20,12,17,8
!
alls
csys,1
lsel,s,loc,y,45
lplot
lesize,all,,,4,
!
alls
csys,1
lsel,s,loc,y,135
lplot
lesize,all,,,4,
!
alls
csys,1
lsel,s,loc,y,225
lplot
lesize,all,,,4,
!
alls
csys,1
lsel,s,loc,y,315
lplot
lesize,all,,,4,
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
LSEL,S,LOC,Y,0
LSEL,A,LOC,Y,90
LSEL,A,LOC,Y,180
LSEL,A,LOC,Y,270
LPLOT
LSEL,U,LOC,X,R1,R2
LPLOT
lesize,all,,,2,

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
LSEL,S,LOC,Y,0
LSEL,A,LOC,Y,90
LSEL,A,LOC,Y,180
LSEL,A,LOC,Y,270
LPLOT
LSEL,U,LOC,X,0,R1
LPLOT
lesize,all,,,1
!
alls
csys,1
asel,s,loc,x,0,r1-7
aplot
type,  1
amesh ,all
alls
asel,u,loc,x,0,r1-7
aplot
type,1
amesh,all
alls

type,2
mat,2
real,2
asel,s,loc,x,0,r1-7

extopt,esize,zzz/20
EXTOPT,ACLEAR,1
EXTOPT,ATTR,0,0,0
csys,0

VEXT,all, , ,0,0,zzz,,,,
eplot

alls
csys,1

type,3
mat,3
real,3
asel,u,loc,x,0,r1-7
asel,u,loc,z,1,zzz

extopt,esize,zzz/20
EXTOPT,ACLEAR,1
EXTOPT,ATTR,0,0,0
csys,0
VEXT,all, , ,0,0,zzz,,,,
alls
eplot

   !接触单元设置
et,4,170
et,5,173
keyopt,5,9,0
keyopt,5,12,2

r,5
real,5
mp,mu,2,0.2

csys,1
vsel,s,type,,3
asel,r,loc,x,r1
aplot
cm,target,area
type,4
                           !设置目标单元
nsla,s,1
esurf,all
esel,s,type,,4
eplot
/psyms,esys,1
allsel
vsel,s,type,,2
asel,r,loc,x,r1
cm,contact,area
type,5
nsla,s,1
esurf,all
esel,s,type,,5
eplot
/psyms,esys,1
allsel
gplot

csys,1
type,1
mat,1
real,1
*do,i,0,15
*do,z,0,zzz/20-1
node1=node(r2,i*90/4,z*20)
node2=node(r2,i*90/4,(z+1)*20)
node4=node(r2,(i+1)*90/4,z*20)
node3=node(r2,(i+1)*90/4,(z+1)*20)
e,node1,node2,node3,node4
*enddo
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!! 求解

NSEL,S,LOC,Z,0
NPLOT
D,all, , , , , ,ALL, , , , ,
ALLSEL,ALL
EPLOT
!————————
/solu
antype,static
outres,basic,all       !输出选项
outpr,basic,all       !窗口输出，这个没什么影响
autots,on             !自动计算时间步长
nlgeom,on          !打开是大变形，非线性
outres,all,all    !输出所有子步结果
time,1
allsel

NSEL,S,LOC,Z,zzz
NPLOT
D,all, ,-20, , , ,UZ, , , , ,
nsubst,100    !荷载子步
neqit,20
eqslv,-1
ncnv,1             !终止分析选项，终止分析且退出程序
lnsrch,on          !线性搜索，使计算稳定
CNVTOL,F, ,0.05,2, ,
allsel
solve

/POST26
numvar,10
csys,1
nsel,s,loc,z,zzz
*get,ncount,node,,count    !获取被选择节点的数量定义到变量ncount中
nnow=0
nsol,6,ndnext(nnow),u,z,Zdisplancement
add,6,6,,,,,,-1
*do,i,1,ncount
nnow=ndnext(nnow)
rforce,8,nnow,f,z,zforce
add,7,7,8,,,,,,1
*enddo
xvar,6
/axlab,x,deflection(mm)
/axlab,y,totle load(N)
prod,7,,,,,,,-1
plvar,7

回复此楼