求解结构动应力的振型函数差分法Ξ.PDFVIP

第 15 卷第 2 期 计 算 力 学 学 报 V o l. 15 N o. 2 1998 年 5 月 CH IN ESE JOU RNAL O F COM PU TA T IONAL M ECHAN ICS M ay 1997 求解结构动应力的振型函数差分法 　　　　赖永星　　王　伟　　　　　　张汴生　　　　　　 ( ) ( ) 郑州工业大学数力系, 郑州, 450002 　　　 河南省科技干部学校, 郑州, 450003 摘　要　本文利用常规有限元方法的计算结果, 结合数值计算方法对振型函数进行 [L ] 算子的微分 计 算, 从而可方便迅速获得到复杂结构动应力响应, 并对梁和板进行了计算, 计算结果表明该方法具 有较高的精度, 较一般的动态有限元具有通用性强, 计算简单等特点。 关键词　动应力; 有限元; 有限差分 分类号　 3471 O 1　引　　言 对于工程设计人员而言, 不仅关心结构在动荷作用下的位移响应, 更为关心它的动应力响 应, 如高层建筑受到地震的作用, 海洋平台受到海浪的冲击, 桥梁受到车辆的振动等等, 因此动 应力的分析计算一直是结构分析中的一个重要课题, 然而至今结构动应力还是一个有待进一 ( 步研究的问题, 有很多不完善的地方。由于结构与方程的复杂性 如复杂的边界形状、边界条件 ) 等 一般很难求得这些动力问题的解析解, 为此人们提出了各种近似计算方法, 如动态有限元, [ 1 ] 但动态有限元局限性较大, 而且计算量宏大 , 很难推广应用。另外在通常的有限元方法中, 求 解结构动态特性与动力位移响应时, 都是直接采用有限元离散建立系统的刚度矩阵 [K ] 和质 量矩阵[M ], 其单元形函数是逼近结点位移条件下的静态位移场的幂函数, 这种单元的静态位 移场是唯一的, 然而当结构振动时其动态位移场是由结构的各阶主振动线性迭加而得出, 与单 元静态位移场有很大差别, 用一般的有限元方法来计算结构的动特性和动位移响应, 还可得到 [ 2 ] 工程要求的精度, 但是用于求结构动应力其精度很低 。本文对此进行了理论分析与计算, 找 出了其造成误差的原因, 并根据有限元方法计算动特性具有较高精度的特点, 利用差分法求出 了具有较高精度的动应力。 2　用有限元方法计算结构动应力及分析 与静力分析相同, 场变量采用有限元插值 ( ) ( ) ( ) {u } = [N x , y , z ]{u e t } 1 其中, [N ] 为形函数矩阵。 对离散单元而言, 由有限元方法可知 ( ) ( ) ( ) {} = [L ][N ]{u e t } = [B ]{ue t } 2 其中 [B ] = [L ][N ] 为几何矩阵, [L ] 为微分算子矩阵。 河南省自然科学基金资助项目 　收稿日期: 修改稿收到日期: 　赖永星: 男, 1960 年生, 硕士, 讲