一个简单的屈曲分析的实例，比较有代表性

zhafeiyang

这个实例虽然比较简单，但是整个的分析过程很具有代表性
对与掌握分析屈曲过程的步骤很有帮助，前两天刚做的
贴出来共享

!悬臂梁受端部轴向压力作用的屈曲分析
!先进行静力分析，在进行特征值屈曲分析，最后进行非线性分析

!静力分析
/prep7                                                 !
*set,f1,-1e6                                           !设置轴向压力荷载参数
et,1,beam189                                           !
sectype,1,beam,I,,2                                    !设置截面参数
secoffset,cent                                         !
secdata,0.15,0.15,0.25,0.015,0.015,0.015,0,0,0,0       !
mp,dens,1,7.85e3                                       !设置材料参数
mp,ex,1,2.06e11                                        !
mp,nuxy,1,0.2                                          !
k,1,0                                                  !
k,2,2.5,0                                              !
k,3,1.25,1                                             !
lstr,1,2                                               !
latt,1,,1,,3,,1                                        !
lesize,1,,,10                                          !
lmesh,1                                                !
/view,1,1,1,1                                          !
/eshape,1.0                                            !
dk,1,,,,0,all,                                         !
fk,2,fx,f1                                             !施加关键点压力
finish                                                 !
                                                       !
/solu                                                  !
antype,0                                               !
eqslv,spar                                             !求解器设置  稀疏矩阵直接法
postres,on                                             !打开预应力开关
solve                                                  !
finish                                                 !

!特征值屈曲分析                                                   
/solu                                                  !
antype,1                                               !
bucopt,lanb,6,0                                        !取前六阶模态分析
mxpand,6,0,0,1,0.001                                   !
solve                                                  !
finish                                                 !
                                                       !
/post1                                                 !
set,first                                              !
pldisp,2                                               !
set,next                                               !
pldisp,2                                               !
set,next                                               !
pldisp,2                                               !
set,next                                               !
pldisp,2                                               !
set,next                                               !
pldisp,2                                               !
set,next                                               !
pldisp,2                                               !
*get,freq1,mode,1,freq                                 !
finish                                                 !

!非线性屈曲分析
/config,nres,200                                       !只记录两百步的结果                                                       !
/prep7                                                 !
tb,biso,1,1,2                                          !定义材料非线性
tbtemp,0                                               !
tbdata,,2.0e8,0                                        !
upgeom,0.01,1,1,'beam-flexure','rst'                   !对有限元模型进行一阶模态的位移结果0.01倍的修改
finish                                                 !
                                                       !
/solu                                                  !
antype,0                                               !
nlgeom,1                                               !打开大变形
outres,all,all                                         !
arclen,1,0                                             !弧长法设置
arctrm,l                                               !弧长法终止准则 达到第一个峰值是终止计算
nsubst,200,,,1                                         !
fk,2,fx,f1*freq1                                       !将轴向压力值放大，放大系数为第一阶模态的主频
solve                                                  !
finish                                                 !
                                                       !
/post26                                                !
nsol,2,2,u,y,deflection                                !提取自由端y方向的位移为变量deflection
rforce,3,1,f,x,reactionf                               !提取固定端x方向的支座反力为变量reactionf
/axlab,x,deflection                                    !
/axlab,y,reactionf                                     !
xvar,2                                                 !
plvar,3                                                !
finish                                                 !

zhafeiyang

模型图

zhafeiyang

当计算达到第一个峰值时停止计算
（计算达到第一个峰值时，因为应变较小，支座反力也达到最大值）

WANGBOWEI2007

不错，搂住有心阿

dengbin_ren

!将轴向压力值放大，放大系数为第一阶模态的主频
为什么这样做？

zhafeiyang

因为第一阶模态的阶特征值
也就是结构屈曲时初始荷载的比例系数
而结构到达第一阶屈曲时特征值相对其他各阶是最大的
因此用它乘以初始荷载就得到了结构的临界荷载

dengbin_ren

谢谢

xiaowei331

楼主  不知道你有没有《ansys结构有限元高级分析方法与范例应用》这本书没？
就是水利水电出版社出的   我想问一下   如果按你这种流程   这书最后面的那个单层网壳的例子算出来的非线性屈曲荷载要比特征值屈曲荷载肯定要大 而且我都难以猜测出究竟大到什么程度   能指点高招吗  怎么求单层网桥的非线性屈曲荷载   谢谢了

liusha1057

不打开PSTRES是无法得到特征值的

太感人

zhafeiyang 发表于 2007-6-19 15:50
当计算达到第一个峰值时停止计算
（计算达到第一个峰值时，因为应变较小，支座反力也达到最大值）
...

为啥没有下降段呢