土木在线论坛 \ 建筑结构 \ ANSYS \ 结构非线性屈曲分析的问题

结构非线性屈曲分析的问题

发布于:2011-04-05 08:47:05 来自:建筑结构/ANSYS [复制转发]
在做加了预应力的非线性分析的时候会出现一些问题
在后处理/post26中,以位移为X轴,以时间为Y轴显示图形,但是我只有一个荷载步,也就是说时间是从0—1之间,但是最后显示出来的时间Y轴,却到了-50多,跪求高手解答下!!!






/filname,,120
/prep7



csys,1
*do,i,1,48,1$k,i,20,360/48*I,0$*enddo
*do,i,1,47,1$lstr,i,i+1$*enddo
lstr,48,1
*do,i,49,88,1$k,i,17,360/40*(I-48),3.330$*enddo
*do,i,49,87,1$lstr,i,I+1$*enddo
lstr,88,49
*do,i,89,120,1$k,i,14,360/32*(I-88),5.71$*enddo
*do,i,89,119,1$lstr,i,I+1$*enddo
lstr,120,89
*do,i,121,144,1$k,i,11,360/24*(I-120),7.45$*enddo
*do,i,121,143,1$lstr,i,I+1$*enddo
lstr,144,121
*do,i,145,160,1$k,i,8,360/16*(I-144),8.69$*enddo
*do,i,145,159,1$lstr,i,I+1$*enddo
l,160,145
*do,i,161,168,1$k,i,4,360/8*(I-160),9.68$*enddo
*do,i,161,167,1$lstr,i,I+1$*enddo
lstr,168,161
k,169,0,0,10
csys,0
l,48,88
l,120,88
l,144,120
l,144,160
l,160,168
l,168,169
l,1,49
l,49,89
l,89,121
l,121,145
l,145,161
l,2,50
l,50,90
l,122,90
l,122,146
l,3,51
l,51,91
l,91,123
l,4,52
l,52,92
l,5,53
l,1,88
l,49,120
l,89,144
l,121,160
l,145,168
l,2,49
l,50,89
l,90,121
l,122,145
l,3,50
l,51,90
l,91,122
l,4,51
l,52,91
l,5,52
csys,1
lgen,8,all,,,,45 x
nummrg,all
ldele,1588
ldele,1384
ldele,1180
ldele,976
ldele,772
ldele,568
ldele,364
ldele,160!!!!!!!!!!!!!!!!!!!!!!!表面网壳单元!!!!!!!!都用beam188!!!!!!!


numstr,line,2001
*do,i,170,209,1$k,i,17,360/40*(I-169),-3.25$*enddo
*do,i,49,88$lstr,i,i+121$*enddo

*do,i,210,241,1$k,i,14,360/32*(I-209),-1.25$*enddo
*do,i,89,120$lstr,i,i+121$*enddo

*do,i,242,265,1$k,i,11,360/24*(I-241),3.5$*enddo
*do,i,121,144,1$lstr,i,i+121$*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!都用beam188!!!!!!!!下部撑杆!!!!!



numstr,line,3001
*do,i,170,208,1$lstr,i,i+1$*enddo
lstr,209,170

*do,i,210,240,1$lstr,i,i+1$*enddo
lstr,241,210

*do,i,242,264,1$lstr,i,i+1$*enddo
lstr,265,242
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!use link10!!!!!!!!!!!!!下部环索!!!!!




csys,0
numstr,line,4001
l,48,209
l,88,241
l,120,265


l,1,209
l,49,241
l,89,265


l,1,170
l,49,210
l,89,242


l,2,170
l,50,210
l,90,242


l,2,171
l,50,211
l,90,243


l,3,172
l,51,212
l,91,244

l,3,171
l,51,211
l,91,243

l,4,172
l,52,212

l,4,173
l,52,213

l,5,173
l,5,174
csys,1
lgen,8,4001,4030,1,,45!基本模型建立
nummrg,all
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!11下部拉杆!!!!用beam188!!!!



!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!暂时不用!!numstr,line,5001
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!暂时不用!!csys,1
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!暂时不用!!*do,i,501,548,1$k,i,20,360/48*(I-500),0$*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!暂时不用!!*do,i,401,448,1$k,i,20,360/48*(I-400),-10$*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!暂时不用!!*do,i,501,548,1$lstr,i,i-100$*enddo



!R1=35/1000$R2=76/1000!赋予属性
!PI=acos(-1)
!A1=PI*R1*R1/4$A2=PI*R2*R2/4


et,1,link10$et,2,beam188

A1=2809e-6

sectype,10,beam,ctube$secdate,0.085,0.09,8

mp,ex,1,1.95e8
mp,prxy,1,0.3
mp,dens,1,7850



r,1,A1,2.98e-3


lsel,s,line,,1,1632,1$$latt,1,,2,,,,10

lsel,s,line,,2001,3000,1$latt,1,,2,,,,10

lsel,s,line,,4001,5000,1$$latt,1,,2,,,,10

lsel,s,line,,3001,4000,1$latt,1,1,1!下部拉索单元



allsel,all
fkdele,all,all
ksel,s,kp,,49,169,1
fk,all,fz,-1625.704


ksel,s,kp,,1,48,1
dk,all,all








allsel,all
lsel,all
lesize,all,,,1
lmesh,all!划分网格













/solu
antype,0
eqslv,spar
pstres,on
solve
finish





/solu
antype,1!求屈曲荷载
bucopt,lanb,1

solve
finish

/post1
set,list
finish


/prep7
upgeom,0.1,1,1,120,rst!初始缺陷



!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!非线性屈曲!!
/solu
antype,0
nlgeom,1
arclen,1

nsubst,100
outres,all,all




fkdele,all,all
ksel,s,kp,,49,169,1
fk,all,fz,-1625.704*1.2





solve
finish




/post26
nsol,2,169,u,z
prod,4,2,,,,,,-1,1,1
prod,5,1,,,,,,1,1,1




xvar,4
plvar,5

全部回复(2 )

只看楼主 我来说两句抢地板
  • whumaomao
    whumaomao 沙发
    没什么奇怪的,你只设置了一个荷载步,但是程序默认有一个按钮会在一些情况下自动计算荷载步,所以出现了50.
    2011-05-15 20:48:15

    回复 举报
    赞同0
  • xuehailangzi
    xuehailangzi 板凳
    这个还是有点难度的,要好好学习才是
    2011-05-12 13:31:12

    回复 举报
    赞同0
这个家伙什么也没有留下。。。

ANSYS

返回版块

4.4 万条内容 · 144 人订阅

猜你喜欢

阅读下一篇

预应力混凝土箱梁应力集中问题

我做的模型,在预应力钢筋锚固点处总是有应力集中,请问各位应该怎么处理?

回帖成功

经验值 +10