埋入水中旋转薄壳谐耦振分析的传递矩阵法.pdfVIP

第28卷第4期 固体力学学报 V01．28No．4 2007年1．2月 ACTAMECHANICASoLIDASINICA December2007 埋入水中旋转薄壳谐耦振分析的传递矩阵法。 邹时智1 苏海东2 向 宇3 黄玉盈p+ (1华中科技大学土木工程与力学学院，武汉，430074) (2长江科学院材料与结构研究所，武汉，430010)(3广西工学院汽车工程系，柳州，545006) 摘要基于一般线弹性薄壳和势流理论，导出了旋转壳状态向量的一阶常微分矩阵方程和水动压力表达 式，再借助齐次扩容技术和精细积分法，应用推广的传递矩阵法对埋入水中旋转薄壳的流固耦振进行了数值求解， 并研究了一些因素对精度的影响．算例表明传递矩阵法和有限元方法相比不仅精度良好，而且有较高的计算效率． 关键词 旋转薄壳，传递矩阵法，齐次扩容精细积分法，流固耦合振动 0 引言 基于上述目的，本文试图首先推导旋转薄壳状 Tottenham和Shimizu[11提出分析圆柱壳自由态向量的一阶常微分矩阵方程，为传递矩阵法的应 振动的传递矩阵法(或称传递函数法)，后经Irie等 用提供重要条件．然后利用新型的齐次扩容精细积 人不断发展已成功解决了变厚度锥壳、锥柱结合壳、 分技术，研究埋入水中旋转薄壳在谐激励下的耦合 非圆柱壳口“3等各类复杂结构在空气中的谐振问 场的响应问题． 题．它属于一种半解析半数值方法，和数值方法相 1旋转薄壳的基本方程 比，其最大优点是每一壳元的状态向量能够用它的 旋转壳体是由一根平面曲线(称为子午线)绕位 传递矩阵完全地表示，从而使得壳体两端支承处的 于同一平面内的一条直线(轴线)旋转而成，如图1 状态向量问有一个简单关系，求解十分方便，而且计 所示．壳体中曲面上任意点的曲线坐标用(妒，口)表 算精度很高．然而这种方法的前提是壳体的控制方 示，其中，妒代表子午线方向，0代表周向(环向)，相 程必须能写成一阶常微分矩阵方程的形式，这常常 应的主曲率半径分别为R。R。．旋转壳的拉梅系数 限制了该方法的应用． 为R。、Ro—R目sin9． 1989年，黄玉盈等[33导出了圆柱壳元流固耦振 设材料的弹性模量为E，泊松比为卢，密度为lD， 状态向量的一阶常微分矩阵方程，进一步用传递矩阵 壳体厚度为h，它们均可以随子午线坐标∞而变化． 法计算了变厚度圆柱蓄水池的自由振动问题，但由于 以下给出旋转薄壳中面的动力学方程[8]． 水动压力不能事先给出，因此产生了反复迭代的麻 1．1几何方程 烦．之后，他们[5、61借助钟万勰提出的精细时程积分 法口1和扩容技术处理非齐次项，用推广的传递矩阵 如图1所示，设“、口、W分别为叭0及法线}方 向的位移．旋转壳中面应变与位移的关系为 法，又解决了水中任意外形弹性环与纵向加肋无限长 非圆柱壳的声弹耦合问题，这时无需任何迭代，获得 e，一