【分享】悬索结构的分析的例题

ouyangll

下面是一道用ANSYS进行悬索结构的找形和计算的例题  
取自<悬索结构设计>沈世钊，中国建筑工业出版设，1997.8  
p149例题3-1  
进行了找形计算和荷载作用下的静力计算  
找形时取较小的弹性模量，静力计算时采用弹性模量的真值  
   
这个例题说明，ANSYS还是可以进行悬索结构的分析的  
  
/prep7  
ET,1,LINK10  
!找形时采用很小的弹性模量  
R,1,0.001468,0.9,  
MP,EX,1,0.60551e9  
MP,NUXY,1,0.3  
!定义节点  
*do,i,1,5  
*do,j,1,2*i-1  
k=(i-1)*(i-1)+j  
n,k,(j-i)*9.15,-36.6+(i-1)*9.15,0  
*enddo  
*enddo  
*do,i,1,4  
*do,j,1,2*i-1  
k=(i-1)*(i-1)+j+25  
n,k,(j-i)*9.15,36.6-(i-1)*9.15,0  
*enddo  
*enddo  
!定义单元  
*do,i,1,4  
*do,j,1,2*i  
k=i*i+j  
e,k,k+1  
*enddo  
*enddo  
*do,i,1,3  
*do,j,1,2*i  
k=i*i+j+25  
e,k,k+1  
*enddo  
*enddo  
*do,i,1,4  
*do,j,1,2*i-1  
k=(i-1)*(i-1)+j  
e,k,k+2*i  
*enddo  
*enddo  
*do,i,1,3  
*do,j,1,2*i-1  
k=(i-1)*(i-1)+j+25  
e,k,k+2*i  
*enddo  
*enddo  
*do,i,18,24  
e,i,i+17  
*enddo  
!施加位移约束  
d,1,ux,0  
d,1,uy,0  
d,1,uz,3.66  
d,2,ux,0  
d,2,uy,0  
d,2,uz,1.83  
d,4,ux,0  
d,4,uy,0  
d,4,uz,1.83  
d,5,ux,0  
d,5,uy,0  
d,5,uz,0  
d,9,ux,0  
d,9,uy,0  
d,9,uz,0  
d,10,ux,0  
d,10,uy,0  
d,10,uz,-1.83  
d,16,ux,0  
d,16,uy,0  
d,16,uz,-1.83  
d,17,ux,0  
d,17,uy,0  
d,17,uz,-3.66  
d,25,ux,0  
d,25,uy,0  
d,25,uz,-3.66  
d,35,ux,0  
d,35,uy,0  
d,35,uz,-1.83  
d,41,ux,0  
d,41,uy,0  
d,41,uz,-1.83  
d,30,ux,0  
d,30,uy,0  
d,30,uz,0  
d,34,ux,0  
d,34,uy,0  
d,34,uz,0  
d,27,ux,0  
d,27,uy,0  
d,27,uz,1.83  
d,29,ux,0  
d,29,uy,0  
d,29,uz,1.83  
d,26,ux,0  
d,26,uy,0  
d,26,uz,3.66  
!设定时间步  
TIME,1  
AUTOTS,0  
NSUBST,20, , ,1  
KBC,0  
!输出每个子步的结果  
OUTPR,BASIC,ALL,  
OUTRES,ALL,ALL,  
!设定大变形，应力刚化  
NLGEOM,1  
SSTIF,ON  
FINISH  
!求截  
/SOLU  
/STATUS,SOLU  
SOLVE  
FINISH  
!后处理  
/POST1  
PRNSOL,DOF,  
PRESOL,SMISC,1  
PLESOL,SMISC,1  
PLNSOL,U,X,0,1  
finish  
/PREP7  
!更新几何形状  
!a.rst为计算结果文件名，最后一个为目录  
!这两个参数应根据你的计算情况定  
UPGEOM,1,LAST,LAST,a,rst,E:\JZD\1\  
!弹性模量恢复为真值  
R,1,0.001468,0.0027248,  
MP,EX,1,2.0e11  
!重新施加位移约束  
d,1,uz,0  
d,2,uz,0  
d,4,uz,0  
d,5,uz,0  
d,9,uz,0  
d,10,uz,0  
d,16,uz,0  
d,17,uz,0  
d,25,uz,0  
d,35,uz,0  
d,41,uz,0  
d,30,uz,0  
d,34,uz,0  
d,27,uz,0  
d,29,uz,0  
d,26,uz,0  
!求截  
/SOLU  
/STATUS,SOLU  
SOLVE  
FINISH  
/POST1  
!观察施加荷载之前结构的位移  
!可见位移基本为零,预张力基本不变  
PLNSOL,U,Z,0,1  
PLESOL,SMISC,1  
FINISH  
!写荷载工况文件01  
/PREP7  
LSWRITE,01,  
!施加节点荷载  
*do,i,1,41  
f,i,fz,-167445  
*enddo  
!设置第2荷载步  
TIME,2  
AUTOTS,0  
NSUBST,20, , ,1  
KBC,0  
LSWRITE,02,  
!求解  
FINISH  
/SOLU  
LSSOLVE,1,2,1,  
后处理  
/POST26  
!可以看到，该结果和书中结果一致  
NSOL,2,21,U,Z,  
PLVAR,2, , , , , , , , , ,  
ESOL,3,33,1,F,Y,  
PLVAR,3, , , , , , , , , ,

FTYL

不错，支持！

foransys

good, thank you!!

cchin

good!

menyueqi

很好呀

dsy3227

上述命令流，实际上进行了悬索的非线性分析，通过此分析，结构在重力状态下达到平衡，在此状态下，如何进行模态分析？