资源描述：
《beam188不能用于计算单轴对称截面梁的弯扭失稳问题》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、beam188不能用于计算单轴对称截面梁的弯扭失稳问题？题：用beam188单元求单轴对称H型截面梁在纯弯作用下的线性屈曲特征值。题目条件：截面高度：300mm上翼缘：150*12mm下翼缘：80*12mm腹板厚度：10mm构件长度：3000mm弹性模量：E=68000MPa泊松比：0.315两端铰接，简支，端部可自由翘曲打开Beam188的翘曲自由度计算结果：Mcr=50.17kN*m根据经典弹性理论，βy=97.25mm（正值，由于上翼缘较大，受压），Mcr=85.93kN*m结果明显错误。于是将上下翼缘颠倒
2、，再计算之。即：上翼缘：80*12mm下翼缘：150*12mm计算结果仍然是：Mcr=50.17kN*m根据经典弹性理论，βy=-97.25mm（负值，由于上翼缘较小，受压），Mcr=29.22kN*m最后，经典弹性理论的计算公式中，直接取βy=0.0mm（不考虑Wagnar效应），可得Mcr=50.11kN*m，这样才和Ansys计算结果相近。结论：beam188不能用于计算单轴对称截面梁的弯扭失稳问题。（也许我还没有找到某个开关，先暂时下此结论。望有高手指教）再以板单元建立模型验证：单元采用SHELL63，上
3、翼缘大时，得Mcr=50.12kN*m颠倒过来，上翼缘小时，得Mcr=29.20kN*m结论：板单元可以用于计算单轴对称截面梁的弯扭失稳问题。也再次验证了beam188的计算错误。2007-12-1804:30#1warsheep助理工程师精华0积分70帖子33水位70技术分0忘了说明：Ansys版本为8.0支座位置位于截面形心上。2007-12-1810:12#2wilsonweic助理工程师精华0积分70帖子35水位70技术分0用beam188/189单元进行线性弯扭屈曲分析，结果不可靠。除了你所说的单轴对称
4、截面外，事实上，双轴对称截面梁的线性弯扭屈曲分析结果也不准确。虽然所涉及的只是线性屈曲分析，但是，我们有理由怀疑，beam188/189非线性弯扭屈曲分析的结果也不一定准确。条件允许的话，最好用壳单元分析结果进行校核。我分别用beam188单元和shell181单元对两端简支和两端固支的双轴对称截面纯弯梁进行了特征值屈曲分析，分析结果与理论屈曲荷载进行了对比。用beam188分析时，分别考虑了KEYOPT1=0和KEYOPT1=1两种情况。结果显示shell181分析的结果与理论值吻合良好，而beam188分析的
5、结果误差很大（见附件）。相关的内容还可参见http://baisi.net/viewthread...mp;bpg=1&age=302007-12-1818:26#3mqiao助理工程师精华0积分69帖子34水位69技术分0您说的有些道理，但部分观点不敢苟同。经过对ansys的beam188单元刚度矩阵的研究，发现原因可能是beam188的几何刚度矩阵一项中只考虑了轴力的影响，而没有考虑和转角φ以及φ'所对应的弯矩和剪力的影响，梁单元的完整的几何刚度矩阵可参考《钢结构构件稳定理论于设计》（陈骥）p445。让人郁闷
6、的是，Ansys手册中没有给出beam188的几何刚度矩阵公式，只给出了一个beam4的几何刚度矩阵公式。（AnsystheoryReference,式（3-60）），所以上述也仅仅是猜测而已。但是，可以确定的是：beam188不能用于计算梁（主要是非双轴对称截面梁）的弯扭失稳，包括线性特征值分析和几何非线性分析。另外经过我对大量算例的分析表明：对于双轴对称截面梁的弯扭屈曲，其计算结果和经典解接近。（但长度很小时差别略大）对于轴心受压构件的弯曲屈曲，其计算结果和经典解接近。对于轴压构件的弯扭屈曲，其计算结果和经典
7、解接近。其他有限元软件也有类似问题。很不幸，我要算的梁太多了，如果全部用板单元恐怕要算到明年了。只能自己编程序了。郁闷中啊。jxln兄，我的以上认识可能也有不妥，有空可相互探讨之，我的qq是1788382。2007-12-1900:03#4710476703工程师精华0积分136帖子38水位79技术分0jxln兄：我粗略研究了你的附件中的命令流，发现您的命令流其实是错误的。您的命令流为：两端简支beam188.c修改后的命令流为（修改处增加了注解）：/prep7et,1,beam188,,,2,keyopt,1
8、,1,1!*******考虑翘曲自由度********mp,ex,1,2.06e11mp,prxy,1,0.3SECTYPE,1,BEAM,I,,2!*******原有命令流截面划分过于粗糙********SECOFFSET,CENTSECDATA,0.1,0.1,0.6,0.01,0.01,0.008,0,0,0,0secplot,1,1k,1,k,2,5,k,3,,-
