24小时热门版块排行榜    

查看: 2630  |  回复: 4

feng135yun

新虫 (初入文坛)

[交流] 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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mygt_hit

专家顾问 (职业作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
最好把报错信息也贴出来,有助于帮助找到问题在哪。
知其然,知其所以然。
2楼2015-11-02 10:43:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1983年xiao巷

银虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
我试了一下,提示单元高度扭曲。意思是说某些单元的节点力过大导致了单元高度变形或奇异,试着增加子步数或增加材料的刚度直至能得出结果,然后再检查一下单元的形状
士不可以不弘毅,任重而道远。
3楼2015-11-05 15:03:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
4楼2015-11-06 11:03:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guo624351005

木虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
我也遇到同样问题,楼主怎么解决的?
5楼2016-02-23 11:06:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 feng135yun 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +4 majunge000 2026-08-11 6/300 2026-08-13 16:50 by 四季常青藤
[基金申请] filecode +14 documentary 2026-08-10 16/800 2026-08-13 16:39 by 润生71
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 36/1800 2026-08-12 22:50 by sdfapple719
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 10/500 2026-08-12 09:12 by 布布和一二
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
信息提示
请填处理意见