旋转柔性梁动态特性的影响分析.docxVIP

旋转柔性梁动态特性的影响分析 1 旋转柔性梁的动力学分析方法 在旋转柔性梁的动力学设计中，基于传统传统刚性假设的设计方法无法正确预测振动性能，因此需要计算纵向变形与旋转运动之间的耦合容量。旋转柔性梁的应力刚化与弹性变形的耦合, 使系统的动态特性更加复杂, 呈现了明显的非线性。文献[1]采用哈密顿原理推导了旋转梁的非线性纵向和弯曲耦合方程, 运用Galerkin方法和多尺度法求得了方程的一阶解析解。文献[2]利用数值迭代方法求解非线性方程, 分析了柔性梁在旋转运动和弹性变形耦合状态下的振动特性。文献[3]基于Kane式建立了旋转梁一般形式的有限元动力学方程, 利用梁截面参数导出了旋转梁的显式三维梁单元矩阵, 最后计算了旋转梁的固有频率并与常规计算进行了比较。除此之外, 很多学者对旋转柔性梁的动态特性展开研究, 做出了显著贡献 在结构动力学中, 商用软件梁单元模型 (如ANSYS, NASTRAN等) 无法同时设置应力刚化及旋转软化效应, 不能用来分析旋转柔性梁的动态特性, 必须采用实体单元模型进行分析 本文采用非线性应变位移关系表达式, 通过变分公式推导出简化的梁单元模型;并编制有限元程序并分析了柔性梁的动力学特性。 2 基本原则 2.1 u梁应力张量的u3000方程 根据文献[6]建立旋转柔性梁的模型, 见图1。在图示平面 (x-z平面) , 梁绕着过端点且垂直该平面的y轴转动, 转动角ψ (弧度) ;u 梁上给定点在变形前后的位置矢量X、X分别为 式中:0≤x≤L, -h/2≤z≤h/2;φ=v′;L为梁的长度;h为梁的高度; 位置列矢量为 总位移矢量为 应变张量为 设定轴向位移u 梁的势能为 式中A、I分别为梁横截面的面积、惯性矩。 按旋转柔性梁在图示平面内的纵向位移u和横向位移v, 将动能T表示为 式中 根据哈密顿原理, 有 可得转动梁的非线性运动方程为 式中表明了惯性力、离心力、旋转软化、加速度、横向位移、弯曲刚度等作用效应。 在推导有限元运动方程的过程中, 对上式简化得运动方程为 在有限元模型中, 忽略式中对振动模态影响较小的陀螺耦合项 2.2 旋转力学性能 由杆及梁单元形函数、势能及动能, 根据哈密顿原理, 可导出旋转梁有限元运动方程为 梁单元轴向、横向对称刚度矩阵K 梁单元科氏力矩阵C 梁单元反对称加速度刚度矩阵K 梁单元旋转离心力刚度矩阵K 旋转软化矩阵K 由式 (17) 可见, 随着旋转速度增大, 旋转软化矩阵呈现非线性特性。参考文献[17-19], 可得M 式中:M 3 旋转软化效应 某悬臂梁长1m, 宽0.1m, 高0.05m (见图2) 。材料参数见表1。试按两种情况分析梁的弯曲振动特性: (1) 梁在x-y平面绕z轴转动; (2) 梁在z-x平面绕y轴转动。 采用有限元梁单元模型, 将绕中性轴z轴﹑y轴的截面惯性参数分别代入式 (12) , 在Matlab平台上进行数值模拟。为便于比较梁在发生旋转软化效应时的弯曲振动特性, 计算分两步: (1) 仅设置应力刚化效应计算弯振特性; (2) 同时设置应力刚化及旋转软化效应计算弯振特性 (最终结果) 。参考图2所示坐标系, 通过数值模拟得到梁在x-y平面内绕z轴转动时的Campbell图, 见图3。图3中仅列出梁的前四阶弯曲振动模态, 其中第一阶﹑第三阶为梁在z-x平面内的弯振模态, 第二阶﹑第四阶为x-y平面内的弯振模态。 从图3可以看出:旋转软化效应对梁振动模态的影响与其旋转平面方位有关。当梁在x-y平面绕z轴转动时, 在旋转软化效应作用下, 位于x-y平面的弯振频率 (图中第二阶及第四阶摆振) 随转速增大而下降 (相对于仅设置应力刚化效应) 。从图中可以看出:仅设置应力刚化效应与同时设置应力刚化及旋转软化效应得到同一阶弯振曲线, 随转速增大而分开 (见图中“分开”) ;x-y平面内第一阶 (图中第二阶) 频率曲线下降幅度大于第二阶 (图中第四阶) 频率曲线的下降幅度, 表明旋转软化效应对低价模态频率的影响较大;而位于z-x平面的第一阶、第二阶 (图中第一阶、第三阶) 频率曲线没有分开, 图中用“重合”表示, 可见:仅设置应力刚化效应与同时设置应力刚化及旋转软化效应得到的频率曲线随转速增大保持重合, 这表明旋转软化效应对其振动模态没有影响。梁在x-y平面绕z轴转动, 则位于z-x平面内的弯曲振动模态不受旋转软化效应的影响。 为了验证采用梁单元模型所得计算结果的准确性, 将本文梁单元模型与Ansys10.0实体单元模型计算结果作了比较。计算中, 同时设置应力刚化及旋转软化效应, 并设定梁在x-y平面绕z轴转动。图4计算结果表明:梁单元模型与实体单元模型计算结果误差不大于5%, 验证了本文梁单元模型求解技术的准确性。 为进一步表明旋转软化效应与梁旋转