用高阶梁单元计算结构附带集中质量的灵敏度.docVIP

振　动　与　冲　击 第 ３２卷第 １５期 ＪＵＲＮＡＬＯＦＶＩＲＡＴＩＮＡＮＤＳＨＯＣＫ Ｖｏｌ３２ Ｎｏ．１５２０１３　  用高阶梁单元计算结构附带集中质量的灵敏度 王　栋，马建军 （西北工业大学 航空学院，西安　７１００７２） 　　摘　要：采用高阶有限梁单元，研究梁或刚架结构附带集中质量的固有频率及其灵敏度计算问题。通过在节点上 增加曲率和曲率的一阶导数两个独立的自由度，可以直接获得频率灵敏度计算所需的基本量。计算结果表明：用高阶梁 单元，可以较精确地获得固有频率及其对集中质量位置的灵敏度值。此外，经过比较发现，在单元网格划分比较密集的情 形下，采用通常的 ２节点 ４自由度梁单元，也能获得比较满意的频率灵敏度计算结果。 关键词：高阶梁单元；梁 －质量系统；集中质量位置；频率灵敏度 中图分类号：Ｏ３２６　　　文献标识码：Ａ Ｓｅｎｓｔｖｉｙｃａｌｕｌｔｏｎｏｆｂｅａｍ ａｓ ｓｒｕｃｔｒｅｕｓｎｇｈｉｈｅｒｒｄｅｒｂｅａｍ ｅｌｍｅｎｔ ＷＡＮＧＤｏｎｇ，ＭＡＪｉｎｕｎ （ＳｃｈｏｏｌｏｆＡｅｒｎａｕｔｃｓ Ｎｏｒｈｗｅｓｅｒ ＰｏｌｔｃｈｎｉａｌＵｎｉｅｒｉｙ Ｘｉａｎ７１００７２，Ｃｈｉａ 　　Ａｂｓｒａｃｔ　Ｔｈｅｃａｌｕｌｔｏｎｏｆｄｙｎａｍｉｂｅｈａｖｉｒｏｆａｂｅａｍｓｒｃｔｒ ｗｉｈａｌｍｐｅｄｍａｓ ｗｉｈｂｏｔ ｔａｎｓａｔｏｎａｌａｎｄ ｒｔｒ ｉｅｒｉｓｗａｓｓｕｄｉｄｕｓｎｇｔｏｔｐｅｓｏｆｈｉｈｅｒｏｒｅｒｂｅａｍｅｌｍｅｎｔ．Ｂｙａｄｄｉｇｒｓｅｃｔｖｅｌｉｄｅｐｅｎｄｅｎｔｖａｒａｂｌｓｏｆ ｔｅｃｕｒａｔｒ ａｎｄｉｓｆｒｔｏｒｅｒｄｅｒｖａｔｖｅａｔａｎｏｄｅ ｔｅｆｅｑｕｅｎｃｙａｎｄｉｓｓｎｓｔｖｉｙｗｉｈｒｓｅｃｔｔ ｔｅｐｏｓｔｏｎｏｆｔｅ ｌｍｐｅｄｍａｓ ａｔａｃｈｅｄｗｅｒ ｃａｌｕｌｔｄａｃｃｕｒｔｌ ｕｓｎｇｔｅｈｉｈｅｒｏｒｅｒｂｅａｍ ｅｌｍｅｎｔ．Ｎｕｍｅｒｃａｌｒｓｌｓｓｏｗｔａｔｔｅ ｕｓ ｏｆｇｅｎｅｒｃｂｅａｍ ｅｌｍｅｎｔｏｆｔｅｃｕｂｉ ｐｏｌｎｏｍｉｌｄｉｐｌｃｅｍｅｎｔｃａｎａｌｏａｃｈｉｖｅｔｅｆｅｑｕｅｎｃｙｓｎｓｔｖｉｙｓｔｓａｃｔｒｌ ｗｉｈａｄｅｎｓ ｆｎｉｅｅｌｍｅｎｔｍｅｓ． Ｋｅｙｗｏｒｄｓ ｈｉｈｅｒｏｒｅｒｂｅａｍ ｅｌｍｅｎｔ ｌｍｐｅｄ ｍａｓ ｐｏｓｔｏｎ； ｆｅｑｕｅｎｃｙ ｓｎｓｔｖｉｙ ｒｔｒ ｉｅｒｉ； ｂｅａｍ ｍａｓ ｓｓｅｍ  　　工程中的许多结构，可以等效成 梁 附 带 集 中 质 量 的力学模型。在结构振动分析、优化设计、模型修改以 及参数识别过程中，结构的固有频率 相 对 于 集 中 质 量 位置的灵敏度（即一阶导数）值，能提供非常重要的指  此，只能通过近似估算，获得节点上的曲率值。 在分析复杂梁结 构 的 固 有 频 率 时，Ｔｏ［５］曾 经 运 用 了两种高阶梁单元：２节点 ６自由度单元和 ２节点 ８自 由度单元。即每个节点的自由度是：横向位移、转角和 导性 信 息 ［１－３］ 曲率，甚至是曲率的 一 阶 导 数。与 传 统 的 ２节 点 ４自 ｍｂｅｒ’ｓｐｒｎｃｉｌ），由惯性力与约束反力等效的 原 则， 得到了梁结构的固有频率，相对于集 中 质 量 位 置 的 灵 敏度计算公式，且同时考虑了集中质量的平动 （移）惯 性和转动惯性。数值计算 结 果 表 明，集 中 质 量 的 转 动 惯性，对结构的高阶频率及其灵敏度值影响很大，应该 引起高度重视。然而，要计 算 梁 －质 量 系 统 的 固 有 频 率对集中质量位置的灵敏度，必须首 先 获 得 集 中 质 量 所在处（点）的曲率值。在简单梁结构的理论分析时， 只要对振型函数两次求导即可获得 ［４］。但对于复杂结 构，求解频率和振型都是采用数值方法，如有限元方法 等。而通常的 ２节点 ４自由度单元 （每个节点只有横 向位移 ｗ和转角 θ二个自由度），并无曲率自由度。因 基金项目：２０１１年西北工业大学本科毕业论文重点扶持项目资助 收稿日期：２０１２－０１－３０　修改稿收到日期：２０１２－０６－２１ 第一作者 王　栋 男，教授，博士生导师，１９６２年生 由度单元相比 ［６］，高阶梁单元位移模式阶次更高，能更 加准确地描述单元内的位移状况。因此在计算结构的 静、动态特性时，具有精度高和收敛速度快的特点。现 已广泛应用于复合材料梯形梁的固有频率和失稳载荷 计算问题 ［７－８］。本 文 探 索 利 用 有 限 元 方 法，计 算 梁 － 质量系统的固有频率，相对于集中质 量 位 置 的 灵 敏 度 的可行性和精度问题。采 用 以 上 两 种 高 阶 梁 单 元，分 析、计算附带有集中质量的梁、或刚（框）架结构的固有 频率及其一阶灵敏度值