我做一个悬索的找形分析的练习时候总提示我自由度约束不足。练习先是按照悬列线建模,再迭代。我是把悬索用很多根link10单元连起来进行分析,还有一种是用bsplin函数连成一条线分析,但是连成一条样条曲线时候划分单元时候变成一条直线,最后还是报错。我参考的ansys工程结构数值分析这本书里面找形,都是先建一link单元再找形,我想问一下如果建了悬索以及吊杆情况下,最后找形完成后吊杆下端标高还是保持原来的标高么?以及可以帮忙看一下错误的原因么?十分感谢!
代码如下:
finish
/clear
/filename,xuansuo
/prep7
antype,0
et,1,link10
et,2,beam3
r,1,1E-5
r,2,0.02
r,3,0.5,1/24,1
mp,ex,1,25e10
mp,dens,1,1000
mp,prxy,1,0.1
mp,ex,2,30e10
mp,dens,2,1000
mp,prxy,2,0.1
*DO,i,1,11,1
x=(i-1)*10
Y=250*COSH(X/250)-250
K,I,X,Y
*ENDDO
*DO,I,1,11,1
X=(I-1)*10
Y=0
K,I+11,X,Y
*ENDDO
*DO,I,1,10,1
L,I,I+1
*ENDDO
*DO,I,1,9,1
L,I+1,I+12
*ENDDO
*DO,I,1,10,1
L,I+11,I+12
*ENDDO
LSEL,S,LINE,,1,11
LATT,1,2,1
LSEL,S,LINE,,11,19
LATT,1,2,1
LSEL,S,LINE,,9,29
LATT,2,3,2
lsel,all
lesize,all,,,1
lmesh,all
NSYM,X,22,all
ESYM,,22,all
NUMMRG,ALL
NUMCMP,ALL
finish
/SOLU
nsel,s,loc,x,-100
nsel,r,loc,y,20.268
d,all,Ux,,,,,Uy
nsel,s,loc,x,-100
nsel,r,loc,y,0
d,all,Uy
nsel,s,loc,x,100
nsel,r,loc,y,20.268
d,all,Ux,,,,,Uy
nsel,s,loc,x,100
nsel,r,loc,y,0
d,all,Ux,,,,,Uy
nsel,s,loc,y,0
nsel,r,loc,x,100
nsel,a,loc,x,-100
nsel,r,loc,y,0
nsel,inve
d,all,Uz
allsel,all
acel,0,-10,0
nropt,auto,,on
esel,s,elem,,1,10
esel,a,elem,,30,39
esel,inve
ekill,all
allsel,all
nlgeom,1
nsubst,10,20,5
time,1
sstif,on
solve
save
finish
/post1
ESEL,S,LIVE
pldisp,1
finish
/prep7
upgeom,1,last,last,xuansuo,rst
finish
save
代码如下:
finish
/clear
/filename,xuansuo
/prep7
antype,0
et,1,link10
et,2,beam3
r,1,1E-5
r,2,0.02
r,3,0.5,1/24,1
mp,ex,1,25e10
mp,dens,1,1000
mp,prxy,1,0.1
mp,ex,2,30e10
mp,dens,2,1000
mp,prxy,2,0.1
*DO,i,1,11,1
x=(i-1)*10
Y=250*COSH(X/250)-250
K,I,X,Y
*ENDDO
*DO,I,1,11,1
X=(I-1)*10
Y=0
K,I+11,X,Y
*ENDDO
*DO,I,1,10,1
L,I,I+1
*ENDDO
*DO,I,1,9,1
L,I+1,I+12
*ENDDO
*DO,I,1,10,1
L,I+11,I+12
*ENDDO
LSEL,S,LINE,,1,11
LATT,1,2,1
LSEL,S,LINE,,11,19
LATT,1,2,1
LSEL,S,LINE,,9,29
LATT,2,3,2
lsel,all
lesize,all,,,1
lmesh,all
NSYM,X,22,all
ESYM,,22,all
NUMMRG,ALL
NUMCMP,ALL
finish
/SOLU
nsel,s,loc,x,-100
nsel,r,loc,y,20.268
d,all,Ux,,,,,Uy
nsel,s,loc,x,-100
nsel,r,loc,y,0
d,all,Uy
nsel,s,loc,x,100
nsel,r,loc,y,20.268
d,all,Ux,,,,,Uy
nsel,s,loc,x,100
nsel,r,loc,y,0
d,all,Ux,,,,,Uy
nsel,s,loc,y,0
nsel,r,loc,x,100
nsel,a,loc,x,-100
nsel,r,loc,y,0
nsel,inve
d,all,Uz
allsel,all
acel,0,-10,0
nropt,auto,,on
esel,s,elem,,1,10
esel,a,elem,,30,39
esel,inve
ekill,all
allsel,all
nlgeom,1
nsubst,10,20,5
time,1
sstif,on
solve
save
finish
/post1
ESEL,S,LIVE
pldisp,1
finish
/prep7
upgeom,1,last,last,xuansuo,rst
finish
save