屈曲分析全过程

1、屈曲分析的过程说明:屈曲分析是一种用于确定结构开始变得不稳定时的临介荷载和屈曲结构发生屈曲响应时的模态形状的技术。ANSYS提供两种结构屈曲荷载和屈曲模态分析方法：非线性屈曲分析和特征值屈曲分析。非线性屈曲分析是在大变形效应开关打开的情况下的一种非线性静力学分析，该分析过程一直进行到结构的极限荷载或最大荷载。非线性屈曲分析的方法是，逐步地施加一个恒定的荷载增量，直到解开始发散为止.尤其重要的是，要一个足够小的荷载增量，来使荷载达到预期的临界屈曲荷载。若荷载增量太大，则屈曲分析所得到的屈曲荷载就可能不准确，在这种情况下打开自动时间步长功能，有助于避免这类问题，打开自动时间步长功能，ANSYS程序

2、将自动寻找屈曲荷载。特征值屈曲分析步骤为：1。建模2 .获得静力解：与一般静力学分析过程一致，但必须激活预应力影响，通常只施加一个单位荷载就行了3 .获得特征屈曲解：Ao进入求解B。定义分析类型Co定义分析选项Do定义荷载步选项E.求解4。扩展解之后就可以察看结果了示例1:!ansys7。0有限元分析实用教程!3.命令流求解!ANSYS命令流：!EigenvalueBucklingFINISH!这两行命令清除当前数据/CLEAR/TITLE,EigenvalueBucklingAnalysis/PREP7!进入前处理器ET,1,BEAM3!选择单元R,1,100,833。333,10!定义实常

3、数MP,EX,1,200000!弹性模量MP,PRXY,1,0.3!泊松比K, 1,0,0!创建梁实体模型L, 2,0,100M, 1,2!创建直线ESIZE,10!单元边长为1mmLMESH,ALL,ALL!划分网格FINISH!退出前处理!屈曲特征值部分/SOLU!进入求解ANTYPE,STATIC从静态分析提取数据PSTRES,ONDK,1,ALLFK,2,FY,-1SOLVEFINISH/SOLUANTYPE,BUCKLEBUCOPT,LANB,1SOLVEFINISH/SOLUEXPASS,ONMXPAND,1SOLVEFINISH!/POST1!在进行屈服分析之前，ANSYS需要!

4、屈服分析中采用预应力!定义约束!顶部施加载荷!求解!退出求解!重新进入求解模型进行屈服分析!屈服分析类型!1阶模态，子空间法!求解!退出求解!重新进入求解展开模态!模态展开打开!定义需要展开的阶数!求解退出求解进入通用后处理SET,LIST!列出特征值求解结果SET,LAST!人感兴趣阶数模态结果PLDISP!显示变形后图形!NonLinearBuckling!非线性分析部分FINISH!这两行命令清除当前数据/CLEAR/TITLE,NonlinearBucklingAnalysis/PREP7!进入前处理ET,1,BEAM3!选择单元MP,EX,1,200000!弹性模量MP,PRXY,1

5、,0.3!泊松比R,1,100,833。333,10!定义实常数K,1,0,0,0!底端节点L,1,2ESIZE,1LMESH,ALLFINISHK,2,0,100,0!顶点连成线!网格尺寸参数设定!划分网格退出前处理/SOLU!进入求解ANTYPE,STATIC!静态分析类型（非屈服分析）NLGEOM,ON!打开非线性大变形设置NSUBST,20NEQIT,1000AUTOTS,ONLNSRCH,ON/ESHAPE,1DK,1,ALL,0OUTRES,ALL,ALL!选择输出数据!5个子步加载!20步迭代!自动时间步长激活线搜索选项!显示二维状态下变形图约束底部节点FK,2,FY,-5000

6、0顶部载荷稍微比特征值分析结果大FK,2,FX,-250!施加水平扰动载荷SOLVE!求解FINISH!退出求解/POST26!进入时间一历程后处理器RFORCE,2,1,F,Y!2#变量表示力NSOL,3,2,U,Y!3#变量表示y方向位移XVAR,2!将x轴显示2#变量PLVAR,3!y轴显示3#变量数据/AXLAB,Y,DEFLECTION!修改y轴标签/AXLAB,X,LOAD!修改x轴标签/REPLOT!重新显示图形示例2:!悬臂梁受端部轴向压力作用的屈曲分析!先进行静力分析，在进行特征值屈曲分析，最后进行非线性分析!静力分析fini/cle/filname,beam-flexure

7、/tittle,beamflexure/prep7*set,f1,-1e6!设置轴向压力荷载参数et,1,beam189mp,dens,1,7,85e3!设置材料参数mp,ex,1,2.06e11mp,nuxy,1,0.2sectype,1,beam,I,2!设置截面参数secoffset,centsecdataQ15,0。15,0。25,0。015,0。015,0.015,0,0,0,0k,1,0k,2,2,5,0k,3,1,25,1lstr,1,2latt,1,1,3,1lesize,1,10lmesh,1/view,1,1,1,1/eshape1.0dk,1,0,all,fk,2,fx,

8、f1finish/soluantype,0eqslv,sparpstres,onsolvefinish!施加关键点压力!求解器设置稀疏矩阵直接法!打开预应力开关!特征值屈曲分析/soluantype,1bucopt,lanb,6,0!取前六阶模态分析mxpand,6,0,0,1,0。001solvefinish/post1*do,i,1,6_set,ipldisp,i*enddo*get,freq1,mode,1,freqfinish!非线性屈曲分析/config,nres,200!只记录两百步的结果/prep7tb,biso,1,1,2!定义材料非线性tbtemp,0tbdata,2.0e8,0upgeom,0。01,1,1,'beam-flexure','rst'!对有限元模型进行一阶模态的位移结果0。01倍的修改save,beam-flexure,dbfinishresu,beamflexure,db/soluantype,0nlgeom,1!打开大变形outres