在做加了预应力的非线性分析的时候会出现一些问题在后处理/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