结构非线性屈曲分析的问题
zenfanqiu
zenfanqiu Lv.2
2011年04月05日 08:47:32
来自于ANSYS
只看楼主

在做加了预应力的非线性分析的时候会出现一些问题在后处理/post26中,以位移为X轴,以时间为Y轴显示图形,但是我只有一个荷载步,也就是说时间是从0—1之间,但是最后显示出来的时间Y轴,却到了-50多,跪求高手解答下!!!/filname,,120/prep7csys,1*do,i,1,48,1$k,i,20,360/48*I,0$*enddo*do,i,1,47,1$lstr,i,i+1$*enddo

在做加了预应力的非线性分析的时候会出现一些问题
在后处理/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
免费打赏
xuehailangzi
2011年05月12日 13:31:01
2楼
这个还是有点难度的,要好好学习才是
回复
whumaomao
2011年05月15日 20:48:51
3楼
没什么奇怪的,你只设置了一个荷载步,但是程序默认有一个按钮会在一些情况下自动计算荷载步,所以出现了50.
回复

相关推荐

APP内打开