应用ANSYS对圆薄板的屈曲响应理论计算结果的验证.pdfVIP

应用ANSYS对圆薄板的屈曲响应理论计算结果的验证 蒋 泉 丁华建 南通大学 建筑工程学院 （226007 ） 摘要： 本文利用有限元软件对圆形薄板的屈曲问题进行了计算，在不同的计算单元模式下，计算 结果差异很大。在选取合适的单元形式时，计算结果与理论分析一致。结果表明，常用的板和 壳单元对于薄板的屈曲问题有着一定的局限性。在进行在有限元软件计算时，选择合适的单元 非常重要。 关键词：有限元 薄板 屈曲问题 Abstract: The finite element software is used to solve the buckling problem of circular thin plate. Various results arise under the different element modes. The numerical results will match the theoretical ones when the proper element mode is selected. It is concluded that the plate and shell mode is fiendish to solve the buckling problems. Keywords: Finite element; Thin plate; Buckling problem 一、前言 随着科学技术的发展，板壳已广泛应用于各种工程结构中。因此，其稳定性的 研究已成为一个重要课题。屈曲作为板壳失稳的一种表现形式，已经为很多科学家 ［ － ］ 及工程师所关注 1 2 。关于薄圆板在面内压缩力作用下的轴对称屈曲问题早已被不 少人所研究，近年来由于对镀层及层合材料脱层问题力学机理的研究，屈曲及其后 ［ ］ 屈曲理论又有了新的应用 3 。有些文献[4-6]利用摄动—有限元方法研究了椭圆板的 弹性失稳。或者利用势能原理和非线性有限元法, 分析了具有中心弹簧支承的极正 交各向异性圆板的轴对称稳定性。 有限元软件由于集成化和易用性受到了人们越来越多的关注，在工程应用和验 证理论结果方面有着无可替代的地位。本文对圆板轴对称后屈曲问题利用不同单元 进行了有限元计算，并且将结果与理论解相比较。说明了计算单元对薄板屈曲问题 的影响，对其他类型的稳定性问题有着一定的借鉴意义。 二、一般圆板轴对称后屈曲的理论推导 设圆形板的圆周半径为 ，厚度为 。用r 表示板中面上的点离开板中心的距 R h 离，ϖ 表示板屈曲以后中面的挠度。 圆板在面内力作用下发生轴对称屈曲以后，随边界径向力的增加，后屈曲变形 增大，同时在环向产生很高的压应力σθ ，在这个力作用下，圆板会在环向产生非轴 1 对称第二次分叉屈曲。以下，由挠度ϖ 和应力函数F 表示的极坐标系中的一对 Karman非线性板方程导出控制圆板二次分叉屈曲和后屈曲的方程，将挠度w 和应力 函数F 写成无量纲形式 ϖ w ϖ +ϖ （1a） h 0 1 f F 2 f 0 +f 1 （1b） Eh 其中，带有下标0的量对应二次分叉屈曲前轴对称屈曲变形的量，带有下标1的量表 示二次分叉屈曲后的增量。二次分叉屈曲后变形的非线性方程为 ［7］ 2 2 i