径向变厚度圆板在均布压力作用下的屈曲.pdfVIP

·688- ‘工程力学)增刊1999年 径向变厚度圆板在均布压力作用下的屈曲 向天宇郑建军周欣竹 (北方交通大学土木毒筑工程学院。北京1000“) 候典杰 邵林 (山东建筑工程学院教务赴．滔博255012)(浙江工业大学土木工程学院。抗州3l∞14) 曩一 本文讨论了变厚度圆板的屈曲．用Bessel函数理论导出了等厚度圆扳单元和环板单元传 递矩阵的正确公式，再应用传递矩阵原理建立了整体变厚度既板届曲计算的特征方程．该法具有计算简 单、节约内存的优点· 关■曩 变厚度四板，屈曲，传递矩阵法 1前言 在计算变厚度圆板屈曲荷载时，一般采用有限元法。为了得到较为精确的结果，必须 细分单元，这导致计算代价成倍增长。传递矩阵法属于一种半解析数值方法，它具有计算 精确、未知数少的优点[1翻。本文应用传递矩阵法研究了变厚度圆板屈曲问题并给出了相应 数值结果，表明了该方法的有效性，具有一定的应用价值。 、 2传递矩阵的推导 对于如图1所示的半径为n的受径向压力只作用的径向变厚度圆板，我们可以证明径 向压力只(，)和环向压力弓(r)满足关系 ’ 尸，(，)=P．(，)=Pf (1) 为了求得临界压力￡，我们可以将圆板划分成如图2所示的月一1个环板单元和1个圆 板单元，对于每个单元我们可以视为等厚度板．根据文献【3】的结果，最低阶失稳为轴对称 失稳，失稳时板的位移满足微分方程 日v4∥+只v2∥=O (2) 上式中，V2=参+；菩，皿是第f个板单元的弯曲刚度·令 q=丽 (3) 向天宇．男。1972．6出生，硕士研究生 ‘工程力学)增刊1999年 ·689‘ 根据Bessel函数理论，方程(2)的解为 ∥=G+Qlllr+GJo(口，)+G虼(口，) (4) 对于圆板单元而言，考虑板中心处挠度和弯矩应为有限值，其解应为 w1=cl。+G。山■。r) (5) 下面给出了这两种情况下相应的径向转角日，、径向弯矩珥和径向剪力一的表达式。 圈l径向变厚度圈扳受径向压力作用 圈2单元划分 (1)圆板单兀的啡、M和I的表达式，即 ’ 钟=一G-^忙1，) 叫=一c3。D。a。3-，。忙，，) 卅喝州如卜等砸。叫 @ (2)环板单元的啡、^‘和一的表达式，即 、 只=一q，诉^忙i，)+争一q≯，墨忙。r)_=一c。qa，’‘(a，)一G。日a。3《扛，) 设第一个圆板单元的半径等于，o-单元周边上的状态向量{f1)为 {r1)={w：口：。叫．心)’ (8) ‘工程力学》增刊1999年 将，=，0代入式(5)、(6)得 ㈦=【，。】{c} (9) 其