ANSYS柱子稳定分析算例

!ANSYS柱子稳定分析算例!ExampleofgeometrynonlinearanalysisinANSYS!要点：预应力，特征值屈曲，添加初始缺陷，几何非线性分析!作者:陆新征，清华大学土木系!Author:LuXinzhengDept.CivilEngrg.ofTsinghuaUniversity!lastrevised:2003.2.!finish/CLEAR/UNITS,SI/PREP7!*FORCE=100OFFSET=0.1!初始缺陷为0.1!建立模型ET,1,BEAM4ET,2,LINK10R,1,0.1*0.12,0.12*0.1**3/12,0.1*0.12**3/12,0.12,0.1,,R,2,0.002*0.002,2e-3,!预应力MPTEMP,1,0MPDATA,EX,1,,200e9MPDATA,PRXY,1,,0.27MPDATA,DENS,1,,7800k,1,k,2,0,0,5k,3,0,0,-5k,11,0.2,0,0k,12,-0.2,0,0k,13,0,0.2k,15,0,-0.2l,1,2l,1,3l,1,11l,1,12l,1,13l,1,15l,2,11l,2,12l,2,13l,2,15l,3,11l,3,12l,3,13l,3,15lsel,,,,1,6latt,1,1,1ALLSEL,ALLlsel,,,,7,14latt,1,2,2ALLSEL,ALLlsel,,,,1,6LESIZE,all,0.3,,,,,,,1lsel,,,,7,14LESIZE,all,,,1,,,,,1ALLSEL,ALLLMESH,ALLFINISH/ESHAPE,1.0!*!求解特征值屈曲荷载*DO,I,1,100FINISH/SOLUDK,3,,,,0,UX,UY,UZ,,,ROTZ,DK,2,,,,0,UX,UY,,,,,FK,2,FZ,-FORCEANTYPE,0!设定时间步TIME,1AUTOTS,0NSUBST,1,,,1!NLGEOM,1SSTIF,ONSOLVEFINISH/SOLUANTYPE,BUCKLE!BucklinganalysisBUCOPT,LANB,1!UseBlockLanczossolutionmethod,extract1modeMXPAND,1!Expand1modeshapePSTRES,ON!INCLUDEPRESTRESSEFFECTSSOLVEFINISH*GET,FREQ1,MODE,1,FREQ*IF,FREQ1,LT,0.01,THENFREQ1=0.01*ENDIF*IF,ABS(FREQ1-1),LT,0.01,THEN!如果频率误差小于1%*EXIT*ENDIFFORCE=FORCE*FREQ1*ENDDO!添加初始缺陷FINISH/POST1SET,FIRSTNSORT,U,SUM*GET,DMAX,SORT,,MAX!*FINISH/prep7upgeom,OFFSET/DMAX,1,1,beammodel,rst!非线性求解FINISH/SOLU!*FORCE=FORCE*3.DK,3,,,,0,UX,UY,UZ,,,ROTZ,DK,2,,,,0,UX,UY,,,,,FK,2,FZ,-FORCEANTYPE,0NLGEOM,1NSUBST,20,0,0OUTRES,ALL,1ARCLEN,1,0,0AUTOTS,-1.0SOLVE!绘制荷载位移曲线FINISH/POST26FILEFILE,'beammodel','rst','.'/UI,COLL,1NUMVAR,200FILLDATA,191,,,,1,1REALVAR,191,191NUMVAR,200FILLDATA,191,,,,1,1REALVAR,191,191!*NSOL,2,2,U,Z,UZ_2STORE,MERGE!*RFORCE,3,19,F,Z,FZ_3STORE,MERGEXVAR,2PLVAR,3,